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.

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.

Developmental fate and cellular maturity encoded in human regulatory DNA landscapes.

Cellular-state information between generations of developing cells may be propagated via regulatory regions. We report consistent patterns of gain and loss of DNase I-hypersensitive sites (DHSs) as cells progress from embryonic stem cells (ESCs) to terminal fates. DHS patterns alone convey rich information about cell fate and lineage relationships distinct from information conveyed by gene expression. Developing cells share a proportion of their DHS landscapes with ESCs; that proportion decreases continuously in each cell type as differentiation progresses, providing a quantitative benchmark of developmental maturity. Developmentally stable DHSs densely encode binding sites for transcription factors involved in autoregulatory feedback circuits. In contrast to normal cells, cancer cells extensively reactivate silenced ESC DHSs and those from developmental programs external to the cell lineage from which the malignancy derives. Our results point to changes in regulatory DNA landscapes as quantitative indicators of cell-fate transitions, lineage relationships, and dysfunction.

Evolutionary history and adaptation from high-coverage whole-genome sequences of diverse African hunter-gatherers.

To reconstruct modern human evolutionary history and identify loci that have shaped hunter-gatherer adaptation, we sequenced the whole genomes of five individuals in each of three different hunter-gatherer populations at > 60× coverage: Pygmies from Cameroon and Khoesan-speaking Hadza and Sandawe from Tanzania. We identify 13.4 million variants, substantially increasing the set of known human variation. We found evidence of archaic introgression in all three populations, and the distribution of time to most recent common ancestors from these regions is similar to that observed for introgressed regions in Europeans. Additionally, we identify numerous loci that harbor signatures of local adaptation, including genes involved in immunity, metabolism, olfactory and taste perception, reproduction, and wound healing. Within the Pygmy population, we identify multiple highly differentiated loci that play a role in growth and anterior pituitary function and are associated with height.

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.

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.

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.

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.

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.

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.

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.

Complex history of admixture between modern humans and Neandertals.

Recent analyses have found that a substantial amount of the Neandertal genome persists in the genomes of contemporary non-African individuals. East Asians have, on average, higher levels of Neandertal ancestry than do Europeans, which might be due to differences in the efficiency of purifying selection, an additional pulse of introgression into East Asians, or other unexplored scenarios. To better define the scope of plausible models of archaic admixture between Neandertals and anatomically modern humans, we analyzed patterns of introgressed sequence in whole-genome data of 379 Europeans and 286 East Asians. We found that inferences of demographic history restricted to neutrally evolving genomic regions allowed a simple one-pulse model to be robustly rejected, suggesting that differences in selection cannot explain the differences in Neandertal ancestry. We show that two additional demographic models, involving either a second pulse of Neandertal gene flow into the ancestors of East Asians or a dilution of Neandertal lineages in Europeans by admixture with an unknown ancestral population, are consistent with the data. Thus, the history of admixture between modern humans and Neandertals is most likely more complex than previously thought.

Sporadic autism exomes reveal a highly interconnected protein network of de novo mutations.

It is well established that autism spectrum disorders (ASD) have a strong genetic component; however, for at least 70% of cases, the underlying genetic cause is unknown. Under the hypothesis that de novo mutations underlie a substantial fraction of the risk for developing ASD in families with no previous history of ASD or related phenotypes--so-called sporadic or simplex families--we sequenced all coding regions of the genome (the exome) for parent-child trios exhibiting sporadic ASD, including 189 new trios and 20 that were previously reported. Additionally, we also sequenced the exomes of 50 unaffected siblings corresponding to these new (n = 31) and previously reported trios (n = 19), for a total of 677 individual exomes from 209 families. Here we show that de novo point mutations are overwhelmingly paternal in origin (4:1 bias) and positively correlated with paternal age, consistent with the modest increased risk for children of older fathers to develop ASD. Moreover, 39% (49 of 126) of the most severe or disruptive de novo mutations map to a highly interconnected ß-catenin/chromatin remodelling protein network ranked significantly for autism candidate genes. In proband exomes, recurrent protein-altering mutations were observed in two genes: CHD8 and NTNG1. Mutation screening of six candidate genes in 1,703 ASD probands identified additional de novo, protein-altering mutations in GRIN2B, LAMC3 and SCN1A. Combined with copy number variant (CNV) data, these results indicate extreme locus heterogeneity but also provide a target for future discovery, diagnostics and therapeutics.

The accessible chromatin landscape of the human genome.

DNase I hypersensitive sites (DHSs) are markers of regulatory DNA and have underpinned the discovery of all classes of cis-regulatory elements including enhancers, promoters, insulators, silencers and locus control regions. Here we present the first extensive map of human DHSs identified through genome-wide profiling in 125 diverse cell and tissue types. We identify ∼2.9 million DHSs that encompass virtually all known experimentally validated cis-regulatory sequences and expose a vast trove of novel elements, most with highly cell-selective regulation. Annotating these elements using ENCODE data reveals novel relationships between chromatin accessibility, transcription, DNA methylation and regulatory factor occupancy patterns. We connect ∼580,000 distal DHSs with their target promoters, revealing systematic pairing of different classes of distal DHSs and specific promoter types. Patterning of chromatin accessibility at many regulatory regions is organized with dozens to hundreds of co-activated elements, and the transcellular DNase I sensitivity pattern at a given region can predict cell-type-specific functional behaviours. The DHS landscape shows signatures of recent functional evolutionary constraint. However, the DHS compartment in pluripotent and immortalized cells exhibits higher mutation rates than that in highly differentiated cells, exposing an unexpected link between chromatin accessibility, proliferative potential and patterns of human variation.

An expansive human regulatory lexicon encoded in transcription factor footprints.

Regulatory factor binding to genomic DNA protects the underlying sequence from cleavage by DNase I, leaving nucleotide-resolution footprints. Using genomic DNase I footprinting across 41 diverse cell and tissue types, we detected 45 million transcription factor occupancy events within regulatory regions, representing differential binding to 8.4 million distinct short sequence elements. Here we show that this small genomic sequence compartment, roughly twice the size of the exome, encodes an expansive repertoire of conserved recognition sequences for DNA-binding proteins that nearly doubles the size of the human cis-regulatory lexicon. We find that genetic variants affecting allelic chromatin states are concentrated in footprints, and that these elements are preferentially sheltered from DNA methylation. High-resolution DNase I cleavage patterns mirror nucleotide-level evolutionary conservation and track the crystallographic topography of protein-DNA interfaces, indicating that transcription factor structure has been evolutionarily imprinted on the human genome sequence. We identify a stereotyped 50-base-pair footprint that precisely defines the site of transcript origination within thousands of human promoters. Finally, we describe a large collection of novel regulatory factor recognition motifs that are highly conserved in both sequence and function, and exhibit cell-selective occupancy patterns that closely parallel major regulators of development, differentiation and pluripotency.

Personal and population genomics of human regulatory variation.

The characteristics and evolutionary forces acting on regulatory variation in humans remains elusive because of the difficulty in defining functionally important noncoding DNA. Here, we combine genome-scale maps of regulatory DNA marked by DNase I hypersensitive sites (DHSs) from 138 cell and tissue types with whole-genome sequences of 53 geographically diverse individuals in order to better delimit the patterns of regulatory variation in humans. We estimate that individuals likely harbor many more functionally important variants in regulatory DNA compared with protein-coding regions, although they are likely to have, on average, smaller effect sizes. Moreover, we demonstrate that there is significant heterogeneity in the level of functional constraint in regulatory DNA among different cell types. We also find marked variability in functional constraint among transcription factor motifs in regulatory DNA, with sequence motifs for major developmental regulators, such as HOX proteins, exhibiting levels of constraint comparable to protein-coding regions. Finally, we perform a genome-wide scan of recent positive selection and identify hundreds of novel substrates of adaptive regulatory evolution that are enriched for biologically interesting pathways such as melanogenesis and adipocytokine signaling. These data and results provide new insights into patterns of regulatory variation in individuals and populations and demonstrate that a large proportion of functionally important variation lies beyond the exome.

Inferring duplications, losses, transfers and incomplete lineage sorting with nonbinary species trees.

MOTIVATION: Gene duplication (D), transfer (T), loss (L) and incomplete lineage sorting (I) are crucial to the evolution of gene families and the emergence of novel functions. The history of these events can be inferred via comparison of gene and species trees, a process called reconciliation, yet current reconciliation algorithms model only a subset of these evolutionary processes. RESULTS: We present an algorithm to reconcile a binary gene tree with a nonbinary species tree under a DTLI parsimony criterion. This is the first reconciliation algorithm to capture all four evolutionary processes driving tree incongruence and the first to reconcile non-binary species trees with a transfer model. Our algorithm infers all optimal solutions and reports complete, temporally feasible event histories, giving the gene and species lineages in which each event occurred. It is fixed-parameter tractable, with polytime complexity when the maximum species outdegree is fixed. Application of our algorithms to prokaryotic and eukaryotic data show that use of an incomplete event model has substantial impact on the events inferred and resulting biological conclusions. AVAILABILITY: Our algorithms have been implemented in Notung, a freely available phylogenetic reconciliation software package, available at http://www.cs.cmu.edu/~durand/Notung. CONTACT: mstolzer@andrew.cmu.edu.

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.

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.

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.