Philip L. F. Johnson

Genetic history of an archaic hominin group from Denisova Cave in Siberia

Using DNA extracted from a finger bone found in Denisova Cave in southern Siberia, we have sequenced the genome of an archaic hominin to about 1.9-fold coverage. This individual is from a group that shares a common origin with Neanderthals. This population was not involved in the putative gene flow from Neanderthals into Eurasians; however, the data suggest that it contributed 4–6% of its genetic material to the genomes of present-day Melanesians. We designate this hominin population ‘Denisovans’ and suggest that it may have been widespread in Asia during the Late Pleistocene epoch. A tooth found in Denisova Cave carries a mitochondrial genome highly similar to that of the finger bone. This tooth shares no derived morphological features with Neanderthals or modern humans, further indicating that Denisovans have an evolutionary history distinct from Neanderthals and modern humans.

The complete genome sequence of a Neanderthal from the Altai Mountains

We present a high-quality genome sequence of a Neanderthal woman from Siberia. We show that her parents were related at the level of half-siblings and that mating among close relatives was common among her recent ancestors. We also sequenced the genome of a Neanderthal from the Caucasus to low coverage. An analysis of the relationships and population history of available archaic genomes and 25 present-day human genomes shows that several gene flow events occurred among Neanderthals, Denisovans and early modern humans, possibly including gene flow into Denisovans from an unknown archaic group. Thus, interbreeding, albeit of low magnitude, occurred among many hominin groups in the Late Pleistocene. In addition, the high-quality Neanderthal genome allows us to establish a definitive list of substitutions that became fixed in modern humans after their separation from the ancestors of Neanderthals and Denisovans.

Patterns of damage in genomic DNA sequences from a Neandertal

High-throughput direct sequencing techniques have recently opened the possibility to sequence genomes from Pleistocene organisms. Here we analyze DNA sequences determined from a Neandertal, a mammoth, and a cave bear. We show that purines are overrepresented at positions adjacent to the breaks in the ancient DNA, suggesting that depurination has contributed to its degradation. We furthermore show that substitutions resulting from miscoding cytosine residues are vastly overrepresented in the DNA sequences and drastically clustered in the ends of the molecules, whereas other substitutions are rare. We present a model where the observed substitution patterns are used to estimate the rate of deamination of cytosine residues in single- and double-stranded portions of the DNA, the length of single-stranded ends, and the frequency of nicks. The results suggest that reliable genome sequences can be obtained from Pleistocene organisms.

A Complete Neandertal Mitochondrial Genome Sequence Determined by High-Throughput Sequencing

A complete mitochondrial (mt) genome sequence was reconstructed from a 38,000 year-old Neandertal individual with 8341 mtDNA sequences identified among 4.8 Gb of DNA generated from approximately 0.3 g of bone. Analysis of the assembled sequence unequivocally establishes that the Neandertal mtDNA falls outside the variation of extant human mtDNAs, and allows an estimate of the divergence date between the two mtDNA lineages of 660,000 +/- 140,000 years. Of the 13 proteins encoded in the mtDNA, subunit 2 of cytochrome c oxidase of the mitochondrial electron transport chain has experienced the largest number of amino acid substitutions in human ancestors since the separation from Neandertals. There is evidence that purifying selection in the Neandertal mtDNA was reduced compared with other primate lineages, suggesting that the effective population size of Neandertals was small.

Recalibrating Equus evolution using the genome sequence of an early Middle Pleistocene horse

The rich fossil record of equids has made them a model for evolutionary processes. Here we present a 1.12-times coverage draft genome from a horse bone recovered from permafrost dated to approximately 560–780 thousand years before present (kyr bp). Our data represent the oldest full genome sequence determined so far by almost an order of magnitude. For comparison, we sequenced the genome of a Late Pleistocene horse (43 kyr bp), and modern genomes of five domestic horse breeds (Equus ferus caballus), a Przewalski’s horse (E. f. przewalskii) and a donkey (E. asinus). Our analyses suggest that the Equus lineage giving rise to all contemporary horses, zebras and donkeys originated 4.0–4.5 million years before present (Myr bp), twice the conventionally accepted time to the most recent common ancestor of the genus Equus. We also find that horse population size fluctuated multiple times over the past 2 Myr, particularly during periods of severe climatic changes. We estimate that the Przewalski’s and domestic horse populations diverged 38–72 kyr bp, and find no evidence of recent admixture between the domestic horse breeds and the Przewalski’s horse investigated. This supports the contention that Przewalski’s horses represent the last surviving wild horse population. We find similar levels of genetic variation among Przewalski’s and domestic populations, indicating that the former are genetically viable and worthy of conservation efforts. We also find evidence for continuous selection on the immune system and olfaction throughout horse evolution. Finally, we identify 29 genomic regions among horse breeds that deviate from neutrality and show low levels of genetic variation compared to the Przewalski’s horse. Such regions could correspond to loci selected early during domestication.

Assuring the quality of next-generation sequencing in clinical laboratory practice

Amy S Gargis, Centers for Disease Control and Prevention Lisa Kalman, Centers for Disease Control and Prevention Meredith W Berry, SeqWright Inc David P Bick, Medical College of Wisconsin David P Dimmock, Medical College of Wisconsin Tina Hambuch, Illumina Clinical Services Fei Lu, SeqWright Inc Elaine Lyon, University of Utah Karl V Voelkerding, University of Utah Barbara Zehnbauer, Emory University

mapDamage2.0: fast approximate Bayesian estimates of ancient DNA damage parameters

Motivation: Ancient DNA (aDNA) molecules in fossilized bones and teeth, coprolites, sediments, mummified specimens and museum collections represent fantastic sources of information for evolutionary biologists, revealing the agents of past epidemics and the dynamics of past populations. However, the analysis of aDNA generally faces two major issues. Firstly, sequences consist of a mixture of endogenous and various exogenous backgrounds, mostly microbial. Secondly, high nucleotide misincorporation rates can be observed as a result of severe post-mortem DNA damage. Such misincorporation patterns are instrumental to authenticate ancient sequences versus modern contaminants. We recently developed the user-friendly mapDamage package that identifies such patterns from next-generation sequencing (NGS) sequence datasets. The absence of formal statistical modeling of the DNA damage process, however, precluded rigorous quantitative comparisons across samples. Results: Here, we describe mapDamage 2.0 that extends the original features of mapDamage by incorporating a statistical model of DNA damage. Assuming that damage events depend only on sequencing position and post-mortem deamination, our Bayesian statistical framework provides estimates of four key features of aDNA molecules: the average length of overhangs (λ), nick frequency (ν) and cytosine deamination rates in both double-stranded regions () and overhangs (). Our model enables rescaling base quality scores according to their probability of being damaged. mapDamage 2.0 handles NGS datasets with ease and is compatible with a wide range of DNA library protocols. Availability: mapDamage 2.0 is available at ginolhac.github.io/mapDamage/ as a Python package and documentation is maintained at the Centre for GeoGenetics Web site (geogenetics.ku.dk/publications/mapdamage2.0/). Contact: jonsson.hakon@gmail.com Supplementary information: Supplementary data are available at Bioinformatics online.

A Revised Timescale for Human Evolution Based on Ancient Mitochondrial Genomes

BACKGROUND Recent analyses of de novo DNA mutations in modern humans have suggested a nuclear substitution rate that is approximately half that of previous estimates based on fossil calibration. This result has led to suggestions that major events in human evolution occurred far earlier than previously thought. RESULTS Here, we use mitochondrial genome sequences from ten securely dated ancient modern humans spanning 40,000 years as calibration points for the mitochondrial clock, thus yielding a direct estimate of the mitochondrial substitution rate. Our clock yields mitochondrial divergence times that are in agreement with earlier estimates based on calibration points derived from either fossils or archaeological material. In particular, our results imply a separation of non-Africans from the most closely related sub-Saharan African mitochondrial DNAs (haplogroup L3) that occurred less than 62-95 kya. CONCLUSIONS Though single loci like mitochondrial DNA (mtDNA) can only provide biased estimates of population divergence times, they can provide valid upper bounds. Our results exclude most of the older dates for African and non-African population divergences recently suggested by de novo mutation rate estimates in the nuclear genome.

Anatomic and functional topography of the dorsal raphe nucleus.

Abstract: Serotonergic systems play an important and generalized role in regulation of sleep‐wake states and behavioral arousal. Recent in vivo electrophysiologic recording studies in animals suggest that several different subtypes of serotonergic neurons with unique behavioral correlates exist within the brainstem raphe nuclei, raising the possibility that topographically organized subpopulations of serotonergic neurons may have unique behavioral or physiologic correlates and unique functional properties. We have shown that the stress‐related and anxiogenic neuropeptide corticotropin‐releasing factor can stimulate the in vitro neuronal firing rates of topographically organized subpopulations of serotonergic neurons within the dorsal raphe nucleus (DR). These findings are consistent with a wealth of behavioral studies suggesting that serotonergic systems within the DR are involved in the modulation of ongoing anxiety‐related behavior and in behavioral sensitization, a process whereby anxiety‐ and fear‐related behavioral responses are sensitized for a period of up to 24 to 48 h. The dorsomedial subdivision of the DR, particularly its middle and caudal aspects, has attracted considerable attention as a region that may play a critical role in the regulation of acute and chronic anxiety states. Future studies aimed at characterization of the molecular and cellular properties of topographically organized subpopulations of serotonergic neurons are likely to lead to major advances in our understanding of the role of serotonergic systems in stress‐related physiology and behavior.

Serotonergic systems associated with arousal and vigilance behaviors following administration of anxiogenic drugs

Serotonergic systems play important roles in modulating behavioral arousal, including behavioral arousal and vigilance associated with anxiety states. To further our understanding of the neural systems associated with increases in anxiety states, we investigated the effects of multiple anxiogenic drugs on topographically organized subpopulations of serotonergic neurons using double immunohistochemical staining for c-Fos and tryptophan hydroxylase combined with topographical analysis of the rat dorsal raphe nucleus (DR). Anxiogenic drugs with diverse pharmacological properties including the adenosine receptor antagonist caffeine, the serotonin 5-HT2A/2C receptor agonist m-chlorophenyl piperazine (mCPP), the alpha2-adrenoreceptor antagonist yohimbine, and the benzodiazepine receptor partial inverse agonist N-methyl-beta-carboline-3-carboxamide (FG-7142) induced increases in behavioral arousal and vigilance behaviors consistent with an increase in anxiety state. In addition, these anxiogenic drugs, excluding yohimbine, had convergent actions on an anatomically-defined subset of serotonergic neurons within the middle and caudal, dorsal subdivision of the DR. High resolution topographical analysis revealed that at the mid-rostrocaudal level, caffeine and FG-7142 had convergent effects on c-Fos expression in serotonergic neurons that were restricted to a previously undefined region, which we have named the shell region of the dorsal part of the dorsal raphe nucleus (DRDSh), that overlaps the anatomical border between the dorsal part of the dorsal raphe nucleus, the ventral part of the dorsal raphe nucleus (DRV), and the ventrolateral part of the dorsal raphe nucleus (DRVL). Retrograde tracing methods revealed that DRDSh contains large numbers of neurons projecting to the basolateral amygdaloid nucleus, a forebrain structure important for emotional appraisal and modulation of anxiety-related physiological and behavioral responses. Together these findings support the hypothesis that there is a functional topographical organization in the DR and are consistent with the hypothesis that anxiogenic drugs have selective actions on a subpopulation of serotonergic neurons projecting to a distributed central autonomic and emotional motor control system regulating anxiety states and anxiety-related physiological and behavioral responses.