Molecular Diversity and Specializations among the Cells of the Adult Mouse Brain.

The mammalian brain is composed of diverse, specialized cell populations. To systematically ascertain and learn from these cellular specializations, we used Drop-seq to profile RNA expression in 690,000 individual cells sampled from 9 regions of the adult mouse brain. We identified 565 transcriptionally distinct groups of cells using computational approaches developed to distinguish biological from technical signals. Cross-region analysis of these 565 cell populations revealed features of brain organization, including a gene-expression module for synthesizing axonal and presynaptic components, patterns in the co-deployment of voltage-gated ion channels, functional distinctions among the cells of the vasculature and specialization of glutamatergic neurons across cortical regions. Systematic neuronal classifications for two complex basal ganglia nuclei and the striatum revealed a rare population of spiny projection neurons. This adult mouse brain cell atlas, accessible through interactive online software (DropViz), serves as a reference for development, disease, and evolution.

Comprehensive Classification of Retinal Bipolar Neurons by Single-Cell Transcriptomics.

Patterns of gene expression can be used to characterize and classify neuronal types. It is challenging, however, to generate taxonomies that fulfill the essential criteria of being comprehensive, harmonizing with conventional classification schemes, and lacking superfluous subdivisions of genuine types. To address these challenges, we used massively parallel single-cell RNA profiling and optimized computational methods on a heterogeneous class of neurons, mouse retinal bipolar cells (BCs). From a population of ∼25,000 BCs, we derived a molecular classification that identified 15 types, including all types observed previously and two novel types, one of which has a non-canonical morphology and position. We validated the classification scheme and identified dozens of novel markers using methods that match molecular expression to cell morphology. This work provides a systematic methodology for achieving comprehensive molecular classification of neurons, identifies novel neuronal types, and uncovers transcriptional differences that distinguish types within a class.

Highly Parallel Genome-wide Expression Profiling of Individual Cells Using Nanoliter Droplets.

Cells, the basic units of biological structure and function, vary broadly in type and state. Single-cell genomics can characterize cell identity and function, but limitations of ease and scale have prevented its broad application. Here we describe Drop-seq, a strategy for quickly profiling thousands of individual cells by separating them into nanoliter-sized aqueous droplets, associating a different barcode with each cell's RNAs, and sequencing them all together. Drop-seq analyzes mRNA transcripts from thousands of individual cells simultaneously while remembering transcripts' cell of origin. We analyzed transcriptomes from 44,808 mouse retinal cells and identified 39 transcriptionally distinct cell populations, creating a molecular atlas of gene expression for known retinal cell classes and novel candidate cell subtypes. Drop-seq will accelerate biological discovery by enabling routine transcriptional profiling at single-cell resolution. VIDEO ABSTRACT.

A concerted neuron-astrocyte program declines in ageing and schizophrenia.

Human brains vary across people and over time; such variation is not yet understood in cellular terms. Here we describe a relationship between people's cortical neurons and cortical astrocytes. We used single-nucleus RNA sequencing to analyse the prefrontal cortex of 191 human donors aged 22-97 years, including healthy individuals and people with schizophrenia. Latent-factor analysis of these data revealed that, in people whose cortical neurons more strongly expressed genes encoding synaptic components, cortical astrocytes more strongly expressed distinct genes with synaptic functions and genes for synthesizing cholesterol, an astrocyte-supplied component of synaptic membranes. We call this relationship the synaptic neuron and astrocyte program (SNAP). In schizophrenia and ageing-two conditions that involve declines in cognitive flexibility and plasticity1,2-cells divested from SNAP: astrocytes, glutamatergic (excitatory) neurons and GABAergic (inhibitory) neurons all showed reduced SNAP expression to corresponding degrees. The distinct astrocytic and neuronal components of SNAP both involved genes in which genetic risk factors for schizophrenia were strongly concentrated. SNAP, which varies quantitatively even among healthy people of similar age, may underlie many aspects of normal human interindividual differences and may be an important point of convergence for multiple kinds of pathophysiology.

Comparative cellular analysis of motor cortex in human, marmoset and mouse.

The primary motor cortex (M1) is essential for voluntary fine-motor control and is functionally conserved across mammals1. Here, using high-throughput transcriptomic and epigenomic profiling of more than 450,000 single nuclei in humans, marmoset monkeys and mice, we demonstrate a broadly conserved cellular makeup of this region, with similarities that mirror evolutionary distance and are consistent between the transcriptome and epigenome. The core conserved molecular identities of neuronal and non-neuronal cell types allow us to generate a cross-species consensus classification of cell types, and to infer conserved properties of cell types across species. Despite the overall conservation, however, many species-dependent specializations are apparent, including differences in cell-type proportions, gene expression, DNA methylation and chromatin state. Few cell-type marker genes are conserved across species, revealing a short list of candidate genes and regulatory mechanisms that are responsible for conserved features of homologous cell types, such as the GABAergic chandelier cells. This consensus transcriptomic classification allows us to use patch-seq (a combination of whole-cell patch-clamp recordings, RNA sequencing and morphological characterization) to identify corticospinal Betz cells from layer 5 in non-human primates and humans, and to characterize their highly specialized physiology and anatomy. These findings highlight the robust molecular underpinnings of cell-type diversity in M1 across mammals, and point to the genes and regulatory pathways responsible for the functional identity of cell types and their species-specific adaptations.

Innovations present in the primate interneuron repertoire.

Primates and rodents, which descended from a common ancestor around 90 million years ago1, exhibit profound differences in behaviour and cognitive capacity; the cellular basis for these differences is unknown. Here we use single-nucleus RNA sequencing to profile RNA expression in 188,776 individual interneurons across homologous brain regions from three primates (human, macaque and marmoset), a rodent (mouse) and a weasel (ferret). Homologous interneuron types-which were readily identified by their RNA-expression patterns-varied in abundance and RNA expression among ferrets, mice and primates, but varied less among primates. Only a modest fraction of the genes identified as 'markers' of specific interneuron subtypes in any one species had this property in another species. In the primate neocortex, dozens of genes showed spatial expression gradients among interneurons of the same type, which suggests that regional variation in cortical contexts shapes the RNA expression patterns of adult neocortical interneurons. We found that an interneuron type that was previously associated with the mouse hippocampus-the 'ivy cell', which has neurogliaform characteristics-has become abundant across the neocortex of humans, macaques and marmosets but not mice or ferrets. We also found a notable subcortical innovation: an abundant striatal interneuron type in primates that had no molecularly homologous counterpart in mice or ferrets. These interneurons expressed a unique combination of genes that encode transcription factors, receptors and neuropeptides and constituted around 30% of striatal interneurons in marmosets and humans.

Author Correction: Innovations present in the primate interneuron repertoire.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Early role for a Na+,K+-ATPase (ATP1A3) in brain development.

Osmotic equilibrium and membrane potential in animal cells depend on concentration gradients of sodium (Na+) and potassium (K+) ions across the plasma membrane, a function catalyzed by the Na+,K+-ATPase α-subunit. Here, we describe ATP1A3 variants encoding dysfunctional α3-subunits in children affected by polymicrogyria, a developmental malformation of the cerebral cortex characterized by abnormal folding and laminar organization. To gain cell-biological insights into the spatiotemporal dynamics of prenatal ATP1A3 expression, we built an ATP1A3 transcriptional atlas of fetal cortical development using mRNA in situ hybridization and transcriptomic profiling of ∼125,000 individual cells with single-cell RNA sequencing (Drop-seq) from 11 areas of the midgestational human neocortex. We found that fetal expression of ATP1A3 is most abundant to a subset of excitatory neurons carrying transcriptional signatures of the developing subplate, yet also maintains expression in nonneuronal cell populations. Moving forward a year in human development, we profiled ∼52,000 nuclei from four areas of an infant neocortex and show that ATP1A3 expression persists throughout early postnatal development, most predominantly in inhibitory neurons, including parvalbumin interneurons in the frontal cortex. Finally, we discovered the heteromeric Na+,K+-ATPase pump complex may form nonredundant cell-type-specific α-ß isoform combinations, including α3-ß1 in excitatory neurons and α3-ß2 in inhibitory neurons. Together, the developmental malformation phenotype of affected individuals and single-cell ATP1A3 expression patterns point to a key role for α3 in human cortex development, as well as a cell-type basis for pre- and postnatal ATP1A3-associated diseases.

Cell diversity and network dynamics in photosensitive human brain organoids.

In vitro models of the developing brain such as three-dimensional brain organoids offer an unprecedented opportunity to study aspects of human brain development and disease. However, the cells generated within organoids and the extent to which they recapitulate the regional complexity, cellular diversity and circuit functionality of the brain remain undefined. Here we analyse gene expression in over 80,000 individual cells isolated from 31 human brain organoids. We find that organoids can generate a broad diversity of cells, which are related to endogenous classes, including cells from the cerebral cortex and the retina. Organoids could be developed over extended periods (more than 9 months), allowing for the establishment of relatively mature features, including the formation of dendritic spines and spontaneously active neuronal networks. Finally, neuronal activity within organoids could be controlled using light stimulation of photosensitive cells, which may offer a way to probe the functionality of human neuronal circuits using physiological sensory stimuli.

Single-nucleus sequencing reveals enriched expression of genetic risk factors in extratelencephalic neurons sensitive to degeneration in ALS.

Amyotrophic lateral sclerosis (ALS) is a neurodegenerative disorder characterized by a progressive loss of motor function linked to degenerating extratelencephalic neurons/Betz cells (ETNs). The reasons why these neurons are selectively affected remain unclear. Here, to understand the unique molecular properties that may sensitize ETNs to ALS, we performed RNA sequencing of 79,169 single nuclei from cortices of patients and controls. In both patients and unaffected individuals, we found significantly higher expression of ALS risk genes in THY1+ ETNs, regardless of diagnosis. In patients, this was accompanied by the induction of genes involved in protein homeostasis and stress responses that were significantly induced in a wide collection of ETNs. Examination of oligodendroglial and microglial nuclei revealed patient-specific downregulation of myelinating genes in oligodendrocytes and upregulation of an endolysosomal reactive state in microglia. Our findings suggest that selective vulnerability of extratelencephalic neurons is partly connected to their intrinsic molecular properties sensitizing them to genetics and mechanisms of degeneration.

Concerted neuron-astrocyte gene expression declines in aging and schizophrenia.

Human brains vary across people and over time; such variation is not yet understood in cellular terms. Here we describe a striking relationship between people's cortical neurons and cortical astrocytes. We used single-nucleus RNA-seq to analyze the prefrontal cortex of 191 human donors ages 22-97 years, including healthy individuals and persons with schizophrenia. Latent-factor analysis of these data revealed that in persons whose cortical neurons more strongly expressed genes for synaptic components, cortical astrocytes more strongly expressed distinct genes with synaptic functions and genes for synthesizing cholesterol, an astrocyte-supplied component of synaptic membranes. We call this relationship the Synaptic Neuron-and-Astrocyte Program (SNAP). In schizophrenia and aging - two conditions that involve declines in cognitive flexibility and plasticity 1,2 - cells had divested from SNAP: astrocytes, glutamatergic (excitatory) neurons, and GABAergic (inhibitory) neurons all reduced SNAP expression to corresponding degrees. The distinct astrocytic and neuronal components of SNAP both involved genes in which genetic risk factors for schizophrenia were strongly concentrated. SNAP, which varies quantitatively even among healthy persons of similar age, may underlie many aspects of normal human interindividual differences and be an important point of convergence for multiple kinds of pathophysiology.

Key Roles of CACNA1C/Cav1.2 and CALB1/Calbindin in Prefrontal Neurons Altered in Cognitive Disorders.

Importance: The risk of mental disorders is consistently associated with variants in CACNA1C (L-type calcium channel Cav1.2) but it is not known why these channels are critical to cognition, and whether they affect the layer III pyramidal cells in the dorsolateral prefrontal cortex that are especially vulnerable in cognitive disorders. Objective: To examine the molecular mechanisms expressed in layer III pyramidal cells in primate dorsolateral prefrontal cortices. Design, Setting, and Participants: The design included transcriptomic analyses from human and macaque dorsolateral prefrontal cortex, and connectivity, protein expression, physiology, and cognitive behavior in macaques. The research was performed in academic laboratories at Yale, Harvard, Princeton, and the University of Pittsburgh. As dorsolateral prefrontal cortex only exists in primates, the work evaluated humans and macaques. Main Outcomes and Measures: Outcome measures included transcriptomic signatures of human and macaque pyramidal cells, protein expression and interactions in layer III macaque pyramidal cells using light and electron microscopy, changes in neuronal firing during spatial working memory, and working memory performance following pharmacological treatments. Results: Layer III pyramidal cells in dorsolateral prefrontal cortex coexpress a constellation of calcium-related proteins, delineated by CALB1 (calbindin), and high levels of CACNA1C (Cav1.2), GRIN2B (NMDA receptor GluN2B), and KCNN3 (SK3 potassium channel), concentrated in dendritic spines near the calcium-storing smooth endoplasmic reticulum. L-type calcium channels influenced neuronal firing needed for working memory, where either blockade or increased drive by ß1-adrenoceptors, reduced neuronal firing by a mean (SD) 37.3% (5.5%) or 40% (6.3%), respectively, the latter via SK potassium channel opening. An L-type calcium channel blocker or ß1-adrenoceptor antagonist protected working memory from stress. Conclusions and Relevance: The layer III pyramidal cells in the dorsolateral prefrontal cortex especially vulnerable in cognitive disorders differentially express calbindin and a constellation of calcium-related proteins including L-type calcium channels Cav1.2 (CACNA1C), GluN2B-NMDA receptors (GRIN2B), and SK3 potassium channels (KCNN3), which influence memory-related neuronal firing. The finding that either inadequate or excessive L-type calcium channel activation reduced neuronal firing explains why either loss- or gain-of-function variants in CACNA1C were associated with increased risk of cognitive disorders. The selective expression of calbindin in these pyramidal cells highlights the importance of regulatory mechanisms in neurons with high calcium signaling, consistent with Alzheimer tau pathology emerging when calbindin is lost with age and/or inflammation.

Sibling chimerism among microglia in marmosets.

Chimerism happens rarely among most mammals but is common in marmosets and tamarins, a result of fraternal twin or triplet birth patterns in which in utero connected circulatory systems (through which stem cells transit) lead to persistent blood chimerism (12-80%) throughout life. The presence of Y-chromosome DNA sequences in other organs of female marmosets has long suggested that chimerism might also affect these organs. However, a longstanding question is whether this chimerism is driven by blood-derived cells or involves contributions from other cell types. To address this question, we analyzed single-cell RNA-seq data from blood, liver, kidney and multiple brain regions across a number of marmosets, using transcribed single nucleotide polymorphisms (SNPs) to identify cells with the sibling's genome in various cell types within these tissues. Sibling-derived chimerism in all tissues arose entirely from cells of hematopoietic origin (i.e., myeloid and lymphoid lineages). In brain tissue this was reflected as sibling-derived chimerism among microglia (20-52%) and macrophages (18-64%) but not among other resident cell types (i.e., neurons, glia or ependymal cells). The percentage of microglia that were sibling-derived showed significant variation across brain regions, even within individual animals, likely reflecting distinct responses by siblings' microglia to local recruitment or proliferation cues or, potentially, distinct clonal expansion histories in different brain areas. In the animals and tissues we analyzed, microglial gene expression profiles bore a much stronger relationship to local/host context than to sibling genetic differences. Naturally occurring marmoset chimerism will provide new ways to understand the effects of genes, mutations and brain contexts on microglial biology and to distinguish between effects of microglia and other cell types on brain phenotypes.

A marmoset brain cell census reveals regional specialization of cellular identities.

The mammalian brain is composed of many brain structures, each with its own ontogenetic and developmental history. We used single-nucleus RNA sequencing to sample over 2.4 million brain cells across 18 locations in the common marmoset, a New World monkey primed for genetic engineering, and examined gene expression patterns of cell types within and across brain structures. The adult transcriptomic identity of most neuronal types is shaped more by developmental origin than by neurotransmitter signaling repertoire. Quantitative mapping of GABAergic types with single-molecule FISH (smFISH) reveals that interneurons in the striatum and neocortex follow distinct spatial principles, and that lateral prefrontal and other higher-order cortical association areas are distinguished by high proportions of VIP+ neurons. We use cell type-specific enhancers to drive AAV-GFP and reconstruct the morphologies of molecularly resolved interneuron types in neocortex and striatum. Our analyses highlight how lineage, local context, and functional class contribute to the transcriptional identity and biodistribution of primate brain cell types.

Brain single cell transcriptomic profiles in episodic memory phenotypes associated with temporal lobe epilepsy.

Memory dysfunction is prevalent in temporal lobe epilepsy (TLE), but little is known about the underlying molecular etiologies. Single-nucleus RNA sequencing technology was used to examine differences in cellular heterogeneity among left (language-dominant) temporal neocortical tissues from patients with TLE with (n = 4) or without (n = 2) impairment in verbal episodic memory. We observed marked cell heterogeneity between memory phenotypes and identified numerous differentially expressed genes across all brain cell types. The most notable differences were observed in glutamatergic (excitatory) and GABAergic (inhibitory) neurons with an overrepresentation of genes associated with long-term potentiation, long-term depression, and MAPK signaling, processes known to be essential for episodic memory formation.

Single-cell RNA sequencing reveals compromised immune microenvironment in precursor stages of multiple myeloma.

Precursor states of Multiple Myeloma (MM) and its native tumor microenvironment need in-depth molecular characterization to better stratify and treat patients at risk. Using single-cell RNA sequencing of bone marrow cells from precursor stages, MGUS and smoldering myeloma (SMM), to full-blown MM alongside healthy donors, we demonstrate early immune changes during patient progression. We find NK cell abundance is frequently increased in early stages, and associated with altered chemokine receptor expression. As early as SMM, we show loss of GrK+ memory cytotoxic T-cells, and show their critical role in MM immunosurveillance in mouse models. Finally, we report MHC class II dysregulation in CD14+ monocytes, which results in T cell suppression in vitro. These results provide a comprehensive map of immune changes at play over the evolution of pre-malignant MM, which will help develop strategies for immune-based patient stratification.

Genetically Distinct Parallel Pathways in the Entopeduncular Nucleus for Limbic and Sensorimotor Output of the Basal Ganglia.

The basal ganglia (BG) integrate inputs from diverse sensorimotor, limbic, and associative regions to guide action-selection and goal-directed behaviors. The entopeduncular nucleus (EP) is a major BG output nucleus and has been suggested to channel signals from distinct BG nuclei to target regions involved in diverse functions. Here we use single-cell transcriptional and molecular analyses to demonstrate that the EP contains at least three classes of projection neurons-glutamate/GABA co-releasing somatostatin neurons, glutamatergic parvalbumin neurons, and GABAergic parvalbumin neurons. These classes comprise functionally and anatomically distinct output pathways that differentially affect EP target regions, such as the lateral habenula (LHb) and thalamus. Furthermore, LHb- and thalamic-projecting EP neurons are differentially innervated by subclasses of striatal and pallidal neurons. Therefore, we identify previously unknown subdivisions within the EP and reveal the existence of cascading, molecularly distinct projections through striatum and globus pallidus to EP targets within epithalamus and thalamus.

A molecular census of arcuate hypothalamus and median eminence cell types.

The hypothalamic arcuate-median eminence complex (Arc-ME) controls energy balance, fertility and growth through molecularly distinct cell types, many of which remain unknown. To catalog cell types in an unbiased way, we profiled gene expression in 20,921 individual cells in and around the adult mouse Arc-ME using Drop-seq. We identify 50 transcriptionally distinct Arc-ME cell populations, including a rare tanycyte population at the Arc-ME diffusion barrier, a new leptin-sensing neuron population, multiple agouti-related peptide (AgRP) and pro-opiomelanocortin (POMC) subtypes, and an orexigenic somatostatin neuron population. We extended Drop-seq to detect dynamic expression changes across relevant physiological perturbations, revealing cell type-specific responses to energy status, including distinct responses in AgRP and POMC neuron subtypes. Finally, integrating our data with human genome-wide association study data implicates two previously unknown neuron populations in the genetic control of obesity. This resource will accelerate biological discovery by providing insights into molecular and cell type diversity from which function can be inferred.

Reduction of 9-nitrocamptothecin-triggered apoptosis in DU-145 human prostate cancer cells by ectopic expression of 14-3-3zeta.

A promising family of anticancer agents, the camptothecins, is noted for their ability to induce apoptosis specifically in malignant cells. However, a major obstacle for successful cancer treatment by these and other chemotherapeutic agents is the intrinsic or acquired resistance to drug treatment. Resistance to 9NC6, a camptothecin derivative, has been modeled in vitro using a human prostate cancer cell line. To elucidate the mechanism for acquired 9NC resistance, we have used a subtractive cloning approach to identify genes whose altered expression level is reflective of 9NC resistance or susceptibility. Differential gene expression was compared between wild-type human prostate cancer cell line, DU-145, and a 9NC-resistant subline, RC1. Results were confirmed by Northern and Western blot analyses. In this report, we focus on one gene, 14-3-3zeta. An expression vector of a full-length myc-epitope-tagged 14-3-3zeta cDNA was constructed and used for transfection into DU-145 cells. We consistently observed that 14-3-3zeta message and protein levels were dramatically increased in 9NC resistant cells. The expression levels of other 14-3-3 family members were unaffected. Strikingly, ectopic overexpression of 14-3-3zeta in wild-type 9NC-susceptible prostate cancer cells decreased 9NC-induced apoptosis. Our results suggest a novel direct or indirect role of 14-3-3zeta in mediating resistance of DU-145 cells to the topoisomerase I inhibitor, 9NC. We are currently exploring whether this represents a more general pathway for drug resistance as well.