Non-hypermutator cancers access driver mutations through reversals in germline mutational bias.

Cancer is an evolutionary disease driven by mutations in asexually-reproducing somatic cells. In asexual microbes, bias reversals in the mutation spectrum can speed adaptation by increasing access to previously undersampled beneficial mutations. By analyzing tumors from 20 tissues, along with normal tissue and the germline, we demonstrate this effect in cancer. Non-hypermutated tumors reverse the germline mutation bias and have consistent spectra across tissues. These spectra changes carry the signature of hypoxia, and they facilitate positive selection in cancer genes. Hypermutated and non-hypermutated tumors thus acquire driver mutations differently: hypermutated tumors by higher mutation rates and non-hypermutated tumors by changing the mutation spectrum to reverse the germline mutation bias.

Shifts in Mutation Bias Promote Mutators by Altering the Distribution of Fitness Effects.

AbstractRecent experimental evidence demonstrates that shifts in mutational biases-for example, increases in transversion frequency-can change the distribution of fitness effects of mutations (DFE). In particular, reducing or reversing a prevailing bias can increase the probability that a de novo mutation is beneficial. It has also been shown that mutator bacteria are more likely to emerge if the beneficial mutations they generate have a larger effect size than observed in the wild type. Here, we connect these two results, demonstrating that mutator strains that reduce or reverse a prevailing bias have a positively shifted DFE, which in turn can dramatically increase their emergence probability. Since changes in mutation rate and bias are often coupled through the gain and loss of DNA repair enzymes, our results predict that the invasion of mutator strains will be facilitated by shifts in mutation bias that offer improved access to previously undersampled beneficial mutations.

Lipids regulate peripheral serotonin release via gut CD1d.

The crosstalk between the immune and neuroendocrine systems is critical for intestinal homeostasis and gut-brain communications. However, it remains unclear how immune cells participate in gut sensation of hormones and neurotransmitters release in response to environmental cues, such as self-lipids and microbial lipids. We show here that lipid-mediated engagement of invariant natural killer T (iNKT) cells with enterochromaffin (EC) cells, a subset of intestinal epithelial cells, promoted peripheral serotonin (5-HT) release via a CD1d-dependent manner, regulating gut motility and hemostasis. We also demonstrated that inhibitory sphingolipids from symbiotic microbe Bacteroides fragilis represses 5-HT release. Mechanistically, CD1d ligation on EC cells transduced a signal and restrained potassium conductance through activation of protein tyrosine kinase Pyk2, leading to calcium influx and 5-HT secretion. Together, our data reveal that by engaging with iNKT cells, gut chemosensory cells selectively perceive lipid antigens via CD1d to control 5-HT release, modulating intestinal and systemic homeostasis.

Extraordinary selection on the human X chromosome associated with archaic admixture.

The X chromosome in non-African humans shows less diversity and less Neanderthal introgression than expected from neutral evolution. Analyzing 162 human male X chromosomes worldwide, we identified fourteen chromosomal regions where nearly identical haplotypes spanning several hundred kilobases are found at high frequencies in non-Africans. Genetic drift alone cannot explain the existence of these haplotypes, which must have been associated with strong positive selection in partial selective sweeps. Moreover, the swept haplotypes are entirely devoid of archaic ancestry as opposed to the non-swept haplotypes in the same genomic regions. The ancient Ust'-Ishim male dated at 45,000 before the present (BP) also carries the swept haplotypes, implying that selection on the haplotypes must have occurred between 45,000 and 55,000 years ago. Finally, we find that the chromosomal positions of sweeps overlap previously reported hotspots of selective sweeps in great ape evolution, suggesting a mechanism of selection unique to X chromosomes.

Shisa7-Dependent Regulation of GABAA Receptor Single-Channel Gating Kinetics.

GABAA receptors (GABAARs) mediate the majority of fast inhibitory transmission throughout the brain. Although it is widely known that pore-forming subunits critically determine receptor function, it is unclear whether their single-channel properties are modulated by GABAAR-associated transmembrane proteins. We previously identified Shisa7 as a GABAAR auxiliary subunit that modulates the trafficking, pharmacology, and deactivation properties of these receptors. However, whether Shisa7 also regulates GABAAR single-channel properties has yet to be determined. Here, we performed single-channel recordings of α2ß3γ2L GABAARs cotransfected with Shisa7 in HEK293T cells and found that while Shisa7 does not change channel slope conductance, it reduced the frequency of receptor openings. Importantly, Shisa7 modulates GABAAR gating by decreasing the duration and open probability within bursts. Through kinetic analysis of individual dwell time components, activation modeling, and macroscopic simulations, we demonstrate that Shisa7 accelerates GABAAR deactivation by governing the time spent between close and open states during gating. Together, our data provide a mechanistic basis for how Shisa7 controls GABAAR gating and reveal for the first time that GABAAR single-channel properties can be modulated by an auxiliary subunit. These findings shed light on processes that shape the temporal dynamics of GABAergic transmission.SIGNIFICANCE STATEMENT Although GABAA receptor (GABAAR) single-channel properties are largely determined by pore-forming subunits, it remains unknown whether they are also controlled by GABAAR-associated transmembrane proteins. Here, we show that Shisa7, a recently identified GABAAR auxiliary subunit, modulates GABAAR activation by altering single-channel burst kinetics. These results reveal that Shisa7 primarily decreases the duration and open probability of receptor burst activity during gating, leading to accelerated GABAAR deactivation. These experiments are the first to assess the gating properties of GABAARs in the presence of an auxiliary subunit and provides a kinetic basis for how Shisa7 modifies temporal attributes of GABAergic transmission at the single-channel level.

Shisa7 phosphorylation regulates GABAergic transmission and neurodevelopmental behaviors.

GABA-A receptors (GABAARs) are crucial for development and function of the brain. Altered GABAergic transmission is hypothesized to be involved in neurodevelopmental disorders. Recently, we identified Shisa7 as a GABAAR auxiliary subunit that modulates GABAAR trafficking and GABAergic transmission. However, the underlying molecular mechanisms remain elusive. Here we generated a knock-in (KI) mouse line that is phospho-deficient at a phosphorylation site in Shisa7 (S405) and combined with electrophysiology, imaging and behavioral assays to illustrate the role of this site in GABAergic transmission and plasticity as well as behaviors. We found that expression of phospho-deficient mutants diminished α2-GABAAR trafficking in heterologous cells. Additionally, α1/α2/α5-GABAAR surface expression and GABAergic inhibition were decreased in hippocampal neurons in KI mice. Moreover, chemically induced inhibitory long-term potentiation was abolished in KI mice. Lastly, KI mice exhibited hyperactivity, increased grooming and impaired sleep homeostasis. Collectively, our study reveals a phosphorylation site critical for Shisa7-dependent GABAARs trafficking which contributes to behavioral endophenotypes displayed in neurodevelopmental disorders.

Interactions between Brainstem Neurons That Regulate the Motility to the Stomach.

Activity in the dorsal vagal complex (DVC) is essential to gastric motility regulation. We and others have previously shown that this activity is greatly influenced by local GABAergic signaling, primarily because of somatostatin (SST)-expressing GABAergic neurons. To further understand the network dynamics associated with gastric motility control in the DVC, we focused on another neuron prominently distributed in this complex, neuropeptide-Y (NPY) neurons. However, the effect of these neurons on gastric motility remains unknown. Here, we investigate the anatomic and functional characteristics of the NPY neurons in the nucleus tractus solitarius (NTS) and their interactions with SST neurons using transgenic mice of both sexes. We sought to determine whether NPY neurons influence the activity of gastric-projecting neurons, synaptically interact with SST neurons, and affect end-organ function. Our results using combined neuroanatomy and optogenetic in vitro and in vivo show that NPY neurons are part of the gastric vagal circuit as they are trans-synaptically labeled by a viral tracer from the gastric antrum, are primarily excitatory as optogenetic activation of these neurons evoke EPSCs in gastric-antrum-projecting neurons, are functionally coupled to each other and reciprocally connected to SST neurons, whose stimulation has a potent inhibitory effect on the action potential firing of the NPY neurons, and affect gastric tone and motility as reflected by their robust optogenetic response in vivo. These findings indicate that interacting NPY and SST neurons are integral to the network that controls vagal transmission to the stomach.SIGNIFICANCE STATEMENT The brainstem neurons in the dorsal nuclear complex are essential for regulating vagus nerve activity that affects the stomach via tone and motility. Two distinct nonoverlapping populations of predominantly excitatory NPY neurons and predominantly inhibitory SST neurons form reciprocal connections with each other in the NTS and with premotor neurons in the dorsal motor nucleus of the vagus to control gastric mechanics. Light activation and inhibition of NTS NPY neurons increased and decreased gastric motility, respectively, whereas both activation and inhibition of NTS SST neurons enhanced gastric motility.

Distinct regulation of tonic GABAergic inhibition by NMDA receptor subtypes.

Tonic inhibition mediated by extrasynaptic GABAARs regulates various brain functions. However, the mechanisms that regulate tonic inhibition remain largely unclear. Here, we report distinct actions of GluN2A- and GluN2B-NMDA receptors (NMDARs) on tonic inhibition in hippocampal neurons under basal and high activity conditions. Specifically, overexpression of GluN2B, but not GluN2A, reduces α5-GABAAR surface expression and tonic currents. Additionally, knockout of GluN2A and GluN2B decreases and increases tonic currents, respectively. Mechanistically, GluN2A-NMDARs inhibit and GluN2B-NMDARs promote α5-GABAAR internalization, resulting in increased and decreased surface α5-GABAAR expression, respectively. Furthermore, GluN2A-NMDARs, but not GluN2B-NMDARs, are required for homeostatic potentiation of tonic inhibition induced by prolonged increase of neuronal activity. Last, tonic inhibition decreases during acute seizures, whereas it increases 24 h later, involving GluN2-NMDAR-dependent signaling. Collectively, these data reveal an NMDAR subunit-specific regulation of tonic inhibition in physiological and pathological conditions and provide mechanistic insight into activity-dependent modulation of tonic inhibition.

Population Genomics in the Great Apes.

The great apes play an important role as model organisms. They are our closest living relatives, allowing us to identify the genetic basis of phenotypic traits that we think of as characteristically human. However, the most significant asset of great apes as model organisms is that they share with humans most of their genetic makeup. This means that we can extend our vast knowledge of the human genome, its genes, and the associated phenotypes to these species. Comparative genomic studies of humans and apes thus reveal how very similar genomes react when exposed to different population genetic regimes. In this way, each species represents a natural experiment, where a genome highly similar to the human one, is differently exposed to the evolutionary forces of demography, population structure, selection, recombination, and admixture/hybridization. The initial sequencing of reference genomes for chimpanzee, orangutan, gorilla, the bonobo, each provided new insights and a second generation of sequencing projects has provided diversity data for all the great apes. In this chapter, we will outline some of the findings that population genomic analysis of great apes has provided, and how comparative studies have helped us understand how the fundamental forces in evolution have contributed to shaping the genomes and the genetic diversity of the great apes.

Looking for Novelty in an "Old" Receptor: Recent Advances Toward Our Understanding of GABAARs and Their Implications in Receptor Pharmacology.

Diverse populations of GABAA receptors (GABAARs) throughout the brain mediate fast inhibitory transmission and are modulated by various endogenous ligands and therapeutic drugs. Deficits in GABAAR signaling underlie the pathophysiology behind neurological and neuropsychiatric disorders such as epilepsy, anxiety, and depression. Pharmacological intervention for these disorders relies on several drug classes that target GABAARs, such as benzodiazepines and more recently neurosteroids. It has been widely demonstrated that subunit composition and receptor stoichiometry impact the biophysical and pharmacological properties of GABAARs. However, current GABAAR-targeting drugs have limited subunit selectivity and produce their therapeutic effects concomitantly with undesired side effects. Therefore, there is still a need to develop more selective GABAAR pharmaceuticals, as well as evaluate the potential for developing next-generation drugs that can target accessory proteins associated with native GABAARs. In this review, we briefly discuss the effects of benzodiazepines and neurosteroids on GABAARs, their use as therapeutics, and some of the pitfalls associated with their adverse side effects. We also discuss recent advances toward understanding the structure, function, and pharmacology of GABAARs with a focus on benzodiazepines and neurosteroids, as well as newly identified transmembrane proteins that modulate GABAARs.

Impact of Mutation Rate and Selection at Linked Sites on DNA Variation across the Genomes of Humans and Other Homininae.

DNA diversity varies across the genome of many species. Variation in diversity across a genome might arise from regional variation in the mutation rate, variation in the intensity and mode of natural selection, and regional variation in the recombination rate. We show that both noncoding and nonsynonymous diversity are positively correlated to a measure of the mutation rate and the recombination rate and negatively correlated to the density of conserved sequences in 50 kb windows across the genomes of humans and nonhuman homininae. Interestingly, we find that although noncoding diversity is equally affected by these three genomic variables, nonsynonymous diversity is mostly dominated by the density of conserved sequences. The positive correlation between diversity and our measure of the mutation rate seems to be largely a direct consequence of regions with higher mutation rates having more diversity. However, the positive correlation with recombination rate and the negative correlation with the density of conserved sequences suggest that selection at linked sites also affect levels of diversity. This is supported by the observation that the ratio of the number of nonsynonymous to noncoding polymorphisms is negatively correlated to a measure of the effective population size across the genome. We show these patterns persist even when we restrict our analysis to GC-conservative mutations, demonstrating that the patterns are not driven by GC biased gene conversion. In conclusion, our comparative analyses describe how recombination rate, gene density, and mutation rate interact to produce the patterns of DNA diversity that we observe along the hominine genomes.

Shisa7 is a GABAA receptor auxiliary subunit controlling benzodiazepine actions.

The function and pharmacology of γ-aminobutyric acid type A receptors (GABAARs) are of great physiological and clinical importance and have long been thought to be determined by the channel pore-forming subunits. We discovered that Shisa7, a single-passing transmembrane protein, localizes at GABAergic inhibitory synapses and interacts with GABAARs. Shisa7 controls receptor abundance at synapses and speeds up the channel deactivation kinetics. Shisa7 also potently enhances the action of diazepam, a classic benzodiazepine, on GABAARs. Genetic deletion of Shisa7 selectively impairs GABAergic transmission and diminishes the effects of diazepam in mice. Our data indicate that Shisa7 regulates GABAAR trafficking, function, and pharmacology and reveal a previously unknown molecular interaction that modulates benzodiazepine action in the brain.

Evolutionary and functional impact of common polymorphic inversions in the human genome.

Inversions are one type of structural variants linked to phenotypic differences and adaptation in multiple organisms. However, there is still very little information about polymorphic inversions in the human genome due to the difficulty of their detection. Here, we develop a new high-throughput genotyping method based on probe hybridization and amplification, and we perform a complete study of 45 common human inversions of 0.1-415 kb. Most inversions promoted by homologous recombination occur recurrently in humans and great apes and they are not tagged by SNPs. Furthermore, there is an enrichment of inversions showing signatures of positive or balancing selection, diverse functional effects, such as gene disruption and gene-expression changes, or association with phenotypic traits. Therefore, our results indicate that the genome is more dynamic than previously thought and that human inversions have important functional and evolutionary consequences, making possible to determine for the first time their contribution to complex traits.

Comparison of the Full Distribution of Fitness Effects of New Amino Acid Mutations Across Great Apes.

The distribution of fitness effects (DFE) is central to many questions in evolutionary biology. However, little is known about the differences in DFE between closely related species. We use >9000 coding genes orthologous one-to-one across great apes, gibbons, and macaques to assess the stability of the DFE across great apes. We use the unfolded site frequency spectrum of polymorphic mutations (n = 8 haploid chromosomes per population) to estimate the DFE. We find that the shape of the deleterious DFE is strikingly similar across great apes. We confirm that effective population size (Ne ) is a strong predictor of the strength of negative selection, consistent with the nearly neutral theory. However, we also find that the strength of negative selection varies more than expected given the differences in Ne between species. Across species, mean fitness effects of new deleterious mutations covaries with Ne , consistent with positive epistasis among deleterious mutations. We find that the strength of negative selection for the smallest populations, bonobos and western chimpanzees, is higher than expected given their Ne This may result from a more efficient purging of strongly deleterious recessive variants in these populations. Forward simulations confirm that these findings are not artifacts of the way we are inferring Ne and DFE parameters. All findings are replicated using only GC-conservative mutations, thereby confirming that GC-biased gene conversion is not affecting our conclusions.

Genetic Deletion of GABAA Receptors Reveals Distinct Requirements of Neurotransmitter Receptors for GABAergic and Glutamatergic Synapse Development.

In the adult brain GABAA receptors (GABAARs) mediate the majority of synaptic inhibition that provides inhibitory balance to excitatory drive and controls neuronal output. In the immature brain GABAAR signaling is critical for neuronal development. However, the cell-autonomous role of GABAARs in synapse development remains largely unknown. We have employed the CRISPR-CAS9 technology to genetically eliminate GABAARs in individual hippocampal neurons and examined GABAergic and glutamatergic synapses. We found that development of GABAergic synapses, but not glutamatergic synapses, critically depends on GABAARs. By combining different genetic approaches, we have also removed GABAARs and two ionotropic glutamate receptors, AMPA receptors (AMPARs) and NMDA receptors (NMDARs), in single neurons and discovered a striking dichotomy. Indeed, while development of glutamatergic synapses and spines does not require signaling mediated by these receptors, inhibitory synapse formation is crucially dependent on them. Our data reveal a critical cell-autonomous role of GABAARs in inhibitory synaptogenesis and demonstrate distinct molecular mechanisms for development of inhibitory and excitatory synapses.

Adaptation and Conservation throughout the Drosophila melanogaster Life-Cycle.

Previous studies of the evolution of genes expressed at different life-cycle stages of Drosophila melanogaster have not been able to disentangle adaptive from nonadaptive substitutions when using nonsynonymous sites. Here, we overcome this limitation by combining whole-genome polymorphism data from D. melanogaster and divergence data between D. melanogaster and Drosophila yakuba. For the set of genes expressed at different life-cycle stages of D. melanogaster, as reported in modENCODE, we estimate the ratio of substitutions relative to polymorphism between nonsynonymous and synonymous sites (α) and then α is discomposed into the ratio of adaptive (ωa) and nonadaptive (ωna) substitutions to synonymous substitutions. We find that the genes expressed in mid- and late-embryonic development are the most conserved, whereas those expressed in early development and postembryonic stages are the least conserved. Importantly, we found that low conservation in early development is due to high rates of nonadaptive substitutions (high ωna), whereas in postembryonic stages it is due, instead, to high rates of adaptive substitutions (high ωa). By using estimates of different genomic features (codon bias, average intron length, exon number, recombination rate, among others), we also find that genes expressed in mid- and late-embryonic development show the most complex architecture: they are larger, have more exons, more transcripts, and longer introns. In addition, these genes are broadly expressed among all stages. We suggest that all these genomic features are related to the conservation of mid- and late-embryonic development. Globally, our study supports the hourglass pattern of conservation and adaptation over the life-cycle.

Nearly Neutral Evolution across the Drosophila melanogaster Genome.

Under the nearly neutral theory of molecular evolution, the proportion of effectively neutral mutations is expected to depend upon the effective population size (Ne). Here, we investigate whether this is the case across the genome of Drosophila melanogaster using polymorphism data from North American and African lines. We show that the ratio of the number of nonsynonymous and synonymous polymorphisms is negatively correlated to the number of synonymous polymorphisms, even when the nonindependence is accounted for. The relationship is such that the proportion of effectively neutral nonsynonymous mutations increases by ∼45% as Ne is halved. However, we also show that this relationship is steeper than expected from an independent estimate of the distribution of fitness effects from the site frequency spectrum. We investigate a number of potential explanations for this and show, using simulation, that this is consistent with a model of genetic hitchhiking: Genetic hitchhiking depresses diversity at neutral and weakly selected sites, but has little effect on the diversity of strongly selected sites.

Mapping Selection within Drosophila melanogaster Embryo's Anatomy.

We present a survey of selection across Drosophila melanogaster embryonic anatomy. Our approach integrates genomic variation, spatial gene expression patterns, and development with the aim of mapping adaptation over the entire embryo's anatomy. Our adaptation map is based on analyzing spatial gene expression information for 5,969 genes (from text-based annotations of in situ hybridization data directly from the BDGP database, Tomancak et al. 2007) and the polymorphism and divergence in these genes (from the project DGRP, Mackay et al. 2012).The proportion of nonsynonymous substitutions that are adaptive, neutral, or slightly deleterious are estimated for the set of genes expressed in each embryonic anatomical structure using the distribution of fitness effects-alpha method (Eyre-Walker and Keightley 2009). This method is a robust derivative of the McDonald and Kreitman test (McDonald and Kreitman 1991). We also explore whether different anatomical structures differ in the phylogenetic age, codon usage, or expression bias of the genes they express and whether genes expressed in many anatomical structures show more adaptive substitutions than other genes.We found that: 1) most of the digestive system and ectoderm-derived structures are under selective constraint, 2) the germ line and some specific mesoderm-derived structures show high rates of adaptive substitution, and 3) the genes that are expressed in a small number of anatomical structures show higher expression bias, lower phylogenetic ages, and less constraint.

Adaptive Evolution Is Substantially Impeded by Hill-Robertson Interference in Drosophila.

Hill-Robertson interference (HRi) is expected to reduce the efficiency of natural selection when two or more linked selected sites do not segregate freely, but no attempt has been done so far to quantify the overall impact of HRi on the rate of adaptive evolution for any given genome. In this work, we estimate how much HRi impedes the rate of adaptive evolution in the coding genome of Drosophila melanogaster. We compiled a data set of 6,141 autosomal protein-coding genes from Drosophila, from which polymorphism levels in D. melanogaster and divergence out to D. yakuba were estimated. The rate of adaptive evolution was calculated using a derivative of the McDonald-Kreitman test that controls for slightly deleterious mutations. We find that the rate of adaptive amino acid substitution at a given position of the genome is positively correlated to both the rate of recombination and the mutation rate, and negatively correlated to the gene density of the region. These correlations are robust to controlling for each other, for synonymous codon bias and for gene functions related to immune response and testes. We show that HRi diminishes the rate of adaptive evolution by approximately 27%. Interestingly, genes with low mutation rates embedded in gene poor regions lose approximately 17% of their adaptive substitutions whereas genes with high mutation rates embedded in gene rich regions lose approximately 60%. We conclude that HRi hampers the rate of adaptive evolution in Drosophila and that the variation in recombination, mutation, and gene density along the genome affects the HRi effect.

Functional Impact and Evolution of a Novel Human Polymorphic Inversion That Disrupts a Gene and Creates a Fusion Transcript.

Despite many years of study into inversions, very little is known about their functional consequences, especially in humans. A common hypothesis is that the selective value of inversions stems in part from their effects on nearby genes, although evidence of this in natural populations is almost nonexistent. Here we present a global analysis of a new 415-kb polymorphic inversion that is among the longest ones found in humans and is the first with clear position effects. This inversion is located in chromosome 19 and has been generated by non-homologous end joining between blocks of transposable elements with low identity. PCR genotyping in 541 individuals from eight different human populations allowed the detection of tag SNPs and inversion genotyping in multiple populations worldwide, showing that the inverted allele is mainly found in East Asia with an average frequency of 4.7%. Interestingly, one of the breakpoints disrupts the transcription factor gene ZNF257, causing a significant reduction in the total expression level of this gene in lymphoblastoid cell lines. RNA-Seq analysis of the effects of this expression change in standard homozygotes and inversion heterozygotes revealed distinct expression patterns that were validated by quantitative RT-PCR. Moreover, we have found a new fusion transcript that is generated exclusively from inverted chromosomes around one of the breakpoints. Finally, by the analysis of the associated nucleotide variation, we have estimated that the inversion was generated ~40,000-50,000 years ago and, while a neutral evolution cannot be ruled out, its current frequencies are more consistent with those expected for a deleterious variant, although no significant association with phenotypic traits has been found so far.