Molecular Architecture of the Mouse Nervous System.

The mammalian nervous system executes complex behaviors controlled by specialized, precisely positioned, and interacting cell types. Here, we used RNA sequencing of half a million single cells to create a detailed census of cell types in the mouse nervous system. We mapped cell types spatially and derived a hierarchical, data-driven taxonomy. Neurons were the most diverse and were grouped by developmental anatomical units and by the expression of neurotransmitters and neuropeptides. Neuronal diversity was driven by genes encoding cell identity, synaptic connectivity, neurotransmission, and membrane conductance. We discovered seven distinct, regionally restricted astrocyte types that obeyed developmental boundaries and correlated with the spatial distribution of key glutamate and glycine neurotransmitters. In contrast, oligodendrocytes showed a loss of regional identity followed by a secondary diversification. The resource presented here lays a solid foundation for understanding the molecular architecture of the mammalian nervous system and enables genetic manipulation of specific cell types.

Schizophrenia genomics: genetic complexity and functional insights.

Determining the causes of schizophrenia has been a notoriously intractable problem, resistant to a multitude of investigative approaches over centuries. In recent decades, genomic studies have delivered hundreds of robust findings that implicate nearly 300 common genetic variants (via genome-wide association studies) and more than 20 rare variants (via whole-exome sequencing and copy number variant studies) as risk factors for schizophrenia. In parallel, functional genomic and neurobiological studies have provided exceptionally detailed information about the cellular composition of the brain and its interconnections in neurotypical individuals and, increasingly, in those with schizophrenia. Taken together, these results suggest unexpected complexity in the mechanisms that drive schizophrenia, pointing to the involvement of ensembles of genes (polygenicity) rather than single-gene causation. In this Review, we describe what we now know about the genetics of schizophrenia and consider the neurobiological implications of this information.

Genome-wide association meta-analysis identifies 48 risk variants and highlights the role of the stria vascularis in hearing loss.

Hearing loss is one of the top contributors to years lived with disability and is a risk factor for dementia. Molecular evidence on the cellular origins of hearing loss in humans is growing. Here, we performed a genome-wide association meta-analysis of clinically diagnosed and self-reported hearing impairment on 723,266 individuals and identified 48 significant loci, 10 of which are novel. A large proportion of associations comprised missense variants, half of which lie within known familial hearing loss loci. We used single-cell RNA-sequencing data from mouse cochlea and brain and mapped common-variant genomic results to spindle, root, and basal cells from the stria vascularis, a structure in the cochlea necessary for normal hearing. Our findings indicate the importance of the stria vascularis in the mechanism of hearing impairment, providing future paths for developing targets for therapeutic intervention in hearing loss.

Probabilistic cell typing enables fine mapping of closely related cell types in situ.

Understanding the function of a tissue requires knowing the spatial organization of its constituent cell types. In the cerebral cortex, single-cell RNA sequencing (scRNA-seq) has revealed the genome-wide expression patterns that define its many, closely related neuronal types, but cannot reveal their spatial arrangement. Here we introduce probabilistic cell typing by in situ sequencing (pciSeq), an approach that leverages previous scRNA-seq classification to identify cell types using multiplexed in situ RNA detection. We applied this method by mapping the inhibitory neurons of mouse hippocampal area CA1, for which ground truth is available from extensive previous work identifying their laminar organization. Our method identified these neuronal classes in a spatial arrangement matching ground truth, and further identified multiple classes of isocortical pyramidal cell in a pattern matching their known organization. This method will allow identifying the spatial organization of closely related cell types across the brain and other tissues.

Postnatal Sox6 Regulates Synaptic Function of Cortical Parvalbumin-Expressing Neurons.

Cortical parvalbumin-expressing (Pvalb+) neurons provide robust inhibition to neighboring pyramidal neurons, crucial for the proper functioning of cortical networks. This class of inhibitory neurons undergoes extensive synaptic formation and maturation during the first weeks after birth and continue to dynamically maintain their synaptic output throughout adulthood. While several transcription factors, such as Nkx2-1, Lhx6, and Sox6, are known to be necessary for the differentiation of progenitors into Pvalb+ neurons, which transcriptional programs underlie the postnatal maturation and maintenance of Pvalb+ neurons' innervation and synaptic function remains largely unknown. Because Sox6 is continuously expressed in Pvalb+ neurons until adulthood, we used conditional knock-out strategies to investigate its putative role in the postnatal maturation and synaptic function of cortical Pvalb+ neurons in mice of both sexes. We found that early postnatal loss of Sox6 in Pvalb+ neurons leads to failure of synaptic bouton growth, whereas later removal in mature Pvalb+ neurons in the adult causes shrinkage of already established synaptic boutons. Paired recordings between Pvalb+ neurons and pyramidal neurons revealed reduced release probability and increased failure rate of Pvalb+ neurons' synaptic output. Furthermore, Pvalb+ neurons lacking Sox6 display reduced expression of full-length tropomyosin-receptor kinase B (TrkB), a key modulator of GABAergic transmission. Once re-expressed in neurons lacking Sox6, TrkB was sufficient to rescue the morphologic synaptic phenotype. Finally, we showed that Sox6 mRNA levels were increased by motor training. Our data thus suggest a constitutive role for Sox6 in the maintenance of synaptic output from Pvalb+ neurons into adulthood.SIGNIFICANCE STATEMENT Cortical parvalbumin-expressing (Pvalb+) inhibitory neurons provide robust inhibition to neighboring pyramidal neurons, crucial for the proper functioning of cortical networks. These inhibitory neurons undergo extensive synaptic formation and maturation during the first weeks after birth and continue to dynamically maintain their synaptic output throughout adulthood. However, it remains largely unknown which transcriptional programs underlie the postnatal maturation and maintenance of Pvalb+ neurons. Here, we show that the transcription factor Sox6 cell-autonomously regulates the synaptic maintenance and output of Pvalb+ neurons until adulthood, leaving unaffected other maturational features of this neuronal population.

Conditional GWAS analysis to identify disorder-specific SNPs for psychiatric disorders.

Substantial genetic liability is shared across psychiatric disorders but less is known about risk variants that are specific to a given disorder. We used multi-trait conditional and joint analysis (mtCOJO) to adjust GWAS summary statistics of one disorder for the effects of genetically correlated traits to identify putative disorder-specific SNP associations. We applied mtCOJO to summary statistics for five psychiatric disorders from the Psychiatric Genomics Consortium-schizophrenia (SCZ), bipolar disorder (BIP), major depression (MD), attention-deficit hyperactivity disorder (ADHD) and autism (AUT). Most genome-wide significant variants for these disorders had evidence of pleiotropy (i.e., impact on multiple psychiatric disorders) and hence have reduced mtCOJO conditional effect sizes. However, subsets of genome-wide significant variants had larger conditional effect sizes consistent with disorder-specific effects: 15 of 130 genome-wide significant variants for schizophrenia, 5 of 40 for major depression, 3 of 11 for ADHD and 1 of 2 for autism. We show that decreased expression of VPS29 in the brain may increase risk to SCZ only and increased expression of CSE1L is associated with SCZ and MD, but not with BIP. Likewise, decreased expression of PCDHA7 in the brain is linked to increased risk of MD but decreased risk of SCZ and BIP.

Classes and continua of hippocampal CA1 inhibitory neurons revealed by single-cell transcriptomics.

Understanding any brain circuit will require a categorization of its constituent neurons. In hippocampal area CA1, at least 23 classes of GABAergic neuron have been proposed to date. However, this list may be incomplete; additionally, it is unclear whether discrete classes are sufficient to describe the diversity of cortical inhibitory neurons or whether continuous modes of variability are also required. We studied the transcriptomes of 3,663 CA1 inhibitory cells, revealing 10 major GABAergic groups that divided into 49 fine-scale clusters. All previously described and several novel cell classes were identified, with three previously described classes unexpectedly found to be identical. A division into discrete classes, however, was not sufficient to describe the diversity of these cells, as continuous variation also occurred between and within classes. Latent factor analysis revealed that a single continuous variable could predict the expression levels of several genes, which correlated similarly with it across multiple cell types. Analysis of the genes correlating with this variable suggested it reflects a range from metabolically highly active faster-spiking cells that proximally target pyramidal cells to slower-spiking cells targeting distal dendrites or interneurons. These results elucidate the complexity of inhibitory neurons in one of the simplest cortical structures and show that characterizing these cells requires continuous modes of variation as well as discrete cell classes.

Spatiotemporal mapping of RNA editing in the developing mouse brain using in situ sequencing reveals regional and cell-type-specific regulation.

BACKGROUND: Adenosine-to-inosine (A-to-I) RNA editing is a process that contributes to the diversification of proteins that has been shown to be essential for neurotransmission and other neuronal functions. However, the spatiotemporal and diversification properties of RNA editing in the brain are largely unknown. Here, we applied in situ sequencing to distinguish between edited and unedited transcripts in distinct regions of the mouse brain at four developmental stages, and investigate the diversity of the RNA landscape. RESULTS: We analyzed RNA editing at codon-altering sites using in situ sequencing at single-cell resolution, in combination with the detection of individual ADAR enzymes and specific cell type marker transcripts. This approach revealed cell-type-specific regulation of RNA editing of a set of transcripts, and developmental and regional variation in editing levels for many of the targeted sites. We found increasing editing diversity throughout development, which arises through regional- and cell type-specific regulation of ADAR enzymes and target transcripts. CONCLUSIONS: Our single-cell in situ sequencing method has proved useful to study the complex landscape of RNA editing and our results indicate that this complexity arises due to distinct mechanisms of regulating individual RNA editing sites, acting both regionally and in specific cell types.

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.

Transcriptomic correlates of electrophysiological and morphological diversity within and across excitatory and inhibitory neuron classes.

In order to further our understanding of how gene expression contributes to key functional properties of neurons, we combined publicly accessible gene expression, electrophysiology, and morphology measurements to identify cross-cell type correlations between these data modalities. Building on our previous work using a similar approach, we distinguished between correlations which were "class-driven," meaning those that could be explained by differences between excitatory and inhibitory cell classes, and those that reflected graded phenotypic differences within classes. Taking cell class identity into account increased the degree to which our results replicated in an independent dataset as well as their correspondence with known modes of ion channel function based on the literature. We also found a smaller set of genes whose relationships to electrophysiological or morphological properties appear to be specific to either excitatory or inhibitory cell types. Next, using data from PatchSeq experiments, allowing simultaneous single-cell characterization of gene expression and electrophysiology, we found that some of the gene-property correlations observed across cell types were further predictive of within-cell type heterogeneity. In summary, we have identified a number of relationships between gene expression, electrophysiology, and morphology that provide testable hypotheses for future studies.

Transcription and Signaling Regulators in Developing Neuronal Subtypes of Mouse and Human Enteric Nervous System.

BACKGROUND & AIMS: The enteric nervous system (ENS) regulates gastrointestinal function via different subtypes of neurons, organized into fine-tuned neural circuits. It is not clear how cell diversity is created within the embryonic ENS; information required for development of cell-based therapies and models of enteric neuropathies. We aimed to identify proteins that regulate ENS differentiation and network formation. METHODS: We generated and compared RNA expression profiles of the entire ENS, ENS progenitor cells, and non-ENS gut cells of mice, collected at embryonic days 11.5 and 15.5, when different subtypes of neurons are formed. Gastrointestinal tissues from R26ReYFP reporter mice crossed to Sox10-CreERT2 or Wnt1-Cre mice were dissected and the 6 populations of cells were isolated by flow cytometry. We used histochemistry to map differentially expressed proteins in mouse and human gut tissues at different stages of development, in different regions. We examined enteric neuronal diversity and gastric function in Wnt1-Cre x Sox6fl/fl mice, which do not express the Sox6 gene in the ENS. RESULTS: We identified 147 transcription and signaling factors that varied in spatial and temporal expression during development of the mouse ENS. Of the factors also analyzed in human ENS, most were conserved. We uncovered 16 signaling pathways (such as fibroblast growth factor and Eph/ephrin pathways). Transcription factors were grouped according to their specific expression in enteric progenitor cells (such as MEF2C), enteric neurons (such as SOX4), or neuron subpopulations (such as SATB1 and SOX6). Lack of SOX6 in the ENS reduced the numbers of gastric dopamine neurons and delayed gastric emptying. CONCLUSIONS: Using transcriptome and histochemical analyses of the developing mouse and human ENS, we mapped expression patterns of transcription and signaling factors. Further studies of these candidate determinants might elucidate the mechanisms by which enteric stem cells differentiate into neuronal subtypes and form distinct connectivity patterns during ENS development. We found expression of SOX6 to be required for development of gastric dopamine neurons.

Histone H2AX-dependent GABA(A) receptor regulation of stem cell proliferation.

Stem cell self-renewal implies proliferation under continued maintenance of multipotency. Small changes in numbers of stem cells may lead to large differences in differentiated cell numbers, resulting in significant physiological consequences. Proliferation is typically regulated in the G1 phase, which is associated with differentiation and cell cycle arrest. However, embryonic stem (ES) cells may lack a G1 checkpoint. Regulation of proliferation in the 'DNA damage' S/G2 cell cycle checkpoint pathway is known for its role in the maintenance of chromatin structural integrity. Here we show that autocrine/paracrine gamma-aminobutyric acid (GABA) signalling by means of GABA(A) receptors negatively controls ES cell and peripheral neural crest stem (NCS) cell proliferation, preimplantation embryonic growth and proliferation in the boundary-cap stem cell niche, resulting in an attenuation of neuronal progenies from this stem cell niche. Activation of GABA(A) receptors leads to hyperpolarization, increased cell volume and accumulation of stem cells in S phase, thereby causing a rapid decrease in cell proliferation. GABA(A) receptors signal through S-phase checkpoint kinases of the phosphatidylinositol-3-OH kinase-related kinase family and the histone variant H2AX. This signalling pathway critically regulates proliferation independently of differentiation, apoptosis and overt damage to DNA. These results indicate the presence of a fundamentally different mechanism of proliferation control in these stem cells, in comparison with most somatic cells, involving proteins in the DNA damage checkpoint pathway.

Distinct genomic signatures and modifiable risk factors underly the comorbidity between major depressive disorder and cardiovascular disease.

Major depressive disorder (MDD) and cardiovascular disease (CVD) are often comorbid, resulting in excess morbidity and mortality. Using genomic data, this study elucidates biological mechanisms, key risk factors, and causal pathways underlying their comorbidity. We show that CVDs share a large proportion of their genetic risk factors with MDD. Multivariate genome-wide association analysis of the shared genetic liability between MDD and atherosclerotic CVD (ASCVD) revealed seven novel loci and distinct patterns of tissue and brain cell-type enrichments, suggesting a role for the thalamus. Part of the genetic overlap was explained by shared inflammatory, metabolic, and psychosocial/lifestyle risk factors. Finally, we found support for causal effects of genetic liability to MDD on CVD risk, but not from most CVDs to MDD, and demonstrated that the causal effects were partly explained by metabolic and psychosocial/lifestyle factors. The distinct signature of MDD-ASCVD comorbidity aligns with the idea of an immunometabolic sub-type of MDD more strongly associated with CVD than overall MDD. In summary, we identify plausible biological mechanisms underlying MDD-CVD comorbidity, as well as key modifiable risk factors for prevention of CVD in individuals with MDD.

Distinct biological signature and modifiable risk factors underlie the comorbidity between major depressive disorder and cardiovascular disease.

Major depressive disorder (MDD) and cardiovascular disease (CVD) are often comorbid, resulting in excess morbidity and mortality. Here we show that CVDs share most of their genetic risk factors with MDD. Multivariate genome-wide association analysis of shared genetic liability between MDD and atherosclerotic CVD revealed seven loci and distinct patterns of tissue and brain cell-type enrichments, suggesting the involvement of the thalamus. Part of the genetic overlap was explained by shared inflammatory, metabolic and psychosocial or lifestyle risk factors. Our data indicated causal effects of genetic liability to MDD on CVD risk, but not from most CVDs to MDD, and showed that the causal effects were partly explained by metabolic and psychosocial or lifestyle factors. The distinct signature of MDD-atherosclerotic CVD comorbidity suggests an immunometabolic subtype of MDD that is more strongly associated with CVD than overall MDD. In summary, we identified biological mechanisms underlying MDD-CVD comorbidity and modifiable risk factors for prevention of CVD in individuals with MDD.

Efficient expansion and dopaminergic differentiation of human fetal ventral midbrain neural stem cells by midbrain morphogens.

Human fetal midbrain tissue grafting has provided proof-of-concept for dopamine cell replacement therapy (CRT) in Parkinson's disease (PD). However, limited tissue availability has hindered the development and widespread use of this experimental therapy. Here we present a method for generating large numbers of midbrain dopaminergic (DA) neurons based on expanding and differentiating neural stem/progenitor cells present in the human ventral midbrain (hVM) tissue. Our results show that hVM neurospheres (hVMN) with low cell numbers, unlike their rodent counterparts, expand the total number of cells 3-fold, whilst retaining their capacity to differentiate into midbrain DA neurons. Moreover, Wnt5a promoted DA differentiation of expanded cells resulting in improved morphological maturation, midbrain DA marker expression, DA release and electrophysiological properties. This method results in cell preparations that, after expansion and differentiation, can contain 6-fold more midbrain DA neurons than the starting VM preparation. Thus, our results provide evidence that by improving expansion and differentiation of progenitors present in the hVM it is possible to greatly enrich cell preparations for DA neurons. This method could substantially reduce the amount of human fetal midbrain tissue necessary for CRT in patients with PD, which could have major implications for the widespread adoption of this approach.

En masse in vitro functional profiling of the axonal mechanosensitivity of sensory neurons.

Perception of the environment relies on somatosensory neurons. Mechanosensory, proprioceptor and many nociceptor subtypes of these neurons have specific mechanosensitivity profiles to adequately differentiate stimulus patterns. Nevertheless, the cellular basis of differential mechanosensation remains largely elusive. Successful transduction of sensory information relies on the recruitment of sensory neurons and mechanosensation occurring at their peripheral axonal endings in vivo. Conspicuously, existing in vitro models aimed to decipher molecular mechanisms of mechanosensation test single sensory neuron somata at any one time. Here, we introduce a compartmental in vitro chamber design to deliver precisely controlled mechanical stimulation of sensory axons with synchronous real-time imaging of Ca(2+) transients in neuronal somata that reliably reflect action potential firing patterns. We report of three previously not characterized types of mechanosensitive neuron subpopulations with distinct intrinsic axonal properties tuned specifically to static indentation or vibration stimuli, showing that different classes of sensory neurons are tuned to specific types of mechanical stimuli. Primary receptor currents of vibration neurons display rapidly adapting conductance reliably detected for every single stimulus during vibration and are consistently converted into action potentials. This result allows for the characterization of two critical steps of mechanosensation in vivo: primary signal detection and signal conversion into specific action potential firing patterns in axons.

Transcriptional diversity in specific synaptic gene sets discriminates cortical neuronal identity.

Synapse diversity has been described from different perspectives, ranging from the specific neurotransmitters released, to their diverse biophysical properties and proteome profiles. However, synapse diversity at the transcriptional level has not been systematically identified across all synapse populations in the brain. To quantify and identify specific synaptic features of neuronal cell types we combined the SynGO (Synaptic Gene Ontology) database with single-cell RNA sequencing data of the mouse neocortex. We show that cell types can be discriminated by synaptic genes alone with the same power as all genes. The cell type discriminatory power is not equally distributed across synaptic genes as we could identify functional categories and synaptic compartments with greater cell type specific expression. Synaptic genes, and specific SynGO categories, belonged to three different types of gene modules: gradient expression over all cell types, gradient expression in selected cell types and cell class- or type-specific profiles. This data provides a deeper understanding of synapse diversity in the neocortex and identifies potential markers to selectively identify synapses from specific neuronal populations.

Transcriptional maintenance of cortical somatostatin interneuron subtype identity during migration.

Although cardinal cortical interneuron identity is established upon cell-cycle exit, it remains unclear whether specific interneuron subtypes are pre-established, and if so, how their identity is maintained prior to circuit integration. We conditionally removed Sox6 (Sox6-cKO) in migrating somatostatin (Sst+) interneurons and assessed the effects on their mature identity. In adolescent mice, five of eight molecular Sst+ subtypes were nearly absent in the Sox6-cKO cortex without a reduction in cell number. Sox6-cKO cells displayed electrophysiological maturity and expressed genes enriched within the broad class of Sst+ interneurons. Furthermore, we could infer subtype identity prior to cortical integration (embryonic day 18.5), suggesting that the loss in subtype was due to disrupted subtype maintenance. Conversely, Sox6 removal at postnatal day 7 did not disrupt marker expression in the mature cortex. Therefore, Sox6 is necessary during migration for maintenance of Sst+ subtype identity, indicating that subtype maintenance requires active transcriptional programs.

Single-cell transcriptomics reveals correct developmental dynamics and high-quality midbrain cell types by improved hESC differentiation.

Stem cell technologies provide new opportunities for modeling cells in health and disease and for regenerative medicine. In both cases, developmental knowledge and defining the molecular properties and quality of the cell types is essential. In this study, we identify developmental factors important for the differentiation of human embryonic stem cells (hESCs) into functional midbrain dopaminergic (mDA) neurons. We found that laminin-511, and dual canonical and non-canonical WNT activation followed by GSK3ß inhibition plus FGF8b, improved midbrain patterning. In addition, neurogenesis and differentiation were enhanced by activation of liver X receptors and inhibition of fibroblast growth factor signaling. Moreover, single-cell RNA-sequencing analysis revealed a developmental dynamics similar to that of the endogenous human ventral midbrain and the emergence of high-quality molecularly defined midbrain cell types, including mDA neurons. Our study identifies novel factors important for human midbrain development and opens the door for a future application of molecularly defined hESC-derived cell types in Parkinson disease.

Sticking out of the crowd: the molecular identity and development of cholecystokinin-containing basket cells.

Certain essential cognitive processes require the precise temporal interplay between glutamatergic (excitatory) pyramidal cells and γ-aminobutyric acid (GABA)-releasing inhibitory interneurons in the hippocampus. Basket cells, the main class of interneurons, target pyramidal cell somata and proximal dendrites and thus are poised to modify network oscillations. Though only present in limited numbers, the impaired development of basket cells can result in changes in the hippocampal circuitry leading to neurological disorders, such as schizophrenia. The diversity of the spatial origins, neurochemical make-up, cytoarchitecture and network contributions amongst basket cells is a provocative example of interneuron heterogeneity in the hippocampus. This review discusses recent data concerned with the developmental trajectories of one subclass, the cholecystokinin-containing basket cell, and emphasizes the significance of the short-range intercellular guidance cues that have recently emerged to impact the formation and function of their inhibitory synapses.