TY - JOUR A1 - Yuan, Junxia A1 - Sheng, Guilian A1 - Preick, Michaela A1 - Sun, Boyang A1 - Hou, Xindong A1 - Chen, Shungang A1 - Taron, Ulrike Helene A1 - Barlow, Axel A1 - Wang, Linying A1 - Hu, Jiaming A1 - Deng, Tao A1 - Lai, Xulong A1 - Hofreiter, Michael T1 - Mitochondrial genomes of Late Pleistocene caballine horses from China belong to a separate clade JF - Quaternary science reviews : the international multidisciplinary research and review journal N2 - 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. KW - Equus dalianensis KW - Equus przewalskii KW - Pleistocene caballine horses KW - ancient DNA KW - phylogenetic relationship KW - divergence time Y1 - 2020 U6 - https://doi.org/10.1016/j.quascirev.2020.106691 SN - 0277-3791 VL - 250 PB - Elsevier CY - Amsterdam [u.a.] ER - TY - GEN A1 - Taron, Ulrike H. A1 - Lell, Moritz A1 - Barlow, Axel A1 - Paijmans, Johanna L. A. T1 - Testing of Alignment Parameters for Ancient Samples BT - Evaluating and Optimizing Mapping Parameters for Ancient Samples Using the TAPAS Tool T2 - Genes N2 - 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. T3 - Zweitveröffentlichungen der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe - 415 KW - ancient DNA KW - short-read mapping KW - palaeogenomics KW - alignment sensitivity / specificity Y1 - 2018 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-409683 ER - TY - JOUR A1 - Taron, Ulrike H. A1 - Lell, Moritz A1 - Barlow, Axel A1 - Paijmans, Johanna L. A. T1 - Testing of Alignment Parameters for Ancient Samples BT - Evaluating and Optimizing Mapping Parameters for Ancient Samples Using the TAPAS Tool JF - Genes N2 - 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. KW - ancient DNA KW - short-read mapping KW - palaeogenomics KW - alignment sensitivity / specificity Y1 - 2018 U6 - https://doi.org/10.3390/genes9030157 SN - 2073-4425 VL - 9 IS - 3 SP - 1 EP - 12 PB - Molecular Diversity Preservation International CY - Basel ER - TY - JOUR A1 - Taron, Ulrike H. A1 - Lell, Moritz A1 - Barlow, Axel A1 - Paijmans, Johanna L. A. T1 - Testing of Alignment Parameters for Ancient Samples BT - Evaluating and Optimizing Mapping Parameters for Ancient Samples Using the TAPAS Tool JF - Genese N2 - 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. KW - ancient DNA KW - short-read mapping KW - palaeogenomics KW - paleogenomics KW - alignment sensitivity/specificity Y1 - 2018 U6 - https://doi.org/10.3390/genes9030157 SN - 2073-4425 VL - 9 IS - 3 PB - MDPI CY - Basel ER - TY - JOUR A1 - Hofreiter, Michael A1 - Paijmans, Johanna L. A. A1 - Goodchild, Helen A1 - Speller, Camilla F. A1 - Barlow, Axel A1 - González-Fortes, Gloria M. A1 - Thomas, Jessica A. A1 - Ludwig, Arne A1 - Collins, Matthew J. T1 - The future of ancient DNA: Technical advances and conceptual shifts JF - Bioessays : ideas that push the boundaries N2 - 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. KW - ancient DNA KW - hybridisation capture KW - multi-locus data KW - next generation sequencing (NGS) KW - palaeogenomics KW - population genomics Y1 - 2015 U6 - https://doi.org/10.1002/bies.201400160 SN - 0265-9247 SN - 1521-1878 VL - 37 IS - 3 SP - 284 EP - 293 PB - Wiley-Blackwell CY - Hoboken ER - TY - JOUR A1 - González-Fortes, Gloria M. A1 - Kolbe, Ben A1 - Fernandes, Daniel A1 - Meleg, Ioana N. A1 - Garcia-Vazquez, Ana A1 - Pinto-Llona, Ana C. A1 - Constantin, Silviu A1 - de Torres, Trino J. A1 - Ortiz, Jose E. A1 - Frischauf, Christine A1 - Rabeder, Gernot A1 - Hofreiter, Michael A1 - Barlow, Axel T1 - Ancient DNA reveals differences in behaviour and sociality between brown bears and extinct cave bears JF - Molecular ecology N2 - 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. KW - ancient DNA KW - extinction KW - homing KW - sociality KW - Ursus arctos KW - Ursus spelaeus Y1 - 2016 U6 - https://doi.org/10.1111/mec.13800 SN - 0962-1083 SN - 1365-294X VL - 25 SP - 4907 EP - 4918 PB - Wiley-Blackwell CY - Hoboken ER - TY - GEN A1 - Barlow, Axel A1 - Hartmann, Stefanie A1 - Gonzalez, Javier A1 - Hofreiter, Michael A1 - Paijmans, Johanna L. A. T1 - Consensify BT - a method for generating pseudohaploid genome sequences from palaeogenomic datasets with reduced error rates T2 - Postprints der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe N2 - 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. T3 - Zweitveröffentlichungen der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe - 1033 KW - palaeogenomics KW - ancient DNA KW - sequencing error KW - error reduction KW - D statistics KW - bioinformatics Y1 - 2020 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-472521 SN - 1866-8372 IS - 1033 ER - TY - JOUR A1 - Barlow, Axel A1 - Hartmann, Stefanie A1 - Gonzalez, Javier A1 - Hofreiter, Michael A1 - Paijmans, Johanna L. A. T1 - Consensify BT - a method for generating pseudohaploid genome sequences from palaeogenomic datasets with reduced error rates JF - Genes / Molecular Diversity Preservation International N2 - 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. KW - palaeogenomics KW - ancient DNA KW - sequencing error KW - error reduction KW - D statistics KW - bioinformatics Y1 - 2020 U6 - https://doi.org/10.3390/genes11010050 SN - 2073-4425 VL - 11 IS - 1 PB - MDPI CY - Basel ER - TY - JOUR A1 - Alberti, Federica A1 - Gonzalez, Javier A1 - Paijmans, Johanna L. A. A1 - Basler, Nikolas A1 - Preick, Michaela A1 - Henneberger, Kirstin A1 - Trinks, Alexandra A1 - Rabeder, Gernot A1 - Conard, Nicholas J. A1 - Muenzel, Susanne C. A1 - Joger, Ulrich A1 - Fritsch, Guido A1 - Hildebrandt, Thomas A1 - Hofreiter, Michael A1 - Barlow, Axel T1 - Optimized DNA sampling of ancient bones using Computed Tomography scans JF - Molecular ecology resources N2 - 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. KW - ancient DNA KW - computer tomography KW - palaeogenomics KW - paleogenetics KW - petrous bone Y1 - 2018 U6 - https://doi.org/10.1111/1755-0998.12911 SN - 1755-098X SN - 1755-0998 VL - 18 IS - 6 SP - 1196 EP - 1208 PB - Wiley CY - Hoboken ER -