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.

Reversing anxiety by targeting a stress-responsive signaling pathway.

Current treatments of anxiety and depressive disorders are plagued by considerable side effects and limited efficacies, underscoring the need for additional molecular targets that can be leveraged to improve medications. Here, we have identified a molecular cascade triggered by chronic stress that exacerbates anxiety- and depressive-like behaviors. Specifically, chronic stress enhances Src kinase activity and tyrosine phosphorylation of calmodulin, which diminishes MyosinVa (MyoVa) interaction with Neuroligin2 (NL2), resulting in decreased inhibitory transmission and heightened anxiety-like behaviors. Importantly, pharmacological inhibition of Src reinstates inhibitory synaptic deficits and effectively reverses heightened anxiety-like behaviors in chronically stressed mice, a process requiring the MyoVa-NL2 interaction. These data demonstrate the reversibility of anxiety- and depressive-like phenotypes at both molecular and behavioral levels and uncover a therapeutic target for anxiety and depressive disorders.

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.

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.

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.

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.

The Drosophila melanogaster Genetic Reference Panel.

A major challenge of biology is understanding the relationship between molecular genetic variation and variation in quantitative traits, including fitness. This relationship determines our ability to predict phenotypes from genotypes and to understand how evolutionary forces shape variation within and between species. Previous efforts to dissect the genotype-phenotype map were based on incomplete genotypic information. Here, we describe the Drosophila melanogaster Genetic Reference Panel (DGRP), a community resource for analysis of population genomics and quantitative traits. The DGRP consists of fully sequenced inbred lines derived from a natural population. Population genomic analyses reveal reduced polymorphism in centromeric autosomal regions and the X chromosome, evidence for positive and negative selection, and rapid evolution of the X chromosome. Many variants in novel genes, most at low frequency, are associated with quantitative traits and explain a large fraction of the phenotypic variance. The DGRP facilitates genotype-phenotype mapping using the power of Drosophila genetics.

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.

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.

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.

Satellite glial contact enhances differentiation and maturation of human iPSC-derived sensory neurons.

Sensory neurons generated from induced pluripotent stem cells (iSNs) are used to model human peripheral neuropathies, however current differentiation protocols produce sensory neurons with an embryonic phenotype. Peripheral glial cells contact sensory neurons early in development and contribute to formation of the canonical pseudounipolar morphology, but these signals are not encompassed in current iSN differentiation protocols. Here, we show that terminal differentiation of iSNs in co-culture with rodent Dorsal Root Ganglion satellite glia (rSG) advances their differentiation and maturation. Co-cultured iSNs develop a pseudounipolar morphology through contact with rSGs. This transition depends on semaphorin-plexin guidance cues and on glial gap junction signaling. In addition to morphological changes, iSNs terminally differentiated in co-culture exhibit enhanced spontaneous action potential firing, more mature gene expression, and increased susceptibility to paclitaxel induced axonal degeneration. Thus, iSNs differentiated in coculture with rSGs provide a better model for investigating human peripheral neuropathies.

Scalable Hypothalamic Arcuate Neuron Differentiation from Human Pluripotent Stem Cells Suitable for Modeling Metabolic and Reproductive Disorders.

The hypothalamus, composed of several nuclei, is essential for maintaining our body's homeostasis. The arcuate nucleus (ARC), located in the mediobasal hypothalamus, contains neuronal populations with eminent roles in energy and glucose homeostasis as well as reproduction. These neuronal populations are of great interest for translational research. To fulfill this promise, we used a robotic cell culture platform to provide a scalable and chemically defined approach for differentiating human pluripotent stem cells (hPSCs) into pro-opiomelanocortin (POMC), somatostatin (SST), tyrosine hydroxylase (TH) and gonadotropin-releasing hormone (GnRH) neuronal subpopulations with an ARC-like signature. This robust approach is reproducible across several distinct hPSC lines and exhibits a stepwise induction of key ventral diencephalon and ARC markers in transcriptomic profiling experiments. This is further corroborated by direct comparison to human fetal hypothalamus, and the enriched expression of genes implicated in obesity and type 2 diabetes (T2D). Genome-wide chromatin accessibility profiling by ATAC-seq identified accessible regulatory regions that can be utilized to predict candidate enhancers related to metabolic disorders and hypothalamic development. In depth molecular, cellular, and functional experiments unveiled the responsiveness of the hPSC-derived hypothalamic neurons to hormonal stimuli, such as insulin, neuropeptides including kisspeptin, and incretin mimetic drugs such as Exendin-4, highlighting their potential utility as physiologically relevant cellular models for disease studies. In addition, differential glucose and insulin treatments uncovered adaptability within the generated ARC neurons in the dynamic regulation of POMC and insulin receptors. In summary, the establishment of this model represents a novel, chemically defined, and scalable platform for manufacturing large numbers of hypothalamic arcuate neurons and serves as a valuable resource for modeling metabolic and reproductive disorders.

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 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.

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.

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.

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.

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.