Ancient DNA Sequences Research Articles

MapDamage: testing for damage patterns in ancient DNA sequences

Ancient DNA extracts consist of a mixture of contaminant DNA molecules, most often originating from environmental microbes, and endogenous fragments exhibiting substantial levels of DNA damage. The latter introduce specific nucleotide misincorporations and DNA fragmentation signatures in sequencing reads that could be advantageously used to argue for sequence validity. mapDamage is a Perl script that computes nucleotide misincorporation and fragmentation patterns using next-generation sequencing reads mapped against a reference genome. The Perl script outputs are further automatically processed in embedded R script in order to detect typical patterns of genuine ancient DNA sequences. The Perl script mapDamage is freely available with documentation and example files at http://geogenetics.ku.dk/all_literature/mapdamage/. The script requires prior installation of the SAMtools suite and R environment and has been validated on both GNU/Linux and MacOSX operating systems.

Development of ancient DNA techniques and its application in evolutionary study of domestic animals

Domestic animals have accompanied the origin and development of human civilization; therefore, the study of their domestication history is important. However, the occurrence of gene flow, or even replacement, among wild or domestic populations has complicated genetic studies tracing the initial domestication and dispersal of modern species. Fortunately, the investigation of DNA from archeological remains (ancient DNA) provides a unique opportunity to unravel these past events. In the present study, the progress of ancient DNA techniques is summarized, including retrieval and evaluation methods, technological innovations in methods of ancient DNA extraction, amplification, and sequencing. As well, the implications of ancient DNA technology for the evolutionary study of four major domestic animals (pig, horse, sheep, and dog) are reviewed. This review will facilitate the development of ancient DNA techniques and their further applications.

Bayesian Estimation of Substitution Rates from Ancient DNA Sequences with Low Information Content

Bayesian Estimation of Substitution Rates from Ancient DNA Sequences with Low Information Content

Late-Quaternary biogeographic scenarios for the brown bear ( Ursus arctos), a wild mammal model species

This review provides an up-to-date synthesis of the matrilineal phylogeography of a uniquely well-studied Holarctic mammal, the brown bear. We extend current knowledge by presenting a DNA sequence derived from one of the earliest known fossils of a polar bear (dated to 115 000 years before present), a species that shares a paraphyletic mitochondrial association with brown bears. A molecular clock analysis of 140 mitochondrial DNA sequences, including our new polar bear sequence, provides novel insights into the times of origin for different brown bear clades. We propose a number of regional biogeographic scenarios based on genetic data, divergence time estimates and paleontological records. The case of the brown bear provides an example for researchers working with less well-studied taxa: it shows clearly that phylogeographic models based on patterns of modern genetic variation alone can be substantially improved by including data on historical patterns of genetic diversity in the form of ancient DNA sequences derived from accurately dated samples and by using an approach to divergence-time estimation that suits the data under analysis. Using such approaches it has been possible to (i) establish that the processes shaping modern genetic diversity in brown bears acted recently, within the last three glacial cycles; (ii) distinguish among hypotheses concerning species’ responses to climatic oscillations in accordance with the lack of phylogeographic structure that existed in brown bears prior to the last glacial maximum (LGM); (iii) reassess theories linking monophyletic brown bear populations to particular LGM refuge areas; and (iv) identify vicariance events and track analogous patterns of migration by brown bears out of Eurasia to North America and Japan.

Ancient DNA sequences point to a large loss of mitochondrial genetic diversity in the saiga antelope (Saiga tatarica) since the Pleistocene

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.

New evidence for the establishment and management of the European fallow deer ( Dama dama dama) in Roman Britain

The semi-domestic status of the European fallow deer ( Dama dama dama) renders its ancient biogeography a reflection of human activity with the potential to provide important insights into the movement, trade patterns and ideology of past societies. Given this potential, fallow deer have received surprisingly little attention from scientists within the fields of archaeology, biology and zoology. Here we present new AMS radiocarbon dates, stable carbon and nitrogen isotope data and genetic evidence (the first ancient DNA sequences for the species) resulting from the analysis of a set of remains recovered from the Roman settlement at Monkton on the Isle of Thanet, Kent, UK. By viewing our results against the very limited available comparative data, this paper provides new information for the establishment and management of fallow deer in Britain. We argue that much more could be achieved with even a slight increase in sample sizes and a plea is made for greater research into this culturally significant species.

A rapid column‐based ancient DNA extraction method for increased sample throughput

Genetic analyses using museum specimens and ancient DNA from fossil samples are becoming increasingly important in phylogenetic and especially population genetic studies. Recent progress in ancient DNA sequencing technologies has substantially increased DNA sequence yields and, in combination with barcoding methods, has enabled large-scale studies using any type of DNA. Moreover, more and more studies now use nuclear DNA sequences in addition to mitochondrial ones. Unfortunately, nuclear DNA is, due to its much lower copy number in living cells compared to mitochondrial DNA, much more difficult to obtain from low-quality samples. Therefore, a DNA extraction method that optimizes DNA yields from low-quality samples and at the same time allows processing many samples within a short time frame is immediately required. In fact, the major bottleneck in the analysis process using samples containing low amounts of degraded DNA now lies in the extraction of samples, as column-based methods using commercial kits are fast but have proven to give very low yields, while more efficient methods are generally very time-consuming. Here, we present a method that combines the high DNA yield of batch-based silica extraction with the time-efficiency of column-based methods. Our results on Pleistocene cave bear samples show that DNA yields are quantitatively comparable, and in fact even slightly better than with silica batch extraction, while at the same time the number of samples that can conveniently be processed in parallel increases and both bench time and costs decrease using this method. Thus, this method is suited for harvesting the power of high-throughput sequencing using the DNA preserved in the millions of paleontological and museums specimens.

Tuberculosis on the north coast of Peru: skeletal and molecular paleopathology of late pre-Hispanic and postcontact mycobacterial disease

This paper examines skeletal and ancient DNA evidence in the study of suspected tuberculosis infection in the late pre-Hispanic and Colonial-era Lambayeque Valley Complex, north coast Peru (A.D. 900–1750). We integrate information on macroscopic lesion characteristics and distribution, radiographic and CT scan imagery, and analysis of Mycobacterium tuberculosis complex rpoB and IS6110 ancient DNA (aDNA) sequences. Destructive lesions were observed in the vertebral bodies of three precontact indigenous adult males, one colonial adolescent female, and in the cranium of a Colonial-period subadult. Assessment of lesion morphology and distribution led us to consider multiple diseases, but tuberculosis represents the most likely diagnostic option in all individuals. DNA was poorly preserved in all samples, but an IS6110 sequence was amplified in one precontact individual consistent with macroscopic diagnosis. These findings expand the geographic and temporal extent of tuberculosis to the late pre-Hispanic and Colonial north coast of Peru to highlight potential synergisms between diet, settlement patterns, and the evolution of Andean tuberculosis before and after European conquest. Moreover, this study helps focus several key questions in Andean tuberculosis research, including possible reassessment of the presence of the IS6110 sequence in the pre-Columbian Americas. Methodological considerations include differential diagnosis – especially with incomplete skeletons – and limitations of aDNA studies underscoring an approach integrating macroscopic, radiographic, and molecular lines of evidence in the paleopathological investigation of one of humankind’s most devastating and destabilizing diseases.

Mitochondrial DNA from 3000-year old chickens at the Teouma site, Vanuatu

Chickens were part of the Lapita cultural complex, transported into and through the Pacific by prehistoric colonists; as such they can be used as a proxy for tracking prehistoric migration and interaction. The Lapita site of Teouma in Vanuatu is well known for the recovery of complete dentate stamped pots and a cemetery containing the largest collection of Lapita period skeletons ever found. Chicken bones recovered from these excavations provide the first ancient DNA sequences from any commensal organism directly associated with a Lapita context. The ancient mtDNA sequences obtained from two Teouma chicken bones are compared with previously published archaeologically derived ancient DNA sequences to extend our understanding of the spread of chickens in Pacific prehistory. The results also show that the haplogroup E signature was present in very early populations of chickens transported into Remote Oceania. This study also adds to the suite of available data relating to isotopic signatures for commensal animals during the early settlement of Vanuatu and may reveal a different diet for Teouma chickens than those from other early prehistoric assemblages in the Pacific.

Human evolution: Stranger from Siberia.

The sequencing of ancient DNA is generating dramatic results. The sequence from a bone fragment has revealed the existence of an unknown type of extinct human ancestor that lived in Asia 40,000 years ago.

Ancient DNA analyses exclude humans as the driving force behind late Pleistocene musk ox ( Ovibos moschatus ) population dynamics

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.

Computational challenges in the analysis of ancient DNA

A new method of next-generation sequencing analysis is presented which takes into account the biases characteristic of ancient, including Neandertal, DNA samples.

A Complete mtDNA Genome of an Early Modern Human from Kostenki, Russia

The recovery of DNA sequences from early modern humans (EMHs) could shed light on their interactions with archaic groups such as Neandertals and their relationships to current human populations. However, such experiments are highly problematic because present-day human DNA frequently contaminates bones [1, 2]. For example, in a recent study of mitochondrial (mt) DNA from Neolithic European skeletons, sequence variants were only taken as authentic if they were absent or rare in the present population, whereas others had to be discounted as possible contamination [3, 4]. This limits analysis to EMH individuals carrying rare sequences and thus yields a biased view of the ancient gene pool. Other approaches of identifying contaminating DNA, such as genotyping all individuals who have come into contact with a sample, restrict analyses to specimens where this is possible [5, 6] and do not exclude all possible sources of contamination. By studying mtDNA in Neandertal remains, where contamination and endogenous DNA can be distinguished by sequence, we show that fragmentation patterns and nucleotide misincorporations can be used to gauge authenticity of ancient DNA sequences. We use these features to determine a complete mtDNA sequence from a approximately 30,000-year-old EMH from the Kostenki 14 site in Russia.

A Practical and Efficient Method for The Retrieval of Ancient DNA Sequence*

Retrieval of ancient DNA (aDNA) sequences from organism remains provide direct view of their evolutionary history.However,researches on aDNA have suffered from lots of technical problems.Specifically,discredited sequences were generated from damaged aDNA templates,and expensive and time-consuming methods were employed.Here,a method which could recover the endogenous aDNA as well as to reduce the cost and research period is described.This is achieved by improving the ancient DNA extraction method of isopropanol precipitation,and reevaluating the method of PCR after N-glycosylase (UNG) treatment,which could remove the damaged DNA from the aDNA extract.The efficiency of these methods were tested by comparing with traditional methods using ancient specimens of pig teeth aged between 4 300 years before present (BP) and 3 900 BP.The results showed that:firstly,the extraction efficiency of the improved method of isopropanol precipitation and current method with silica-based spin column were all 60%.Furthermore,the research period at least could be reduced by half with the application of the improved methods and the cost to 1/10 of the current method.Secondly,sequences obtained through the method of PCR after UNG treatment were 100% authentic.In contrast,66%～88% sequences were authentic based on the results obtained with the method of multiple PCRs without UNG treatment.And the research cost and period needed by the method with UNG treatment were only half of the later one.These results demonstrate that the improved extraction method of isopropanol precipitation combined with the method of PCR after UNG treatment could increase the success rate of authentic DNA amplified and at least reduce the research cost and period by half.Therefore,this method can be applied in the large-scale detection of ancient specimens.

Processed pseudogenes: the ‘fossilized footprints’ of past gene expression

Although our knowledge of the genes and genomes of extinct organisms is improving as a result of progress in sequencing ancient DNA, the transcriptomes of extinct organisms remain inaccessible, owing to the rapid degradation of messenger RNA after death. We provide empirical evidence that gene expression levels in the reproductive tissues of mice and during early mouse development correlate highly with the rate of inherited retroposition: the source of processed pseudogenes in the genome. Thus, processed pseudogenes might serve as fossilized footprints of the expression of their parent genes, shedding light on ancient transcriptomes that could provide significant insights into the evolution of gene expression.

Targeted Retrieval and Analysis of Five Neandertal mtDNA Genomes

Analysis of Neandertal DNA holds great potential for investigating the population history of this group of hominins, but progress has been limited due to the rarity of samples and damaged state of the DNA. We present a method of targeted ancient DNA sequence retrieval that greatly reduces sample destruction and sequencing demands and use this method to reconstruct the complete mitochondrial DNA (mtDNA) genomes of five Neandertals from across their geographic range. We find that mtDNA genetic diversity in Neandertals that lived 38,000 to 70,000 years ago was approximately one-third of that in contemporary modern humans. Together with analyses of mtDNA protein evolution, these data suggest that the long-term effective population size of Neandertals was smaller than that of modern humans and extant great apes.

An algebraic hypothesis about the primeval genetic code architecture

A plausible architecture of an ancient genetic code is derived from an extended base triplet vector space over the Galois field of the extended base alphabet {D, A, C, G, U}, where symbol D represents one or more hypothetical bases with unspecific pairings. We hypothesized that the high degeneration of a primeval genetic code with five bases and the gradual origin and improvement of a primeval DNA repair system could make possible the transition from ancient to modern genetic codes. Our results suggest that the Watson–Crick base pairing G ≡ C and A = U and the non-specific base pairing of the hypothetical ancestral base D used to define the sum and product operations are enough features to determine the coding constraints of the primeval and the modern genetic code, as well as, the transition from the former to the latter. Geometrical and algebraic properties of this vector space reveal that the present codon assignment of the standard genetic code could be induced from a primeval codon assignment. Besides, the Fourier spectrum of the extended DNA genome sequences derived from the multiple sequence alignment suggests that the called period-3 property of the present coding DNA sequences could also exist in the ancient coding DNA sequences. The phylogenetic analyses achieved with metrics defined in the N-dimensional vector space ( B 3 ) N of DNA sequences and with the new evolutionary model presented here also suggest that an ancient DNA coding sequence with five or more bases does not contradict the expected evolutionary history.

DNA content and distribution in ancient feathers and potential to reconstruct the plumage of extinct avian taxa

Feathers are known to contain amplifiable DNA at their base (calamus) and have provided an important genetic source from museum specimens. However, feathers in subfossil deposits generally only preserve the upper shaft and feather 'vane' which are thought to be unsuitable for DNA analysis. We analyse subfossil moa feathers from Holocene New Zealand rockshelter sites and demonstrate that both ancient DNA and plumage information can be recovered from their upper portion, allowing species identification and a means to reconstruct the appearance of extinct taxa. These ancient DNA sequences indicate that the distal portions of feathers are an untapped resource for studies of museum, palaeontological and modern specimens. We investigate the potential to reconstruct the plumage of pre-historically extinct avian taxa using subfossil remains, rather than assuming morphological uniformity with closely related extant taxa. To test the notion of colour persistence in subfossil feathers, we perform digital comparisons of feathers of the red-crowned parakeet (Cyanoramphus novaezelandiae novaezelandiae) excavated from the same horizons as the moa feathers, with modern samples. The results suggest that the coloration of the moa feathers is authentic, and computer software is used to perform plumage reconstructions of moa based on subfossil remains.

Genotyping human ancient mtDNA control and coding region polymorphisms with a multiplexed Single-Base-Extension assay: the singular maternal history of the Tyrolean Iceman.

BackgroundProgress in the field of human ancient DNA studies has been severely restricted due to the myriad sources of potential contamination, and because of the pronounced difficulty in identifying authentic results. Improving the robustness of human aDNA results is a necessary pre-requisite to vigorously testing hypotheses about human evolution in Europe, including possible admixture with Neanderthals. This study approaches the problem of distinguishing between authentic and contaminating sequences from common European mtDNA haplogroups by applying a multiplexed Single-Base-Extension assay, containing both control and coding region sites, to DNA extracted from the Tyrolean Iceman.ResultsThe multiplex assay developed for this study was able to confirm that the Iceman's mtDNA belongs to a new European mtDNA clade with a very limited distribution amongst modern data sets. Controlled contamination experiments show that the correct results are returned by the multiplex assay even in the presence of substantial amounts of exogenous DNA. The overall level of discrimination achieved by targeting both control and coding region polymorphisms in a single reaction provides a methodology capable of dealing with most cases of homoplasy prevalent in European haplogroups.ConclusionThe new genotyping results for the Iceman confirm the extreme fallibility of human aDNA studies in general, even when authenticated by independent replication. The sensitivity and accuracy of the multiplex Single-Base-Extension methodology forms part of an emerging suite of alternative techniques for the accurate retrieval of ancient DNA sequences from both anatomically modern humans and Neanderthals. The contamination of laboratories remains a pressing concern in aDNA studies, both in the pre and post-PCR environments, and the adoption of a forensic style assessment of a priori risks would significantly improve the credibility of results.

Using Classical Population Genetics Tools with Heterochroneous Data: Time Matters!

BackgroundNew polymorphism datasets from heterochroneous data have arisen thanks to recent advances in experimental and microbial molecular evolution, and the sequencing of ancient DNA (aDNA). However, classical tools for population genetics analyses do not take into account heterochrony between subsets, despite potential bias on neutrality and population structure tests. Here, we characterize the extent of such possible biases using serial coalescent simulations.Methodology/Principal FindingsWe first use a coalescent framework to generate datasets assuming no or different levels of heterochrony and contrast most classical population genetic statistics. We show that even weak levels of heterochrony (∼10% of the average depth of a standard population tree) affect the distribution of polymorphism substantially, leading to overestimate the level of polymorphism θ, to star like trees, with an excess of rare mutations and a deficit of linkage disequilibrium, which are the hallmark of e.g. population expansion (possibly after a drastic bottleneck). Substantial departures of the tests are detected in the opposite direction for more heterochroneous and equilibrated datasets, with balanced trees mimicking in particular population contraction, balancing selection, and population differentiation. We therefore introduce simple corrections to classical estimators of polymorphism and of the genetic distance between populations, in order to remove heterochrony-driven bias. Finally, we show that these effects do occur on real aDNA datasets, taking advantage of the currently available sequence data for Cave Bears (Ursus spelaeus), for which large mtDNA haplotypes have been reported over a substantial time period (22–130 thousand years ago (KYA)).Conclusions/SignificanceConsidering serial sampling changed the conclusion of several tests, indicating that neglecting heterochrony could provide significant support for false past history of populations and inappropriate conservation decisions. We therefore argue for systematically considering heterochroneous models when analyzing heterochroneous samples covering a large time scale.