Nadin Rohland

A High-Coverage Genome Sequence from an Archaic Denisovan Individual

Ancient Genomics The Denisovans were archaic humans closely related to Neandertals, whose populations overlapped with the ancestors of modern-day humans. Using a single-stranded library preparation method, Meyer et al. (p. 222, published online 30 August) provide a detailed analysis of a high-quality Denisovan genome. The genomic sequence provides evidence for very low rates of heterozygosity in the Denisova, probably not because of recent inbreeding, but instead because of a small population size. The genome sequence also illuminates the relationships between humans and archaics, including Neandertals, and establishes a catalog of genetic changes within the human lineage. A close-up look provides clues to the relationships between modern humans, Denisovans, and Neandertals. We present a DNA library preparation method that has allowed us to reconstruct a high-coverage (30×) genome sequence of a Denisovan, an extinct relative of Neandertals. The quality of this genome allows a direct estimation of Denisovan heterozygosity indicating that genetic diversity in these archaic hominins was extremely low. It also allows tentative dating of the specimen on the basis of “missing evolution” in its genome, detailed measurements of Denisovan and Neandertal admixture into present-day human populations, and the generation of a near-complete catalog of genetic changes that swept to high frequency in modern humans since their divergence from Denisovans.

Ancient Admixture in Human History

Population mixture is an important process in biology. We present a suite of methods for learning about population mixtures, implemented in a software package called ADMIXTOOLS, that support formal tests for whether mixture occurred and make it possible to infer proportions and dates of mixture. We also describe the development of a new single nucleotide polymorphism (SNP) array consisting of 629,433 sites with clearly documented ascertainment that was specifically designed for population genetic analyses and that we genotyped in 934 individuals from 53 diverse populations. To illustrate the methods, we give a number of examples that provide new insights about the history of human admixture. The most striking finding is a clear signal of admixture into northern Europe, with one ancestral population related to present-day Basques and Sardinians and the other related to present-day populations of northeast Asia and the Americas. This likely reflects a history of admixture between Neolithic migrants and the indigenous Mesolithic population of Europe, consistent with recent analyses of ancient bones from Sweden and the sequencing of the genome of the Tyrolean “Iceman.”

Massive migration from the steppe was a source for Indo-European languages in Europe

We generated genome-wide data from 69 Europeans who lived between 8,000–3,000 years ago by enriching ancient DNA libraries for a target set of almost 400,000 polymorphisms. Enrichment of these positions decreases the sequencing required for genome-wide ancient DNA analysis by a median of around 250-fold, allowing us to study an order of magnitude more individuals than previous studies and to obtain new insights about the past. We show that the populations of Western and Far Eastern Europe followed opposite trajectories between 8,000–5,000 years ago. At the beginning of the Neolithic period in Europe, ∼8,000–7,000 years ago, closely related groups of early farmers appeared in Germany, Hungary and Spain, different from indigenous hunter-gatherers, whereas Russia was inhabited by a distinctive population of hunter-gatherers with high affinity to a ∼24,000-year-old Siberian. By ∼6,000–5,000 years ago, farmers throughout much of Europe had more hunter-gatherer ancestry than their predecessors, but in Russia, the Yamnaya steppe herders of this time were descended not only from the preceding eastern European hunter-gatherers, but also from a population of Near Eastern ancestry. Western and Eastern Europe came into contact ∼4,500 years ago, as the Late Neolithic Corded Ware people from Germany traced ∼75% of their ancestry to the Yamnaya, documenting a massive migration into the heartland of Europe from its eastern periphery. This steppe ancestry persisted in all sampled central Europeans until at least ∼3,000 years ago, and is ubiquitous in present-day Europeans. These results provide support for a steppe origin of at least some of the Indo-European languages of Europe.

Genomically recoded organisms expand biological functions.

Changing the Code Easily and efficiently expanding the genetic code could provide tools to genome engineers with broad applications in medicine, energy, agriculture, and environmental safety. Lajoie et al. (p. 357) replaced all known UAG stop codons with synonymous UAA stop codons in Escherichia coli MG1655, as well as release factor 1 (RF1; terminates translation at UAG), thereby eliminating natural UAG translation function without impairing fitness. This made it possible to reassign UAG as a dedicated codon to genetically encode nonstandard amino acids while avoiding deleterious incorporation at native UAG positions. The engineered E. coli incorporated nonstandard amino acids into its proteins and showed enhanced resistance to bacteriophage T7. In a second paper, Lajoie et al. (p. 361) demonstrated the recoding of 13 codons in 42 highly expressed essential genes in E. coli. Codon usage was malleable, but synonymous codons occasionally were nonequivalent in unpredictable ways. Bacteria engineered to use nonstandard amino acids show increased resistance to bacteriophage attack. We describe the construction and characterization of a genomically recoded organism (GRO). We replaced all known UAG stop codons in Escherichia coli MG1655 with synonymous UAA codons, which permitted the deletion of release factor 1 and reassignment of UAG translation function. This GRO exhibited improved properties for incorporation of nonstandard amino acids that expand the chemical diversity of proteins in vivo. The GRO also exhibited increased resistance to T7 bacteriophage, demonstrating that new genetic codes could enable increased viral resistance.

Ancient DNA extraction from bones and teeth

This method is designed to maximize recovery of PCR-amplifiable DNA from ancient bone and teeth specimens and at the same time to minimize co-extraction of substances that inhibit PCR. This is achieved by a combination of DNA extraction from bone powder using a buffer consisting solely of EDTA and proteinase K, and purification of the DNA by binding to silica in the presence of high concentrations of guanidinium thiocyanate. All steps are performed at room temperature (20–23 °C), thereby reducing further degradation of the already damaged and fragile ancient DNA and providing an optimal trade-off between DNA release and degradation. Furthermore, the purification step removes most of the various types of PCR inhibitors present in ancient bone samples, thereby optimizing the amount of ancient DNA available for subsequent enzymatic manipulation, such as PCR amplification. The protocol presented here allows DNA extraction from ancient bone and teeth with a minimum of working steps and equipment and yields DNA extracts within 2 working days.

Cost-effective, high-throughput DNA sequencing libraries for multiplexed target capture

Improvements in technology have reduced the cost of DNA sequencing to the point that the limiting factor for many experiments is the time and reagent cost of sample preparation. We present an approach in which 192 sequencing libraries can be produced in a single day of technician time at a cost of about

Genome-wide patterns of selection in 230 ancient Eurasians

15 per sample. These libraries are effective not only for low-pass whole-genome sequencing, but also for simultaneously enriching them in pools of approximately 100 individually barcoded samples for a subset of the genome without substantial loss in efficiency of target capture. We illustrate the power and effectiveness of this approach on about 2000 samples from a prostate cancer study.

Comparison and optimization of ancient DNA extraction

Ancient DNA makes it possible to observe natural selection directly by analysing samples from populations before, during and after adaptation events. Here we report a genome-wide scan for selection using ancient DNA, capitalizing on the largest ancient DNA data set yet assembled: 230 West Eurasians who lived between 6500 and 300 bc, including 163 with newly reported data. The new samples include, to our knowledge, the first genome-wide ancient DNA from Anatolian Neolithic farmers, whose genetic material we obtained by extracting from petrous bones, and who we show were members of the population that was the source of Europe’s first farmers. We also report a transect of the steppe region in Samara between 5600 and 300 bc, which allows us to identify admixture into the steppe from at least two external sources. We detect selection at loci associated with diet, pigmentation and immunity, and two independent episodes of selection on height.

A melanocortin 1 receptor allele suggests varying pigmentation among Neanderthals

Ancient DNA analyses rely on the extraction of the tiny amounts of DNA remaining in samples that are hundreds to tens of thousands of years old. Despite the critical role extraction efficiency plays in this field of research, no study has comprehensively compared ancient DNA extraction techniques to date. There are a wide range of methods currently in use, which rely on such disparate principles as spin columns, alcohol precipitation, or binding to silica. We have compared a number of these methods using quantitative PCR and then optimized each step of the most promising method. We found that most chemicals routinely added to ancient DNA extraction buffers do not increase, and sometimes even decrease, DNA yields. Consequently, our optimized method uses a buffer consisting solely of EDTA and proteinase K for bone digestion and binding DNA to silica via guanidinium thiocyanate for DNA purification. In a comparison with published methods, this minimalist approach, on average, outperforms all other methods in terms of DNA yields as measured using quantitative PCR. We also found that the addition of bovine serum albumin (BSA) to the PCR helps to overcome inhibitors in ancient DNA extracts. Finally, we observed a marked difference in the performance between different types of DNA polymerases, as measured by amplification success.

The landscape of recombination in African Americans

The melanocortin 1 receptor (MC1R) regulates pigmentation in humans and other vertebrates. Variants of MC1R with reduced function are associated with pale skin color and red hair in humans of primarily European origin. We amplified and sequenced a fragment of the MC1R gene (mc1r) from two Neanderthal remains. Both specimens have a mutation that was not found in ∼3700 modern humans analyzed. Functional analyses show that this variant reduces MC1R activity to a level that alters hair and/or skin pigmentation in humans. The impaired activity of this variant suggests that Neanderthals varied in pigmentation levels, potentially on the scale observed in modern humans. Our data suggest that inactive MC1R variants evolved independently in both modern humans and Neanderthals.