Search

The neocortex is disproportionately expanded in human compared with mouse1,2, both in its total volume relative to subcortical structures and in the proportion occupied by supragranular layers composed of neurons that selectively make connections within the neocortex and with other telencephalic structures. Single-cell transcriptomic analyses of human and mouse neocortex show an increased diversity of glutamatergic neuron types in supragranular layers in human neocortex and pronounced gradients as a function of cortical depth3. Here, to probe the functional and anatomical correlates of this transcriptomic diversity, we developed a robust platform combining patch clamp recording, biocytin staining and single-cell RNA-sequencing (Patch-seq) to examine neurosurgically resected human tissues. We demonstrate a strong correspondence between morphological, physiological and transcriptomic phenotypes of five human glutamatergic supragranular neuron types. These were enriched in but not restricted to layers, with one type varying continuously in all phenotypes across layers 2 and 3. The deep portion of layer 3 contained highly distinctive cell types, two of which express a neurofilament protein that labels long-range projection neurons in primates that are selectively depleted in Alzheimer’s disease4,5. Together, these results demonstrate the explanatory power of transcriptomic cell-type classification, provide a structural underpinning for increased complexity of cortical function in humans, and implicate discrete transcriptomic neuron types as selectively vulnerable in disease., Combined patch clamp recording, biocytin staining and single-cell RNA-sequencing of human neurocortical neurons shows an expansion of glutamatergic neuron types relative to mouse that characterizes the greater complexity of the human neocortex.

The ability to spatially resolve the cellular architecture of human cortical cell types over informative areas is essential to understanding brain function. We combined in situ sequencing gene expression data and single-nucleus RNA-sequencing cell type definitions to spatially map cells in sections of the human cortex via probabilistic cell typing. We mapped and classified a total of 59,816 cells into all 75 previously defined subtypes to create a first spatial atlas of human cortical cells in their native position, their abundances and genetic signatures. We also examined the precise within- and across-layer distributions of all the cell types and provide a resource for the cell atlas community. The abundances and locations presented here could serve as a reference for further studies, that include human brain tissues and disease applications at the cell type level., Langseth et al. present cell type maps of human cortical tissue sections and show an efficient and robust workflow to accurately resolve anatomical organization of human brain tissue.

The left temporal lobe is an integral part of the language system and its cortical structure and function associate with general intelligence. However, whether cortical laminar architecture and cellular properties of this brain area relate to verbal intelligence is unknown. Here, we addressed this using histological analysis and cellular recordings of neurosurgically resected temporal cortex in combination with presurgical IQ scores. We find that subjects with higher general and verbal IQ scores have thicker left (but not right) temporal cortex (Brodmann area 21, BA21). The increased thickness is due to the selective increase in layers 2 and 3 thickness, accompanied by lower neuron densities, and larger dendrites and cell body size of pyramidal neurons in these layers. Furthermore, these neurons sustain faster action potential kinetics, which improves information processing. Our results indicate that verbal mental ability associates with selective adaptations of supragranular layers and their cellular micro-architecture and function in left, but not right temporal cortex.

Single-cell genomics is rapidly advancing our knowledge of the diversity of cell phenotypes, including both cell types and cell states. Driven by single-cell/-nucleus RNA sequencing (scRNA-seq), comprehensive cell atlas projects characterizing a wide range of organisms and tissues are currently underway. As a result, it is critical that the transcriptional phenotypes discovered are defined and disseminated in a consistent and concise manner. Molecular biomarkers have historically played an important role in biological research, from defining immune cell types by surface protein expression to defining diseases by their molecular drivers. Here, we describe a machine learning-based marker gene selection algorithm, NS-Forest version 2.0, which leverages the nonlinear attributes of random forest feature selection and a binary expression scoring approach to discover the minimal marker gene expression combinations that optimally capture the cell type identity represented in complete scRNA-seq transcriptional profiles. The marker genes selected provide an expression barcode that serves as both a useful tool for downstream biological investigation and the necessary and sufficient characteristics for semantic cell type definition. The use of NS-Forest to identify marker genes for human brain middle temporal gyrus cell types reveals the importance of cell signaling and noncoding RNAs in neuronal cell type identity.

Increasing interest in studies of prenatal human brain development, particularly using new single-cell genomics and anatomical technologies to create cell atlases, creates a strong need for accurate and detailed anatomical reference atlases. In this study, we present two cellular-resolution digital anatomical atlases for prenatal human brain at post-conceptional weeks (PCW) 15 and 21. Both atlases were annotated on sequential Nissl-stained sections covering brain-wide structures on the basis of combined analysis of cytoarchitecture, acetylcholinesterase staining and an extensive marker gene expression dataset. This high information content dataset allowed reliable and accurate demarcation of developing cortical and subcortical structures and their subdivisions. Furthermore, using the anatomical atlases as a guide, spatial expression of 37 and 5 genes from the brains respectively at PCW 15 and 21 was annotated, illustrating reliable marker genes for many developing brain structures. Finally, the present study uncovered several novel developmental features, such as the lack of an outer subventricular zone in the hippocampal formation and entorhinal cortex, and the apparent extension of both cortical (excitatory) and subcortical (inhibitory) progenitors into the prenatal olfactory bulb. These comprehensive atlases provide useful tools for visualization, segmentation, targeting, imaging and interpretation of brain structures of prenatal human brain, and for guiding and interpreting the next generation of cell census and connectome studies. This article is protected by copyright. All rights reserved.

Single cell/nucleus RNA sequencing (scRNAseq) is emerging as an essential tool to unravel the phenotypic heterogeneity of cells in complex biological systems. While computational methods for scRNAseq cell type clustering have advanced, the ability to integrate datasets to identify common and novel cell types across experiments remains a challenge. Here, we introduce a cluster-to-cluster cell type matching method—FR-Match—that utilizes supervised feature selection for dimensionality reduction and incorporates shared information among cells to determine whether two cell type clusters share the same underlying multivariate gene expression distribution. FR-Match is benchmarked with existing cell-to-cell and cell-to-cluster cell type matching methods using both simulated and real scRNAseq data. FR-Match proved to be a stringent method that produced fewer erroneous matches of distinct cell subtypes and had the unique ability to identify novel cell phenotypes in new datasets. In silico validation demonstrated that the proposed workflow is the only self-contained algorithm that was robust to increasing numbers of true negatives (i.e. non-represented cell types). FR-Match was applied to two human brain scRNAseq datasets sampled from cortical layer 1 and full thickness middle temporal gyrus. When mapping cell types identified in specimens isolated from these overlapping human brain regions, FR-Match precisely recapitulated the laminar characteristics of matched cell type clusters, reflecting their distinct neuroanatomical distributions. An R package and Shiny application are provided at https://github.com/JCVenterInstitute/FRmatch for users to interactively explore and match scRNAseq cell type clusters with complementary visualization tools.

Neurons in the developing brain undergo extensive structural refinement as nascent circuits adopt their mature form1. This transformation is facilitated by the engulfment and degradation of excess axonal branches and inappropriate synapses by surrounding glial cells, including microglia and astrocytes2,3. However, the small size of phagocytic organelles and the complex, highly ramified morphology of glia has made it difficult to determine the contribution of these and other glial cell types to this process. Here, we used large scale, serial electron microscopy (ssEM) with computational volume segmentation to reconstruct the complete 3D morphologies of distinct glial types in the mouse visual cortex. Unexpectedly, we discovered that the fine processes of oligodendrocyte precursor cells (OPCs), a population of abundant, highly dynamic glial progenitors4, frequently surrounded terminal axon branches and included numerous phagolysosomes (PLs) containing fragments of axons and presynaptic terminals. Single- nucleus RNA sequencing indicated that cortical OPCs express key phagocytic genes, as well as neuronal transcripts, consistent with active axonal engulfment. PLs were ten times more abundant in OPCs than in microglia in P36 mice, and declined with age and lineage progression, suggesting that OPCs contribute very substantially to refinement of neuronal circuits during later phases of cortical development.

To elucidate cortical microcircuit structure and synaptic properties we present a unique, extensive, and public synaptic physiology dataset and analysis platform. Through its application, we reveal principles that relate cell type to synapse properties and intralaminar circuit organization in the mouse and human cortex. The dynamics of excitatory synapses align with the postsynaptic cell subclass, whereas inhibitory synapse dynamics partly align with presynaptic cell subclass but with considerable overlap. Despite these associations, synaptic properties are heterogeneous in most subclass to subclass connections. The two main axes of heterogeneity are strength and variability. Cell subclasses divide along the variability axis, while the strength axis accounts for significant heterogeneity within the subclass. In human cortex, excitatory to excitatory synapse dynamics are distinct from those in mouse and short-term plasticity varies with depth across layers 2 and 3. With a novel connectivity analysis that enables fair comparisons between circuit elements, we find that intralaminar connection probability among cell subclasses exhibits a strong layer dependence.These and other findings combined with the analysis platform create new opportunities for the neuroscience community to advance our understanding of cortical microcircuits.

SUMMARY Viral genetic tools that target specific brain cell types could transform basic neuroscience and targeted gene therapy. Here, we use comparative open chromatin analysis to identify thousands of human-neocortical-sub-class-specific putative enhancers from across the genome to control gene expression in adeno-associated virus (AAV) vectors. The cellular specificity of reporter expression from enhancer-AAVs is established by molecular profiling after systemic AAV delivery in mouse. Over 30% of enhancer-AAVs produce specific expression in the targeted subclass, including both excitatory and inhibitory subclasses. We present a collection of Parvalbumin (PVALB) enhancer-AAVs that show highly enriched expression not only in cortical PVALB cells but also in some subcortical PVALB populations. Five vectors maintain PVALB-enriched expression in primate neocortex. These results demonstrate how genome-wide open chromatin data mining and cross-species AAV validation can be used to create the next generation of non-species-restricted viral genetic tools., Graphical Abstract, In brief Viral genetic tools that target specific brain cell types could transform basic neuroscience and targeted gene therapy. Mich et al. use comparative open chromatin analysis to identify human neocortical enhancers that can drive gene expression from AAV vectors in a subclass-specific fashion in multiple species.

The expansion of supragranular cortical layers is thought to have enabled evolutionary development of human cognition and language. However, whether increased volume of supragranular cortical layers can actually support greater cognitive and language abilities in humans has not been demonstrated. Here, we find that subjects with higher general and verbal intelligence test (VIQ) scores have selectively expanded layers 2 and 3 only in the left temporal cortex, an area associated with language and IQ-test performance. This expansion is accompanied by lower neuron densities and larger cell-body size. Furthermore, individuals with higher VIQ scores had neurons with larger dendritic trees in left temporal cortex, potentially impacting their function. Indeed, neurons of subjects with higher VIQ scores had faster action potential upstroke kinetics, which improves information processing. These data show that expansion of supragranular layer volume, cortical and cellular micro-architecture and function are associated with improved verbal mental ability in human subjects.

Abundant evidence supports the presence of at least three distinct types of thalamocortical (TC) neurons in the primate dorsal lateral geniculate nucleus (dLGN) of the thalamus, the brain region that conveys visual information from the retina to the primary visual cortex (V1). Different types of TC neurons in mice, humans, and macaques have distinct morphologies, distinct connectivity patterns, and convey different aspects of visual information to the cortex. To investigate the molecular underpinnings of these cell types, and how these relate to differences in dLGN between human, macaque, and mice, we profiled gene expression in single nuclei and cells using RNA-sequencing. These efforts identified four distinct types of TC neurons in the primate dLGN: magnocellular (M) neurons, parvocellular (P) neurons, and two types of koniocellular (K) neurons. Despite extensively documented morphological and physiological differences between M and P neurons, we identified few genes with significant differential expression between transcriptomic cell types corresponding to these two neuronal populations. Likewise, the dominant feature of TC neurons of the adult mouse dLGN is high transcriptomic similarity, with an axis of heterogeneity that aligns with core vs. shell portions of mouse dLGN. Together, these data show that transcriptomic differences between principal cell types in the mature mammalian dLGN are subtle relative to the observed differences in morphology and cortical projection targets. Finally, alignment of transcriptome profiles across species highlights expanded diversity of GABAergic neurons in primate versus mouse dLGN and homologous types of TC neurons in primates that are distinct from TC neurons in mouse.

The Patch-seq approach is a powerful variation of the standard patch clamp technique that allows for the combined electrophysiological, morphological, and transcriptomic characterization of individual neurons. To generate Patch-seq datasets at a scale and quality that can be integrated with high-throughput dissociated cell transcriptomic data, we have optimized the technique by identifying and refining key factors that contribute to the efficient collection of high-quality data. To rapidly generate high-quality electrophysiology data, we developed patch clamp electrophysiology software with analysis functions specifically designed to automate acquisition with online quality control. We recognized a substantial improvement in transcriptomic data quality when the nucleus was extracted following the recording. For morphology success, the importance of maximizing the neuron’s membrane integrity during the extraction of the nucleus was much more critical to success than varying the duration of the electrophysiology recording. We compiled the lab protocol with the analysis and acquisition software athttps://github.com/AllenInstitute/patchseqtools. This resource can be used by individual labs to generate Patch-seq data across diverse mammalian species and that is compatible with recent large-scale publicly available Allen Institute Patch-seq datasets.

In the neocortex, subcerebral axonal projections originate largely from layer 5 (L5) extratelencephalic-projecting (ET) neurons. The highly distinctive morpho-electric properties of these neurons have mainly been described in rodents, where ET neurons can be labeled by retrograde tracers or transgenic lines. Similar labeling strategies are not possible in the human neocortex, rendering the translational relevance of findings in rodents unclear. We leveraged the recent discovery of a transcriptomically-defined L5 ET neuron type to study the properties of human L5 ET neurons in neocortical brain slices derived from neurosurgeries. Patch-seq recordings, where transcriptome, physiology and morphology are assayed from the same cell, revealed many conserved morpho-electric properties of human and rodent L5 ET neurons. Divergent properties were also apparent but were often smaller than differences between cell types within these two species. These data suggest a conserved function of L5 ET neurons in the neocortical hierarchy, but also highlight marked phenotypic divergence possibly related to functional specialization of human neocortex.

Single cell genomics is rapidly advancing our knowledge of cell phenotypic types and states. Driven by single cell/nucleus RNA sequencing (scRNA-seq) data, comprehensive atlas projects covering a wide range of organisms and tissues are currently underway. As a result, it is critical that the cell transcriptional phenotypes discovered are defined and disseminated in a consistent and concise manner. Molecular biomarkers have historically played an important role in biological research, from defining immune cell-types by surface protein expression to defining diseases by molecular drivers. Here we describe a machine learning-based marker gene selection algorithm, NS-Forest version 2.0, which leverages the non-linear attributes of random forest feature selection and a binary expression scoring approach to discover the minimal marker gene expression combinations that precisely captures the cell type identity represented in the complete scRNA-seq transcriptional profiles. The marker genes selected provide a barcode of the necessary and sufficient characteristics for semantic cell type definition and serve as useful tools for downstream biological investigation. The use of NS-Forest to identify marker genes for human brain middle temporal gyrus cell types reveals the importance of cell signaling and non-coding RNAs in neuronal cell type identity.

Single cell/nucleus RNA sequencing (scRNAseq) is emerging as an essential tool to unravel the phenotypic heterogeneity of cells in complex biological systems. While computational methods for scRNAseq cell type clustering have advanced, the ability to integrate datasets to identify common and novel cell types across experiments remains a challenge. Here, we introduce a cluster-to-cluster cell type matching method – FR-Match – that utilizes supervised feature selection for dimensionality reduction and incorporates shared information among cells to determine whether two cell type clusters share the same underlying multivariate gene expression distribution. FR-Match is benchmarked with existing cell-to-cell and cell-to-cluster cell type matching methods using both simulated and real scRNAseq data. FR-Match proved to be a stringent method that produced fewer erroneous matches of distinct cell subtypes and had the unique ability to identify novel cell phenotypes in new datasets. In silico validation demonstrated that the proposed workflow is the only self-contained algorithm that was robust to increasing numbers of true negatives (i.e. non-represented cell types). FR-Match was applied to two human brain scRNAseq datasets sampled from cortical layer 1 and full thickness middle temporal gyrus. When mapping cell types identified in specimens isolated from these overlapping human brain regions, FR-Match precisely recapitulated the laminar characteristics of matched cell type clusters, reflecting their distinct neuroanatomical distributions. An R package and Shiny application are provided at https://github.com/JCVenterInstitute/FRmatch for users to interactively explore and match scRNAseq cell type clusters with complementary visualization tools.

The neocortex is disproportionately expanded in human compared to mouse, both in its total volume relative to subcortical structures and in the proportion occupied by supragranular layers that selectively make connections within the cortex and other telencephalic structures. Single-cell transcriptomic analyses of human and mouse cortex show an increased diversity of glutamatergic neuron types in supragranular cortex in human and pronounced gradients as a function of cortical depth. To probe the functional and anatomical correlates of this transcriptomic diversity, we describe a robust Patch-seq platform using neurosurgically-resected human tissues. We characterize the morphological and physiological properties of five transcriptomically defined human glutamatergic supragranular neuron types. Three of these types have properties that are specialized compared to the more homogeneous properties of transcriptomically defined homologous mouse neuron types. The two remaining supragranular neuron types, located exclusively in deep layer 3, do not have clear mouse homologues in supragranular cortex but are transcriptionally most similar to deep layer mouse intratelencephalic-projecting neuron types. Furthermore, we reveal the transcriptomic types in deep layer 3 that express high levels of non-phosphorylated heavy chain neurofilament protein that label long-range neurons known to be selectively depleted in Alzheimer’s disease. Together, these results demonstrate the power of transcriptomic cell type classification, provide a mechanistic underpinning for increased complexity of cortical function in human cortical evolution, and implicate discrete transcriptomic cell types as selectively vulnerable in disease.

Neurons are frequently classified into distinct groups or cell types on the basis of structural, physiological, or genetic attributes. To better constrain the definition of neuronal cell types, we characterized the transcriptomes and intrinsic physiological properties of over 3,700 GABAergic mouse visual cortical neurons and reconstructed the local morphologies of 350 of those neurons. We found that most transcriptomic types (t-types) occupy specific laminar positions within mouse visual cortex, and many of those t-types exhibit consistent electrophysiological and morphological features. We observed that these properties could vary continuously between t-types, which limited the ability to predict specific t-types from other data modalities. Despite that, the data support the presence of at least 20 interneuron met-types that have congruent morphological, electrophysiological, and transcriptomic properties.HighlightsPatch-seq data obtained from >3,700 GABAergic cortical interneuronsComprehensive characterization of morpho-electric features of transcriptomic types20 interneuron met-types that have congruent properties across data modalitiesDifferent Sst met-types preferentially innervate different cortical layers

Neurons are frequently classified into distinct groups or cell types on the basis of structural, physiological, or genetic attributes. To better constrain the definition of neuronal cell types, we characterized the transcriptomes and intrinsic physiological properties of over 3,700 GABAergic mouse visual cortical neurons and reconstructed the local morphologies of 350 of those neurons. We found that most transcriptomic types (t-types) occupy specific laminar positions within mouse visual cortex, and many of those t-types exhibit consistent electrophysiological and morphological features. We observed that these properties could vary continuously between t- types, which limited the ability to predict specific t-types from other data modalities. Despite that, the data support the presence of at least 20 interneuron met-types that have congruent morphological, electrophysiological, and transcriptomic properties.

In the neocortex, subcerebral axonal projections originate largely from layer 5 (L5) extratelencephalic-projecting (ET) neurons. The unique morpho-electric properties of these neurons have been mainly described in rodents, where retrograde tracers or transgenic lines can label them. Similar labeling strategies are infeasible in the human neocortex, rendering the translational relevance of findings in rodents unclear. We leveraged the recent discovery of a transcriptomically defined L5 ET neuron type to study the properties of human L5 ET neurons in neocortical brain slices derived from neurosurgeries. Patch-seq recordings, where transcriptome, physiology, and morphology were assayed from the same cell, revealed many conserved morpho-electric properties of human and rodent L5 ET neurons. Divergent properties were often subtler than differences between L5 cell types within these two species. These data suggest a conserved function of L5 ET neurons in the neocortical hierarchy but also highlight phenotypic divergence possibly related to functional specialization of human neocortex.

To better understand the molecular and cellular differences in brain organization between human and nonhuman primates, we performed transcriptome sequencing of 16 regions of adult human, chimpanzee, and macaque brains. Integration with human single-cell transcriptomic data revealed global, regional, and cell-type–specific species expression differences in genes representing distinct functional categories. We validated and further characterized the human specificity of genes enriched in distinct cell types through histological and functional analyses, including rare subpallial-derived interneurons expressing dopamine biosynthesis genes enriched in the human striatum and absent in the nonhuman African ape neocortex. Our integrated analysis of the generated data revealed diverse molecular and cellular features of the phylogenetic reorganization of the human brain across multiple levels, with relevance for brain function and disease., Data was generated as part of the PsychENCODE Consortium, supported by MH103339, MH106934, and MH110926. Additional support was provided by NIH grants MH109904, MH106874, AG048918, DK111178, and NS092988 (National Chimpanzee Brain Resource), the Kavli Foundation, the James S. McDonnell Foundation, NSF grant BCS-1316829, MINECO BFU2014-55090-P (FEDER), Howard Hughes International Early Career, and Secretaria d’Universitats i Recerca del Departament d’Economia i Coneixement de la Generalitat de Catalunya.

SummaryViral genetic tools to target specific brain cell types in humans and non-genetic model organisms will transform basic neuroscience and targeted gene therapy. Here we used comparative epigenetics to identify thousands of human neuronal subclass-specific putative enhancers to regulate viral tools, and 34% of these were conserved in mouse. We established an AAV platform to evaluate cellular specificity of functional enhancers by multiplexed fluorescent in situ hybridization (FISH) and single cell RNA sequencing. Initial testing in mouse neocortex yields a functional enhancer discovery success rate of over 30%. We identify enhancers with specificity for excitatory and inhibitory classes and subclasses including PVALB, LAMP5, and VIP/LAMP5 cells, some of which maintain specificityin vivoorex vivoin monkey and human neocortex. Finally, functional enhancers can be proximal or distal to cellular marker genes, conserved or divergent across species, and could yield brain-wide specificity greater than the most selective marker genes.

Defining a complete set of cell types within the cortex requires reconciling disparate results achieved through diverging methodologies. To address this correspondence problem, multiple methodologies must be applied to the same cells across multiple single-cell experiments. Here we present a new approach applying spatial transcriptomics using multiplexed fluorescencein situhybridization, (mFISH) to brain tissue previously interrogated through two photon optogenetic mapping of synaptic connectivity. This approach can resolve the anatomical, transcriptomic, connectomic, electrophysiological, and morphological characteristics of single cells within the mouse cortex.

INTRODUCTION The brain is responsible for cognition, behavior, and much of what makes us uniquely human. The development of the brain is a highly complex process, and this process is reliant on precise regulation of molecular and cellular events grounded in the spatiotemporal regulation of the transcriptome. Disruption of this regulation can lead to neuropsychiatric disorders. RATIONALE The regulatory, epigenomic, and transcriptomic features of the human brain have not been comprehensively compiled across time, regions, or cell types. Understanding the etiology of neuropsychiatric disorders requires knowledge not just of endpoint differences between healthy and diseased brains but also of the developmental and cellular contexts in which these differences arise. Moreover, an emerging body of research indicates that many aspects of the development and physiology of the human brain are not well recapitulated in model organisms, and therefore it is necessary that neuropsychiatric disorders be understood in the broader context of the developing and adult human brain. RESULTS Here we describe the generation and analysis of a variety of genomic data modalities at the tissue and single-cell levels, including transcriptome, DNA methylation, and histone modifications across multiple brain regions ranging in age from embryonic development through adulthood. We observed a widespread transcriptomic transition beginning during late fetal development and consisting of sharply decreased regional differences. This reduction coincided with increases in the transcriptional signatures of mature neurons and the expression of genes associated with dendrite development, synapse development, and neuronal activity, all of which were temporally synchronous across neocortical areas, as well as myelination and oligodendrocytes, which were asynchronous. Moreover, genes including MEF2C , SATB2 , and TCF4 , with genetic associations to multiple brain-related traits and disorders, converged in a small number of modules exhibiting spatial or spatiotemporal specificity. CONCLUSION We generated and applied our dataset to document transcriptomic and epigenetic changes across human development and then related those changes to major neuropsychiatric disorders. These data allowed us to identify genes, cell types, gene coexpression modules, and spatiotemporal loci where disease risk might converge, demonstrating the utility of the dataset and providing new insights into human development and disease. Spatiotemporal dynamics of human brain development and neuropsychiatric risks. Human brain development begins during embryonic development and continues through adulthood (top). Integrating data modalities (bottom left) revealed age- and cell type–specific properties and global patterns of transcriptional dynamics, including a late fetal transition (bottom middle). We related the variation in gene expression (brown, high; purple, low) to regulatory elements in the fetal and adult brains, cell type–specific signatures, and genetic loci associated with neuropsychiatric disorders (bottom right; gray circles indicate enrichment for corresponding features among module genes). Relationships depicted in this panel do not correspond to specific observations. CBC, cerebellar cortex; STR, striatum; HIP, hippocampus; MD, mediodorsal nucleus of thalamus; AMY, amygdala.

Neurons are frequently classified into distinct types on the basis of structural, physiological, or genetic attributes. To better constrain the definition of neuronal cell types, we characterized the transcriptomes and intrinsic physiological properties of over 4,200 mouse visual cortical GABAergic interneurons and reconstructed the local morphologies of 517 of those neurons. We find that most transcriptomic types (t-types) occupy specific laminar positions within visual cortex, and, for most types, the cells mapping to a t-type exhibit consistent electrophysiological and morphological properties. These properties display both discrete and continuous variation among t-types. Through multimodal integrated analysis, we define 28 met-types that have congruent morphological, electrophysiological, and transcriptomic properties and robust mutual predictability. We identify layer-specific axon innervation pattern as a defining feature distinguishing different met-types. These met-types represent a unified definition of cortical GABAergic interneuron types, providing a systematic framework to capture existing knowledge and bridge future analyses across different modalities.

SummaryGene expression studies suggest that differential ion channel expression contributes to differences in rodent versus human neuronal physiology. We tested whether h-channels more prominently contribute to the physiological properties of human compared to mouse supragranular pyramidal neurons. Single cell/nucleus RNA sequencing revealed ubiquitous HCN1-subunit expression in excitatory neurons in human, but not mouse supragranular layers. Using patch-clamp recordings, we found stronger h-channel-related membrane properties in supragranular pyramidal neurons in human temporal cortex, compared to mouse supragranular pyramidal neurons in temporal association area. The magnitude of these differences depended upon cortical depth and was largest in pyramidal neurons in deep L3. Additionally, pharmacologically blocking h-channels produced a larger change in membrane properties in human compared to mouse neurons. Finally, using biophysical modeling, we provided evidence that h-channels promote the transfer of theta frequencies from dendrite-to-soma in human L3 pyramidal neurons. Thus, h-channels contribute to between-species differences in a fundamental neuronal property.

This protocol is a practical guide to the N-methyl-D-glucamine (NMDG) protective recovery method of brain slice preparation. Numerous recent studies have validated the utility of this method for enhancing neuronal preservation and overall brain slice viability. The implementation of this technique by early adopters has facilitated detailed investigations into brain function using diverse experimental applications and spanning a wide range of animal ages, brain regions, and cell types. Steps are outlined for carrying out the protective recovery brain slice technique using an optimized NMDG artificial cerebrospinal fluid (aCSF) media formulation and enhanced procedure to reliably obtain healthy brain slices for patch clamp electrophysiology. With this updated approach, a substantial improvement is observed in the speed and reliability of gigaohm seal formation during targeted patch clamp recording experiments while maintaining excellent neuronal preservation, thereby facilitating challenging experimental applications. Representative results are provided from multi-neuron patch clamp recording experiments to assay synaptic connectivity in neocortical brain slices prepared from young adult transgenic mice and mature adult human neurosurgical specimens. Furthermore, the optimized NMDG protective recovery method of brain slicing is compatible with both juvenile and adult animals, thus resolving a limitation of the original methodology. In summary, a single media formulation and brain slicing procedure can be implemented across various species and ages to achieve excellent viability and tissue preservation.

Cells are fundamental functional units of multicellular organisms, with different cell types playing distinct physiological roles in the body. The recent advent of single cell transcriptional profiling using RNA sequencing is producing “big data”, enabling the identification of novel human cell types at an unprecedented rate. In this review, we summarize recent work characterizing cell types in the human central nervous and immune systems using single cell and single nuclei RNA sequencing, and discuss the implications that these discoveries are having on the representation of cell types in the reference Cell Ontology (CL). We propose a method based on random forest machine learning for identifying sets of necessary and sufficient marker genes that can be used to assemble consistent and reproducible cell type definitions for incorporation into the CL. The representation of defined cell type classes and their relationships in the CL using this strategy will make the cell type classes findable, accessible, interoperable, and reusable (FAIR), allowing the CL to serve as a reference knowledgebase of information about the role that distinct cellular phenotypes play in human health and disease.

Gene expression studies suggest that differential ion channel expression contributes to differences in rodent versus human neuronal physiology. We tested whether h-channels more prominently contribute to the physiological properties of human compared to mouse supragranular pyramidal neurons. Single cell/nucleus RNA sequencing revealed ubiquitous HCN1-subunit expression in excitatory neurons in human, but not mouse supragranular layers. Using patch-clamp recordings, we found stronger h-channel-related membrane properties in supragranular pyramidal neurons in human temporal cortex, compared to mouse supragranular pyramidal neurons in temporal association area. The magnitude of these differences depended upon cortical depth and was largest in pyramidal neurons in deep L3. Additionally, pharmacologically blocking h-channels produced a larger change in membrane properties in human compared to mouse neurons. Finally, using biophysical modeling, we provided evidence that h-channels promote the transfer of theta frequencies from dendrite-to-soma in human L3 pyramidal neurons. Thus, h-channels contribute to between-species differences in a fundamental neuronal property.

Transcriptional profiling of complex tissues by RNA-sequencing of single nuclei presents some advantages over whole cell analysis. It enables unbiased cellular coverage, lack of cell isolation-based transcriptional effects, and application to archived frozen specimens. Using a well-matched pair of single-nucleus RNA-seq (snRNA-seq) and single-cell RNA-seq (scRNA-seq) SMART-Seq v4 datasets from mouse visual cortex, we demonstrate that similarly high-resolution clustering of closely related neuronal types can be achieved with both methods if intronic sequences are included in nuclear RNA-seq analysis. More transcripts are detected in individual whole cells (∼11,000 genes) than nuclei (∼7,000 genes), but the majority of genes have similar detection across cells and nuclei. We estimate that the nuclear proportion of total cellular mRNA varies from 20% to over 50% for large and small pyramidal neurons, respectively. Together, these results illustrate the high information content of nuclear RNA for characterization of cellular diversity in brain tissues.

Neocortex contains a multitude of cell types segregated into layers and functionally distinct regions. To investigate the diversity of cell types across the mouse neocortex, we analyzed 12,714 cells from the primary visual cortex (VISp), and 9,035 cells from the anterior lateral motor cortex (ALM) by deep single-cell RNA-sequencing (scRNA-seq), identifying 116 transcriptomic cell types. These two regions represent distant poles of the neocortex and perform distinct functions. We define 50 inhibitory transcriptomic cell types, all of which are shared across both cortical regions. In contrast, 49 of 52 excitatory transcriptomic types were found in either VISp or ALM, with only three present in both. By combining single cell RNA-seq and retrograde labeling, we demonstrate correspondence between excitatory transcriptomic types and their region-specific long-range target specificity. This study establishes a combined transcriptomic and projectional taxonomy of cortical cell types from functionally distinct regions of the mouse cortex.

The neocortex contains a multitude of cell types that are segregated into layers and functionally distinct areas. To investigate the diversity of cell types across the mouse neocortex, here we analysed 23,822 cells from two areas at distant poles of the mouse neocortex: the primary visual cortex and the anterior lateral motor cortex. We define 133 transcriptomic cell types by deep, single-cell RNA sequencing. Nearly all types of GABA (γ-aminobutyric acid)-containing neurons are shared across both areas, whereas most types of glutamatergic neurons were found in one of the two areas. By combining single-cell RNA sequencing and retrograde labelling, we match transcriptomic types of glutamatergic neurons to their long-range projection specificity. Our study establishes a combined transcriptomic and projectional taxonomy of cortical cell types from functionally distinct areas of the adult mouse cortex.

The structure and function of the human brain are highly stereotyped, implying a conserved molecular program responsible for its development, cellular structure and function. We applied a correlation-based metric called differential stability to assess reproducibility of gene expression patterning across 132 structures in six individual brains, revealing mesoscale genetic organization. The genes with the highest differential stability are highly biologically relevant, with enrichment for brain-related annotations, disease associations, drug targets and literature citations. Using genes with high differential stability, we identified 32 anatomically diverse and reproducible gene expression signatures, which represent distinct cell types, intracellular components and/or associations with neurodevelopmental and neurodegenerative disorders. Genes in neuron-associated compared to non-neuronal networks showed higher preservation between human and mouse; however, many diversely patterned genes displayed marked shifts in regulation between species. Finally, highly consistent transcriptional architecture in neocortex is correlated with resting state functional connectivity, suggesting a link between conserved gene expression and functionally relevant circuitry.

Developmental changes in the temporal and spatial regulation of gene expression drive the emergence of normal mature brain function, while disruptions in these processes underlie many neurodevelopmental abnormalities. To solidify our foundational knowledge of such changes in a primate brain with an extended period of postnatal maturation like in human, we investigated the whole-genome transcriptional profiles of rhesus monkey brains from birth to adulthood. We found that gene expression dynamics are largest from birth through infancy, after which gene expression profiles transition to a relatively stable state by young adulthood. Biological pathway enrichment analysis revealed that genes more highly expressed at birth are associated with cell adhesion and neuron differentiation, while genes more highly expressed in juveniles and adults are associated with cell death. Neocortex showed significantly greater differential expression over time than subcortical structures, and this trend likely reflects the protracted postnatal development of the cortex. Using network analysis, we identified 27 co-expression modules containing genes with highly correlated expression patterns that are associated with specific brain regions, ages or both. In particular, one module with high expression in neonatal cortex and striatum that decreases during infancy and juvenile development was significantly enriched for autism spectrum disorder (ASD)-related genes. This network was enriched for genes associated with axon guidance and interneuron differentiation, consistent with a disruption in the formation of functional cortical circuitry in ASD.

The stereotyped spatial architecture of the brain is both beautiful and fundamentally related to its function, extending from gross morphology to individual neuron types, where soma position, dendritic architecture, and axonal projections determine their roles in functional circuitry. Our understanding of the cell types that make up the brain is rapidly accelerating, driven in particular by recent advances in single-cell transcriptomics. However, understanding brain function, development, and disease will require linking molecular cell types to morphological, physiological, and behavioral correlates. Emerging spatially resolved transcriptomic methods promise to fill this gap by localizing molecularly defined cell types in tissues, with simultaneous detection of morphology, activity, or connectivity. Here, we review the requirements for spatial transcriptomic methods toward these goals, consider the challenges ahead, and describe promising applications.

Summary Gene expression studies suggest that differential ion channel expression contributes to differences in rodent versus human neuronal physiology. We tested whether h-channels more prominently contribute to the physiological properties of human compared to mouse supragranular pyramidal neurons. Single-cell/nucleus RNA sequencing revealed ubiquitous HCN1-subunit expression in excitatory neurons in human, but not mouse, supragranular layers. Using patch-clamp recordings, we found stronger h-channel-related membrane properties in supragranular pyramidal neurons in human temporal cortex, compared to mouse supragranular pyramidal neurons in temporal association area. The magnitude of these differences depended upon cortical depth and was largest in pyramidal neurons in deep L3. Additionally, pharmacologically blocking h-channels produced a larger change in membrane properties in human compared to mouse neurons. Finally, using biophysical modeling, we provide evidence that h-channels promote the transfer of theta frequencies from dendrite-to-soma in human L3 pyramidal neurons. Thus, h-channels contribute to between-species differences in a fundamental neuronal property.

The neurogenic potential of the subgranular zone (SGZ) of the hippocampal dentate gyrus is likely to be regulated by molecular cues arising from its complex heterogeneous cellular environment. Through transcriptome analysis using laser microdissection coupled with DNA microarrays, in combination with analysis of genome-wide in situ hybridization data, we identified 363 genes selectively enriched in adult mouse SGZ. These genes reflect expression in the different constituent cell types, including progenitor and dividing cells, immature granule cells, astrocytes, oligodendrocytes and GABAergic interneurons. Similar transcriptional profiling in the rhesus monkey dentate gyrus across postnatal development identified a highly overlapping set of SGZ-enriched genes, which can be divided based on temporal profiles to reflect maturation of glia versus granule neurons. Furthermore, we identified a neurogenesis-related gene network with decreasing postnatal expression that is highly correlated with the declining number of proliferating cells in dentate gyrus over postnatal development. Many of the genes in this network showed similar postnatal downregulation in mouse, suggesting a conservation of molecular mechanisms underlying developmental and adult neurogenesis in rodents and primates. Conditional deletion of Sox4 and Sox11, encoding two neurogenesis-related transcription factors central in this network, produces a mouse with no hippocampus, confirming the crucial role for these genes in regulating hippocampal neurogenesis.

The transcriptional underpinnings of brain development remain poorly understood, particularly in humans and closely related non-human primates. We describe a high-resolution transcriptional atlas of rhesus monkey (Macaca mulatta) brain development that combines dense temporal sampling of prenatal and postnatal periods with fine anatomical division of cortical and subcortical regions associated with human neuropsychiatric disease. Gene expression changes more rapidly before birth, both in progenitor cells and maturing neurons. Cortical layers and areas acquire adult-like molecular profiles surprisingly late in postnatal development. Disparate cell populations exhibit distinct developmental timing of gene expression, but also unexpected synchrony of processes underlying neural circuit construction including cell projection and adhesion. Candidate risk genes for neurodevelopmental disorders including primary microcephaly, autism spectrum disorder, intellectual disability, and schizophrenia show disease-specific spatiotemporal enrichment within developing neocortex. Human developmental expression trajectories are more similar to monkey than rodent, although approximately 9% of genes show human-specific regulation with evidence for prolonged maturation or neoteny compared to monkey.

Neuroanatomically precise, genome-wide maps of transcript distributions are critical resources to complement genomic sequence data and to correlate functional and genetic brain architecture. Here we describe the generation and analysis of a transcriptional atlas of the adult human brain, comprising extensive histological analysis and comprehensive microarray profiling of ~900 neuroanatomically precise subdivisions in two individuals. Transcriptional regulation varies enormously by anatomical location, with different regions and their constituent cell types displaying robust molecular signatures that are highly conserved between individuals. Analysis of differential gene expression and gene co-expression relationships demonstrates that brain-wide variation strongly reflects the distributions of major cell classes such as neurons, oligodendrocytes, astrocytes and microglia. Local neighbourhood relationships between fine anatomical subdivisions are associated with discrete neuronal subtypes and genes involved with synaptic transmission. The neocortex displays a relatively homogeneous transcriptional pattern, but with distinct features associated selectively with primary sensorimotor cortices and with enriched frontal lobe expression. Notably, the spatial topography of the neocortex is strongly reflected in its molecular topography— the closer two cortical regions, the more similar their transcriptomes. This freely accessible online data resource forms a high-resolution transcriptional baseline for neurogenetic studies of normal and abnormal human brain function.

Neuronal responses during sensory processing are influenced by both the organization of intracortical connections and the statistical features of sensory stimuli. How these intrinsic and extrinsic factors govern the activity of excitatory and inhibitory populations is unclear. Using two-photon calcium imaging in vivo and intracellular recordings in vitro, we investigated the dependencies between synaptic connectivity, feature selectivity and network activity in pyramidal cells and fast-spiking parvalbumin-expressing (PV) interneurons in mouse visual cortex. In pyramidal cell populations, patterns of neuronal correlations were largely stimulus-dependent, indicating that their responses were not strongly dominated by functionally biased recurrent connectivity. By contrast, visual stimulation only weakly modified co-activation patterns of fast-spiking PV cells, consistent with the observation that these broadly tuned interneurons received very dense and strong synaptic input from nearby pyramidal cells with diverse feature selectivities. Therefore, feedforward and recurrent network influences determine the activity of excitatory and inhibitory ensembles in fundamentally different ways.

Considerable progress has been made in understanding variations in gene sequence and expression level associated with phenotype, yet how genetic diversity translates into complex phenotypic differences remains poorly understood. Here, we examine the relationship between genetic background and spatial patterns of gene expression across seven strains of mice, providing the most extensive cellular-resolution comparative analysis of gene expression in the mammalian brain to date. Using comprehensive brainwide anatomic coverage (more than 200 brain regions), we applied in situ hybridization to analyze the spatial expression patterns of 49 genes encoding well-known pharmaceutical drug targets. Remarkably, over 50% of the genes examined showed interstrain expression variation. In addition, the variability was nonuniformly distributed across strain and neuroanatomic region, suggesting certain organizing principles. First, the degree of expression variance among strains mirrors genealogic relationships. Second, expression pattern differences were concentrated in higher-order brain regions such as the cortex and hippocampus. Divergence in gene expression patterns across the brain could contribute significantly to variations in behavior and responses to neuroactive drugs in laboratory mouse strains and may help to explain individual differences in human responsiveness to neuroactive drugs.

Although cytoarchitectonic organization of the mammalian cortex into different lamina has been well-studied, identifying the architectural differences that distinguish cortical areas from one another is more challenging. Localization of large anatomical structures is possible using magnetic resonance imaging or invasive techniques (such as anterograde or retrograde tracing), but identifying patterns in gene expression architecture is limited as gene products do not necessarily identify an immediate functional consequence of a specialized area. Expression of specific genes in the mouse and human cortex is most often identified across entire lamina, and areal patterning of expression (when it exists) is most easily differentiated on a layer-by-layer basis. Since cortical organization is defined by the expression of large sets of genes, the task of identifying individual (or groups of structures) cannot be done using individual areal markers. In this manuscript we describe a methodology for clustering gene expression correlation profiles in the C57Bl/6J mouse cortex to identify large-scale genetic relationships between layers and areas. By using the Anatomic Gene Expression Atlas (http://mouse.brain-map.org/agea/) derived from in situ hybridization data in the Allen Brain Atlas, we show that a consistent expression based organization of areal patterning in the mouse cortex exists when clustered on a laminar basis. Surface-based mapping and visualization techniques are used as a representation to clarify these relationships.

The Allen Brain Atlas (ABA, www.brain-map.org) is a genome wide, spatially registered collection of cellular resolution in situ hybridization gene expression image data of the C57Bl/6J mouse brain. Derived from the ABA, the Anatomic Gene Expression Atlas (AGEA, http://mouse.brain-map.org/agea) has demonstrated both laminar and areal spatial gene expression correlations in the mouse cortex. While the mouse cortex is lissencephalic, its curvature and substantial bending in boundary areas renders it difficult to visualize and analyze laminar versus areal effects in a rectilinear coordinate framework. In context of human and non-human primate cortex, surface-based representation has proven useful for understanding relative locations of laminar, columnar, and areal features. In this paper, we describe a methodology for constructing surface-based flatmaps of the mouse cortex that enables mapping of gene expression data from individual genes in the ABA, or probabilistic expression maps from the AGEA, to identify and visualize genetic relationships between layers and areas.

The molecular and cellular landscape of glioblastoma is highly complex and its relationship to histologic features routinely used for diagnosis is unclear. To investigate this relationship, we generated an anatomic transcriptional atlas of human glioblastoma, adopting a highly-systematized, large-scale, histology-driven approach to the characterization of anatomic features and cancer stem cell populations. The atlas of 42 tumors consists of several data modalities, including 270 transcriptomes, ∼11,500 semi-annotated pathology images registered to ∼23,000 in situ hybridization gene expression images, ∼400 MRI scans, tumor-derived cell lines and xenografts, and supporting longitudinal clinical information. We show that gene expression varies considerably by anatomic feature and cancer stem cell population, exhibiting molecular signatures that are more highly conserved within tumors than between tumors, and reflecting the cell types and microenvironment of each feature and population. These freely-accessible online data resources of the Ivy Glioblastoma Atlas Project (Ivy GAP), one for the fully-annotated anatomic transcriptional atlas (http://glioblastoma.alleninstitute.org/), and one for the detailed clinical and genomic data (http://ivygap.org/), constitute a unique platform for exploring the anatomic and genetic basis of glioblastoma at the cellular and molecular levels. This project was supported by grants from the Ben and Catherine Ivy Foundation to the Allen Institute (PI, RBP) and Swedish Neuroscience Institute (PIs, GDF and NS).

Song-learning birds and humans share independently evolved similarities in brain pathways for vocal learning that are essential for song and speech and are not found in most other species. Comparisons of brain transcriptomes of song-learning birds and humans relative to vocal nonlearners identified convergent gene expression specializations in specific song and speech brain regions of avian vocal learners and humans. The strongest shared profiles relate bird motor and striatal song-learning nuclei, respectively, with human laryngeal motor cortex and parts of the striatum that control speech production and learning. Most of the associated genes function in motor control and brain connectivity. Thus, convergent behavior and neural connectivity for a complex trait are associated with convergent specialized expression of multiple genes.

GABAergic interneurons are essential for neural circuit function, and their loss or dysfunction is implicated in human neuropsychiatric disease. In vitro methods for interneuron generation hold promise for studying human cellular and functional properties and, ultimately, for therapeutic cell replacement. Here we describe a protocol for generating cortical interneurons from hESCs and analyze the properties and maturation time course of cell types using single-cell RNA-seq. We find that the cell types produced mimic in vivo temporal patterns of neuron and glial production, with immature progenitors and neurons observed early and mature cortical neurons and glial cell types produced late. By comparing the transcriptomes of immature interneurons to those of more mature neurons, we identified genes important for human interneuron differentiation. Many of these genes were previously implicated in neurodevelopmental and neuropsychiatric disorders.