Extensive human DNA contamination in extracts from ancient dog bones and teeth

Ancient DNA (aDNA) sequences, especially those of human origin, are notoriously difficult to analyze due to molecular damage and exogenous DNA contamination. Relatively few systematic studies have focused on this problem. Here we investigate the extent and origin of human DNA contamination in the most frequently used sources for aDNA studies, that is, bones and teeth from museum collections. To distinguish contaminant DNA from authentic DNA we extracted DNA from dog (Canis familiaris) specimens. We monitored the presence of a 148-bp human-specific and a 152-bp dog-specific mitochondrial DNA (mtDNA) fragment in DNA extracts as well as in negative controls. The total number of human and dog template molecules were quantified using real-time polymerase chain reaction (PCR), and the sequences were characterized by amplicon cloning and sequencing. Although standard precautions to avoid contamination were taken, we found that all samples from the 29 dog specimens contained human DNA, often at levels exceeding the amount of authentic ancient dog DNA. The level of contaminating human DNA was also significantly higher in the dog extracts than in the negative controls, and an experimental setup indicated that this was not caused by the carrier effect. This suggests that the contaminating human DNA mainly originated from the dog bones rather than from laboratory procedures. When cloned, fragments within a contaminated PCR product generally displayed several different sequences, although one haplotype was often found in majority. This leads us to believe that recognized criteria for authenticating aDNA cannot separate contamination from ancient human DNA the way they are presently used.     

Fragmentation of contaminant and endogenous DNA in ancient samples determined by shotgun sequencing; prospects for human palaeogenomics

BACKGROUND: Despite the successful retrieval of genomes from past remains, the prospects for human palaeogenomics remain unclear because of the difficulty of distinguishing contaminant from endogenous DNA sequences. Previous sequence data generated on high-throughput sequencing platforms indicate that fragmentation of ancient DNA sequences is a characteristic trait primarily arising due to depurination processes that create abasic sites leading to DNA breaks. METHODOLOGY/PRINCIPALS FINDINGS: To investigate whether this pattern is present in ancient remains from a temperate environment, we have 454-FLX pyrosequenced different samples dated between 5,500 and 49,000 years ago: a bone from an extinct goat (Myotragus balearicus) that was treated with a depurinating agent (bleach), an Iberian lynx bone not subjected to any treatment, a human Neolithic sample from Barcelona (Spain), and a Neandertal sample from the El Sidrón site (Asturias, Spain). The efficiency of retrieval of endogenous sequences is below 1% in all cases. We have used the non-human samples to identify human sequences (0.35 and 1.4%, respectively), that we positively know are contaminants. CONCLUSIONS: We observed that bleach treatment appears to create a depurination-associated fragmentation pattern in resulting contaminant sequences that is indistinguishable from previously described endogenous sequences. Furthermore, the nucleotide composition pattern observed in 5' and 3' ends of contaminant sequences is much more complex than the flat pattern previously described in some Neandertal contaminants. Although much research on samples with known contaminant histories is needed, our results suggest that endogenous and contaminant sequences cannot be distinguished by the fragmentation pattern alone.     

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

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

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.     

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.     

Technical note: improved ancient DNA purification for PCR using ion-exchange columns

A novel method of ancient DNA (aDNA) purification was developed using ion-exchange columns to improve PCR-amplifiable DNA extraction from ancient bone samples. Thirteen PCR-resistant ancient bone samples aged 500-3,300 years were tested to extract aDNA using a recently reported, silica-based aDNA extraction method and an ion-exchange column method for the further purification. The PCR success rates of the aDNA extracts were evaluated for the amplification ability of the fragments of mitochondrial DNA, a high-copy DNA, and amelogenin, a low-copy DNA. The results demonstrate that the further purification of silica-based aDNA extracts using ion-exchange columns considerably improved PCR amplification. We suggest that the ion-exchange column-based method will be useful for the improvement of PCR-amplifiable aDNA extraction, particularly from the poorly preserved, PCR-resistant, ancient samples.     

Hypersensitive PCR, ancient human mtDNA, and contamination

When highly efficient polymerase was used with high cycle numbers (50-60), strong amplifications were observed, but negative controls were also unexpectedly amplified in a study of ancient human mtDNA from 2000-year-old skeletons. The results of a series of tests revealed that the hypersensitive polymerase chain reaction (PCR) generated by higher cycles and the presence of contaminant DNA (though at extremely low levels) should be responsible for the amplification of negative controls. We suggest that PCR sensitivity be optimized to take advantage of highly efficient polymerase and at the same time prevent "background DNA" from becoming "contaminant DNA" and obscuring the analysis of authentic ancient DNA. We propose the use of multiple positive controls when amplifying ancient human mtDNA samples to indicate the sensitivity of individual PCR amplifications and to monitor the contamination levels of modern human DNA. This study provides some suggestions as to how to amplify and analyze ancient human mtDNA when unavoidable and extremely tiny amounts of modern human DNA exist.     

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.     

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.     

Novel high-resolution characterization of ancient DNA reveals C > U-type base modification events as the sole cause of post mortem miscoding lesions

Ancient DNA (aDNA) research has long depended on the power of PCR to amplify trace amounts of surviving genetic material from preserved specimens. While PCR permits specific loci to be targeted and amplified, in many ways it can be intrinsically unsuited to damaged and degraded aDNA templates. PCR amplification of aDNA can produce highly-skewed distributions with significant contributions from miscoding lesion damage and non-authentic sequence artefacts. As traditional PCR-based approaches have been unable to fully resolve the molecular nature of aDNA damage over many years, we have developed a novel single primer extension (SPEX)-based approach to generate more accurate sequence information. SPEX targets selected template strands at defined loci and can generate a quantifiable redundancy of coverage; providing new insights into the molecular nature of aDNA damage and fragmentation. SPEX sequence data reveals inherent limitations in both traditional and metagenomic PCR-based approaches to aDNA, which can make current damage analyses and correct genotyping of ancient specimens problematic. In contrast to previous aDNA studies, SPEX provides strong quantitative evidence that C > U-type base modifications are the sole cause of authentic endogenous damage-derived miscoding lesions. This new approach could allow ancient specimens to be genotyped with unprecedented accuracy.     

Ancient DNA sheds new light on the Svalbard foraminiferal fossil record of the last millennium

Recent palaeogenetic studies have demonstrated the occurrence of preserved ancient DNA (aDNA) in various types of fossilised material. Environmental aDNA sequences assigned to modern species have been recovered from marine sediments dating to the Pleistocene. However, the match between the aDNA and the fossil record still needs to be evaluated for the environmental DNA approaches to be fully exploited. Here, we focus on foraminifera in sediments up to one thousand years old retrieved from the Hornsund fjord (Svalbard). We compared the diversity of foraminiferal microfossil assemblages with the diversity of aDNA sequenced from subsurface sediment samples using both cloning and high‐throughput sequencing (HTS). Our study shows that 57% of the species archived in the fossil record were also detected in the aDNA data. However, the relative abundance of aDNA sequence reads and fossil specimens differed considerably. We also found a limited match between the stratigraphic occurrence of some fossil species and their aDNA sequences, especially in the case of rare taxa. The aDNA data comprised a high proportion of non‐fossilised monothalamous species, which are known to dominate in modern foraminiferal communities of the Svalbard region. Our results confirm the relevance of HTS for studying past micro‐eukaryotic diversity and provide insight into its ability to reflect fossil assemblages. Palaeogenetic studies including aDNA analyses of non‐fossilised groups expand the range of palaeoceanographical proxies and therefore may increase the accuracy of palaeoenvironmental reconstructions.     

Brief communication: identification of the authentic ancient DNA sequence in a human bone contaminated with modern DNA

We present a method to distinguish authentic ancient DNA from contaminating DNA in a human bone. This is achieved by taking account of the spatial distribution of the various sequence families within the bone and the extent of degradation of the template DNAs, as revealed by the error content of the sequences. To demonstrate the veracity of the method, we handled two ancient human tibiae in order to contaminate them with modern DNA, and then subjected segments of the bones to various decontaminating treatments, including removal of the outer 1-2 mm, before extracting DNA, cloning, and obtaining a total of 107 mitochondrial DNA sequences. Sequences resulting from the deliberate contamination were located exclusively in the outer 1-2 mm of the bones, and only one of these 27 sequences contained an error that could be ascribed to DNA degradation. A second, much smaller set of relatively error-free sequences, which we ascribe to contamination during excavation or curation, was also located exclusively in the outer 1-2 mm. In contrast, a family of 72 sequences, displaying extensive degradation products but identifiable as haplogroup U5a1a, was distributed throughout one of the bones and represents the authentic ancient DNA content of this specimen.     

Damage and repair of ancient DNA

Under certain conditions small amounts of DNA can survive for long periods of time and can be used as polymerase chain reaction (PCR) substrates for the study of phylogenetic relationships and population genetics of extinct plants and animals, including hominids. Because of extensive DNA degradation, these studies are limited to species that lived within the past 10(4)-10(5) years (Late Pleistocene), although DNA sequences from 10(6) years have been reported. Ancient DNA (aDNA) has been used to study phylogenetic relationships of protists, fungi, algae, plants, and higher eukaryotes such as extinct horses, cave bears, the marsupial wolf, the moa, and Neanderthal. In the past few years, this technology has been extended to the study of infectious disease in ancient Egyptian and South American mummies, the dietary habits of ancient animals, and agricultural practices and population dynamics of early native Americans. Hence, ancient DNA contains information pertinent to numerous fields of study including evolution, population genetics, ecology, climatology, medicine, archeology, and behavior. The major obstacles to the study of aDNA are its extremely low yield, contamination with modern DNA, and extensive degradation. In the course of this review, we will discuss the current aDNA literature describing the importance of aDNA studies as they relate to important biological questions and the difficulties associated with extracting useful information from highly degraded and damaged substrates derived from limited sources. In addition, we will present some of our own preliminary and published data on mechanisms of DNA degradation and some speculative thoughts on strategies for repair and restoration of aDNA.     

Combining bleach and mild predigestion improves ancient DNA recovery from bones

The feasibility of genome‐scale studies from archaeological material remains critically dependent on the ability to access endogenous, authentic DNA. In the majority of cases, this represents a few per cent of the DNA extract, at most. A number of specific pre‐extraction protocols for bone powder aimed to improve ancient DNA recovery before library amplification have recently been developed. Here, we test the effects of combining two of such protocols, a bleach wash and a predigestion step, on 12 bone samples of Atlantic cod and domestic horse aged 750–1350 cal. years before present. Using high‐throughput sequencing, we show that combined together, bleach wash and predigestion consistently yield DNA libraries with higher endogenous content than either of these methods alone. Additionally, the molecular complexity of these libraries is improved and endogenous DNA templates show larger size distributions. Other library characteristics, such as DNA damage profiles or the composition of microbial communities, are little affected by the pre‐extraction protocols. Application of the combined protocol presented in this study will facilitate the genetic analysis of an increasing number of ancient remains and will reduce the cost of whole‐genome sequencing.     

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.     

Blocking human contaminant DNA during PCR allows amplification of rare mammal species from sedimentary ancient DNA

Analyses of degraded DNA are typically hampered by contamination, especially when employing universal primers such as commonly used in environmental DNA studies. In addition to false-positive results, the amplification of contaminant DNA may cause false-negative results because of competition, or bias, during the PCR. In this study, we test the utility of human-specific blocking primers in mammal diversity analyses of ancient permafrost samples from Siberia. Using quantitative PCR (qPCR) on human and mammoth DNA, we first optimized the design and concentration of blocking primer in the PCR. Subsequently, 454 pyrosequencing of ancient permafrost samples amplified with and without the addition of blocking primer revealed that DNA sequences from a diversity of mammalian representatives of the Beringian megafauna were retrieved only when the blocking primer was added to the PCR. Notably, we observe the first retrieval of woolly rhinoceros (Coelodonta antiquitatis) DNA from ancient permafrost cores. In contrast, reactions without blocking primer resulted in complete dominance by human DNA sequences. These results demonstrate that in ancient environmental analyses, the PCR can be biased towards the amplification of contaminant sequences to such an extent that retrieval of the endogenous DNA is severely restricted. The application of blocking primers is a promising tool to avoid this bias and can greatly enhance the quantity and the diversity of the endogenous DNA sequences that are amplified.     

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.     

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.     

Beyond ancient microbial DNA: nonnucleotidic biomolecules for paleomicrobiology

Identifying the causes of past epidemics depends on the specific detection of pathogens in buried individuals; this field of research is known as paleomicrobiology, an emerging field that has benefited from technological advances in microbiology. For almost 15 years, the detection, identification, and characterization of microbes in ancient environmental and human specimens emerged on the basis of ancient DNA (aDNA) analyses. aDNA limitations due to potential contamination by modern DNA and altered aDNA led to the development of alternative methods for the detection and characterization of nonnucleotidic biomolecules, including mycolic acids (of ancient mycobacteria) and proteins. Accordingly, immunohistochemistry, immunochromatography, and enzyme-linked immunosorbent assay techniques have been developed for the specific detection of microbes from ancient human and environmental specimens. Protein analysis by mass spectrometry, a standard for ancient animal identification, has also recently emerged as a technique for ancient mycobacteria detection, while immuno-PCR is yet another promising technique. As with aDNA, strict protocols must be enforced to ensure authenticity of the data. Here we review the analysis of nonnucleotidic biomolecules from ancient microbes and the ability of these analyses to complement aDNA analyses, which opens new opportunities for identification of ancient microbes as well as new avenues to potentially resolve controversies regarding the cause of some historical pandemics and study the coevolution of microbes and hosts.