Ana B. Muñoz-Manchado

Cell types in the mouse cortex and hippocampus revealed by single-cell RNA-seq

Cellular diversity in the brain revealed The mammalian brain has an extraordinarily large number of cells. Although there are quite a few different cell types, many cells in any one category tend to look alike. Zeisel et al. analyzed the transcriptomes of mouse brain cells to reveal more than meets the eye. Interneurons of similar type were found in dissimilar regions of the brain. Oligodendrocytes that seemed to be all of one class were differentiated by their molecular signatures into a half-dozen classes. Microglia associated with blood vessels were distinguished from look-alike perivascular macrophages. Thus, the complex microanatomy of the brain can be revealed by the RNAs expressed in its cells. Science, this issue p. 1138 A close look at the genes expressed by cells in the brain reveals hidden and coordinated cellular complexity. The mammalian cerebral cortex supports cognitive functions such as sensorimotor integration, memory, and social behaviors. Normal brain function relies on a diverse set of differentiated cell types, including neurons, glia, and vasculature. Here, we have used large-scale single-cell RNA sequencing (RNA-seq) to classify cells in the mouse somatosensory cortex and hippocampal CA1 region. We found 47 molecularly distinct subclasses, comprising all known major cell types in the cortex. We identified numerous marker genes, which allowed alignment with known cell types, morphology, and location. We found a layer I interneuron expressing Pax6 and a distinct postmitotic oligodendrocyte subclass marked by Itpr2. Across the diversity of cortical cell types, transcription factors formed a complex, layered regulatory code, suggesting a mechanism for the maintenance of adult cell type identity.

Oligodendrocyte heterogeneity in the mouse juvenile and adult central nervous system

One size does not fit all Oligodendrocytes are best known for their ability to myelinate brain neurons, thus increasing the speed of signal transmission. Marques et al. surveyed oligodendrocytes of developing mice and found unexpected heterogeneity. Transcriptional analysis identified 12 populations, ranging from precursors to mature oligodendrocytes. Transcriptional profiles diverged as the oligodendrocytes matured, building distinct populations. One population was responsive to motor learning, and another, with a different transcriptome, traveled along blood vessels. Science, this issue p. 1326 Brain oligodendrocytes express transcriptional heterogeneity between brain regions and age of development. Oligodendrocytes have been considered as a functionally homogeneous population in the central nervous system (CNS). We performed single-cell RNA sequencing on 5072 cells of the oligodendrocyte lineage from 10 regions of the mouse juvenile and adult CNS. Thirteen distinct populations were identified, 12 of which represent a continuum from Pdgfra+ oligodendrocyte precursor cells (OPCs) to distinct mature oligodendrocytes. Initial stages of differentiation were similar across the juvenile CNS, whereas subsets of mature oligodendrocytes were enriched in specific regions in the adult brain. Newly formed oligodendrocytes were detected in the adult CNS and were responsive to complex motor learning. A second Pdgfra+ population, distinct from OPCs, was found along vessels. Our study reveals the dynamics of oligodendrocyte differentiation and maturation, uncoupling them at a transcriptional level and highlighting oligodendrocyte heterogeneity in the CNS.

Novel Striatal GABAergic Interneuron Populations Labeled in the 5HT3aEGFP Mouse

Histological and morphological studies indicate that approximately 5% of striatal neurons are cholinergic or γ-aminobutyric acidergic (GABAergic) interneurons (gINs). However, the number of striatal neurons expressing known interneuron markers is too small to account for the entire interneuron population. We therefore studied the serotonin (5HT) receptor 3a-enhanced green fluorescent protein (5HT3a(EGFP)) mouse, in which we found that a large number of striatal gINs are labeled. Roughly 20% of 5HT3a(EGFP)-positive cells co-express parvalbumin and exhibit fast-spiking (FS) electrophysiological properties. However, the majority of labeled neurons do not overlap with known molecular interneuron markers. Intrinsic electrical properties reveal at least 2 distinct novel subtypes: a late-spiking (LS) neuropeptide-Y (NPY)-negative neurogliaform (NGF) interneuron, and a large heterogeneous population with several features resembling low-threshold-spiking (LTS) interneurons that do not express somatostatin, NPY, or neuronal nitric oxide synthase. Although the 5HT3a(EGFP) NGF and LTS-like interneurons have electrophysiological properties similar to previously described populations, they are pharmacologically distinct. In direct contrast to previously described NPY(+) LTS and NGF cells, LTS-like 5HT3a(EGFP) cells show robust responses to nicotine administration, while the 5HT3a(EGFP) NGF cell type shows little or no response. By constructing a molecular map of the overlap between these novel populations and existing interneuron populations, we are able to reconcile the morphological and molecular estimates of striatal interneuron numbers.

Genetic identification of brain cell types underlying schizophrenia

With few exceptions, the marked advances in knowledge about the genetic basis of schizophrenia have not converged on findings that can be confidently used for precise experimental modeling. By applying knowledge of the cellular taxonomy of the brain from single-cell RNA sequencing, we evaluated whether the genomic loci implicated in schizophrenia map onto specific brain cell types. We found that the common-variant genomic results consistently mapped to pyramidal cells, medium spiny neurons (MSNs) and certain interneurons, but far less consistently to embryonic, progenitor or glial cells. These enrichments were due to sets of genes that were specifically expressed in each of these cell types. We also found that many of the diverse gene sets previously associated with schizophrenia (genes involved in synaptic function, those encoding mRNAs that interact with FMRP, antipsychotic targets, etc.) generally implicated the same brain cell types. Our results suggest a parsimonious explanation: the common-variant genetic results for schizophrenia point at a limited set of neurons, and the gene sets point to the same cells. The genetic risk associated with MSNs did not overlap with that of glutamatergic pyramidal cells and interneurons, suggesting that different cell types have biologically distinct roles in schizophrenia.Integration of single-cell RNA sequencing with genome-wide association data implicates specific brain cell types in schizophrenia. Gene sets previously associated with schizophrenia implicate the same cell types, which include pyramidal cells and medium spiny neurons.

Integrated Bayesian analysis of rare exonic variants to identify risk genes for schizophrenia and neurodevelopmental disorders

BackgroundIntegrating rare variation from trio family and case–control studies has successfully implicated specific genes contributing to risk of neurodevelopmental disorders (NDDs) including autism spectrum disorders (ASD), intellectual disability (ID), developmental disorders (DDs), and epilepsy (EPI). For schizophrenia (SCZ), however, while sets of genes have been implicated through the study of rare variation, only two risk genes have been identified.MethodsWe used hierarchical Bayesian modeling of rare-variant genetic architecture to estimate mean effect sizes and risk-gene proportions, analyzing the largest available collection of whole exome sequence data for SCZ (1,077 trios, 6,699 cases, and 13,028 controls), and data for four NDDs (ASD, ID, DD, and EPI; total 10,792 trios, and 4,058 cases and controls).ResultsFor SCZ, we estimate there are 1,551 risk genes. There are more risk genes and they have weaker effects than for NDDs. We provide power analyses to predict the number of risk-gene discoveries as more data become available. We confirm and augment prior risk gene and gene set enrichment results for SCZ and NDDs. In particular, we detected 98 new DD risk genes at FDR < 0.05. Correlations of risk-gene posterior probabilities are high across four NDDs (ρ>0.55), but low between SCZ and the NDDs (ρ<0.3). An in-depth analysis of 288 NDD genes shows there is highly significant protein–protein interaction (PPI) network connectivity, and functionally distinct PPI subnetworks based on pathway enrichment, single-cell RNA-seq cell types, and multi-region developmental brain RNA-seq.ConclusionsWe have extended a pipeline used in ASD studies and applied it to infer rare genetic parameters for SCZ and four NDDs (https://github.com/hoangtn/extTADA). We find many new DD risk genes, supported by gene set enrichment and PPI network connectivity analyses. We find greater similarity among NDDs than between NDDs and SCZ. NDD gene subnetworks are implicated in postnatally expressed presynaptic and postsynaptic genes, and for transcriptional and post-transcriptional gene regulation in prenatal neural progenitor and stem cells.

BCL11B/CTIP2 is highly expressed in GABAergic interneurons of the mouse somatosensory cortex

In the nervous system, BCL11B is crucial for the development of deep layer corticospinal projection neurons and striatal medium spiny neurons and is often used as a marker for the aforementioned cell types. However, the expression of BCL11B in subtypes of non-excitatory neurons in the primary somatosensory cortex (S1) has not been reported in the mouse. In this study we show that BCL11B is extensively expressed in S1 GABAergic interneurons, throughout the three main subgroups (somatostatin-, parvalbumin- and 5HT3a-expresssing). Almost all BCL11B positive cells in the upper S1 layers were GABAergic interneurons and surprisingly, almost 40% of the BCL11B positive neurons in layer V were GABAergic interneurons. Single cell mRNA sequencing data revealed higher Bcl11b expression in S1 interneurons compared to deep layer pyramidal neurons. The highest levels of Bcl11b expression were found within the 5HT3a population, specifically in putative neurogliaform interneuron subclasses (5HT3a-positive but not expressing vasoactive intestinal peptide). In the light of our findings we suggest caution using BCL11B as a single marker to identify neurons.

Genome-wide Analysis of Insomnia (N=1,331,010) Identifies Novel Loci and Functional Pathways

Insomnia is the second-most prevalent mental disorder, with no sufficient treatment available. Despite a substantial role of genetic factors, only a handful of genes have been implicated and insight into the associated neurobiological pathways remains limited. Here, we use an unprecedented large genetic association sample (N=1,331,010) to allow detection of a substantial number of genetic variants and gain insight into biological functions, cell types and tissues involved in insomnia. We identify 202 genome-wide significant loci implicating 956 genes through positional, eQTL and chromatin interaction mapping. We show involvement of the axonal part of neurons, of specific cortical and subcortical tissues, and of two specific cell-types in insomnia: striatal medium spiny neurons and hypothalamic neurons. These cell-types have been implicated previously in the regulation of reward processing, sleep and arousal in animal studies, but have never been genetically linked to insomnia in humans. We found weak genetic correlations with other sleep-related traits, but strong genetic correlations with psychiatric and metabolic traits. Mendelian randomization identified causal effects of insomnia on specific psychiatric and metabolic traits. Our findings reveal key brain areas and cells implicated in the neurobiology of insomnia and its related disorders, and provide novel targets for treatment.

Meta-analysis of genome-wide association studies for neuroticism in 449,484 individuals identifies novel genetic loci and pathways

Neuroticism is an important risk factor for psychiatric traits, including depression1, anxiety2,3, and schizophrenia4–6. At the time of analysis, previous genome-wide association studies7–12 (GWAS) reported 16 genomic loci associated to neuroticism10–12. Here we conducted a large GWAS meta-analysis (n = 449,484) of neuroticism and identified 136 independent genome-wide significant loci (124 new at the time of analysis), which implicate 599 genes. Functional follow-up analyses showed enrichment in several brain regions and involvement of specific cell types, including dopaminergic neuroblasts (P = 3.49 × 10−8), medium spiny neurons (P = 4.23 × 10−8), and serotonergic neurons (P = 1.37 × 10−7). Gene set analyses implicated three specific pathways: neurogenesis (P = 4.43 × 10−9), behavioral response to cocaine processes (P = 1.84 × 10−7), and axon part (P = 5.26 × 10−8). We show that neuroticism’s genetic signal partly originates in two genetically distinguishable subclusters13 (‘depressed affect’ and ‘worry’), suggesting distinct causal mechanisms for subtypes of individuals. Mendelian randomization analysis showed unidirectional and bidirectional effects between neuroticism and multiple psychiatric traits. These results enhance neurobiological understanding of neuroticism and provide specific leads for functional follow-up experiments.A meta-analysis of genome-wide association studies for neuroticism identifies novel loci, pathways and potential drug targets. Further analysis implicates specific brain regions and evaluates genetic overlap with other neuropsychiatric traits.

Biological annotation of genetic loci associated with intelligence in a meta-analysis of 87,740 individuals

Variance in IQ is associated with a wide range of health outcomes, and 1% of the population are affected by intellectual disability. Despite a century of research, the fundamental neural underpinnings of intelligence remain unclear. We integrate results from genome-wide association studies (GWAS) of intelligence with brain tissue and single cell gene expression data to identify tissues and cell types associated with intelligence. GWAS data for IQ (N = 78,308) were meta-analyzed with a study comparing 1247 individuals with mean IQ ~170 to 8185 controls. Genes associated with intelligence implicate pyramidal neurons of the somatosensory cortex and CA1 region of the hippocampus, and midbrain embryonic GABAergic neurons. Tissue-specific analyses find the most significant enrichment for frontal cortex brain expressed genes. These results suggest specific neuronal cell types and genes may be involved in intelligence and provide new hypotheses for neuroscience experiments using model systems.

Diversity of Interneurons in the Dorsal Striatum Revealed by Single-Cell RNA Sequencing and PatchSeq

Summary Striatal locally projecting neurons, or interneurons, act on nearby circuits and shape functional output to the rest of the basal ganglia. We performed single-cell RNA sequencing of striatal cells enriching for interneurons. We find seven discrete interneuron types, six of which are GABAergic. In addition to providing specific markers for the populations previously described, including those expressing Sst/Npy, Th, Npy without Sst, and Chat, we identify two small populations of cells expressing Cck with or without Vip. Surprisingly, the Pvalb-expressing cells do not constitute a discrete cluster but rather are part of a larger group of cells expressing Pthlh with a spatial gradient of Pvalb expression. Using PatchSeq, we show that Pthlh cells exhibit a continuum of electrophysiological properties correlated with expression of Pvalb. Furthermore, we find significant molecular differences that correlate with differences in electrophysiological properties between Pvalb-expressing cells of the striatum and those of the cortex.