The cell type composition of the adult mouse brain revealed by single cell and spatial genomics

Building and analyzing metacells in single-cell genomics data

The advent of high-throughput single-cell genomics technologies has fundamentally transformed biological sciences. Currently, millions of cells from complex biological tissues can be phenotypically profiled across multiple modalities. The scaling of computational methods to analyze and visualize such data is a constant challenge, and tools need to be regularly updated, if not redesigned, to cope with ever-growing numbers of cells. Over the last few years, metacells have been introduced to reduce the size and complexity of single-cell genomics data while preserving biologically relevant information and improving interpretability. Here, we review recent studies that capitalize on the concept of metacells-and the many variants in nomenclature that have been used. We further outline how and when metacells should (or should not) be used to analyze single-cell genomics data and what should be considered when analyzing such data at the metacell level. To facilitate the exploration of metacells, we provide a comprehensive tutorial on the construction and analysis of metacells from single-cell RNA-seq data ( https://github.com/GfellerLab/MetacellAnalysisTutorial ) as well as a fully integrated pipeline to rapidly build, visualize and evaluate metacells with different methods ( https://github.com/GfellerLab/MetacellAnalysisToolkit ).

Strain in the Midbrain: Impact of Traumatic Brain Injury on the Central Serotonin System

Traumatic brain injury (TBI) is a pervasive public health crisis that severely impacts the quality of life of affected individuals. Like peripheral forms of trauma, TBI results from extraordinarily heterogeneous environmental forces being imparted on the cranial space, resulting in heterogeneous disease pathologies. This has made therapies for TBI notoriously difficult to develop, and currently, there are no FDA-approved pharmacotherapies specifically for the acute or chronic treatment of TBI. TBI is associated with changes in cognition and can precipitate the onset of debilitating psychiatric disorders like major depressive disorder (MDD), generalized anxiety disorder (GAD), and post-traumatic stress disorder (PTSD). Complicating these effects of TBI, FDA-approved pharmacotherapies utilized to treat these disorders often fail to reach the desired level of efficacy in the context of neurotrauma. Although a complicated association, decades of work have linked central serotonin (5-HT) neurotransmission as being involved in the etiology of a myriad of neuropsychiatric disorders, including MDD and GAD. 5-HT is a biogenic monoamine neurotransmitter that is highly conserved across scales of biology. Though the majority of 5-HT is isolated to peripheral sites such as the gastrointestinal (GI) tract, 5-HT neurotransmission within the CNS exerts exquisite control over diverse biological functions, including sleep, appetite and respiration, while simultaneously establishing normal mood, perception, and attention. Although several key studies have begun to elucidate how various forms of neurotrauma impact central 5-HT neurotransmission, a full determination of precisely how TBI disrupts the highly regulated dynamics of 5-HT neuron function and/or 5-HT neurotransmission has yet to be conceptually or experimentally resolved. The purpose of the current review is, therefore, to integrate the disparate bodies of 5-HT and TBI research and synthesize insight into how new combinatorial research regarding 5-HT neurotransmission and TBI may offer an informed perspective into the nature of TBI-induced neuropsychiatric complications.

Single-cell analysis of chromatin accessibility in the adult mouse brain

Recent advances in single-cell technologies have led to the discovery of thousands of brain cell types; however, our understanding of the gene regulatory programs in these cell types is far from complete1-4. Here we report a comprehensive atlas of candidate cis-regulatory DNA elements (cCREs) in the adult mouse brain, generated by analysing chromatin accessibility in 2.3 million individual brain cells from 117 anatomical dissections. The atlas includes approximately 1 million cCREs and their chromatin accessibility across 1,482 distinct brain cell populations, adding over 446,000 cCREs to the most recent such annotation in the mouse genome. The mouse brain cCREs are moderately conserved in the human brain. The mouse-specific cCREs-specifically, those identified from a subset of cortical excitatory neurons-are strongly enriched for transposable elements, suggesting a potential role for transposable elements in the emergence of new regulatory programs and neuronal diversity. Finally, we infer the gene regulatory networks in over 260 subclasses of mouse brain cells and develop deep-learning models to predict the activities of gene regulatory elements in different brain cell types from the DNA sequence alone. Our results provide a resource for the analysis of cell-type-specific gene regulation programs in both mouse and human brains.

The molecular landscape of neurological disorders: insights from single-cell RNA sequencing in neurology and neurosurgery

Single-cell ribonucleic acid sequencing (scRNA-seq) has emerged as a transformative technology in neurological and neurosurgical research, revolutionising our comprehension of complex neurological disorders. In brain tumours, scRNA-seq has provided valuable insights into cancer heterogeneity, the tumour microenvironment, treatment resistance, and invasion patterns. It has also elucidated the brain tri-lineage cancer hierarchy and addressed limitations of current models. Neurodegenerative diseases such as Alzheimer's disease, Parkinson's disease, and amyotrophic lateral sclerosis have been molecularly subtyped, dysregulated pathways have been identified, and potential therapeutic targets have been revealed using scRNA-seq. In epilepsy, scRNA-seq has explored the cellular and molecular heterogeneity underlying the condition, uncovering unique glial subpopulations and dysregulation of the immune system. ScRNA-seq has characterised distinct cellular constituents and responses to spinal cord injury in spinal cord diseases, as well as provided molecular signatures of various cell types and identified interactions involved in vascular remodelling. Furthermore, scRNA-seq has shed light on the molecular complexities of cerebrovascular diseases, such as stroke, providing insights into specific genes, cell-specific expression patterns, and potential therapeutic interventions. This review highlights the potential of scRNA-seq in guiding precision medicine approaches, identifying clinical biomarkers, and facilitating therapeutic discovery. However, challenges related to data analysis, standardisation, sample acquisition, scalability, and cost-effectiveness need to be addressed. Despite these challenges, scRNA-seq has the potential to transform clinical practice in neurological and neurosurgical research by providing personalised insights and improving patient outcomes.

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.

Thyroid hormone rewires cortical circuits to coordinate body-wide metabolism and exploratory drive

Animals adapt to varying environmental conditions by modifying the function of their internal organs, including the brain. To be adaptive, alterations in behavior must be coordinated with the functional state of organs throughout the body. Here we find that thyroid hormone- a prominent regulator of metabolism in many peripheral organs- activates cell-type specific transcriptional programs in anterior regions of cortex of adult mice via direct activation of thyroid hormone receptors. These programs are enriched for axon-guidance genes in glutamatergic projection neurons, synaptic regulators across both astrocytes and neurons, and pro-myelination factors in oligodendrocytes, suggesting widespread remodeling of cortical circuits. Indeed, whole-cell electrophysiology recordings revealed that thyroid hormone induces local transcriptional programs that rewire cortical neural circuits via pre-synaptic mechanisms, resulting in increased excitatory drive with a concomitant sensitization of recruited inhibition. We find that thyroid hormone bidirectionally regulates innate exploratory behaviors and that the transcriptionally mediated circuit changes in anterior cortex causally promote exploratory decision-making. Thus, thyroid hormone acts directly on adult cerebral cortex to coordinate exploratory behaviors with whole-body metabolic state.

Cell-type specific molecular signatures of aging revealed in a brain-wide transcriptomic cell-type atlas

Biological aging can be defined as a gradual loss of homeostasis across various aspects of molecular and cellular function. Aging is a complex and dynamic process which influences distinct cell types in a myriad of ways. The cellular architecture of the mammalian brain is heterogeneous and diverse, making it challenging to identify precise areas and cell types of the brain that are more susceptible to aging than others. Here, we present a high-resolution single-cell RNA sequencing dataset containing ~1.2 million high-quality single-cell transcriptomic profiles of brain cells from young adult and aged mice across both sexes, including areas spanning the forebrain, midbrain, and hindbrain. We find age-associated gene expression signatures across nearly all 130+ neuronal and non-neuronal cell subclasses we identified. We detect the greatest gene expression changes in non-neuronal cell types, suggesting that different cell types in the brain vary in their susceptibility to aging. We identify specific, age-enriched clusters within specific glial, vascular, and immune cell types from both cortical and subcortical regions of the brain, and specific gene expression changes associated with cell senescence, inflammation, decrease in new myelination, and decreased vasculature integrity. We also identify genes with expression changes across multiple cell subclasses, pointing to certain mechanisms of aging that may occur across wide regions or broad cell types of the brain. Finally, we discover the greatest gene expression changes in cell types localized to the third ventricle of the hypothalamus, including tanycytes, ependymal cells, and Tbx3+ neurons found in the arcuate nucleus that are part of the neuronal circuits regulating food intake and energy homeostasis. These findings suggest that the area surrounding the third ventricle in the hypothalamus may be a hub for aging in the mouse brain. Overall, we reveal a dynamic landscape of cell-type-specific transcriptomic changes in the brain associated with normal aging that will serve as a foundation for the investigation of functional changes in the aging process and the interaction of aging and diseases.

Single-cell analysis of the nervous system at small and large scales with instant partitions

Single-cell RNA sequencing is a new frontier across all biology, particularly in neuroscience. While powerful for answering numerous neuroscience questions, limitations in sample input size, and initial capital outlay can exclude some researchers from its application. Here, we tested a recently introduced method for scRNAseq across diverse scales and neuroscience experiments. We benchmarked against a major current scRNAseq technology and found that PIPseq performed similarly, in line with earlier benchmarking data. Across dozens of samples, PIPseq recovered many brain cell types at small and large scales (1,000-100,000 cells/sample) and was able to detect differentially expressed genes in an inflammation paradigm. Similarly, PIPseq could detect expected and new differentially expressed genes in a brain single cell suspension from a knockout mouse model; it could also detect rare, virally-la-belled cells following lentiviral targeting and gene knockdown. Finally, we used PIPseq to investigate gene expression in a nontraditional model species, the little skate (Leucoraja erinacea). In total, PIPSeq was able to detect single-cell gene expression changes across models and species, with an added benefit of large scale capture and sequencing of each sample.

Whole-brain annotation and multi-connectome cell typing quantifies circuit stereotypy in Drosophila

The fruit fly Drosophila melanogaster combines surprisingly sophisticated behaviour with a highly tractable nervous system. A large part of the fly's success as a model organism in modern neuroscience stems from the concentration of collaboratively generated molecular genetic and digital resources. As presented in our FlyWire companion paper 1 , this now includes the first full brain connectome of an adult animal. Here we report the systematic and hierarchical annotation of this ~130,000-neuron connectome including neuronal classes, cell types and developmental units (hemilineages). This enables any researcher to navigate this huge dataset and find systems and neurons of interest, linked to the literature through the Virtual Fly Brain database 2 . Crucially, this resource includes 4,552 cell types. 3,094 are rigorous consensus validations of cell types previously proposed in the hemibrain connectome 3 . In addition, we propose 1,458 new cell types, arising mostly from the fact that the FlyWire connectome spans the whole brain, whereas the hemibrain derives from a subvolume. Comparison of FlyWire and the hemibrain showed that cell type counts and strong connections were largely stable, but connection weights were surprisingly variable within and across animals. Further analysis defined simple heuristics for connectome interpretation: connections stronger than 10 unitary synapses or providing >1% of the input to a target cell are highly conserved. Some cell types showed increased variability across connectomes: the most common cell type in the mushroom body, required for learning and memory, is almost twice as numerous in FlyWire as the hemibrain. We find evidence for functional homeostasis through adjustments of the absolute amount of excitatory input while maintaining the excitation-inhibition ratio. Finally, and surprisingly, about one third of the cell types proposed in the hemibrain connectome could not yet be reliably identified in the FlyWire connectome. We therefore suggest that cell types should be defined to be robust to inter-individual variation, namely as groups of cells that are quantitatively more similar to cells in a different brain than to any other cell in the same brain. Joint analysis of the FlyWire and hemibrain connectomes demonstrates the viability and utility of this new definition. Our work defines a consensus cell type atlas for the fly brain and provides both an intellectual framework and open source toolchain for brain-scale comparative connectomics.

Structural alterations in the amygdala and impaired social incentive learning in a mouse model of a genetic variant associated with neurodevelopmental disorders

Copy number variants (CNVs) are robustly associated with psychiatric disorders and their dimensions and changes in brain structures and behavior. However, as CNVs contain many genes, the precise gene-phenotype relationship remains unclear. Although various volumetric alterations in the brains of 22q11.2 CNV carriers have been identified in humans and mouse models, it is unknown how the genes in the 22q11.2 region individually contribute to structural alterations and associated mental illnesses and their dimensions. Our previous studies have identified Tbx1 , a T-box family transcription factor encoded in 22q11.2 CNV, as a driver gene for social interaction and communication, spatial and working memory, and cognitive flexibility. However, it remains unclear how TBX1 impacts the volumes of various brain regions and their functionally linked behavioral dimensions. In this study, we used volumetric magnetic resonance imaging analysis to comprehensively evaluate brain region volumes in congenic Tbx1 heterozygous mice. Our data show that the volumes of anterior and posterior portions of the amygdaloid complex and its surrounding cortical regions were reduced in Tbx1 heterozygous mice. Moreover, we examined the behavioral consequences of an altered volume of the amygdala. Tbx1 heterozygous mice were impaired for their ability to detect the incentive value of a social partner in a task that depends on the amygdala. Our findings identify the structural basis for a specific social dimension associated with loss-of-function variants of TBX1 and 22q11.2 CNV.

Single-cell DNA Methylome and 3D Multi-omic Atlas of the Adult Mouse Brain

Cytosine DNA methylation is essential in brain development and has been implicated in various neurological disorders. A comprehensive understanding of DNA methylation diversity across the entire brain in the context of the brain's 3D spatial organization is essential for building a complete molecular atlas of brain cell types and understanding their gene regulatory landscapes. To this end, we employed optimized single-nucleus methylome (snmC-seq3) and multi-omic (snm3C-seq1) sequencing technologies to generate 301,626 methylomes and 176,003 chromatin conformation/methylome joint profiles from 117 dissected regions throughout the adult mouse brain. Using iterative clustering and integrating with companion whole-brain transcriptome and chromatin accessibility datasets, we constructed a methylation-based cell type taxonomy that contains 4,673 cell groups and 261 cross-modality-annotated subclasses. We identified millions of differentially methylated regions (DMRs) across the genome, representing potential gene regulation elements. Notably, we observed spatial cytosine methylation patterns on both genes and regulatory elements in cell types within and across brain regions. Brain-wide multiplexed error-robust fluorescence in situ hybridization (MERFISH2) data validated the association of this spatial epigenetic diversity with transcription and allowed the mapping of the DNA methylation and topology information into anatomical structures more precisely than our dissections. Furthermore, multi-scale chromatin conformation diversities occur in important neuronal genes, highly associated with DNA methylation and transcription changes. Brain-wide cell type comparison allowed us to build a regulatory model for each gene, linking transcription factors, DMRs, chromatin contacts, and downstream genes to establish regulatory networks. Finally, intragenic DNA methylation and chromatin conformation patterns predicted alternative gene isoform expression observed in a companion whole-brain SMART-seq3 dataset. Our study establishes the first brain-wide, single-cell resolution DNA methylome and 3D multi-omic atlas, providing an unparalleled resource for comprehending the mouse brain's cellular-spatial and regulatory genome diversity.