Search

Most ancient specimens contain very low levels of endogenous DNA, precluding the shotgun sequencing of many interesting samples because of cost. Ancient DNA (aDNA) libraries often contain <1% endogenous DNA, with the majority of sequencing capacity taken up by environmental DNA. Here we present a capture-based method for enriching the endogenous component of aDNA sequencing libraries. By using biotinylated RNA baits transcribed from genomic DNA libraries, we are able to capture DNA fragments from across the human genome. We demonstrate this method on libraries created from four Iron Age and Bronze Age human teeth from Bulgaria, as well as bone samples from seven Peruvian mummies and a Bronze Age hair sample from Denmark. Prior to capture, shotgun sequencing of these libraries yielded an average of 1.2% of reads mapping to the human genome (including duplicates). After capture, this fraction increased substantially, with up to 59% of reads mapped to human and enrichment ranging from 6- to 159-fold. Furthermore, we maintained coverage of the majority of regions sequenced in the precapture library. Intersection with the 1000 Genomes Project reference panel yielded an average of 50,723 SNPs (range 3,062-147,243) for the postcapture libraries sequenced with 1 million reads, compared with 13,280 SNPs (range 217-73,266) for the precapture libraries, increasing resolution in population genetic analyses. Our whole-genome capture approach makes it less costly to sequence aDNA from specimens containing very low levels of endogenous DNA, enabling the analysis of larger numbers of samples., (Copyright © 2013 The American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)

Claims of extreme survival of DNA have emphasized the need for reliable models of DNA degradation through time. By analysing mitochondrial DNA (mtDNA) from 158 radiocarbon-dated bones of the extinct New Zealand moa, we confirm empirically a long-hypothesized exponential decay relationship. The average DNA half-life within this geographically constrained fossil assemblage was estimated to be 521 years for a 242 bp mtDNA sequence, corresponding to a per nucleotide fragmentation rate (k) of 5.50 × 10(-6) per year. With an effective burial temperature of 13.1°C, the rate is almost 400 times slower than predicted from published kinetic data of in vitro DNA depurination at pH 5. Although best described by an exponential model (R(2) = 0.39), considerable sample-to-sample variance in DNA preservation could not be accounted for by geologic age. This variation likely derives from differences in taphonomy and bone diagenesis, which have confounded previous, less spatially constrained attempts to study DNA decay kinetics. Lastly, by calculating DNA fragmentation rates on Illumina HiSeq data, we show that nuclear DNA has degraded at least twice as fast as mtDNA. These results provide a baseline for predicting long-term DNA survival in bone.

The genetic background of the European Mesolithic and the extent of population replacement during the Neolithic is poorly understood, both due to the scarcity of human remains from that period and the inherent methodological difficulties of ancient DNA research. However, advances in sequencing technologies are both increasing data yields and providing supporting evidence for data authenticity, such as nucleotide misincorporation patterns. We use these methods to characterize both the mitochondrial DNA genome and generate shotgun genomic data from two exceptionally well-preserved 7,000-year-old Mesolithic individuals from La Braña-Arintero site in León (Northwestern Spain). The mitochondria of both individuals are assigned to U5b2c1, a haplotype common among the small number of other previously studied Mesolithic individuals from Northern and Central Europe. This suggests a remarkable genetic uniformity and little phylogeographic structure over a large geographic area of the pre-Neolithic populations. Using Approximate Bayesian Computation, a model of genetic continuity from Mesolithic to Neolithic populations is poorly supported. Furthermore, analyses of 1.34% and 0.53% of their nuclear genomes, containing about 50,000 and 20,000 ancestry informative SNPs, respectively, show that these two Mesolithic individuals are not related to current populations from either the Iberian Peninsula or Southern Europe., (Copyright © 2012 Elsevier Ltd. All rights reserved.)

The Paisley Caves in Oregon record the oldest directly dated human remains (DNA) in the Western Hemisphere. More than 100 high-precision radiocarbon dates show that deposits containing artifacts and coprolites ranging in age from 12,450 to 2295 (14)C years ago are well stratified. Western Stemmed projectile points were recovered in deposits dated to 11,070 to 11,340 (14)C years ago, a time contemporaneous with or preceding the Clovis technology. There is no evidence of diagnostic Clovis technology at the site. These two distinct technologies were parallel developments, not the product of a unilinear technological evolution. "Blind testing" analysis of coprolites by an independent laboratory confirms the presence of human DNA in specimens of pre-Clovis age. The colonization of the Americas involved multiple technologically divergent, and possibly genetically divergent, founding groups.

Nunataks are isolated bedrocks protruding through ice sheets. They vary in age, but represent island environments in 'oceans' of ice through which organism dispersals and replacements can be studied over time. The J.A.D. Jensen's Nunataks at the southern Greenland ice sheet are the most isolated nunataks on the northern hemisphere - some 30 km from the nearest biological source. They constitute around 2 km(2) of ice-free land that was established in the early Holocene. We have investigated the changes in plant composition at these nunataks using both the results of surveys of the flora over the last 130 years and through reconstruction of the vegetation from the end of the Holocene Thermal Maximum (5528 ± 75 cal year BP) using meta-barcoding of plant DNA recovered from the nunatak sediments (sedaDNA). Our results show that several of the plant species detected with sedaDNA are described from earlier vegetation surveys on the nunataks (in 1878, 1967 and 2009). In 1967, a much higher biodiversity was detected than from any other of the studied periods. While this may be related to differences in sampling efforts for the oldest period, it is not the case when comparing the 1967 and 2009 levels where the botanical survey was exhaustive. As no animals and humans are found on the nunataks, this change in diversity over a period of just 42 years must relate to environmental changes probably being climate-driven. This suggests that even the flora of fairly small and isolated ice-free areas reacts quickly to a changing climate., (© 2011 Blackwell Publishing Ltd.)

Although ancient DNA from sediments (sedaDNA) has been used to investigate past ecosystems, the approach has never been directly compared with the traditional methods of pollen and macrofossil analysis. We conducted a comparative survey of 18 ancient permafrost samples spanning the Late Pleistocene (46-12.5 thousand years ago), from the Taymyr Peninsula in northern Siberia. The results show that pollen, macrofossils and sedaDNA are complementary rather than overlapping and, in combination, reveal more detailed information on plant palaeocommunities than can be achieved by each individual approach. SedaDNA and macrofossils share greater overlap in plant identifications than with pollen, suggesting that sedaDNA is local in origin. These two proxies also permit identification to lower taxonomic levels than pollen, enabling investigation into temporal changes in species composition and the determination of indicator species to describe environmental changes. Combining data from all three proxies reveals an area continually dominated by a mosaic vegetation of tundra-steppe, pioneer and wet-indicator plants. Such vegetational stability is unexpected, given the severe climate changes taking place in the Northern Hemisphere during this time, with changes in average annual temperatures of >22 °C. This may explain the abundance of ice-age mammals such as horse and bison in Taymyr Peninsula during the Pleistocene and why it acted as a refugium for the last mainland woolly mammoth. Our finding reveals the benefits of combining sedaDNA, pollen and macrofossil for palaeovegetational reconstruction and adds to the increasing evidence suggesting large areas of the Northern Hemisphere remained ecologically stable during the Late Pleistocene.

It is commonly believed that trees were absent in Scandinavia during the last glaciation and first recolonized the Scandinavian Peninsula with the retreat of its ice sheet some 9000 years ago. Here, we show the presence of a rare mitochondrial DNA haplotype of spruce that appears unique to Scandinavia and with its highest frequency to the west-an area believed to sustain ice-free refugia during most of the last ice age. We further show the survival of DNA from this haplotype in lake sediments and pollen of Trøndelag in central Norway dating back ~10,300 years and chloroplast DNA of pine and spruce in lake sediments adjacent to the ice-free Andøya refugium in northwestern Norway as early as ~22,000 and 17,700 years ago, respectively. Our findings imply that conifer trees survived in ice-free refugia of Scandinavia during the last glaciation, challenging current views on survival and spread of trees as a response to climate changes.

We used high-sensitivity, high-resolution tandem mass spectrometry to shotgun sequence ancient protein remains extracted from a 43 000 year old woolly mammoth ( Mammuthus primigenius ) bone preserved in the Siberian permafrost. For the first time, 126 unique protein accessions, mostly low-abundance extracellular matrix and plasma proteins, were confidently identified by solid molecular evidence. Among the best characterized was the carrier protein serum albumin, presenting two single amino acid substitutions compared to extant African ( Loxodonta africana ) and Indian ( Elephas maximus ) elephants. Strong evidence was observed of amino acid modifications due to post-mortem hydrolytic and oxidative damage. A consistent subset of this permafrost bone proteome was also identified in more recent Columbian mammoth ( Mammuthus columbi ) samples from temperate latitudes, extending the potential of the approach described beyond subpolar environments. Mass spectrometry-based ancient protein sequencing offers new perspectives for future molecular phylogenetic inference and physiological studies on samples not amenable to ancient DNA investigation. This approach therefore represents a further step into the ongoing integration of different high-throughput technologies for identification of ancient biomolecules, unleashing the field of paleoproteomics.

We present the first set of microsatellite markers developed exclusively for an extinct taxon. Microsatellite data have been analysed in thousands of genetic studies on extant species but the technology can be problematic when applied to low copy number (LCN) DNA. It is therefore rarely used on substrates more than a few decades old. Now, with the primers and protocols presented here, microsatellite markers are available to study the extinct New Zealand moa (Aves: Dinornithiformes) and, as with single nucleotide polymorphism (SNP) technology, the markers represent a means by which the field of ancient DNA can (preservation allowing) move on from its reliance on mitochondrial DNA. Candidate markers were identified using high throughput sequencing technology (GS-FLX) on DNA extracted from fossil moa bone and eggshell. From the 'shotgun' reads, >60 primer pairs were designed and tested on DNA from bones of the South Island giant moa (Dinornis robustus). Six polymorphic loci were characterised and used to assess measures of genetic diversity. Because of low template numbers, typical of ancient DNA, allelic dropout was observed in 36-70% of the PCR reactions at each microsatellite marker. However, a comprehensive survey of allelic dropout, combined with supporting quantitative PCR data, allowed us to establish a set of criteria that maximised data fidelity. Finally, we demonstrated the viability of the primers and the protocols, by compiling a full Dinornis microsatellite dataset representing fossils of c. 600-5000 years of age. A multi-locus genotype was obtained from 74 individuals (84% success rate), and the data showed no signs of being compromised by allelic dropout. The methodology presented here provides a framework by which to generate and evaluate microsatellite data from samples of much greater antiquity than attempted before, and opens new opportunities for ancient DNA research.

The remains of 12 Neandertal individuals have been found at the El Sidrón site (Asturias, Spain), consisting of six adults, three adolescents, two juveniles, and one infant. Archaeological, paleontological, and geological evidence indicates that these individuals represent all or part of a contemporaneous social group of Neandertals, who died at around the same time and later were buried together as a result of a collapse of an underground karst. We sequenced phylogenetically informative positions of mtDNA hypervariable regions 1 and 2 from each of the remains. Our results show that the 12 individuals stem from three different maternal lineages, accounting for seven, four, and one individual(s), respectively. Using a Y-chromosome assay to confirm the morphological determination of sex for each individual, we found that, although the three adult males carried the same mtDNA lineage, each of the three adult females carried different mtDNA lineages. These findings provide evidence to indicate that Neandertal groups not only were small and characterized by low genetic diversity but also were likely to have practiced patrilocal mating behavior.

Prior to the Holocene, the range of the saiga antelope (Saiga tatarica) spanned from France to the Northwest Territories of Canada. Although its distribution subsequently contracted to the steppes of Central Asia, historical records indicate that it remained extremely abundant until the end of the Soviet Union, after which its populations were reduced by over 95%. We have analysed the mitochondrial control region sequence variation of 27 ancient and 38 modern specimens, to assay how the species' genetic diversity has changed since the Pleistocene. Phylogenetic analyses reveal the existence of two well-supported, and clearly distinct, clades of saiga. The first, spanning a time range from >49,500 (14) C ybp to the present, comprises all the modern specimens and ancient samples from the Northern Urals, Middle Urals and Northeast Yakutia. The second clade is exclusive to the Northern Urals and includes samples dating from between 40,400 to 10,250 (14) C ybp. Current genetic diversity is much lower than that present during the Pleistocene, an observation that data modelling using serial coalescent indicates cannot be explained by genetic drift in a population of constant size. Approximate Bayesian Computation analyses show the observed data is more compatible with a drastic population size reduction (c. 66-77%) following either a demographic bottleneck in the course of the Holocene or late Pleistocene, or a geographic fragmentation (followed by local extinction of one subpopulation) at the Holocene/Pleistocene transition., (© 2010 Blackwell Publishing Ltd.)

Owing to exceptional biomolecule preservation, fossil avian eggshell has been used extensively in geochronology and palaeodietary studies. Here, we show, to our knowledge, for the first time that fossil eggshell is a previously unrecognized source of ancient DNA (aDNA). We describe the successful isolation and amplification of DNA from fossil eggshell up to 19 ka old. aDNA was successfully characterized from eggshell obtained from New Zealand (extinct moa and ducks), Madagascar (extinct elephant birds) and Australia (emu and owl). Our data demonstrate excellent preservation of the nucleic acids, evidenced by retrieval of both mitochondrial and nuclear DNA from many of the samples. Using confocal microscopy and quantitative PCR, this study critically evaluates approaches to maximize DNA recovery from powdered eggshell. Our quantitative PCR experiments also demonstrate that moa eggshell has approximately 125 times lower bacterial load than bone, making it a highly suitable substrate for high-throughput sequencing approaches. Importantly, the preservation of DNA in Pleistocene eggshell from Australia and Holocene deposits from Madagascar indicates that eggshell is an excellent substrate for the long-term preservation of DNA in warmer climates. The successful recovery of DNA from this substrate has implications in a number of scientific disciplines; most notably archaeology and palaeontology, where genotypes and/or DNA-based species identifications can add significantly to our understanding of diets, environments, past biodiversity and evolutionary processes.

The causes of the late Pleistocene megafaunal extinctions are poorly understood. Different lines of evidence point to climate change, the arrival of humans, or a combination of these events as the trigger. Although many species went extinct, others, such as caribou and bison, survived to the present. The musk ox has an intermediate story: relatively abundant during the Pleistocene, it is now restricted to Greenland and the Arctic Archipelago. In this study, we use ancient DNA sequences, temporally unbiased summary statistics, and Bayesian analytical techniques to infer musk ox population dynamics throughout the late Pleistocene and Holocene. Our results reveal that musk ox genetic diversity was much higher during the Pleistocene than at present, and has undergone several expansions and contractions over the past 60,000 years. Northeast Siberia was of key importance, as it was the geographic origin of all samples studied and held a large diverse population until local extinction at approximately 45,000 radiocarbon years before present ((14)C YBP). Subsequently, musk ox genetic diversity reincreased at ca. 30,000 (14)C YBP, recontracted at ca. 18,000 (14)C YBP, and finally recovered in the middle Holocene. The arrival of humans into relevant areas of the musk ox range did not affect their mitochondrial diversity, and both musk ox and humans expanded into Greenland concomitantly. Thus, their population dynamics are better explained by a nonanthropogenic cause (for example, environmental change), a hypothesis supported by historic observations on the sensitivity of the species to both climatic warming and fluctuations.

The mitochondrial genome has been the traditional focus of most research into ancient DNA, owing to its high copy number and population-level variability. Despite this long-standing interest in mitochondrial DNA, it was only in 2001 that the first complete ancient mitogenomic sequences were obtained. As a result of various methodological developments, including the introduction of high-throughput sequencing techniques, the total number of ancient mitogenome sequences has increased rapidly over the past few years. In this review, we present a brief history of ancient mitogenomics and describe the technical challenges that face researchers in the field. We catalogue the diverse sequencing methods and source materials used to obtain ancient mitogenomic sequences, summarise the associated genetic and phylogenetic studies that have been conducted, and evaluate the future prospects of the field.

Genetic variation in microsatellites is rarely examined in the field of ancient DNA (aDNA) due to the low quantity of nuclear DNA in the fossil record together with the lack of characterized nuclear markers in extinct species. 454 sequencing platforms provide a new high-throughput technology capable of generating up to 1 gigabases per run as short (200-400-bp) read lengths. 454 data were generated from the fossil bone of an extinct New Zealand moa (Aves: Dinornithiformes). We identified numerous short tandem repeat (STR) motifs, and here present the successful isolation and characterization of one polymorphic microsatellite (Moa&#95;MS2). Primers designed to flank this locus amplified all three moa species tested here. The presented method proved to be a fast and efficient way of identifying microsatellite markers in ancient DNA templates and, depending on biomolecule preservation, has the potential of enabling high-resolution population genetic studies of extinct taxa. As sequence read lengths of the 454 platforms and its competitors (e.g., the SOLEXA and SOLiD platforms) increase, this approach will become increasingly powerful in identifying microsatellites in extinct (and extant) organisms, and will afford new opportunities to study past biodiversity and extinction processes.

The timing of the first human migration into the Americas and its relation to the appearance of the Clovis technological complex in North America at about 11,000 to 10,800 radiocarbon years before the present (14C years B.P.) remains contentious. We establish that humans were present at Paisley 5 Mile Point Caves, in south-central Oregon, by 12,300 14C years B.P., through the recovery of human mitochondrial DNA (mtDNA) from coprolites, directly dated by accelerator mass spectrometry. The mtDNA corresponds to Native American founding haplogroups A2 and B2. The dates of the coprolites are >1000 14C years earlier than currently accepted dates for the Clovis complex.

We used authentication tests developed for ancient DNA to evaluate claims by Asara et al. (Reports, 13 April 2007, p. 280) of collagen peptide sequences recovered from mastodon and Tyrannosaurus rex fossils. Although the mastodon samples pass these tests, absence of amino acid composition data, lack of evidence for peptide deamidation, and association of alpha1(I) collagen sequences with amphibians rather than birds suggest that T. rex does not.

Poulakakis and colleagues (Poulakakis et al. 2006: Biol. Lett. 2, 451-454), report the recovery of 'authentic' mammoth DNA from an 800,000-year-old fragment of bone excavated on the island of Crete. In light of results from other ancient DNA studies that indicate how DNA survival is unlikely in samples, which are recovered from warm environments and are relatively old (e.g. more than 100,000 years), these findings come as a great surprise. Here, we show that problems exist with the methodological approaches used in the study. First, the nested PCR technique as reported is nonsensical--one of the second round 'nested' primers falls outside the amplicon of the first round PCR. More worryingly, the binding region of one of the first round primers (Elcytb320R) falls within the short 43 base pair reported mammoth sequence, specifically covering two of the three reportedly diagnostic Elephas polymorphisms. Finally, we demonstrate using a simple BLAST search in GenBank that the claimed 'uniquely derived character state' for mammoths is in fact also found within modern elephants.

Although ancient DNA (aDNA) miscoding lesions have been studied since the earliest days of the field, their nature remains a source of debate. A variety of conflicting hypotheses exist about which miscoding lesions constitute true aDNA damage as opposed to PCR polymerase amplification error. Furthermore, considerable disagreement and speculation exists on which specific damage events underlie observed miscoding lesions. The root of the problem is that it has previously been difficult to assemble sufficient data to test the hypotheses, and near-impossible to accurately determine the specific strand of origin of observed damage events. With the advent of emulsion-based clonal amplification (emPCR) and the sequencing-by-synthesis technology this has changed. In this paper we demonstrate how data produced on the Roche GS20 genome sequencer can determine miscoding lesion strands of origin, and subsequently be interpreted to enable characterization of the aDNA damage behind the observed phenotypes. Through comparative analyses on 390,965 bp of modern chloroplast and 131,474 bp of ancient woolly mammoth GS20 sequence data we conclusively demonstrate that in this sample at least, a permafrost preserved specimen, Type 2 (cytosine-->thymine/guanine-->adenine) miscoding lesions represent the overwhelming majority of damage-derived miscoding lesions. Additionally, we show that an as yet unidentified guanine-->adenine analogue modification, not the conventionally argued cytosine-->uracil deamination, underpins a significant proportion of Type 2 damage. How widespread these implications are for aDNA will become apparent as future studies analyse data recovered from a wider range of substrates.

To date, the field of ancient DNA has relied almost exclusively on mitochondrial DNA (mtDNA) sequences. However, a number of recent studies have reported the successful recovery of ancient nuclear DNA (nuDNA) sequences, thereby allowing the characterization of genetic loci directly involved in phenotypic traits of extinct taxa. It is well documented that postmortem damage in ancient mtDNA can lead to the generation of artifactual sequences. However, as yet no one has thoroughly investigated the damage spectrum in ancient nuDNA. By comparing clone sequences from 23 fossil specimens, recovered from environments ranging from permafrost to desert, we demonstrate the presence of miscoding lesion damage in both the mtDNA and nuDNA, resulting in insertion of erroneous bases during amplification. Interestingly, no significant differences in the frequency of miscoding lesion damage are recorded between mtDNA and nuDNA despite great differences in cellular copy numbers. For both mtDNA and nuDNA, we find significant positive correlations between total sequence heterogeneity and the rates of type 1 transitions (adenine --> guanine and thymine --> cytosine) and type 2 transitions (cytosine --> thymine and guanine --> adenine), respectively. Type 2 transitions are by far the most dominant and increase relative to those of type 1 with damage load. The results suggest that the deamination of cytosine (and 5-methyl cytosine) to uracil (and thymine) is the main cause of miscoding lesions in both ancient mtDNA and nuDNA sequences. We argue that the problems presented by postmortem damage, as well as problems with contamination from exogenous sources of conserved nuclear genes, allelic variation, and the reliance on single nucleotide polymorphisms, call for great caution in studies relying on ancient nuDNA sequences.

Ancient DNA studies represent a powerful tool that can be used to obtain genetic insights into the past. However, despite the publication of large numbers of apparently successful ancient DNA studies, a number of problems exist with the field that are often ignored. Therefore, questions exist as to how reliable the conclusions of many of the published studies are. In this paper we outline first the problems associated with aDNA studies, and secondly present potential guidelines designed so as to enable non-specialist readers the opportunity to critically assess the quality of aDNA publications.

The diversity of fungi in permanently frozen soil from northeastern Siberia was studied by culture-independent PCR amplification of diverse environmental 18S rRNA genes. Elaborate protocols to avoid contamination during drilling, sampling, and amplification were used. A broad diversity of eukaryotic DNA sequences that were 510 bp long, including sequences of various fungi, plants, and invertebrates, could be obtained reproducibly from samples that were up to 300,000 to 400,000 years old. The sequences revealed that ancient fungal communities included a diversity of cold-adapted yeasts, dark-pigmented fungi, plant-parasitic fungi, and lichen mycobionts. DNA traces of tree-associated macrofungi in a modern tundra sample indicated that there was a shift in fungal diversity following the last ice age and supported recent results showing that there was a severe change in the plant composition in northeastern Siberia during this period. Interestingly, DNA sequences with high homology to sequences of coprophilic and keratinophilic fungi indicated that feces, hair, skin, and nails could have been sources of ancient megafauna DNA recently reported to be present in small amounts of Siberian permafrost sediments.

The widespread extinctions of large mammals at the end of the Pleistocene epoch have often been attributed to the depredations of humans; here we present genetic evidence that questions this assumption. We used ancient DNA and Bayesian techniques to reconstruct a detailed genetic history of bison throughout the late Pleistocene and Holocene epochs. Our analyses depict a large diverse population living throughout Beringia until around 37,000 years before the present, when the population's genetic diversity began to decline dramatically. The timing of this decline correlates with environmental changes associated with the onset of the last glacial cycle, whereas archaeological evidence does not support the presence of large populations of humans in Eastern Beringia until more than 15,000 years later.

This study helps to clarify the debate on the Western and Eastern genetic influences in Central Asia. Thirty-six skeletal remains from Kazakhstan (Central Asia), excavated from different sites dating between the fifteenth century BC to the fifth century AD, have been analysed for the hypervariable control region (HVR-I) and haplogroup diagnostic single nucleotide polymorphisms (SNPs) of the mitochondrial DNA genome. Standard authentication criteria for ancient DNA studies, including multiple extractions, cloning of PCR products and independent replication, have been followed. The distribution of east and west Eurasian lineages through time in the region is concordant with the available archaeological information: prior to the thirteenth-seventh century BC, all Kazakh samples belong to European lineages; while later an arrival of east Eurasian sequences that coexisted with the previous west Eurasian genetic substratum can be detected. The presence of an ancient genetic substratum of European origin in West Asia may be related to the discovery of ancient mummies with European features in Xinjiang and to the existence of an extinct Indo-European language, Tocharian. This study demonstrates the usefulness of the ancient DNA in unravelling complex patterns of past human migrations so as to help decipher the origin of present-day admixed populations.