Nuclear DNA from two early Neandertals reveals 80,000 years of genetic continuity in Europe

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.

Functional characterization of archaic-specific variants in mitonuclear genes: insights from comparative analysis in S. cerevisiae

Neanderthal and Denisovan hybridisation with modern humans has generated a non-random genomic distribution of introgressed regions, the result of drift and selection dynamics. Cross-species genomic incompatibility and more efficient removal of slightly deleterious archaic variants have been proposed as selection-based processes involved in the post-hybridisation purge of archaic introgressed regions. Both scenarios require the presence of functionally different alleles across Homo species onto which selection operated differently according to which populations hosted them, but only a few of these variants have been pinpointed so far. In order to identify functionally divergent archaic variants removed in humans, we focused on mitonuclear genes, which are underrepresented in the genomic landscape of archaic humans. We searched for non-synonymous, fixed, archaic-derived variants present in mitonuclear genes, rare or absent in human populations. We then compared the functional impact of archaic and human variants in the model organism Saccharomyces cerevisiae. Notably, a variant within the mitochondrial tyrosyl-tRNA synthetase 2 (YARS2) gene exhibited a significant decrease in respiratory activity and a substantial reduction of Cox2 levels, a proxy for mitochondrial protein biosynthesis, coupled with the accumulation of the YARS2 protein precursor and a lower amount of mature enzyme. Our work suggests that this variant is associated with mitochondrial functionality impairment, thus contributing to the purging of archaic introgression in YARS2. While different molecular mechanisms may have impacted other mitonuclear genes, our approach can be extended to the functional screening of mitonuclear genetic variants present across species and populations.

The Neanderthal niche space of Western Eurasia 145 ka to 30 ka ago

Neanderthals occupied Western Eurasia between 350 ka and 40 ka ago, during the climatically volatile Pleistocene. A key issue is to what extent Neanderthal populations expanded into areas of Western Eurasia and what conditions facilitated such range expansions. The range extent of Neanderthals is generally based on the distribution of Neanderthal material, but the land-altering nature of glacial periods has erased much of the already sparse material evidence of Neanderthals, particularly in the northern latitudes. To overcome this obstacle species distribution models can estimate past distributions of Neanderthals, however, most implementations are generally constrained spatially and temporally and may be artificially truncating the Neanderthal niche space. Using dated contexts from Neanderthal sites from across Western Eurasia, millennial-scale paleoclimate reconstructions, and a spatiotemporal species distribution model, we infer the fundamental climatic niche space of Neanderthals and estimate the extent of Neanderthal occupation. We find that (a.) despite the long timeframe, Neanderthals occupy a relatively narrow fundamental climatic niche space, (b.) the estimated projected potential Neanderthal niche space suggests a larger geographic range than the material record suggests, and (c.) that there was a general decline in the size of the projected potential Neanderthal niche from 145 ka ago onward, possibly contributing to their extinction.

Sex estimation of the adult Neandertal Regourdou 1 (Montignac, France): Implications for sexing human fossil remains

Sex is a biological trait fundamental to the study of hominin fossils. Among the many questions that can be addressed are those related to taxonomy, biological variability, sexual dimorphism, paleoobstetrics, funerary selection, and paleodemography. While new methodologies such as paleogenomics or paleoproteomics can be used to determine sex, they have not been systematically applied to Pleistocene human remains due to their destructive nature. Therefore, we estimated sex from the coxal bone of the newly discovered pelvic remains of the Regourdou 1 Neandertal (Southwest France, MIS 5) based on morphological and metric data employing two methods that have been recently revised and shown to be reliable in multiple studies. Both methods calculate posterior probabilities of the estimate. The right coxal bone of Regourdou 1 was partially reconstructed providing additional traits for sex estimation. These methods were cross validated on 14 sufficiently preserved coxal bones of specimens from the Neandertal lineage. Our results show that the Regourdou 1 individual, whose postcranial skeleton is not robust, is a male, and that previous sex attributions of comparative Neandertal specimens are largely in agreement with those obtained here. Our results encourage additional morphological research of fossil hominins in order to develop a set of methods that are applicable, reliable, and reproducible.

More than a decade of genetic research on the Denisovans

Denisovans, a group of now extinct humans who lived in Eastern Eurasia in the Middle and Late Pleistocene, were first identified from DNA sequences just over a decade ago. Only ten fragmentary remains from two sites have been attributed to Denisovans based entirely on molecular information. Nevertheless, there has been great interest in using genetic data to understand Denisovans and their place in human history. From the reconstruction of a single high-quality genome, it has been possible to infer their population history, including events of admixture with other human groups. Additionally, the identification of Denisovan DNA in the genomes of present-day individuals has provided insights into the timing and routes of dispersal of ancient modern humans into Asia and Oceania, as well as the contributions of archaic DNA to the physiology of present-day people. In this Review, we synthesize more than a decade of research on Denisovans, reconcile controversies and summarize insights into their population history and phenotype. We also highlight how our growing knowledge about Denisovans has provided insights into our own evolutionary history.

Deep-time paleogenomics and the limits of DNA survival

Although most ancient DNA studies have focused on the last 50,000 years, paleogenomic approaches can now reach into the early Pleistocene, an epoch of repeated environmental changes that shaped present-day biodiversity. Emerging deep-time genomic transects, including from DNA preserved in sediments, will enable inference of adaptive evolution, discovery of unrecognized species, and exploration of how glaciations, volcanism, and paleomagnetic reversals shaped demography and community composition. In this Review, we explore the state-of-the-art in paleogenomics and discuss key challenges, including technical limitations, evolutionary divergence and associated biases, and the need for more precise dating of remains and sediments. We conclude that with improvements in laboratory and computational methods, the emerging field of deep-time paleogenomics will expand the range of questions addressable using ancient DNA.

Anatomically modern human in the Châtelperronian hominin collection from the Grotte du Renne (Arcy-sur-Cure, Northeast France)

Around 42,000 years ago, anatomically modern humans appeared in Western Europe to the detriment of indigenous Neanderthal groups. It is during this period that new techno-cultural complexes appear, such as the Châtelperronian that extends from northern Spain to the Paris Basin. The Grotte du Renne (Arcy-sur-Cure) is a key site for discussing the biological identity of its makers. This deposit has yielded several Neanderthal human remains in its Châtelperronian levels. However, the last inventory of the paleoanthropological collection attributed to this techno-complex allowed the identification of an ilium belonging to a neonate (AR-63) whose morphology required a thorough analysis to assess its taxonomic attribution. Using geometric morphometrics, we quantified its morphology and compared it to that of 2 Neanderthals and 32 recent individuals deceased during the perinatal period to explore their morphological variation. Our results indicate a morphological distinction between the ilia of Neanderthals and anatomically modern neonates. Although AR-63 is slightly outside recent variability, it clearly differs from the Neanderthals. We propose that this is due to its belonging to an early modern human lineage whose morphology differs slightly from present-day humans. We also explore different hypotheses about the presence of this anatomically modern neonate ilium among Neanderthal remains.

Nutrigenomics in the context of evolution

Nutrigenomics describes the interaction between nutrients and our genome. Since the origin of our species most of these nutrient-gene communication pathways have not changed. However, our genome experienced over the past 50,000 years a number of evolutionary pressures, which are based on the migration to new environments concerning geography and climate, the transition from hunter-gatherers to farmers including the zoonotic transfer of many pathogenic microbes and the rather recent change of societies to a preferentially sedentary lifestyle and the dominance of Western diet. Human populations responded to these challenges not only by specific anthropometric adaptations, such as skin color and body stature, but also through diversity in dietary intake and different resistance to complex diseases like the metabolic syndrome, cancer and immune disorders. The genetic basis of this adaptation process has been investigated by whole genome genotyping and sequencing including that of DNA extracted from ancient bones. In addition to genomic changes, also the programming of epigenomes in pre- and postnatal phases of life has an important contribution to the response to environmental changes. Thus, insight into the variation of our (epi)genome in the context of our individual's risk for developing complex diseases, helps to understand the evolutionary basis how and why we become ill. This review will discuss the relation of diet, modern environment and our (epi)genome including aspects of redox biology. This has numerous implications for the interpretation of the risks for disease and their prevention.

Modelling Neanderthals' dispersal routes from Caucasus towards east

The study of the cultural materials associated with the Neanderthal physical remains from the sites in the Caucasus, Central Asia and Siberian Altai and adjacent areas documents two distinct techno-complexes of Micoquian and Mousterian. These findings potentially outline two dispersal routes for the Neanderthals out of Europe. Using data on topography and Palaeoclimate, we generated computer-based least-cost-path modelling for the Neanderthal dispersal routes from Caucasus towards the east. In this regard, two dispersal routes have been identified: A northern route from Greater Caucasus associated with Micoquian techno-complex towards Siberian Altai and a southern route from Lesser Caucasus associated with Mousterian towards Siberian Altai via the Southern Caspian Corridor. Based on archaeological, bio- and physio-geographical data, our model hypothesises that during climatic deterioration phases (e.g. MIS 4) the connection between Greater and Lesser Caucasus was limited. This issue perhaps resulted in the separate development and spread of two cultural groups of Micoquian and Mousterian with an input from two different population sources of Neanderthal influxes: eastern and southern Europe refugia for these two northern and southern dispersal routes respectively. Of these two, we focus on the southern dispersal route, for it comprises a 'rapid dispersal route' towards east. The significant location of the Southern Caspian corridor between high mountains of Alborz and the Caspian Sea, provided a special biogeographical zone and a refugium. This exceptional physio-geographic condition brings forward the Southern Caspian corridor as a potential place of admixture of different hominin species including Neanderthals and homo sapiens.

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.

Oxytocin and vasotocin receptor variation and the evolution of human prosociality

Modern human lifestyle strongly depends on complex social traits like empathy, tolerance and cooperation. These diverse facets of social cognition have been associated with variation in the oxytocin receptor (OTR) and its sister genes, the vasotocin/vasopressin receptors (VTR1A/AVPR1A and AVPR1B/VTR1B). Here, we compared the available genomic sequences of these receptors between modern humans, archaic humans, and 12 non-human primate species, and identified sites that show heterozygous variation in modern humans and archaic humans distinct from variation in other primates, and for which we could find association studies with clinical implications. On these sites, we performed a range of analyses (variant clustering, pathogenicity prediction, regulation, linkage disequilibrium frequency), and reviewed the literature on selection data in different modern-human populations. We found five sites with modern human specific variation, where the modern human allele is the major allele in the global population (OTR: rs1042778, rs237885, rs6770632; VTR1A: rs10877969; VTR1B: rs33985287). Among them, variation in the OTR-rs6770632 site was predicted to be the most functional. Two alleles (OTR: rs59190448 and rs237888) present only in modern humans and archaic humans were putatively under positive selection in modern humans, with rs237888 predicted to be a highly functional site. Three sites showed convergent evolution between modern humans and bonobos (OTR: rs2228485 and rs237897; VTR1A: rs1042615), with OTR-rs2228485 ranking highly in terms of functionality and reported to be under balancing selection in modern humans (Schaschl, 2015) [1]. Our findings have implications for understanding hominid prosociality, as well as the similarities between modern human and bonobo social behavior.

•

We compared the oxytocin/vasotocin receptors in modern humans, archaic humans, and non-human primates.

•

We found 5 sites with modern human specific variation.

•

In these sites, the modern human allele is the major allele in the global population.

•

Several sites were predicted to be functional and with selection signatures in modern humans.

•

We also identified 3 sites of convergent evolution in modern humans and bonobos.

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.

Vitamin D in the Context of Evolution

For at least 1.2 billion years, eukaryotes have been able to synthesize sterols and, therefore, can produce vitamin D when exposed to UV-B. Vitamin D endocrinology was established some 550 million years ago in animals, when the high-affinity nuclear receptor VDR (vitamin D receptor), transport proteins and enzymes for vitamin D metabolism evolved. This enabled vitamin D to regulate, via its target genes, physiological process, the first of which were detoxification and energy metabolism. In this way, vitamin D was enabled to modulate the energy-consuming processes of the innate immune system in its fight against microbes. In the evolving adaptive immune system, vitamin D started to act as a negative regulator of growth, which prevents overboarding reactions of T cells in the context of autoimmune diseases. When, some 400 million years ago, species left the ocean and were exposed to gravitation, vitamin D endocrinology took over the additional role as a major regulator of calcium homeostasis, being important for a stable skeleton. Homo sapiens evolved approximately 300,000 years ago in East Africa and had adapted vitamin D endocrinology to the intensive exposure of the equatorial sun. However, when some 75,000 years ago, when anatomically modern humans started to populate all continents, they also reached regions with seasonally low or no UV-B, i.e., and under these conditions vitamin D became a vitamin.

The evolutionary history of human spindle genes includes back-and-forth gene flow with Neandertals

Proteins associated with the spindle apparatus, a cytoskeletal structure that ensures the proper segregation of chromosomes during cell division, experienced an unusual number of amino acid substitutions in modern humans after the split from the ancestors of Neandertals and Denisovans. Here, we analyze the history of these substitutions and show that some of the genes in which they occur may have been targets of positive selection. We also find that the two changes in the kinetochore scaffold 1 (KNL1) protein, previously believed to be specific to modern humans, were present in some Neandertals. We show that the KNL1 gene of these Neandertals shared a common ancestor with present-day Africans about 200,000 years ago due to gene flow from the ancestors (or relatives) of modern humans into Neandertals. Subsequently, some non-Africans inherited this modern human-like gene variant from Neandertals, but none inherited the ancestral gene variants. These results add to the growing evidence of early contacts between modern humans and archaic groups in Eurasia and illustrate the intricate relationships among these groups.

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 earliest Denisovans and their cultural adaptation

Since the initial identification of the Denisovans a decade ago, only a handful of their physical remains have been discovered. Here we analysed ~3,800 non-diagnostic bone fragments using collagen peptide mass fingerprinting to locate new hominin remains from Denisova Cave (Siberia, Russia). We identified five new hominin bones, four of which contained sufficient DNA for mitochondrial analysis. Three carry mitochondrial DNA of the Denisovan type and one was found to carry mtDNA of the Neanderthal type. The former come from the same archaeological layer near the base of the cave's sequence and are the oldest securely dated evidence of Denisovans at 200 ka (thousand years ago) (205-192 ka at 68.2% or 217-187 ka at 95% probability). The stratigraphic context in which they were located contains a wealth of archaeological material in the form of lithics and faunal remains, allowing us to determine the material culture associated with these early hominins and explore their behavioural and environmental adaptations. The combination of bone collagen fingerprinting and genetic analyses has so far more-than-doubled the number of hominin bones at Denisova Cave and has expanded our understanding of Denisovan and Neanderthal interactions, as well as their archaeological signatures.

Paleogenomics of the prehistory of Europe: human migrations, domestication and disease

A substantial portion of ancient DNA research has been centred on understanding European populations' origin and evolution. A rchaeological evidence has already shown that the peopling of Europe involved an intricate pattern of demic and/or cultural diffusion since the Upper Palaeolithic, which became more evident during the Neolithic and Bronze Age periods. However, ancient DNA data has been crucial in determining if cultural changes occurred due to the movement of ideas or people. With the advent of next-generation sequencing and population-based paleogenomic research, ancient DNA studies have been directed not only at the study of continental human migrations, but also to the detailed analysis of particular archaeological sites, the processes of domestication, or the spread of disease during prehistoric times. With this vast paleogenomic effort added to a proper archaeological contextualisation of results, a deeper understanding of Europe's peopling is starting to emanate.

PyDamage: automated ancient damage identification and estimation for contigs in ancient DNA de novo assembly

DNA de novo assembly can be used to reconstruct longer stretches of DNA (contigs), including genes and even genomes, from short DNA sequencing reads. Applying this technique to metagenomic data derived from archaeological remains, such as paleofeces and dental calculus, we can investigate past microbiome functional diversity that may be absent or underrepresented in the modern microbiome gene catalogue. However, compared to modern samples, ancient samples are often burdened with environmental contamination, resulting in metagenomic datasets that represent mixtures of ancient and modern DNA. The ability to rapidly and reliably establish the authenticity and integrity of ancient samples is essential for ancient DNA studies, and the ability to distinguish between ancient and modern sequences is particularly important for ancient microbiome studies. Characteristic patterns of ancient DNA damage, namely DNA fragmentation and cytosine deamination (observed as C-to-T transitions) are typically used to authenticate ancient samples and sequences, but existing tools for inspecting and filtering aDNA damage either compute it at the read level, which leads to high data loss and lower quality when used in combination with de novo assembly, or require manual inspection, which is impractical for ancient assemblies that typically contain tens to hundreds of thousands of contigs. To address these challenges, we designed PyDamage, a robust, automated approach for aDNA damage estimation and authentication of de novo assembled aDNA. PyDamage uses a likelihood ratio based approach to discriminate between truly ancient contigs and contigs originating from modern contamination. We test PyDamage on both on simulated aDNA data and archaeological paleofeces, and we demonstrate its ability to reliably and automatically identify contigs bearing DNA damage characteristic of aDNA. Coupled with aDNA de novo assembly, Pydamage opens up new doors to explore functional diversity in ancient metagenomic datasets.

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.

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.

Systematic benchmark of ancient DNA read mapping

The current standard practice for assembling individual genomes involves mapping millions of short DNA sequences (also known as DNA 'reads') against a pre-constructed reference genome. Mapping vast amounts of short reads in a timely manner is a computationally challenging task that inevitably produces artefacts, including biases against alleles not found in the reference genome. This reference bias and other mapping artefacts are expected to be exacerbated in ancient DNA (aDNA) studies, which rely on the analysis of low quantities of damaged and very short DNA fragments (~30-80 bp). Nevertheless, the current gold-standard mapping strategies for aDNA studies have effectively remained unchanged for nearly a decade, during which time new software has emerged. In this study, we used simulated aDNA reads from three different human populations to benchmark the performance of 30 distinct mapping strategies implemented across four different read mapping software-BWA-aln, BWA-mem, NovoAlign and Bowtie2-and quantified the impact of reference bias in downstream population genetic analyses. We show that specific NovoAlign, BWA-aln and BWA-mem parameterizations achieve high mapping precision with low levels of reference bias, particularly after filtering out reads with low mapping qualities. However, unbiased NovoAlign results required the use of an IUPAC reference genome. While relevant only to aDNA projects where reference population data are available, the benefit of using an IUPAC reference demonstrates the value of incorporating population genetic information into the aDNA mapping process, echoing recent results based on graph genome representations.

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.

Distance metrics for ranked evolutionary trees

Genealogical tree modeling is essential for estimating evolutionary parameters in population genetics and phylogenetics. Recent mathematical results concerning ranked genealogies without leaf labels unlock opportunities in the analysis of evolutionary trees. In particular, comparisons between ranked genealogies facilitate the study of evolutionary processes of different organisms sampled at multiple time periods. We propose metrics on ranked tree shapes and ranked genealogies for lineages isochronously and heterochronously sampled. Our proposed tree metrics make it possible to conduct statistical analyses of ranked tree shapes and timed ranked tree shapes or ranked genealogies. Such analyses allow us to assess differences in tree distributions, quantify estimation uncertainty, and summarize tree distributions. We show the utility of our metrics via simulations and an application in infectious diseases.

Reconstructing double-stranded DNA fragments on a single-molecule level reveals patterns of degradation in ancient samples

Extensive manipulations involved in the preparation of DNA samples for sequencing have hitherto made it impossible to determine the precise structure of double-stranded DNA fragments being sequenced, such as the presence of blunt ends, single-stranded overhangs, or single-strand breaks. We here describe MatchSeq, a method that combines single-stranded DNA library preparation from diluted DNA samples with computational sequence matching, allowing the reconstruction of double-stranded DNA fragments on a single-molecule level. The application of MatchSeq to Neanderthal DNA, a particularly complex source of degraded DNA, reveals that 1- or 2-nt overhangs and blunt ends dominate the ends of ancient DNA molecules and that short gaps exist, which are predominantly caused by the loss of individual purines. We further show that deamination of cytosine to uracil occurs in both single- and double-stranded contexts close to the ends of molecules, and that single-stranded parts of DNA fragments are enriched in pyrimidines. MatchSeq provides unprecedented resolution for interrogating the structures of fragmented double-stranded DNA and can be applied to fragmented double-stranded DNA isolated from any biological source. The method relies on well-established laboratory techniques and can easily be integrated into routine data generation. This possibility is shown by the successful reconstruction of double-stranded DNA fragments from previously published single-stranded sequence data, allowing a more comprehensive characterization of the biochemical properties not only of ancient DNA but also of cell-free DNA from human blood plasma, a clinically relevant marker for the diagnosis and monitoring of disease.

Removing reference bias and improving indel calling in ancient DNA data analysis by mapping to a sequence variation graph

BACKGROUND

During the last decade, the analysis of ancient DNA (aDNA) sequence has become a powerful tool for the study of past human populations. However, the degraded nature of aDNA means that aDNA molecules are short and frequently mutated by post-mortem chemical modifications. These features decrease read mapping accuracy and increase reference bias, in which reads containing non-reference alleles are less likely to be mapped than those containing reference alleles. Alternative approaches have been developed to replace the linear reference with a variation graph which includes known alternative variants at each genetic locus. Here, we evaluate the use of variation graph software vg to avoid reference bias for aDNA and compare with existing methods.

RESULTS

We use vg to align simulated and real aDNA samples to a variation graph containing 1000 Genome Project variants and compare with the same data aligned with bwa to the human linear reference genome. Using vg leads to a balanced allelic representation at polymorphic sites, effectively removing reference bias, and more sensitive variant detection in comparison with bwa, especially for insertions and deletions (indels). Alternative approaches that use relaxed bwa parameter settings or filter bwa alignments can also reduce bias but can have lower sensitivity than vg, particularly for indels.

CONCLUSIONS

Our findings demonstrate that aligning aDNA sequences to variation graphs effectively mitigates the impact of reference bias when analyzing aDNA, while retaining mapping sensitivity and allowing detection of variation, in particular indel variation, that was previously missed.

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

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

Triangulating Neanderthal cognition: A tale of not seeing the forest for the trees

The inference of Neanderthal cognition, including their cultural and linguistic capabilities, has persisted as a fiercely debated research topic for decades. This lack of consensus is substantially based on inherent uncertainties in reconstructing prehistory out of indirect evidence as well as other methodological limitations. Further factors include systemic difficulties within interdisciplinary discourse, data artifacts, historic research biases, and the sheer scope of the relevant research. Given the degrees of freedom in interpretation ensuing from these complications, any attempt to find approximate answers to the yet unsettled pertinent discourse may not rest on single studies, but instead a careful and comprehensive interdisciplinary synthesis of findings. Triangulating Neanderthals' cognition by considering the plethora of data, diverse perspectives and aforementioned complexities present within the literature constitutes the currently most reliable pathway to tentative conclusions. While some uncertainties remain, such an approach paints the picture of an extensive shared humanity between anatomically modern humans and Neanderthals. This article is categorized under: Cognitive Biology > Evolutionary Roots of Cognition Linguistics > Evolution of Language.

New perspectives on Neanderthal dispersal and turnover from Stajnia Cave (Poland)

The Micoquian is the broadest and longest enduring cultural facies of the Late Middle Palaeolithic that spread across the periglacial and boreal environments of Europe between Eastern France, Poland, and Northern Caucasus. Here, we present new data from the archaeological record of Stajnia Cave (Poland) and the paleogenetic analysis of a Neanderthal molar S5000, found in a Micoquian context. Our results demonstrate that the mtDNA genome of Stajnia S5000 dates to MIS 5a making the tooth the oldest Neanderthal specimen from Central-Eastern Europe. Furthermore, S5000 mtDNA has the fewest number of differences to mtDNA of Mezmaiskaya 1 Neanderthal from Northern Caucasus, and is more distant from almost contemporaneous Neanderthals of Scladina and Hohlenstein-Stadel. This observation and the technological affinity between Poland and the Northern Caucasus could be the result of increased mobility of Neanderthals that changed their subsistence strategy for coping with the new low biomass environments and the increased foraging radius of gregarious animals. The Prut and Dniester rivers were probably used as the main corridors of dispersal. The persistence of the Micoquian techno-complex in South-Eastern Europe infers that this axis of mobility was also used at the beginning of MIS 3 when a Neanderthal population turnover occurred in the Northern Caucasus.

Skin colour and vitamin D: An update

Homo sapiens evolved in East Africa and had dark skin, hair, and eyes, in order to protect against deleterious consequences of intensive UV radiation at equatorial latitudes. Intensive skin pigmentation was thought to bear the risk of inefficient vitamin D3 synthesis in the skin. This initiated the hypothesis that within the past 75 000 years, in which humans migrated to higher latitudes in Asia and Europe, the need for vitamin D3 synthesis served as an evolutionary driver for skin lightening. In this review, we summarize the recent archeogenomic reconstruction of population admixture in Europe and demonstrate that skin lightening happened as late as 5000 years ago through immigration of lighter pigmented populations from western Anatolia and the Russian steppe but not primarily via evolutionary pressure for vitamin D3 synthesis. We show that variations in genes encoding for proteins being responsible for the transport, metabolism and signalling of vitamin D provide alternative mechanisms of adaptation to a life in northern latitudes without suffering from consequences of vitamin D deficiency. This includes hypotheses explaining differences in the vitamin D status and response index of European populations.

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.

A high-coverage Neandertal genome from Chagyrskaya Cave

We sequenced the genome of a Neandertal from Chagyrskaya Cave in the Altai Mountains, Russia, to 27-fold genomic coverage. We show that this Neandertal was a female and that she was more related to Neandertals in western Eurasia [Prüfer et al., Science 358, 655-658 (2017); Hajdinjak et al., Nature 555, 652-656 (2018)] than to Neandertals who lived earlier in Denisova Cave [Prüfer et al., Nature 505, 43-49 (2014)], which is located about 100 km away. About 12.9% of the Chagyrskaya genome is spanned by homozygous regions that are between 2.5 and 10 centiMorgans (cM) long. This is consistent with the fact that Siberian Neandertals lived in relatively isolated populations of less than 60 individuals. In contrast, a Neandertal from Europe, a Denisovan from the Altai Mountains, and ancient modern humans seem to have lived in populations of larger sizes. The availability of three Neandertal genomes of high quality allows a view of genetic features that were unique to Neandertals and that are likely to have been at high frequency among them. We find that genes highly expressed in the striatum in the basal ganglia of the brain carry more amino-acid-changing substitutions than genes expressed elsewhere in the brain, suggesting that the striatum may have evolved unique functions in Neandertals.

What have the revelations about Neanderthal DNA revealed about Homo sapiens?

Genetic studies have presented increasing indications about the complexity of the interactions between Homo sapiens, Neanderthals and Denisovans, during Pleistocene. The results indicate potential replacement or admixture of the groups of hominins that lived in the same region at different times. Recently, the time of separation among these hominins in relation to the Last Common Ancestor – LCA has been reasonably well established. Events of mixing with emphasis on the Neanderthal gene flow into H. sapiens outside Africa, Denisovans into H. sapiens ancestors in Oceania and continental Asia, Neanderthals into Denisovans, as well as the origin of some phenotypic features in specific populations such as the color of the skin, eyes, hair and predisposition to develop certain kinds of diseases have also been found. The current information supports the existence of both replacement and interbreeding events, and indicates the need to revise the two main explanatory models, the Multiregional and the Out-of-Africa hypotheses, about the origin and evolution of H. sapiens and its co-relatives. There is definitely no longer the possibility of justifying only one model over the other. This paper aims to provide a brief review and update on the debate around this issue, considering the advances brought about by the recent genetic as well as morphological traits analyses.

Archaic hominin introgression into modern human genomes

Ancient genomes from multiple Neanderthal and the Denisovan individuals, along with DNA sequence data from diverse contemporary human populations strongly support the prevalence of gene flow among different hominins. Recent studies now provide evidence for multiple gene flow events that leave genetic signatures in extant and ancient human populations. These events include older gene flow from an unknown hominin in Africa predating out-of-Africa migrations, and in the last 50,000-100,000 years, multiple gene flow events from Neanderthals into ancestral Eurasian human populations, and at least three distinct introgression events from a lineage close to Denisovans into ancestors of extant Southeast Asian and Oceanic populations. Some of these introgression events may have happened as late as 20,000 years before present and reshaped the way in which we think about human evolution. In this review, I aim to answer anthropologically relevant questions with regard to recent research on ancient hominin introgression in the human lineage. How have genomic data from archaic hominins changed our view of human evolution? Is there any doubt about whether introgression from ancient hominins to the ancestors of present-day humans occurred? What is the current view of human evolutionary history from the genomics perspective? What is the impact of introgression on human phenotypes?

A genetic analysis of the Gibraltar Neanderthals

The Forbes' Quarry and Devil's Tower partial crania from Gibraltar are among the first Neanderthal remains ever found. Here, we show that small amounts of ancient DNA are preserved in the petrous bones of the 2 individuals despite unfavorable climatic conditions. However, the endogenous Neanderthal DNA is present among an overwhelming excess of recent human DNA. Using improved DNA library construction methods that enrich for DNA fragments carrying deaminated cytosine residues, we were able to sequence 70 and 0.4 megabase pairs (Mbp) nuclear DNA of the Forbes' Quarry and Devil's Tower specimens, respectively, as well as large parts of the mitochondrial genome of the Forbes' Quarry individual. We confirm that the Forbes' Quarry individual was a female and the Devil's Tower individual a male. We also show that the Forbes' Quarry individual is genetically more similar to the ∼120,000-y-old Neanderthals from Scladina Cave in Belgium (Scladina I-4A) and Hohlenstein-Stadel Cave in Germany, as well as to a ∼60,000- to 70,000-y-old Neanderthal from Russia (Mezmaiskaya 1), than to a ∼49,000-y-old Neanderthal from El Sidrón (El Sidrón 1253) in northern Spain and other younger Neanderthals from Europe and western Asia. This suggests that the Forbes' Quarry fossil predates the latter Neanderthals. The preservation of archaic human DNA in the warm coastal climate of Gibraltar, close to the shores of Africa, raises hopes for the future recovery of archaic human DNA from regions in which climatic conditions are less than optimal for DNA preservation.