Recharacterization of ancient DNA miscoding lesions: insights in the era of sequencing-by-synthesis

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.     

Preferential access to genetic information from endogenous hominin ancient DNA and accurate quantitative SNP-typing via SPEX

The analysis of targeted genetic loci from ancient, forensic and clinical samples is usually built upon polymerase chain reaction (PCR)-generated sequence data. However, many studies have shown that PCR amplification from poor-quality DNA templates can create sequence artefacts at significant levels. With hominin (human and other hominid) samples, the pervasive presence of highly PCR-amplifiable human DNA contaminants in the vast majority of samples can lead to the creation of recombinant hybrids and other non-authentic artefacts. The resulting PCR-generated sequences can then be difficult, if not impossible, to authenticate. In contrast, single primer extension (SPEX)-based approaches can genotype single nucleotide polymorphisms from ancient fragments of DNA as accurately as modern DNA. A single SPEX-type assay can amplify just one of the duplex DNA strands at target loci and generate a multi-fold depth-of-coverage, with non-authentic recombinant hybrids reduced to undetectable levels. Crucially, SPEX-type approaches can preferentially access genetic information from damaged and degraded endogenous ancient DNA templates over modern human DNA contaminants. The development of SPEX-type assays offers the potential for highly accurate, quantitative genotyping from ancient hominin samples.     

Improving ancient DNA read mapping against modern reference genomes

ABSTRACT: BACKGROUND: Next-Generation Sequencing has revolutionized our approach to ancient DNA (aDNA) research, by providing complete genomic sequences of ancient individuals and extinct species. However, the recovery of genetic material from long-dead organisms is still complicated by a number of issues, including post-mortem DNA damage and high levels of environmental contamination. Together with error profiles specific to the type of sequencing platforms used, these specificities could limit our ability to map sequencing reads against modern reference genomes and therefore limit our ability to identify endogenous ancient reads, reducing the efficiency of shotgun sequencing aDNA. RESULTS: In this study, we compare different computational methods for improving the accuracy and sensitivity of aDNA sequence identification, based on shotgun sequencing reads recovered from Pleistocene horse extracts using Illumina GAIIx and Helicos Heliscope platforms. We show that the performance of the Burrows Wheeler Aligner (BWA), that has been developed for mapping of undamaged sequencing reads using platforms with low rates of indel-types of sequencing errors, can be employed at acceptable run-times by modifying default parameters in a platform-specific manner. We also examine if trimming likely damaged positions at read ends can increase the recovery of genuine aDNA fragments and if accurate identification of human contamination can be achieved using a strategy previously suggested based on best hit filtering. We show that combining our different mapping and filtering approaches can increase the number of high-quality endogenous hits recovered by up to 33%. CONCLUSIONS: We have shown that Illumina and Helicos sequences recovered from aDNA extracts could not be aligned to modern reference genomes with the same efficiency unless mapping parameters are optimized for the specific types of errors generated by these platforms and by post-mortem DNA damage. Our findings have important implications for future aDNA research, as we define mapping guidelines that improve our ability to identify genuine aDNA sequences, which in turn could improve the genotyping accuracy of ancient specimens. Our framework provides a significant improvement to the standard procedures used for characterizing ancient genomes, which is challenged by contamination and often low amounts of DNA material.     

Tracking down human contamination in ancient human teeth

DNA contamination arising from the manipulation of ancient calcified tissue samples is a poorly understood, yet fundamental, problem that affects the reliability of ancient DNA (aDNA) studies. We have typed the mitochondrial DNA hypervariable region I of the only 6 people involved in the excavation, washing, and subsequent anthropological and genetic study of 23 Neolithic remains excavated from Granollers (Barcelona, Spain) and searched for their presence among the 572 clones generated during the aDNA analyses of teeth from these samples. Of the cloned sequences, 17.13% could be unambiguously identified as contaminants, with those derived from the people involved in the retrieval and washing of the remains present in higher frequencies than those of the anthropologist and genetic researchers. This finding confirms, for the first time, previous hypotheses that teeth samples are most susceptible to contamination at their initial excavation. More worrying, the cloned contaminant sequences exhibit substitutions that can be attributed to DNA damage after the contamination event, and we demonstrate that the level of such damage increases with time: contaminants that are >10 years old have approximately 5 times more damage than those that are recent. Furthermore, we demonstrate that in this data set, the damage rate of the old contaminant sequences is indistinguishable from that of the endogenous DNA sequences. As such, the commonly used argument that miscoding lesions observed among cloned aDNA sequences can be used to support data authenticity is misleading in scenarios where the presence of old contaminant sequences is possible. We argue therefore that the typing of those involved in the manipulation of the ancient human specimens is critical in order to ensure that generated results are accurate.     

Bayesian estimation of sequence damage in ancient DNA

DNA extracted from archaeological and paleontological remains is usually damaged by biochemical processes postmortem. Some of these processes lead to changes in the structure of the DNA molecule, which can result in the incorporation of incorrect nucleotides during polymerase chain reaction. These base misincorporations, or miscoding lesions, can lead to the inclusion of spurious additional mutations in ancient DNA (aDNA) data sets. This has the potential to affect the outcome of phylogenetic and population genetic analyses, including estimates of mutation rates and genetic diversity. We present a novel model, termed the delta model, which estimates the amount of damage in DNA data and accounts for its effects in a Bayesian phylogenetic framework. The ability of the delta model to estimate damage is first investigated using a simulation study. The model is then applied to 13 aDNA data sets. The amount of damage in these data sets is shown to be significant but low (about 1 damaged base per 750 nt), suggesting that precautions for limiting the influence of damaged sites, such as cloning and enzymatic treatment, are worthwhile. The results also suggest that relatively high rates of mutation previously estimated from aDNA data are not entirely an artifact of sequence damage and are likely to be due to other factors such as the persistence of transient polymorphisms. The delta model appears to be particularly useful for placing upper credibility limits on the amount of sequence damage in an alignment, and this capacity might be beneficial for future aDNA studies or for the estimation of sequencing errors in modern DNA.     

Application de la technique de PCR en temps réel à l’étude de l’ADN ancien

Real-time Polymerase Chain Reaction (PCR) for the study of ancient DNA. The properties of ancient DNA (aDNA) make difficult the retrieval of DNA sequence. The advantage of Real-Time PCR was exploited, for the first time, in the study of aDNA. We determined the optimal condition to amplify, in one round of PCR, aDNA, which should be directly sequenced. Beside the verification of aDNA authenticity, we compared two cleaning bone methods: scalpel and ethanol. The ethanol specimens showed the best DNA yield. The aDNA was extracted and amplified (mitochondrial hypervariable region I) from five skeletons exhumed from the archaeological site of Notre-Dame-du-Bourg (France), dated from 3rd to 17th century. To cite this article: R. Kefi et al., C. R. Palevol 2 (2003) 125–132.     

A Statistical Approach to Identify Ancient Template DNA

One of the key problems in the study of ancient DNA is that of authenticating sequences obtained from PCR amplifications of highly degraded samples. Contamination of ancient samples and postmortem damage to endogenous DNA templates are the major obstacles facing researchers in this task. In particular, the authentication of sequences obtained from ancient human remains is thought by many to be rather challenging. We propose a novel approach, based on the c statistic, that can be employed to help identify the sequence motif of an endogenous template, based on a sample of sequences that reflect the nucleotide composition of individual template molecules obtained from ancient tissues (such as cloned products from a PCR amplification). The c statistic exploits as information the most common form of postmortem damage observed among clone sequences in ancient DNA studies, namely, lesion-induced substitutions caused by cytosine deamination events. Analyses of simulated sets of templates with miscoding lesions and real sets of clone sequences from the literature indicate that the c-based approach is highly effective in identifying endogenous sequence motifs, even when they are not present among the sampled clones. The proposed approach is likely to be of general use to researchers working with DNA from ancient tissues, particularly from human remains, where authentication of results has been most challenging. [Reviewing Editor: Dr. Magnus Nordborg]     

Assessing the fidelity of ancient DNA sequences amplified from nuclear genes

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.     

The De Novo Assembly of Mitochondrial Genomes of the Extinct Passenger Pigeon (Ectopistes migratorius) with Next Generation Sequencing

The information from ancient DNA (aDNA) provides an unparalleled opportunity to infer phylogenetic relationships and population history of extinct species and to investigate genetic evolution directly. However, the degraded and fragmented nature of aDNA has posed technical challenges for studies based on conventional PCR amplification. In this study, we present an approach based on next generation sequencing to efficiently sequence the complete mitochondrial genome (mitogenome) of two extinct passenger pigeons (Ectopistes migratorius) using de novo assembly of massive short (90 bp), paired-end or single-end reads. Although varying levels of human contamination and low levels of postmortem nucleotide lesion were observed, they did not impact sequencing accuracy. Our results demonstrated that the de novo assembly of shotgun sequence reads could be a potent approach to sequence mitogenomes, and offered an efficient way to infer evolutionary history of extinct species.     

To clone or not to clone: method analysis for retrieving consensus sequences in ancient DNA samples

The challenges associated with the retrieval and authentication of ancient DNA (aDNA) evidence are principally due to post-mortem damage which makes ancient samples particularly prone to contamination from "modern" DNA sources. The necessity for authentication of results has led many aDNA researchers to adopt methods considered to be "gold standards" in the field, including cloning aDNA amplicons as opposed to directly sequencing them. However, no standardized protocol has emerged regarding the necessary number of clones to sequence, how a consensus sequence is most appropriately derived, or how results should be reported in the literature. In addition, there has been no systematic demonstration of the degree to which direct sequences are affected by damage or whether direct sequencing would provide disparate results from a consensus of clones.To address this issue, a comparative study was designed to examine both cloned and direct sequences amplified from ～3,500 year-old ancient northern fur seal DNA extracts. Majority rules and the Consensus Confidence Program were used to generate consensus sequences for each individual from the cloned sequences, which exhibited damage at 31 of 139 base pairs across all clones. In no instance did the consensus of clones differ from the direct sequence. This study demonstrates that, when appropriate, cloning need not be the default method, but instead, should be used as a measure of authentication on a case-by-case basis, especially when this practice adds time and cost to studies where it may be superfluous.     

一种实用可靠的古代 DNA 获取方法

古代DNA序列信息能够为物种演化研究提供最直接的分子证据,但获取古代DNA的技术仍存在诸多瓶颈,尤其是扩增中存在受损伤DNA模板的干扰、获取成本高和实验周期长等问题．改进了异丙醇沉淀提取法,并采用了尿嘧啶糖苷酶(UNG)去除受损伤DNA模板后进行扩增的方法,最终可以高效地获取真实的古代DNA序列．实验利用距今4 300～3 900年前的猪牙样本,将改进的古 DNA 获取方法与常规方法进行比较研究,结果表明,改进的异丙醇沉淀法提取结合UNG处理后进行PCR扩增的方法,可以在保证古代DNA获取成功率并提高获得的DNA序列可靠性的前提下,将经费投入和实验周期都各减少至常规方法的50%以下．这可以为开展大规模古代样本检测提供一种切实可行的 DNA 获取方法．     

Ancient genomics

The past decade has witnessed a revolution in ancient DNA (aDNA) research. Although the field''s focus was previously limited to mitochondrial DNA and a few nuclear markers, whole genome sequences from the deep past can now be retrieved. This breakthrough is tightly connected to the massive sequence throughput of next generation sequencing platforms and the ability to target short and degraded DNA molecules. Many ancient specimens previously unsuitable for DNA analyses because of extensive degradation can now successfully be used as source materials. Additionally, the analytical power obtained by increasing the number of sequence reads to billions effectively means that contamination issues that have haunted aDNA research for decades, particularly in human studies, can now be efficiently and confidently quantified. At present, whole genomes have been sequenced from ancient anatomically modern humans, archaic hominins, ancient pathogens and megafaunal species. Those have revealed important functional and phenotypic information, as well as unexpected adaptation, migration and admixture patterns. As such, the field of aDNA has entered the new era of genomics and has provided valuable information when testing specific hypotheses related to the past.     

Optimization of DNA Recovery and Amplification from Non-Carbonized Archaeobotanical Remains

Ancient DNA (aDNA) recovered from archaeobotanical remains can provide key insights into many prominent archaeological research questions, including processes of domestication, past subsistence strategies, and human interactions with the environment. However, it is often difficult to isolate aDNA from ancient plant materials, and furthermore, such DNA extracts frequently contain inhibitory substances that preclude successful PCR amplification. In the age of high-throughput sequencing, this problem is even more significant because each additional endogenous aDNA molecule improves analytical resolution. Therefore, in this paper, we compare a variety of DNA extraction techniques on primarily desiccated archaeobotanical remains and identify which method consistently yields the greatest amount of purified DNA. In addition, we test five DNA polymerases to determine how well they replicate DNA extracted from non-charred ancient plant remains. Based upon the criteria of resistance to enzymatic inhibition, behavior in quantitative real-time PCR, replication fidelity, and compatibility with aDNA damage, we conclude these polymerases have nuanced properties, requiring researchers to make educated decisions as to which one to use for a given task. The experimental findings should prove useful to the aDNA and archaeological communities by guiding future research methodologies and ensuring precious archaeobotanical remains are studied in optimal ways, and may thereby yield important new perspectives on the interactions between humans and past plant communities.     

Characterization of Nucleotide Misincorporation Patterns in the Iceman's Mitochondrial DNA

Background

The degradation of DNA represents one of the main issues in the genetic analysis of archeological specimens. In the recent years, a particular kind of post-mortem DNA modification giving rise to nucleotide misincorporation (“miscoding lesions”) has been the object of extensive investigations.

Methodology/Principal Findings

To improve our knowledge regarding the nature and incidence of ancient DNA nucleotide misincorporations, we have utilized 6,859 (629,975 bp) mitochondrial (mt) DNA sequences obtained from the 5,350–5,100-years-old, freeze-desiccated human mummy popularly known as the Tyrolean Iceman or Ötzi. To generate the sequences, we have applied a mixed PCR/pyrosequencing procedure allowing one to obtain a particularly high sequence coverage. As a control, we have produced further 8,982 (805,155 bp) mtDNA sequences from a contemporary specimen using the same system and starting from the same template copy number of the ancient sample. From the analysis of the nucleotide misincorporation rate in ancient, modern, and putative contaminant sequences, we observed that the rate of misincorporation is significantly lower in modern and putative contaminant sequence datasets than in ancient sequences. In contrast, type 2 transitions represent the vast majority (85%) of the observed nucleotide misincorporations in ancient sequences.

Conclusions/Significance

This study provides a further contribution to the knowledge of nucleotide misincorporation patterns in DNA sequences obtained from freeze-preserved archeological specimens. In the Iceman system, ancient sequences can be clearly distinguished from contaminants on the basis of nucleotide misincorporation rates. This observation confirms a previous identification of the ancient mummy sequences made on a purely phylogenetical basis. The present investigation provides further indication that the majority of ancient DNA damage is reflected by type 2 (cytosine→thymine/guanine→adenine) transitions and that type 1 transitions are essentially PCR artifacts.     

Poly(A) tailing of ancient DNA: a method for reproducible microsatellite genotyping

Microsatellites could be of great potential use in the analysis of ancient remains, but so far such analyses have failed to be reproducible mainly because of the high degree of ancient DNA (aDNA) degradation. During PCR, annealing of the primers to the complementary sequences of microsatellites occurs together with cross-annealing of partially degraded repeated sequences. This could create chimeric alleles that do not correspond to the authentic ones. Here we report a simple method for processing aDNA fragments prior to PCR that greatly reduces the production of chimeric alleles. This approach eliminates aDNA molecules broken within the repeats as targets for Taq polymerase by adding poly(A) tails at the 3(') ends of the DNA fragments, which disrupts the homology in the region and thus prevents annealing out of register. We have analyzed one dinucleotide- (D6S337) and two trinucleotide-containing loci (IT15 and SCA1) using poly(A)-tailed and the same untreated aDNA as template. aDNAs were isolated from 28 human remains, 600 and 7000 years of age. In repeated experiments with untreated aDNAs we obtained three to five times more alleles compared to poly(A)-tailed aDNAs. According to our results, modification of aDNA by poly(A) tailing is an efficient pretreatment for accurate genotyping.     

Estimating temporally variable selection intensity from ancient DNA data with the flexibility of modelling linkage and epistasis

Innovations in ancient DNA (aDNA) preparation and sequencing technologies have exponentially increased the quality and quantity of aDNA data extracted from ancient biological materials. The additional temporal component from the incoming aDNA data can provide improved power to address fundamental evolutionary questions like characterizing selection processes that shape the phenotypes and genotypes of contemporary populations or species. However, utilizing aDNA to study past selection processes still involves considerable hurdles like how to eliminate the confounding factor of genetic interactions in the inference of selection. To address this issue, we extend the approach of He et al., 2023 to infer temporally variable selection from the aDNA data in the form of genotype likelihoods with the flexibility of modelling linkage and epistasis in this work. Our posterior computation is carried out by a robust adaptive version of the particle marginal Metropolis-Hastings algorithm with a coerced acceptance rate. Our extension inherits the desirable features of He et al., 2023 such as modelling sample uncertainty resulting from the damage and fragmentation of aDNA molecules and reconstructing underlying gamete frequency trajectories of the population. We evaluate its performance through extensive simulations and show its utility with an application to the aDNA data from pigmentation loci in horses.     

Analysis of ancient DNA from coprolites: a perspective with random amplified polymorphic DNA-polymerase chain reaction approach

The aim of this work was to determine approaches that would improve the quality of ancient DNA (aDNA) present in coprolites to enhance the possibility of success in retrieving specific sequence targets. We worked with coprolites from South American archaeological sites in Brazil and Chile dating up to 7,000 years ago. Using established protocols for aDNA extraction we obtained samples showing high degradation as usually happens with this kind of material. The reconstructive polymerization pretreatment was essential to overcome the DNA degradation and the serial dilutions helped with to prevent polymerase chain reaction (PCR) inhibitors. Moreover, the random amplified polymorphic DNA-PCR has been shown to be a reliable technique for further experiments to recover specific aDNA sequences.     

Identification of ancient Olea europaea L. and Cornus mas L. seeds by DNA barcoding

The analysis of ancient DNA (aDNA) provides archaeologists and anthropologists with innovative, scientific and accurate data to study and understand the past. In this work, ancient seeds, found in the "Mora Cavorso" archaeological site (Latium, Central Italy), were analyzed to increase information about Italian Neolithic populations (plant use, agriculture, diet, trades, customs and ecology). We performed morphological and genetic techniques to identify fossil botanical species. In particular, this study also suggests and emphasizes the use of DNA barcode method for ancient plant sample analysis. Scanning electron microscope (SEM) observations showed seed compact structure and irregular surface but they did not permit a precise nor empirical classification: so, a molecular approach was necessary. DNA was extracted from ancient seeds and then it was used, as template, for PCR amplifications of standardized barcode genes. Although aDNA could be highly degraded by the time, successful PCR products were obtained, sequenced and compared to nucleotide sequence databases. Positive outcomes (supported by morphological comparison with modern seeds, geographical distribution and historical data) indicated that seeds could be identified as belonging to two plant species: Olea europaea L. and Cornus mas L.     

Adaptor-mediated amplification of minute amounts of severely fragmented ancient nucleic acids

It has been repeatedly shown that high copy number mitochondrial DNA sequences can be recovered from ancient samples. A significant increase in the volume of information available to researchers will be observed when the amplification of nuclear DNA becomes commonplace and reproducible. To this end we established a modification of the Rapid Amplification of cDNA Ends (RACE) procedure normally used for the generation of cDNA ends from adaptor-ligated expressed sequence tag libraries. The modifications were designed to specifically address the problems associated with the highly damaged nucleic acids extracted from palaeontological specimens. For this study we used 6 human samples dating to 450 AD and approximately 6.500 BP that were refractory to reliable amplification of single copy loci by PCR. Racemate contents (ratio of D/L enantiomers) of aspartic acid, alanine, and leucine also indicated that no amplifiable DNA is present in 5 of the 6 samples. The proposed technique allowed us (i) to amplify four X-chromosomal loci from 5 human specimens, and (ii) to correct allelic drop-out phenomena at the amelogenin locus in one individual; thus showing that the threshold of 80 x 10-3 for D/Lasp as a borderline for the presence/absence of amplifiable aDNA requires reassessment. Reliability of the proposed technique (i.e. amplification of DNA sequences endogenous to the find) was validated by the application of "ancient RACE" (aRACE) to prehistoric animal samples.     

How many clones need to be sequenced from a single forensic or ancient DNA sample in order to determine a reliable consensus sequence?

Forensic and ancient DNA (aDNA) extracts are mixtures of endogenous aDNA, existing in more or less damaged state, and contaminant DNA. To obtain the true aDNA sequence, it is not sufficient to generate a single direct sequence of the mixture, even where the authentic aDNA is the most abundant (e.g. 25% or more) in the component mixture. Only bacterial cloning can elucidate the components of this mixture. We calculate the number of clones that need to be sampled (for various mixture ratios) in order to be confident (at various levels of confidence) to have identified the major component. We demonstrate that to be >95% confident of identifying the most abundant sequence present at 70% in the ancient sample, 20 clones must be sampled. We make recommendations and offer a free-access web-based program, which constructs the most reliable consensus sequence from the user's input clone sequences and analyses the confidence limits for each nucleotide position and for the whole consensus sequence. Accepted authentication methods must be employed in order to assess the authenticity and endogeneity of the resulting consensus sequences (e.g. quantification and replication by another laboratory, blind testing, amelogenin sex versus morphological sex, the effective use of controls, etc.) and determine whether they are indeed aDNA.