@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} } @article{FoersterBullLenzetal.2018, author = {F{\"o}rster, Daniel W. and Bull, James K. and Lenz, Dorina and Autenrieth, Marijke and Paijmans, Johanna L. A. and Kraus, Robert H. S. and Nowak, Carsten and Bayerl, Helmut and K{\"u}hn, Ralph and Saveljev, Alexander P. and Sindicic, Magda and Hofreiter, Michael and Schmidt, Krzysztof and Fickel, J{\"o}rns}, title = {Targeted resequencing of coding DNA sequences for SNP discovery in nonmodel species}, series = {Molecular ecology resources}, volume = {18}, journal = {Molecular ecology resources}, number = {6}, publisher = {Wiley}, address = {Hoboken}, issn = {1755-098X}, doi = {10.1111/1755-0998.12924}, pages = {1356 -- 1373}, year = {2018}, abstract = {Targeted capture coupled with high-throughput sequencing can be used to gain information about nuclear sequence variation at hundreds to thousands of loci. Divergent reference capture makes use of molecular data of one species to enrich target loci in other (related) species. This is particularly valuable for nonmodel organisms, for which often no a priori knowledge exists regarding these loci. Here, we have used targeted capture to obtain data for 809 nuclear coding DNA sequences (CDS) in a nonmodel organism, the Eurasian lynx Lynx lynx, using baits designed with the help of the published genome of a related model organism (the domestic cat Felis catus). Using this approach, we were able to survey intraspecific variation at hundreds of nuclear loci in L. lynx across the species' European range. A large set of biallelic candidate SNPs was then evaluated using a high-throughput SNP genotyping platform (Fluidigm), which we then reduced to a final 96 SNP-panel based on assay performance and reliability; validation was carried out with 100 additional Eurasian lynx samples not included in the SNP discovery phase. The 96 SNP-panel developed from CDS performed very successfully in the identification of individuals and in population genetic structure inference (including the assignment of individuals to their source population). In keeping with recent studies, our results show that genic SNPs can be valuable for genetic monitoring of wildlife species.}, language = {en} } @article{BarlowCahillHartmannetal.2018, author = {Barlow, Axel and Cahill, James A. and Hartmann, Stefanie and Theunert, Christoph and Xenikoudakis, Georgios and Gonzalez-Fortes, Gloria M. and Paijmans, Johanna L. A. and Rabeder, Gernot and Frischauf, Christine and Garcia-Vazquez, Ana and Murtskhvaladze, Marine and Saarma, Urmas and Anijalg, Peeter and Skrbinsek, Tomaz and Bertorelle, Giorgio and Gasparian, Boris and Bar-Oz, Guy and Pinhasi, Ron and Slatkin, Montgomery and Dalen, Love and Shapiro, Beth and Hofreiter, Michael}, title = {Partial genomic survival of cave bears in living brown bears}, series = {Nature Ecology \& Evolution}, volume = {2}, journal = {Nature Ecology \& Evolution}, number = {10}, publisher = {Nature Publ. Group}, address = {London}, issn = {2397-334X}, doi = {10.1038/s41559-018-0654-8}, pages = {1563 -- 1570}, year = {2018}, abstract = {Although many large mammal species went extinct at the end of the Pleistocene epoch, their DNA may persist due to past episodes of interspecies admixture. However, direct empirical evidence of the persistence of ancient alleles remains scarce. Here, we present multifold coverage genomic data from four Late Pleistocene cave bears (Ursus spelaeus complex) and show that cave bears hybridized with brown bears (Ursus arctos) during the Pleistocene. We develop an approach to assess both the directionality and relative timing of gene flow. We find that segments of cave bear DNA still persist in the genomes of living brown bears, with cave bears contributing 0.9 to 2.4\% of the genomes of all brown bears investigated. Our results show that even though extinction is typically considered as absolute, following admixture, fragments of the gene pool of extinct species can survive for tens of thousands of years in the genomes of extant recipient species.}, 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} } @article{PaijmansBarlowFoersteretal.2018, author = {Paijmans, Johanna L. A. and Barlow, Axel and F{\"o}rster, Daniel W. and Henneberger, Kirstin and Meyer, Matthias and Nickel, Birgit and Nagel, Doris and Wors{\o}e Havm{\o}ller, Rasmus and Baryshnikov, Gennady F. and Joger, Ulrich and Rosendahl, Wilfried and Hofreiter, Michael}, title = {Historical biogeography of the leopard (Panthera pardus) and its extinct Eurasian populations}, series = {BMC Evolutionary Biology}, volume = {18}, journal = {BMC Evolutionary Biology}, number = {156}, publisher = {BioMed Central und Springer}, address = {London, Berlin und Heidelberg}, issn = {1471-2148}, doi = {10.1186/s12862-018-1268-0}, pages = {12}, year = {2018}, abstract = {Background Resolving the historical biogeography of the leopard (Panthera pardus) is a complex issue, because patterns inferred from fossils and from molecular data lack congruence. Fossil evidence supports an African origin, and suggests that leopards were already present in Eurasia during the Early Pleistocene. Analysis of DNA sequences however, suggests a more recent, Middle Pleistocene shared ancestry of Asian and African leopards. These contrasting patterns led researchers to propose a two-stage hypothesis of leopard dispersal out of Africa: an initial Early Pleistocene colonisation of Asia and a subsequent replacement by a second colonisation wave during the Middle Pleistocene. The status of Late Pleistocene European leopards within this scenario is unclear: were these populations remnants of the first dispersal, or do the last surviving European leopards share more recent ancestry with their African counterparts? Results In this study, we generate and analyse mitogenome sequences from historical samples that span the entire modern leopard distribution, as well as from Late Pleistocene remains. We find a deep bifurcation between African and Eurasian mitochondrial lineages (~ 710 Ka), with the European ancient samples as sister to all Asian lineages (~ 483 Ka). The modern and historical mainland Asian lineages share a relatively recent common ancestor (~ 122 Ka), and we find one Javan sample nested within these. Conclusions The phylogenetic placement of the ancient European leopard as sister group to Asian leopards suggests that these populations originate from the same out-of-Africa dispersal which founded the Asian lineages. The coalescence time found for the mitochondrial lineages aligns well with the earliest undisputed fossils in Eurasia, and thus encourages a re-evaluation of the identification of the much older putative leopard fossils from the region. The relatively recent ancestry of all mainland Asian leopard lineages suggests that these populations underwent a severe population bottleneck during the Pleistocene. Finally, although only based on a single sample, the unexpected phylogenetic placement of the Javan leopard could be interpreted as evidence for exchange of mitochondrial lineages between Java and mainland Asia, calling for further investigation into the evolutionary history of this subspecies.}, language = {en} }