Genome-scale sequencing and analysis of human, wolf, and bison DNA from 25,000-year-old sediment

Cave sediments have been shown to preserve ancient DNA but so far have not yielded the genome-scale information of skeletal remains. We retrieved and analyzed human and mammalian nuclear and mitochondrial environmental "shotgun" genomes from a single 25,000-year-old Upper Paleolithic sediment sample from Satsurblia cave, western Georgia:first, a human environmental genome with substantial basal Eurasian ancestry, which was an ancestral component of the majority of post-Ice Age people in the Near East, North Africa, and parts of Europe; second, a wolf environmental genome that is basal to extant Eurasian wolves and dogs and represents a previously unknown, likely extinct, Caucasian lineage; and third, a European bison environmental genome that is basal to present-day populations, suggesting that population structure has been substantially reshaped since the Last Glacial Maximum. Our results provide new insights into the Late Pleistocene genetic histories of these three species and demonstrate that direct shotgun sequencing of sediment DNA, without target enrichment methods, can yield genome-wide data informative of ancestry and phylogenetic relationships.

Pleistocene sediment DNA reveals hominin and faunal turnovers at Denisova Cave

Denisova Cave in southern Siberia is the type locality of the Denisovans, an archaic hominin group who were related to Neanderthals^1–4. The dozen hominin remains recovered from the deposits also include Neanderthals^5,6 and the child of a Neanderthal and a Denisovan⁷, which suggests that Denisova Cave was a contact zone between these archaic hominins. However, uncertainties persist about the order in which these groups appeared at the site, the timing and environmental context of hominin occupation, and the association of particular hominin groups with archaeological assemblages^5,8–11. Here we report the analysis of DNA from 728 sediment samples that were collected in a grid-like manner from layers dating to the Pleistocene epoch. We retrieved ancient faunal and hominin mitochondrial (mt)DNA from 685 and 175 samples, respectively. The earliest evidence for hominin mtDNA is of Denisovans, and is associated with early Middle Palaeolithic stone tools that were deposited approximately 250,000 to 170,000 years ago; Neanderthal mtDNA first appears towards the end of this period. We detect a turnover in the mtDNA of Denisovans that coincides with changes in the composition of faunal mtDNA, and evidence that Denisovans and Neanderthals occupied the site repeatedly—possibly until, or after, the onset of the Initial Upper Palaeolithic at least 45,000 years ago, when modern human mtDNA is first recorded in the sediments.

Ancient mitochondrial DNA from sediments reveals the sequence of Denisovan, Neanderthal and faunal occupation of Denisova Cave, and evidence for the appearance of modern humans at least 45,000 years ago.

Environmental genomics of Late Pleistocene black bears and giant short-faced bears

Analysis of ancient environmental DNA (eDNA) has revolutionized our ability to describe biological communities in space and time,1-3 by allowing for parallel sequencing of DNA from all trophic levels.4-8 However, because environmental samples contain sparse and fragmented data from multiple individuals, and often contain closely related species,9 the field of ancient eDNA has so far been limited to organellar genomes in its contribution to population and phylogenetic studies.5,6,10,11 This is in contrast to data from fossils12,13 where full-genome studies are routine, despite these being rare and their destruction for sequencing undesirable.14-16 Here, we report the retrieval of three low-coverage (0.03×) environmental genomes from American black bear (Ursus americanus) and a 0.04× environmental genome of the extinct giant short-faced bear (Arctodus simus) from cave sediment samples from northern Mexico dated to 16-14 thousand calibrated years before present (cal kyr BP), which we contextualize with a new high-coverage (26×) and two lower-coverage giant short-faced bear genomes obtained from fossils recovered from Yukon Territory, Canada, which date to ∼22-50 cal kyr BP. We show that the Late Pleistocene black bear population in Mexico is ancestrally related to the present-day Eastern American black bear population, and that the extinct giant short-faced bears present in Mexico were deeply divergent from the earlier Beringian population. Our findings demonstrate the ability to separately analyze genomic-scale DNA sequences of closely related species co-preserved in environmental samples, which brings the use of ancient eDNA into the era of population genomics and phylogenetics.

Peopling of the Americas as inferred from ancient genomics

In less than a decade, analyses of ancient genomes have transformed our understanding of the Indigenous peopling and population history of the Americas. These studies have shown that this history, which began in the late Pleistocene epoch and continued episodically into the Holocene epoch, was far more complex than previously thought. It is now evident that the initial dispersal involved the movement from northeast Asia of distinct and previously unknown populations, including some for whom there are no currently known descendants. The first peoples, once south of the continental ice sheets, spread widely, expanded rapidly and branched into multiple populations. Their descendants-over the next fifteen millennia-experienced varying degrees of isolation, admixture, continuity and replacement, and their genomes help to illuminate the relationships among major subgroups of Native American populations. Notably, all ancient individuals in the Americas, save for later-arriving Arctic peoples, are more closely related to contemporary Indigenous American individuals than to any other population elsewhere, which challenges the claim-which is based on anatomical evidence-that there was an early, non-Native American population in the Americas. Here we review the patterns revealed by ancient genomics that help to shed light on the past peoples who created the archaeological landscape, and together lead to deeper insights into the population and cultural history of the Americas.

Reconstruction of ancient microbial genomes from the human gut

Loss of gut microbial diversity^1–6 in industrial populations is associated with chronic diseases⁷, underscoring the importance of studying our ancestral gut microbiome. However, relatively little is known about the composition of pre-industrial gut microbiomes. Here we performed a large-scale de novo assembly of microbial genomes from palaeofaeces. From eight authenticated human palaeofaeces samples (1,000–2,000 years old) with well-preserved DNA from southwestern USA and Mexico, we reconstructed 498 medium- and high-quality microbial genomes. Among the 181 genomes with the strongest evidence of being ancient and of human gut origin, 39% represent previously undescribed species-level genome bins. Tip dating suggests an approximate diversification timeline for the key human symbiont Methanobrevibacter smithii. In comparison to 789 present-day human gut microbiome samples from eight countries, the palaeofaeces samples are more similar to non-industrialized than industrialized human gut microbiomes. Functional profiling of the palaeofaeces samples reveals a markedly lower abundance of antibiotic-resistance and mucin-degrading genes, as well as enrichment of mobile genetic elements relative to industrial gut microbiomes. This study facilitates the discovery and characterization of previously undescribed gut microorganisms from ancient microbiomes and the investigation of the evolutionary history of the human gut microbiota through genome reconstruction from palaeofaeces.

Ancient microbiomes from palaeofaeces are more similar to non-industrialized than industrialized human gut microbiomes regardless of geography, but 39% of their de novo reconstructed genomes represent previously undescribed microbial species.

Earliest known human burial in Africa

The origin and evolution of hominin mortuary practices are topics of intense interest and debate^1-3. Human burials dated to the Middle Stone Age (MSA) are exceedingly rare in Africa and unknown in East Africa^1-6. Here we describe the partial skeleton of a roughly 2.5- to 3.0-year-old child dating to 78.3 ± 4.1 thousand years ago, which was recovered in the MSA layers of Panga ya Saidi (PYS), a cave site in the tropical upland coast of Kenya^7,8. Recent excavations have revealed a pit feature containing a child in a flexed position. Geochemical, granulometric and micromorphological analyses of the burial pit content and encasing archaeological layers indicate that the pit was deliberately excavated. Taphonomical evidence, such as the strict articulation or good anatomical association of the skeletal elements and histological evidence of putrefaction, support the in-place decomposition of the fresh body. The presence of little or no displacement of the unstable joints during decomposition points to an interment in a filled space (grave earth), making the PYS finding the oldest known human burial in Africa. The morphological assessment of the partial skeleton is consistent with its assignment to Homo sapiens, although the preservation of some primitive features in the dentition supports increasing evidence for non-gradual assembly of modern traits during the emergence of our species. The PYS burial sheds light on how MSA populations interacted with the dead.

Bioinformatic pipelines combining denoising and clustering tools allow for more comprehensive prokaryotic and eukaryotic metabarcoding

Environmental DNA metabarcoding is a powerful tool for studying biodiversity. However, bioinformatic approaches need to adjust to the diversity of taxonomic compartments targeted as well as to each barcode gene specificities. We built and tested a pipeline based on read correction with DADA2 allowing analysing metabarcoding data from prokaryotic (16S) and eukaryotic (18S, COI) life compartments. We implemented the option to cluster amplicon sequence variants (ASVs) into operational taxonomic units (OTUs) with swarm, a network-based clustering algorithm, and the option to curate ASVs/OTUs using LULU. Finally, taxonomic assignment was implemented via the Ribosomal Database Project Bayesian classifier (RDP) and BLAST. We validated this pipeline with ribosomal and mitochondrial markers using metazoan mock communities and 42 deep-sea sediment samples. The results show that ASVs and OTUs describe different levels of biotic diversity, the choice of which depends on the research questions. They underline the advantages and complementarity of clustering and LULU-curation for producing metazoan biodiversity inventories at a level approaching the one obtained using morphological criteria. While clustering removes intraspecific variation, LULU effectively removes spurious clusters, originating from errors or intragenomic variability. Swarm clustering affected alpha and beta diversity differently depending on genetic marker. Specifically, d-values > 1 appeared to be less appropriate with 18S for metazoans. Similarly, increasing LULU's minimum ratio level proved essential to avoid losing species in sample-poor data sets. Comparing BLAST and RDP underlined that accurate assignments of deep-sea species can be obtained with RDP, but highlighted the need for a concerted effort to build comprehensive, ecosystem-specific databases.

Towards the Optimization of eDNA/eRNA Sampling Technologies for Marine Biosecurity Surveillance

The field of eDNA is growing exponentially in response to the need for detecting rare and invasive species for management and conservation decisions. Developing technologies and standard protocols within the biosecurity sector must address myriad challenges associated with marine environments, including salinity, temperature, advective and deposition processes, hydrochemistry and pH, and contaminating agents. These approaches must also provide a robust framework that meets the need for biosecurity management decisions regarding threats to human health, environmental resources, and economic interests, especially in areas with limited clean-laboratory resources and experienced personnel. This contribution aims to facilitate dialogue and innovation within this sector by reviewing current approaches for sample collection, post-sampling capture and concentration of eDNA, preservation, and extraction, all through a biosecurity monitoring lens.

Unearthing Neanderthal population history using nuclear and mitochondrial DNA from cave sediments

Bones and teeth are important sources of Pleistocene hominin DNA, but are rarely recovered at archaeological sites. Mitochondrial DNA (mtDNA) has been retrieved from cave sediments but provides limited value for studying population relationships. We therefore developed methods for the enrichment and analysis of nuclear DNA from sediments and applied them to cave deposits in western Europe and southern Siberia dated to between 200,000 and 50,000 years ago. We detected a population replacement in northern Spain about 100,000 years ago, which was accompanied by a turnover of mtDNA. We also identified two radiation events in Neanderthal history during the early part of the Late Pleistocene. Our work lays the ground for studying the population history of ancient hominins from trace amounts of nuclear DNA in sediments.

Successful extraction of insect DNA from recent copal inclusions: limits and perspectives

Insects entombed in copal, the sub-fossilized resin precursor of amber, represent a potential source of genetic data for extinct and extant, but endangered or elusive, species. Despite several studies demonstrated that it is not possible to recover endogenous DNA from insect inclusions, the preservation of biomolecules in fossilized resins samples is still under debate. In this study, we tested the possibility of obtaining endogenous ancient DNA (aDNA) molecules from insects preserved in copal, applying experimental protocols specifically designed for aDNA recovery. We were able to extract endogenous DNA molecules from one of the two samples analyzed, and to identify the taxonomic status of the specimen. Even if the sample was found well protected from external contaminants, the recovered DNA was low concentrated and extremely degraded, compared to the sample age. We conclude that it is possible to obtain genomic data from resin-entombed organisms, although we discourage aDNA analysis because of the destructive method of extraction protocols and the non-reproducibility of the results.

Reproducible, portable, and efficient ancient genome reconstruction with nf-core/eager

The broadening utilisation of ancient DNA to address archaeological, palaeontological, and biological questions is resulting in a rising diversity in the size of laboratories and scale of analyses being performed. In the context of this heterogeneous landscape, we present an advanced, and entirely redesigned and extended version of the EAGER pipeline for the analysis of ancient genomic data. This Nextflow pipeline aims to address three main themes: accessibility and adaptability to different computing configurations, reproducibility to ensure robust analytical standards, and updating the pipeline to the latest routine ancient genomic practices. The new version of EAGER has been developed within the nf-core initiative to ensure high-quality software development and maintenance support; contributing to a long-term life-cycle for the pipeline. nf-core/eager will assist in ensuring that a wider range of ancient DNA analyses can be applied by a diverse range of research groups and fields.

Reevaluating the timing of Neanderthal disappearance in Northwest Europe

Elucidating when Neanderthal populations disappeared from Eurasia is a key question in paleoanthropology, and Belgium is one of the key regions for studying the Middle to Upper Paleolithic transition. Previous radiocarbon dating placed the Spy Neanderthals among the latest surviving Neanderthals in Northwest Europe with reported dates as young as 23,880 ± 240 B.P. (OxA-8912). Questions were raised, however, regarding the reliability of these dates. Soil contamination and carbon-based conservation products are known to cause problems during the radiocarbon dating of bulk collagen samples. Employing a compound-specific approach that is today the most efficient in removing contamination and ancient genomic analysis, we demonstrate here that previous dates produced on Neanderthal specimens from Spy were inaccurately young by up to 10,000 y due to the presence of unremoved contamination. Our compound-specific radiocarbon dates on the Neanderthals from Spy and those from Engis and Fonds-de-Forêt demonstrate that they disappeared from Northwest Europe at 44,200 to 40,600 cal B.P. (at 95.4% probability), much earlier than previously suggested. Our data contribute significantly to refining models for Neanderthal disappearance in Europe and, more broadly, show that chronometric models regarding the appearance or disappearance of animal or hominin groups should be based only on radiocarbon dates obtained using robust pretreatment methods.

Environmental palaeogenomic reconstruction of an Ice Age algal population

Palaeogenomics has greatly increased our knowledge of past evolutionary and ecological change, but has been restricted to the study of species that preserve either as or within fossils. Here we show the potential of shotgun metagenomics to reveal population genomic information for a taxon that does not preserve in the body fossil record, the algae Nannochloropsis. We shotgun sequenced two lake sediment samples dated to the Last Glacial Maximum and reconstructed full chloroplast and mitochondrial genomes to explore within-lake population genomic variation. This revealed two major haplogroups for each organellar genome, which could be assigned to known varieties of N. limnetica, although we show that at least three haplotypes were present using our minimum haplotype diversity estimation method. These approaches demonstrate the utility of lake sedimentary ancient DNA (sedaDNA) for population genomic analysis, thereby opening the door to environmental palaeogenomics, which will unlock the full potential of sedaDNA.

Lammers et al. use sedimentary ancient DNA to reconstruct palaeogenomes of Nannochloropsis. This study demonstrates the value of sedaDNA for palaeogenomic reconstructions and population genomic analysis.

Hybridization capture of larch (Larix Mill.) chloroplast genomes from sedimentary ancient DNA reveals past changes of Siberian forest

Siberian larch (Larix Mill.) forests dominate vast areas of northern Russia and contribute important ecosystem services to the world. It is important to understand the past dynamics of larches in order to predict their likely response to a changing climate in the future. Sedimentary ancient DNA extracted from lake sediment cores can serve as archives to study past vegetation. However, the traditional method of studying sedimentary ancient DNA-metabarcoding-focuses on small fragments, which cannot resolve Larix to species level nor allow a detailed study of population dynamics. Here, we use shotgun sequencing and hybridization capture with long-range PCR-generated baits covering the complete Larix chloroplast genome to study Larix populations from a sediment core reaching back to 6700 years from the Taymyr region in northern Siberia. In comparison with shotgun sequencing, hybridization capture results in an increase in taxonomically classified reads by several orders of magnitude and the recovery of complete chloroplast genomes of Larix. Variation in the chloroplast reads corroborates an invasion of Larix gmelinii into the range of Larix sibirica before 6700 years ago. Since then, both species have been present at the site, although larch populations have decreased with only a few trees remaining in what was once a forested area. This study demonstrates for the first time that hybridization capture applied directly to ancient DNA of plants extracted from lake sediments can provide genome-scale information and is a viable tool for studying past genomic changes in populations of single species, irrespective of a preservation as macrofossil.

Lake Sedimentary DNA Research on Past Terrestrial and Aquatic Biodiversity: Overview and Recommendations

The use of lake sedimentary DNA to track the long-term changes in both terrestrial and aquatic biota is a rapidly advancing field in paleoecological research. Although largely applied nowadays, knowledge gaps remain in this field and there is therefore still research to be conducted to ensure the reliability of the sedimentary DNA signal. Building on the most recent literature and seven original case studies, we synthesize the state-of-the-art analytical procedures for effective sampling, extraction, amplification, quantification and/or generation of DNA inventories from sedimentary ancient DNA (sedaDNA) via high-throughput sequencing technologies. We provide recommendations based on current knowledge and best practises.

Origins of modern human ancestry

New finds in the palaeoanthropological and genomic records have changed our view of the origins of modern human ancestry. Here we review our current understanding of how the ancestry of modern humans around the globe can be traced into the deep past, and which ancestors it passes through during our journey back in time. We identify three key phases that are surrounded by major questions, and which will be at the frontiers of future research. The most recent phase comprises the worldwide expansion of modern humans between 40 and 60 thousand years ago (ka) and their last known contacts with archaic groups such as Neanderthals and Denisovans. The second phase is associated with a broadly construed African origin of modern human diversity between 60 and 300 ka. The oldest phase comprises the complex separation of modern human ancestors from archaic human groups from 0.3 to 1 million years ago. We argue that no specific point in time can currently be identified at which modern human ancestry was confined to a limited birthplace, and that patterns of the first appearance of anatomical or behavioural traits that are used to define Homo sapiens are consistent with a range of evolutionary histories.

Evaluation of DNA Extraction Methods Developed for Forensic and Ancient DNA Applications Using Bone Samples of Different Age

The efficient extraction of DNA from challenging samples, such as bones, is critical for the success of downstream genotyping analysis in molecular genetic disciplines. Even though the ancient DNA community has developed several protocols targeting small DNA fragments that are typically present in decomposed or old specimens, only recently forensic geneticists have started to adopt those protocols. Here, we compare an ancient DNA extraction protocol (Dabney) with a bone extraction method (Loreille) typically used in forensics. Real-time quantitative PCR and forensically representative typing methods including fragment size analysis and sequencing were used to assess protocol performance. We used four bone samples of different age in replicates to study the effects of both extraction methods. Our results confirm Loreille’s overall increased gain of DNA when enough tissue is available and Dabney’s improved efficiency for retrieving shorter DNA fragments that is beneficial when highly degraded DNA is present. The results suggest that the choice of extraction method needs to be based on available sample, degradation state, and targeted genotyping method. We modified the Dabney protocol by pooling parallel lysates prior to purification to study gain and performance in single tube typing assays and found that up to six parallel lysates lead to an almost linear gain of extracted DNA. These data are promising for further forensic investigations as the adapted Dabney protocol combines increased sensitivity for degraded DNA with necessary total DNA amount for forensic applications.

Finding karstic caves and rockshelters in the Inner Asian mountain corridor using predictive modelling and field survey

The area of the Inner Asian Mountain Corridor (IAMC) follows the foothills and piedmont zones around the northern limits of Asia's interior mountains, connecting two important areas for human evolution: the Fergana valley and the Siberian Altai. Prior research has suggested the IAMC may have provided an area of connected refugia from harsh climates during the Pleistocene. To date, this region contains very few secure, dateable Pleistocene sites, but its widely available carbonate units present an opportunity for discovering cave sites, which generally preserve longer sequences and organic remains. Here we present two models for predicting karstic cave and rockshelter features in the Kazakh portion of the IAMC. The 2018 model used a combination of lithological data and unsupervised landform classification, while the 2019 model used feature locations from the results of our 2017-2018 field surveys in a supervised classification using a minimum-distance classifier and morphometric features derived from the ASTER digital elevation model (DEM). We present the results of two seasons of survey using two iterations of the karstic cave models (2018 and 2019), and evaluate their performance during survey. In total, we identified 105 cave and rockshelter features from 2017-2019. We conclude that this model-led approach significantly reduces the target area for foot survey.

Examining Natural History through the Lens of Palaeogenomics

The many high-resolution tools that are uniquely applicable to specimens from the Quaternary period (the past ~2.5 Ma) provide an opportunity to cross-validate data and test hypotheses based on the morphology and distribution of fossils. Among these tools is palaeogenomics - the genome-scale sequencing of genetic material from ancient specimens - that can provide direct insight into ecology and evolution, potentially improving the accuracy of inferences about past ecological communities over longer timescales. Palaeogenomics has revealed instances of over- and underestimation of extinct diversity, detected cryptic faunal migration and turnover, allowed quantification of widespread sex biases and sexual dimorphism in the fossil record, revealed past hybridisation events and hybrid individuals, and has highlighted previously unrecognised routes of zoonotic disease transfer.

Ancient DNA from Guam and the peopling of the Pacific

Humans reached the Mariana Islands in the western Pacific by ∼3,500 y ago, contemporaneous with or even earlier than the initial peopling of Polynesia. They crossed more than 2,000 km of open ocean to get there, whereas voyages of similar length did not occur anywhere else until more than 2,000 y later. Yet, the settlement of Polynesia has received far more attention than the settlement of the Marianas. There is uncertainty over both the origin of the first colonizers of the Marianas (with different lines of evidence suggesting variously the Philippines, Indonesia, New Guinea, or the Bismarck Archipelago) as well as what, if any, relationship they might have had with the first colonizers of Polynesia. To address these questions, we obtained ancient DNA data from two skeletons from the Ritidian Beach Cave Site in northern Guam, dating to ∼2,200 y ago. Analyses of complete mitochondrial DNA genome sequences and genome-wide SNP data strongly support ancestry from the Philippines, in agreement with some interpretations of the linguistic and archaeological evidence, but in contradiction to results based on computer simulations of sea voyaging. We also find a close link between the ancient Guam skeletons and early Lapita individuals from Vanuatu and Tonga, suggesting that the Marianas and Polynesia were colonized from the same source population, and raising the possibility that the Marianas played a role in the eventual settlement of Polynesia.

A revised AMS and tephra chronology for the Late Middle to Early Upper Paleolithic occupations of Ortvale Klde, Republic of Georgia

The nature and timing of the shift from the Late Middle Paleolithic (LMP) to the Early Upper Paleolithic (EUP) varied geographically, temporally, and substantively across the Near East and Eurasia; however, the result of this process was the archaeological disappearance of Middle Paleolithic technologies across the length and breadth of their geographic distribution. Ortvale Klde rockshelter (Republic of Georgia) contains the most detailed LMP-EUP archaeological sequence in the Caucasus, an environmentally and topographically diverse region situated between southwest Asia and Europe. Tephrochronological investigations at the site reveal volcanic ash (tephra) from various volcanic sources and provide a tephrostratigraphy for the site that will facilitate future correlations in the region. We correlate one of the cryptotephra layers to the large, caldera-forming Nemrut Formation eruption (30,000 years ago) from Nemrut volcano in Turkey. We integrate this tephrochronological constraint with new radiocarbon dates and published ages in an OxCal Bayesian age model to produce a revised chronology for the site. This model increases the ages for the end of the LMP (∼47.5-44.2 ka cal BP) and appearance of the EUP (∼46.7-43.6 ka cal BP) at Ortvale Klde, which are earlier than those currently reported for other sites in the Caucasus but similar to estimates for specific sites in southwest Asia and eastern Europe. These data, coupled with archaeological, stratigraphic, and taphonomic observations, suggest that at Ortvale Klde, (1) the appearance of EUP technologies of bone and stone has no technological roots in the preceding LMP, (2) a LMP population vacuum likely preceded the appearance of these EUP technologies, and (3) the systematic combination of tephra correlations and absolute dating chronologies promises to substantially improve our inter-regional understanding of this critical time interval of human evolution and the potential interconnectedness of hominins at different sites.

Genome-scale target capture of mitochondrial and nuclear environmental DNA from water samples

Environmental DNA (eDNA) provides a promising supplement to traditional sampling methods for population genetic inferences, but current studies have almost entirely focused on short mitochondrial markers. Here, we develop one mitochondrial and one nuclear set of target capture probes for the whale shark (Rhincodon typus) and test them on seawater samples collected in Qatar to investigate the potential of target capture for eDNA‐based population studies. The mitochondrial target capture successfully retrieved ~235× (90× − 352× per base position) coverage of the whale shark mitogenome. Using a minor allele frequency of 5%, we find 29 variable sites throughout the mitogenome, indicative of at least five contributing individuals. We also retrieved numerous mitochondrial reads from an abundant nontarget species, mackerel tuna (Euthynnus affinis), showing a clear relationship between sequence similarity to the capture probes and the number of captured reads. The nuclear target capture probes retrieved only a few reads and polymorphic variants from the whale shark, but we successfully obtained millions of reads and thousands of polymorphic variants with different allele frequencies from E. affinis. We demonstrate that target capture of complete mitochondrial genomes and thousands of nuclear loci is possible from aquatic eDNA samples. Our results highlight that careful probe design, taking into account the range of divergence between target and nontarget sequences as well as presence of nontarget species at the sampling site, is crucial to consider. eDNA sampling coupled with target capture approaches provide an efficient means with which to retrieve population genomic data from aggregating and spawning aquatic species.

Micro-contextual identification of archaeological lipid biomarkers using resin-impregnated sediment slabs

Characterizing organic matter preserved in archaeological sediment is crucial to behavioral and paleoenvironmental investigations. This task becomes particularly challenging when considering microstratigraphic complexity. Most of the current analytical methods rely on loose sediment samples lacking spatial and temporal resolution at a microstratigraphic scale, adding uncertainty to the results. Here, we explore the potential of targeted molecular and isotopic biomarker analysis on polyester resin-impregnated sediment slabs from archaeological micromorphology, a technique that provides microstratigraphic control. We performed gas chromatography-mass spectrometry (GC-MS) and gas chromatography-isotope ratio mass spectromety (GC-IRMS) analyses on a set of samples including drill dust from resin-impregnated experimental and archaeological samples, loose samples from the same locations and resin control samples to assess the degree of interference of polyester resin in the GC-MS and Carbon-IRMS signals of different lipid fractions (n-alkanes, aromatics, n-ketones, alcohols, fatty acids and other high polarity lipids). The results show that biomarkers within the n-alkane, aromatic, n-ketone, and alcohol fractions can be identified. Further work is needed to expand the range of identifiable lipid biomarkers. This study represents the first micro-contextual approach to archaeological lipid biomarkers and contributes to the advance of archaeological science by adding a new method to obtain behavioral or paleoenvironmental proxies.

A systematic investigation of human DNA preservation in medieval skeletons

Ancient DNA (aDNA) analyses necessitate the destructive sampling of archaeological material. Currently, the cochlea, part of the osseous inner ear located inside the petrous pyramid, is the most sought after skeletal element for molecular analyses of ancient humans as it has been shown to yield high amounts of endogenous DNA. However, destructive sampling of the petrous pyramid may not always be possible, particularly in cases where preservation of skeletal morphology is of top priority. To investigate alternatives, we present a survey of human aDNA preservation for each of ten skeletal elements in a skeletal collection from Medieval Germany. Through comparison of human DNA content and quality we confirm best performance of the petrous pyramid and identify seven additional sampling locations across four skeletal elements that yield adequate aDNA for most applications in human palaeogenetics. Our study provides a better perspective on DNA preservation across the human skeleton and takes a further step toward the more responsible use of ancient materials in human aDNA studies.

Estimating molecular preservation of the intestinal microbiome via metagenomic analyses of latrine sediments from two medieval cities

Ancient latrine sediments, which contain the concentrated collective biological waste of past whole human communities, have the potential to be excellent proxies for human gastrointestinal health on the population level. A rich body of literature explores their use to detect the presence of gut-associated eukaryotic parasites through microscopy, immunoassays and genetics. Despite this interest, a lack of studies have explored the whole genetic content of ancient latrine sediments through consideration not only of gut-associated parasites, but also of core community gut microbiome signals that remain from the group that used the latrine. Here, we present a metagenomic analysis of bulk sediment from medieval latrines in Riga (Latvia) and Jerusalem. Our analyses reveal survival of microbial DNA representative of intestinal flora as well as numerous parasites. These data are compared against parasite taxon identifications obtained via microscopy and ELISA techniques. Together, these findings provide a first glimpse into the rich prokaryotic and eukaryotic intestinal flora of pre-industrial agricultural populations, which may give a better context for interpreting the health of modern microbiomes.

This article is part of the theme issue ‘Insights into health and disease from ancient biomolecules’.

Manual and automated preparation of single-stranded DNA libraries for the sequencing of DNA from ancient biological remains and other sources of highly degraded DNA

It has been shown that highly fragmented DNA is most efficiently converted into DNA libraries for sequencing if both strands of the DNA fragments are processed independently. We present an updated protocol for library preparation from single-stranded DNA, which is based on the splinted ligation of an adapter oligonucleotide to the 3' ends of single DNA strands, the synthesis of a complementary strand using a DNA polymerase and the addition of a 5' adapter via blunt-end ligation. The efficiency of library preparation is determined individually for each sample using a spike-in oligonucleotide. The whole workflow, including library preparation, quantification and amplification, requires two work days for up to 16 libraries. Alternatively, we provide documentation and electronic protocols enabling automated library preparation of 96 samples in parallel on a Bravo NGS Workstation (Agilent Technologies). After library preparation, molecules with uninformative short inserts (shorter than ~30-35 base pairs) can be removed by polyacrylamide gel electrophoresis if desired.

Evidence of human occupation in Mexico around the Last Glacial Maximum

The initial colonization of the Americas remains a highly debated topic¹, and the exact timing of the first arrivals is unknown. The earliest archaeological record of Mexico-which holds a key geographical position in the Americas-is poorly known and understudied. Historically, the region has remained on the periphery of research focused on the first American populations². However, recent investigations provide reliable evidence of a human presence in the northwest region of Mexico^3,4, the Chiapas Highlands⁵, Central Mexico⁶ and the Caribbean coast^7-9 during the Late Pleistocene and Early Holocene epochs. Here we present results of recent excavations at Chiquihuite Cave-a high-altitude site in central-northern Mexico-that corroborate previous findings in the Americas^10-17of cultural evidence that dates to the Last Glacial Maximum (26,500-19,000 years ago)¹⁸, and which push back dates for human dispersal to the region possibly as early as 33,000-31,000 years ago. The site yielded about 1,900 stone artefacts within a 3-m-deep stratified sequence, revealing a previously unknown lithic industry that underwent only minor changes over millennia. More than 50 radiocarbon and luminescence dates provide chronological control, and genetic, palaeoenvironmental and chemical data document the changing environments in which the occupants lived. Our results provide new evidence for the antiquity of humans in the Americas, illustrate the cultural diversity of the earliest dispersal groups (which predate those of the Clovis culture) and open new directions of research.

The timing and effect of the earliest human arrivals in North America

The peopling of the Americas marks a major expansion of humans across the planet. However, questions regarding the timing and mechanisms of this dispersal remain, and the previously accepted model (termed 'Clovis-first')-suggesting that the first inhabitants of the Americas were linked with the Clovis tradition, a complex marked by distinctive fluted lithic points¹-has been effectively refuted. Here we analyse chronometric data from 42 North American and Beringian archaeological sites using a Bayesian age modelling approach, and use the resulting chronological framework to elucidate spatiotemporal patterns of human dispersal. We then integrate these patterns with the available genetic and climatic evidence. The data obtained show that humans were probably present before, during and immediately after the Last Glacial Maximum (about 26.5-19 thousand years ago)^2,3 but that more widespread occupation began during a period of abrupt warming, Greenland Interstadial 1 (about 14.7-12.9 thousand years before AD 2000)⁴. We also identify the near-synchronous commencement of Beringian, Clovis and Western Stemmed cultural traditions, and an overlap of each with the last dates for the appearance of 18 now-extinct faunal genera. Our analysis suggests that the widespread expansion of humans through North America was a key factor in the extinction of large terrestrial mammals.

Human evolutionary history in Eastern Eurasia using insights from ancient DNA

Advances in ancient genomics are providing unprecedented insight into modern human history. Here, we review recent progress uncovering prehistoric populations in Eastern Eurasia based on ancient DNA studies from the Upper Pleistocene to the Holocene. Many ancient populations existed during the Upper Pleistocene of Eastern Eurasia-some with no substantial ancestry related to present-day populations, some with an affinity to East Asians, and some who contributed to Native Americans. By the Holocene, the genetic composition across East Asia greatly shifted, with several substantial migrations. Three are southward: an increase in northern East Asian-related ancestry in southern East Asia; movement of East Asian-related ancestry into Southeast Asia, mixing with Basal Asian ancestry; and movement of southern East Asian ancestry to islands of Southeast Asia and the Southwest Pacific through the expansion of Austronesians. We anticipate that additional ancient DNA will magnify our understanding of the genetic history in Eastern Eurasia.

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.

Persistence of environmental DNA in cultivated soils: implication of this memory effect for reconstructing the dynamics of land use and cover changes

eDNA refers to DNA extracted from an environmental sample with the goal of identifying the occurrence of past or current biological communities in aquatic and terrestrial environments. However, there is currently a lack of knowledge regarding the soil memory effect and its potential impact on lake sediment eDNA records. To investigate this issue, two contrasted sites located in cultivated environments in France were studied. In the first site, soil samples were collected (n = 30) in plots for which the crop rotation history was documented since 1975. In the second site, samples were collected (n = 40) to compare the abundance of currently observed taxa versus detected taxa in cropland and other land uses. The results showed that the last cultivated crop was detected in 100% of the samples as the most abundant. In addition, weeds were the most abundant taxa identified in both sites. Overall, these results illustrate the potential of eDNA analyses for identifying the recent (< 10 years) land cover history of soils and outline the detection of different taxa in cultivated plots. The capacity of detection of plant species grown on soils delivering sediments to lacustrine systems is promising to improve our understanding of sediment transfer processes over short timescales.

Unroll Please: Deciphering the Genetic Code in Scrolls and Other Ancient Materials

The unrelenting development of ancient DNA methods now allows researchers to obtain archaeogenetic data from increasingly diverse sources. In a new study in this issue of Cell, researchers apply the latest DNA technologies to unravel the mysteries of the Dead Sea Scrolls, one of the world's most famous and influential sets of ancient parchments.

Generations of genomes: advances in paleogenomics technology and engagement for Indigenous people of the Americas

For decades, scientists have collected genomic information from Indigenous peoples and their ancestors with the goal of elucidating human migration events, understanding ancestral origins, and identifying ancestral variants contributing to disease. However, such studies may not have offered much benefit to the Indigenous groups who contributed DNA, and many have instead perpetuated stereotypes and other harms. With recent advances in genomic technology facilitating the study of both ancient and present-day DNA, researchers and Indigenous communities have new opportunities to begin collaboratively addressing important questions about human health and history. Yet, while there are increased efforts to ethically engage Indigenous communities, more work is still needed as the discipline struggles to absolve itself of the racialized science and extractive biocolonialism that defined its past.

An optimized method for the extraction of ancient eukaryote DNA from marine sediments

Marine sedimentary ancient DNA (sedaDNA) provides a powerful means to reconstruct marine palaeo-communities across the food web. However, currently there are few optimized sedaDNA extraction protocols available to maximize the yield of small DNA fragments typical of ancient DNA (aDNA) across a broad diversity of eukaryotes. We compared seven combinations of sedaDNA extraction treatments and sequencing library preparations using marine sediments collected at a water depth of 104 m off Maria Island, Tasmania, in 2018. These seven methods contrasted frozen versus refrigerated sediment, bead-beating induced cell lysis versus ethylenediaminetetraacetic acid (EDTA) incubation, DNA binding in silica spin columns versus in silica-solution, diluted versus undiluted DNA in shotgun library preparations to test potential inhibition issues during amplification steps, and size-selection of low molecular-weight (LMW) DNA to increase the extraction efficiency of sedaDNA. Maximum efficiency was obtained from frozen sediments subjected to a combination of EDTA incubation and bead-beating, DNA binding in silica-solution, and undiluted DNA in shotgun libraries, across 45 marine eukaryotic taxa. We present an optimized extraction protocol integrating these steps, with an optional post-library LMW size-selection step to retain DNA fragments of ≤500 base pairs. We also describe a stringent bioinformatic filtering approach for metagenomic data and provide a comprehensive list of contaminants as a reference for future sedaDNA studies. The new extraction and data-processing protocol should improve quantitative paleo-monitoring of eukaryotes from marine sediments, as well as other studies relying on the detection of highly fragmented and degraded eukaryote DNA in sediments.

Hybridization ddRAD-sequencing for population genomics of nonmodel plants using highly degraded historical specimen DNA

Species' responses at the genetic level are key to understanding the long-term consequences of anthropogenic global change. Herbaria document such responses, and, with contemporary sampling, provide high-resolution time-series of plant evolutionary change. Characterizing genetic diversity is straightforward for model species with small genomes and a reference sequence. For nonmodel species-with small or large genomes-diversity is traditionally assessed using restriction-enzyme-based sequencing. However, age-related DNA damage and fragmentation preclude the use of this approach for ancient herbarium DNA. Here, we combine reduced-representation sequencing and hybridization-capture to overcome this challenge and efficiently compare contemporary and historical specimens. Specifically, we describe how homemade DNA baits can be produced from reduced-representation libraries of fresh samples, and used to efficiently enrich historical libraries for the same fraction of the genome to produce compatible sets of sequence data from both types of material. Applying this approach to both Arabidopsis thaliana and the nonmodel plant Cardamine bulbifera, we discovered polymorphisms de novo in an unbiased, reference-free manner. We show that the recovered genetic variation recapitulates known genetic diversity in A. thaliana, and recovers geographical origin in both species and over time, independent of bait diversity. Hence, our method enables fast, cost-efficient, large-scale integration of contemporary and historical specimens for assessment of genome-wide genetic trends over time, independent of genome size and presence of a reference genome.

Non-destructive ZooMS identification reveals strategic bone tool raw material selection by Neandertals

Five nearly identical fragments of specialized bone tools, interpreted as lissoirs (French for “smoothers”), have been found at two Middle Paleolithic sites in southwest France. The finds span three separate archaeological deposits, suggesting continuity in the behavior of late Neandertals. Using standard morphological assessments, we determined that the lissoirs were produced on ribs of medium-sized ungulates. However, since these bones are highly fragmented and anthropogenically modified, species determinations were challenging. Also, conservative curation policy recommends minimizing destructive sampling of rare, fragile, or small artifacts for molecular identification methods. To better understand raw material selection for these five lissoirs, we reassess their taxonomy using a non-destructive ZooMS methodology based on triboelectric capture of collagen. We sampled four storage containers and obtained identifiable MALDI-TOF MS collagen fingerprints, all indicative of the same taxonomic clade, which includes aurochs and bison (Bos sp. and Bison sp.). The fifth specimen, which was stored in a plastic bag, provided no useful MALDI-TOF MS spectra. We show that the choice of large bovid ribs in an archaeological layer dominated by reindeer (Rangifer tarandus) demonstrates strategic selection by these Neandertals. Furthermore, our results highlight the value of a promising technique for the non-destructive analysis of bone artifacts.

Ancient Plant Genomics in Archaeology, Herbaria, and the Environment

The ancient DNA revolution of the past 35 years has driven an explosion in the breadth, nuance, and diversity of questions that are approachable using ancient biomolecules, and plant research has been a constant, indispensable facet of these developments. Using archaeological, paleontological, and herbarium plant tissues, researchers have probed plant domestication and dispersal, plant evolution and ecology, paleoenvironmental composition and dynamics, and other topics across related disciplines. Here, we review the development of the ancient DNA discipline and the role of plant research in its progress and refinement. We summarize our understanding of long-term plant DNA preservation and the characteristics of degraded DNA. In addition, we discuss challenges in ancient DNA recovery and analysis and the laboratory and bioinformatic strategies used to mitigate them. Finally, we review recent applications of ancient plant genomic research.

The Treasure Vault Can be Opened: Large-Scale Genome Skimming Works Well Using Herbarium and Silica Gel Dried Material

Genome skimming has the potential for generating large data sets for DNA barcoding and wider biodiversity genomic studies, particularly via the assembly and annotation of full chloroplast (cpDNA) and nuclear ribosomal DNA (nrDNA) sequences. We compare the success of genome skims of 2051 herbarium specimens from Norway/Polar regions with 4604 freshly collected, silica gel dried specimens mainly from the European Alps and the Carpathians. Overall, we were able to assemble the full chloroplast genome for 67% of the samples and the full nrDNA cluster for 86%. Average insert length, cover and full cpDNA and rDNA assembly were considerably higher for silica gel dried than herbarium-preserved material. However, complete plastid genomes were still assembled for 54% of herbarium samples compared to 70% of silica dried samples. Moreover, there was comparable recovery of coding genes from both tissue sources (121 for silica gel dried and 118 for herbarium material) and only minor differences in assembly success of standard barcodes between silica dried (89% ITS2, 96% matK and rbcL) and herbarium material (87% ITS2, 98% matK and rbcL). The success rate was > 90% for all three markers in 1034 of 1036 genera in 160 families, and only Boraginaceae worked poorly, with 7 genera failing. Our study shows that large-scale genome skims are feasible and work well across most of the land plant families and genera we tested, independently of material type. It is therefore an efficient method for increasing the availability of plant biodiversity genomic data to support a multitude of downstream applications.

A mobile laboratory for ancient DNA analysis

Mobile devices for on-field DNA analysis have been used for medical diagnostics at the point-of-care, forensic investigations and environmental surveys, but still have to be validated for ancient DNA studies. We report here on a mobile laboratory that we setup using commercially available devices, including a compact real-time PCR machine, and describe procedures to perform DNA extraction and analysis from a variety of archeological samples within 4 hours. The process is carried out on 50 mg samples that are identified at the species level using custom TaqMan real-time PCR assays for mitochondrial DNA fragments. We evaluated the potential of this approach in museums lacking facilities for DNA studies by analyzing samples from the Enlène (MIS 2 layer) and the Portel-Ouest cave (MIS 3 deposits), and also performed experiments during an excavation campaign at the Roc-en-Pail (MIS 5) open-air site. Enlène Bovinae bone samples only yielded DNA for the extinct steppe bison (Bison priscus), whereas Portel-Ouest cave coprolites contained cave hyena (Crocuta crocuta spelaea) DNA together, for some of them, with DNA for the European bison sister species/subspecies (Bison schoetensacki/Bb1-X), thus highlighting the cave hyena diet. Roc-en-Pail Bovinae bone and tooth samples also contained DNA for the Bison schoetensacki/Bb1-X clade, and Cervidae bone samples only yielded reindeer (Rangifer tarandus) DNA. Subsequent DNA sequencing analyses confirmed that correct species identification had been achieved using our TaqMan assays, hence validating these assays for future studies. We conclude that our approach enables the rapid genetic characterization of tens of millennia-old archeological samples and is expected to be useful for the on-site screening of museums and freshly excavated samples for DNA content. Because our mobile laboratory is made up of commercially available instruments, this approach is easily accessible to other investigators.

Paleo-metagenomics of North American fossil packrat middens: Past biodiversity revealed by ancient DNA

Fossil rodent middens are powerful tools in paleoecology. In arid parts of western North America, packrat (Neotoma spp.) middens preserve plant and animal remains for tens of thousands of years. Midden contents are so well preserved that fragments of endogenous ancient DNA (aDNA) can be extracted and analyzed across millennia. Here, we explore the use of shotgun metagenomics to study the aDNA obtained from packrat middens up to 32,000 C14 years old. Eleven Illumina HiSeq 2500 libraries were successfully sequenced, and between 0.11% and 6.7% of reads were classified using Centrifuge against the NCBI "nt" database. Eukaryotic taxa identified belonged primarily to vascular plants with smaller proportions mapping to ascomycete fungi, arthropods, chordates, and nematodes. Plant taxonomic diversity in the middens is shown to change through time and tracks changes in assemblages determined by morphological examination of the plant remains. Amplicon sequencing of ITS2 and rbcL provided minimal data for some middens, but failed at amplifying the highly fragmented DNA present in others. With repeated sampling and deep sequencing, analysis of packrat midden aDNA from well-preserved midden material can provide highly detailed characterizations of past communities of plants, animals, bacteria, and fungi present as trace DNA fossils. The prospects for gaining more paleoecological insights from aDNA for rodent middens will continue to improve with optimization of laboratory methods, decreasing sequencing costs, and increasing computational power.

Effects of sample age on data quality from targeted sequencing of museum specimens: what are we capturing in time?

Background

Next generation sequencing (NGS) can recover DNA data from valuable extant and extinct museum specimens. However, archived or preserved DNA is difficult to sequence because of its fragmented, damaged nature, such that the most successful NGS methods for preserved specimens remain sub-optimal. Improving wet-lab protocols and comprehensively determining the effects of sample age on NGS library quality are therefore of vital importance. Here, I examine the relationship between sample age and several indicators of library quality following targeted NGS sequencing of ~ 1300 loci using 271 samples of pinned moth specimens (Helicoverpa armigera) ranging in age from 5 to 117 years.

Results

I find that older samples have lower DNA concentrations following extraction and thus require a higher number of indexing PCR cycles during library preparation. When sequenced reads are aligned to a reference genome or to only the targeted region, older samples have a lower number of sequenced and mapped reads, lower mean coverage, and lower estimated library sizes, while the percentage of adapters in sequenced reads increases significantly as samples become older. Older samples also show the poorest capture success, with lower enrichment and a higher improved coverage anticipated from further sequencing.

Conclusions

Sample age has significant, measurable impacts on the quality of NGS data following targeted enrichment. However, incorporating a uracil-removing enzyme into the blunt end-repair step during library preparation could help to repair DNA damage, and using a method that prevents adapter-dimer formation may result in improved data yields.

Population-level inferences from environmental DNA-Current status and future perspectives

Environmental DNA (eDNA) extracted from water samples has recently shown potential as a valuable source of population genetic information for aquatic macroorganisms. This approach offers several potential advantages compared with conventional tissue-based methods, including the fact that eDNA sampling is noninvasive and generally more cost-efficient. Currently, eDNA approaches have been limited to single-marker studies of mitochondrial DNA (mtDNA), and the relationship between eDNA haplotype composition and true haplotype composition still needs to be thoroughly verified. This will require testing of bioinformatic and statistical software to correct for erroneous sequences, as well as biases and random variation in relative sequence abundances. However, eDNA-based population genetic methods have far-reaching potential for both basic and applied research. In this paper, we present a brief overview of the achievements of eDNA-based population genetics to date, and outline the prospects for future developments in the field, including the estimation of nuclear DNA (nuDNA) variation and epigenetic information. We discuss the challenges associated with eDNA samples as opposed to those of individual tissue samples and assess whether eDNA might offer additional types of information unobtainable with tissue samples. Lastly, we provide recommendations for determining whether an eDNA approach would be a useful and suitable choice in different research settings. We limit our discussion largely to contemporary aquatic systems, but the advantages, challenges, and perspectives can to a large degree be generalized to eDNA studies with a different spatial and temporal focus.

Archaeological evidence for two separate dispersals of Neanderthals into southern Siberia

Neanderthals were once widespread across Europe and western Asia. They also penetrated into the Altai Mountains of southern Siberia, but the geographical origin of these populations and the timing of their dispersal have remained elusive. Here we describe an archaeological assemblage from Chagyrskaya Cave, situated in the Altai foothills, where around 90,000 Middle Paleolithic artifacts and 74 Neanderthal remains have been recovered from deposits dating to between 59 and 49 thousand years ago (age range at 95.4% probability). Environmental reconstructions suggest that the Chagyrskaya hominins were adapted to the dry steppe and hunted bison. Their distinctive toolkit closely resembles Micoquian assemblages from central and eastern Europe, including the northern Caucasus, more than 3,000 kilometers to the west of Chagyrskaya Cave. At other Altai sites, evidence of earlier Neanderthal populations lacking associated Micoquian-like artifacts implies two or more Neanderthal incursions into this region. We identify eastern Europe as the most probable ancestral source region for the Chagyrskaya toolmakers, supported by DNA results linking the Neanderthal remains with populations in northern Croatia and the northern Caucasus, and providing a rare example of a long-distance, intercontinental population movement associated with a distinctive Paleolithic toolkit.

Quaternary DNA: A Multidisciplinary Research Field

The purpose of this Milankovitch review is to explain the significance of Quaternary DNA studies and the importance of the recent methodological advances that have enabled the study of late Quaternary remains in more detail, and the testing of new assumptions in evolutionary biology and phylogeography to reconstruct the past. The topic is wide, and this review is not intended to be an exhaustive account of all the aDNA work performed in the last three decades on late-Quaternary remains. Instead, it is a selection of relevant studies aimed at illustrating how aDNA has been used to reconstruct not only environments of the past, but also the history of many species including our own.