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.

Dynamics and turnover of memory CD8 T cell responses following yellow fever vaccination.

Understanding how immunological memory lasts a lifetime requires quantifying changes in the number of memory cells as well as how their division and death rates change over time. We address these questions by using a statistically powerful mixed-effects differential equations framework to analyze data from two human studies that follow CD8 T cell responses to the yellow fever vaccine (YFV-17D). Models were first fit to the frequency of YFV-specific memory CD8 T cells and deuterium enrichment in those cells 42 days to 1 year post-vaccination. A different dataset, on the loss of YFV-specific CD8 T cells over three decades, was used to assess out of sample predictions of our models. The commonly used exponential and bi-exponential decline models performed relatively poorly. Models with the cell loss following a power law (exactly or approximately) were most predictive. Notably, using only the first year of data, these models accurately predicted T cell frequencies up to 30 years post-vaccination. Our analyses suggest that division rates of these cells drop and plateau at a low level (0.1% per day, ∼ double the estimated values for naive T cells) within one year following vaccination, whereas death rates continue to decline for much longer. Our results show that power laws can be predictive for T cell memory, a finding that may be useful for vaccine evaluation and epidemiological modeling. Moreover, since power laws asymptotically decline more slowly than any exponential decline, our results help explain the longevity of immune memory phenomenologically.

Population structure, intergroup interaction, and human contact govern infectious disease impacts in mountain gorilla populations.

Infectious zoonotic diseases are a threat to wildlife conservation and global health. They are especially a concern for wild apes, which are vulnerable to many human infectious diseases. As ecotourism, deforestation, and great ape field research increase, the threat of human-sourced infections to wild populations becomes more substantial and could result in devastating population declines. The endangered mountain gorillas (Gorilla beringei beringei) of the Virunga Massif in east-central Africa suffer periodic disease outbreaks and are exposed to infections from human-sourced pathogens. It is important to understand the possible risks of disease introduction and spread in this population and how human contact may facilitate disease transmission. Here we present and evaluate an individual-based, stochastic, discrete-time disease transmission model to predict epidemic outcomes and better understand health risks to the Virunga mountain gorilla population. To model disease transmission we have derived estimates for gorilla contact, interaction, and migration rates. The model shows that the social structure of gorilla populations plays a profound role in governing disease impacts with subdivided populations experiencing less than 25% of the outbreak levels of a single homogeneous population. It predicts that gorilla group dispersal and limited group interactions are strong factors in preventing widespread population-level outbreaks of infectious disease after such diseases have been introduced into the population. However, even a moderate amount of human contact increases disease spread and can lead to population-level outbreaks.

Linking high GC content to the repair of double strand breaks in prokaryotic genomes.

Genomic GC content varies widely among microbes for reasons unknown. While mutation bias partially explains this variation, prokaryotes near-universally have a higher GC content than predicted solely by this bias. Debate surrounds the relative importance of the remaining explanations of selection versus biased gene conversion favoring GC alleles. Some environments (e.g. soils) are associated with a high genomic GC content of their inhabitants, which implies that either high GC content is a selective adaptation to particular habitats, or that certain habitats favor increased rates of gene conversion. Here, we report a novel association between the presence of the non-homologous end joining DNA double-strand break repair pathway and GC content; this observation suggests that DNA damage may be a fundamental driver of GC content, leading in part to the many environmental patterns observed to-date. We discuss potential mechanisms accounting for the observed association, and provide preliminary evidence that sites experiencing higher rates of double-strand breaks are under selection for increased GC content relative to the genomic background.

Exploring the functional composition of the human microbiome using a hand-curated microbial trait database.

BACKGROUND: Even when microbial communities vary wildly in their taxonomic composition, their functional composition is often surprisingly stable. This suggests that a functional perspective could provide much deeper insight into the principles governing microbiome assembly. Much work to date analyzing the functional composition of microbial communities, however, relies heavily on inference from genomic features. Unfortunately, output from these methods can be hard to interpret and often suffers from relatively high error rates. RESULTS: We built and analyzed a domain-specific microbial trait database from known microbe-trait pairs recorded in the literature to better understand the functional composition of the human microbiome. Using a combination of phylogentically conscious machine learning tools and a network science approach, we were able to link particular traits to areas of the human body, discover traits that determine the range of body areas a microbe can inhabit, and uncover drivers of metabolic breadth. CONCLUSIONS: Domain-specific trait databases are an effective compromise between noisy methods to infer complex traits from genomic data and exhaustive, expensive attempts at database curation from the literature that do not focus on any one subset of taxa. They provide an accurate account of microbial traits and, by limiting the number of taxa considered, are feasible to build within a reasonable time-frame. We present a database specific for the human microbiome, in the hopes that this will prove useful for research into the functional composition of human-associated microbial communities.

Selection influences naive CD8+ TCR-ß repertoire sharing.

Within each individual, the adaptive immune system generates a repertoire of cells expressing receptors capable of recognizing diverse potential pathogens. The theoretical diversity of the T-cell receptor (TCR) repertoire exceeds the actual size of the T-cell population in an individual by several orders of magnitude - making the observation of identical TCRs in different individuals extremely improbable if all receptors were equally likely. Despite this disparity between the theoretical and the realized diversity of the repertoire, these 'public' receptor sequences have been identified in autoimmune, cancer and pathogen interaction contexts. Biased generation processes explain the presence of public TCRs in the naive repertoire, but do not adequately explain the different abundances of these public TCRs. We investigate and characterize the distribution of genomic TCR-ß sequences of naive CD8+ T cells from three genetically identical mice, comparing non-productive (non-functional sequences) and productive sequences. We find public TCR-ß sequences at higher abundances compared with unshared sequences in the productive, but not in the non-productive, repertoire. We show that neutral processes such as recombination biases, codon degeneracy and generation probability do not fully account for these differences, and conclude that thymic or peripheral selection plays an important role in increasing the abundances of public TCR-ß sequences.

The Rewarding and Anxiolytic Properties of Ethanol within the Central Nucleus of the Amygdala: Mediated by Genetic Background and Nociceptin.

In humans, alcohol is consumed for its rewarding and anxiolytic effects. The central nucleus of the amygdala (CeA) is considered a neuronal nexus that regulates fear, anxiety, and drug self-administration. Manipulations of the CeA alter ethanol (EtOH) consumption under numerous EtOH self-administration models. The experiments determined whether EtOH is reinforcing/anxiolytic within the CeA, whether selective breeding for high alcohol consumption alters the rewarding properties of EtOH in the CeA, and whether the reinforcing/anxiolytic effects of EtOH in the CeA are mediated by the neuropeptides corticotropin-releasing factor (CRF) and nociceptin. The reinforcing properties of EtOH were determined by having male Wistar and Taconic alcohol-preferring (tP) rats self-administer EtOH directly into the CeA. The expression of anxiety-like behaviors was assessed through multiple behavioral models (social interaction, acoustic startle, and open field). Coadministration of EtOH and a CRF1 antagonist (NBI35965) or nociceptin on self-administration into the CeA and anxiety-like behaviors was determined. EtOH was self-administered directly into the lateral CeA, and tP rats self-administered a lower concentration of EtOH than Wistar rats. EtOH microinjected into the lateral CeA reduced the expression of anxiety-like behaviors, indicating an anxiolytic effect. Coadministration of NBI35965 failed to alter the rewarding/anxiolytic properties of EtOH in the CeA. In contrast, coadministration of the nociceptin enhanced both EtOH reward and anxiolysis in the CeA. Overall, the data indicate that the lateral CeA is a key anatomic location that mediates the rewarding and anxiolytic effects of EtOH, and local nociceptin receptors, but not local CRF1 receptors, are involved in these behaviors. SIGNIFICANCE STATEMENT: Alcohol is consumed for the stimulatory, rewarding, and anxiolytic properties of the drug of abuse. The current data are the first to establish that alcohol is reinforcing and anxiolytic within the lateral central nucleus of the amygdala (CeA) and that the nociceptin system regulates these effects of alcohol within the CeA.

Genome sequence of a 45,000-year-old modern human from western Siberia.

We present the high-quality genome sequence of a ∼45,000-year-old modern human male from Siberia. This individual derives from a population that lived before-or simultaneously with-the separation of the populations in western and eastern Eurasia and carries a similar amount of Neanderthal ancestry as present-day Eurasians. However, the genomic segments of Neanderthal ancestry are substantially longer than those observed in present-day individuals, indicating that Neanderthal gene flow into the ancestors of this individual occurred 7,000-13,000 years before he lived. We estimate an autosomal mutation rate of 0.4 × 10(-9) to 0.6 × 10(-9) per site per year, a Y chromosomal mutation rate of 0.7 × 10(-9) to 0.9 × 10(-9) per site per year based on the additional substitutions that have occurred in present-day non-Africans compared to this genome, and a mitochondrial mutation rate of 1.8 × 10(-8) to 3.2 × 10(-8) per site per year based on the age of the bone.

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.

Select panicogenic drugs and stimuli induce consistent increases in tail skin flushes and decreases in core body temperature.

Panic attacks (PAs) are episodes of intense fear or discomfort that are accompanied by a variety of both psychological and somatic symptoms. Panic induction in preclinical models (e.g. rats) has largely been assayed through flight and avoidance behavioral tests and cardiorespiratory activity. Yet, the literature pertaining to PAs shows that thermal sensations (hot flushes/heat sensations and chills) are also a common symptom during PAs in humans. Considering that temperature alterations are objectively measurable in rodents, we hypothesized that select panicogenic drugs and stimuli induce consistent changes in thermoregulation related to hot flushes and chills. Specifically, we challenged male rats with intraperitoneal injections of the GABAergic inverse agonist FG-7142; the α2 adrenoceptor antagonist yohimbine; the serotonin agonist D-fenfluramine, and 20% CO2 (an interoceptive homeostatic challenge). We assayed core body temperature and tail skin temperature using implanted radiotelemetry probes and tail thermistors/thermal imaging camera, respectively, and found that all challenges elicited rapid, high-amplitude (~7-9°C) increase in tail skin temperature and delayed decreases (~1-3°C) in core body temperature. We propose that thermal sensations such as these may be an additional indicator of a panic response in rodents and humans, as these panicogenic compounds or stimuli are known to precipitate PAs in persons with panic disorder.

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.

Initial viral load determines the magnitude of the human CD8 T cell response to yellow fever vaccination.

CD8 T cells are a potent tool for eliminating intracellular pathogens and tumor cells. Thus, eliciting robust CD8 T-cell immunity is the basis for many vaccines under development. However, the relationship between antigen load and the magnitude of the CD8 T-cell response is not well-described in a human immune response. Here we address this issue by quantifying viral load and the CD8 T-cell response in a cohort of 80 individuals immunized with the live attenuated yellow fever vaccine (YFV-17D) by sampling peripheral blood at days 0, 1, 2, 3, 5, 7, 9, 11, 14, 30, and 90. When the virus load was below a threshold (peak virus load < 225 genomes per mL, or integrated virus load < 400 genome days per mL), the magnitude of the CD8 T-cell response correlated strongly with the virus load (R(2) ∼ 0.63). As the virus load increased above this threshold, the magnitude of the CD8 T-cell responses saturated. Recent advances in CD8 T-cell-based vaccines have focused on replication-incompetent or single-cycle vectors. However, these approaches deliver relatively limited amounts of antigen after immunization. Our results highlight the requirement that T-cell-based vaccines should deliver sufficient antigen during the initial period of the immune response to elicit a large number of CD8 T cells that may be needed for protection.

Genome-wide ancestry of 17th-century enslaved Africans from the Caribbean.

Between 1500 and 1850, more than 12 million enslaved Africans were transported to the New World. The vast majority were shipped from West and West-Central Africa, but their precise origins are largely unknown. We used genome-wide ancient DNA analyses to investigate the genetic origins of three enslaved Africans whose remains were recovered on the Caribbean island of Saint Martin. We trace their origins to distinct subcontinental source populations within Africa, including Bantu-speaking groups from northern Cameroon and non-Bantu speakers living in present-day Nigeria and Ghana. To our knowledge, these findings provide the first direct evidence for the ethnic origins of enslaved Africans, at a time for which historical records are scarce, and demonstrate that genomic data provide another type of record that can shed new light on long-standing historical questions.

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.

A selective orexin-1 receptor antagonist attenuates stress-induced hyperarousal without hypnotic effects.

Orexins (OXs) are peptides produced by perifornical (PeF) and lateral hypothalamic neurons that exert a prominent role in arousal-related processes, including stress. A critical role for the orexin-1 receptor (OX1R) in complex emotional behavior is emerging, such as overactivation of the OX1R pathway being associated with panic or anxiety states. Here we characterize a brain-penetrant, selective, and high-affinity OX1R antagonist, compound 56 [N-({3-[(3-ethoxy-6-methylpyridin-2-yl)carbonyl]-3-azabicyclo[4.1.0]hept-4-yl}methyl)-5-(trifluoromethyl)pyrimidin-2-amine]. Ex vivo receptor binding studies demonstrated that, after subcutaneous administration, compound 56 crossed the blood-brain barrier and occupied OX1Rs in the rat brain at lower doses than standard OX1R antagonists GSK-1059865 [5-bromo-N-({1-[(3-fluoro-2-methoxyphenyl)carbonyl]-5-methylpiperidin-2-yl}methyl)pyridin-2-amine], SB-334867 [1-(2-methyl-1,3-benzoxazol-6-yl)-3-(1,5-naphthyridin-4-yl)urea], and SB-408124 [1-(6,8-difluoro-2-methylquinolin-4-yl)-3-[4-(dimethylamino)phenyl]urea]. Although compound 56 did not alter spontaneous sleep in rats and in wild-type mice, its administration in orexin-2 receptor knockout mice selectively promoted rapid eye movement sleep, demonstrating target engagement and specific OX1R blockade. In a rat model of psychological stress induced by cage exchange, the OX1R antagonist prevented the prolongation of sleep onset without affecting sleep duration. In a rat model of panic vulnerability (involving disinhibition of the PeF OX region) to threatening internal state changes (i.e., intravenous sodium lactate infusion), compound 56 attenuated sodium lactate-induced panic-like behaviors and cardiovascular responses without altering baseline locomotor or autonomic activity. In conclusion, OX1R antagonism represents a novel therapeutic strategy for the treatment of various psychiatric disorders associated with stress or hyperarousal states.

OREXIN 1 AND 2 RECEPTOR INVOLVEMENT IN CO2 -INDUCED PANIC-ASSOCIATED BEHAVIOR AND AUTONOMIC RESPONSES.

BACKGROUND: The neuropeptides orexin A and B play a role in reward and feeding and are critical for arousal. However, it was not initially appreciated that most prepro-orexin synthesizing neurons are almost exclusively concentrated in the perifornical hypothalamus, which when stimulated elicits panic-associated behavior and cardiovascular responses in rodents and self-reported "panic attacks" and "fear of dying" in humans. More recent studies support a role for the orexin system in coordinating an integrative stress response. For instance, orexin neurons are highly reactive to anxiogenic stimuli, are hyperactive in anxiety pathology, and have strong projections to anxiety and panic-associated circuitry. Although the two cognate orexin receptors are colocalized in many brain regions, the orexin 2 receptor (OX2R) most robustly maps to the histaminergic wake-promoting region, while the orexin 1 receptor (OX1R) distribution is more exclusive and dense in anxiety and panic circuitry regions, such as the locus ceruleus. Overall, this suggests that OX1Rs play a critical role in mobilizing anxiety and panic responses. METHODS: Here, we used a CO2 -panic provocation model to screen a dual OX1/2R antagonist (DORA-12) to globally inhibit orexin activity, then a highly selective OX1R antagonist (SORA1, Compound 56) or OX2R antagonist (SORA2, JnJ10397049) to assess OX1R and OX2R involvement. RESULTS: All compounds except the SORA2 attenuated CO2 -induced anxiety-like behaviors, and all but the SORA2 and DORA attenuated CO2 -induced cardiovascular responses. CONCLUSIONS: SORA1s may represent a novel method of treating anxiety disorders, with no apparent sedative effects that were present with a benzodiazepine.

Peripheral selection rather than thymic involution explains sudden contraction in naive CD4 T-cell diversity with age.

A diverse array of T cells is required for defense against pathogens. The naive CD4 T-cell repertoire reaches its peak diversity by early human adulthood and is maintained until older age. Surprisingly, around age 70, this diversity appears to plummet abruptly. A similar qualitative pattern holds for the CD4 T memory-cell population. We used mathematical models to explore different hypotheses for how such a loss of diversity might occur. The prevailing hypotheses suggest that the loss of diversity is due to a decline in emigration of cells from the thymus or a contraction in total number of cells. Our models reject these mechanisms because they yield only a gradual and minimal decline in the repertoire instead of the observed sudden and profound decrease later in life. We propose that an abrupt decline in the repertoire could be caused by the accumulation of mutations (defined here as any cell-intrinsic heritable event) that provide a short-term fitness advantage to a small number of T-cell clones (e.g., by an increased division rate or decreased death rate), with the person as a whole incurring the long-term cost of a decreased ability to fight infections.

A population biological approach to understanding the maintenance and loss of the T-cell repertoire during aging.

The adaptive immune system requires a diverse T-cell repertoire to be able to respond to a wide variety of pathogens. Worryingly, the repertoire diversity declines dramatically in old age. As thymic output generates novel T cells, the conventional view holds that a decrease in this output with age is responsible for the loss in the repertoire. However, many additional factors affect the repertoire such as homeostatic turnover and antigen-dependent expansion in response to infection. Mathematical models taking a population biology perspective are important tools for understanding how the interplay between these factors affects the immune repertoire. These models suggest that thymic decline is not a major factor but rather that some combination of virus-induced proliferation and T-cell-intrinsic genetic or epigenetic changes gives rise to the oligoclonal expansions that cause the decline in T-cell diversity. We also discuss consequences for strategies to rejuvenate the immune repertoire in old age.

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/). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.