Conserved and divergent gene regulatory programs of the mammalian neocortex.

Divergence of cis-regulatory elements drives species-specific traits1, but how this manifests in the evolution of the neocortex at the molecular and cellular level remains unclear. Here we investigated the gene regulatory programs in the primary motor cortex of human, macaque, marmoset and mouse using single-cell multiomics assays, generating gene expression, chromatin accessibility, DNA methylome and chromosomal conformation profiles from a total of over 200,000 cells. From these data, we show evidence that divergence of transcription factor expression corresponds to species-specific epigenome landscapes. We find that conserved and divergent gene regulatory features are reflected in the evolution of the three-dimensional genome. Transposable elements contribute to nearly 80% of the human-specific candidate cis-regulatory elements in cortical cells. Through machine learning, we develop sequence-based predictors of candidate cis-regulatory elements in different species and demonstrate that the genomic regulatory syntax is highly preserved from rodents to primates. Finally, we show that epigenetic conservation combined with sequence similarity helps to uncover functional cis-regulatory elements and enhances our ability to interpret genetic variants contributing to neurological disease and traits.

Coexpression networks implicate human midfetal deep cortical projection neurons in the pathogenesis of autism.

Autism spectrum disorder (ASD) is a complex developmental syndrome of unknown etiology. Recent studies employing exome- and genome-wide sequencing have identified nine high-confidence ASD (hcASD) genes. Working from the hypothesis that ASD-associated mutations in these biologically pleiotropic genes will disrupt intersecting developmental processes to contribute to a common phenotype, we have attempted to identify time periods, brain regions, and cell types in which these genes converge. We have constructed coexpression networks based on the hcASD "seed" genes, leveraging a rich expression data set encompassing multiple human brain regions across human development and into adulthood. By assessing enrichment of an independent set of probable ASD (pASD) genes, derived from the same sequencing studies, we demonstrate a key point of convergence in midfetal layer 5/6 cortical projection neurons. This approach informs when, where, and in what cell types mutations in these specific genes may be productively studied to clarify ASD pathophysiology.

Large-scale cellular-resolution gene profiling in human neocortex reveals species-specific molecular signatures.

Although there have been major advances in elucidating the functional biology of the human brain, relatively little is known of its cellular and molecular organization. Here we report a large-scale characterization of the expression of ∼1,000 genes important for neural functions by in situ hybridization at a cellular resolution in visual and temporal cortices of adult human brains. These data reveal diverse gene expression patterns and remarkable conservation of each individual gene's expression among individuals (95%), cortical areas (84%), and between human and mouse (79%). A small but substantial number of genes (21%) exhibited species-differential expression. Distinct molecular signatures, comprised of genes both common between species and unique to each, were identified for each major cortical cell type. The data suggest that gene expression profile changes may contribute to differential cortical function across species, and in particular, a shift from corticosubcortical to more predominant corticocortical communications in the human brain.

A comparison of anatomic and cellular transcriptome structures across 40 human brain diseases.

Genes associated with risk for brain disease exhibit characteristic expression patterns that reflect both anatomical and cell type relationships. Brain-wide transcriptomic patterns of disease risk genes provide a molecular-based signature, based on differential co-expression, that is often unique to that disease. Brain diseases can be compared and aggregated based on the similarity of their signatures which often associates diseases from diverse phenotypic classes. Analysis of 40 common human brain diseases identifies 5 major transcriptional patterns, representing tumor-related, neurodegenerative, psychiatric and substance abuse, and 2 mixed groups of diseases affecting basal ganglia and hypothalamus. Further, for diseases with enriched expression in cortex, single-nucleus data in the middle temporal gyrus (MTG) exhibits a cell type expression gradient separating neurodegenerative, psychiatric, and substance abuse diseases, with unique excitatory cell type expression differentiating psychiatric diseases. Through mapping of homologous cell types between mouse and human, most disease risk genes are found to act in common cell types, while having species-specific expression in those types and preserving similar phenotypic classification within species. These results describe structural and cellular transcriptomic relationships of disease risk genes in the adult brain and provide a molecular-based strategy for classifying and comparing diseases, potentially identifying novel disease relationships.

Transcriptomic Perspectives on Neocortical Structure, Development, Evolution, and Disease.

The cerebral cortex is the source of our most complex cognitive capabilities and a vulnerable target of many neurological and neuropsychiatric disorders. Transcriptomics offers a new approach to understanding the cortex at the level of its underlying genetic code, and rapid technological advances have propelled this field to the high-throughput study of the complete set of transcribed genes at increasingly fine resolution to the level of individual cells. These tools have revealed features of the genetic architecture of adult cortical areas, layers, and cell types, as well as spatiotemporal patterning during development. This has allowed a fresh look at comparative anatomy as well, illustrating surprisingly large differences between mammals while at the same time revealing conservation of some features from avians to mammals. Finally, transcriptomics is fueling progress in understanding the causes of neurodevelopmental diseases such as autism, linking genetic association studies to specific molecular pathways and affected brain regions.

Oligodendrocyte precursor cells ingest axons in the mouse neocortex.

Neurons in the developing brain undergo extensive structural refinement as nascent circuits adopt their mature form. This physical transformation of neurons is facilitated by the engulfment and degradation of axonal branches and synapses by surrounding glial cells, including microglia and astrocytes. However, the small size of phagocytic organelles and the complex, highly ramified morphology of glia have made it difficult to define the contribution of these and other glial cell types to this crucial process. Here, we used large-scale, serial section transmission electron microscopy (TEM) with computational volume segmentation to reconstruct the complete 3D morphologies of distinct glial types in the mouse visual cortex, providing unprecedented resolution of their morphology and composition. Unexpectedly, we discovered that the fine processes of oligodendrocyte precursor cells (OPCs), a population of abundant, highly dynamic glial progenitors, frequently surrounded small branches of axons. Numerous phagosomes and phagolysosomes (PLs) containing fragments of axons and vesicular structures were present inside their processes, suggesting that OPCs engage in axon pruning. Single-nucleus RNA sequencing from the developing mouse cortex revealed that OPCs express key phagocytic genes at this stage, as well as neuronal transcripts, consistent with active axon engulfment. Although microglia are thought to be responsible for the majority of synaptic pruning and structural refinement, PLs were ten times more abundant in OPCs than in microglia at this stage, and these structures were markedly less abundant in newly generated oligodendrocytes, suggesting that OPCs contribute substantially to the refinement of neuronal circuits during cortical development.

Strong and reliable synaptic communication between pyramidal neurons in adult human cerebral cortex.

Synaptic transmission constitutes the primary mode of communication between neurons. It is extensively studied in rodent but not human neocortex. We characterized synaptic transmission between pyramidal neurons in layers 2 and 3 using neurosurgically resected human middle temporal gyrus (MTG, Brodmann area 21), which is part of the distributed language circuitry. We find that local connectivity is comparable with mouse layer 2/3 connections in the anatomical homologue (temporal association area), but synaptic connections in human are 3-fold stronger and more reliable (0% vs 25% failure rates, respectively). We developed a theoretical approach to quantify properties of spinous synapses showing that synaptic conductance and voltage change in human dendritic spines are 3-4-folds larger compared with mouse, leading to significant NMDA receptor activation in human unitary connections. This model prediction was validated experimentally by showing that NMDA receptor activation increases the amplitude and prolongs decay of unitary excitatory postsynaptic potentials in human but not in mouse connections. Since NMDA-dependent recurrent excitation facilitates persistent activity (supporting working memory), our data uncovers cortical microcircuit properties in human that may contribute to language processing in MTG.

FR-Match: robust matching of cell type clusters from single cell RNA sequencing data using the Friedman-Rafsky non-parametric test.

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.

Parallel RNA and DNA analysis after deep sequencing (PRDD-seq) reveals cell type-specific lineage patterns in human brain.

Elucidating the lineage relationships among different cell types is key to understanding human brain development. Here we developed parallel RNA and DNA analysis after deep sequencing (PRDD-seq), which combines RNA analysis of neuronal cell types with analysis of nested spontaneous DNA somatic mutations as cell lineage markers, identified from joint analysis of single-cell and bulk DNA sequencing by single-cell MosaicHunter (scMH). PRDD-seq enables simultaneous reconstruction of neuronal cell type, cell lineage, and sequential neuronal formation ("birthdate") in postmortem human cerebral cortex. Analysis of two human brains showed remarkable quantitative details that relate mutation mosaic frequency to clonal patterns, confirming an early divergence of precursors for excitatory and inhibitory neurons, and an "inside-out" layer formation of excitatory neurons as seen in other species. In addition our analysis allows an estimate of excitatory neuron-restricted precursors (about 10) that generate the excitatory neurons within a cortical column. Inhibitory neurons showed complex, subtype-specific patterns of neurogenesis, including some patterns of development conserved relative to mouse, but also some aspects of primate cortical interneuron development not seen in mouse. PRDD-seq can be broadly applied to characterize cell identity and lineage from diverse archival samples with single-cell resolution and in potentially any developmental or disease condition.

Hybridization-based in situ sequencing (HybISS) for spatially resolved transcriptomics in human and mouse brain tissue.

Visualization of the transcriptome in situ has proven to be a valuable tool in exploring single-cell RNA-sequencing data, providing an additional spatial dimension to investigate multiplexed gene expression, cell types, disease architecture or even data driven discoveries. In situ sequencing (ISS) method based on padlock probes and rolling circle amplification has been used to spatially resolve gene transcripts in tissue sections of various origins. Here, we describe the next iteration of ISS, HybISS, hybridization-based in situ sequencing. Modifications in probe design allows for a new barcoding system via sequence-by-hybridization chemistry for improved spatial detection of RNA transcripts. Due to the amplification of probes, amplicons can be visualized with standard epifluorescence microscopes for high-throughput efficiency and the new sequencing chemistry removes limitations bound by sequence-by-ligation chemistry of ISS. HybISS design allows for increased flexibility and multiplexing, increased signal-to-noise, all without compromising throughput efficiency of imaging large fields of view. Moreover, the current protocol is demonstrated to work on human brain tissue samples, a source that has proven to be difficult to work with image-based spatial analysis techniques. Overall, HybISS technology works as a targeted amplification detection method for improved spatial transcriptomic visualization, and importantly, with an ease of implementation.

Cell type discovery using single-cell transcriptomics: implications for ontological representation.

Cells are fundamental function 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, which 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 being identified by high-throughput/high-content technologies 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.

The complete genome sequence of a Neanderthal from the Altai Mountains.

We present a high-quality genome sequence of a Neanderthal woman from Siberia. We show that her parents were related at the level of half-siblings and that mating among close relatives was common among her recent ancestors. We also sequenced the genome of a Neanderthal from the Caucasus to low coverage. An analysis of the relationships and population history of available archaic genomes and 25 present-day human genomes shows that several gene flow events occurred among Neanderthals, Denisovans and early modern humans, possibly including gene flow into Denisovans from an unknown archaic group. Thus, interbreeding, albeit of low magnitude, occurred among many hominin groups in the Late Pleistocene. In addition, the high-quality Neanderthal genome allows us to establish a definitive list of substitutions that became fixed in modern humans after their separation from the ancestors of Neanderthals and Denisovans.

Functional organization of the hippocampal longitudinal axis.

The precise functional role of the hippocampus remains a topic of much debate. The dominant view is that the dorsal (or posterior) hippocampus is implicated in memory and spatial navigation and the ventral (or anterior) hippocampus mediates anxiety-related behaviours. However, this 'dichotomy view' may need revision. Gene expression studies demonstrate multiple functional domains along the hippocampal long axis, which often exhibit sharply demarcated borders. By contrast, anatomical studies and electrophysiological recordings in rodents suggest that the long axis is organized along a gradient. Together, these observations suggest a model in which functional long-axis gradients are superimposed on discrete functional domains. This model provides a potential framework to explain and test the multiple functions ascribed to the hippocampus.

An anatomically comprehensive atlas of the adult human brain transcriptome.

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.

Cell type discovery and representation in the era of high-content single cell phenotyping.

BACKGROUND: A fundamental characteristic of multicellular organisms is the specialization of functional cell types through the process of differentiation. These specialized cell types not only characterize the normal functioning of different organs and tissues, they can also be used as cellular biomarkers of a variety of different disease states and therapeutic/vaccine responses. In order to serve as a reference for cell type representation, the Cell Ontology has been developed to provide a standard nomenclature of defined cell types for comparative analysis and biomarker discovery. Historically, these cell types have been defined based on unique cellular shapes and structures, anatomic locations, and marker protein expression. However, we are now experiencing a revolution in cellular characterization resulting from the application of new high-throughput, high-content cytometry and sequencing technologies. The resulting explosion in the number of distinct cell types being identified is challenging the current paradigm for cell type definition in the Cell Ontology. RESULTS: In this paper, we provide examples of state-of-the-art cellular biomarker characterization using high-content cytometry and single cell RNA sequencing, and present strategies for standardized cell type representations based on the data outputs from these cutting-edge technologies, including "context annotations" in the form of standardized experiment metadata about the specimen source analyzed and marker genes that serve as the most useful features in machine learning-based cell type classification models. We also propose a statistical strategy for comparing new experiment data to these standardized cell type representations. CONCLUSION: The advent of high-throughput/high-content single cell technologies is leading to an explosion in the number of distinct cell types being identified. It will be critical for the bioinformatics community to develop and adopt data standard conventions that will be compatible with these new technologies and support the data representation needs of the research community. The proposals enumerated here will serve as a useful starting point to address these challenges.

Spatiotemporal dynamics of the postnatal developing primate brain transcriptome.

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.

Patches of disorganization in the neocortex of children with autism.

BACKGROUND: Autism involves early brain overgrowth and dysfunction, which is most strongly evident in the prefrontal cortex. As assessed on pathological analysis, an excess of neurons in the prefrontal cortex among children with autism signals a disturbance in prenatal development and may be concomitant with abnormal cell type and laminar development. METHODS: To systematically examine neocortical architecture during the early years after the onset of autism, we used RNA in situ hybridization with a panel of layer- and cell-type-specific molecular markers to phenotype cortical microstructure. We assayed markers for neurons and glia, along with genes that have been implicated in the risk of autism, in prefrontal, temporal, and occipital neocortical tissue from postmortem samples obtained from children with autism and unaffected children between the ages of 2 and 15 years. RESULTS: We observed focal patches of abnormal laminar cytoarchitecture and cortical disorganization of neurons, but not glia, in prefrontal and temporal cortical tissue from 10 of 11 children with autism and from 1 of 11 unaffected children. We observed heterogeneity between cases with respect to cell types that were most abnormal in the patches and the layers that were most affected by the pathological features. No cortical layer was uniformly spared, with the clearest signs of abnormal expression in layers 4 and 5. Three-dimensional reconstruction of layer markers confirmed the focal geometry and size of patches. CONCLUSIONS: In this small, explorative study, we found focal disruption of cortical laminar architecture in the cortexes of a majority of young children with autism. Our data support a probable dysregulation of layer formation and layer-specific neuronal differentiation at prenatal developmental stages. (Funded by the Simons Foundation and others.).

Conserved molecular signatures of neurogenesis in the hippocampal subgranular zone of rodents and primates.

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.