Present-Day DNA Contamination in Ancient DNA Datasets

Inferring diet, disease and antibiotic resistance from ancient human oral microbiomes

The interaction between a host and its microbiome is an area of intense study. For the human host, it is known that the various body-site-associated microbiomes impact heavily on health and disease states. For instance, the oral microbiome is a source of various pathogens and potential antibiotic resistance gene pools. The effect of historical changes to the human host and environment to the associated microbiome, however, has been less well explored. In this review, we characterize several historical and prehistoric events which are considered to have impacted the oral environment and therefore the bacterial communities residing within it. The link between evolutionary changes to the oral microbiota and the significant societal and behavioural changes occurring during the pre-Neolithic, Agricultural Revolution, Industrial Revolution and Antibiotic Era is outlined. While previous studies suggest the functional profile of these communities may have shifted over the centuries, there is currently a gap in knowledge that needs to be filled. Biomolecular archaeological evidence of innate antimicrobial resistance within the oral microbiome shows an increase in the abundance of antimicrobial resistance genes since the advent and widespread use of antibiotics in the modern era. Nevertheless, a lack of research into the prevalence and evolution of antimicrobial resistance within the oral microbiome throughout history hinders our ability to combat antimicrobial resistance in the modern era.

Assessing the impact of post-mortem damage and contamination on imputation performance in ancient DNA

Low-coverage imputation is becoming ever more present in ancient DNA (aDNA) studies. Imputation pipelines commonly used for present-day genomes have been shown to yield accurate results when applied to ancient genomes. However, post-mortem damage (PMD), in the form of C-to-T substitutions at the reads termini, and contamination with DNA from closely related species can potentially affect imputation performance in aDNA. In this study, we evaluated imputation performance (i) when using a genotype caller designed for aDNA, ATLAS, compared to bcftools, and (ii) when contamination is present. We evaluated imputation performance with principal component analyses and by calculating imputation error rates. With a particular focus on differently imputed sites, we found that using ATLAS prior to imputation substantially improved imputed genotypes for a very damaged ancient genome (42% PMD). Trimming the ends of the sequencing reads led to similar improvements in imputation accuracy. For the remaining genomes, ATLAS brought limited gains. Finally, to examine the effect of contamination on imputation, we added various amounts of reads from two present-day genomes to a previously downsampled high-coverage ancient genome. We observed that imputation accuracy drastically decreased for contamination rates above 5%. In conclusion, we recommend (i) accounting for PMD by either trimming sequencing reads or using a genotype caller such as ATLAS before imputing highly damaged genomes and (ii) only imputing genomes containing up to 5% of contamination.

decOM: similarity-based microbial source tracking of ancient oral samples using k-mer-based methods

BACKGROUND

The analysis of ancient oral metagenomes from archaeological human and animal samples is largely confounded by contaminant DNA sequences from modern and environmental sources. Existing methods for Microbial Source Tracking (MST) estimate the proportions of environmental sources, but do not perform well on ancient metagenomes. We developed a novel method called decOM for Microbial Source Tracking and classification of ancient and modern metagenomic samples using k-mer matrices.

RESULTS

We analysed a collection of 360 ancient oral, modern oral, sediment/soil and skin metagenomes, using stratified five-fold cross-validation. decOM estimates the contributions of these source environments in ancient oral metagenomic samples with high accuracy, outperforming two state-of-the-art methods for source tracking, FEAST and mSourceTracker.

CONCLUSIONS

decOM is a high-accuracy microbial source tracking method, suitable for ancient oral metagenomic data sets. The decOM method is generic and could also be adapted for MST of other ancient and modern types of metagenomes. We anticipate that decOM will be a valuable tool for MST of ancient metagenomic studies. Video Abstract.

Developing the Protocol Infrastructure for DNA Sequencing Natural History Collections

Intentionally preserved biological material in natural history collections represents a vast repository of biodiversity. Advances in laboratory and sequencing technologies have made these specimens increasingly accessible for genomic analyses, offering a window into the genetic past of species and often permitting access to information that can no longer be sampled in the wild. Due to their age, preparation and storage conditions, DNA retrieved from museum and herbarium specimens is often poor in yield, heavily fragmented and biochemically modified. This not only poses methodological challenges in recovering nucleotide sequences, but also makes such investigations susceptible to environmental and laboratory contamination. In this paper, we review the practical challenges associated with making the recovery of DNA sequence data from museum collections more routine. We first review key operational principles and issues to address, to guide the decision-making process and dialogue between researchers and curators about when and how to sample museum specimens for genomic analyses. We then outline the range of steps that can be taken to reduce the likelihood of contamination including laboratory set-ups, workflows and working practices. We finish by presenting a series of case studies, each focusing on protocol practicalities for the application of different mainstream methodologies to museum specimens including: (i) shotgun sequencing of insect mitogenomes, (ii) whole genome sequencing of insects, (iii) genome skimming to recover plant plastid genomes from herbarium specimens, (iv) target capture of multi-locus nuclear sequences from herbarium specimens, (v) RAD-sequencing of bird specimens and (vi) shotgun sequencing of ancient bovid bone samples.

Imputation of ancient human genomes

Due to postmortem DNA degradation and microbial colonization, most ancient genomes have low depth of coverage, hindering genotype calling. Genotype imputation can improve genotyping accuracy for low-coverage genomes. However, it is unknown how accurate ancient DNA imputation is and whether imputation introduces bias to downstream analyses. Here we re-sequence an ancient trio (mother, father, son) and downsample and impute a total of 43 ancient genomes, including 42 high-coverage (above 10x) genomes. We assess imputation accuracy across ancestries, time, depth of coverage, and sequencing technology. We find that ancient and modern DNA imputation accuracies are comparable. When downsampled at 1x, 36 of the 42 genomes are imputed with low error rates (below 5%) while African genomes have higher error rates. We validate imputation and phasing results using the ancient trio data and an orthogonal approach based on Mendel's rules of inheritance. We further compare the downstream analysis results between imputed and high-coverage genomes, notably principal component analysis, genetic clustering, and runs of homozygosity, observing similar results starting from 0.5x coverage, except for the African genomes. These results suggest that, for most populations and depths of coverage as low as 0.5x, imputation is a reliable method that can improve ancient DNA studies.

Artefact Profiling: Panomics Approaches for Understanding the Materiality of Written Artefacts

This review explains the strategies behind genomics, proteomics, metabolomics, metallomics and isotopolomics approaches and their applicability to written artefacts. The respective sub-chapters give an insight into the analytical procedure and the conclusions drawn from such analyses. A distinction is made between information that can be obtained from the materials used in the respective manuscript and meta-information that cannot be obtained from the manuscript itself, but from residues of organisms such as bacteria or the authors and readers. In addition, various sampling techniques are discussed in particular, which pose a special challenge in manuscripts. The focus is on high-resolution, non-targeted strategies that can be used to extract the maximum amount of information about ancient objects. The combination of the various omics disciplines (panomics) especially offers potential added value in terms of the best possible interpretations of the data received. The information obtained can be used to understand the production of ancient artefacts, to gain impressions of former living conditions, to prove their authenticity, to assess whether there is a toxic hazard in handling the manuscripts, and to be able to determine appropriate measures for their conservation and restoration.

SNP analysis of challenging bone DNA samples using the HID-Ion AmpliSeq™ Identity Panel: facts and artefacts

PCR-MPS is an emerging tool for the analysis of low-quality DNA samples. In this study, we used PCR-MPS to analyse 32 challenging bone DNA samples from three Second World War victims, which previously yielded no results in conventional STR PCR-CE typing. The Identity Panel was used with 27 cycles of PCR. Despite that we only had an average of 6.8 pg of degraded DNA as template, 30 out of 32 libraries (93.8%) produced sequencing data for about 63/90 autosomal markers per sample. Out of the 30 libraries, 14 (46.7%) yielded single source genetic profiles in agreement with the biological identity of the donor, whereas 12 cases (40.0%) resulted in SNP profiles that did not match or were mixed. The misleading outcomes for those 12 cases were likely due to hidden exogenous human contamination, as shown by the higher frequencies of allelic imbalance, unusual high frequencies of allelic drop-ins, high heterozygosity levels in the consensus profiles generated from challenging samples, and traces of amplified molecular products in four out of eight extraction negative controls. Even if the source and the time of the contamination were not identified, it is likely that it occurred along the multi-step bone processing workflow. Our results suggest that only positive identification by statistical tools (e.g. likelihood ratio) should be accepted as reliable; oppositely, the results leading to exclusion should be treated as inconclusive because of potential contamination issues. Finally, strategies are discussed for monitoring the workflow of extremely challenging bone samples in PCR-MPS experiments with an increased number of PCR cycles.

KIN: a method to infer relatedness from low-coverage ancient DNA

Genetic kinship of ancient individuals can provide insights into their culture and social hierarchy, and is relevant for downstream genetic analyses. However, estimating relatedness from ancient DNA is difficult due to low-coverage, ascertainment bias, or contamination from various sources. Here, we present KIN, a method to estimate the relatedness of a pair of individuals from the identical-by-descent segments they share. KIN accurately classifies up to 3rd-degree relatives using at least 0.05x sequence coverage and differentiates siblings from parent-child pairs. It incorporates additional models to adjust for contamination and detect inbreeding, which improves classification accuracy.

Mitochondrial DNA Consensus Calling and Quality Filtering for Constructing Ancient Human Mitogenomes: Comparison of Two Widely Applied Methods

Retrieving high-quality endogenous ancient DNA (aDNA) poses several challenges, including low molecular copy number, high rates of fragmentation, damage at read termini, and potential presence of exogenous contaminant DNA. All these factors complicate a reliable reconstruction of consensus aDNA sequences in reads from high-throughput sequencing platforms. Here, we report findings from a thorough evaluation of two alternative tools (ANGSD and schmutzi) aimed at overcoming these issues and constructing high-quality ancient mitogenomes. Raw genomic data (BAM/FASTQ) from a total of 17 previously published whole ancient human genomes ranging from the 14th to the 7th millennium BCE were retrieved and mitochondrial consensus sequences were reconstructed using different quality filters, with their accuracy measured and compared. Moreover, the influence of different sequence parameters (number of reads, sequenced bases, mean coverage, and rate of deamination and contamination) as predictors of derived sequence quality was evaluated. Complete mitogenomes were successfully reconstructed for all ancient samples, and for the majority of them, filtering substantially improved mtDNA consensus calling and haplogroup prediction. Overall, the schmutzi pipeline, which estimates and takes into consideration exogenous contamination, appeared to have the edge over the much faster and user-friendly alternative method (ANGSD) in moderate to high coverage samples (>1,000,000 reads). ANGSD, however, through its read termini trimming filter, showed better capabilities in calling the consensus sequence from low-quality samples. Among all the predictors of overall sample quality examined, the strongest correlation was found for the available number of sequence reads and bases. In the process, we report a previously unassigned haplogroup (U3b) for an Early Chalcolithic individual from Southern Anatolia/Northern Levant.

Ancient Metagenomic Studies: Considerations for the Wider Scientific Community

Like modern metagenomics, ancient metagenomics is a highly data-rich discipline, with the added challenge that the DNA of interest is degraded and, depending on the sample type, in low abundance. This requires the application of specialized measures during molecular experiments and computational analyses. Furthermore, researchers often work with finite sample sizes, which impedes optimal experimental design and control of confounding factors, and with ethically sensitive samples necessitating the consideration of additional guidelines. In September 2020, early career researchers in the field of ancient metagenomics met (Standards, Precautions & Advances in Ancient Metagenomics 2 [SPAAM2] community meeting) to discuss the state of the field and how to address current challenges. Here, in an effort to bridge the gap between ancient and modern metagenomics, we highlight and reflect upon some common misconceptions, provide a brief overview of the challenges in our field, and point toward useful resources for potential reviewers and newcomers to the field.

Inline index helped in cleaning up data contamination generated during library preparation and the subsequent steps

BACKGROUND

High-throughput sequencing involves library preparation and amplification steps, which may induce contamination across samples or between samples and the environment.

METHODS

We tested the effect of applying an inline-index strategy, in which DNA indices of 6 bp were added to both ends of the inserts at the ligation step of library prep for resolving the data contamination problem.

RESULTS

Our results showed that the contamination ranged from 0.29 to 1.25% in one experiment and from 0.83 to 27.01% in the other. We also found that contamination could be environmental or from reagents besides cross-contamination between samples.

CONCLUSIONS

Inline-index method is a useful experimental design to clean up the data and address the contamination problem which has been plaguing high-throughput sequencing data in many applications.

Badania kopalnego DNA – możliwości i ograniczenia

Ostatnie cztery dekady przyniosły znaczący rozwój archeologii molekularnej i badania nad kopalnym DNA (aDNA). Nowatorskie metody uwzględniają szeroki zakres badań, począwszy od sekwencjonowania niewielkich fragmentów mitochondrialnego DNA po wielkoskalowe badania całych populacji, łączące sekwencjonowanie genomów mitochondrialnych, genów podlegających doborowi naturalnemu, jak i całych genomów jądrowych. Postęp, zwłaszcza w dziedzinie technologii sekwencjonowania DNA, umożliwił pozyskanie informacji ze szczątków paleontologicznych i materiału archeologicznego, umożliwiając zbadanie związków filogenetycznych między wymarłymi i współczesnymi gatunkami. Dzięki zastosowaniu technologii sekwencjonowania nowej generacji możliwe stało się poznanie sekwencji DNA nie tylko bezpośrednio ze szczątków ludzkich lub zwierzęcych, ale także z osadów sedymentacyjnych z głębin jezior oraz jaskiń. W artykule przedstawiono możliwości i ograniczenia występujące w badaniach nad kopalnym DNA ludzi, zwierząt czy bakterii z podkreśleniem wkładu polskich badaczy w rozwój tej dziedziny nauki.

Integrating multi-taxon palaeogenomes and sedimentary ancient DNA to study past ecosystem dynamics

Ancient DNA (aDNA) has played a major role in our understanding of the past. Important advances in the sequencing and analysis of aDNA from a range of organisms have enabled a detailed understanding of processes such as past demography, introgression, domestication, adaptation and speciation. However, to date and with the notable exception of microbiomes and sediments, most aDNA studies have focused on single taxa or taxonomic groups, making the study of changes at the community level challenging. This is rather surprising because current sequencing and analytical approaches allow us to obtain and analyse aDNA from multiple source materials. When combined, these data can enable the simultaneous study of multiple taxa through space and time, and could thus provide a more comprehensive understanding of ecosystem-wide changes. It is therefore timely to develop an integrative approach to aDNA studies by combining data from multiple taxa and substrates. In this review, we discuss the various applications, associated challenges and future prospects of such an approach.

Quantitative Human Paleogenetics: What can Ancient DNA Tell us About Complex Trait Evolution?

Genetic association data from national biobanks and large-scale association studies have provided new prospects for understanding the genetic evolution of complex traits and diseases in humans. In turn, genomes from ancient human archaeological remains are now easier than ever to obtain, and provide a direct window into changes in frequencies of trait-associated alleles in the past. This has generated a new wave of studies aiming to analyse the genetic component of traits in historic and prehistoric times using ancient DNA, and to determine whether any such traits were subject to natural selection. In humans, however, issues about the portability and robustness of complex trait inference across different populations are particularly concerning when predictions are extended to individuals that died thousands of years ago, and for which little, if any, phenotypic validation is possible. In this review, we discuss the advantages of incorporating ancient genomes into studies of trait-associated variants, the need for models that can better accommodate ancient genomes into quantitative genetic frameworks, and the existing limits to inferences about complex trait evolution, particularly with respect to past populations.

Successful application of ancient DNA extraction and library construction protocols to museum wet collection specimens

Millions of scientific specimens are housed in museum collections, a large part of which are fluid preserved. The use of formaldehyde as fixative and subsequent storage in ethanol is especially common in ichthyology and herpetology. This type of preservation damages DNA and reduces the chance of successful retrieval of genetic data. We applied ancient DNA extraction and single stranded library construction protocols to a variety of vertebrate samples obtained from wet collections and of different ages. Our results show that almost all samples tested yielded endogenous DNA. Archival DNA extraction was successful across different tissue types as well as using small amounts of tissue. Conversion of archival DNA fragments into single-stranded libraries resulted in usable data even for samples with initially undetectable DNA amounts. Subsequent target capture approaches for mitochondrial DNA using homemade baits on a subset of 30 samples resulted in almost complete mitochondrial genome sequences in several instances. Thus, application of ancient DNA methodology makes wet collection specimens, including type material as well as rare, old or extinct species, accessible for genetic and genomic analyses. Our results, accompanied by detailed step-by-step protocols, are a large step forward to open the DNA archive of museum wet collections for scientific studies.

HuConTest: Testing Human Contamination in Great Ape Samples

Modern human contamination is a common problem in ancient DNA studies. We provide evidence that this issue is also present in studies in great apes, which are our closest living relatives, for example in noninvasive samples. Here, we present a simple method to detect human contamination in short-read sequencing data from different species: HuConTest. We demonstrate its feasibility using blood and tissue samples from these species. This test is particularly useful for more complex samples (such as museum and noninvasive samples) which have smaller amounts of endogenous DNA, as we show here.

Ancient and modern stickleback genomes reveal the demographic constraints on adaptation

Adaptation is typically studied by comparing modern populations with contrasting environments. Individuals persisting in the ancestral habitat are typically used to represent the ancestral founding population; however, it has been questioned whether these individuals are good proxies for the actual ancestors.¹ To address this, we applied a paleogenomics approach² to directly access the ancestral genepool: partially sequencing the genomes of two 11- to 13,000-year-old stickleback recovered from the transitionary layer between marine and freshwater sediments of two Norwegian isolation lakes³ and comparing them with 30 modern stickleback genomes from the same lakes and adjacent marine fjord, in addition to a global dataset of 20 genomes.⁴ The ancient stickleback shared genome-wide ancestry with the modern fjord population, whereas modern lake populations have lost substantial ancestral variation following founder effects, and subsequent drift and selection. Freshwater-adaptive alleles found in one ancient stickleback genome have not risen to high frequency in the present-day population from the same lake. Comparison to the global dataset suggested incomplete adaptation to freshwater in our modern lake populations. Our findings reveal the impact of population bottlenecks in constraining adaptation due to reduced efficacy of selection on standing variation present in founder populations.

Community-curated and standardised metadata of published ancient metagenomic samples with AncientMetagenomeDir

Ancient DNA and RNA are valuable data sources for a wide range of disciplines. Within the field of ancient metagenomics, the number of published genetic datasets has risen dramatically in recent years, and tracking this data for reuse is particularly important for large-scale ecological and evolutionary studies of individual taxa and communities of both microbes and eukaryotes. AncientMetagenomeDir (archived at 10.5281/zenodo.3980833) is a collection of annotated metagenomic sample lists derived from published studies that provide basic, standardised metadata and accession numbers to allow rapid data retrieval from online repositories. These tables are community-curated and span multiple sub-disciplines to ensure adequate breadth and consensus in metadata definitions, as well as longevity of the database. Internal guidelines and automated checks facilitate compatibility with established sequence-read archives and term-ontologies, and ensure consistency and interoperability for future meta-analyses. This collection will also assist in standardising metadata reporting for future ancient metagenomic studies.

Measurement(s)	genome • Metagenome • Metadata • Ancient DNA
Technology Type(s)	digital curation
Factor Type(s)	geographic location • sample age

Machine-accessible metadata file describing the reported data: 10.6084/m9.figshare.13241537

Competitive mapping allows for the identification and exclusion of human DNA contamination in ancient faunal genomic datasets

BACKGROUND

After over a decade of developments in field collection, laboratory methods and advances in high-throughput sequencing, contamination remains a key issue in ancient DNA research. Currently, human and microbial contaminant DNA still impose challenges on cost-effective sequencing and accurate interpretation of ancient DNA data.

RESULTS

Here we investigate whether human contaminating DNA can be found in ancient faunal sequencing datasets. We identify variable levels of human contamination, which persists even after the sequence reads have been mapped to the faunal reference genomes. This contamination has the potential to affect a range of downstream analyses.

CONCLUSIONS

We propose a fast and simple method, based on competitive mapping, which allows identifying and removing human contamination from ancient faunal DNA datasets with limited losses of true ancient data. This method could represent an important tool for the ancient DNA field.