The Neandertal genome and ancient DNA authenticity

Evaluation of the usefulness of insertion-null markers in critical skeletal remains

Forensic DNA analysis in compromised skeletal remains may pose challenges due to DNA degradation, often resulting in partial or negative autosomal STRs profiles. To address this issue, alternative approaches such as mitochondrial DNA or SNPs typing may be employed; however, they are labour-intensive and costly. Insertion-null alleles (INNULs), short interspersed nuclear elements, have been suggested as a valuable tool for human identification in challenging samples due to their small amplicon size. A commercial kit including 20 INNULs markers along with amelogenin (InnoTyper® 21) has been developed. This study assesses its utility using degraded skeletal remains, comparing the results obtained (the number of detected alleles, RFU values, PHR, and the number of reportable markers) to those obtained using GlobalFiler™. Subsequently, the random match probability of the two profiles for each sample was determined using Familias version 3 to evaluate the power of discrimination of the results obtained from each kit. In every sample, InnoTyper® 21 yielded more alleles, higher RFU values, and a greater number of reportable loci. However, in most cases, both profiles were similarly informative. In conclusion, InnoTyper® 21 serves as a valuable complement to the analysis of challenging samples in cases where a poor or negative profile was obtained.

Performance comparison of four qPCR and three autosomal STR commercial kits from degraded skeletal remains

This research evaluates the current DNA quantification (Quantifiler™ Trio, PowerQuant®, Investigator® Quantiplex® Pro and InnoQuant® HY Fast) and autosomal STRs amplification kits (GlobalFiler™, PowerPlex® Fusion 6 C, Investigator® 24Plex QS) using 62 degraded skeletal remains from armed conflicts (petrous bone, femur, tibia, and tooth) with several parameters (autosomal small, large, and male target, degradation index, probability of degradation, number of alleles above analytical threshold, number of alleles above stochastic threshold, RFU, peak height ratio, number of reportable loci). The best qPCR/autosomal STRs amplification tandem was determined by comparing quantification results by a DNA quantity estimation based on sample average RFU. InnoQuant® HY Fast was the most sensitive kit, and no significative differences were observed among amplification kits; however, Investigator® 24 Plex QS was found to be the most sensitive in our samples. That is why InnoQuant™ and Investigator® 24Plex QS were determined to be the best tandem.

Review: Computational analysis of human skeletal remains in ancient DNA and forensic genetics

Degraded DNA is used to answer questions in the fields of ancient DNA (aDNA) and forensic genetics. While aDNA studies typically center around human evolution and past history, and forensic genetics is often more concerned with identifying a specific individual, scientists in both fields face similar challenges. The overlap in source material has prompted periodic discussions and studies on the advantages of collaboration between fields toward mutually beneficial methodological advancements. However, most have been centered around wet laboratory methods (sampling, DNA extraction, library preparation, etc.). In this review, we focus on the computational side of the analytical workflow. We discuss limitations and considerations to consider when working with degraded DNA. We hope this review provides a framework to researchers new to computational workflows for how to think about analyzing highly degraded DNA and prompts an increase of collaboration between the forensic genetics and aDNA fields.

Demographic reconstruction of the Western sheep expansion from whole-genome sequences

As one of the earliest livestock, sheep (Ovis aries) were domesticated in the Fertile Crescent about 12,000-10,000 years ago and have a nearly worldwide distribution today. Most of our knowledge about the timing of their expansions stems from archaeological data but it is unclear how the genetic diversity of modern sheep fits with these dates. We used whole-genome sequencing data of 63 domestic breeds and their wild relatives, the Asiatic mouflon (O. gmelini, previously known as O. orientalis), to explore the demographic history of sheep. On the global scale, our analysis revealed geographic structuring among breeds with unidirectional recent gene flow from domestics into Asiatic mouflons. We then selected 4 representative breeds from Spain, Morocco, the United Kingdom, and Iran to build a comprehensive demographic model of the Western sheep expansion. We inferred a single domestication event around 11,000 years ago. The subsequent westward expansion is dated to approximately 7,000 years ago, later than the original Neolithic expansion of sheep and slightly predating the Secondary Product Revolution associated with wooly sheep. We see some signals of recent gene flow from an ancestral population into Southern European breeds which could reflect admixture with feral European mouflon. Furthermore, our results indicate that many breeds experienced a reduction of their effective population size during the last centuries, probably associated with modern breed development. Our study provides insights into the complex demographic history of Western Eurasian sheep, highlighting interactions between breeds and their wild counterparts.

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.

A comparison between petrous bone and tooth, femur and tibia DNA analysis from degraded skeletal remains

Skeletal remains are the only biological material that remains after long periods; however, environmental conditions such as temperature, humidity, and pH affect DNA preservation, turning skeletal remains into a challenging sample for DNA laboratories. Sample selection is a key factor, and femur and tooth have been traditionally recommended as the best substrate of genetic material. Recently, petrous bone (cochlear area) has been suggested as a better option due to its DNA yield. This research aims to evaluate the efficiency of petrous bone compared to other cranium samples (tooth) and postcranial long bones (femur and tibia). A total amount of 88 samples were selected from 38 different individuals. The samples were extracted by using an organic extraction protocol, DNA quantification by Quantifiler Trio kit and amplified with GlobalFiler kit. Results show that petrous bone outperforms other bone remains in quantification data, yielding 15-30 times more DNA than the others. DNA profile data presented likeness between petrous bone and tooth regarding detected alleles; however, the amount of DNA extracted in petrous bones allowed us to obtain more informative DNA profiles with superior quality. In conclusion, petrous bone or teeth sampling is recommended if DNA typing is going to be performed with environmentally degraded skeletal remains.

Density separation of petrous bone powders for optimized ancient DNA yields

Density separation is a process routinely used to segregate minerals, organic matter, and even microplastics, from soils and sediments. Here we apply density separation to archaeological bone powders before DNA extraction to increase endogenous DNA recovery relative to a standard control extraction of the same powders. Using nontoxic heavy liquid solutions, we separated powders from the petrous bones of 10 individuals of similar archaeological preservation into eight density intervals (2.15 to 2.45 g/cm3, in 0.05 increments). We found that the 2.30 to 2.35 g/cm3 and 2.35 to 2.40 g/cm3 intervals yielded up to 5.28-fold more endogenous unique DNA than the corresponding standard extraction (and up to 8.53-fold before duplicate read removal), while maintaining signals of ancient DNA authenticity and not reducing library complexity. Although small 0.05 g/cm3 intervals may maximally optimize yields, a single separation to remove materials with a density above 2.40 g/cm3 yielded up to 2.57-fold more endogenous DNA on average, which enables the simultaneous separation of samples that vary in preservation or in the type of material analyzed. While requiring no new ancient DNA laboratory equipment and fewer than 30 min of extra laboratory work, the implementation of density separation before DNA extraction can substantially boost endogenous DNA yields without decreasing library complexity. Although subsequent studies are required, we present theoretical and practical foundations that may prove useful when applied to other ancient DNA substrates such as teeth, other bones, and sediments.

Advancements and Challenges in Ancient DNA Research: Bridging the Global North-South Divide

Ancient DNA (aDNA) research first began in 1984 and ever since has greatly expanded our understanding of evolution and migration. Today, aDNA analysis is used to solve various puzzles about the origin of mankind, migration patterns, and the spread of infectious diseases. The incredible findings ranging from identifying the new branches within the human family to studying the genomes of extinct flora and fauna have caught the world by surprise in recent times. However, a closer look at these published results points out a clear Global North and Global South divide. Therefore, through this research, we aim to emphasize encouraging better collaborative opportunities and technology transfer to support researchers in the Global South. Further, the present research also focuses on expanding the scope of the ongoing conversation in the field of aDNA by reporting relevant literature published around the world and discussing the advancements and challenges in the field.

A mobile lab for ancient DNA extraction in Perug

We report the use of a mobile laboratory set up to extract ancient DNA (aDNA) from 34 human coprolites (fossilized faeces) samples. Our approach enabled the rapid genetic characterization of 5,000 years old archeological samples. It is useful for the on-site screening of museums and freshly excavated samples for DNA. This approach is accessible to other investigators as the mobile laboratory was set up using commercially available instruments.

Fine-scale population genetic structure and barriers to gene flow in a widespread seabird (Ardenna pacifica)

Pelagic seabirds are highly mobile, reducing opportunities for population isolation that might promote differentiation and speciation. At the same time, many seabirds are philopatric, and their tendency to return to their natal islands to breed might reduce gene flow sufficiently to permit local adaptation and differentiation. To test the net impact of these competing processes, estimates of differentiation and gene flow based on comprehensive geographical sampling are required. We leveraged diverse source material to achieve comprehensive geographical sampling in a widespread seabird, the Wedge-tailed Shearwater (Ardenna pacifica). Using data from sequence capture and high-throughput sequencing of 2402 loci containing 20 780 single nucleotide polymorphisms, we tested for population differentiation and gene flow among breeding areas. We found little evidence of deep divergences within A. pacifica but were able to resolve fine-scale differentiation across island groups. This differentiation was sufficient to assign individuals sampled away from breeding areas to their likely source populations. Estimated effective migration surfaces revealed reduced migration between the Indian Ocean and Pacific Ocean, presumably owing to land barriers, and across the equatorial Pacific Ocean, perhaps associated with differences in breeding schedule. Our results reveal that, despite their mobility, pelagic seabirds can exhibit fine-scale population differentiation and reduced gene flow among ocean basins.

Bioarchaeological evidence of one of the earliest Islamic burials in the Levant

The Middle East plays a central role in human history harbouring a vast diversity of ethnic, cultural and religious groups. However, much remains to be understood about past and present genomic diversity in this region. Here we present a multidisciplinary bioarchaeological analysis of two individuals dated to the late 7th and early 8th centuries, the Umayyad Era, from Tell Qarassa, an open-air site in modern-day Syria. Radiocarbon dates and burial type are consistent with one of the earliest Islamic Arab burials in the Levant. Interestingly, we found genomic similarity to a genotyped group of modern-day Bedouins and Saudi rather than to most neighbouring Levantine groups. This study represents the genomic analysis of a secondary use site with characteristics consistent with an early Islamic burial in the Levant. We discuss our findings and possible historic scenarios in the light of forces such as genetic drift and their possible interaction with religious and cultural processes (including diet and subsistence practices).

Ancient genomic and archaeological data combine to identify a surprisingly early Islamic burial in modern day Syria.

A guide to avian museomics: Insights gained from resequencing hundreds of avian study skins

Biological specimens in natural history collections constitute a massive repository of genetic information. Many specimens have been collected in areas in which they no longer exist or in areas where present‐day collecting is not possible. There are also specimens in collections representing populations or species that have gone extinct. Furthermore, species or populations may have been sampled throughout an extensive time period, which is particularly valuable for studies of genetic change through time. With the advent of high‐throughput sequencing, natural history museum resources have become accessible for genomic research. Consequently, these unique resources are increasingly being used across many fields of natural history. In this paper, we summarize our experiences of resequencing hundreds of genomes from historical avian museum specimens. We publish the protocols we have used and discuss the entire workflow from sampling and laboratory procedures, to the bioinformatic processing of historical specimen data.

The high diversity of SARS-CoV-2-related coronaviruses in pangolins alerts potential ecological risks

Understanding the zoonotic origin and evolution history of SARS-CoV-2 will provide critical insights for alerting and preventing future outbreaks. A significant gap remains for the possible role of pangolins as a reservoir of SARS-CoV-2 related coronaviruses (SC2r-CoVs). Here, we screened SC2r-CoVs in 172 samples from 163 pangolin individuals of four species, and detected positive signals in muscles of four Manis javanica and, for the first time, one M. pentadactyla. Phylogeographic analysis of pangolin mitochondrial DNA traced their origins from Southeast Asia. Using in-solution hybridization capture sequencing, we assembled a partial pangolin SC2r-CoV (pangolin-CoV) genome sequence of 22 895 bp (MP20) from the M. pentadactyla sample. Phylogenetic analyses revealed MP20 was very closely related to pangolin-CoVs that were identified in M. javanica seized by Guangxi Customs. A genetic contribution of bat coronavirus to pangolin-CoVs via recombination was indicated. Our analysis revealed that the genetic diversity of pangolin-CoVs is substantially higher than previously anticipated. Given the potential infectivity of pangolin-CoVs, the high genetic diversity of pangolin-CoVs alerts the ecological risk of zoonotic evolution and transmission of pathogenic SC2r-CoVs.

Genome-wide autosomal, mtDNA, and Y chromosome analysis of King Bela III of the Hungarian Arpad dynasty

The ancient Hungarians, "Madzsars", established their control of the Carpathian Basin in the late ninth century and founded the Hungarian Kingdom around 1000AD. The origin of the Magyars as a tribal federation has been much debated in the past. From the time of the conquest to the early fourteenth century they were ruled by descendants of the Arpad family. In order to learn more about the genetic origin of this family, we here analyzed the genome of Bela III one of the most prominent members of the early Hungarian dynasty that ruled the Hungarian Kingdom from 1172 to 1196. The Y-Chromosome of Bela III belongs to haplogroup R1a-Z2123 that is today found in highest frequency in Central Asia, supporting a Central Asian origin for the ruling lineage of the Hungarian kingdom. The autosomal DNA profile of Bela III, however, falls within the genetic variation of present-day east European populations. This is further supported through his mtDNA genome that belongs to haplogroup H, the most common European maternal lineage, but also found in Central Asia. However, we didn't find an exact haplotype match for Bela III. The typical autosomal and maternal Central Eastern European ancestry among Bela III autosomes might be best explained by consecutive intermarriage with local European ruling families.

Inferring archaic introgression from hominin genetic data

Questions surrounding the timing, extent, and evolutionary consequences of archaic admixture into human populations have a long history in evolutionary anthropology. More recently, advances in human genetics, particularly in the field of ancient DNA, have shed new light on the question of whether or not Homo sapiens interbred with other hominin groups. By the late 1990s, published genetic work had largely concluded that archaic groups made no lasting genetic contribution to modern humans; less than a decade later, this conclusion was reversed following the successful DNA sequencing of an ancient Neanderthal. This reversal of consensus is noteworthy, but the reasoning behind it is not widely understood across all academic communities. There remains a communication gap between population geneticists and paleoanthropologists. In this review, we endeavor to bridge this gap by outlining how technological advancements, new statistical methods, and notable controversies ultimately led to the current consensus.

Heterogeneous Hunter-Gatherer and Steppe-Related Ancestries in Late Neolithic and Bell Beaker Genomes from Present-Day France

The transition from the Late Neolithic to the Bronze Age has witnessed important population and societal changes in western Europe.¹ These include massive genomic contributions of pastoralist herders originating from the Pontic-Caspian steppes^2,3 into local populations, resulting from complex interactions between collapsing hunter-gatherers and expanding farmers of Anatolian ancestry.^4-8 This transition is documented through extensive ancient genomic data from present-day Britain,^9,10 Ireland,^11,12 Iberia,¹³ Mediterranean islands,^14,15 and Germany.⁸ It remains, however, largely overlooked in France, where most focus has been on the Middle Neolithic (n = 63),^8,9,16 with the exception of one Late Neolithic genome sequenced at 0.05× coverage.¹⁶ This leaves the key transitional period covering ∼3,400-2,700 cal. years (calibrated years) BCE genetically unsampled and thus the exact time frame of hunter-gatherer persistence and arrival of steppe migrations unknown. To remediate this, we sequenced 24 ancient human genomes from France spanning ∼3,400-1,600 cal. years BCE. This reveals Late Neolithic populations that are genetically diverse and include individuals with dark skin, hair, and eyes. We detect heterogeneous hunter-gatherer ancestries within Late Neolithic communities, reaching up to ∼63.3% in some individuals, and variable genetic contributions of steppe herders in Bell Beaker populations. We provide an estimate as late as ∼3,800 years BCE for the admixture between Neolithic and Mesolithic populations and as early as ∼2,650 years BCE for the arrival of steppe-related ancestry. The genomic heterogeneity characterized underlines the complex history of human interactions even at the local scale.

An efficient pipeline for ancient DNA mapping and recovery of endogenous ancient DNA from whole-genome sequencing data

Ancient DNA research has developed rapidly over the past few decades due to improvements in PCR and next-generation sequencing (NGS) technologies, but challenges still exist. One major challenge in relation to ancient DNA research is to recover genuine endogenous ancient DNA sequences from raw sequencing data. This is often difficult due to degradation of ancient DNA and high levels of contamination, especially homologous contamination that has extremely similar genetic background with that of the real ancient DNA. In this study, we collected whole-genome sequencing (WGS) data from 6 ancient samples to compare different mapping algorithms. To further explore more effective methods to separate endogenous DNA from homologous contaminations, we attempted to recover reads based on ancient DNA specific characteristics of deamination, depurination, and DNA fragmentation with different parameters. We propose a quick and improved pipeline for separating endogenous ancient DNA while simultaneously decreasing homologous contaminations to very low proportions. Our goal in this research was to develop useful recommendations for ancient DNA mapping and for separation of endogenous DNA to facilitate future studies of ancient DNA.

AuthentiCT: a model of ancient DNA damage to estimate the proportion of present-day DNA contamination

Contamination from present-day DNA is a fundamental issue when studying ancient DNA from historical or archaeological material, and quantifying the amount of contamination is essential for downstream analyses. We present AuthentiCT, a command-line tool to estimate the proportion of present-day DNA contamination in ancient DNA datasets generated from single-stranded DNA libraries. The prediction is based solely on the patterns of post-mortem damage observed on ancient DNA sequences. The method has the power to quantify contamination from as few as 10,000 mapped sequences, making it particularly useful for analysing specimens that are poorly preserved or for which little data is available.

Present-Day DNA Contamination in Ancient DNA Datasets

Present-day contamination can lead to false conclusions in ancient DNA studies. A number of methods are available to estimate contamination, which use a variety of signals and are appropriate for different types of data. Here an overview of currently available methods highlighting their strengths and weaknesses is provided, and a classification based on the signals used to estimate contamination is proposed. This overview aims at enabling researchers to choose the most appropriate methods for their dataset. Based on this classification, potential avenues for the further development of methods are discussed.

The past, present and future of ancient bacterial DNA

Groundbreaking studies conducted in the mid-1980s demonstrated the possibility of sequencing ancient DNA (aDNA), which has allowed us to answer fundamental questions about the human past. Microbiologists were thus given a powerful tool to glimpse directly into inscrutable bacterial history, hitherto inaccessible due to a poor fossil record. Initially plagued by concerns regarding contamination, the field has grown alongside technical progress, with the advent of high-throughput sequencing being a breakthrough in sequence output and authentication. Albeit burdened with challenges unique to the analysis of bacteria, a growing number of viable sources for aDNA has opened multiple avenues of microbial research. Ancient pathogens have been extracted from bones, dental pulp, mummies and historical medical specimens and have answered focal historical questions such as identifying the aetiological agent of the black death as Yersinia pestis. Furthermore, ancient human microbiomes from fossilized faeces, mummies and dental plaque have shown shifts in human commensals through the Neolithic demographic transition and industrial revolution, whereas environmental isolates stemming from permafrost samples have revealed signs of ancient antimicrobial resistance. Culminating in an ever-growing repertoire of ancient genomes, the quickly expanding body of bacterial aDNA studies has also enabled comparisons of ancient genomes to their extant counterparts, illuminating the evolutionary history of bacteria. In this review we summarize the present avenues of research and contextualize them in the past of the field whilst also pointing towards questions still to be answered.

Analysis of oral microbiome from fossil human remains revealed the significant differences in virulence factors of modern and ancient Tannerella forsythia

BACKGROUND

Recent advances in the next-generation sequencing (NGS) allowed the metagenomic analyses of DNA from many different environments and sources, including thousands of years old skeletal remains. It has been shown that most of the DNA extracted from ancient samples is microbial. There are several reports demonstrating that the considerable fraction of extracted DNA belonged to the bacteria accompanying the studied individuals before their death.

RESULTS

In this study we scanned 344 microbiomes from 1000- and 2000- year-old human teeth. The datasets originated from our previous studies on human ancient DNA (aDNA) and on microbial DNA accompanying human remains. We previously noticed that in many samples infection-related species have been identified, among them Tannerella forsythia, one of the most prevalent oral human pathogens. Samples containing sufficient amount of T. forsythia aDNA for a complete genome assembly were selected for thorough analyses. We confirmed that the T. forsythia-containing samples have higher amounts of the periodontitis-associated species than the control samples. Despites, other pathogens-derived aDNA was found in the tested samples it was too fragmented and damaged to allow any reasonable reconstruction of these bacteria genomes. The anthropological examination of ancient skulls from which the T. forsythia-containing samples were obtained revealed the pathogenic alveolar bone loss in tooth areas characteristic for advanced periodontitis. Finally, we analyzed the genetic material of ancient T. forsythia strains. As a result, we assembled four ancient T. forsythia genomes - one 2000- and three 1000- year-old. Their comparison with contemporary T. forsythia genomes revealed a lower genetic diversity within the four ancient strains than within contemporary strains. We also investigated the genes of T. forsythia virulence factors and found that several of them (KLIKK protease and bspA genes) differ significantly between ancient and modern bacteria.

CONCLUSIONS

In summary, we showed that NGS screening of the ancient human microbiome is a valid approach for the identification of disease-associated microbes. Following this protocol, we provided a new set of information on the emergence, evolution and virulence factors of T. forsythia, the member of the oral dysbiotic microbiome.

Human origins in a southern African palaeo-wetland and first migrations

Anatomically modern humans originated in Africa around 200 thousand years ago (ka)^1-4. Although some of the oldest skeletal remains suggest an eastern African origin², southern Africa is home to contemporary populations that represent the earliest branch of human genetic phylogeny^5,6. Here we generate, to our knowledge, the largest resource for the poorly represented and deepest-rooting maternal L0 mitochondrial DNA branch (198 new mitogenomes for a total of 1,217 mitogenomes) from contemporary southern Africans and show the geographical isolation of L0d1'2, L0k and L0g KhoeSan descendants south of the Zambezi river in Africa. By establishing mitogenomic timelines, frequencies and dispersals, we show that the L0 lineage emerged within the residual Makgadikgadi-Okavango palaeo-wetland of southern Africa⁷, approximately 200 ka (95% confidence interval, 240-165 ka). Genetic divergence points to a sustained 70,000-year-long existence of the L0 lineage before an out-of-homeland northeast-southwest dispersal between 130 and 110 ka. Palaeo-climate proxy and model data suggest that increased humidity opened green corridors, first to the northeast then to the southwest. Subsequent drying of the homeland corresponds to a sustained effective population size (L0k), whereas wet-dry cycles and probable adaptation to marine foraging allowed the southwestern migrants to achieve population growth (L0d1'2), as supported by extensive south-coastal archaeological evidence^8-10. Taken together, we propose a southern African origin of anatomically modern humans with sustained homeland occupation before the first migrations of people that appear to have been driven by regional climate changes.

Ancient pathogen genomics as an emerging tool for infectious disease research

Over the past decade, a genomics revolution, made possible through the development of high-throughput sequencing, has triggered considerable progress in the study of ancient DNA, enabling complete genomes of past organisms to be reconstructed. A newly established branch of this field, ancient pathogen genomics, affords an in-depth view of microbial evolution by providing a molecular fossil record for a number of human-associated pathogens. Recent accomplishments include the confident identification of causative agents from past pandemics, the discovery of microbial lineages that are now extinct, the extrapolation of past emergence events on a chronological scale and the characterization of long-term evolutionary history of microorganisms that remain relevant to public health today. In this Review, we discuss methodological advancements, persistent challenges and novel revelations gained through the study of ancient pathogen genomes.

Authentication and Assessment of Contamination in Ancient DNA

Contamination from both present-day humans and postmortem microbial sources is a common challenge in ancient DNA studies. Here we present a suite of tools to assist in the assessment of contamination in ancient DNA data sets. These tools perform standard tests of authenticity of ancient DNA data including detecting the presence of postmortem damage signatures in sequence alignments and quantifying the amount of present-day human contamination.

ARIADNA: machine learning method for ancient DNA variant discovery

Ancient DNA (aDNA) studies often rely on standard methods of mutation calling, optimized for high-quality contemporary DNA but not for excessive contamination, time- or environment-related damage of aDNA. In the absence of validated datasets and despite showing extreme sensitivity to aDNA quality, these methods have been used in many published studies, sometimes with additions of arbitrary filters or modifications, designed to overcome aDNA degradation and contamination problems. The general lack of best practices for aDNA mutation calling may lead to inaccurate results. To address these problems, we present ARIADNA (ARtificial Intelligence for Ancient DNA), a novel approach based on machine learning techniques, using specific aDNA characteristics as features to yield improved mutation calls. In our comparisons of variant callers across several ancient genomes, ARIADNA consistently detected higher-quality genome variants with fast runtimes, while reducing the false positive rate compared with other approaches.

Quantifying and reducing spurious alignments for the analysis of ultra-short ancient DNA sequences

BACKGROUND

The study of ancient DNA is hampered by degradation, resulting in short DNA fragments. Advances in laboratory methods have made it possible to retrieve short DNA fragments, thereby improving access to DNA preserved in highly degraded, ancient material. However, such material contains large amounts of microbial contamination in addition to DNA fragments from the ancient organism. The resulting mixture of sequences constitutes a challenge for computational analysis, since microbial sequences are hard to distinguish from the ancient sequences of interest, especially when they are short.

RESULTS

Here, we develop a method to quantify spurious alignments based on the presence or absence of rare variants. We find that spurious alignments are enriched for mismatches and insertion/deletion differences and lack substitution patterns typical of ancient DNA. The impact of spurious alignments can be reduced by filtering on these features and by imposing a sample-specific minimum length cutoff. We apply this approach to sequences from four ~ 430,000-year-old Sima de los Huesos hominin remains, which contain particularly short DNA fragments, and increase the amount of usable sequence data by 17-150%. This allows us to place a third specimen from the site on the Neandertal lineage.

CONCLUSIONS

Our method maximizes the sequence data amenable to genetic analysis from highly degraded ancient material and avoids pitfalls that are associated with the analysis of ultra-short DNA sequences.

Ratio of mitochondrial to nuclear DNA affects contamination estimates in ancient DNA analysis

In the last decade, ancient DNA research has grown rapidly and started to overcome several of its earlier limitations through Next-Generation-Sequencing (NGS). Among other advances, NGS allows direct estimation of sample contamination from modern DNA sources. First NGS-based approaches of estimating contamination measured heterozygosity. These measurements, however, could only be performed on haploid genomic regions, i.e. the mitochondrial genome or male X chromosomes, but provided no measures of contamination in the nuclear genome of females with their two X chromosomes. Instead, female nuclear contamination is routinely extrapolated from mitochondrial contamination estimates, but it remains unclear if this extrapolation is reliable and to what degree variation in mitochondrial to nuclear DNA ratios affects this extrapolation. We therefore analyzed ancient DNA from 317 samples of different skeletal elements from multiple sites, spanning a temporal range from 7,000 BP to 386 AD. We found that the mitochondrial to nuclear DNA (mt/nc) ratio negatively correlates with an increase in endogenous DNA content and strongly influenced mitochondrial and nuclear contamination estimates in males. The ratio of mt to nc contamination estimates remained stable for overall mt/nc ratios below 200, as found particularly often in petrous bones but less in other skeletal elements and became more variable above that ratio.

Metagenomic Composition Analysis of an Ancient Sequenced Polar Bear Jawbone from Svalbard

The sequencing of ancient DNA samples provides a novel way to find, characterize, and distinguish exogenous genomes of endogenous targets. After sequencing, computational composition analysis enables filtering of undesired sources in the focal organism, with the purpose of improving the quality of assemblies and subsequent data analysis. More importantly, such analysis allows extinct and extant species to be identified without requiring a specific or new sequencing run. However, the identification of exogenous organisms is a complex task, given the nature and degradation of the samples, and the evident necessity of using efficient computational tools, which rely on algorithms that are both fast and highly sensitive. In this work, we relied on a fast and highly sensitive tool, FALCON-meta, which measures similarity against whole-genome reference databases, to analyse the metagenomic composition of an ancient polar bear (Ursus maritimus) jawbone fossil. The fossil was collected in Svalbard, Norway, and has an estimated age of 110,000 to 130,000 years. The FASTQ samples contained 349 GB of nonamplified shotgun sequencing data. We identified and localized, relative to the FASTQ samples, the genomes with significant similarities to reference microbial genomes, including those of viruses, bacteria, and archaea, and to fungal, mitochondrial, and plastidial sequences. Among other striking features, we found significant similarities between modern-human, some bacterial and viral sequences (contamination) and the organelle sequences of wild carrot and tomato relative to the whole samples. For each exogenous candidate, we ran a damage pattern analysis, which in addition to revealing shallow levels of damage in the plant candidates, identified the source as contamination.

Testing of Alignment Parameters for Ancient Samples: Evaluating and Optimizing Mapping Parameters for Ancient Samples Using the TAPAS Tool

High-throughput sequence data retrieved from ancient or other degraded samples has led to unprecedented insights into the evolutionary history of many species, but the analysis of such sequences also poses specific computational challenges. The most commonly used approach involves mapping sequence reads to a reference genome. However, this process becomes increasingly challenging with an elevated genetic distance between target and reference or with the presence of contaminant sequences with high sequence similarity to the target species. The evaluation and testing of mapping efficiency and stringency are thus paramount for the reliable identification and analysis of ancient sequences. In this paper, we present 'TAPAS', (Testing of Alignment Parameters for Ancient Samples), a computational tool that enables the systematic testing of mapping tools for ancient data by simulating sequence data reflecting the properties of an ancient dataset and performing test runs using the mapping software and parameter settings of interest. We showcase TAPAS by using it to assess and improve mapping strategy for a degraded sample from a banded linsang (Prionodon linsang), for which no closely related reference is currently available. This enables a 1.8-fold increase of the number of mapped reads without sacrificing mapping specificity. The increase of mapped reads effectively reduces the need for additional sequencing, thus making more economical use of time, resources, and sample material.

Evolutionary Patterns and Processes: Lessons from Ancient DNA

Ever since its emergence in 1984, the field of ancient DNA has struggled to overcome the challenges related to the decay of DNA molecules in the fossil record. With the recent development of high-throughput DNA sequencing technologies and molecular techniques tailored to ultra-damaged templates, it has now come of age, merging together approaches in phylogenomics, population genomics, epigenomics, and metagenomics. Leveraging on complete temporal sample series, ancient DNA provides direct access to the most important dimension in evolution—time, allowing a wealth of fundamental evolutionary processes to be addressed at unprecedented resolution. This review taps into the most recent findings in ancient DNA research to present analyses of ancient genomic and metagenomic data.

Applications of Probe Capture Enrichment Next Generation Sequencing for Whole Mitochondrial Genome and 426 Nuclear SNPs for Forensically Challenging Samples

The application of next generation sequencing (NGS) for the analysis of mitochondrial (mt) DNA, short tandem repeats (STRs), and single nucleotide polymorphism (SNPs) has demonstrated great promise for challenging forensic specimens, such as degraded, limited, and mixed samples. Target enrichment using probe capture rather than PCR amplification offers advantages for analysis of degraded DNA since two intact PCR primer sites in the template DNA molecule are not required. Furthermore, NGS software programs can help remove PCR duplicates to determine initial template copy numbers of a shotgun library. Moreover, the same shotgun library prepared from a limited DNA source can be enriched for mtDNA as well as nuclear markers by hybrid capture with the relevant probe panels. Here, we demonstrate the use of this strategy in the analysis of limited and mock degraded samples using our custom probe capture panels for massively parallel sequencing of the whole mtgenome and 426 SNP markers. We also applied the mtgenome capture panel in a mixed sample and analyzed using both phylogenetic and variant frequency based bioinformatics tools to resolve the minor and major contributors. Finally, the results obtained on individual telogen hairs demonstrate the potential of probe capture NGS analysis for both mtDNA and nuclear SNPs for challenging forensic specimens.

Genome of the Tasmanian tiger provides insights into the evolution and demography of an extinct marsupial carnivore

The Tasmanian tiger or thylacine (Thylacinus cynocephalus) was the largest carnivorous Australian marsupial to survive into the modern era. Despite last sharing a common ancestor with the eutherian canids ~160 million years ago, their phenotypic resemblance is considered the most striking example of convergent evolution in mammals. The last known thylacine died in captivity in 1936 and many aspects of the evolutionary history of this unique marsupial apex predator remain unknown. Here we have sequenced the genome of a preserved thylacine pouch young specimen to clarify the phylogenetic position of the thylacine within the carnivorous marsupials, reconstruct its historical demography and examine the genetic basis of its convergence with canids. Retroposon insertion patterns placed the thylacine as the basal lineage in Dasyuromorphia and suggest incomplete lineage sorting in early dasyuromorphs. Demographic analysis indicated a long-term decline in genetic diversity starting well before the arrival of humans in Australia. In spite of their extraordinary phenotypic convergence, comparative genomic analyses demonstrated that amino acid homoplasies between the thylacine and canids are largely consistent with neutral evolution. Furthermore, the genes and pathways targeted by positive selection differ markedly between these species. Together, these findings support models of adaptive convergence driven primarily by cis-regulatory evolution.

Human evolution: a tale from ancient genomes

The field of human ancient DNA (aDNA) has moved from mitochondrial sequencing that suffered from contamination and provided limited biological insights, to become a fully genomic discipline that is changing our conception of human history. Recent successes include the sequencing of extinct hominins, and true population genomic studies of Bronze Age populations. Among the emerging areas of aDNA research, the analysis of past epigenomes is set to provide more new insights into human adaptation and disease susceptibility through time. Starting as a mere curiosity, ancient human genetics has become a major player in the understanding of our evolutionary history.This article is part of the themed issue 'Evo-devo in the genomics era, and the origins of morphological diversity'.

gargammel: a sequence simulator for ancient DNA

Summary

Ancient DNA has emerged as a remarkable tool to infer the history of extinct species and past populations. However, many of its characteristics, such as extensive fragmentation, damage and contamination, can influence downstream analyses. To help investigators measure how these could impact their analyses in silico, we have developed gargammel, a package that simulates ancient DNA fragments given a set of known reference genomes. Our package simulates the entire molecular process from post-mortem DNA fragmentation and DNA damage to experimental sequencing errors, and reproduces most common bias observed in ancient DNA datasets.

Availability and Implementation

The package is publicly available on github: https://grenaud.github.io/gargammel/ and released under the GPL.

Supplementary information

Supplementary data are available at Bioinformatics online.

Optimized mtDNA Control Region Primer Extension Capture Analysis for Forensically Relevant Samples and Highly Compromised mtDNA of Different Age and Origin

The analysis of mitochondrial DNA (mtDNA) has proven useful in forensic genetics and ancient DNA (aDNA) studies, where specimens are often highly compromised and DNA quality and quantity are low. In forensic genetics, the mtDNA control region (CR) is commonly sequenced using established Sanger-type Sequencing (STS) protocols involving fragment sizes down to approximately 150 base pairs (bp). Recent developments include Massively Parallel Sequencing (MPS) of (multiplex) PCR-generated libraries using the same amplicon sizes. Molecular genetic studies on archaeological remains that harbor more degraded aDNA have pioneered alternative approaches to target mtDNA, such as capture hybridization and primer extension capture (PEC) methods followed by MPS. These assays target smaller mtDNA fragment sizes (down to 50 bp or less), and have proven to be substantially more successful in obtaining useful mtDNA sequences from these samples compared to electrophoretic methods. Here, we present the modification and optimization of a PEC method, earlier developed for sequencing the Neanderthal mitochondrial genome, with forensic applications in mind. Our approach was designed for a more sensitive enrichment of the mtDNA CR in a single tube assay and short laboratory turnaround times, thus complying with forensic practices. We characterized the method using sheared, high quantity mtDNA (six samples), and tested challenging forensic samples (n = 2) as well as compromised solid tissue samples (n = 15) up to 8 kyrs of age. The PEC MPS method produced reliable and plausible mtDNA haplotypes that were useful in the forensic context. It yielded plausible data in samples that did not provide results with STS and other MPS techniques. We addressed the issue of contamination by including four generations of negative controls, and discuss the results in the forensic context. We finally offer perspectives for future research to enable the validation and accreditation of the PEC MPS method for final implementation in forensic genetic laboratories.

A Robust Framework for Microbial Archaeology

Microbial archaeology is flourishing in the era of high-throughput sequencing, revealing the agents behind devastating historical plagues, identifying the cryptic movements of pathogens in prehistory, and reconstructing the ancestral microbiota of humans. Here, we introduce the fundamental concepts and theoretical framework of the discipline, then discuss applied methodologies for pathogen identification and microbiome characterization from archaeological samples. We give special attention to the process of identifying, validating, and authenticating ancient microbes using high-throughput DNA sequencing data. Finally, we outline standards and precautions to guide future research in the field.

Ancient DNA reveals the Arctic origin of Viking Age cod from Haithabu, Germany

Knowledge of the range and chronology of historic trade and long-distance transport of natural resources is essential for determining the impacts of past human activities on marine environments. However, the specific biological sources of imported fauna are often difficult to identify, in particular if species have a wide spatial distribution and lack clear osteological or isotopic differentiation between populations. Here, we report that ancient fish-bone remains, despite being porous, brittle, and light, provide an excellent source of endogenous DNA (15-46%) of sufficient quality for whole-genome reconstruction. By comparing ancient sequence data to that of modern specimens, we determine the biological origin of 15 Viking Age (800-1066 CE) and subsequent medieval (1066-1280 CE) Atlantic cod (Gadus morhua) specimens from excavation sites in Germany, Norway, and the United Kingdom. Archaeological context indicates that one of these sites was a fishing settlement for the procurement of local catches, whereas the other localities were centers of trade. Fish from the trade sites show a mixed ancestry and are statistically differentiated from local fish populations. Moreover, Viking Age samples from Haithabu, Germany, are traced back to the North East Arctic Atlantic cod population that has supported the Lofoten fisheries of Norway for centuries. Our results resolve a long-standing controversial hypothesis and indicate that the marine resources of the North Atlantic Ocean were used to sustain an international demand for protein as far back as the Viking Age.

Mining Metagenomic Data Sets for Ancient DNA: Recommended Protocols for Authentication

While a comparatively young area of research, investigations relying on ancient DNA data have been highly valuable in revealing snapshots of genetic variation in both the recent and the not-so-recent past. Born out of a tradition of single-locus PCR-based approaches that often target individual species, stringent criteria for both data acquisition and analysis were introduced early to establish high standards of data quality. Today, the immense volume of data made available through next-generation sequencing has significantly increased the analytical resolution offered by processing ancient tissues and permits parallel analyses of host and microbial communities. The adoption of this new approach to data acquisition, however, requires an accompanying update on methods of DNA authentication, especially given that ancient molecules are expected to exist in low proportions in archaeological material, where an environmental signal is likely to dominate. In this review, we provide a summary of recent data authentication approaches that have been successfully used to distinguish between endogenous and nonendogenous DNA sequences in metagenomic data sets. While our discussion mostly centers on the detection of ancient human and ancient bacterial pathogen DNA, their applicability is far wider.

Paleogenomic Evidence for Multi-generational Mixing between Neolithic Farmers and Mesolithic Hunter-Gatherers in the Lower Danube Basin

The transition from hunting and gathering to farming involved profound cultural and technological changes. In Western and Central Europe, these changes occurred rapidly and synchronously after the arrival of early farmers of Anatolian origin [1-3], who largely replaced the local Mesolithic hunter-gatherers [1, 4-6]. Further east, in the Baltic region, the transition was gradual, with little or no genetic input from incoming farmers [7]. Here we use ancient DNA to investigate the relationship between hunter-gatherers and farmers in the Lower Danube basin, a geographically intermediate area that is characterized by a rapid Neolithic transition but also by the presence of archaeological evidence that points to cultural exchange, and thus possible admixture, between hunter-gatherers and farmers. We recovered four human paleogenomes (1.1× to 4.1× coverage) from Romania spanning a time transect between 8.8 thousand years ago (kya) and 5.4 kya and supplemented them with two Mesolithic genomes (1.7× and 5.3×) from Spain to provide further context on the genetic background of Mesolithic Europe. Our results show major Western hunter-gatherer (WHG) ancestry in a Romanian Eneolithic sample with a minor, but sizeable, contribution from Anatolian farmers, suggesting multiple admixture events between hunter-gatherers and farmers. Dietary stable-isotope analysis of this sample suggests a mixed terrestrial/aquatic diet. Our results provide support for complex interactions among hunter-gatherers and farmers in the Danube basin, demonstrating that in some regions, demic and cultural diffusion were not mutually exclusive, but merely the ends of a continuum for the process of Neolithization.

Genetic Time Travel

At its core, genetics is a historical discipline. Mutations are passed on from generation to generation and accumulate as a result of chance as well as of selection within and between populations and species. However, until recently, geneticists were confined to the study of present-day genetic variation and could only indirectly make inferences about the historical processes that resulted in the variation in present-day gene pools. This "time trap" has now been overcome thanks to the ability to analyze DNA extracted from ancient remains, and this is about to revolutionize several aspects of genetics.

Polymerase chain reaction with nearby primers

DNA analysis of biological specimens containing degraded nucleic acids such as mortal remains, archaeological artefacts, forensic samples etc. has gained more attention in recent years. DNA extracted from these samples is often inapplicable for conventional polymerase chain reaction (PCR), so for its amplification the nearby primers are commonly used. Here we report the data that clarify the features of PCR with nearby and abutting primers. We have shown that the proximity of primers leads to significant reduction of the reaction time and ensures the successful performance of DNA amplification even in the presence of PCR inhibitors. The PCR with abutting primers is usually characterized by the absence of nonspecific amplification products that causes extreme sensitivity with limit of detection on single copy level. The feasibility of PCR with abutting primers was demonstrated on species identification of 100 years old rotten wood.

A Molecular Approach to the Sexing of the Triple Burial at the Upper Paleolithic Site of Dolní Věstonice

In the past decades ancient DNA research has brought numerous insights to archaeological research where traditional approaches were limited. The determination of sex in human skeletal remains is often challenging for physical anthropologists when dealing with incomplete, juvenile or pathological specimens. Molecular approaches allow sexing on the basis of sex-specific markers or by calculating the ratio of DNA derived from different chromosomes. Here we propose a novel approach that relies on the ratio of X chromosome-derived shotgun sequencing data to the autosomal coverage, thus establishing the probability of an XX or XY karyotype. Applying this approach to the individuals of the Upper Paleolithic triple burial of Dolní Věstonice reveals that all three skeletons, including the individual DV 15, whose sex has long been debated due to a pathological condition, were male.

DNA Sequencing in Cultural Heritage

During the last three decades, DNA analysis on degraded samples revealed itself as an important research tool in anthropology, archaeozoology, molecular evolution, and population genetics. Application on topics such as determination of species origin of prehistoric and historic objects, individual identification of famous personalities, characterization of particular samples important for historical, archeological, or evolutionary reconstructions, confers to the paleogenetics an important role also for the enhancement of cultural heritage. A really fast improvement in methodologies in recent years led to a revolution that permitted recovering even complete genomes from highly degraded samples with the possibility to go back in time 400,000 years for samples from temperate regions and 700,000 years for permafrozen remains and to analyze even more recent material that has been subjected to hard biochemical treatments. Here we propose a review on the different methodological approaches used so far for the molecular analysis of degraded samples and their application on some case studies.

Toward high-resolution population genomics using archaeological samples

The term ‘ancient DNA’ (aDNA) is coming of age, with over 1,200 hits in the PubMed database, beginning in the early 1980s with the studies of ‘molecular paleontology’. Rooted in cloning and limited sequencing of DNA from ancient remains during the pre-PCR era, the field has made incredible progress since the introduction of PCR and next-generation sequencing. Over the last decade, aDNA analysis ushered in a new era in genomics and became the method of choice for reconstructing the history of organisms, their biogeography, and migration routes, with applications in evolutionary biology, population genetics, archaeogenetics, paleo-epidemiology, and many other areas. This change was brought by development of new strategies for coping with the challenges in studying aDNA due to damage and fragmentation, scarce samples, significant historical gaps, and limited applicability of population genetics methods. In this review, we describe the state-of-the-art achievements in aDNA studies, with particular focus on human evolution and demographic history. We present the current experimental and theoretical procedures for handling and analysing highly degraded aDNA. We also review the challenges in the rapidly growing field of ancient epigenomics. Advancement of aDNA tools and methods signifies a new era in population genetics and evolutionary medicine research.

Ancient DNA and human history

We review studies of genomic data obtained by sequencing hominin fossils with particular emphasis on the unique information that ancient DNA (aDNA) can provide about the demographic history of humans and our closest relatives. We concentrate on nuclear genomic sequences that have been published in the past few years. In many cases, particularly in the Arctic, the Americas, and Europe, aDNA has revealed historical demographic patterns in a way that could not be resolved by analyzing present-day genomes alone. Ancient DNA from archaic hominins has revealed a rich history of admixture between early modern humans, Neanderthals, and Denisovans, and has allowed us to disentangle complex selective processes. Information from aDNA studies is nowhere near saturation, and we believe that future aDNA sequences will continue to change our understanding of hominin history.