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 individuals1 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 limited2-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.

Characterization of RNA content in individual phase-separated coacervate microdroplets.

Condensates formed by complex coacervation are hypothesized to have played a crucial part during the origin-of-life. In living cells, condensation organizes biomolecules into a wide range of membraneless compartments. Although RNA is a key component of biological condensates and the central component of the RNA world hypothesis, little is known about what determines RNA accumulation in condensates and to which extend single condensates differ in their RNA composition. To address this, we developed an approach to read the RNA content from single synthetic and protein-based condensates using high-throughput sequencing. We find that certain RNAs efficiently accumulate in condensates. These RNAs are strongly enriched in sequence motifs which show high sequence similarity to short interspersed elements (SINEs). We observe similar results for protein-derived condensates, demonstrating applicability across different in vitro reconstituted membraneless organelles. Thus, our results provide a new inroad to explore the RNA content of phase-separated droplets at single condensate resolution.

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 or sediments 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 explored the possibility of improving the economics of sample preparation via multiplex capture, that is, 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 reamplification 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.

Microstratigraphic preservation of ancient faunal and hominin DNA in Pleistocene cave 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.

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 Neanderthals1-4. The dozen hominin remains recovered from the deposits also include Neanderthals5,6 and the child of a Neanderthal and a Denisovan7, 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 assemblages5,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.

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.

Initial Upper Palaeolithic humans in Europe had recent Neanderthal ancestry.

Modern humans appeared in Europe by at least 45,000 years ago1-5, but the extent of their interactions with Neanderthals, who disappeared by about 40,000 years ago6, and their relationship to the broader expansion of modern humans outside Africa are poorly understood. Here we present genome-wide data from three individuals dated to between 45,930 and 42,580 years ago from Bacho Kiro Cave, Bulgaria1,2. They are the earliest Late Pleistocene modern humans known to have been recovered in Europe so far, and were found in association with an Initial Upper Palaeolithic artefact assemblage. Unlike two previously studied individuals of similar ages from Romania7 and Siberia8 who did not contribute detectably to later populations, these individuals are more closely related to present-day and ancient populations in East Asia and the Americas than to later west Eurasian populations. This indicates that they belonged to a modern human migration into Europe that was not previously known from the genetic record, and provides evidence that there was at least some continuity between the earliest modern humans in Europe and later people in Eurasia. Moreover, we find that all three individuals had Neanderthal ancestors a few generations back in their family history, confirming that the first European modern humans mixed with Neanderthals and suggesting that such mixing could have been common.

The evolutionary history of Neanderthal and Denisovan Y chromosomes.

Ancient DNA has provided new insights into many aspects of human history. However, we lack comprehensive studies of the Y chromosomes of Denisovans and Neanderthals because the majority of specimens that have been sequenced to sufficient coverage are female. Sequencing Y chromosomes from two Denisovans and three Neanderthals shows that the Y chromosomes of Denisovans split around 700 thousand years ago from a lineage shared by Neanderthals and modern human Y chromosomes, which diverged from each other around 370 thousand years ago. The phylogenetic relationships of archaic and modern human Y chromosomes differ from the population relationships inferred from the autosomal genomes and mirror mitochondrial DNA phylogenies, indicating replacement of both the mitochondrial and Y chromosomal gene pools in late Neanderthals. This replacement is plausible if the low effective population size of Neanderthals resulted in an increased genetic load in Neanderthals relative to modern humans.

Human Stem Cell Resources Are an Inroad to Neandertal DNA Functions.

Induced pluripotent stem cells (iPSCs) from diverse humans offer the potential to study human functional variation in controlled culture environments. A portion of this variation originates from an ancient admixture between modern humans and Neandertals, which introduced alleles that left a phenotypic legacy on individual humans today. Here, we show that a large iPSC repository harbors extensive Neandertal DNA, including alleles that contribute to human phenotypes and diseases, encode hundreds of amino acid changes, and alter gene expression in specific tissues. We provide a database of the inferred introgressed Neandertal alleles for each individual iPSC line, together with the annotation of the predicted functional variants. We also show that transcriptomic data from organoids generated from iPSCs can be used to track Neandertal-derived RNA over developmental processes. Human iPSC resources provide an opportunity to experimentally explore Neandertal DNA function and its contribution to present-day phenotypes, and potentially study Neandertal traits.

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.

Limits of long-term selection against Neandertal introgression.

Several studies have suggested that introgressed Neandertal DNA was subjected to negative selection in modern humans. A striking observation in support of this is an apparent monotonic decline in Neandertal ancestry observed in modern humans in Europe over the past 45,000 years. Here, we show that this decline is an artifact likely caused by gene flow between modern human populations, which is not taken into account by statistics previously used to estimate Neandertal ancestry. When we apply a statistic that avoids assumptions about modern human demography by taking advantage of two high-coverage Neandertal genomes, we find no evidence for a change in Neandertal ancestry in Europe over the past 45,000 years. We use whole-genome simulations of selection and introgression to investigate a wide range of model parameters and find that negative selection is not expected to cause a significant long-term decline in genome-wide Neandertal ancestry. Nevertheless, these models recapitulate previously observed signals of selection against Neandertal alleles, in particular the depletion of Neandertal ancestry in conserved genomic regions. Surprisingly, we find that this depletion is strongest in regulatory and conserved noncoding regions and in the most conserved portion of protein-coding sequences.

admixr-R package for reproducible analyses using ADMIXTOOLS.

SUMMARY: We present a new R package admixr, which provides a convenient interface for performing reproducible population genetic analyses (f3, D, f4, f4-ratio, qpWave and qpAdm), as implemented by command-line programs in the ADMIXTOOLS software suite. In a traditional ADMIXTOOLS workflow, the user must first generate a set of text configuration files tailored to each individual analysis, often using a combination of shell scripting and manual text editing. The non-tabular output files then need to be parsed to extract values of interest prior to further analyses. Our package simplifies this process by automating all low-level configuration and parsing steps, making analyses as simple as running a single R command. Furthermore, we provide a set of R functions for processing, filtering and manipulating datasets in the EIGENSTRAT format. By unifying all steps of the workflow under a single R framework, this package enables the automation of analytic pipelines, significantly improving the reproducibility of population genetic studies. AVAILABILITY AND IMPLEMENTATION: The source code of the R package is available under the MIT license. Installation instructions, reference manual and a tutorial can be found on the package website at https://bioinf.eva.mpg.de/admixr. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Neandertal Introgression Sheds Light on Modern Human Endocranial Globularity.

One of the features that distinguishes modern humans from our extinct relatives and ancestors is a globular shape of the braincase [1-4]. As the endocranium closely mirrors the outer shape of the brain, these differences might reflect altered neural architecture [4, 5]. However, in the absence of fossil brain tissue, the underlying neuroanatomical changes as well as their genetic bases remain elusive. To better understand the biological foundations of modern human endocranial shape, we turn to our closest extinct relatives: the Neandertals. Interbreeding between modern humans and Neandertals has resulted in introgressed fragments of Neandertal DNA in the genomes of present-day non-Africans [6, 7]. Based on shape analyses of fossil skull endocasts, we derive a measure of endocranial globularity from structural MRI scans of thousands of modern humans and study the effects of introgressed fragments of Neandertal DNA on this phenotype. We find that Neandertal alleles on chromosomes 1 and 18 are associated with reduced endocranial globularity. These alleles influence expression of two nearby genes, UBR4 and PHLPP1, which are involved in neurogenesis and myelination, respectively. Our findings show how integration of fossil skull data with archaic genomics and neuroimaging can suggest developmental mechanisms that may contribute to the unique modern human endocranial shape.

The genome of the offspring of a Neanderthal mother and a Denisovan father.

Neanderthals and Denisovans are extinct groups of hominins that separated from each other more than 390,000 years ago1,2. Here we present the genome of 'Denisova 11', a bone fragment from Denisova Cave (Russia)3 and show that it comes from an individual who had a Neanderthal mother and a Denisovan father. The father, whose genome bears traces of Neanderthal ancestry, came from a population related to a later Denisovan found in the cave4-6. The mother came from a population more closely related to Neanderthals who lived later in Europe2,7 than to an earlier Neanderthal found in Denisova Cave8, suggesting that migrations of Neanderthals between eastern and western Eurasia occurred sometime after 120,000 years ago. The finding of a first-generation Neanderthal-Denisovan offspring among the small number of archaic specimens sequenced to date suggests that mixing between Late Pleistocene hominin groups was common when they met.

Evolutionary history and adaptation of a human pygmy population of Flores Island, Indonesia.

Flores Island, Indonesia, was inhabited by the small-bodied hominin species Homo floresiensis, which has an unknown evolutionary relationship to modern humans. This island is also home to an extant human pygmy population. Here we describe genome-scale single-nucleotide polymorphism data and whole-genome sequences from a contemporary human pygmy population living on Flores near the cave where H. floresiensis was found. The genomes of Flores pygmies reveal a complex history of admixture with Denisovans and Neanderthals but no evidence for gene flow with other archaic hominins. Modern individuals bear the signatures of recent positive selection encompassing the FADS (fatty acid desaturase) gene cluster, likely related to diet, and polygenic selection acting on standing variation that contributed to their short-stature phenotype. Thus, multiple independent instances of hominin insular dwarfism occurred on Flores.

The Predecessors Within . . .

By examining the genomes of present-day people from Asia, researchers show that modern humans met and interbred with Denisovans, distant relatives to Neanderthals, on at least two occasions. As a result, people today carry DNA from two different Denisovan populations.

A high-coverage Neandertal genome from Vindija Cave in Croatia.

To date, the only Neandertal genome that has been sequenced to high quality is from an individual found in Southern Siberia. We sequenced the genome of a female Neandertal from ~50,000 years ago from Vindija Cave, Croatia, to ~30-fold genomic coverage. She carried 1.6 differences per 10,000 base pairs between the two copies of her genome, fewer than present-day humans, suggesting that Neandertal populations were of small size. Our analyses indicate that she was more closely related to the Neandertals that mixed with the ancestors of present-day humans living outside of sub-Saharan Africa than the previously sequenced Neandertal from Siberia, allowing 10 to 20% more Neandertal DNA to be identified in present-day humans, including variants involved in low-density lipoprotein cholesterol concentrations, schizophrenia, and other diseases.

Archaic Hominin Admixture Facilitated Adaptation to Out-of-Africa Environments.

As modern humans dispersed from Africa throughout the world, they encountered and interbred with archaic hominins, including Neanderthals and Denisovans [1, 2]. Although genome-scale maps of introgressed sequences have been constructed [3-6], considerable gaps in knowledge remain about the functional, phenotypic, and evolutionary significance of archaic hominin DNA that persists in present-day individuals. Here, we describe a comprehensive set of analyses that identified 126 high-frequency archaic haplotypes as putative targets of adaptive introgression in geographically diverse populations. These loci are enriched for immune-related genes (such as OAS1/2/3, TLR1/6/10, and TNFAIP3) and also encompass genes (including OCA2 and BNC2) that influence skin pigmentation phenotypes. Furthermore, we leveraged existing and novel large-scale gene expression datasets to show many positively selected archaic haplotypes act as expression quantitative trait loci (eQTLs), suggesting that modulation of transcript abundance was a common mechanism facilitating adaptive introgression. Our results demonstrate that hybridization between modern and archaic hominins provided an important reservoir of advantageous alleles that enabled adaptation to out-of-Africa environments.

Differences and similarities between human and chimpanzee neural progenitors during cerebral cortex development.

Human neocortex expansion likely contributed to the remarkable cognitive abilities of humans. This expansion is thought to primarily reflect differences in proliferation versus differentiation of neural progenitors during cortical development. Here, we have searched for such differences by analysing cerebral organoids from human and chimpanzees using immunohistofluorescence, live imaging, and single-cell transcriptomics. We find that the cytoarchitecture, cell type composition, and neurogenic gene expression programs of humans and chimpanzees are remarkably similar. Notably, however, live imaging of apical progenitor mitosis uncovered a lengthening of prometaphase-metaphase in humans compared to chimpanzees that is specific to proliferating progenitors and not observed in non-neural cells. Consistent with this, the small set of genes more highly expressed in human apical progenitors points to increased proliferative capacity, and the proportion of neurogenic basal progenitors is lower in humans. These subtle differences in cortical progenitors between humans and chimpanzees may have consequences for human neocortex evolution.