@article{YuanShengPreicketal.2020,
  author    = {Yuan, Junxia and Sheng, Guilian and Preick, Michaela and Sun, Boyang and Hou, Xindong and Chen, Shungang and Taron, Ulrike Helene and Barlow, Axel and Wang, Linying and Hu, Jiaming and Deng, Tao and Lai, Xulong and Hofreiter, Michael},
  title     = {Mitochondrial genomes of Late Pleistocene caballine horses from China belong to a separate clade},
  series = {Quaternary science reviews : the international multidisciplinary research and review journal},
  volume    = {250},
  journal   = {Quaternary science reviews : the international multidisciplinary research and review journal},
  publisher = {Elsevier},
  address   = {Amsterdam [u.a.]},
  issn      = {0277-3791},
  doi       = {10.1016/j.quascirev.2020.106691},
  pages     = {8},
  year      = {2020},
  abstract  = {There were several species of Equus in northern China during the Late Pleistocene, including Equus przewalskii and Equus dalianensis. A number of morphological studies have been carried out on E. przewalskii and E. dalianensis, but their evolutionary history is still unresolved. In this study, we retrieved near-complete mitochondrial genomes from E. dalianensis and E. przewalskii specimens excavated from Late Pleistocene strata in northeastern China. Phylogenetic analyses revealed that caballoid horses were divided into two subclades: the New World and the Old World caballine horse subclades. The Old World caballine horses comprise of two deep phylogenetic lineages, with modern and ancient Equus caballus and modern E. przewalskii forming lineage I, and the individuals in this study together with one Yakut specimen forming lineage II. Our results indicate that Chinese Late Pleistocene caballoid horses showed a closer relationship to other Eurasian caballine horses than that to Pleistocene horses from North America. In addition, phylogenetic analyses suggested a close relationship between E. dalianensis and the Chinese fossil E. przewalskii, in agreement with previous researches based on morphological analyses. Interestingly, E. dalianensis and the fossil E. przewalskii were intermixed rather than split into distinct lineages, suggesting either that gene flow existed between these two species or that morphology-based species assignment of palaeontological specimens is not always correct. Moreover, Bayesian analysis showed that the divergence time between the New World and the Old World caballoid horses was at 1.02 Ma (95\% CI: 0.86-1.24 Ma), and the two Old World lineages (I \& II) split at 0.88 Ma (95\% CI: 0.69-1.13 Ma), which indicates that caballoid horses seem to have evolved into different populations in the Old World soon after they migrated from North America via the Bering Land Bridge. Finally, the TMRCA of E. dalianensis was estimated at 0.20 Ma (95\% CI: 0.15-0.28 Ma), and it showed a relative low genetic diversity compared with other Equus species.},
  language  = {en}
}
@article{XiangGaoYuetal.2014,
  author    = {Xiang, Hai and Gao, Jianqiang and Yu, Baoquan and Zhou, Hui and Cai, Dawei and Zhang, Youwen and Chen, Xiaoyong and Wang, Xi and Hofreiter, Michael and Zhao, Xingbo},
  title     = {Early Holocene chicken domestication in northern China},
  series = {Proceedings of the National Academy of Sciences of the United States of America},
  volume    = {111},
  journal   = {Proceedings of the National Academy of Sciences of the United States of America},
  number    = {49},
  publisher = {National Acad. of Sciences},
  address   = {Washington},
  issn      = {0027-8424},
  doi       = {10.1073/pnas.1411882111},
  pages     = {17564 -- 17569},
  year      = {2014},
  abstract  = {Chickens represent by far the most important poultry species, yet the number, locations, and timings of their domestication have remained controversial for more than a century. Here we report ancient mitochondrial DNA sequences from the earliest archaeological chicken bones from China, dating back to similar to 10,000 B.P. The results clearly show that all investigated bones, including the oldest from the Nanzhuangtou site, are derived from the genus Gallus, rather than any other related genus, such as Phasianus. Our analyses also suggest that northern China represents one region of the earliest chicken domestication, possibly dating as early as 10,000 y B.P. Similar to the evidence from pig domestication, our results suggest that these early domesticated chickens contributed to the gene pool of modern chicken populations. Moreover, our results support the idea that multiple members of the genus Gallus, specifically Gallus gallus and Gallus sonneratii contributed to the gene pool of the modern domestic chicken. Our results provide further support for the growing evidence of an early mixed agricultural complex in northern China.},
  language  = {en}
}
@phdthesis{Wutke2016,
  author    = {Wutke, Saskia},
  title     = {Tracing Changes in Space and Time},
  school      = {Universit{\"a}t Potsdam},
  pages     = {x, 84},
  year      = {2016},
  abstract  = {The horse is a fascinating animal symbolizing power, beauty, strength and grace. Among all the animal species domesticated the horse had the largest impact on the course of human history due to its importance for warfare and transportation. Studying the process of horse domestication contributes to the knowledge about the history of horses and even of our own species. Research based on molecular methods has increasingly focused on the genetic basis of horse domestication. Mitochondrial DNA (mtDNA) analyses of modern and ancient horses detected immense maternal diversity, probably due to many mares that contributed to the domestic population. However, mtDNA does not provide an informative phylogeographic structure. In contrast, Y chromosome analyses displayed almost complete uniformity in modern stallions but relatively high diversity in a few ancient horses. Further molecular markers that seem to be well suited to infer the domestication history of horses or genetic and phenotypic changes during this process are loci associated with phenotypic traits. This doctoral thesis consists of three different parts for which I analyzed various single nucleotide polymorphisms (SNPs) associated with coat color, locomotion or Y chromosomal variation of horses. These SNPs were genotyped in 350 ancient horses from the Chalcolithic (5,000 BC) to the Middle Ages (11th century). The distribution of the samples ranges from China to the Iberian Peninsula and Iceland. By applying multiplexed next-generation sequencing (NGS) I sequenced short amplicons covering the relevant positions: i) eight coat-color-associated mutations in six genes to deduce the coat color phenotype; ii) the so-called 'Gait-keeper' SNP in the DMRT3 gene to screen for the ability to amble; iii) 16 SNPs previously detected in ancient horses to infer the corresponding haplotype. Based on these data I investigated the occurrence and frequencies of alleles underlying the respective phenotypes as well as Y chromosome haplotypes at different times and regions. Also, selection coefficients for several Y chromosome lineages or phenotypes were estimated. Concerning coat color differences in ancient horses my work constitutes the most comprehensive study to date. I detected an increase of chestnut horses in the Middle Ages as well as differential selection for spotted and solid phenotypes over time which reflects changing human preferences. With regard to ambling horses, the corresponding allele was present in medieval English and Icelandic horses. Based on these results I argue that Norse settlers, who frequently invaded parts of Britain, brought ambling individuals to Iceland from the British Isles which can be regarded the origin of this trait. Moreover, these settlers appear to have selected for ambling in Icelandic horses. Relating to the third trait, the paternal diversity, these findings represent the largest ancient dataset of Y chromosome variation in non-humans. I proved the existence of several Y chromosome haplotypes in early domestic horses. The decline of Y chromosome variation coincides with the movement of nomadic peoples from the Eurasian steppes and later with different breeding practices in the Roman period. In conclusion, positive selection was estimated for several phenotypes/lineages in different regions or times which indicates that these were preferred by humans. Furthermore, I could successfully infer the distribution and dispersal of horses in association with human movements and actions. Thereby, a better understanding of the influence of people on the changing appearance and genetic diversity of domestic horses could be gained. My results also emphasize the close relationship of ancient genetics and archeology or history and that only in combination well-founded conclusions can be reached.},
  language  = {en}
}
@misc{WestburyBalekaBarlowetal.2017,
  author    = {Westbury, Michael V. and Baleka, Sina Isabelle and Barlow, Axel and Hartmann, Stefanie and Paijmans, Johanna L. A. and Kramarz, Alejandro and Forasiepi, Anal{\´i}a M. and Bond, Mariano and Gelfo, Javier N. and Reguero, Marcelo A. and L{\´o}pez-Mendoza, Patricio and Taglioretti, Matias and Scaglia, Fernando and Rinderknecht, Andr{\´e}s and Jones, Washington and Mena, Francisco and Billet, Guillaume and de Muizon, Christian and Aguilar, Jos{\´e} Luis and MacPhee, Ross D.E. and Hofreiter, Michael},
  title     = {A mitogenomic timetree for Darwin's enigmatic South American mammal Macrauchenia patachonica},
  series = {Postprints der Universit{\"a}t Potsdam Mathematisch-Naturwissenschaftliche Reihe},
  journal   = {Postprints der Universit{\"a}t Potsdam Mathematisch-Naturwissenschaftliche Reihe},
  number    = {793},
  issn      = {1866-8372},
  doi       = {10.25932/publishup-44080},
  url       = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-440801},
  pages     = {8},
  year      = {2017},
  abstract  = {The unusual mix of morphological traits displayed by extinct South American native ungulates (SANUs) confounded both Charles Darwin, who first discovered them, and Richard Owen, who tried to resolve their relationships. Here we report an almost complete mitochondrial genome for the litoptern Macrauchenia. Our dated phylogenetic tree places Macrauchenia as sister to Perissodactyla, but close to the radiation of major lineages within Laurasiatheria. This position is consistent with a divergence estimate of B66Ma (95\% credibility interval, 56.64-77.83 Ma) obtained for the split between Macrauchenia and other Panperissodactyla. Combined with their morphological distinctiveness, this evidence supports the positioning of Litopterna (possibly in company with other SANU groups) as a separate order within Laurasiatheria. We also show that, when using strict criteria, extinct taxa marked by deep divergence times and a lack of close living relatives may still be amenable to palaeogenomic analysis through iterative mapping against more distant relatives.},
  language  = {en}
}
@article{ThomasCarvalhoHaileetal.2017,
  author    = {Thomas, Jessica E. and Carvalho, Gary R. and Haile, James and Martin, Michael D. and Castruita, Jose A. Samaniego and Niemann, Jonas and Sinding, Mikkel-Holger S. and Sandoval-Velasco, Marcela and Rawlence, Nicolas J. and Fuller, Errol and Fjeldsa, Jon and Hofreiter, Michael and Stewart, John R. and Gilbert, M. Thomas P. and Knapp, Michael},
  title     = {An ‛Aukward' tale},
  series = {Genes},
  volume    = {8},
  journal   = {Genes},
  number    = {6},
  publisher = {MDPI},
  address   = {Basel},
  issn      = {2073-4425},
  doi       = {10.3390/genes8060164},
  pages     = {164},
  year      = {2017},
  abstract  = {One hundred and seventy-three years ago, the last two Great Auks, Pinguinus impennis, ever reliably seen were killed. Their internal organs can be found in the collections of the Natural History Museum of Denmark, but the location of their skins has remained a mystery. In 1999, Great Auk expert Errol Fuller proposed a list of five potential candidate skins in museums around the world. Here we take a palaeogenomic approach to test which—if any—of Fuller's candidate skins likely belong to either of the two birds. Using mitochondrial genomes from the five candidate birds (housed in museums in Bremen, Brussels, Kiel, Los Angeles, and Oldenburg) and the organs of the last two known individuals, we partially solve the mystery that has been on Great Auk scholars' minds for generations and make new suggestions as to the whereabouts of the still-missing skin from these two birds.},
  language  = {en}
}
@article{TeasdalevanDoornFiddymentetal.2015,
  author    = {Teasdale, Matthew David and van Doorn, N. L. and Fiddyment, S. and Webb, C. C. and Hofreiter, Michael and Collins, Matthew J. and Bradley, Daniel G.},
  title     = {Paging through history: parchment as a reservoir of ancient DNA for next generation sequencing},
  series = {Philosophical transactions of the Royal Society of London : B, Biological sciences},
  volume    = {370},
  journal   = {Philosophical transactions of the Royal Society of London : B, Biological sciences},
  number    = {1660},
  publisher = {Royal Society},
  address   = {London},
  issn      = {0962-8436},
  doi       = {10.1098/rstb.2013.0379},
  pages     = {7},
  year      = {2015},
  abstract  = {Parchment represents an invaluable cultural reservoir. Retrieving an additional layer of information from these abundant, dated livestock-skins via the use of ancient DNA (aDNA) sequencing has been mooted by a number of researchers. However, prior PCR-based work has indicated that this may be challenged by cross-individual and cross-species contamination, perhaps from the bulk parchment preparation process. Here we apply next generation sequencing to two parchments of seventeenth and eighteenth century northern English provenance. Following alignment to the published sheep, goat, cow and human genomes, it is clear that the only genome displaying substantial unique homology is sheep and this species identification is confirmed by collagen peptide mass spectrometry. Only 4\% of sequence reads align preferentially to a different species indicating low contamination across species. Moreover, mitochondrial DNA sequences suggest an upper bound of contamination at 5\%. Over 45\% of reads aligned to the sheep genome, and even this limited sequencing exercise yield 9 and 7\% of each sampled sheep genome post filtering, allowing the mapping of genetic affinity to modern British sheep breeds. We conclude that parchment represents an excellent substrate for genomic analyses of historical livestock.},
  language  = {en}
}
@article{TassiVaiGhirottoetal.2017,
  author    = {Tassi, Francesca and Vai, Stefania and Ghirotto, Silvia and Lari, Martina and Modi, Alessandra and Pilli, Elena and Brunelli, Andrea and Susca, Roberta Rosa and Budnik, Alicja and Labuda, Damian and Alberti, Federica and Lalueza-Fox, Carles and Reich, David and Caramelli, David and Barbujani, Guido},
  title     = {Genome diversity in the Neolithic Globular Amphorae culture and the spread of Indo-European languages},
  series = {Proceedings of the Royal Society of London : B, Biological sciences},
  volume    = {284},
  journal   = {Proceedings of the Royal Society of London : B, Biological sciences},
  publisher = {Royal Society},
  address   = {London},
  issn      = {0962-8452},
  doi       = {10.1098/rspb.2017.1540},
  pages     = {9},
  year      = {2017},
  abstract  = {It is unclear whether Indo-European languages in Europe spread from the Pontic steppes in the late Neolithic, or from Anatolia in the Early Neolithic. Under the former hypothesis, people of the Globular Amphorae culture (GAC) would be descended from Eastern ancestors, likely representing the Yamnaya culture. However, nuclear (six individuals typed for 597 573 SNPs) and mitochondrial (11 complete sequences) DNA from the GAC appear closer to those of earlier Neolithic groups than to the DNA of all other populations related to the Pontic steppe migration. Explicit comparisons of alternative demographic models via approximate Bayesian computation confirmed this pattern. These results are not in contrast to Late Neolithic gene flow from the Pontic steppes into Central Europe. However, they add nuance to this model, showing that the eastern affinities of the GAC in the archaeological record reflect cultural influences from other groups from the East, rather than the movement of people.},
  language  = {en}
}
@misc{TaronLellBarlowetal.2018,
  author    = {Taron, Ulrike H. and Lell, Moritz and Barlow, Axel and Paijmans, Johanna L. A.},
  title     = {Testing of Alignment Parameters for Ancient Samples},
  series = {Genes},
  journal   = {Genes},
  url       = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-409683},
  pages     = {12},
  year      = {2018},
  abstract  = {High-throughput sequence data retrieved from ancient or other degraded samples has led to unprecedented insights into the evolutionary history of many species, but the analysis of such sequences also poses specific computational challenges. The most commonly used approach involves mapping sequence reads to a reference genome. However, this process becomes increasingly challenging with an elevated genetic distance between target and reference or with the presence of contaminant sequences with high sequence similarity to the target species. The evaluation and testing of mapping efficiency and stringency are thus paramount for the reliable identification and analysis of ancient sequences. In this paper, we present 'TAPAS', (Testing of Alignment Parameters for Ancient Samples), a computational tool that enables the systematic testing of mapping tools for ancient data by simulating sequence data reflecting the properties of an ancient dataset and performing test runs using the mapping software and parameter settings of interest. We showcase TAPAS by using it to assess and improve mapping strategy for a degraded sample from a banded linsang (Prionodon linsang), for which no closely related reference is currently available. This enables a 1.8-fold increase of the number of mapped reads without sacrificing mapping specificity. The increase of mapped reads effectively reduces the need for additional sequencing, thus making more economical use of time, resources, and sample material.},
  language  = {en}
}
@article{TaronLellBarlowetal.2018,
  author    = {Taron, Ulrike H. and Lell, Moritz and Barlow, Axel and Paijmans, Johanna L. A.},
  title     = {Testing of Alignment Parameters for Ancient Samples},
  series = {Genes},
  volume    = {9},
  journal   = {Genes},
  number    = {3},
  publisher = {Molecular Diversity Preservation International},
  address   = {Basel},
  issn      = {2073-4425},
  doi       = {10.3390/genes9030157},
  pages     = {1 -- 12},
  year      = {2018},
  abstract  = {High-throughput sequence data retrieved from ancient or other degraded samples has led to unprecedented insights into the evolutionary history of many species, but the analysis of such sequences also poses specific computational challenges. The most commonly used approach involves mapping sequence reads to a reference genome. However, this process becomes increasingly challenging with an elevated genetic distance between target and reference or with the presence of contaminant sequences with high sequence similarity to the target species. The evaluation and testing of mapping efficiency and stringency are thus paramount for the reliable identification and analysis of ancient sequences. In this paper, we present 'TAPAS', (Testing of Alignment Parameters for Ancient Samples), a computational tool that enables the systematic testing of mapping tools for ancient data by simulating sequence data reflecting the properties of an ancient dataset and performing test runs using the mapping software and parameter settings of interest. We showcase TAPAS by using it to assess and improve mapping strategy for a degraded sample from a banded linsang (Prionodon linsang), for which no closely related reference is currently available. This enables a 1.8-fold increase of the number of mapped reads without sacrificing mapping specificity. The increase of mapped reads effectively reduces the need for additional sequencing, thus making more economical use of time, resources, and sample material.},
  language  = {en}
}
@article{TaronLellBarlowetal.2018,
  author    = {Taron, Ulrike H. and Lell, Moritz and Barlow, Axel and Paijmans, Johanna L. A.},
  title     = {Testing of Alignment Parameters for Ancient Samples},
  series = {Genese},
  volume    = {9},
  journal   = {Genese},
  number    = {3},
  publisher = {MDPI},
  address   = {Basel},
  issn      = {2073-4425},
  doi       = {10.3390/genes9030157},
  pages     = {12},
  year      = {2018},
  abstract  = {High-throughput sequence data retrieved from ancient or other degraded samples has led to unprecedented insights into the evolutionary history of many species, but the analysis of such sequences also poses specific computational challenges. The most commonly used approach involves mapping sequence reads to a reference genome. However, this process becomes increasingly challenging with an elevated genetic distance between target and reference or with the presence of contaminant sequences with high sequence similarity to the target species. The evaluation and testing of mapping efficiency and stringency are thus paramount for the reliable identification and analysis of ancient sequences. In this paper, we present 'TAPAS', (Testing of Alignment Parameters for Ancient Samples), a computational tool that enables the systematic testing of mapping tools for ancient data by simulating sequence data reflecting the properties of an ancient dataset and performing test runs using the mapping software and parameter settings of interest. We showcase TAPAS by using it to assess and improve mapping strategy for a degraded sample from a banded linsang (Prionodon linsang), for which no closely related reference is currently available. This enables a 1.8-fold increase of the number of mapped reads without sacrificing mapping specificity. The increase of mapped reads effectively reduces the need for additional sequencing, thus making more economical use of time, resources, and sample material.},
  language  = {en}
}
@misc{SpellervandenHurkCharpentieretal.2016,
  author    = {Speller, Camilla and van den Hurk, Youri and Charpentier, Anne and Rodrigues, Ana and Gardeisen, Armelle and Wilkens, Barbara and McGrath, Krista and Rowsell, Keri and Spindler, Luke and Collins, Matthew J. and Hofreiter, Michael},
  title     = {Barcoding the largest animals on Earth: ongoing challenges and molecular solutions in the taxonomic identification of ancient cetaceans},
  series = {Philosophical transactions of the Royal Society of London : B, Biological sciences},
  volume    = {371},
  journal   = {Philosophical transactions of the Royal Society of London : B, Biological sciences},
  publisher = {Royal Society},
  address   = {London},
  issn      = {0962-8436},
  doi       = {10.1098/rstb.2015.0332},
  pages     = {11},
  year      = {2016},
  language  = {en}
}
@misc{SiskaJonesJeonetal.2017,
  author    = {Siska, Veronika and Jones, Eppie Ruth and Jeon, Sungwon and Bhak, Youngjune and Kim, Hak-Min and Cho, Yun Sung and Kim, Hyunho and Lee, Kyusang and Veselovskaya, Elizaveta and Balueva, Tatiana and Gallego-Llorente, Marcos and Hofreiter, Michael and Bradley, Daniel G. and Eriksson, Anders and Pinhasi, Ron and Bhak, Jong and Manica, Andrea},
  title     = {Genome-wide data from two early Neolithic East Asian individuals dating to 7700 years ago},
  series = {Postprints der Universit{\"a}t Potsdam Mathematisch-Naturwissenschaftliche Reihe},
  journal   = {Postprints der Universit{\"a}t Potsdam Mathematisch-Naturwissenschaftliche Reihe},
  number    = {791},
  issn      = {1866-8372},
  doi       = {10.25932/publishup-43997},
  url       = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-439977},
  pages     = {11},
  year      = {2017},
  abstract  = {Ancient genomes have revolutionized our understanding of Holocene prehistory and, particularly, the Neolithic transition in western Eurasia. In contrast, East Asia has so far received little attention, despite representing a core region at which the Neolithic transition took place independently similar to 3 millennia after its onset in the Near East. We report genome-wide data from two hunter-gatherers from Devil's Gate, an early Neolithic cave site (dated to similar to 7.7 thousand years ago) located in East Asia, on the border between Russia and Korea. Both of these individuals are genetically most similar to geographically close modern populations from the Amur Basin, all speaking Tungusic languages, and, in particular, to the Ulchi. The similarity to nearby modern populations and the low levels of additional genetic material in the Ulchi imply a high level of genetic continuity in this region during the Holocene, a pattern that markedly contrasts with that reported for Europe.},
  language  = {en}
}
@article{SignorePaijmansHofreiteretal.2019,
  author    = {Signore, Anthony V. and Paijmans, Johanna L. A. and Hofreiter, Michael and Fago, Angela and Weber, Roy E. and Springer, Mark S. and Campbell, Kevin L.},
  title     = {Emergence of a chimeric globin pseudogene and increased Hemoglobin Oxygen Affinity Underlie the evolution of aquatic specializations in Sirenia},
  series = {Molecular biology and evolution},
  volume    = {36},
  journal   = {Molecular biology and evolution},
  number    = {6},
  publisher = {Oxford Univ. Press},
  address   = {Oxford},
  issn      = {0737-4038},
  doi       = {10.1093/molbev/msz044},
  pages     = {1134 -- 1147},
  year      = {2019},
  abstract  = {As limits on O2 availability during submergence impose severe constraints on aerobic respiration, the oxygen binding globin proteins of marine mammals are expected to have evolved under strong evolutionary pressures during their land-to-sea transition. Here, we address this question for the order Sirenia by retrieving, annotating, and performing detailed selection analyses on the globin repertoire of the extinct Steller's sea cow (Hydrodamalis gigas), dugong (Dugong dugon), and Florida manatee (Trichechus manatus latirostris) in relation to their closest living terrestrial relatives (elephants and hyraxes). These analyses indicate most loci experienced elevated nucleotide substitution rates during their transition to a fully aquatic lifestyle. While most of these genes evolved under neutrality or strong purifying selection, the rate of nonsynonymous/synonymous replacements increased in two genes (Hbz-T1 and Hba-T1) that encode the α-type chains of hemoglobin (Hb) during each stage of life. Notably, the relaxed evolution of Hba-T1 is temporally coupled with the emergence of a chimeric pseudogene (Hba-T2/Hbq-ps) that contributed to the tandemly linked Hba-T1 of stem sirenians via interparalog gene conversion. Functional tests on recombinant Hb proteins from extant and ancestral sirenians further revealed that the molecular remodeling of Hba-T1 coincided with increased Hb-O2 affinity in early sirenians. Available evidence suggests that this trait evolved to maximize O2 extraction from finite lung stores and suppress tissue O2 offloading, thereby facilitating the low metabolic intensities of extant sirenians. In contrast, the derived reduction in Hb-O2 affinity in (sub)Arctic Steller's sea cows is consistent with fueling increased thermogenesis by these once colossal marine herbivores.},
  language  = {en}
}
@phdthesis{Schulte2022,
  author    = {Schulte, Luise},
  title     = {Dynamics of Larix (Mill.) species in Siberia during the last 50,000 years inferred from sedimentary ancient DNA},
  doi       = {10.25932/publishup-55878},
  url       = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-558782},
  school      = {Universit{\"a}t Potsdam},
  pages     = {xi, 121},
  year      = {2022},
  abstract  = {The deciduous needle tree larch (Larix Mill.) covers more than 80\% of the Asian boreal forests. Only a few Larix species constitute the vast forests and these species differ markedly in their ecological traits, most importantly in their ability to grow on and stabilize underlying permafrost. The pronounced dominance of the summergreen larches makes the Asian boreal forests unique, as the rest of the northern hemisphere boreal forests is almost exclusively dominated by evergreen needle-leaf forests. Global warming is impacting the whole world but is especially pronounced in the arctic and boreal regions. Although adapted to extreme climatic conditions, larch forests are sensitive to varying climatic conditions. By their sheer size, changes in Asian larch forests as range shifts or changes in species composition and the resulting vegetation-climate feedbacks are of global relevance. It is however still uncertain if larch forests will persist under the ongoing warming climate or if they will be replaced by evergreen forests. It is therefore of great importance to understand how these ecosystems will react to future climate warmings and if they will maintain their dominance. One step in the better understanding of larch dynamics is to study how the vast dominant forests developed and why they only established in northern Asia. A second step is to study how the species reacted to past changes in the climate. The first objective of this thesis was to review and identify factors promoting Asian larch dominance. I achieved this by synthesizing and comparing reported larch occurrences and influencing components on the northern hemisphere continents in the present and in the past. The second objective was to find a possibility to directly study past Larix populations in Siberia and specifically their genetic variation, enabling the study of geographic movements. For this, I established chloroplast enrichment by hybridization capture from sedimentary ancient DNA (sedaDNA) isolated from lake sediment records. The third objective was to use the established method to track past larch populations, their glacial refugia during the Last Glacial Maximum (LGM) around 21,000 years before present (ka BP), and their post-glacial migration patterns. To study larch promoting factors, I compared the present state of larch species ranges, areas of dominance, their bioclimatic niches, and the distribution on different extents and thaw depths of permafrost. The species comparison showed that the bioclimatic niches greatly overlap between the American and Asian species and that it is only in the extremely continental climates in which only the Asian larch species can persist. I revealed that the area of dominance is strongly connected to permafrost extent but less linked to permafrost seasonal thaw depths. Comparisons of the paleorecord of larch between the continents suggest differences in the recolonization history. Outside of northern Asia and Alaska, glacial refugial populations of larch were confined to the southern regions and thus recolonization could only occur as migration from south to north. Alaskan larch populations could not establish wide-range dominant forest which could be related to their own genetically depletion as separated refugial population. In Asia, it is still unclear whether or not the northern refugial populations contributed and enhanced the postglacial colonization or whether they were replaced by populations invading from the south in the course of climate warming. Asian larch dominance is thus promoted partly by adaptions to extremely continental climates and by adaptations to grow on continuous permafrost but could be also connected to differences in glacial survival and recolonization history of Larix species. Except for extremely rare macrofossil findings of fossilized cones, traditional methods to study past vegetation are not able to distinguish between larch species or populations. Within the scope of this thesis, I therefore established a method to retrieve genetic information of past larch populations to distinguish between species. Using the Larix chloroplast genome as target, I successfully applied the method of DNA target enrichment by hybridization capture on sedaDNA samples from lake records and showed that it is able to distinguish between larch species. I then used the method on samples from lake records from across Siberia dating back up to 50 ka BP. The results allowed me to address the question of glacial survival and post-glacial recolonization mode in Siberian larch species. The analyzed pattern showed that LGM refugia were almost exclusively constituted by L. gmelinii, even in sites of current L. sibirica distribution. For included study sites, L. sibirica migrated into its extant northern distribution area only in the Holocene. Consequently, the post-glacial recolonization of L. sibirica was not enhanced by northern glacial refugia. In case of sites in extant distribution area of L. gmelinii, the absence of a genetic turn-over point to a continuous population rather than an invasion of southern refugia. The results suggest that climate has a strong influence on the distribution of Larix species and that species may also respond differently to future climate warming. Because species differ in their ecological characteristics, species distribution is also relevant with respect to further feedbacks between vegetation and climate. With this thesis, I give an overview of present and past larch occurrences and evaluate which factors promote their dominance. Furthermore, I provide the tools to study past Larix species and give first important insights into the glacial history of Larix populations.},
  language  = {en}
}
@misc{MohandesanSpellerPetersetal.2017,
  author    = {Mohandesan, Elmira and Speller, Camilla F. and Peters, Joris and Uerpmann, Hans-Peter and Uerpmann, Margarethe and De Cupere, Bea and Hofreiter, Michael and Burger, Pamela A.},
  title     = {Combined hybridization capture and shotgun sequencing for ancient DNA analysis of extinct wild and domestic dromedary camel},
  series = {Postprints der Universit{\"a}t Potsdam Mathematisch-Naturwissenschaftliche Reihe},
  journal   = {Postprints der Universit{\"a}t Potsdam Mathematisch-Naturwissenschaftliche Reihe},
  number    = {789},
  issn      = {1866-8372},
  doi       = {10.25932/publishup-43995},
  url       = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-439955},
  pages     = {300 -- 313},
  year      = {2017},
  abstract  = {The performance of hybridization capture combined with next-generation sequencing (NGS) has seen limited investigation with samples from hot and arid regions until now. We applied hybridization capture and shotgun sequencing to recover DNA sequences from bone specimens of ancient-domestic dromedary (Camelus dromedarius) and its extinct ancestor, the wild dromedary from Jordan, Syria, Turkey and the Arabian Peninsula, respectively. Our results show that hybridization capture increased the percentage of mitochondrial DNA (mtDNA) recovery by an average 187-fold and in some cases yielded virtually complete mitochondrial (mt) genomes at multifold coverage in a single capture experiment. Furthermore, we tested the effect of hybridization temperature and time by using a touchdown approach on a limited number of samples. We observed no significant difference in the number of unique dromedary mtDNA reads retrieved with the standard capture compared to the touchdown method. In total, we obtained 14 partial mitochondrial genomes from ancient-domestic dromedaries with 17-95\% length coverage and 1.27-47.1-fold read depths for the covered regions. Using whole-genome shotgun sequencing, we successfully recovered endogenous dromedary nuclear DNA (nuDNA) from domestic and wild dromedary specimens with 1-1.06-fold read depths for covered regions. Our results highlight that despite recent methodological advances, obtaining ancient DNA (aDNA) from specimens recovered from hot, arid environments is still problematic. Hybridization protocols require specific optimization, and samples at the limit of DNA preservation need multiple replications of DNA extraction and hybridization capture as has been shown previously for Middle Pleistocene specimens.},
  language  = {en}
}
@article{MohandesanSpellerPetersetal.2017,
  author    = {Mohandesan, Elmira and Speller, Camilla F. and Peters, Joris and Uerpmann, Hans-Peter and Uerpmann, Margarethe and De Cupere, Bea and Hofreiter, Michael and Burger, Pamela A.},
  title     = {Combined hybridization capture and shotgun sequencing for ancient DNA analysis of extinct wild and domestic dromedary camel},
  series = {Molecular ecology resources},
  volume    = {17},
  journal   = {Molecular ecology resources},
  number    = {2},
  publisher = {Wiley},
  address   = {Hoboken},
  issn      = {1755-098X},
  doi       = {10.1111/1755-0998.12551},
  pages     = {300 -- 313},
  year      = {2017},
  abstract  = {The performance of hybridization capture combined with next-generation sequencing (NGS) has seen limited investigation with samples from hot and arid regions until now. We applied hybridization capture and shotgun sequencing to recover DNA sequences from bone specimens of ancient-domestic dromedary (Camelus dromedarius) and its extinct ancestor, the wild dromedary from Jordan, Syria, Turkey and the Arabian Peninsula, respectively. Our results show that hybridization capture increased the percentage of mitochondrial DNA (mtDNA) recovery by an average 187-fold and in some cases yielded virtually complete mitochondrial (mt) genomes at multifold coverage in a single capture experiment. Furthermore, we tested the effect of hybridization temperature and time by using a touchdown approach on a limited number of samples. We observed no significant difference in the number of unique dromedary mtDNA reads retrieved with the standard capture compared to the touchdown method. In total, we obtained 14 partial mitochondrial genomes from ancient-domestic dromedaries with 17-95\% length coverage and 1.27-47.1-fold read depths for the covered regions. Using whole-genome shotgun sequencing, we successfully recovered endogenous dromedary nuclear DNA (nuDNA) from domestic and wild dromedary specimens with 1-1.06-fold read depths for covered regions. Our results highlight that despite recent methodological advances, obtaining ancient DNA (aDNA) from specimens recovered from hot, arid environments is still problematic. Hybridization protocols require specific optimization, and samples at the limit of DNA preservation need multiple replications of DNA extraction and hybridization capture as has been shown previously for Middle Pleistocene specimens.},
  language  = {en}
}
@article{LudwigReissmannBeneckeetal.2015,
  author    = {Ludwig, Arne and Reissmann, Monika and Benecke, Norbert and Bellone, Rebecca and Sandoval-Castellanos, Edson and Cieslak, Michael and Gonz{\´a}lez-Fortes, Gloria M. and Morales-Muniz, Arturo and Hofreiter, Michael and Pruvost, Melanie},
  title     = {Twenty-five thousand years of fluctuating selection on leopard complex spotting and congenital night blindness in horses},
  series = {Philosophical transactions of the Royal Society of London : B, Biological sciences},
  volume    = {370},
  journal   = {Philosophical transactions of the Royal Society of London : B, Biological sciences},
  number    = {1660},
  publisher = {Royal Society},
  address   = {London},
  issn      = {0962-8436},
  doi       = {10.1098/rstb.2013.0386},
  pages     = {7},
  year      = {2015},
  abstract  = {Leopard complex spotting is inherited by the incompletely dominant locus, LP, which also causes congenital stationary night blindness in homozygous horses. We investigated an associated single nucleotide polymorphism in the TRPM1 gene in 96 archaeological bones from 31 localities from Late Pleistocene (approx. 17 000 YBP) to medieval times. The first genetic evidence of LP spotting in Europe dates back to the Pleistocene. We tested for temporal changes in the LP associated allele frequency and estimated coefficients of selection by means of approximate Bayesian computation analyses. Our results show that at least some of the observed frequency changes are congruent with shifts in artificial selection pressure for the leopard complex spotting phenotype. In early domestic horses from Kirklareli-Kanligecit (Turkey) dating to 2700-2200 BC, a remarkably high number of leopard spotted horses (six of 10 individuals) was detected including one adult homozygote. However, LP seems to have largely disappeared during the late Bronze Age, suggesting selection against this phenotype in early domestic horses. During the Iron Age, LP reappeared, probably by reintroduction into the domestic gene pool from wild animals. This picture of alternating selective regimes might explain how genetic diversity was maintained in domestic animals despite selection for specific traits at different times.},
  language  = {en}
}
@misc{HuynenSuzukiOguraetal.2014,
  author    = {Huynen, Leon and Suzuki, Takayuki and Ogura, Toshihiko and Watanabe, Yusuke and Millar, Craig D. and Hofreiter, Michael and Smith, Craig and Mirmoeini, Sara and Lambert, David M.},
  title     = {Reconstruction and in vivo analysis of the extinct tbx5 gene from ancient wingless moa (Aves: Dinornithiformes)},
  series = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe},
  journal   = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe},
  number    = {1117},
  issn      = {1866-8372},
  doi       = {10.25932/publishup-43159},
  url       = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-431599},
  pages     = {10},
  year      = {2014},
  abstract  = {Background The forelimb-specific gene tbx5 is highly conserved and essential for the development of forelimbs in zebrafish, mice, and humans. Amongst birds, a single order, Dinornithiformes, comprising the extinct wingless moa of New Zealand, are unique in having no skeletal evidence of forelimb-like structures. Results To determine the sequence of tbx5 in moa, we used a range of PCR-based techniques on ancient DNA to retrieve all nine tbx5 exons and splice sites from the giant moa, Dinornis. Moa Tbx5 is identical to chicken Tbx5 in being able to activate the downstream promotors of fgf10 and ANF. In addition we show that missexpression of moa tbx5 in the hindlimb of chicken embryos results in the formation of forelimb features, suggesting that Tbx5 was fully functional in wingless moa. An alternatively spliced exon 1 for tbx5 that is expressed specifically in the forelimb region was shown to be almost identical between moa and ostrich, suggesting that, as well as being fully functional, tbx5 is likely to have been expressed normally in moa since divergence from their flighted ancestors, approximately 60 mya. Conclusions The results suggests that, as in mice, moa tbx5 is necessary for the induction of forelimbs, but is not sufficient for their outgrowth. Moa Tbx5 may have played an important role in the development of moa's remnant forelimb girdle, and may be required for the formation of this structure. Our results further show that genetic changes affecting genes other than tbx5 must be responsible for the complete loss of forelimbs in moa.},
  language  = {en}
}
@phdthesis{Huang2021,
  author    = {Huang, Sichao},
  title     = {Past and present biodiversity in northeastern Siberia inferred from sedimentary DNA metabarcoding},
  school      = {Universit{\"a}t Potsdam},
  year      = {2021},
  abstract  = {The arctic-boreal treeline is a transition zone from taiga to tundra covering a vast area in Siberia. It often features large environmental gradients and reacts sensitively to changes in the environment. For example, the expansion of shrubs and a northward movement of the treeline are observable in Siberia as a response to the warming climate. The changes in vegetation across the treeline are known to influence the water chemistry in the lakes. This causes further alteration to the composition and diversity of sensitive aquatic organisms such as diatoms and macrophytes. Despite the rising awareness of the complex climate-feedback mechanisms of terrestrial plants, the understanding of their assembly rules and about responses of aquatic biomes in the surrounding treeline lakes is still limited. The goal of this thesis is to examine the previous and present biodiversity of terrestrial and freshwater biomes from the Siberian treeline ecotone, as well as their reactions to environmental changes. In particular, this thesis attempts to examine the performance of applying sedimentary DNA metabarcoding in terrestrial plants, aquatic macrophytes and diatoms, their spatial patterns along the environmental gradients and their temporal patterns throughout the climate transition from the late Pleistocene to Holocene. Sedimentary DNA metabarcoding combined with next-generation sequencing is applied as a primary tool to explore the composition and diversity of terrestrial plants, diatoms and aquatic macrophytes. The main study area is located in Chukotka of northeastern Siberia in the Arctic, a biodiversity hotspot due to its continental location and the diverse habitats of the glacial refugium. The modern diatom diversity was assessed with a specific diatom metabarcoding marker and morphological identification. Both approaches agree to a dominance of Fragilariaceae and Aulacoseiraceae, as well as on the environmental influential indicators of the diatom community. The high diversity of Fragilariaceae identified in the thermokarst lakes is found to follow the vegetation gradient along the treeline, suggesting that diatom metabarcoding can decipher relationships between diatom assemblage shifts and the relevant environmental changes. In particular, the metabarcoding approach detects diversification of fragilarioids in glacial lakes which is not visible using morphology. Sedimentary ancient DNA records indicate a vegetation mosaic of forb-dominated steppe-tundra during 28-19 ka, followed by a shift to dwarf-shrub tundra during 19-14 ka. During the most recent 14 thousand years, the vegetation consists of deciduous shrublands, then a change to boreal forest is observed. Investigations on the alpha diversity of the vegetation show that species richness is unexpectedly highest during pre-LGM, which is likely related to the extensive area that allows for more taxa. The optimum Holocene warming during 9-6 ka is not accompanied by a high richness as widely believed, but with an evenly distributed community by the fulfilment of erect shrubs. Furthermore, changes in taxonomic and phylogenetic diversity show complementary results in understanding community diversity. The composition and richness in the modern macrophytes community from Siberian Arctic and Chinese alpine are best co-influenced by July temperature and electrical conductivity.. Past macrophyte turnover during the late Pleistocene-Holocene is less noticeable in Siberia, whereas a pronounced community change from emergent to submerged plants is detected from Chinese alpine regions at about 14 ka due to increasing temperature and varying water conductivity. Finally, sedimentary DNA metabarcoding is a cost-effective and powerful proxy for ecological application, whereas completeness of the reference library, coverage and resolution of the metabarcoding marker are the major limitations of sedimentary DNA based diversity monitoring. The composition and richness in modern vegetation and macrophytes across broad spatial gradients is constrained by environmental variables, suggesting a potential usage for environmental monitoring. Diatom distributions are driven by different water variables along the treeline. Past records indicate that the shrub coverage has a noticeable influence on the assemblies of both terrestrial plants and aquatic macrophytes, though the shift in macrophyte community is relatively minor in the past 28 thousand years. In the long-term, the shrub expansion may eventually result in a genetically more diverse vegetation community but reduced species richness. When exceeding the optimal temperatures, further warming may lead to a decrease and putative loss of macrophytes and diatoms.},
  language  = {en}
}
@article{HofreiterPaijmansGoodchildetal.2015,
  author    = {Hofreiter, Michael and Paijmans, Johanna L. A. and Goodchild, Helen and Speller, Camilla F. and Barlow, Axel and Gonz{\´a}lez-Fortes, Gloria M. and Thomas, Jessica A. and Ludwig, Arne and Collins, Matthew J.},
  title     = {The future of ancient DNA: Technical advances and conceptual shifts},
  series = {Bioessays : ideas that push the boundaries},
  volume    = {37},
  journal   = {Bioessays : ideas that push the boundaries},
  number    = {3},
  publisher = {Wiley-Blackwell},
  address   = {Hoboken},
  issn      = {0265-9247},
  doi       = {10.1002/bies.201400160},
  pages     = {284 -- 293},
  year      = {2015},
  abstract  = {Technological innovations such as next generation sequencing and DNA hybridisation enrichment have resulted in multi-fold increases in both the quantity of ancient DNA sequence data and the time depth for DNA retrieval. To date, over 30 ancient genomes have been sequenced, moving from 0.7x coverage (mammoth) in 2008 to more than 50x coverage (Neanderthal) in 2014. Studies of rapid evolutionary changes, such as the evolution and spread of pathogens and the genetic responses of hosts, or the genetics of domestication and climatic adaptation, are developing swiftly and the importance of palaeogenomics for investigating evolutionary processes during the last million years is likely to increase considerably. However, these new datasets require new methods of data processing and analysis, as well as conceptual changes in interpreting the results. In this review we highlight important areas of future technical and conceptual progress and discuss research topics in the rapidly growing field of palaeogenomics.},
  language  = {en}
}
@misc{HofreiterPaijmansGoodchildetal.2015,
  author    = {Hofreiter, Michael and Paijmans, Johanna L. A. and Goodchild, Helen and Speller, Camilla F. and Barlow, Axel and Gonzalez-Fortes, Gloria M. and Thomas, Jessica A. and Ludwig, Arne and Collins, Matthew J.},
  title     = {The future of ancient DNA},
  series = {Postprints der Universit{\"a}t Potsdam : Mathematisch Naturwissenschaftliche Reihe},
  journal   = {Postprints der Universit{\"a}t Potsdam : Mathematisch Naturwissenschaftliche Reihe},
  number    = {908},
  issn      = {1866-8372},
  doi       = {10.25932/publishup-43881},
  url       = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-438816},
  pages     = {284 -- 295},
  year      = {2015},
  abstract  = {Technological innovations such as next generation sequencing and DNA hybridisation enrichment have resulted in multi-fold increases in both the quantity of ancient DNA sequence data and the time depth for DNA retrieval. To date, over 30 ancient genomes have been sequenced, moving from 0.7x coverage (mammoth) in 2008 to more than 50x coverage (Neanderthal) in 2014. Studies of rapid evolutionary changes, such as the evolution and spread of pathogens and the genetic responses of hosts, or the genetics of domestication and climatic adaptation, are developing swiftly and the importance of palaeogenomics for investigating evolutionary processes during the last million years is likely to increase considerably. However, these new datasets require new methods of data processing and analysis, as well as conceptual changes in interpreting the results. In this review we highlight important areas of future technical and conceptual progress and discuss research topics in the rapidly growing field of palaeogenomics.},
  language  = {en}
}
@article{HofreiterHartmann2020,
  author    = {Hofreiter, Michael and Hartmann, Stefanie},
  title     = {Reconstructing protein-coding sequences from ancient DNA},
  series = {Odorant binding and chemosensory proteins},
  volume    = {642},
  journal   = {Odorant binding and chemosensory proteins},
  publisher = {Academic Press, an imprint of Elsevier},
  address   = {Cambridge, MA.},
  isbn      = {978-0-12-821157-1},
  issn      = {0076-6879},
  doi       = {10.1016/bs.mie.2020.05.008},
  pages     = {21 -- 33},
  year      = {2020},
  abstract  = {Obtaining information about functional details of proteins of extinct species is of critical importance for a better understanding of the real-life appearance, behavior and ecology of these lost entries in the book of life. In this chapter, we discuss the possibilities to retrieve the necessary DNA sequence information from paleogenomic data obtained from fossil specimens, which can then be used to express and subsequently analyze the protein of interest. We discuss the problems specific to ancient DNA, including mis-coding lesions, short read length and incomplete paleogenome assemblies. Finally, we discuss an alternative, but currently rarely used approach, direct PCR amplification, which is especially useful for comparatively short proteins.},
  language  = {en}
}
@unpublished{HagelbergHofreiterKeyser2015,
  author    = {Hagelberg, Erika and Hofreiter, Michael and Keyser, Christine},
  title     = {Ancient DNA: the first three decades},
  series = {Philosophical transactions of the Royal Society of London : B, Biological sciences},
  volume    = {370},
  journal   = {Philosophical transactions of the Royal Society of London : B, Biological sciences},
  number    = {1660},
  publisher = {Royal Society},
  address   = {London},
  issn      = {0962-8436},
  doi       = {10.1098/rstb.2013.0371},
  pages     = {6},
  year      = {2015},
  language  = {en}
}
@article{GonzalezFortesKolbeFernandesetal.2016,
  author    = {Gonz{\´a}lez-Fortes, Gloria M. and Kolbe, Ben and Fernandes, Daniel and Meleg, Ioana N. and Garcia-Vazquez, Ana and Pinto-Llona, Ana C. and Constantin, Silviu and de Torres, Trino J. and Ortiz, Jose E. and Frischauf, Christine and Rabeder, Gernot and Hofreiter, Michael and Barlow, Axel},
  title     = {Ancient DNA reveals differences in behaviour and sociality between brown bears and extinct cave bears},
  series = {Molecular ecology},
  volume    = {25},
  journal   = {Molecular ecology},
  publisher = {Wiley-Blackwell},
  address   = {Hoboken},
  issn      = {0962-1083},
  doi       = {10.1111/mec.13800},
  pages     = {4907 -- 4918},
  year      = {2016},
  abstract  = {Ancient DNA studies have revolutionized the study of extinct species and populations, providing insights on phylogeny, phylogeography, admixture and demographic history. However, inferences on behaviour and sociality have been far less frequent. Here, we investigate the complete mitochondrial genomes of extinct Late Pleistocene cave bears and middle Holocene brown bears that each inhabited multiple geographically proximate caves in northern Spain. In cave bears, we find that, although most caves were occupied simultaneously, each cave almost exclusively contains a unique lineage of closely related haplotypes. This remarkable pattern suggests extreme fidelity to their birth site in cave bears, best described as homing behaviour, and that cave bears formed stable maternal social groups at least for hibernation. In contrast, brown bears do not show any strong association of mitochondrial lineage and cave, suggesting that these two closely related species differed in aspects of their behaviour and sociality. This difference is likely to have contributed to cave bear extinction, which occurred at a time in which competition for caves between bears and humans was likely intense and the ability to rapidly colonize new hibernation sites would have been crucial for the survival of a species so dependent on caves for hibernation as cave bears. Our study demonstrates the potential of ancient DNA to uncover patterns of behaviour and sociality in ancient species and populations, even those that went extinct many tens of thousands of years ago.},
  language  = {en}
}
@misc{GonzalezFortesJonesLightfootetal.2017,
  author    = {Gonz{\´a}lez-Fortes, Gloria M. and Jones, Eppie R. and Lightfoot, Emma and Bonsall, Clive and Lazar, Catalin and Grandal-d'Anglade, Aurora and Garralda, Mar{\´i}a Dolores and Drak, Labib and Siska, Veronika and Simalcsik, Angela and Boroneant, Adina and Roman{\´i}, Juan Ram{\´o}n Vidal and Vaqueiro Rodr{\´i}guez, Marcos and Arias, Pablo and Pinhasi, Ron and Manica, Andrea and Hofreiter, Michael},
  title     = {Paleogenomic evidence for multi-generational mixing between Neolithic Farmers and Mesolithic Hunter-Gatherers in the lower Danube Basin},
  series = {Postprints der Universit{\"a}t Potsdam Mathematisch-Naturwissenschaftliche Reihe},
  journal   = {Postprints der Universit{\"a}t Potsdam Mathematisch-Naturwissenschaftliche Reihe},
  number    = {794},
  issn      = {1866-8372},
  doi       = {10.25932/publishup-44011},
  url       = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-440115},
  pages     = {1801 -- 1820},
  year      = {2017},
  abstract  = {The transition from hunting and gathering to farming involved profound cultural and technological changes. In Western and Central Europe, these changes occurred rapidly and synchronously after the arrival of early farmers of Anatolian origin [1-3], who largely replaced the local Mesolithic hunter-gatherers [1, 4-6]. Further east, in the Baltic region, the transition was gradual, with little or no genetic input from incoming farmers [7]. Here we use ancient DNA to investigate the relationship between hunter-gatherers and farmers in the Lower Danube basin, a geographically intermediate area that is characterized by a rapid Neolithic transition but also by the presence of archaeological evidence that points to cultural exchange, and thus possible admixture, between hunter-gatherers and farmers. We recovered four human paleogenomes (1.13 to 4.13 coverage) from Romania spanning a time transect between 8.8 thousand years ago (kya) and 5.4 kya and supplemented them with two Mesolithic genomes (1.73- and 5.33) from Spain to provide further context on the genetic background of Mesolithic Europe. Our results show major Western hunter-gatherer (WHG) ancestry in a Romanian Eneolithic sample with a minor, but sizeable, contribution from Anatolian farmers, suggesting multiple admixture events between hunter-gatherers and farmers. Dietary stableisotope analysis of this sample suggests a mixed terrestrial/ aquatic diet. Our results provide support for complex interactions among hunter-gatherers and farmers in the Danube basin, demonstrating that in some regions, demic and cultural diffusion were not mutually exclusive, but merely the ends of a continuum for the process of Neolithization.},
  language  = {en}
}
@article{GluecklerGengGrimmetal.2022,
  author    = {Gl{\"u}ckler, Ramesh and Geng, Rongwei and Grimm, Lennart and Baisheva, Izabella and Herzschuh, Ulrike and Stoof-Leichsenring, Kathleen R. and Kruse, Stefan and Andreev, Andrej Aleksandrovic and Pestryakova, Luidmila and Dietze, Elisabeth},
  title     = {Holocene wildfire and vegetation dynamics in Central Yakutia, Siberia, reconstructed from lake-sediment proxies},
  series = {Frontiers in Ecology and Evolution},
  volume    = {10},
  journal   = {Frontiers in Ecology and Evolution},
  publisher = {Frontiers Media},
  address   = {Lausanne},
  issn      = {2296-701X},
  doi       = {10.3389/fevo.2022.962906},
  pages     = {19},
  year      = {2022},
  abstract  = {Wildfires play an essential role in the ecology of boreal forests. In eastern Siberia, fire activity has been increasing in recent years, challenging the livelihoods of local communities. Intensifying fire regimes also increase disturbance pressure on the boreal forests, which currently protect the permafrost beneath from accelerated degradation. However, long-term relationships between changes in fire regime and forest structure remain largely unknown. We assess past fire-vegetation feedbacks using sedimentary proxy records from Lake Satagay, Central Yakutia, Siberia, covering the past c. 10,800 years. Results from macroscopic and microscopic charcoal analyses indicate high amounts of burnt biomass during the Early Holocene, and that the present-day, low-severity surface fire regime has been in place since c. 4,500 years before present. A pollen-based quantitative reconstruction of vegetation cover and a terrestrial plant record based on sedimentary ancient DNA metabarcoding suggest a pronounced shift in forest structure toward the Late Holocene. Whereas the Early Holocene was characterized by postglacial open larch-birch woodlands, forest structure changed toward the modern, mixed larch-dominated closed-canopy forest during the Mid-Holocene. We propose a potential relationship between open woodlands and high amounts of burnt biomass, as well as a mediating effect of dense larch forest on the climate-driven intensification of fire regimes. Considering the anticipated increase in forest disturbances (droughts, insect invasions, and wildfires), higher tree mortality may force the modern state of the forest to shift toward an open woodland state comparable to the Early Holocene. Such a shift in forest structure may result in a positive feedback on currently intensifying wildfires. These new long-term data improve our understanding of millennial-scale fire regime changes and their relationships to changes of vegetation in Central Yakutia, where the local population is already being confronted with intensifying wildfire seasons.},
  language  = {en}
}
@misc{GambaJonesTeasdaleetal.2014,
  author    = {Gamba, Cristina and Jones, Eppie R. and Teasdale, Matthew D. and McLaughlin, Russell L. and Gonz{\´a}lez-Fortes, Gloria M. and Mattiangeli, Valeria and Dombor{\´o}czki, L{\´a}szl{\´o} and Kőv{\´a}ri, Ivett and Pap, Ildik{\´o} and Anders, Alexandra and Whittle, Alasdair and Dani, J{\´a}nos and Raczky, P{\´a}l and Higham, Thomas F. G. and Hofreiter, Michael and Bradley, Daniel G. and Pinhasi, Ron},
  title     = {Genome flux and stasis in a five millennium transect of European prehistory},
  series = {Zweitver{\"o}ffentlichungen der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe},
  volume    = {5},
  journal   = {Zweitver{\"o}ffentlichungen der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe},
  number    = {1332},
  issn      = {1866-8372},
  doi       = {10.25932/publishup-43799},
  url       = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-437999},
  pages     = {9},
  year      = {2014},
  abstract  = {The Great Hungarian Plain was a crossroads of cultural transformations that have shaped European prehistory. Here we analyse a 5,000-year transect of human genomes, sampled from petrous bones giving consistently excellent endogenous DNA yields, from 13 Hungarian Neolithic, Copper, Bronze and Iron Age burials including two to high (similar to 22x) and seven to similar to 1x coverage, to investigate the impact of these on Europe's genetic landscape. These data suggest genomic shifts with the advent of the Neolithic, Bronze and Iron Ages, with interleaved periods of genome stability. The earliest Neolithic context genome shows a European hunter-gatherer genetic signature and a restricted ancestral population size, suggesting direct contact between cultures after the arrival of the first farmers into Europe. The latest, Iron Age, sample reveals an eastern genomic influence concordant with introduced Steppe burial rites. We observe transition towards lighter pigmentation and surprisingly, no Neolithic presence of lactase persistence.},
  language  = {en}
}
@misc{ElsnerSchiblerHofreiteretal.2015,
  author    = {Elsner, Julia and Schibler, J{\"o}rg and Hofreiter, Michael and Schlumbaum, Angela},
  title     = {Burial condition is the most important factor for mtDNA PCR amplification success in Palaeolithic equid remains from the Alpine foreland},
  series = {Postprints der Universit{\"a}t Potsdam Mathematisch-Naturwissenschaftliche Reihe},
  journal   = {Postprints der Universit{\"a}t Potsdam Mathematisch-Naturwissenschaftliche Reihe},
  number    = {727},
  issn      = {1866-8372},
  doi       = {10.25932/publishup-42976},
  url       = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-429763},
  pages     = {505 -- 515},
  year      = {2015},
  abstract  = {Faunal remains from Palaeolithic sites are important genetic sources to study preglacial and postglacial populations and to investigate the effect of climate change and human impact. Post mortem decay, resulting in fragmented and chemically modified DNA, is a key obstacle in ancient DNA analyses. In the absence of reliable methods to determine the presence of endogenous DNA in sub-fossil samples, temporal and spatial surveys of DNA survival on a regional scale may help to estimate the potential of faunal remains from a given time period and region. We therefore investigated PCR amplification success, PCR performance and post mortem damage in c. 47,000 to c. 12,000-year-old horse remains from 14 Palaeolithic sites along the Swiss Jura Mountains in relation to depositional context, tissue type, storage time and age, potentially influencing DNA preservation. The targeted 75 base pair mitochondrial DNA fragment could be amplified solely from equid remains from caves and not from any of the open dry and (temporary) wetland sites. Whether teeth are better than bones cannot be ultimately decided; however, both storage time after excavation and age significantly affect PCR amplification and performance, albeit not in a linear way. This is best explained by the—inevitable—heterogeneity of the data set. The extent of post mortem damage is not related to any of the potential impact factors. The results encourage comprehensive investigations of Palaeolithic cave sites, even from temperate regions.},
  language  = {en}
}
@article{ChenLiZhangetal.2019,
  author    = {Chen, Shun-Gang and Li, Ji and Zhang, Fan and Xiao, Bo and Hu, Jia-Ming and Cui, Yin-Qiu and Hofreiter, Michael and Hou, Xin-Dong and Sheng, Gui-Lian and Lai, Xu-Long and Yuan, Jun-Xia},
  title     = {Different maternal lineages revealed by ancient mitochondrial genome of Camelus bactrianus from China},
  series = {Mitochondrial DNA Part A},
  volume    = {30},
  journal   = {Mitochondrial DNA Part A},
  number    = {7},
  publisher = {Routledge, Taylor \& Francis Group},
  address   = {Abingdon},
  issn      = {2470-1394},
  doi       = {10.1080/24701394.2019.1659250},
  pages     = {786 -- 793},
  year      = {2019},
  abstract  = {Domestic Bactrian camel (Camelus bactrianus) used to be one of the most important livestock species in Chinese history, as well as the major transport carrier on the ancient Silk Road. However, archeological studies on Chinese C. bactrianus are still limited, and molecular biology research on this species is mainly focused on modern specimens. In this study, we retrieved the complete mitochondrial genome from a C. bactrianus specimen, which was excavated from northwestern China and dated at 1290-1180 cal. Phylogenetic analyses using 18 mitochondrial genomes indicated that the C. bactrianus clade was divided into two maternal lineages. The majority of samples originating from Iran to Japan and Mongolia belong to subclade A1, while our sample together with two Mongolian individuals formed the much smaller subclade A2. Furthermore, the divergence time of these two maternal lineages was estimated as 165 Kya (95\% credibility interval 117-222 Kya), this might indicate that several different evolutionary lineages were incorporated into the domestic gene pool during the initial domestication process. Bayesian skyline plot (BSP) analysis a slow increase in female effective population size of C. bactrianus from 5000 years ago, which to the beginning of domestication of C. bactrianus. The present study also revealed that there were extensive exchanges of genetic information among C. bactrianus populations in regions along the Silk Road.},
  language  = {en}
}
@article{CasasMarceMarmesatSorianoetal.2017,
  author    = {Casas-Marce, Mireia and Marmesat, Elena and Soriano, Laura and Martinez-Cruz, Begona and Lucena-Perez, Maria and Nocete, Francisco and Rodriguez-Hidalgo, Antonio and Canals, Antoni and Nadal, Jordi and Detry, Cleia and Bernaldez-Sanchez, Eloisa and Fernandez-Rodriguez, Carlos and Perez-Ripoll, Manuel and Stiller, Mathias and Hofreiter, Michael and Rodriguez, Alejandro and Revilla, Eloy and Delibes, Miguel and Godoy, Jose A.},
  title     = {Spatiotemporal Dynamics of Genetic Variation in the Iberian Lynx along Its Path to Extinction Reconstructed with Ancient DNA},
  series = {Molecular biology and evolution},
  volume    = {34},
  journal   = {Molecular biology and evolution},
  publisher = {Oxford Univ. Press},
  address   = {Oxford},
  issn      = {0737-4038},
  doi       = {10.1093/molbev/msx222},
  pages     = {2893 -- 2907},
  year      = {2017},
  abstract  = {There is the tendency to assume that endangered species have been both genetically and demographically healthier in the past, so that any genetic erosion observed today was caused by their recent decline. The Iberian lynx (Lynx pardinus) suffered a dramatic and continuous decline during the 20th century, and now shows extremely low genome- and species-wide genetic diversity among other signs of genomic erosion. We analyze ancient (N\&\#8201;=\&\#8201;10), historical (N\&\#8201;=\&\#8201;245), and contemporary (N\&\#8201;=\&\#8201;172) samples with microsatellite and mitogenome data to reconstruct the species' demography and investigate patterns of genetic variation across space and time. Iberian lynx populations transitioned from low but significantly higher genetic diversity than today and shallow geographical differentiation millennia ago, through a structured metapopulation with varying levels of diversity during the last centuries, to two extremely genetically depauperate and differentiated remnant populations by 2002. The historical subpopulations show varying extents of genetic drift in relation to their recent size and time in isolation, but these do not predict whether the populations persisted or went finally extinct. In conclusion, current genetic patterns were mainly shaped by genetic drift, supporting the current admixture of the two genetic pools and calling for a comprehensive genetic management of the ongoing conservation program. This study illustrates how a retrospective analysis of demographic and genetic patterns of endangered species can shed light onto their evolutionary history and this, in turn, can inform conservation actions.},
  language  = {en}
}
@misc{CampbellHofreiter2015,
  author    = {Campbell, Kevin L. and Hofreiter, Michael},
  title     = {Resurrecting phenotypes from ancient DNA sequences: promises and perspectives},
  series = {Canadian journal of zoology = Revue canadienne de zoologie},
  volume    = {93},
  journal   = {Canadian journal of zoology = Revue canadienne de zoologie},
  number    = {9},
  publisher = {NRC Research Press},
  address   = {Ottawa},
  issn      = {0008-4301},
  doi       = {10.1139/cjz-2014-0337},
  pages     = {701 -- 710},
  year      = {2015},
  abstract  = {Anatomical changes in extinct mammalian lineages over evolutionary time, such as the loss of fingers and teeth and the rapid increase in body size that accompanied the late Miocene dispersal of the progenitors of Steller's sea cows (Hydrodamalis gigas (Zimmermann, 1780)) into North Pacific waters and the convergent development of a thick pelage and accompanying reductions in ear and tail surface area of woolly mammoths (Mammuthus primigenius (Blumenbach, 1799)) and woolly rhinoceros (Coelodonta antiquitatis (Blumenbach, 1799)), are prime examples of adaptive evolution underlying the exploitation of new habitats. It is likely, however, that biochemical specializations adopted during these evolutionary transitions were of similar or even greater biological importance. As these "living" processes do not fossilize, direct information regarding the physiological attributes of extinct species has largely remained beyond the range of scientific inquiry. However, the ability to retrieve genomic sequences from ancient DNA samples, combined with ectopic expression systems, now permit the evolutionary origins and structural and functional properties of authentic prehistoric proteins to be examined in great detail. Exponential technical advances in ancient DNA retrieval, enrichment, and sequencing will soon permit targeted generation of complete genomes from hundreds of extinct species across the last one million years that, in combination with emerging in vitro expression, genome engineering, and cell differentiation techniques, promises to herald an exciting new trajectory of evolutionary research at the interface of biochemistry, genomics, palaeontology, and cell biology.},
  language  = {en}
}
@misc{BarlowHartmannGonzalezetal.2020,
  author    = {Barlow, Axel and Hartmann, Stefanie and Gonzalez, Javier and Hofreiter, Michael and Paijmans, Johanna L. A.},
  title     = {Consensify},
  series = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe},
  journal   = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe},
  number    = {1033},
  issn      = {1866-8372},
  doi       = {10.25932/publishup-47252},
  url       = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-472521},
  pages     = {24},
  year      = {2020},
  abstract  = {A standard practise in palaeogenome analysis is the conversion of mapped short read data into pseudohaploid sequences, frequently by selecting a single high-quality nucleotide at random from the stack of mapped reads. This controls for biases due to differential sequencing coverage, but it does not control for differential rates and types of sequencing error, which are frequently large and variable in datasets obtained from ancient samples. These errors have the potential to distort phylogenetic and population clustering analyses, and to mislead tests of admixture using D statistics. We introduce Consensify, a method for generating pseudohaploid sequences, which controls for biases resulting from differential sequencing coverage while greatly reducing error rates. The error correction is derived directly from the data itself, without the requirement for additional genomic resources or simplifying assumptions such as contemporaneous sampling. For phylogenetic and population clustering analysis, we find that Consensify is less affected by artefacts than methods based on single read sampling. For D statistics, Consensify is more resistant to false positives and appears to be less affected by biases resulting from different laboratory protocols than other frequently used methods. Although Consensify is developed with palaeogenomic data in mind, it is applicable for any low to medium coverage short read datasets. We predict that Consensify will be a useful tool for future studies of palaeogenomes.},
  language  = {en}
}
@article{BarlowHartmannGonzalezetal.2020,
  author    = {Barlow, Axel and Hartmann, Stefanie and Gonzalez, Javier and Hofreiter, Michael and Paijmans, Johanna L. A.},
  title     = {Consensify},
  series = {Genes / Molecular Diversity Preservation International},
  volume    = {11},
  journal   = {Genes / Molecular Diversity Preservation International},
  number    = {1},
  publisher = {MDPI},
  address   = {Basel},
  issn      = {2073-4425},
  doi       = {10.3390/genes11010050},
  pages     = {22},
  year      = {2020},
  abstract  = {A standard practise in palaeogenome analysis is the conversion of mapped short read data into pseudohaploid sequences, frequently by selecting a single high-quality nucleotide at random from the stack of mapped reads. This controls for biases due to differential sequencing coverage, but it does not control for differential rates and types of sequencing error, which are frequently large and variable in datasets obtained from ancient samples. These errors have the potential to distort phylogenetic and population clustering analyses, and to mislead tests of admixture using D statistics. We introduce Consensify, a method for generating pseudohaploid sequences, which controls for biases resulting from differential sequencing coverage while greatly reducing error rates. The error correction is derived directly from the data itself, without the requirement for additional genomic resources or simplifying assumptions such as contemporaneous sampling. For phylogenetic and population clustering analysis, we find that Consensify is less affected by artefacts than methods based on single read sampling. For D statistics, Consensify is more resistant to false positives and appears to be less affected by biases resulting from different laboratory protocols than other frequently used methods. Although Consensify is developed with palaeogenomic data in mind, it is applicable for any low to medium coverage short read datasets. We predict that Consensify will be a useful tool for future studies of palaeogenomes.},
  language  = {en}
}
@misc{AlterMeyerPostetal.2015,
  author    = {Alter, S. Elizabeth and Meyer, Matthias and Post, Klaas and Czechowski, Paul and Gravlund, Peter and Gaines, Cork and Rosenbaum, Howard C. and Kaschner, Kristin and Turvey, Samuel T. and van der Plicht, Johannes and Shapiro, Beth and Hofreiter, Michael},
  title     = {Climate impacts on transocean dispersal and habitat in gray whales from the Pleistocene to 2100},
  series = {Postprints der Universit{\"a}t Potsdam : Mathematisch Naturwissenschaftliche Reihe},
  journal   = {Postprints der Universit{\"a}t Potsdam : Mathematisch Naturwissenschaftliche Reihe},
  number    = {965},
  issn      = {1866-8372},
  doi       = {10.25932/publishup-43892},
  url       = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-438920},
  pages     = {1510 -- 1522},
  year      = {2015},
  abstract  = {Arctic animals face dramatic habitat alteration due to ongoing climate change. Understanding how such species have responded to past glacial cycles can help us forecast their response to today's changing climate. Gray whales are among those marine species likely to be strongly affected by Arctic climate change, but a thorough analysis of past climate impacts on this species has been complicated by lack of information about an extinct population in the Atlantic. While little is known about the history of Atlantic gray whales or their relationship to the extant Pacific population, the extirpation of the Atlantic population during historical times has been attributed to whaling. We used a combination of ancient and modern DNA, radiocarbon dating and predictive habitat modelling to better understand the distribution of gray whales during the Pleistocene and Holocene. Our results reveal that dispersal between the Pacific and Atlantic was climate dependent and occurred both during the Pleistocene prior to the last glacial period and the early Holocene immediately following the opening of the Bering Strait. Genetic diversity in the Atlantic declined over an extended interval that predates the period of intensive commercial whaling, indicating this decline may have been precipitated by Holocene climate or other ecological causes. These first genetic data for Atlantic gray whales, particularly when combined with predictive habitat models for the year 2100, suggest that two recent sightings of gray whales in the Atlantic may represent the beginning of the expansion of this species' habitat beyond its currently realized range.},
  language  = {en}
}
@article{AlterMeyerPostetal.2015,
  author    = {Alter, S. Elizabeth and Meyer, Matthias and Post, Klaas and Czechowski, Paul and Gravlund, Peter and Gaines, Cork and Rosenbaum, Howard C. and Kaschner, Kristin and Turvey, Samuel T. and van der Plicht, Johannes and Shapiro, Beth and Hofreiter, Michael},
  title     = {Climate impacts on transocean dispersal and habitat in gray whales from the Pleistocene to 2100},
  series = {Molecular ecology},
  volume    = {24},
  journal   = {Molecular ecology},
  number    = {7},
  publisher = {Wiley-Blackwell},
  address   = {Hoboken},
  issn      = {0962-1083},
  doi       = {10.1111/mec.13121},
  pages     = {1510 -- 1522},
  year      = {2015},
  abstract  = {Arctic animals face dramatic habitat alteration due to ongoing climate change. Understanding how such species have responded to past glacial cycles can help us forecast their response to today's changing climate. Gray whales are among those marine species likely to be strongly affected by Arctic climate change, but a thorough analysis of past climate impacts on this species has been complicated by lack of information about an extinct population in the Atlantic. While little is known about the history of Atlantic gray whales or their relationship to the extant Pacific population, the extirpation of the Atlantic population during historical times has been attributed to whaling. We used a combination of ancient and modern DNA, radiocarbon dating and predictive habitat modelling to better understand the distribution of gray whales during the Pleistocene and Holocene. Our results reveal that dispersal between the Pacific and Atlantic was climate dependent and occurred both during the Pleistocene prior to the last glacial period and the early Holocene immediately following the opening of the Bering Strait. Genetic diversity in the Atlantic declined over an extended interval that predates the period of intensive commercial whaling, indicating this decline may have been precipitated by Holocene climate or other ecological causes. These first genetic data for Atlantic gray whales, particularly when combined with predictive habitat models for the year 2100, suggest that two recent sightings of gray whales in the Atlantic may represent the beginning of the expansion of this species' habitat beyond its currently realized range.},
  language  = {en}
}
@article{AlbertiGonzalezPaijmansetal.2018,
  author    = {Alberti, Federica and Gonzalez, Javier and Paijmans, Johanna L. A. and Basler, Nikolas and Preick, Michaela and Henneberger, Kirstin and Trinks, Alexandra and Rabeder, Gernot and Conard, Nicholas J. and Muenzel, Susanne C. and Joger, Ulrich and Fritsch, Guido and Hildebrandt, Thomas and Hofreiter, Michael and Barlow, Axel},
  title     = {Optimized DNA sampling of ancient bones using Computed Tomography scans},
  series = {Molecular ecology resources},
  volume    = {18},
  journal   = {Molecular ecology resources},
  number    = {6},
  publisher = {Wiley},
  address   = {Hoboken},
  issn      = {1755-098X},
  doi       = {10.1111/1755-0998.12911},
  pages     = {1196 -- 1208},
  year      = {2018},
  abstract  = {The prevalence of contaminant microbial DNA in ancient bone samples represents the principal limiting factor for palaeogenomic studies, as it may comprise more than 99\% of DNA molecules obtained. Efforts to exclude or reduce this contaminant fraction have been numerous but also variable in their success. Here, we present a simple but highly effective method to increase the relative proportion of endogenous molecules obtained from ancient bones. Using computed tomography (CT) scanning, we identify the densest region of a bone as optimal for sampling. This approach accurately identifies the densest internal regions of petrous bones, which are known to be a source of high-purity ancient DNA. For ancient long bones, CT scans reveal a high-density outermost layer, which has been routinely removed and discarded prior to DNA extraction. For almost all long bones investigated, we find that targeted sampling of this outermost layer provides an increase in endogenous DNA content over that obtained from softer, trabecular bone. This targeted sampling can produce as much as 50-fold increase in the proportion of endogenous DNA, providing a directly proportional reduction in sequencing costs for shotgun sequencing experiments. The observed increases in endogenous DNA proportion are not associated with any reduction in absolute endogenous molecule recovery. Although sampling the outermost layer can result in higher levels of human contamination, some bones were found to have more contamination associated with the internal bone structures. Our method is highly consistent, reproducible and applicable across a wide range of bone types, ages and species. We predict that this discovery will greatly extend the potential to study ancient populations and species in the genomics era.},
  language  = {en}
}