Transcriptomic diversity of cell types across the adult human brain

Chromatin accessibility during human first-trimester neurodevelopment

The human brain develops through a tightly organized cascade of patterning events, induced by transcription factor expression and changes in chromatin accessibility. Although gene expression across the developing brain has been described at single-cell resolution1, similar atlases of chromatin accessibility have been primarily focused on the forebrain2-4. Here we describe chromatin accessibility and paired gene expression across the entire developing human brain during the first trimester (6-13 weeks after conception). We defined 135 clusters and used multiomic measurements to link candidate cis-regulatory elements to gene expression. The number of accessible regions increased both with age and along neuronal differentiation. Using a convolutional neural network, we identified putative functional transcription factor-binding sites in enhancers characterizing neuronal subtypes. We applied this model to cis-regulatory elements linked to ESRRB to elucidate its activation mechanism in the Purkinje cell lineage. Finally, by linking disease-associated single nucleotide polymorphisms to cis-regulatory elements, we validated putative pathogenic mechanisms in several diseases and identified midbrain-derived GABAergic neurons as being the most vulnerable to major depressive disorder-related mutations. Our findings provide a more detailed view of key gene regulatory mechanisms underlying the emergence of brain cell types during the first trimester and a comprehensive reference for future studies related to human neurodevelopment.

Multimodal cortical neuronal cell type classification

Since more than a century, neuroscientists have distinguished excitatory (glutamatergic) neurons with long-distance projections from inhibitory (GABAergic) neurons with local projections and established layer-dependent schemes for the ~ 80% excitatory (principal) cells as well as the ~ 20% inhibitory neurons. Whereas, in the early days, mainly morphological criteria were used to define cell types, later supplemented by electrophysiological and neurochemical properties, nowadays. single-cell transcriptomics is the method of choice for cell type classification. Bringing recent insight together, we conclude that despite all established layer- and area-dependent differences, there is a set of reliably identifiable cortical cell types that were named (among others) intratelencephalic (IT), extratelencephalic (ET), and corticothalamic (CT) for the excitatory cells, which altogether comprise ~ 56 transcriptomic cell types (t-types). By the same means, inhibitory neurons were subdivided into parvalbumin (PV), somatostatin (SST), vasoactive intestinal polypeptide (VIP), and "other (i.e. Lamp5/Sncg)" subpopulations, which altogether comprise ~ 60 t-types. The coming years will show which t-types actually translate into "real" cell types that show a common set of multimodal features, including not only transcriptome but also physiology and morphology as well as connectivity and ultimately function. Only with the better knowledge of clear-cut cell types and experimental access to them, we will be able to reveal their specific functions, a task which turned out to be difficult in a part of the brain being so much specialized for cognition as the cerebral cortex.

Molecular and cellular mechanisms of selective vulnerability in neurodegenerative diseases

The selective vulnerability of specific neuronal subtypes is a hallmark of neurodegenerative diseases. In this Review, I summarize our current understanding of the brain regions and cell types that are selectively vulnerable in different neurodegenerative diseases and describe the proposed underlying cell-autonomous and non-cell-autonomous mechanisms. I highlight how recent methodological innovations - including single-cell transcriptomics, CRISPR-based screens and human cell-based models of disease - are enabling new breakthroughs in our understanding of selective vulnerability. An understanding of the molecular mechanisms that determine selective vulnerability and resilience would shed light on the key processes that drive neurodegeneration and point to potential therapeutic strategies to protect vulnerable cell populations.

Discovery of cell type classification marker genes from single cell RNA sequencing data using NS-Forest

The use of single-cell transcriptomic technologies that quantitively describe cell transcriptional phenotypes using single cell/nucleus RNA sequencing (scRNA-seq) is revolutionizing our understanding of cell biology, leading to new insights in cell type identification, disease mechanisms, and drug development. The tremendous growth in scRNA-seq data has posed new challenges in eYiciently characterizing data-driven cell types and identifying quantifiable marker genes for cell type classification. The use of machine learning and explainable artificial intelligence has emerged as an eYective approach to study large-scale scRNA-seq data. NS-Forest is a random forest machine learning-based algorithm that aims to provide a scalable data-driven solution to identify minimum combinations of necessary and suYicient marker genes that capture cell type identity with maximum classification accuracy. Here, we describe the latest version, NS-Forest version 4.0 and its companion Python package ( https://github.com/JCVenterInstitute/NSForest ), with several enhancements, to select marker gene combinations that exhibit selective expression patterns among closely related cell types and more eYiciently perform marker gene selection for large-scale scRNA-seq data atlases with millions of cells. By modularizing the final decision tree step, NS-Forest v4.0 can be used to compare the performance of user-defined marker genes with the NS-Forest computationally-derived marker genes based on the decision tree classifiers. To quantify how well the identified markers exhibit the desired pattern of being exclusively expressed at high levels within their target cell types, we introduce the On-Target Fraction metric that ranges from 0 to1, with a metric of 1 given to markers that are only expressed within their target cell types and not in cells of any other cell types. We have applied NS-Forest v4.0 on scRNA-seq datasets from three human organs, including the brain, kidney, and lung. We observe that NS-Forest v4.0 outperforms previous versions on its ability to identify markers with higher On-Target Fraction values for closely related cell types and outperforms other marker gene selection approaches on the classification performance with significantly higher F-beta scores.

Diverse and asymmetric patterns of single-neuron projectome in regulating interhemispheric connectivity

The corpus callosum, historically considered primarily for homotopic connections, supports many heterotopic connections, indicating complex interhemispheric connectivity. Understanding this complexity is crucial yet challenging due to diverse cell-specific wiring patterns. Here, we utilized public AAV bulk tracing and single-neuron tracing data to delineate the anatomical connection patterns of mouse brains and conducted wide-field calcium imaging to assess functional connectivity across various brain states in male mice. The single-neuron data uncovered complex and dense interconnected patterns, particularly for interhemispheric-heterotopic connections. We proposed a metric "heterogeneity" to quantify the complexity of the connection patterns. Computational modeling of these patterns suggested that the heterogeneity of upstream projections impacted downstream homotopic functional connectivity. Furthermore, higher heterogeneity observed in interhemispheric-heterotopic projections would cause lower strength but higher stability in functional connectivity than their intrahemispheric counterparts. These findings were corroborated by our wide-field functional imaging data, underscoring the important role of heterotopic-projection heterogeneity in interhemispheric communication.

Precision and Accuracy of Single-Cell/Nuclei RNA Sequencing Data

Single-cell/nuclei RNA sequencing (sc/snRNA-Seq) is widely used for profiling cell-type gene expressions in biomedical research. An important but underappreciated issue is the quality of sc/snRNA-Seq data that would impact the reliability of downstream analyses. Here we evaluated the precision and accuracy in 18 sc/snRNA-Seq datasets. The precision was assessed on data from human brain studies with a total of 3,483,905 cells from 297 individuals, by utilizing technical replicates. The accuracy was evaluated with sample-matched scRNA-Seq and pooled-cell RNA-Seq data of cultured mononuclear phagocytes from four species. The results revealed low precision and accuracy at the single-cell level across all evaluated data. Cell number and RNA quality were highlighted as two key factors determining the expression precision, accuracy, and reproducibility of differential expression analysis in sc/snRNA-Seq. This study underscores the necessity of sequencing enough high-quality cells per cell type per individual, preferably in the hundreds, to mitigate noise in expression quantification.

Unraveling the enigma: housekeeping gene Ugt1a7c as a universal biomarker for microglia

Background

Microglia, brain resident macrophages, play multiple roles in maintaining homeostasis, including immunity, surveillance, and protecting the central nervous system through their distinct activation processes. Identifying all types of microglia-driven populations is crucial due to the presence of various phenotypes that differ based on developmental stages or activation states. During embryonic development, the E8.5 yolk sac contains erythromyeloid progenitors that go through different growth phases, eventually resulting in the formation of microglia. In addition, microglia are present in neurological diseases as a diverse population. So far, no individual biomarker for microglia has been discovered that can accurately identify and monitor their development and attributes.

Summary

Here, we highlight the newly defined biomarker of mouse microglia, UGT1A7C, which exhibits superior stability in expression during microglia development and activation compared to other known microglia biomarkers. The UGT1A7C sensing chemical probe labels all microglia in the 3xTG AD mouse model. The expression of Ugt1a7c is stable during development, with only a 4-fold variation, while other microglia biomarkers, such as Csf1r and Cx3cr1, exhibit at least a 10-fold difference. The UGT1A7C expression remains constant throughout its lifespan. In addition, the expression and activity of UGT1A7C are the same in response to different types of inflammatory activators' treatment in vitro.

Conclusion

We propose employing UGT1A7C as the representative biomarker for microglia, irrespective of their developmental state, age, or activation status. Using UGT1A7C can reduce the requirement for using multiple biomarkers, enhance the precision of microglia analysis, and even be utilized as a standard for gene/protein expression.

Benchmark of cellular deconvolution methods using a multi-assay reference dataset from postmortem human prefrontal cortex

Background

Cellular deconvolution of bulk RNA-sequencing (RNA-seq) data using single cell or nuclei RNA-seq (sc/snRNA-seq) reference data is an important strategy for estimating cell type composition in heterogeneous tissues, such as human brain. Computational methods for deconvolution have been developed and benchmarked against simulated data, pseudobulked sc/snRNA-seq data, or immunohistochemistry reference data. A major limitation in developing improved deconvolution algorithms has been the lack of integrated datasets with orthogonal measurements of gene expression and estimates of cell type proportions on the same tissue sample. Deconvolution algorithm performance has not yet been evaluated across different RNA extraction methods (cytosolic, nuclear, or whole cell RNA), different library preparation types (mRNA enrichment vs. ribosomal RNA depletion), or with matched single cell reference datasets.

Results

A rich multi-assay dataset was generated in postmortem human dorsolateral prefrontal cortex (DLPFC) from 22 tissue blocks. Assays included spatially-resolved transcriptomics, snRNA-seq, bulk RNA-seq (across six library/extraction RNA-seq combinations), and RNAScope/Immunofluorescence (RNAScope/IF) for six broad cell types. The Mean Ratio method, implemented in the DeconvoBuddies R package, was developed for selecting cell type marker genes. Six computational deconvolution algorithms were evaluated in DLPFC and predicted cell type proportions were compared to orthogonal RNAScope/IF measurements.

Conclusions

Bisque and hspe were the most accurate methods, were robust to differences in RNA library types and extractions. This multi-assay dataset showed that cell size differences, marker genes differentially quantified across RNA libraries, and cell composition variability in reference snRNA-seq impact the accuracy of current deconvolution methods.

Multimodal Nature of the Single-cell Primate Brain Atlas: Morphology, Transcriptome, Electrophysiology, and Connectivity

Primates exhibit complex brain structures that augment cognitive function. The neocortex fulfills high-cognitive functions through billions of connected neurons. These neurons have distinct transcriptomic, morphological, and electrophysiological properties, and their connectivity principles vary. These features endow the primate brain atlas with a multimodal nature. The recent integration of next-generation sequencing with modified patch-clamp techniques is revolutionizing the way to census the primate neocortex, enabling a multimodal neuronal atlas to be established in great detail: (1) single-cell/single-nucleus RNA-seq technology establishes high-throughput transcriptomic references, covering all major transcriptomic cell types; (2) patch-seq links the morphological and electrophysiological features to the transcriptomic reference; (3) multicell patch-clamp delineates the principles of local connectivity. Here, we review the applications of these technologies in the primate neocortex and discuss the current advances and tentative gaps for a comprehensive understanding of the primate neocortex.

Morphological heterogeneity of neurons in the human central amygdaloid nucleus

The central amygdaloid nucleus (CeA) has an ancient phylogenetic development and functions relevant for animal survival. Local cells receive intrinsic amygdaloidal information that codes emotional stimuli of fear, integrate them, and send cortical and subcortical output projections that prompt rapid visceral and social behavior responses. We aimed to describe the morphology of the neurons that compose the human CeA (N = 8 adult men). Cells within CeA coronal borders were identified using the thionine staining and were further analyzed using the "single-section" Golgi method followed by open-source software procedures for two-dimensional and three-dimensional image reconstructions. Our results evidenced varied neuronal cell body features, number and thickness of primary shafts, dendritic branching patterns, and density and shape of dendritic spines. Based on these criteria, we propose the existence of 12 morphologically different spiny neurons in the human CeA and discuss the variability in the dendritic architecture within cellular types, including likely interneurons. Some dendritic shafts were long and straight, displayed few collaterals, and had planar radiation within the coronal neuropil volume. Most of the sampled neurons showed a few to moderate density of small stubby/wide spines. Long spines (thin and mushroom) were observed occasionally. These novel data address the synaptic processing and plasticity in the human CeA. Our morphological description can be combined with further transcriptomic, immunohistochemical, and electrophysiological/connectional approaches. It serves also to investigate how neurons are altered in neurological and psychiatric disorders with hindered emotional perception, in anxiety, following atrophy in schizophrenia, and along different stages of Alzheimer's disease.

Single-cell and spatial transcriptomics: Bridging current technologies with long-read sequencing

Single-cell technologies have transformed biomedical research over the last decade, opening up new possibilities for understanding cellular heterogeneity, both at the genomic and transcriptomic level. In addition, more recent developments of spatial transcriptomics technologies have made it possible to profile cells in their tissue context. In parallel, there have been substantial advances in sequencing technologies, and the third generation of methods are able to produce reads that are tens of kilobases long, with error rates matching the second generation short reads. Long reads technologies make it possible to better map large genome rearrangements and quantify isoform specific abundances. This further improves our ability to characterize functionally relevant heterogeneity. Here, we show how researchers have begun to combine single-cell, spatial transcriptomics, and long-read technologies, and how this is resulting in powerful new approaches to profiling both the genome and the transcriptome. We discuss the achievements so far, and we highlight remaining challenges and opportunities.

Microglial phagolysosome dysfunction and altered neural communication amplify phenotypic severity in Prader-Willi Syndrome with larger deletion

Prader-Willi Syndrome (PWS) is a rare neurodevelopmental disorder of genetic etiology, characterized by paternal deletion of genes located at chromosome 15 in 70% of cases. Two distinct genetic subtypes of PWS deletions are characterized, where type I (PWS T1) carries four extra haploinsufficient genes compared to type II (PWS T2). PWS T1 individuals display more pronounced physiological and cognitive abnormalities than PWS T2, yet the exact neuropathological mechanisms behind these differences remain unclear. Our study employed postmortem hypothalamic tissues from PWS T1 and T2 individuals, conducting transcriptomic analyses and cell-specific protein profiling in white matter, neurons, and glial cells to unravel the cellular and molecular basis of phenotypic severity in PWS sub-genotypes. In PWS T1, key pathways for cell structure, integrity, and neuronal communication are notably diminished, while glymphatic system activity is heightened compared to PWS T2. The microglial defect in PWS T1 appears to stem from gene haploinsufficiency, as global and myeloid-specific Cyfip1 haploinsufficiency in murine models demonstrated. Our findings emphasize microglial phagolysosome dysfunction and altered neural communication as crucial contributors to the severity of PWS T1's phenotype.

ICARUS v3, a massively scalable web server for single-cell RNA-seq analysis of millions of cells

MOTIVATION

In recent years, improvements in throughput of single-cell RNA-seq have resulted in a significant increase in the number of cells profiled. The generation of single-cell RNA-seq datasets comprising >1 million cells is becoming increasingly common, giving rise to demands for more efficient computational workflows.

RESULTS

We present an update to our single-cell RNA-seq analysis web server application, ICARUS (available at https://launch.icarus-scrnaseq.cloud.edu.au) that allows effective analysis of large-scale single-cell RNA-seq datasets. ICARUS v3 utilizes the geometric cell sketching method to subsample cells from the overall dataset for dimensionality reduction and clustering that can be then projected to the large dataset. We then extend this functionality to select a representative subset of cells for downstream data analysis applications including differential expression analysis, gene co-expression network construction, gene regulatory network construction, trajectory analysis, cell-cell communication inference, and cell cluster associations to GWAS traits. We demonstrate analysis of single-cell RNA-seq datasets using ICARUS v3 of 1.3 million cells completed within the hour.

AVAILABILITY AND IMPLEMENTATION

ICARUS is available at https://launch.icarus-scrnaseq.cloud.edu.au.

Temporal patterning of the vertebrate developing neural tube

The chronologically ordered generation of distinct cell types is essential for the establishment of neuronal diversity and the formation of neuronal circuits. Recently, single-cell transcriptomic analyses of various areas of the developing vertebrate nervous system have provided evidence for the existence of a shared temporal patterning program that partitions neurons based on the timing of neurogenesis. In this review, I summarize the findings that lead to the proposal of this shared temporal program before focusing on the developing spinal cord to discuss how temporal patterning in general and this program specifically contributes to the ordered formation of neuronal circuits.

Heterogeneity of layer 4 in visual areas of rhesus macaque cortex

Recently, single-cell RNA-sequencing (scRNA-seq) has enabled unprecedented insights to the cellular landscape of the brains of many different species, among them the rhesus macaque as a key animal model. Building on previous, broader surveys of the macaque brain, we closely examined five immediately neighboring areas within the visual cortex of the rhesus macaque: V1, V2, V4, MT and TEO. To facilitate this, we first devised a novel pipeline for brain spatial archive - the BrainSPACE - which enabled robust archiving and sampling from the whole unfixed brain. SnRNA-sequencing of ~100,000 nuclei from visual areas V1 and V4 revealed conservation within the GABAergic neuron subtypes, while seven and one distinct principle neuron subtypes were detected in V1 and V4, respectively, all most likely located in layer 4. Moreover, using small molecule fluorescence in situ hybridization, we identified cell type density gradients across V1, V2, V4, MT, and TEO appearing to reflect the visual hierarchy. These findings demonstrate an association between the clear areal specializations among neighboring areas with the hierarchical levels within the visual cortex of the rhesus macaque.

Whole organism snRNA-seq reveals systemic peripheral changes in Alzheimer's Disease fly models

Peripheral tissues become disrupted in Alzheimer's Disease (AD). However, a comprehensive understanding of how the expression of AD-associated toxic proteins, Aβ42 and Tau, in neurons impacts the periphery is lacking. Using Drosophila, a prime model organism for studying aging and neurodegeneration, we generated the Alzheimer's Disease Fly Cell Atlas (AD-FCA): whole-organism single-nucleus transcriptomes of 219 cell types from adult flies neuronally expressing human Aβ42 or Tau. In-depth analyses and functional data reveal impacts on peripheral sensory neurons by Aβ42 and on various non-neuronal peripheral tissues by Tau, including the gut, fat body, and reproductive system. This novel AD atlas provides valuable insights into potential biomarkers and the intricate interplay between the nervous system and peripheral tissues in response to AD-associated proteins.

A spatially-resolved transcriptional atlas of the murine dorsal pons at single-cell resolution

The "dorsal pons", or "dorsal pontine tegmentum" (dPnTg), is part of the brainstem. It is a complex, densely packed region whose nuclei are involved in regulating many vital functions. Notable among them are the parabrachial nucleus, the Kölliker Fuse, the Barrington nucleus, the locus coeruleus, and the dorsal, laterodorsal, and ventral tegmental nuclei. In this study, we applied single-nucleus RNA-seq (snRNA-seq) to resolve neuronal subtypes based on their unique transcriptional profiles and then used multiplexed error robust fluorescence in situ hybridization (MERFISH) to map them spatially. We sampled ~1 million cells across the dPnTg and defined the spatial distribution of over 120 neuronal subtypes. Our analysis identified an unpredicted high transcriptional diversity in this region and pinpointed the unique marker genes of many neuronal subtypes. We also demonstrated that many neuronal subtypes are transcriptionally similar between humans and mice, enhancing this study's translational value. Finally, we developed a freely accessible, GPU and CPU-powered dashboard ( http://harvard.heavy.ai:6273/ ) that combines interactive visual analytics and hardware-accelerated SQL into a data science framework to allow the scientific community to query and gain insights into the data.

Drug targeting in psychiatric disorders - how to overcome the loss in translation?

In spite of major efforts and investment in development of psychiatric drugs, many clinical trials have failed in recent decades, and clinicians still prescribe drugs that were discovered many years ago. Although multiple reasons have been discussed for the drug development deadlock, we focus here on one of the major possible biological reasons: differences between the characteristics of drug targets in preclinical models and the corresponding targets in patients. Importantly, based on technological advances in single-cell analysis, we propose here a framework for the use of available and newly emerging knowledge from single-cell and spatial omics studies to evaluate and potentially improve the translational predictivity of preclinical models before commencing preclinical and, in particular, clinical studies. We believe that these recommendations will improve preclinical models and the ability to assess drugs in clinical trials, reducing failure rates in expensive late-stage trials and ultimately benefitting psychiatric drug discovery and development.

Crym-positive striatal astrocytes gate perseverative behaviour

Astrocytes are heterogeneous glial cells of the central nervous system1-3. However, the physiological relevance of astrocyte diversity for neural circuits and behaviour remains unclear. Here we show that a specific population of astrocytes in the central striatum expresses μ-crystallin (encoded by Crym in mice and CRYM in humans) that is associated with several human diseases, including neuropsychiatric disorders4-7. In adult mice, reducing the levels of μ-crystallin in striatal astrocytes through CRISPR-Cas9-mediated knockout of Crym resulted in perseverative behaviours, increased fast synaptic excitation in medium spiny neurons and dysfunctional excitatory-inhibitory synaptic balance. Increased perseveration stemmed from the loss of astrocyte-gated control of neurotransmitter release from presynaptic terminals of orbitofrontal cortex-striatum projections. We found that perseveration could be remedied using presynaptic inhibitory chemogenetics8, and that this treatment also corrected the synaptic deficits. Together, our findings reveal converging molecular, synaptic, circuit and behavioural mechanisms by which a molecularly defined and allocated population of striatal astrocytes gates perseveration phenotypes that accompany neuropsychiatric disorders9-12. Our data show that Crym-positive striatal astrocytes have key biological functions within the central nervous system, and uncover astrocyte-neuron interaction mechanisms that could be targeted in treatments for perseveration.

scGPT: toward building a foundation model for single-cell multi-omics using generative AI

Generative pretrained models have achieved remarkable success in various domains such as language and computer vision. Specifically, the combination of large-scale diverse datasets and pretrained transformers has emerged as a promising approach for developing foundation models. Drawing parallels between language and cellular biology (in which texts comprise words; similarly, cells are defined by genes), our study probes the applicability of foundation models to advance cellular biology and genetic research. Using burgeoning single-cell sequencing data, we have constructed a foundation model for single-cell biology, scGPT, based on a generative pretrained transformer across a repository of over 33 million cells. Our findings illustrate that scGPT effectively distills critical biological insights concerning genes and cells. Through further adaptation of transfer learning, scGPT can be optimized to achieve superior performance across diverse downstream applications. This includes tasks such as cell type annotation, multi-batch integration, multi-omic integration, perturbation response prediction and gene network inference.

Functional 3'-UTR Variants Identify Regulatory Mechanisms Impacting Alcohol Use Disorder and Related Traits

Although genome-wide association studies (GWAS) have identified loci associated with alcohol consumption and alcohol use disorder (AUD), they do not identify which variants are functional. To approach this, we evaluated the impact of variants in 3' untranslated regions (3'-UTRs) of genes in loci associated with substance use and neurological disorders using a massively parallel reporter assay (MPRA) in neuroblastoma and microglia cells. Functionally impactful variants explained a higher proportion of heritability of alcohol traits than non-functional variants. We identified genes whose 3'UTR activities are associated with AUD and alcohol consumption by combining variant effects from MPRA with GWAS results. We examined their effects by evaluating gene expression after CRISPR inhibition of neuronal cells and stratifying brain tissue samples by MPRA-derived 3'-UTR activity. A pathway analysis of differentially expressed genes identified inflammation response pathways. These analyses suggest that variation in response to inflammation contributes to the propensity to increase alcohol consumption.

A deep-learning strategy to identify cell types across species from high-density extracellular recordings

High-density probes allow electrophysiological recordings from many neurons simultaneously across entire brain circuits but fail to determine each recorded neuron's cell type. Here, we develop a strategy to identify cell types from extracellular recordings in awake animals, opening avenues to unveil the computational roles of neurons with distinct functional, molecular, and anatomical properties. We combine optogenetic activation and pharmacology using the cerebellum as a testbed to generate a curated ground-truth library of electrophysiological properties for Purkinje cells, molecular layer interneurons, Golgi cells, and mossy fibers. We train a semi-supervised deep-learning classifier that predicts cell types with greater than 95% accuracy based on waveform, discharge statistics, and layer of the recorded neuron. The classifier's predictions agree with expert classification on recordings using different probes, in different laboratories, from functionally distinct cerebellar regions, and across animal species. Our approach provides a general blueprint for cell-type identification from extracellular recordings across the brain.

Connecting genomic results for psychiatric disorders to human brain cell types and regions reveals convergence with functional connectivity

Understanding the temporal and spatial brain locations etiological for psychiatric disorders is essential for targeted neurobiological research. Integration of genomic insights from genome-wide association studies with single-cell transcriptomics is a powerful approach although past efforts have necessarily relied on mouse atlases. Leveraging a comprehensive atlas of the adult human brain, we prioritized cell types via the enrichment of SNP-heritabilities for brain diseases, disorders, and traits, progressing from individual cell types to brain regions. Our findings highlight specific neuronal clusters significantly enriched for the SNP-heritabilities for schizophrenia, bipolar disorder, and major depressive disorder along with intelligence, education, and neuroticism. Extrapolation of cell-type results to brain regions reveals important patterns for schizophrenia with distinct subregions in the hippocampus and amygdala exhibiting the highest significance. Cerebral cortical regions display similar enrichments despite the known prefrontal dysfunction in those with schizophrenia highlighting the importance of subcortical connectivity. Using functional MRI connectivity from cases with schizophrenia and neurotypical controls, we identified brain networks that distinguished cases from controls that also confirmed involvement of the central and lateral amygdala, hippocampal body, and prefrontal cortex. Our findings underscore the value of single-cell transcriptomics in decoding the polygenicity of psychiatric disorders and offer a promising convergence of genomic, transcriptomic, and brain imaging modalities toward common biological targets.

Multimodal profiling reveals tissue-directed signatures of human immune cells altered with age

The immune system comprises multiple cell lineages and heterogeneous subsets found in blood and tissues throughout the body. While human immune responses differ between sites and over age, the underlying sources of variation remain unclear as most studies are limited to peripheral blood. Here, we took a systems approach to comprehensively profile RNA and surface protein expression of over 1.25 million immune cells isolated from blood, lymphoid organs, and mucosal tissues of 24 organ donors aged 20-75 years. We applied a multimodal classifier to annotate the major immune cell lineages (T cells, B cells, innate lymphoid cells, and myeloid cells) and their corresponding subsets across the body, leveraging probabilistic modeling to define bases for immune variations across donors, tissue, and age. We identified dominant tissue-specific effects on immune cell composition and function across lineages for lymphoid sites, intestines, and blood-rich tissues. Age-associated effects were intrinsic to both lineage and site as manifested by macrophages in mucosal sites, B cells in lymphoid organs, and T and NK cells in blood-rich sites. Our results reveal tissue-specific signatures of immune homeostasis throughout the body and across different ages. This information provides a basis for defining the transcriptional underpinnings of immune variation and potential associations with disease-associated immune pathologies across the human lifespan.

Human cerebellar organoids with functional Purkinje cells

Research on human cerebellar development and disease has been hampered by the need for a human cell-based system that recapitulates the human cerebellum's cellular diversity and functional features. Here, we report a human organoid model (human cerebellar organoids [hCerOs]) capable of developing the complex cellular diversity of the fetal cerebellum, including a human-specific rhombic lip progenitor population that have never been generated in vitro prior to this study. 2-month-old hCerOs form distinct cytoarchitectural features, including laminar organized layering, and create functional connections between inhibitory and excitatory neurons that display coordinated network activity. Long-term culture of hCerOs allows healthy survival and maturation of Purkinje cells that display molecular and electrophysiological hallmarks of their in vivo counterparts, addressing a long-standing challenge in the field. This study therefore provides a physiologically relevant, all-human model system to elucidate the cell-type-specific mechanisms governing cerebellar development and disease.

Alzheimer's Disease: Models and Molecular Mechanisms Informing Disease and Treatments

Alzheimer's Disease (AD) is a complex neurodegenerative disease resulting in progressive loss of memory, language and motor abilities caused by cortical and hippocampal degeneration. This review captures the landscape of understanding of AD pathology, diagnostics, and current therapies. Two major mechanisms direct AD pathology: (1) accumulation of amyloid β (Aβ) plaque and (2) tau-derived neurofibrillary tangles (NFT). The most common variants in the Aβ pathway in APP, PSEN1, and PSEN2 are largely responsible for early-onset AD (EOAD), while MAPT, APOE, TREM2 and ABCA7 have a modifying effect on late-onset AD (LOAD). More recent studies implicate chaperone proteins and Aβ degrading proteins in AD. Several tests, such as cognitive function, brain imaging, and cerebral spinal fluid (CSF) and blood tests, are used for AD diagnosis. Additionally, several biomarkers seem to have a unique AD specific combination of expression and could potentially be used in improved, less invasive diagnostics. In addition to genetic perturbations, environmental influences, such as altered gut microbiome signatures, affect AD. Effective AD treatments have been challenging to develop. Currently, there are several FDA approved drugs (cholinesterase inhibitors, Aß-targeting antibodies and an NMDA antagonist) that could mitigate AD rate of decline and symptoms of distress.

Multiomic spatial landscape of innate immune cells at human central nervous system borders

The innate immune compartment of the human central nervous system (CNS) is highly diverse and includes several immune-cell populations such as macrophages that are frequent in the brain parenchyma (microglia) and less numerous at the brain interfaces as CNS-associated macrophages (CAMs). Due to their scantiness and particular location, little is known about the presence of temporally and spatially restricted CAM subclasses during development, health and perturbation. Here we combined single-cell RNA sequencing, time-of-flight mass cytometry and single-cell spatial transcriptomics with fate mapping and advanced immunohistochemistry to comprehensively characterize the immune system at human CNS interfaces with over 356,000 analyzed transcriptomes from 102 individuals. We also provide a comprehensive analysis of resident and engrafted myeloid cells in the brains of 15 individuals with peripheral blood stem cell transplantation, revealing compartment-specific engraftment rates across different CNS interfaces. Integrated multiomic and high-resolution spatial transcriptome analysis of anatomically dissected glioblastoma samples shows regionally distinct myeloid cell-type distributions driven by hypoxia. Notably, the glioblastoma-associated hypoxia response was distinct from the physiological hypoxia response in fetal microglia and CAMs. Our results highlight myeloid diversity at the interfaces of the human CNS with the periphery and provide insights into the complexities of the human brain's immune system.

Automatic cell-type harmonization and integration across Human Cell Atlas datasets

Harmonizing cell types across the single-cell community and assembling them into a common framework is central to building a standardized Human Cell Atlas. Here, we present CellHint, a predictive clustering tree-based tool to resolve cell-type differences in annotation resolution and technical biases across datasets. CellHint accurately quantifies cell-cell transcriptomic similarities and places cell types into a relationship graph that hierarchically defines shared and unique cell subtypes. Application to multiple immune datasets recapitulates expert-curated annotations. CellHint also reveals underexplored relationships between healthy and diseased lung cell states in eight diseases. Furthermore, we present a workflow for fast cross-dataset integration guided by harmonized cell types and cell hierarchy, which uncovers underappreciated cell types in adult human hippocampus. Finally, we apply CellHint to 12 tissues from 38 datasets, providing a deeply curated cross-tissue database with ∼3.7 million cells and various machine learning models for automatic cell annotation across human tissues.

The molecular cytoarchitecture of the adult mouse brain

The function of the mammalian brain relies upon the specification and spatial positioning of diversely specialized cell types. Yet, the molecular identities of the cell types and their positions within individual anatomical structures remain incompletely known. To construct a comprehensive atlas of cell types in each brain structure, we paired high-throughput single-nucleus RNA sequencing with Slide-seq1,2-a recently developed spatial transcriptomics method with near-cellular resolution-across the entire mouse brain. Integration of these datasets revealed the cell type composition of each neuroanatomical structure. Cell type diversity was found to be remarkably high in the midbrain, hindbrain and hypothalamus, with most clusters requiring a combination of at least three discrete gene expression markers to uniquely define them. Using these data, we developed a framework for genetically accessing each cell type, comprehensively characterized neuropeptide and neurotransmitter signalling, elucidated region-specific specializations in activity-regulated gene expression and ascertained the heritability enrichment of neurological and psychiatric phenotypes. These data, available as an online resource ( www.BrainCellData.org ), should find diverse applications across neuroscience, including the construction of new genetic tools and the prioritization of specific cell types and circuits in the study of brain diseases.

A high-resolution transcriptomic and spatial atlas of cell types in the whole mouse brain

The mammalian brain consists of millions to billions of cells that are organized into many cell types with specific spatial distribution patterns and structural and functional properties1-3. Here we report a comprehensive and high-resolution transcriptomic and spatial cell-type atlas for the whole adult mouse brain. The cell-type atlas was created by combining a single-cell RNA-sequencing (scRNA-seq) dataset of around 7 million cells profiled (approximately 4.0 million cells passing quality control), and a spatial transcriptomic dataset of approximately 4.3 million cells using multiplexed error-robust fluorescence in situ hybridization (MERFISH). The atlas is hierarchically organized into 4 nested levels of classification: 34 classes, 338 subclasses, 1,201 supertypes and 5,322 clusters. We present an online platform, Allen Brain Cell Atlas, to visualize the mouse whole-brain cell-type atlas along with the single-cell RNA-sequencing and MERFISH datasets. We systematically analysed the neuronal and non-neuronal cell types across the brain and identified a high degree of correspondence between transcriptomic identity and spatial specificity for each cell type. The results reveal unique features of cell-type organization in different brain regions-in particular, a dichotomy between the dorsal and ventral parts of the brain. The dorsal part contains relatively fewer yet highly divergent neuronal types, whereas the ventral part contains more numerous neuronal types that are more closely related to each other. Our study also uncovered extraordinary diversity and heterogeneity in neurotransmitter and neuropeptide expression and co-expression patterns in different cell types. Finally, we found that transcription factors are major determinants of cell-type classification and identified a combinatorial transcription factor code that defines cell types across all parts of the brain. The whole mouse brain transcriptomic and spatial cell-type atlas establishes a benchmark reference atlas and a foundational resource for integrative investigations of cellular and circuit function, development and evolution of the mammalian brain.

Cellular control of protein levels: A systems biology perspective

How cells regulate protein levels is a central question of biology. Over the past decades, molecular biology research has provided profound insights into the mechanisms and the molecular machinery governing each step of the gene expression process, from transcription to protein degradation. Recent advances in transcriptomics and proteomics have complemented our understanding of these fundamental cellular processes with a quantitative, systems-level perspective. Multi-omic studies revealed significant quantitative, kinetic and functional differences between the genome, transcriptome and proteome. While protein levels often correlate with mRNA levels, quantitative investigations have demonstrated a substantial impact of translation and protein degradation on protein expression control. In addition, protein-level regulation appears to play a crucial role in buffering protein abundances against undesirable mRNA expression variation. These findings have practical implications for many fields, including gene function prediction and precision medicine.

A spatially-resolved transcriptional atlas of the murine dorsal pons at single-cell resolution

The "dorsal pons", or "dorsal pontine tegmentum" (dPnTg), is part of the brainstem. It is a complex, densely packed region whose nuclei are involved in regulating many vital functions. Notable among them are the parabrachial nucleus, the Kölliker Fuse, the Barrington nucleus, the locus coeruleus, and the dorsal, laterodorsal, and ventral tegmental nuclei. In this study, we applied single-nucleus RNA-seq (snRNA-seq) to resolve neuronal subtypes based on their unique transcriptional profiles and then used multiplexed error robust fluorescence in situ hybridization (MERFISH) to map them spatially. We sampled ~1 million cells across the dPnTg and defined the spatial distribution of over 120 neuronal subtypes. Our analysis identified an unpredicted high transcriptional diversity in this region and pinpointed many neuronal subtypes' unique marker genes. We also demonstrated that many neuronal subtypes are transcriptionally similar between humans and mice, enhancing this study's translational value. Finally, we developed a freely accessible, GPU and CPU-powered dashboard (http://harvard.heavy.ai:6273/) that combines interactive visual analytics and hardware-accelerated SQL into a data science framework to allow the scientific community to query and gain insights into the data.

Single-cell genomics reveals region-specific developmental trajectories underlying neuronal diversity in the human hypothalamus

The development and diversity of neuronal subtypes in the human hypothalamus has been insufficiently characterized. To address this, we integrated transcriptomic data from 241,096 cells (126,840 newly generated) in the prenatal and adult human hypothalamus to reveal a temporal trajectory from proliferative stem cell populations to mature hypothalamic cell types. Iterative clustering of the adult neurons identified 108 robust transcriptionally distinct neuronal subtypes representing 10 hypothalamic nuclei. Pseudotime trajectories provided insights into the genes driving formation of these nuclei. Comparisons to single-cell transcriptomic data from the mouse hypothalamus suggested extensive conservation of neuronal subtypes despite certain differences in species-enriched gene expression. The uniqueness of hypothalamic neuronal lineages was examined developmentally by comparing excitatory lineages present in cortex and inhibitory lineages in ganglionic eminence, revealing both distinct and shared drivers of neuronal maturation across the human forebrain. These results provide a comprehensive transcriptomic view of human hypothalamus development through gestation and adulthood at cellular resolution.

Functional genomics and systems biology in human neuroscience

Neuroscience research has entered a phase of key discoveries in the realm of neurogenomics owing to strong financial and intellectual support for resource building and tool development. The previous challenge of tissue heterogeneity has been met with the application of techniques that can profile individual cells at scale. Moreover, the ability to perturb genes, gene regulatory elements and neuronal activity in a cell-type-specific manner has been integrated with gene expression studies to uncover the functional underpinnings of the genome at a systems level. Although these insights have necessarily been grounded in model systems, we now have the opportunity to apply these approaches in humans and in human tissue, thanks to advances in human genetics, brain imaging and tissue collection. We acknowledge that there will probably always be limits to the extent to which we can apply the genomic tools developed in model systems to human neuroscience; however, as we describe in this Perspective, the neuroscience field is now primed with an optimal foundation for tackling this ambitious challenge. The application of systems-level network analyses to these datasets will facilitate a deeper appreciation of human neurogenomics that cannot otherwise be achieved from directly observable phenomena.

A comparative atlas of single-cell chromatin accessibility in the human brain

Recent advances in single-cell transcriptomics have illuminated the diverse neuronal and glial cell types within the human brain. However, the regulatory programs governing cell identity and function remain unclear. Using a single-nucleus assay for transposase-accessible chromatin using sequencing (snATAC-seq), we explored open chromatin landscapes across 1.1 million cells in 42 brain regions from three adults. Integrating this data unveiled 107 distinct cell types and their specific utilization of 544,735 candidate cis-regulatory DNA elements (cCREs) in the human genome. Nearly a third of the cCREs demonstrated conservation and chromatin accessibility in the mouse brain cells. We reveal strong links between specific brain cell types and neuropsychiatric disorders including schizophrenia, bipolar disorder, Alzheimer's disease (AD), and major depression, and have developed deep learning models to predict the regulatory roles of noncoding risk variants in these disorders.

Single-cell DNA methylation and 3D genome architecture in the human brain

Delineating the gene-regulatory programs underlying complex cell types is fundamental for understanding brain function in health and disease. Here, we comprehensively examined human brain cell epigenomes by probing DNA methylation and chromatin conformation at single-cell resolution in 517 thousand cells (399 thousand neurons and 118 thousand non-neurons) from 46 regions of three adult male brains. We identified 188 cell types and characterized their molecular signatures. Integrative analyses revealed concordant changes in DNA methylation, chromatin accessibility, chromatin organization, and gene expression across cell types, cortical areas, and basal ganglia structures. We further developed single-cell methylation barcodes that reliably predict brain cell types using the methylation status of select genomic sites. This multimodal epigenomic brain cell atlas provides new insights into the complexity of cell-type-specific gene regulation in adult human brains.

Comprehensive cell atlas of the first-trimester developing human brain

The adult human brain comprises more than a thousand distinct neuronal and glial cell types, a diversity that emerges during early brain development. To reveal the precise sequence of events during early brain development, we used single-cell RNA sequencing and spatial transcriptomics and uncovered cell states and trajectories in human brains at 5 to 14 postconceptional weeks (pcw). We identified 12 major classes that are organized as ~600 distinct cell states, which map to precise spatial anatomical domains at 5 pcw. We described detailed differentiation trajectories of the human forebrain and midbrain and found a large number of region-specific glioblasts that mature into distinct pre-astrocytes and pre-oligodendrocyte precursor cells. Our findings reveal the establishment of cell types during the first trimester of human brain development.

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.

Signature morphoelectric properties of diverse GABAergic interneurons in the human neocortex

Human cortex transcriptomic studies have revealed a hierarchical organization of γ-aminobutyric acid-producing (GABAergic) neurons from subclasses to a high diversity of more granular types. Rapid GABAergic neuron viral genetic labeling plus Patch-seq (patch-clamp electrophysiology plus single-cell RNA sequencing) sampling in human brain slices was used to reliably target and analyze GABAergic neuron subclasses and individual transcriptomic types. This characterization elucidated transitions between PVALB and SST subclasses, revealed morphological heterogeneity within an abundant transcriptomic type, identified multiple spatially distinct types of the primate-specialized double bouquet cells (DBCs), and shed light on cellular differences between homologous mouse and human neocortical GABAergic neuron types. These results highlight the importance of multimodal phenotypic characterization for refinement of emerging transcriptomic cell type taxonomies and for understanding conserved and specialized cellular properties of human brain cell types.

A family portrait of human brain cells

A cell census provides information on the source of human brain specialization.