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 - 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 - Hofreiter, Michael A1 - Hartmann, Stefanie T1 - Reconstructing protein-coding sequences from ancient DNA JF - Odorant binding and chemosensory proteins N2 - 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. KW - re-sequencing KW - mapping KW - genome assembly KW - targeted assembly KW - SRAssembler KW - ancient DNA KW - reference sequence KW - paleogenomics Y1 - 2020 SN - 978-0-12-821157-1 U6 - https://doi.org/10.1016/bs.mie.2020.05.008 SN - 0076-6879 VL - 642 SP - 21 EP - 33 PB - Academic Press, an imprint of Elsevier CY - Cambridge, MA. ER - TY - JOUR A1 - Mohandesan, Elmira A1 - Speller, Camilla F. A1 - Peters, Joris A1 - Uerpmann, Hans-Peter A1 - Uerpmann, Margarethe A1 - De Cupere, Bea A1 - Hofreiter, Michael A1 - Burger, Pamela A. T1 - Combined hybridization capture and shotgun sequencing for ancient DNA analysis of extinct wild and domestic dromedary camel JF - Molecular ecology resources N2 - 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. KW - ancient DNA KW - Camelus dromedarius KW - capture enrichment KW - degraded DNA KW - mitochondrial genome (mtDNA) KW - next-generation sequencing Y1 - 2017 U6 - https://doi.org/10.1111/1755-0998.12551 SN - 1755-098X SN - 1755-0998 VL - 17 IS - 2 SP - 300 EP - 313 PB - Wiley CY - Hoboken ER - TY - JOUR A1 - Thomas, Jessica E. A1 - Carvalho, Gary R. A1 - Haile, James A1 - Martin, Michael D. A1 - Castruita, Jose A. Samaniego A1 - Niemann, Jonas A1 - Sinding, Mikkel-Holger S. A1 - Sandoval-Velasco, Marcela A1 - Rawlence, Nicolas J. A1 - Fuller, Errol A1 - Fjeldsa, Jon A1 - Hofreiter, Michael A1 - Stewart, John R. A1 - Gilbert, M. Thomas P. A1 - Knapp, Michael T1 - An ‛Aukward’ tale BT - a genetic approach to discover the whereabouts of the Last Great Auks JF - Genes N2 - 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. KW - ancient DNA KW - extinct birds KW - mitochondrial genome KW - museum specimens KW - palaeogenomics Y1 - 2017 U6 - https://doi.org/10.3390/genes8060164 SN - 2073-4425 VL - 8 IS - 6 SP - 164 PB - MDPI 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 - 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 - TY - JOUR A1 - Signore, Anthony V. A1 - Paijmans, Johanna L. A. A1 - Hofreiter, Michael A1 - Fago, Angela A1 - Weber, Roy E. A1 - Springer, Mark S. A1 - Campbell, Kevin L. T1 - Emergence of a chimeric globin pseudogene and increased Hemoglobin Oxygen Affinity Underlie the evolution of aquatic specializations in Sirenia JF - Molecular biology and evolution N2 - 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. KW - ancient DNA KW - aquatic adaptation KW - gene conversion KW - hemoglobin KW - oxygen affinity KW - molecular evolution KW - myoglobin KW - neuroglobin KW - cytoglobin KW - pseudogene Y1 - 2019 U6 - https://doi.org/10.1093/molbev/msz044 SN - 0737-4038 SN - 1537-1719 VL - 36 IS - 6 SP - 1134 EP - 1147 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Chen, Shun-Gang A1 - Li, Ji A1 - Zhang, Fan A1 - Xiao, Bo A1 - Hu, Jia-Ming A1 - Cui, Yin-Qiu A1 - Hofreiter, Michael A1 - Hou, Xin-Dong A1 - Sheng, Gui-Lian A1 - Lai, Xu-Long A1 - Yuan, Jun-Xia T1 - Different maternal lineages revealed by ancient mitochondrial genome of Camelus bactrianus from China JF - Mitochondrial DNA Part A N2 - 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. KW - Camelus bactrianus KW - mitochondrial genome KW - ancient DNA KW - phylogeny KW - maternal lineages Y1 - 2019 U6 - https://doi.org/10.1080/24701394.2019.1659250 SN - 2470-1394 SN - 2470-1408 VL - 30 IS - 7 SP - 786 EP - 793 PB - Routledge, Taylor & Francis Group CY - Abingdon ER - TY - JOUR A1 - Casas-Marce, Mireia A1 - Marmesat, Elena A1 - Soriano, Laura A1 - Martinez-Cruz, Begona A1 - Lucena-Perez, Maria A1 - Nocete, Francisco A1 - Rodriguez-Hidalgo, Antonio A1 - Canals, Antoni A1 - Nadal, Jordi A1 - Detry, Cleia A1 - Bernaldez-Sanchez, Eloisa A1 - Fernandez-Rodriguez, Carlos A1 - Perez-Ripoll, Manuel A1 - Stiller, Mathias A1 - Hofreiter, Michael A1 - Rodriguez, Alejandro A1 - Revilla, Eloy A1 - Delibes, Miguel A1 - Godoy, Jose A. T1 - Spatiotemporal Dynamics of Genetic Variation in the Iberian Lynx along Its Path to Extinction Reconstructed with Ancient DNA JF - Molecular biology and evolution N2 - 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 = 10), historical (N = 245), and contemporary (N = 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. KW - Iberian lynx KW - ancient DNA KW - paleogenetics KW - genetic erosion KW - endangered species Y1 - 2017 U6 - https://doi.org/10.1093/molbev/msx222 SN - 0737-4038 SN - 1537-1719 VL - 34 SP - 2893 EP - 2907 PB - Oxford Univ. Press CY - Oxford ER -