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

Integrating musealized archaeological sediment collections into current geoarchaeological analytical frameworks for sustainable research practices

We present a review of the latest framework achievements in geoarchaeological sciences applied to microstratigraphic and biomolecular studies of prehistoric archaeological contexts, highlighting the importance of musealized archaeological stratigraphies. We assess how today's scientific and technological accomplishments can be tailored for archaeological human ecology studies with analytical ensembles that provide unprecedented results. Sampling and processing workflows originating from resin consolidation and thin section micromorphology of undisturbed blocks of archaeological soils and sediments, guarantee subsampling accuracy at the micrometre scale granting access to individual components otherwise impossible to target: the achievable information yield makes even the smallest soil samples potential sources of pioneering discoveries. Yet, archaeological excavations are still the primary mode of retrieving new soil samples. We argue that, when dealing with archaeological sites that were excavated and partially musealised in the past, the exploration of related museum collections should be prioritized as perspective source of new study samples. Analysing old and potentially very informative samples, with an approach that we define as "Green Archaeology", may represent a source of well-structured primary data as well as a means for planning new excavations, delivering novel discoveries while safeguarding site integrity and promoting Third Mission valorisation of sites and heritage dormant collections.

Sedimentary DNA: Archaeology reinvented

How can we better understand the human agency in prehistory? A recent study explores a novel way of using sedimentary DNA extracted from sediments of a rock shelter to investigate the human impact on the landscape.

Tracing early pastoralism in Central Europe using sedimentary ancient DNA

Central European forests have been shaped by complex human interactions throughout the Holocene, with significant changes following the introduction of domesticated animals in the Neolithic (∼7.5-6.0 ka before present [BP]). However, understanding early pastoral practices and their impact on forests is limited by methods for detecting animal movement across past landscapes. Here, we examine ancient sedimentary DNA (sedaDNA) preserved at the Velký Mamuťák rock shelter in northern Bohemia (Czech Republic), which has been a forested enclave since the early Holocene. We find that domesticated animals, their associated microbiomes, and plants potentially gathered for fodder have clear representation by the Late Neolithic, around 6.0 ka BP, and persist throughout the Bronze Age into recent times. We identify a change in dominant grazing species from sheep to pigs in the Bronze Age (∼4.1-3.0 ka BP) and interpret the impact this had in the mid-Holocene retrogressions that still define the structure of Central European forests today. This study highlights the ability of ancient metagenomics to bridge archaeological and paleoecological methods and provide an enhanced perspective on the roots of the "Anthropocene."

soibean: High-Resolution Taxonomic Identification of Ancient Environmental DNA Using Mitochondrial Pangenome Graphs

Ancient environmental DNA (aeDNA) is becoming a powerful tool to gain insights about past ecosystems, overcoming the limitations of conventional fossil records. However, several methodological challenges remain, particularly for classifying the DNA to species level and conducting phylogenetic analysis. Current methods, primarily tailored for modern datasets, fail to capture several idiosyncrasies of aeDNA, including species mixtures from closely related species and ancestral divergence. We introduce soibean, a novel tool that utilizes mitochondrial pangenomic graphs for identifying species from aeDNA reads. It outperforms existing methods in accurately identifying species from multiple closely related sources within a sample, enhancing phylogenetic analysis for aeDNA. soibean employs a damage-aware likelihood model for precise identification at low coverage with a high damage rate. Additionally, we reconstructed ancestral sequences for soibean's database to handle aeDNA that is highly diverged from modern references. soibean demonstrates effectiveness through simulated data tests and empirical validation. Notably, our method uncovered new empirical results in published datasets, including using porpoise whales as food in a Mesolithic community in Sweden, demonstrating its potential to reveal previously unrecognized findings in aeDNA studies.

Long genetic and social isolation in Neanderthals before their extinction

Neanderthal genomes have been recovered from sites across Eurasia, painting an increasingly complex picture of their populations' structure that mostly indicates that late European Neanderthals belonged to a single metapopulation with no significant evidence of population structure. Here, we report the discovery of a late Neanderthal individual, nicknamed "Thorin," from Grotte Mandrin in Mediterranean France, and his genome. These dentognathic fossils, including a rare example of distomolars, are associated with a rich archeological record of Neanderthal final technological traditions in this region ∼50-42 thousand years ago. Thorin's genome reveals a relatively early divergence of ∼105 ka with other late Neanderthals. Thorin belonged to a population with a small group size that showed no genetic introgression with other known late European Neanderthals, revealing some 50 ka of genetic isolation of his lineage despite them living in neighboring regions. These results have important implications for resolving competing hypotheses about causes of the disappearance of the Neanderthals.

Maximizing efficiency in sedimentary ancient DNA analysis: a novel extract pooling approach

In the last few decades, the field of ancient DNA has taken a new direction towards using sedimentary ancient DNA (sedaDNA) for studying human and mammalian population dynamics as well as past ecosystems. However, the screening of numerous sediment samples from archaeological sites remains a time-consuming and costly endeavor, particularly when targeting hominin DNA. Here, we present a novel high-throughput method that facilitates the fast and efficient analysis of sediment samples by applying a pooled testing approach. This method combines multiple extracts, enabling early parallelization of laboratory procedures and effective aDNA screening. Pooled samples with detectable aDNA signals undergo detailed analysis, while empty pools are discarded. We have successfully applied our method to multiple sediment samples from Middle and Upper Paleolithic sites in Europe, Asia, and Africa. Notably, our results reveal that an aDNA signal remains discernible even when pooled with four negative samples. We also demonstrate that the DNA yield of double-stranded libraries increases significantly when reducing the extract input, potentially mitigating the effects of inhibition. By embracing this innovative approach, researchers can analyze large numbers of sediment samples for aDNA preservation, achieving significant cost reductions of up to 70% and reducing hands-on laboratory time to one-fifth.

Improving data archiving practices in ancient genomics

Ancient DNA is producing a rich record of past genetic diversity in humans and other species. However, unless the primary data is appropriately archived, its long-term value will not be fully realised. I surveyed publicly archived data from 42 recent ancient genomics studies. Half of the studies archived incomplete datasets, preventing accurate replication and representing a loss of data of potential future use. No studies met all criteria that could be considered best practice. Based on these results, I make six recommendations for data producers: (1) archive all sequencing reads, not just those that aligned to a reference genome, (2) archive read alignments too, but as secondary analysis files, (3) provide correct experiment metadata on samples, libraries and sequencing runs, (4) provide informative sample metadata, (5) archive data from low-coverage and negative experiments, and (6) document archiving choices in papers, and peer review these. Given the reliance on destructive sampling of finite material, ancient genomics studies have a particularly strong responsibility to ensure the longevity and reusability of generated data.

Unravelling reference bias in ancient DNA datasets

MOTIVATION

The alignment of sequencing reads is a critical step in the characterization of ancient genomes. However, reference bias and spurious mappings pose a significant challenge, particularly as cutting-edge wet lab methods generate datasets that push the boundaries of alignment tools. Reference bias occurs when reference alleles are favoured over alternative alleles during mapping, whereas spurious mappings stem from either contamination or when endogenous reads fail to align to their correct position. Previous work has shown that these phenomena are correlated with read length but a more thorough investigation of reference bias and spurious mappings for ancient DNA has been lacking. Here, we use a range of empirical and simulated palaeogenomic datasets to investigate the impacts of mapping tools, quality thresholds, and reference genome on mismatch rates across read lengths.

RESULTS

For these analyses, we introduce AMBER, a new bioinformatics tool for assessing the quality of ancient DNA mapping directly from BAM-files and informing on reference bias, read length cut-offs and reference selection. AMBER rapidly and simultaneously computes the sequence read mapping bias in the form of the mismatch rates per read length, cytosine deamination profiles at both CpG and non-CpG sites, fragment length distributions, and genomic breadth and depth of coverage. Using AMBER, we find that mapping algorithms and quality threshold choices dictate reference bias and rates of spurious alignment at different read lengths in a predictable manner, suggesting that optimized mapping parameters for each read length will be a key step in alleviating reference bias and spurious mappings.

AVAILABILITY AND IMPLEMENTATION

AMBER is available for noncommercial use on GitHub (https://github.com/tvandervalk/AMBER.git). Scripts used to generate and analyse simulated datasets are available on Github (https://github.com/sdolenz/refbias_scripts).

Archaeology meets environmental genomics: implementing sedaDNA in the study of the human past

Sedimentary ancient DNA (sedaDNA) has become one of the standard applications in the field of paleogenomics in recent years. It has been used for paleoenvironmental reconstructions, detecting the presence of prehistoric species in the absence of macro remains and even investigating the evolutionary history of a few species. However, its application in archaeology has been limited and primarily focused on humans. This article argues that sedaDNA holds significant potential in addressing key archaeological questions concerning the origins, lifestyles, and environments of past human populations. Our aim is to facilitate the integration of sedaDNA into the standard workflows in archaeology as a transformative tool, thereby unleashing its full potential for studying the human past. Ultimately, we not only underscore the challenges inherent in the sedaDNA field but also provide a research agenda for essential enhancements needed for implementing sedaDNA into the archaeological workflow.

Utilization of qPCR to Determine Duration and Environmental Drivers Contributing to the Persistence of Human DNA in Soil

Little is known about the underlying mechanisms that contribute to the persistence and degradation of DNA within soil. The goals of this study are to determine the duration of mitochondrial DNA (mtDNA) and nuclear DNA (nuDNA) persistence in soils enriched by surface-level human decomposition and to better understand the contribution of environmental factors. The surface-level decomposition of three human cadavers was documented over 11 weeks. Based on quantitative PCR results, we found nuDNA to persist in soils six weeks post-placement, while mtDNA was recoverable for the entire 11-week decomposition period. Principle components analyses and Spearman's rank correlations revealed that (1) time, (2) total body score, and (3) weekly average air temperature were significantly correlated with concentrations of nuDNA and mtDNA in soil, suggesting these factors play a role in the degradation of DNA in soils.

Late Pleistocene stickleback environmental genomes reveal the chronology of freshwater adaptation

Directly observing the chronology and tempo of adaptation in response to ecological change is rarely possible in natural ecosystems. Sedimentary ancient DNA (sedaDNA) has been shown to be a tractable source of genome-scale data of long-dead organisms1,2,3 and to thereby potentially provide an understanding of the evolutionary histories of past populations.4,5 To date, time series of ecosystem biodiversity have been reconstructed from sedaDNA, typically using DNA metabarcoding or shotgun sequence data generated from less than 1 g of sediment.6,7 Here, we maximize sequence coverage by extracting DNA from ∼50× more sediment per sample than the majority of previous studies1,2,3 to achieve genotype resolution. From a time series of Late Pleistocene sediments spanning from a marine to freshwater ecosystem, we compare adaptive genotypes reconstructed from the environmental genomes of three-spined stickleback at key time points of this transition. We find a staggered temporal dynamic in which freshwater alleles at known loci of large effect in marine-freshwater divergence of three-spined stickleback (e.g., EDA)8 were already established during the brackish phase of the formation of the isolation basin. However, marine alleles were still detected across the majority of marine-freshwater divergence-associated loci, even after the complete isolation of the lake from marine ingression. Our retrospective approach to studying adaptation from environmental genomes of three-spined sticklebacks at the end of the last glacial period complements contemporary experimental approaches9,10,11 and highlights the untapped potential for retrospective "evolve and resequence" natural experiments using sedaDNA.

Exploring different methodological approaches to unlock paleobiodiversity in peat profiles using ancient DNA

Natural and human-induced environmental changes deeply affected terrestrial ecosystems throughout the Holocene. Paleoenvironmental reconstructions provide information about the past and allow us to predict/model future scenarios. Among potential records, peat bogs are widely used because they present a precise stratigraphy and act as natural archives of highly diverse organic remains. Over the decades, several techniques have been developed to identify debris occurring in peat, including their morphological description. However, this is strongly constrained by the researcher's ability to distinguish residues at the species level, which typically requires many years of experience. In addition, potential contamination hampers using these techniques to obtain information from organisms such as fungi or bacteria. Environmental DNA metabarcoding and shotgun metagenome sequencing could represent a solution to detect specific groups of organisms without any a priori knowledge of their characteristics and/or to identify organisms that have rarely been considered in previous investigations. Moreover, shotgun metagenomics may allow the identification of bacteria and fungi (including both yeast and filamentous life forms), ensuring discrimination between ancient and modern organisms through the study of deamination/damage patterns. In the present review, we aim to i) present the state-of-the-art methodologies in paleoecological and paleoclimatic studies focusing on peat core analyses, proposing alternative approaches to the classical morphological identification of plant residues, and ii) suggest biomolecular approaches that will allow the use of proxies such as invertebrates, fungi, and bacteria, which are rarely employed in paleoenvironmental reconstructions.

Binding interaction of environmental DNA with typical emerging perfluoroalkyl acids and its impact on bioavailability

As the replacement compounds of perfluoroalkyl acids (PFAAs), emerging PFAAs generally exhibit equal or more hazardous toxicity than legacy PFAAs. Numerous DNA as environmental organic matters coexists with emerging PFAAs, but their interactions and the resulting interaction impacts on the bioavailability of emerging PFAAs remain insufficiently understood. Here, we studied the binding strength and mechanism between DNA and emerging PFAAs (perfluorobutyric acid, perfluorobutylsulfonic acid, and hexafluoropropylene oxide dimer acid) using perfluorooctanoic acid as the control, and further investigated the impacts of DNA binding on the bioavailability of the emerging PFAAs. Isothermal titration calorimetry and quantum chemical calculation found that the emerging PFAAs could bind with DNA bases (main thymine) by van der Waals force and halogen-bond, showing the binding affinities in the range of 7.87 × 104 to L/mol to 6.54 × 106 L/mol. The PFAAs-DNA binding significantly decreased the bioavailability of the PFAAs in both seedlings and plants of pakchoi (Brassica chinensis L.), with little differences in bioavailability change extent among PFAAs. The findings highlight the universality and similarity of the DNA binding effects on PFAAs bioavailability, which can be the natural detoxification mechanism for response to the PFAAs pollution.

slendr: a framework for spatio-temporal population genomic simulations on geographic landscapes

One of the goals of population genetics is to understand how evolutionary forces shape patterns of genetic variation over time. However, because populations evolve across both time and space, most evolutionary processes also have an important spatial component, acting through phenomena such as isolation by distance, local mate choice, or uneven distribution of resources. This spatial dimension is often neglected, partly due to the lack of tools specifically designed for building and evaluating complex spatio-temporal population genetic models. To address this methodological gap, we present a new framework for simulating spatially-explicit genomic data, implemented in a new R package called slendr (www.slendr.net), which leverages a SLiM simulation back-end script bundled with the package. With this framework, the users can programmatically and visually encode spatial population ranges and their temporal dynamics (i.e., population displacements, expansions, and contractions) either on real Earth landscapes or on abstract custom maps, and schedule splits and gene-flow events between populations using a straightforward declarative language. Additionally, slendr can simulate data from traditional, non-spatial models, either with SLiM or using an alternative built-in coalescent msprime back end. Together with its R-idiomatic interface to the tskit library for tree-sequence processing and analysis, slendr opens up the possibility of performing efficient, reproducible simulations of spatio-temporal genomic data entirely within the R environment, leveraging its wealth of libraries for geospatial data analysis, statistics, and visualization. Here, we present the design of the slendr R package and demonstrate its features on several practical example workflows.

aMeta: an accurate and memory-efficient ancient metagenomic profiling workflow

Analysis of microbial data from archaeological samples is a growing field with great potential for understanding ancient environments, lifestyles, and diseases. However, high error rates have been a challenge in ancient metagenomics, and the availability of computational frameworks that meet the demands of the field is limited. Here, we propose aMeta, an accurate metagenomic profiling workflow for ancient DNA designed to minimize the amount of false discoveries and computer memory requirements. Using simulated data, we benchmark aMeta against a current state-of-the-art workflow and demonstrate its superiority in microbial detection and authentication, as well as substantially lower usage of computer memory.

A glimpse of the paleome in endolithic microbial communities

BACKGROUND

The terrestrial subsurface is home to a significant proportion of the Earth's microbial biomass. Our understanding about terrestrial subsurface microbiomes is almost exclusively derived from groundwater and porous sediments mainly by using 16S rRNA gene surveys. To obtain more insights about biomass of consolidated rocks and the metabolic status of endolithic microbiomes, we investigated interbedded limestone and mudstone from the vadose zone, fractured aquifers, and deep aquitards.

RESULTS

By adapting methods from microbial archaeology and paleogenomics, we could recover sufficient DNA for downstream metagenomic analysis from seven rock specimens independent of porosity, lithology, and depth. Based on the extracted DNA, we estimated between 2.81 and 4.25 × 105 cells × g-1 rock. Analyzing DNA damage patterns revealed paleome signatures (genetic records of past microbial communities) for three rock specimens, all obtained from the vadose zone. DNA obtained from deep aquitards isolated from surface input was not affected by DNA decay indicating that water saturation and not flow is controlling subsurface microbial survival. Decoding the taxonomy and functional potential of paleome communities revealed increased abundances for sequences affiliated with chemolithoautotrophs and taxa such as Cand. Rokubacteria. We also found a broader metabolic potential in terms of aromatic hydrocarbon breakdown, suggesting a preferred utilization of sedimentary organic matter in the past.

CONCLUSIONS

Our study suggests that limestones function as archives for genetic records of past microbial communities including those sensitive to environmental stress at modern times, due to their specific conditions facilitating long-term DNA preservation. Video Abstract.

Neanderthal coexistence with Homo sapiens in Europe was affected by herbivore carrying capacity

It has been proposed that climate change and the arrival of modern humans in Europe affected the disappearance of Neanderthals due to their impact on trophic resources; however, it has remained challenging to quantify the effect of these factors. By using Bayesian age models to derive the chronology of the European Middle to Upper Paleolithic transition, followed by a dynamic vegetation model that provides the Net Primary Productivity, and a macroecological model to compute herbivore abundance, we show that in continental regions where the ecosystem productivity was low or unstable, Neanderthals disappeared before or just after the arrival of Homo sapiens. In contrast, regions with high and stable productivity witnessed a prolonged coexistence between both species. The temporal overlap between Neanderthals and H. sapiens is significantly correlated with the carrying capacity of small- and medium-sized herbivores. These results suggest that herbivore abundance released the trophic pressure of the secondary consumers guild, which affected the coexistence likelihood between both human species.

Climate shifts orchestrated hominin interbreeding events across Eurasia

When, where, and how often hominin interbreeding happened is largely unknown. We study the potential for Neanderthal-Denisovan admixture using species distribution models that integrate extensive fossil, archaeological, and genetic data with transient coupled general circulation model simulations of global climate and biomes. Our Pleistocene hindcast of past hominins' habitat suitability reveals pronounced climate-driven zonal shifts in the main overlap region of Denisovans and Neanderthals in central Eurasia. These shifts, which influenced the timing and intensity of potential interbreeding events, can be attributed to the response of climate and vegetation to past variations in atmospheric carbon dioxide and Northern Hemisphere ice-sheet volume. Therefore, glacial-interglacial climate swings likely played an important role in favoring gene flow between archaic humans.

Releasing secrets bound to ancient remains with modern DNA extraction techniques: an interview with Elena Essel

Elena Essel (Msc) spoke to Ebony Torrington, Managing Editor of BioTechniques. Essel is a molecular biologist in Matthias Meyer's Advanced DNA Sequencing Techniques group at the Max Planck Institute for Evolutionary Anthropology in Leipzig (Germany). Essel studied biology at University of Erlangen-Nuremberg (Erlangen, Germany) for her bachelor's and in Martin-Luther-University Halle-Wittenberg (Halle an der Saale, Germany) for her master's. Essel worked in Meyer's group on DNA extraction of very degraded material for her master's thesis. Meyer is an expert in developing new cutting-edge methods for researching ancient DNA, with a focus on skeletal remains, and more recently on sediment remains. Essel now focusses on DNA sampling and extraction aspects of the pipeline at Meyer's lab for the ancient DNA workflow.

Prospecting for Zoonotic Pathogens by Using Targeted DNA Enrichment

More than 60 zoonoses are linked to small mammals, including some of the most devastating pathogens in human history. Millions of museum-archived tissues are available to understand natural history of those pathogens. Our goal was to maximize the value of museum collections for pathogen-based research by using targeted sequence capture. We generated a probe panel that includes 39,916 80-bp RNA probes targeting 32 pathogen groups, including bacteria, helminths, fungi, and protozoans. Laboratory-generated, mock-control samples showed that we are capable of enriching targeted loci from pathogen DNA 2,882‒6,746-fold. We identified bacterial species in museum-archived samples, including Bartonella, a known human zoonosis. These results showed that probe-based enrichment of pathogens is a highly customizable and efficient method for identifying pathogens from museum-archived tissues.

Ancient human DNA recovered from a Palaeolithic pendant

Artefacts made from stones, bones and teeth are fundamental to our understanding of human subsistence strategies, behaviour and culture in the Pleistocene. Although these resources are plentiful, it is impossible to associate artefacts to specific human individuals¹ who can be morphologically or genetically characterized, unless they are found within burials, which are rare in this time period. Thus, our ability to discern the societal roles of Pleistocene individuals based on their biological sex or genetic ancestry is limited^2-5. Here we report the development of a non-destructive method for the gradual release of DNA trapped in ancient bone and tooth artefacts. Application of the method to an Upper Palaeolithic deer tooth pendant from Denisova Cave, Russia, resulted in the recovery of ancient human and deer mitochondrial genomes, which allowed us to estimate the age of the pendant at approximately 19,000-25,000 years. Nuclear DNA analysis identifies the presumed maker or wearer of the pendant as a female individual with strong genetic affinities to a group of Ancient North Eurasian individuals who lived around the same time but were previously found only further east in Siberia. Our work redefines how cultural and genetic records can be linked in prehistoric archaeology.

Endocranial anatomy of the Guercy 1 early Neanderthal from Baume Moula-Guercy (Soyons, Ardèche, France)

We provide the first comparative description of the endocranium of the Guercy 1 Early Neanderthal and examine its affinities to Preneanderthals, Neanderthals, and Homo sapiens. The Guercy 1 cranium derives from deposits chronostratigraphically and biostratigraphically dated to the Eemian Interglacial (MIS 5e). For comparative purposes, we compiled a sample of European and Southwest Asian subadult and adult Middle-to-Late Pleistocene hominins (≈MIS 12-MIS 1; N = 65). We sampled both a Preneanderthal-Neanderthal group and a Homo sapiens group. The Preneanderthal-Neanderthal group was further divided into three time-successive subgroups defined by associated MIS stages. Metric and morphological observations were made on original fossils and physical and virtual endocranial reconstructions. Guercy 1 and other Early Neanderthals, differ from Preneanderthals by increased development of the prefrontal cortex, precentral and postcentral gyri, inferior parietal lobule, and frontoparietal operculum. Early Neanderthal differ, in general, from Late Neanderthals by exhibiting less development in most of the latter brain structures. The late group additionally differentiates itself from the early group by a greater development of the rostral superior parietal lobule, angular gyrus, superior and middle temporal gyri, and caudal branches of the superior temporal gyrus. Endocranial morphology assessed along the Preneanderthal-Neanderthal sequence show that brain structures prominent in Preneanderthals are accentuated in Early-to-Late Neanderthals. However, both the Early and Late groups differentiate themselves by also showing regionally specific changes in brain development. This pattern of morphological change is consistent with a mosaic pattern of neural evolution in these Middle-to-Late Pleistocene hominins.

Inference of population structure and admixture proportion from Y chromosomal data of Chinese population

In the past two decades, Y chromosome data has been generated for human population genetic studies. These Y chromosome datasets were produced with various testing methods and markers, thus difficult to combine them for a comprehensive analysis. In this study, we combine four human Y chromosomal datasets of Han, Tibetan, Hui, and Li ethnic groups. The dataset contains 27 microsatellites and 137 single nucleotide polymorphisms these populations share in common. We assembled a single dataset containing 2439 individuals from 25 nationwide populations in China. A systematic analysis of genetic distance and clustering was performed. To determine the gene flow of the studied population with worldwide populations, we modeled the ancestry informative markers. The reference panel was regarded as a mixture of South Asian (SAS), East Asian (EAS), European (EUR), African (AFR), and American (AMR) populations from 1000 Genomes data of Y chromosome using nonlinear data-fitting. We then calculated the admixture proportion of these four studied populations with 26 worldwide populations. The results showed that the Han and Hui have great genetic affinity, and Hui is the most admixed ethnic group, with 61.53% EAS, 34.65% SAS, 1.91% AFR, 1.56% AMR, and 0.04% EUR ancestry component (the AMR is highly admixed and thus should be ignored). All the other three ethnic groups contained more than 97% EAS ancestry component. The Li is the least admixed population in this study. The combined dataset in this study is the largest of this kind reported to date and proposes reference population data for use in future paternal genetic studies and forensic genealogical identification.

A Case of Chronic Maxillary Sinusitis in a Late Neanderthal Population of the Altai Mountains

We describe a likely case of chronic maxillary sinusitis (CMS) in a Neanderthal skeletal sample from Chagyrskaya Cave, in the Altai Mountains. Signs of CMS were recorded in the Chagyrskaya 57 specimen, which is a fragment of a left maxilla. Alveoli of the upper fi rst molar are partially preserved, and so are the second and third upper molars, with adjacent parts of the walls, and the fl oor of the maxillary sinus. The fragment was found in layer 6b, dating to 53,100–51,100 BP. We analyze the factors that had caused the development of the disease, and assess its etiology. In the 3D-model, generated by computed microtomography, and in the original specimen, porotic changes were registered, situated at the fracture line of the alveoli of M1, lost post-mortem, and near the vestibular roots of both preserved molars. Also, there were isolated bone spicules, 1.0–2.6 mm in size. These signs indicate incipient CMS, evidently caused by chronic periodontal disease combined with a deep alveolar recess of the maxillary sinus. As the periodontal gap expanded, several small nutrient foramina, piercing the bottom of the sinus, merged. As a result, several oro-antral channels formed, whereupon the infection spread into the maxillary sinus. Since the deep alveolar recess is observed in the vast majority of Neanderthal crania with published images or reconstructed maxillary cavities, it can be assumed that Neanderthals were predisposed to odontogenic CMS.

Genetic insights into the social organization of Neanderthals

Genomic analyses of Neanderthals have previously provided insights into their population history and relationship to modern humans1-8, but the social organization of Neanderthal communities remains poorly understood. Here we present genetic data for 13 Neanderthals from two Middle Palaeolithic sites in the Altai Mountains of southern Siberia: 11 from Chagyrskaya Cave9,10 and 2 from Okladnikov Cave11-making this one of the largest genetic studies of a Neanderthal population to date. We used hybridization capture to obtain genome-wide nuclear data, as well as mitochondrial and Y-chromosome sequences. Some Chagyrskaya individuals were closely related, including a father-daughter pair and a pair of second-degree relatives, indicating that at least some of the individuals lived at the same time. Up to one-third of these individuals' genomes had long segments of homozygosity, suggesting that the Chagyrskaya Neanderthals were part of a small community. In addition, the Y-chromosome diversity is an order of magnitude lower than the mitochondrial diversity, a pattern that we found is best explained by female migration between communities. Thus, the genetic data presented here provide a detailed documentation of the social organization of an isolated Neanderthal community at the easternmost extent of their known range.

The contribution of Neanderthal introgression to modern human traits

Neanderthals, our closest extinct relatives, lived in western Eurasia from 400,000 years ago until they went extinct around 40,000 years ago. DNA retrieved from ancient specimens revealed that Neanderthals mated with modern human contemporaries. As a consequence, introgressed Neanderthal DNA survives scattered across the human genome such that 1-4% of the genome of present-day people outside Africa are inherited from Neanderthal ancestors. Patterns of Neanderthal introgressed genomic sequences suggest that Neanderthal alleles had distinct fates in the modern human genetic background. Some Neanderthal alleles facilitated human adaptation to new environments such as novel climate conditions, UV exposure levels and pathogens, while others had deleterious consequences. Here, we review the body of work on Neanderthal introgression over the past decade. We describe how evolutionary forces shaped the genomic landscape of Neanderthal introgression and highlight the impact of introgressed alleles on human biology and phenotypic variation.

Extended longevity of DNA preservation in Levantine Paleolithic sediments, Sefunim Cave, Israel

Paleogenomic research can elucidate the evolutionary history of human and faunal populations. Although the Levant is a key land-bridge between Africa and Eurasia, thus far, relatively little ancient DNA data has been generated from this region, since DNA degrades faster in warm climates. As sediments can be a source of ancient DNA, we analyzed 33 sediment samples from different sedimentological contexts in the Paleolithic layers of Sefunim Cave (Israel). Four contained traces of ancient Cervidae and Hyaenidae mitochondrial DNA. Dating by optical luminescence and radiocarbon indicates that the DNA comes from layers between 30,000 and 70,000 years old, surpassing theoretical expectations regarding the longevity of DNA deposited in such a warm environment. Both identified taxa are present in the zooarchaeological record of the site but have since gone extinct from the region, and a geoarchaeological study suggests little movement of the sediments after their deposition, lending further support to our findings. We provide details on the local conditions in the cave, which we hypothesize were particularly conducive to the long-term preservation of DNA-information that will be pertinent for future endeavors aimed at recovering ancient DNA from the Levant and other similarly challenging contexts.

Genomic analysis of a novel Neanderthal from Mezmaiskaya Cave provides insights into the genetic relationships of Middle Palaeolithic populations

The Mezmaiskaya cave is located on the North Caucasus near the border that divides Europe and Asia. Previously, fossil remains for two Neanderthals were reported from Mezmaiskaya Cave. A tooth from the third archaic hominin specimen (Mezmaiskaya 3) was retrieved from layer 3 in Mezmaiskaya Cave. We performed genome sequencing of Mezmaiskaya 3. Analysis of partial nuclear genome sequence revealed that it belongs to a Homo sapiens neanderthalensis female. Based on a high-coverage mitochondrial genome sequence, we demonstrated that the relationships of Mezmaiskaya 3 to Mezmaiskaya 1 and Stajnia S5000 individuals were closer than those to other Neanderthals. Our data demonstrate the close genetic connections between the early Middle Palaeolithic Neanderthals that were replaced by genetically distant later group in the same geographic areas. Based on mitochondrial DNA (mtDNA) data, we suggest that Mezmaiskaya 3 was the latest Neanderthal individual from the early Neanderthal’s branches. We proposed a hierarchical nomenclature for the mtDNA haplogroups of Neanderthals. In addition, we retrieved ancestral mtDNA mutations in presumably functional sites fixed in the Neanderthal clades, and also provided the first data showing mtDNA heteroplasmy in Neanderthal specimen.

Evolving ancient DNA techniques and the future of human history

Ancient DNA (aDNA) techniques applied to human genomics have significantly advanced in the past decade, enabling large-scale aDNA research, sometimes independent of human remains. This commentary reviews the major milestones of aDNA techniques and explores future directions to expand the scope of aDNA research and insights into present-day human health.

Predicting Archaic Hominin Phenotypes from Genomic Data

Ancient DNA provides a powerful window into the biology of extant and extinct species, including humans' closest relatives: Denisovans and Neanderthals. Here, we review what is known about archaic hominin phenotypes from genomic data and how those inferences have been made. We contend that understanding the influence of variants on lower-level molecular phenotypes-such as gene expression and protein function-is a promising approach to using ancient DNA to learn about archaic hominin traits. Molecular phenotypes have simpler genetic architectures than organism-level complex phenotypes, and this approach enables moving beyond association studies by proposing hypotheses about the effects of archaic variants that are testable in model systems. The major challenge to understanding archaic hominin phenotypes is broadening our ability to accurately map genotypes to phenotypes, but ongoing advances ensure that there will be much more to learn about archaic hominin phenotypes from their genomes. Expected final online publication date for the Annual Review of Genomics and Human Genetics, Volume 23 is October 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

The intrusive nature of the Châtelperronian in the Iberian Peninsula

Multiple factors have been proposed to explain the disappearance of Neandertals between ca. 50 and 40 kyr BP. Central to these discussions has been the identification of new techno-cultural complexes that overlap with the period of Neandertal demise in Europe. One such complex is the Châtelperronian, which extends from the Paris Basin to the Northern Iberian Peninsula between 43,760–39,220 BP. In this study we present the first open-air Châtelperronian site in the Northern Iberian Peninsula, Aranbaltza II. The technological features of its stone tool assemblage show no links with previous Middle Paleolithic technology in the region, and chronological modeling reveals a gap between the latest Middle Paleolithic and the Châtelperronian in this area. We interpret this as evidence of local Neandertal extinction and replacement by other Neandertal groups coming from southern France, illustrating how local extinction episodes could have played a role in the process of disappearance of Neandertals.

Quantifying and reducing cross-contamination in single- and multiplex hybridization capture of ancient DNA

The use of hybridization capture has enabled a massive upscaling in sample sizes for ancient DNA studies, allowing the analysis of hundreds of skeletal remains (Mathieson et al., 2015; Narasimhan et al., 2019) or sediments (Vernot et al., 2021; Wang et al., 2021; Zavala et al., 2021) in single studies. Nevertheless, demands in throughput continue to grow, and hybridization capture has become a limiting step in sample preparation due to the large consumption of reagents, consumables and time. Here we explore the possibility of improving the economics of sample preparation via multiplex capture, i.e. the hybridization capture of pools of double-indexed ancient DNA libraries. We demonstrate that this strategy is feasible, at least for small genomic targets such as mitochondrial DNA, if the annealing temperature is increased and PCR cycles are limited in post-capture amplification to avoid index swapping by jumping PCR, which manifests as cross-contamination in resulting sequence data. We also show that the re-amplification of double-indexed libraries to PCR plateau before or after hybridization capture can sporadically lead to small, but detectable cross-contamination even if libraries are amplified in separate reactions. We provide protocols for both manual capture and automated capture in 384-well format that are compatible with single- and multiplex capture and effectively suppress cross-contamination and artefact formation. Last, we provide a simple computational method for quantifying cross-contamination due to index swapping in double-indexed libraries, which we recommend using for routine quality checks in studies that are sensitive to cross-contamination.

Genomics of Ancient Pathogens: First Advances and Prospects

Paleogenomics is one of the urgent and promising areas of interdisciplinary research in the today's world science. New genomic methods of ancient DNA (aDNA) analysis, such as next generation sequencing (NGS) technologies, make it possible not only to obtain detailed genetic information about historical and prehistoric human populations, but also to study individual microbial and viral pathogens and microbiomes from different ancient and historical objects. Studies of aDNA of pathogens by reconstructing their genomes have so far yielded complete sequences of the ancient pathogens that played significant role in the history of the world: Yersinia pestis (plague), Variola virus (smallpox), Vibrio cholerae (cholera), HBV (hepatitis B virus), as well as the equally important endemic human infectious agents: Mycobacterium tuberculosis (tuberculosis), Mycobacterium leprae (leprosy), and Treponema pallidum (syphilis). Genomic data from these pathogens complemented the information previously obtained by paleopathologists and allowed not only to identify pathogens from the past pandemics, but also to recognize the pathogen lineages that are now extinct, to refine chronology of the pathogen appearance in human populations, and to reconstruct evolutionary history of the pathogens that are still relevant to public health today. In this review, we describe state-of-the-art genomic research of the origins and evolution of many ancient pathogens and viruses and examine mechanisms of the emergence and spread of the ancient infections in the mankind history.

Ancient DNA Methods Improve Forensic DNA Profiling of Korean War and World War II Unknowns

The integration of massively parallel sequencing (MPS) technology into forensic casework has been of particular benefit to the identification of unknown military service members. However, highly degraded or chemically treated skeletal remains often fail to provide usable DNA profiles, even with sensitive mitochondrial (mt) DNA capture and MPS methods. In parallel, the ancient DNA field has developed workflows specifically for degraded DNA, resulting in the successful recovery of nuclear DNA and mtDNA from skeletal remains as well as sediment over 100,000 years old. In this study we use a set of disinterred skeletal remains from the Korean War and World War II to test if ancient DNA extraction and library preparation methods improve forensic DNA profiling. We identified an ancient DNA extraction protocol that resulted in the recovery of significantly more human mtDNA fragments than protocols previously used in casework. In addition, utilizing single-stranded rather than double-stranded library preparation resulted in increased attainment of reportable mtDNA profiles. This study emphasizes that the combination of ancient DNA extraction and library preparation methods evaluated here increases the success rate of DNA profiling, and likelihood of identifying historical remains.

Microstratigraphic preservation of ancient faunal and hominin DNA in Pleistocene cave sediments

DNA preserved in sediments has emerged as an important source of information about past ecosystems, independent of the discovery of skeletal remains. However, little is known about the sources of sediment DNA, the factors affecting its long-term preservation, and the extent to which it may be translocated after deposition. Here, we show that impregnated blocks of intact sediment are excellent archives of DNA. DNA distribution is highly heterogeneous at the microscale in the cave sediment we studied, suggesting that postdepositional movement of DNA is unlikely to be a common phenomenon in cases where the stratigraphy is undisturbed. Combining micromorphological analysis with microstratigraphic retrieval of ancient DNA therefore allows genetic information to be associated with the detailed archaeological and ecological record preserved in sediments.

Ancient DNA recovered from Pleistocene sediments represents a rich resource for the study of past hominin and environmental diversity. However, little is known about how DNA is preserved in sediments and the extent to which it may be translocated between archaeological strata. Here, we investigate DNA preservation in 47 blocks of resin-impregnated archaeological sediment collected over the last four decades for micromorphological analyses at 13 prehistoric sites in Europe, Asia, Africa, and North America and show that such blocks can preserve DNA of hominins and other mammals. Extensive microsampling of sediment blocks from Denisova Cave in the Altai Mountains reveals that the taxonomic composition of mammalian DNA differs drastically at the millimeter-scale and that DNA is concentrated in small particles, especially in fragments of bone and feces (coprolites), suggesting that these are substantial sources of DNA in sediments. Three microsamples taken in close proximity in one of the blocks yielded Neanderthal DNA from at least two male individuals closely related to Denisova 5, a Neanderthal toe bone previously recovered from the same layer. Our work indicates that DNA can remain stably localized in sediments over time and provides a means of linking genetic information to the archaeological and ecological records on a microstratigraphic scale.

Microvolume DNA extraction methods for microscale amplicon and metagenomic studies

Investigating the composition and metabolic capacity of aquatic microbial assemblages usually requires the filtration of multi-litre samples, which are up to 1 million-fold larger than the microenvironments within which microbes are predicted to be spatially organised. To determine if community profiles can be reliably generated from microlitre volumes, we sampled seawater at a coastal and an oceanic site, filtered and homogenised them, and extracted DNA from bulk samples (2 L) and microvolumes (100, 10 and 1 μL) using two new approaches. These microvolume DNA extraction methods involve either physical or chemical lysis (through pH/thermal shock and lytic enzymes/surfactants, respectively), directly followed by the capture of DNA on magnetic beads. Downstream analysis of extracted DNA using both amplicon sequencing and metagenomics, revealed strong correlation with standard large volume approaches, demonstrating the fidelity of taxonomic and functional profiles of microbial communities in as little as 1 μL of seawater. This volume is six orders of magnitude smaller than most standard operating procedures for marine metagenomics, which will allow precise sampling of the heterogenous landscape that microbes inhabit.

Ancient DNA diffuses from human bones to cave stones

Recent studies have demonstrated the potential to recover ancient human mitochondrial DNA and nuclear DNA from cave sediments. However, the source of such sedimentary ancient DNA is still under discussion. Here we report the case of a Bronze Age human skeleton, found in a limestone cave, which was covered with layers of calcite stone deposits. By analyzing samples representing bones and stone deposits from this cave, we were able to: i) reconstruct the full human mitochondrial genome from the bones and the stones (same haplotype); ii) determine the sex of the individual; iii) reconstruct six ancient bacterial and archaeal genomes; and finally iv) demonstrate better ancient DNA preservation in the stones than in the bones. Thereby, we demonstrate the direct diffusion of human DNA from bones into the surrounding environment and show the potential to reconstruct ancient microbial genomes from such cave deposits, which represent an additional paleoarcheological archive resource.

Response to Comment on "A global environmental crisis 42,000 years ago"

[Figure: see text].

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

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

Palaeogenetics: Dirt, what is it good for? Everything

Researchers in the ancient DNA community have suspected for well over a decade that ancient whole genomes can be found in sediments. Three new studies now provide such evidence and, with it, endless possibilities for future studies of sedimentary ancient DNA.

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.

Our Tangled Family Tree: New Genomic Methods Offer Insight into the Legacy of Archaic Admixture

The archaic ancestry present in the human genome has captured the imagination of both scientists and the wider public in recent years. This excitement is the result of new studies pushing the envelope of what we can learn from the archaic genetic information that has survived for over 50,000 years in the human genome. Here, we review the most recent ten years of literature on the topic of archaic introgression, including the current state of knowledge on Neanderthal and Denisovan introgression, as well as introgression from other as-yet unidentified archaic populations. We focus this review on four topics: 1) a reimagining of human demographic history, including evidence for multiple admixture events between modern humans, Neanderthals, Denisovans, and other archaic populations; 2) state-of-the-art methods for detecting archaic ancestry in population-level genomic data; 3) how these novel methods can detect archaic introgression in modern African populations; and 4) the functional consequences of archaic gene variants, including how those variants were co-opted into novel function in modern human populations. The goal of this review is to provide a simple-to-access reference for the relevant methods and novel data, which has changed our understanding of the relationship between our species and its siblings. This body of literature reveals the large degree to which the genetic legacy of these extinct hominins has been integrated into the human populations of today.