Cited by 1,250Single-cell and spatial transcriptomicsCell Image Analysis TechniquesGenomics and Chromatin Dynamics

The DOI matches its Crossref record; 189 of the 189 references it cites were found in a registry or an index; 97 of the 118 sources its claims cite were checked for support.

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

Open access
Yao, Zizhen, van Velthoven, Cindy T. J., Kunst, Michael, Zhang, Meng, McMillen, Delissa, Lee, Changkyu, Jung, Won, Goldy, Jeff, Abdelhak, Aliya, Aitken, Matthew, Baker, Katherine, Baker, Pamela, Barkan, Eliza, Bertagnolli, Darren, Bhandiwad, Ashwin, Bielstein, Cameron, Bishwakarma, Prajal, Campos, Jazmin, Carey, Daniel, Casper, Tamara, Chakka, Anish Bhaswanth, Chakrabarty, Rushil, Chavan, Sakshi, Chen, Min, Clark, Michael, Close, Jennie, Crichton, Kirsten, Daniel, Scott, DiValentin, Peter, Dolbeare, Tim, Ellingwood, Lauren, Fiabane, Elysha, Fliss, Timothy, Gee, James, Gerstenberger, James, Glandon, Alexandra, Gloe, Jessica, Gould, Joshua, Gray, James, Guilford, Nathan, Guzman, Junitta, Hirschstein, Daniel, Ho, Windy, Hooper, Marcus, Huang, Mike, Hupp, Madie, Jin, Kelly, Kroll, Matthew, Lathia, Kanan, Leon, Arielle, Li, Su, Long, Brian, Madigan, Zach, Malloy, Jessica, Malone, Jocelin, Maltzer, Zoe, Martin, Naomi, McCue, Rachel, McGinty, Ryan, Mei, Nicholas, Melchor, Jose, Meyerdierks, Emma, Mollenkopf, Tyler, Moonsman, Skyler, Nguyen, Thuc Nghi, Otto, Sven, Pham, Trangthanh, Rimorin, Christine, Ruiz, Augustin, Sanchez, Raymond, Sawyer, Lane, Shapovalova, Nadiya, Shepard, Noah, Slaughterbeck, Cliff, Sulc, Josef, Tieu, Michael, Torkelson, Amy, Tung, Herman, Valera Cuevas, Nasmil, Vance, Shane, Wadhwani, Katherine, Ward, Katelyn, Levi, Boaz, Farrell, Colin, Young, Rob, Staats, Brian, Wang, Ming-Qiang Michael, Thompson, Carol L., Mufti, Shoaib, Pagan, Chelsea M., Kruse, Lauren, Dee, Nick, Sunkin, Susan M., Esposito, Luke, Hawrylycz, Michael J., Waters, Jack, Ng, Lydia, Smith, Kimberly, Tasic, Bosiljka, Zhuang, Xiaowei, Zeng, Hongkui·2023·Nature

Published version of a preprint

Abstract

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 properties 1–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.

Overview

AI snapshot

Written from the paper’s full text.

TL;DR

This paper presents a whole adult mouse brain transcriptomic and spatial cell-type atlas combining ~7 million scRNA-seq cells (~4.0 million QC-passed) with ~4.3 million MERFISH cells. The atlas defines 34 classes, 338 subclasses, 1,201 supertypes and 5,322 clusters, and reports high correspondence between transcriptomic identity and spatial specificity.

Key findings

  • The final taxonomy contains 34 classes, 338 subclasses, 1,201 supertypes and 5,322 clusters; neuronal types account for 29 classes (85%), 315 subclasses (93%), 1,156 supertypes (96%) and 5,205 clusters (98%), while non-neuronal types form 5 classes, 23 subclasses, 45 supertypes and 117 clusters.
  • Neuronal subclasses are spatially restricted, with high inequality in spatial distribution by Gini coefficient and Shannon diversity; transcriptomically similar cell types are spatially closer, and mutual prediction between transcriptomic identity and CCFv3 location is high, with confusions mostly among closely related subclasses or neighbouring regions.
  • The brain shows a dorsal–ventral dichotomy: dorsal regions contain fewer yet highly divergent neuronal types, whereas ventral regions contain more numerous, more closely related types; hypothalamus, midbrain and hindbrain have the largest numbers of clusters and smaller spatial spans, while pallium, CNU, thalamus and cerebellum have fewer clusters after volume normalization.
  • The study identifies 62 glutamate–GABA dual-transmitter clusters, widespread modulatory neurotransmitter co-release, neuropeptide expression in about 79% of clusters, and 52 transcription factor co-expression modules; transcription factor-based subclass recall is 0.94 versus 0.98 with all DEGs, 0.90 with functional genes, 0.88 with adhesion genes and 0.75 with random DEGs.

Methods

  • Generated 781 scRNA-seq libraries from CCFv3-guided microdissections, yielding around 7.0 million single-cell transcriptomes; stringent QC retained around 4.0 million high-quality cells, with 10xv2 and 10xv3 datasets jointly clustered and supplemented by 1,687 10xMulti snRNA-seq nuclei across 33 clusters.
  • Used MERFISH with two gene panels (1,147 genes and 500 genes); the AIBS MERFISH dataset comprised 59 serial full coronal sections at 200-µm intervals, around 4.3 million segmented and QC-passed cells, registered to Allen CCFv3 and hierarchically mapped to the scRNA-seq taxonomy.
  • Main analyses included iterative clustering, 8,460 differentially expressed marker genes, fivefold cross-validation, UMAP and constellation plots, Gini coefficient and Shannon diversity, mutual prediction between transcriptome and spatial location, and WGCNA for transcription factor modules.

Limitations

  • Stated: dissociation of adult brain tissue produced many low-quality cells, damaged cells and mixed debris; stringent QC disqualified a high proportion of cells and likely selectively depleted large neurons such as Purkinje cells and large midbrain/hindbrain motor neurons, so midbrain and hindbrain cell types may not be fully represented or resolved.
  • Stated: the complex relationships among cell types cannot be fully captured by a one-dimensional hierarchy or two-dimensional UMAPs; not every cluster may represent a true cell type, and the categorization may need revision with better methods or more experimental evidence, especially developmental data.
  • Evident: spatial quantification used only 51/59 MERFISH sections within CCFv3 boundaries and 3,062,367 cells after filtering; 7 small subregions were completely missed, and the authors note that CCFv3 registration inaccuracies may contribute to confusion among neighbouring regions.

Model-written summary. Check it against the paper before relying on it.

Integrity signals

Registry

Crossref · 10.1038/s41586-023-06812-z

Cites

189 of 189 references verified0 citation errors · 0 not checkable

Statistics

Nothing to recompute.

Figures

Fig. 1 Transcriptomic cell-type taxonomy of the whole mouse brain. a, Left, the transcriptomic taxonomy tree of 338 subclasses organized in a dendrogram (10xv2: n = 1,699,939 cells; 10xv3: n = 2,341,350 cells; 10x Multiome: n = 1,687 nuclei). The neighbourhood and class levels are marked on the taxonomy tree. Classes marked with asterisks are included in the NN–IMN-GC neighbourhood. The IDs of every third subclass are shown to the right of the dendrogram. Full subclass names are provided in Supplementary Table 7. Following subclass IDs, bar plots represent (left to right): major neurotransmitter type, region distribution of profiled cells, number of clusters per subclass, number of RNA-seq cells analysed per subclass, and number of cells analysed by MERFISH per subclass. Subclasses marked with grey dots contain sex-dominant clusters. Sex-dominant clusters within a subclass are identified by calculating the odds and log P value for male/female distribution per cluster. Clusters with odds < 0.2 and log10(P value) < −10 are considered to be sex-dominant. b–e, UMAP representation of all cell types coloured by class (b), subclass (c), brain region (d) and major neurotransmitter type (e). Colour schemes for a–e are shown in the key at the bottom right of the figure. Astro, astrocyte; CB, cerebellum; CGE, caudal ganglionic eminence; CNU, cerebral nuclei; CR, Cajal–Retzius; CT, corticothalamic; CTX, cerebral cortex; CTXsp, cortical subplate; DG, dentate gyrus; EA, extended amygdala; Epen, ependymal; EPI, epithalamus; ET, extratelencephalic; GC, granule cell; HB, hindbrain; HPF, hippocampal formation; HY, hypothalamus; HYa, anterior hypothalamic; IMN, immature neurons; IT, intratelencephalic; L6b, layer 6b; LGE, lateral ganglionic eminence; LH, lateral habenula; LSX, lateral septal complex; MB, midbrain; MGE, medial ganglionic eminence; MH, medial habenula; MM, medial mammillary nucleus; MY, medulla; NN, non-neuronal; NP, near-projecting; OB, olfactory bulb; OEC, olfactory ensheathing cells; OLF, olfactory areas; Oligo, oligodendrocytes; OPC, oligodendrocyte precursor cells; P, pons; PAL, pallidum; STR, striatum; TH, thalamus. Neurotransmitter types: Chol, cholinergic; Dopa, dopaminergic; GABA, GABAergic; Glut, glutamatergic; Glyc, glycinergic; Hist, histaminergic; Nora, noradrenergic; Sero, serotonergic; NA, not applicable (no neurotransmitter detected).
Fig. 1 Transcriptomic cell-type taxonomy of the whole mouse brain. a, Left, the transcriptomic taxonomy tree of 338 subclasses organized in a dendrogram (10xv2: n = 1,699,939 cells; 10xv3: n = 2,341,350 cells; 10x Multiome: n = 1,687 nuclei). The neighbourhood and class levels are marked on the taxonomy tree. Classes marked with asterisks are included in the NN–IMN-GC neighbourhood. The IDs of every third subclass are shown to the right of the dendrogram. Full subclass names are provided in Supplementary Table 7. Following subclass IDs, bar plots represent (left to right): major neurotransmitter type, region distribution of profiled cells, number of clusters per subclass, number of RNA-seq cells analysed per subclass, and number of cells analysed by MERFISH per subclass. Subclasses marked with grey dots contain sex-dominant clusters. Sex-dominant clusters within a subclass are identified by calculating the odds and log P value for male/female distribution per cluster. Clusters with odds < 0.2 and log10(P value) < −10 are considered to be sex-dominant. b–e, UMAP representation of all cell types coloured by class (b), subclass (c), brain region (d) and major neurotransmitter type (e). Colour schemes for a–e are shown in the key at the bottom right of the figure. Astro, astrocyte; CB, cerebellum; CGE, caudal ganglionic eminence; CNU, cerebral nuclei; CR, Cajal–Retzius; CT, corticothalamic; CTX, cerebral cortex; CTXsp, cortical subplate; DG, dentate gyrus; EA, extended amygdala; Epen, ependymal; EPI, epithalamus; ET, extratelencephalic; GC, granule cell; HB, hindbrain; HPF, hippocampal formation; HY, hypothalamus; HYa, anterior hypothalamic; IMN, immature neurons; IT, intratelencephalic; L6b, layer 6b; LGE, lateral ganglionic eminence; LH, lateral habenula; LSX, lateral septal complex; MB, midbrain; MGE, medial ganglionic eminence; MH, medial habenula; MM, medial mammillary nucleus; MY, medulla; NN, non-neuronal; NP, near-projecting; OB, olfactory bulb; OEC, olfactory ensheathing cells; OLF, olfactory areas; Oligo, oligodendrocytes; OPC, oligodendrocyte precursor cells; P, pons; PAL, pallidum; STR, striatum; TH, thalamus. Neurotransmitter types: Chol, cholinergic; Dopa, dopaminergic; GABA, GABAergic; Glut, glutamatergic; Glyc, glycinergic; Hist, histaminergic; Nora, noradrenergic; Sero, serotonergic; NA, not applicable (no neurotransmitter detected).CC BY · PubMed Central
Fig. 2 Neuronal cell-type classification and distribution across the brain. a–l, UMAP representation (a,c,e,g,i,k) and representative MERFISH sections (b,d,f,h,j,l) of Pallium-Glut (a,b), Subpallium-GABA (c,d), HY–EA-Glut–GABA (e,f), TH–EPI-Glut (g,h), MB–HB-Glut–Sero–Dopa (i,j) and MB–HB–CB-GABA (k,l) neighbourhoods coloured by subclass. Each subclass is labelled with its ID and shown in the same colour in UMAP representations and MERFISH sections. a,c,e,g,i,k, Outlines in UMAP representations show cell classes. For full subclass names, see Supplementary Table 7.
Fig. 2 Neuronal cell-type classification and distribution across the brain. a–l, UMAP representation (a,c,e,g,i,k) and representative MERFISH sections (b,d,f,h,j,l) of Pallium-Glut (a,b), Subpallium-GABA (c,d), HY–EA-Glut–GABA (e,f), TH–EPI-Glut (g,h), MB–HB-Glut–Sero–Dopa (i,j) and MB–HB–CB-GABA (k,l) neighbourhoods coloured by subclass. Each subclass is labelled with its ID and shown in the same colour in UMAP representations and MERFISH sections. a,c,e,g,i,k, Outlines in UMAP representations show cell classes. For full subclass names, see Supplementary Table 7.CC BY · PubMed Central
Fig. 3 Modulatory neurotransmitter types and their distribution throughout the brain. a,b, Neuronal subclasses containing clusters that release modulatory neurotransmitters and their various co-release combinations with glutamate and/or GABA. UMAPs are coloured by subclass (a) and neurotransmitter type (b). c, Representative MERFISH sections showing the location of neuronal types expressing modulatory neurotransmitters. Cells are coloured by neurotransmitter type and labelled by subclass ID. See Supplementary Table 7 for detailed neurotransmitter assignment for each cluster. ADP, anterodorsal preoptic nucleus; AHN, anterior hypothalamic nucleus; ARH, arcuate hypothalamic nucleus; CLI, central linear nucleus raphe; CUN, cuneiform nucleus; DMH, dorsomedial nucleus of the hypothalamus; DMX, dorsal motor nucleus of the vagus nerve; IF, interfascicular nucleus raphe; LHA, lateral hypothalamic area; MDRN, medullary reticular nucleus; MPN, medial preoptic nucleus; MPO, medial preoptic area; MV, medial vestibular nucleus; NTS, nucleus of the solitary tract; PAG, periaqueductal grey; PARN, parvicellular reticular nucleus; PB, parabrachial nucleus; PBG, parabigeminal nucleus; PGRN, paragigantocellular reticular nucleus; PGRNd, paragigantocellular reticular nucleus, dorsal part; PH, posterior hypothalamic nucleus; PMv, ventral premammillary nucleus; PPN, pedunculopontine nucleus; PVa, periventricular hypothalamic nucleus, anterior part; PVHd, paraventricular hypothalamic nucleus, descending division; PVi, periventricular hypothalamic nucleus, intermediate part; PVpo, periventricular hypothalamic nucleus, preoptic part; PVR, periventricular region; RAmb, midbrain raphe nuclei; RL, rostral linear nucleus raphe; SBPV, subparaventricular zone; SNc, substantia nigra, compact part; SPIV, spinal vestibular nucleus; TMv, tuberomammillary nucleus, ventral part; VII, facial motor nucleus; VMPO, ventromedial preoptic nucleus; VTA, ventral tegmental area; ZI, zona incerta.
Fig. 3 Modulatory neurotransmitter types and their distribution throughout the brain. a,b, Neuronal subclasses containing clusters that release modulatory neurotransmitters and their various co-release combinations with glutamate and/or GABA. UMAPs are coloured by subclass (a) and neurotransmitter type (b). c, Representative MERFISH sections showing the location of neuronal types expressing modulatory neurotransmitters. Cells are coloured by neurotransmitter type and labelled by subclass ID. See Supplementary Table 7 for detailed neurotransmitter assignment for each cluster. ADP, anterodorsal preoptic nucleus; AHN, anterior hypothalamic nucleus; ARH, arcuate hypothalamic nucleus; CLI, central linear nucleus raphe; CUN, cuneiform nucleus; DMH, dorsomedial nucleus of the hypothalamus; DMX, dorsal motor nucleus of the vagus nerve; IF, interfascicular nucleus raphe; LHA, lateral hypothalamic area; MDRN, medullary reticular nucleus; MPN, medial preoptic nucleus; MPO, medial preoptic area; MV, medial vestibular nucleus; NTS, nucleus of the solitary tract; PAG, periaqueductal grey; PARN, parvicellular reticular nucleus; PB, parabrachial nucleus; PBG, parabigeminal nucleus; PGRN, paragigantocellular reticular nucleus; PGRNd, paragigantocellular reticular nucleus, dorsal part; PH, posterior hypothalamic nucleus; PMv, ventral premammillary nucleus; PPN, pedunculopontine nucleus; PVa, periventricular hypothalamic nucleus, anterior part; PVHd, paraventricular hypothalamic nucleus, descending division; PVi, periventricular hypothalamic nucleus, intermediate part; PVpo, periventricular hypothalamic nucleus, preoptic part; PVR, periventricular region; RAmb, midbrain raphe nuclei; RL, rostral linear nucleus raphe; SBPV, subparaventricular zone; SNc, substantia nigra, compact part; SPIV, spinal vestibular nucleus; TMv, tuberomammillary nucleus, ventral part; VII, facial motor nucleus; VMPO, ventromedial preoptic nucleus; VTA, ventral tegmental area; ZI, zona incerta.CC BY · PubMed Central
Fig. 4 Non-neuronal cell types and immature neuronal types. a, UMAP representation of the NN–IMN–GC neighbourhood coloured by subclass. Outlines show cell classes. b–d, Three subpopulations indicated in a are highlighted and further investigated: astrocytes (b), ependymal cells (c) and VLMCs (d). UMAP representation and representative MERFISH sections of astrocytes (b), ependymal cells (c) and VLMCs (d) are coloured and numbered by cluster. b,c, Outlines in UMAPs show subclasses. e, Colocalization of tanycyte, ependymal cell and VLMC clusters around V3 and ME, as shown in selected MERFISH sections. f, Colocalization of VLMC, CHOR and ependymal cell clusters in various ventricles, as shown in selected MERFISH sections. ABC, arachnoid barrier cells; BAM, border-associated macrophages; CHOR, choroid plexus; DC, dendritic cells; DCO, dorsal cochlear nucleus; Endo, endothelial cells; NT, non-telencephalon; Peri, pericytes; PIR, piriform cortex; SMC, smooth muscle cells; TE, telencephalon; UBC, unipolar brush cells; VLMC, vascular leptomeningeal cells.
Fig. 4 Non-neuronal cell types and immature neuronal types. a, UMAP representation of the NN–IMN–GC neighbourhood coloured by subclass. Outlines show cell classes. b–d, Three subpopulations indicated in a are highlighted and further investigated: astrocytes (b), ependymal cells (c) and VLMCs (d). UMAP representation and representative MERFISH sections of astrocytes (b), ependymal cells (c) and VLMCs (d) are coloured and numbered by cluster. b,c, Outlines in UMAPs show subclasses. e, Colocalization of tanycyte, ependymal cell and VLMC clusters around V3 and ME, as shown in selected MERFISH sections. f, Colocalization of VLMC, CHOR and ependymal cell clusters in various ventricles, as shown in selected MERFISH sections. ABC, arachnoid barrier cells; BAM, border-associated macrophages; CHOR, choroid plexus; DC, dendritic cells; DCO, dorsal cochlear nucleus; Endo, endothelial cells; NT, non-telencephalon; Peri, pericytes; PIR, piriform cortex; SMC, smooth muscle cells; TE, telencephalon; UBC, unipolar brush cells; VLMC, vascular leptomeningeal cells.CC BY · PubMed Central
Fig. 5 Transcription factor modules across the whole mouse brain. a, Distribution of the number of differentially expressed transcription factors (TFs) between neuronal and non-neuronal classes, between classes, between subclasses, and within subclasses. b, Cross-validation accuracy for each cluster (left) or subclass (right) using classifiers built based on all 8,460 marker genes (all), 534 transcription factor marker genes (TF), 541 functional marker genes, 857 marker genes encoding adhesion molecules (adhesion), 534 randomly selected adhesion marker genes (random adhesion), or 534 randomly selected marker genes (random). c, Density plot showing distribution of correlation of marker gene expression between clusters using all markers, adhesion marker genes, functional genes and transcription factors. Correlation values are derived from full correlation matrices shown in Extended Data Fig. 3. d, Expression of key transcription factors for each subclass in the taxonomy tree, organized in transcription factor co-expression modules shown as colour bars on both sides of the heat map. Module IDs are shown on the left, exemplary transcription factor genes are shown on the right. For a full list of transcription factor genes in each module (in the same order as in this heat map), see Supplementary Table 8. Avg, mean. Source Data
Fig. 5 Transcription factor modules across the whole mouse brain. a, Distribution of the number of differentially expressed transcription factors (TFs) between neuronal and non-neuronal classes, between classes, between subclasses, and within subclasses. b, Cross-validation accuracy for each cluster (left) or subclass (right) using classifiers built based on all 8,460 marker genes (all), 534 transcription factor marker genes (TF), 541 functional marker genes, 857 marker genes encoding adhesion molecules (adhesion), 534 randomly selected adhesion marker genes (random adhesion), or 534 randomly selected marker genes (random). c, Density plot showing distribution of correlation of marker gene expression between clusters using all markers, adhesion marker genes, functional genes and transcription factors. Correlation values are derived from full correlation matrices shown in Extended Data Fig. 3. d, Expression of key transcription factors for each subclass in the taxonomy tree, organized in transcription factor co-expression modules shown as colour bars on both sides of the heat map. Module IDs are shown on the left, exemplary transcription factor genes are shown on the right. For a full list of transcription factor genes in each module (in the same order as in this heat map), see Supplementary Table 8. Avg, mean. Source DataCC BY · PubMed Central
Fig. 6 Brain region-specific features of cell types. a, Heat map showing the CCFv3 region distribution (y axis) in each subclass (x axis) for MERFISH cells. Bar graphs on the left show the broad CCFv3 regions, proportion of neuronal versus glial cells per region of interest (ROI), and proportion of neurotransmitter types per ROI. Bar graphs on the right show broad CCFv3 regions, Shannon diversity per subclass and supertype, and number of cells per ROI. Bar graphs on the top show number of cells per subclass, Gini coefficient and class assignment. Bar graphs on the bottom show subclass and class annotations. b, Scatter plot showing the number of neuronal clusters identified per major brain region versus the number of neuronal cells profiled by scRNA-seq in the corresponding region. Each neuronal cluster is assigned to the most dominant region. c, As in b, except numbers of clusters and profiled cells are normalized by the region volume. d, Distribution of the number of DEGs (identified in scRNA-seq data) between every pair of neuronal clusters within each major brain region, split into indicated quantiles. The curves show the spread of the number of DEGs between more similar types at 0.1 quantile versus the more distinct types at 0.9 quantile. e, Scatter plot showing the number of cells mapped to a given neuronal cluster versus the span (as measured by IQR) of their 3D coordinates along the anteroposterior, dorsoventral and mediolateral axes based on the MERFISH dataset, stratified by the major brain regions. Note that both axes are in log scales. The plot shows how localized the clusters are within each region along each spatial axis. IQR, inter-quantile range (the difference between 75% quantile and 25% quantile). Pall, pallium. Source Data
Fig. 6 Brain region-specific features of cell types. a, Heat map showing the CCFv3 region distribution (y axis) in each subclass (x axis) for MERFISH cells. Bar graphs on the left show the broad CCFv3 regions, proportion of neuronal versus glial cells per region of interest (ROI), and proportion of neurotransmitter types per ROI. Bar graphs on the right show broad CCFv3 regions, Shannon diversity per subclass and supertype, and number of cells per ROI. Bar graphs on the top show number of cells per subclass, Gini coefficient and class assignment. Bar graphs on the bottom show subclass and class annotations. b, Scatter plot showing the number of neuronal clusters identified per major brain region versus the number of neuronal cells profiled by scRNA-seq in the corresponding region. Each neuronal cluster is assigned to the most dominant region. c, As in b, except numbers of clusters and profiled cells are normalized by the region volume. d, Distribution of the number of DEGs (identified in scRNA-seq data) between every pair of neuronal clusters within each major brain region, split into indicated quantiles. The curves show the spread of the number of DEGs between more similar types at 0.1 quantile versus the more distinct types at 0.9 quantile. e, Scatter plot showing the number of cells mapped to a given neuronal cluster versus the span (as measured by IQR) of their 3D coordinates along the anteroposterior, dorsoventral and mediolateral axes based on the MERFISH dataset, stratified by the major brain regions. Note that both axes are in log scales. The plot shows how localized the clusters are within each region along each spatial axis. IQR, inter-quantile range (the difference between 75% quantile and 25% quantile). Pall, pallium. Source DataCC BY · PubMed Central
Extended Data Fig. 1 scRNA-seq data analysis workflow. (a) Number of cells at each step in the scRNA-seq data analysis pipeline. The identification of doublets and low-quality clusters is described in more detail in Methods. The 10xv2 and 10xv3 data were first QC-ed and analyzed separately. After initial clustering the datasets were combined and QC-ed again before and after joint clustering. 10x Multiome snRNA-seq data was added to fill in gaps that were identified after joint clustering of 10xv2 and 10xv3 scRNA-seq data. (b-c) Gene count and qc score thresholds used for each of the four major cell populations (neuroglial cells, neurons, immature neurons and granule cells, and other) on the 10xv2 (b) and 10xv3 (c) datasets. (d-e) Number of cells isolated from dissection ROIs (pre-QC) and number of cells passing QC (post-QC) for 10xv2 (d) and 10xv3 (e) datasets. We didn’t profile LSX, STR, sAMY, PAL, Pons, MY, and CB by 10xv2. Some regions were collected using different dissections between 10xv2 and 10xv3, but all regions were covered by 10xv3.
Extended Data Fig. 1 scRNA-seq data analysis workflow. (a) Number of cells at each step in the scRNA-seq data analysis pipeline. The identification of doublets and low-quality clusters is described in more detail in Methods. The 10xv2 and 10xv3 data were first QC-ed and analyzed separately. After initial clustering the datasets were combined and QC-ed again before and after joint clustering. 10x Multiome snRNA-seq data was added to fill in gaps that were identified after joint clustering of 10xv2 and 10xv3 scRNA-seq data. (b-c) Gene count and qc score thresholds used for each of the four major cell populations (neuroglial cells, neurons, immature neurons and granule cells, and other) on the 10xv2 (b) and 10xv3 (c) datasets. (d-e) Number of cells isolated from dissection ROIs (pre-QC) and number of cells passing QC (post-QC) for 10xv2 (d) and 10xv3 (e) datasets. We didn’t profile LSX, STR, sAMY, PAL, Pons, MY, and CB by 10xv2. Some regions were collected using different dissections between 10xv2 and 10xv3, but all regions were covered by 10xv3.CC BY · PubMed Central
Extended Data Fig. 2 MERFISH data generation, data processing and summary of results. (a) Workflow for generating and processing MERFISH data. (b) Correlation of gene detection between MERFISH and bulk RNA-sequencing for four different brain regions. (c) Histogram displaying the distribution of gene detection correlation between adjacent MERFISH sections. (d-f) Violin plots displaying distribution of cell volumes (d), gene detection (e), and mRNA molecule detection (f) for individual sections ordered from anterior to posterior (left panel) or cumulative distribution for the whole brain (right panel). Red dashed lines indicate cutoff for filtering. (g) Cumulative histogram showing the relative contribution of each subclass to each section ordered from anterior to posterior. (h) Heatmap showing the proportion of cells per region for each subclass in the MERFISH data (left) and scRNA-seq data (right). The heatmap in the middle shows the correlation between region distribution of MERFISH and scRNA-seq data.
Extended Data Fig. 2 MERFISH data generation, data processing and summary of results. (a) Workflow for generating and processing MERFISH data. (b) Correlation of gene detection between MERFISH and bulk RNA-sequencing for four different brain regions. (c) Histogram displaying the distribution of gene detection correlation between adjacent MERFISH sections. (d-f) Violin plots displaying distribution of cell volumes (d), gene detection (e), and mRNA molecule detection (f) for individual sections ordered from anterior to posterior (left panel) or cumulative distribution for the whole brain (right panel). Red dashed lines indicate cutoff for filtering. (g) Cumulative histogram showing the relative contribution of each subclass to each section ordered from anterior to posterior. (h) Heatmap showing the proportion of cells per region for each subclass in the MERFISH data (left) and scRNA-seq data (right). The heatmap in the middle shows the correlation between region distribution of MERFISH and scRNA-seq data.CC BY · PubMed Central
Extended Data Fig. 3 Marker gene expression correlation matrices showing relatedness among cell types. (a-d) Heatmaps showing pairwise Pearson correlation of gene expression levels for each pair of clusters using marker gene sets of 534 transcription factors (a), 541 functional genes (including neuropeptides, GPCRs, ion channels, transporters, etc.) (b), 857 adhesion molecules (c), and all 8,460 marker genes (d). Correlations were computed using 10xv3 scRNA-seq data only except for the 31 nuclei-dominated clusters where 10xMulti snRNA-seq data were used. (e-g) All marker gene expression correlation between clusters compared to correlation between clusters of expression of functional marker genes (e), adhesion marker genes (f), and transcription factor (TF) marker genes (g). Correlation values are derived from a-d.
Extended Data Fig. 3 Marker gene expression correlation matrices showing relatedness among cell types. (a-d) Heatmaps showing pairwise Pearson correlation of gene expression levels for each pair of clusters using marker gene sets of 534 transcription factors (a), 541 functional genes (including neuropeptides, GPCRs, ion channels, transporters, etc.) (b), 857 adhesion molecules (c), and all 8,460 marker genes (d). Correlations were computed using 10xv3 scRNA-seq data only except for the 31 nuclei-dominated clusters where 10xMulti snRNA-seq data were used. (e-g) All marker gene expression correlation between clusters compared to correlation between clusters of expression of functional marker genes (e), adhesion marker genes (f), and transcription factor (TF) marker genes (g). Correlation values are derived from a-d.CC BY · PubMed Central
Extended Data Fig. 4 Comparison of the initial automated computational assignment of classes and subclasses with the manually revised, final assignment of classes and subclasses. Size of the dot corresponds to the number of overlapping cells (frequency) in corresponding classes or subclasses, and color represents the Jaccard similarity between corresponding classes or subclasses.
Extended Data Fig. 4 Comparison of the initial automated computational assignment of classes and subclasses with the manually revised, final assignment of classes and subclasses. Size of the dot corresponds to the number of overlapping cells (frequency) in corresponding classes or subclasses, and color represents the Jaccard similarity between corresponding classes or subclasses.CC BY · PubMed Central
Extended Data Fig. 5 Transcriptomic cell type taxonomy of the whole mouse brain with additional metadata information. (a-b) Number of genes (a) or number of UMIs (b) detected per cell in 10xv2 (top) or 10xv3 (bottom) datasets for each cell type neighborhood. The data shown is post-QC. Numbers of cells for 10xv2: Pallium-Glut n = 1,128,664, Subpallium-GABA n = 269,307, HY-EA-Glut-GABA n = 107,706, TH-EPI-Glut n = 73,702, MB-HB-CB-GABA n = 18,590, MB-HB-Glut-Sero-Dopa n = 20,089, IMN-GC n = 123,650, Neuroglial n = 80,959, Vascular n = 6,894, Immune n = 4,941. Numbers of cells for 10xv3: Pallium-Glut n = 366,137, Subpallium-GABA n = 342,116, HY-EA-Glut-GABA n = 187,742, TH-EPI-Glut n = 52,469, MB-HB-CB-GABA n = 167,425, MB-HB-Glut-Sero-Dopa n = 159,653, IMN-GC n = 209,310, Neuroglial n = 774,537, Vascular n = 130,599, Immune n = 87,639. (c-d) Number of genes (c) or number of UMIs (d) detected per nucleus in the 10xMulti dataset (post-QC) for each subclass where 10xMulti nuclei were added. Numbers of nuclei: 157 RN Spp1 Glut n = 48, 233 NLL-SOC Spp1 Glut n = 115, 254 VCO Mafa Meis2 Glut n = 490, 261 HB Calcb Chol n = 274, 278 NLL Gata3 Gly-Gaba n = 339, 279 PSV Pax2 Gly-Gaba n = 43, 280 NLL-po Pax7 Gaba n = 17, 297 CU-ECU Pax2 Gly-Gaba n = 15, 313 CBX Purkinje Gaba n = 346. All box plots include the median line, the box denotes the interquartile range (IQR), whiskers denote the rest of the data distribution, and outliers are denoted by points greater than ± 1.5× IQR. (e) The transcriptomic taxonomy tree of 338 subclasses organized in a dendrogram (same as Fig. 1a). From left to right, the bar plots represent neurotransmitter (NT) type assignment, heat map showing expression of major neurotransmitter marker genes, sex distribution, platform distribution, light-dark distribution of profiled cells, and number of donors that contributed to each subclass.
Extended Data Fig. 5 Transcriptomic cell type taxonomy of the whole mouse brain with additional metadata information. (a-b) Number of genes (a) or number of UMIs (b) detected per cell in 10xv2 (top) or 10xv3 (bottom) datasets for each cell type neighborhood. The data shown is post-QC. Numbers of cells for 10xv2: Pallium-Glut n = 1,128,664, Subpallium-GABA n = 269,307, HY-EA-Glut-GABA n = 107,706, TH-EPI-Glut n = 73,702, MB-HB-CB-GABA n = 18,590, MB-HB-Glut-Sero-Dopa n = 20,089, IMN-GC n = 123,650, Neuroglial n = 80,959, Vascular n = 6,894, Immune n = 4,941. Numbers of cells for 10xv3: Pallium-Glut n = 366,137, Subpallium-GABA n = 342,116, HY-EA-Glut-GABA n = 187,742, TH-EPI-Glut n = 52,469, MB-HB-CB-GABA n = 167,425, MB-HB-Glut-Sero-Dopa n = 159,653, IMN-GC n = 209,310, Neuroglial n = 774,537, Vascular n = 130,599, Immune n = 87,639. (c-d) Number of genes (c) or number of UMIs (d) detected per nucleus in the 10xMulti dataset (post-QC) for each subclass where 10xMulti nuclei were added. Numbers of nuclei: 157 RN Spp1 Glut n = 48, 233 NLL-SOC Spp1 Glut n = 115, 254 VCO Mafa Meis2 Glut n = 490, 261 HB Calcb Chol n = 274, 278 NLL Gata3 Gly-Gaba n = 339, 279 PSV Pax2 Gly-Gaba n = 43, 280 NLL-po Pax7 Gaba n = 17, 297 CU-ECU Pax2 Gly-Gaba n = 15, 313 CBX Purkinje Gaba n = 346. All box plots include the median line, the box denotes the interquartile range (IQR), whiskers denote the rest of the data distribution, and outliers are denoted by points greater than ± 1.5× IQR. (e) The transcriptomic taxonomy tree of 338 subclasses organized in a dendrogram (same as Fig. 1a). From left to right, the bar plots represent neurotransmitter (NT) type assignment, heat map showing expression of major neurotransmitter marker genes, sex distribution, platform distribution, light-dark distribution of profiled cells, and number of donors that contributed to each subclass.CC BY · PubMed Central
Extended Data Fig. 6 Constellation plot of the global relatedness between subclasses. Each subclass is represented by a disk, labeled by the subclass ID and positioned at the subclass centroid in UMAP coordinates shown in Fig. 1c. The size of the disk corresponds to the number of cells within each subclass, and the edge weights correspond to the fraction of shared neighbors (see Methods) between subclasses. Each subclass is colored by the class it belongs to. Curved outlines drawn around subclasses show the major neighborhoods. Source Data
Extended Data Fig. 6 Constellation plot of the global relatedness between subclasses. Each subclass is represented by a disk, labeled by the subclass ID and positioned at the subclass centroid in UMAP coordinates shown in Fig. 1c. The size of the disk corresponds to the number of cells within each subclass, and the edge weights correspond to the fraction of shared neighbors (see Methods) between subclasses. Each subclass is colored by the class it belongs to. Curved outlines drawn around subclasses show the major neighborhoods. Source DataCC BY · PubMed Central
Extended Data Fig. 7 Validation of data integration across 10xv2, 10xv3, and MERFISH datasets. (a-c) UMAP representation of all cell types colored by profiling platform (a), region (b), and subclass (c). Other than the regions only profiled by 10xv3 (LSX, STR, sAMY, PAL, Pons, MY), the cells from both 10xv2 and 10xv3 platforms integrate very well. Cell types in isocortex and HPF have a lot more 10xv2 cells, consistent with our sampling plan. For cell types/clusters containing many cells, we observed separation of 10xv2 and 10xv3 data in the UMAP space, but not at the cluster level. (d) Correlation of gene expression between 10xv2 and 10xv3 and between 10xv3 and MERFISH. For each gene, we computed the Pearson correlation of its average expression in each cluster across clusters between 10xv2 and 10xv3, and the correlation between 10xv3 and MERFISH. For 10xv3 and MERFISH comparison, distribution of the correlation values of all 500 genes in the MERFISH panel is shown. For 10xv3 and 10xv2 comparison, we show the correlation of 5383 marker genes based on 10xv2, and 466 10xv2 marker genes that are also present on the MERFISH gene panel (the other 34 MERFISH genes not shown have low expression in 10xv2 clusters). We manually inspected several genes with poor correlation and found them to have poor gene annotation or show relatively small variations across clusters. Most genes with low correlations between 10xv3 and MERFISH data are *Rik genes that are more likely to be poorly annotated, and the MERFISH probes selected for them might not work well. (e) 2D density plot showing on the X-axis the number of DEGs (based on 10xv3 dataset) present on the MERFISH gene panel between all pairs of clusters, and on the Y-axis the number of such DEGs showing the same direction of changes between corresponding pairs of mapped MERFISH clusters. Almost all the DEGs between all pairs of clusters show the same direction of changes between 10xv3 and MERFISH. (f) 2D density plot showing on the X-axis the number of DEGs (based on 10xv3 dataset) present on the MERFISH gene panel between all pairs of clusters, and on the Y-axis the number of such DEGs showing the same direction of changes, and |log2(FC)| > 1 between corresponding pairs of mapped MERFISH clusters. About 60% of DEGs between all pairs of clusters based on 10xv3 show significant fold change (FC) in MERFISH. (g) Similar analysis as in (f) but shown as violin plot by binning the number of 10xv3 DEGs present on the MERFISH gene panel on the X-axis, with better resolution on closely related pairs with four or fewer DEGs present on MERFISH gene panels. The MERFISH dataset can resolve the vast majority of clusters due to strong correlation of DEG expression between 10xv3 and MERFISH clusters. On the other hand, a few hundred pairs of clusters with fewer than two DEGs on the MERFISH gene panel remain unresolvable in the MERFISH data, and they are usually sibling clusters with indistinguishable spatial distribution. Source Data
Extended Data Fig. 7 Validation of data integration across 10xv2, 10xv3, and MERFISH datasets. (a-c) UMAP representation of all cell types colored by profiling platform (a), region (b), and subclass (c). Other than the regions only profiled by 10xv3 (LSX, STR, sAMY, PAL, Pons, MY), the cells from both 10xv2 and 10xv3 platforms integrate very well. Cell types in isocortex and HPF have a lot more 10xv2 cells, consistent with our sampling plan. For cell types/clusters containing many cells, we observed separation of 10xv2 and 10xv3 data in the UMAP space, but not at the cluster level. (d) Correlation of gene expression between 10xv2 and 10xv3 and between 10xv3 and MERFISH. For each gene, we computed the Pearson correlation of its average expression in each cluster across clusters between 10xv2 and 10xv3, and the correlation between 10xv3 and MERFISH. For 10xv3 and MERFISH comparison, distribution of the correlation values of all 500 genes in the MERFISH panel is shown. For 10xv3 and 10xv2 comparison, we show the correlation of 5383 marker genes based on 10xv2, and 466 10xv2 marker genes that are also present on the MERFISH gene panel (the other 34 MERFISH genes not shown have low expression in 10xv2 clusters). We manually inspected several genes with poor correlation and found them to have poor gene annotation or show relatively small variations across clusters. Most genes with low correlations between 10xv3 and MERFISH data are *Rik genes that are more likely to be poorly annotated, and the MERFISH probes selected for them might not work well. (e) 2D density plot showing on the X-axis the number of DEGs (based on 10xv3 dataset) present on the MERFISH gene panel between all pairs of clusters, and on the Y-axis the number of such DEGs showing the same direction of changes between corresponding pairs of mapped MERFISH clusters. Almost all the DEGs between all pairs of clusters show the same direction of changes between 10xv3 and MERFISH. (f) 2D density plot showing on the X-axis the number of DEGs (based on 10xv3 dataset) present on the MERFISH gene panel between all pairs of clusters, and on the Y-axis the number of such DEGs showing the same direction of changes, and |log2(FC)| > 1 between corresponding pairs of mapped MERFISH clusters. About 60% of DEGs between all pairs of clusters based on 10xv3 show significant fold change (FC) in MERFISH. (g) Similar analysis as in (f) but shown as violin plot by binning the number of 10xv3 DEGs present on the MERFISH gene panel on the X-axis, with better resolution on closely related pairs with four or fewer DEGs present on MERFISH gene panels. The MERFISH dataset can resolve the vast majority of clusters due to strong correlation of DEG expression between 10xv3 and MERFISH clusters. On the other hand, a few hundred pairs of clusters with fewer than two DEGs on the MERFISH gene panel remain unresolvable in the MERFISH data, and they are usually sibling clusters with indistinguishable spatial distribution. Source DataCC BY · PubMed Central
Extended Data Fig. 8 Validation of gene expression patterns of scRNA-seq transcriptomes imputed into MERFISH space. (a) Correlation of expression for all 500 genes in the MERFISH panel between MERFISH and 10xv3 (red), imputed MERFISH and MERFISH (green), and imputed MERFISH and 10xv3 (blue). To test the accuracy of MERFISH imputation, one gene is excluded from the gene panel at a time from KNN computation at all levels and its imputed gene expression is compared with its original gene expression. The distribution of correlations between imputed expression and the original MERFISH expression or the reference 10xv3 expression is shown for each gene at the cluster level. (b) Scatterplots showing the correlation between imputed MERFISH gene expression vs. MERFISH gene expression (left panels) and imputed MERFISH gene expression vs. 10xv3 gene expression (right panels) for selected genes, Calb2 (top row), Baiap3 (middle row), and Lypd1 (bottom row). (c-e) Examples of spatial gene expression patterns from MERFISH data (left panels), imputed MERFISH data (middle panels), and images from the Allen in situ hybridization (ISH) atlas (right panels) for select genes, Calb2 (c), Baiap3 (d), and Lypd1 (e). (f) Representative MERFISH sections showing imputed expression of Foxp2 (which was not directly profiled by MERFISH) and Allen ISH images. ISH image credit: Allen Institute for Brain Science, https://mouse.brain-map.org/. Source Data
Extended Data Fig. 8 Validation of gene expression patterns of scRNA-seq transcriptomes imputed into MERFISH space. (a) Correlation of expression for all 500 genes in the MERFISH panel between MERFISH and 10xv3 (red), imputed MERFISH and MERFISH (green), and imputed MERFISH and 10xv3 (blue). To test the accuracy of MERFISH imputation, one gene is excluded from the gene panel at a time from KNN computation at all levels and its imputed gene expression is compared with its original gene expression. The distribution of correlations between imputed expression and the original MERFISH expression or the reference 10xv3 expression is shown for each gene at the cluster level. (b) Scatterplots showing the correlation between imputed MERFISH gene expression vs. MERFISH gene expression (left panels) and imputed MERFISH gene expression vs. 10xv3 gene expression (right panels) for selected genes, Calb2 (top row), Baiap3 (middle row), and Lypd1 (bottom row). (c-e) Examples of spatial gene expression patterns from MERFISH data (left panels), imputed MERFISH data (middle panels), and images from the Allen in situ hybridization (ISH) atlas (right panels) for select genes, Calb2 (c), Baiap3 (d), and Lypd1 (e). (f) Representative MERFISH sections showing imputed expression of Foxp2 (which was not directly profiled by MERFISH) and Allen ISH images. ISH image credit: Allen Institute for Brain Science, https://mouse.brain-map.org/. Source DataCC BY · PubMed Central
Extended Data Fig. 9 Distribution of glutamate-GABA dual transmitting cell types throughout the brain. (a-b) Neuronal subclasses containing clusters releasing glutamate-GABA dual transmitters. UMAPs are colored by subclass (a) and neurotransmitter type (b). Glutamate-GABA co-releasing clusters are labeled by cluster ID in panel (b). (c) MERFISH sections showing glutamate-GABA co-releasing clusters colored by the subclass to which they belong. See Supplementary Table 7 for detailed neurotransmitter assignment for each cluster.
Extended Data Fig. 9 Distribution of glutamate-GABA dual transmitting cell types throughout the brain. (a-b) Neuronal subclasses containing clusters releasing glutamate-GABA dual transmitters. UMAPs are colored by subclass (a) and neurotransmitter type (b). Glutamate-GABA co-releasing clusters are labeled by cluster ID in panel (b). (c) MERFISH sections showing glutamate-GABA co-releasing clusters colored by the subclass to which they belong. See Supplementary Table 7 for detailed neurotransmitter assignment for each cluster.CC BY · PubMed Central
Extended Data Fig. 10 Neuropeptide distribution across the whole mouse brain. (a) Scatter plot of Tau score over the number of clusters each neuropeptide is expressed in at the level of log2(CPM) > 3. The Tau score is a measurement of cell type specificity, which varies from 0 to 1 where 0 means uniformly expressed and 1 means highly specific to one type. (b) Scatter plot of Tau score over the number of clusters each peptide-liganded G-protein coupled receptor (GPCR) gene is expressed in at the level of log2(CPM) > 3. (c) Expression level of neuropeptide in log2(CPM) per cluster. For each neuropeptide along the Y axis, clusters are sorted from the highest to lowest mean gene expression level along the X axis. (d) Expression level of neuropeptide in log2(CPM) per cluster. For each neuropeptide along the Y axis, clusters are sorted from the highest to lowest mean gene expression level along the X axis. For each gene, only the top 200 highest-expressing clusters out of 5,322 clusters are shown. (e) Representative MERFISH sections highlighting the spatial location of clusters expressing each of the 20 highly cell-type-specific neuropeptide genes (expressed in 8 or fewer clusters). (f) Representative MERFISH sections showing the expression of the neuropeptides present on the MERFISH gene panel that are widely expressed. We also note that the relationships between mRNA levels, the post-translationally processed peptide levels, and the functional levels are unknown for most neuropeptides, thus, it is difficult to predict what mRNA levels would lead to sufficient functional expression of a given neuropeptide.
Extended Data Fig. 10 Neuropeptide distribution across the whole mouse brain. (a) Scatter plot of Tau score over the number of clusters each neuropeptide is expressed in at the level of log2(CPM) > 3. The Tau score is a measurement of cell type specificity, which varies from 0 to 1 where 0 means uniformly expressed and 1 means highly specific to one type. (b) Scatter plot of Tau score over the number of clusters each peptide-liganded G-protein coupled receptor (GPCR) gene is expressed in at the level of log2(CPM) > 3. (c) Expression level of neuropeptide in log2(CPM) per cluster. For each neuropeptide along the Y axis, clusters are sorted from the highest to lowest mean gene expression level along the X axis. (d) Expression level of neuropeptide in log2(CPM) per cluster. For each neuropeptide along the Y axis, clusters are sorted from the highest to lowest mean gene expression level along the X axis. For each gene, only the top 200 highest-expressing clusters out of 5,322 clusters are shown. (e) Representative MERFISH sections highlighting the spatial location of clusters expressing each of the 20 highly cell-type-specific neuropeptide genes (expressed in 8 or fewer clusters). (f) Representative MERFISH sections showing the expression of the neuropeptides present on the MERFISH gene panel that are widely expressed. We also note that the relationships between mRNA levels, the post-translationally processed peptide levels, and the functional levels are unknown for most neuropeptides, thus, it is difficult to predict what mRNA levels would lead to sufficient functional expression of a given neuropeptide.CC BY · PubMed Central
Extended Data Fig. 11 Additional non-neuronal UMAPs and marker genes. (a-c) UMAP representation of the NN-IMN-GC neighborhood colored by subclass (a), region (b), and supertype (c). (d) Dot plot showing marker gene expression in non-neuronal subclasses. Dot size and color indicate proportion of expressing cells and average expression level in each subclass, respectively. (e) Dot plot showing marker gene expression in all clusters in the Astro-Epen class. Dot size and color indicate proportion of expressing cells and average expression level in each cluster, respectively. (f) Dot plot showing marker gene expression in VLMC clusters. Dot size and color indicate proportion of expressing cells and average expression level in each cluster, respectively. (g) Representative MERFISH sections showing the spatial gradient of OEC clusters. (h) Co-localization of VLMCs with interlaminar astrocytes (ILA) as shown in selected MERFISH sections. (i) UMAP representation of OPCs and oligodendrocytes colored and labeled by supertype. (j-k) Representative MERFISH sections showing the spatial distribution of OPC (j) and Oligo (k) supertypes.
Extended Data Fig. 11 Additional non-neuronal UMAPs and marker genes. (a-c) UMAP representation of the NN-IMN-GC neighborhood colored by subclass (a), region (b), and supertype (c). (d) Dot plot showing marker gene expression in non-neuronal subclasses. Dot size and color indicate proportion of expressing cells and average expression level in each subclass, respectively. (e) Dot plot showing marker gene expression in all clusters in the Astro-Epen class. Dot size and color indicate proportion of expressing cells and average expression level in each cluster, respectively. (f) Dot plot showing marker gene expression in VLMC clusters. Dot size and color indicate proportion of expressing cells and average expression level in each cluster, respectively. (g) Representative MERFISH sections showing the spatial gradient of OEC clusters. (h) Co-localization of VLMCs with interlaminar astrocytes (ILA) as shown in selected MERFISH sections. (i) UMAP representation of OPCs and oligodendrocytes colored and labeled by supertype. (j-k) Representative MERFISH sections showing the spatial distribution of OPC (j) and Oligo (k) supertypes.CC BY · PubMed Central
Extended Data Fig. 12 Gene expression patterns in immature neuronal populations and RMS astrocytes. (a) UMAP representation of immature neuron populations colored by supertype. Maturation trajectories in dentate gyrus (DG) (b), inner olfactory bulb (c), and outer olfactory bulb (d) are highlighted. (b-d) Representative MERFISH sections showing location of immature neuronal supertypes from the three trajectories shown in (a). (e-g) Heatmap showing gene expression changes as immature neurons transition to mature cell types, conserved between OB (left) and DG (right) cell type development (e), specific to OB cell types (f), and specific to DG cell types (g). Key marker genes at each stage of development are highlighted. It seems, however, that the scRNA-seq data might not have captured all cell states along the DG maturation trajectory based on the gaps between clusters in the UMAP and absence of expression for genes like Ascl1, Pax6, Top2a, and Mki67 along the DG trajectory. Various studies have tried to capture the transitional states between neural stem and neuronal progenitor cells in the DG with most making use of transgenic mice to isolate specific states142,143. (h) MERFISH sections showing the co-localization of immature neurons and astrocytes in the rostral migratory stream (RMS). The dashed boxes in (d) show the location of the highlighted regions in (h). (i) Heatmap showing gene expression changes in astrocytes associated with the RMS from SVZ to OB. Highlighted genes are conserved between the RMS-associated astrocytes and the OB trajectory.
Extended Data Fig. 12 Gene expression patterns in immature neuronal populations and RMS astrocytes. (a) UMAP representation of immature neuron populations colored by supertype. Maturation trajectories in dentate gyrus (DG) (b), inner olfactory bulb (c), and outer olfactory bulb (d) are highlighted. (b-d) Representative MERFISH sections showing location of immature neuronal supertypes from the three trajectories shown in (a). (e-g) Heatmap showing gene expression changes as immature neurons transition to mature cell types, conserved between OB (left) and DG (right) cell type development (e), specific to OB cell types (f), and specific to DG cell types (g). Key marker genes at each stage of development are highlighted. It seems, however, that the scRNA-seq data might not have captured all cell states along the DG maturation trajectory based on the gaps between clusters in the UMAP and absence of expression for genes like Ascl1, Pax6, Top2a, and Mki67 along the DG trajectory. Various studies have tried to capture the transitional states between neural stem and neuronal progenitor cells in the DG with most making use of transgenic mice to isolate specific states142,143. (h) MERFISH sections showing the co-localization of immature neurons and astrocytes in the rostral migratory stream (RMS). The dashed boxes in (d) show the location of the highlighted regions in (h). (i) Heatmap showing gene expression changes in astrocytes associated with the RMS from SVZ to OB. Highlighted genes are conserved between the RMS-associated astrocytes and the OB trajectory.CC BY · PubMed Central
Extended Data Fig. 13 Transcription factor families. Expression of key transcription factors for each subclass in the taxonomy tree, organized by transcription factor gene families. The lines divide the dendrogram into neighborhoods. Source Data
Extended Data Fig. 13 Transcription factor families. Expression of key transcription factors for each subclass in the taxonomy tree, organized by transcription factor gene families. The lines divide the dendrogram into neighborhoods. Source DataCC BY · PubMed Central
Extended Data Fig. 14 Distribution of Gini coefficient for subclasses. (a) Illustration explaining the concept of the Gini coefficient (GC). For each subclass, brain regions are ordered by number of cells present in x-axis. The y-axis is the cumulative fraction of cells for each subclass. The Gini coefficient is calculated by dividing the area (Area A) between the line of perfect equality and the observed distribution curve (the Lorenz curve) by the total area under the line of perfect equality (Areas A + B). The result is a value between 0 and 1 with 0 representing perfect equality and 1 maximum inequality. (b) Distribution of GCs for all subclasses. Color scheme is the same as used for Fig. 6a. (c) Ridge plot showing the distribution of GCs for subclasses grouped by class. (d) 3D example plots of subclasses for 4 classes (02 NP-CT-L6b, 07 CTX-MGE-GABA, 30 Astro-Epen, and 33 Vascular), illustrating the wide range of GCs present. Within each class, plots are ordered by GC from lowest to highest.
Extended Data Fig. 14 Distribution of Gini coefficient for subclasses. (a) Illustration explaining the concept of the Gini coefficient (GC). For each subclass, brain regions are ordered by number of cells present in x-axis. The y-axis is the cumulative fraction of cells for each subclass. The Gini coefficient is calculated by dividing the area (Area A) between the line of perfect equality and the observed distribution curve (the Lorenz curve) by the total area under the line of perfect equality (Areas A + B). The result is a value between 0 and 1 with 0 representing perfect equality and 1 maximum inequality. (b) Distribution of GCs for all subclasses. Color scheme is the same as used for Fig. 6a. (c) Ridge plot showing the distribution of GCs for subclasses grouped by class. (d) 3D example plots of subclasses for 4 classes (02 NP-CT-L6b, 07 CTX-MGE-GABA, 30 Astro-Epen, and 33 Vascular), illustrating the wide range of GCs present. Within each class, plots are ordered by GC from lowest to highest.CC BY · PubMed Central
Extended Data Fig. 15 Predicting transcriptomic subclasses based on spatial location of MERFISH cells. (a) Correspondence between assigned subclasses for MERFISH cells in glutamatergic, dopaminergic and serotonergic subclasses, and predicted subclasses based on the spatial coordinates of these MERFISH cells. Each row is normalized by dividing by the maximum number. Insert in the lower left corner shows the correspondence between assigned and predicted glutamatergic, dopaminergic, and serotonergic classes. Insert in the upper right corner highlights the correspondence between assigned and predicted subclasses in the MB Glut class. (b) Correspondence between assigned GABAergic subclasses and predicted subclasses based on the spatial coordinates of MERFISH cells. Insert in the lower left corner shows the correspondence between assigned and predicted GABAergic classes. Insert in the upper right corner highlights the correspondence between assigned and predicted subclasses in the MB GABA class. Source Data
Extended Data Fig. 15 Predicting transcriptomic subclasses based on spatial location of MERFISH cells. (a) Correspondence between assigned subclasses for MERFISH cells in glutamatergic, dopaminergic and serotonergic subclasses, and predicted subclasses based on the spatial coordinates of these MERFISH cells. Each row is normalized by dividing by the maximum number. Insert in the lower left corner shows the correspondence between assigned and predicted glutamatergic, dopaminergic, and serotonergic classes. Insert in the upper right corner highlights the correspondence between assigned and predicted subclasses in the MB Glut class. (b) Correspondence between assigned GABAergic subclasses and predicted subclasses based on the spatial coordinates of MERFISH cells. Insert in the lower left corner shows the correspondence between assigned and predicted GABAergic classes. Insert in the upper right corner highlights the correspondence between assigned and predicted subclasses in the MB GABA class. Source DataCC BY · PubMed Central
Extended Data Fig. 16 Predicting anatomical regions based on imputed MERFISH transcriptomes. (a) Correspondence between assigned CCFv3 regions for MERFISH cells in glutamatergic, dopaminergic, and serotonergic cell types to predicted CCFv3 regions based on imputed transcriptomes of these MERFISH cells. Each row is normalized by dividing by the maximum number. Insert in the lower left corner shows the correspondence between assigned and predicted regions for glutamatergic, dopaminergic, and serotonergic cell types. Insert in the upper right corner highlights the correspondence between assigned and predicted subregions in midbrain. (b) Correspondence between assigned CCFv3 regions for MERFISH cells in GABAergic cell types and predicted CCFv3 regions based on imputed transcriptomes of these MERFISH cells. Insert in the lower left corner shows the correspondence between assigned and predicted regions for GABAergic cell types. Insert in the upper right corner highlights the correspondence between assigned and predicted subregions in midbrain. Source Data
Extended Data Fig. 16 Predicting anatomical regions based on imputed MERFISH transcriptomes. (a) Correspondence between assigned CCFv3 regions for MERFISH cells in glutamatergic, dopaminergic, and serotonergic cell types to predicted CCFv3 regions based on imputed transcriptomes of these MERFISH cells. Each row is normalized by dividing by the maximum number. Insert in the lower left corner shows the correspondence between assigned and predicted regions for glutamatergic, dopaminergic, and serotonergic cell types. Insert in the upper right corner highlights the correspondence between assigned and predicted subregions in midbrain. (b) Correspondence between assigned CCFv3 regions for MERFISH cells in GABAergic cell types and predicted CCFv3 regions based on imputed transcriptomes of these MERFISH cells. Insert in the lower left corner shows the correspondence between assigned and predicted regions for GABAergic cell types. Insert in the upper right corner highlights the correspondence between assigned and predicted subregions in midbrain. Source DataCC BY · PubMed Central
Extended Data Fig. 17 Distribution of cluster numbers and cluster sizes across different brain regions. (a) Number of clusters per fine CCFv3 region (Supplementary Table 9) as analyzed using the MERFISH data. Bars are colored by broad CCFv3 regions. (b) Distribution of cluster size (number of cells per cluster) per major brain region in scRNA-seq data and MERFISH data.
Extended Data Fig. 17 Distribution of cluster numbers and cluster sizes across different brain regions. (a) Number of clusters per fine CCFv3 region (Supplementary Table 9) as analyzed using the MERFISH data. Bars are colored by broad CCFv3 regions. (b) Distribution of cluster size (number of cells per cluster) per major brain region in scRNA-seq data and MERFISH data.CC BY · PubMed Central
Extended Data Fig. 18 Circadian cycle associated expression changes in clock genes. (a-b) Dot plots showing the expression of clock genes in light-phase and dark-phase cells within each cell class (a) or selected subclasses that have any clock genes with fold change |log2(FC)| > 1 between light and dark phases (b). Dot size and color indicate proportion of expressing cells and average expression level in each class or subclass, respectively. (c) Heatmap showing the log2(FC) difference between light and dark phases for clock genes in the selected subclasses as in (b).
Extended Data Fig. 18 Circadian cycle associated expression changes in clock genes. (a-b) Dot plots showing the expression of clock genes in light-phase and dark-phase cells within each cell class (a) or selected subclasses that have any clock genes with fold change |log2(FC)| > 1 between light and dark phases (b). Dot size and color indicate proportion of expressing cells and average expression level in each class or subclass, respectively. (c) Heatmap showing the log2(FC) difference between light and dark phases for clock genes in the selected subclasses as in (b).CC BY · PubMed Central

Verification

Claims the paper makes

182 Own result · 45 Own interpretation · 38 Own method · 31 No basis given · 28 Partly supported · 19 Supported · 5 No full text of the sources · 3 Supported so far · 3 Not checked yet · 1 Mixed support · 1 Sources not found · accuracy not yet measured

Show all 170 claims

  • It may state two separate claims as one. 140
  • Probably not a substantial claim. 70
  • The quote only partly states this claim. 3
186 of 186
  • The paper relies on this without support.

    The mouse is a suitable model for a comprehensive definition of mammalian brain composition and architecture.

    “The mouse (Mus musculus) is the most widely used mammalian model organism and is therefore a natural first choice for a comprehensive definition of mammalian brain composition and architecture.”
  • The paper relies on this without support.

    The CCFv3-registered MERFISH dataset accurately represents the spatial distribution of cell types.

    “using the CCFv3-registered whole-brain MERFISH dataset”
  • The paper relies on this without support.

    The donor sex assignment is accurate.

    “cells derived from the two sexes”
  • The paper relies on this without support.

    The number of clusters is a proxy for cell-type complexity.

    “indicating a high degree of cell-type complexity”
  • The paper relies on this without support.

    Transcriptomic similarity reflects cell-type similarity.

    “Transcriptomically more similar cell types are located closer to each other spatially”
  • The paper relies on this without support.

    The marker genes used are specific to the cell types they mark.

    “can be distinguished by highly specific marker genes at all levels of hierarchy”
  • The paper relies on this without support.

    The Tau score is a valid measure of cell-type specificity of gene expression.

    “We measured the cell-type specificity of expression of these genes using the Tau score60”
  • The authors state this as an assumption.

    The QC process effectively removes doublets and low-quality cells.

    “Our QC process excluded the possibility of doublet or low-quality cell contamination”
  • The paper relies on this without support.

    The 094 SCH Six6 Cdc14a Gaba subclass is specific to the suprachiasmatic nucleus, which is the circadian pacemaker of the brain.

    “which is specific to the suprachiasmatic nucleus (SCH), the circadian pacemaker of the brain”
  • The paper relies on this without support.

    Noise clusters can confuse the integration algorithm and reduce cell-type resolution.

    “The presence of such clusters can confuse the integration algorithm and reduce the cell-type resolution.”
  • The paper relies on this without support.

    A cluster with at least 50% fewer UMIs or >100 lower QC score, smaller size, and no more than one up-regulated gene relative to another cluster is a degraded version of that other cluster.

    “In these cases, one cluster was a degraded version of another cluster and therefore removed.”
  • The paper relies on this without support.

    Due to significantly increased dataset sizes, the Jaccard similarity graph can be extremely large, impossible to fit in memory.

    “Due to significantly increased dataset sizes, the Jaccard similarity graph can be extremely large, impossible to fit in memory.”
  • The paper relies on this without support.

    Sampling of some hindbrain cell types was not comprehensive enough; they are likely very rare and difficult to collect due to non-neuronal cells and high myelination.

    “First, sampling of these cell types was not comprehensive enough. They are likely very rare and, given the large fraction of non-neuronal cells present in these areas along with high level of myelination, these cells can be very difficult to collect.”
  • The paper relies on this without support.

    Transcriptomic differences in hindbrain neuronal types are subtler than in other brain regions.

    “Finally, transcriptomic differences in hindbrain neuronal types appear to be subtler compared to neuronal types in other brain regions.”
  • The paper relies on this without support.

    Full-scale integration of the entire 10xMulti snRNA-seq dataset could compromise the many high-resolution clusters already present in the scRNA-seq taxonomy.

    “we did not proceed with full-scale integration of the entire 10xMulti snRNA-seq dataset which could compromise the many high-resolution clusters already present in our scRNA-seq taxonomy.”
  • The paper relies on this without support.

    Any nonlinear curve can be accurately approximated using a series of linear segments.

    “as any nonlinear curve can be accurately approximated using a series of linear segments.”
  • The paper relies on this without support.

    Slc6a1 is widely expressed in all GABAergic neurons, as well as astrocytes, and even many subcortical glutamatergic neurons.

    “because Slc6a1 is widely expressed in all GABAergic neurons, as well as astrocytes, and even many subcortical glutamatergic neurons”
17 of 17

Claims other papers make about this one

Each is checked against this paper’s own text: does it say what the citing paper says it does?

Self-check

Percentages, sums and reported differences re-computed from the paper’s own numbers.

Each contradiction between two passages listed here was confirmed by an independent check.

1 of 1

Citations

Cites

Show all 157 references

Research works printed but not linked: 1

  1. Paxinos, G. The Rat Nervous System 4th edn (Academic Press, 2014).

Cited by

Sorted by citations

Items 801–900 of 1,250

Metadata

Authors
Yao, Z., van Velthoven, C. T. J., Kunst, M., Zhang, M., McMillen, D., Lee, C., Jung, W., Goldy, J., Abdelhak, A., Aitken, M., Baker, K., Baker, P., Barkan, E., Bertagnolli, D., Bhandiwad, A., Bielstein, C., Bishwakarma, P., Campos, J., Carey, D., Casper, T., Chakka, A. B., Chakrabarty, R., Chavan, S., Chen, M., Clark, M., Close, J., Crichton, K., Daniel, S., DiValentin, P., Dolbeare, T., Ellingwood, L., Fiabane, E., Fliss, T., Gee, J., Gerstenberger, J., Glandon, A., Gloe, J., Gould, J., Gray, J., Guilford, N., Guzman, J., Hirschstein, D., Ho, W., Hooper, M., Huang, M., Hupp, M., Jin, K., Kroll, M., Lathia, K., Leon, A., Li, S., Long, B., Madigan, Z., Malloy, J., Malone, J., Maltzer, Z., Martin, N., McCue, R., McGinty, R., Mei, N., Melchor, J., Meyerdierks, E., Mollenkopf, T., Moonsman, S., Nguyen, T. N., Otto, S., Pham, T., Rimorin, C., Ruiz, A., Sanchez, R., Sawyer, L., Shapovalova, N., Shepard, N., Slaughterbeck, C., Sulc, J., Tieu, M., Torkelson, A., Tung, H., Valera Cuevas, N., Vance, S., Wadhwani, K., Ward, K., Levi, B., Farrell, C., Young, R., Staats, B., Wang, M.-Q. M., Thompson, C. L., Mufti, S., Pagan, C. M., Kruse, L., Dee, N., Sunkin, S. M., Esposito, L., Hawrylycz, M. J., Waters, J., Ng, L., Smith, K., Tasic, B., Zhuang, X., & Zeng, H.
Affiliations
  • Allen Institute for Brain Science (US)
  • Howard Hughes Medical Institute (US)
  • University of Pennsylvania (US)
  • Genentech (US)
Corresponding author
  • Yao, Zizhen
  • Zeng, Hongkui
Venue
Nature
Published
2023-12-13
Volume
624
Issue
7991
Pages
317-332
ISSN
  • 0028-0836 (ISSN-L)
  • 1476-4687
OpenAlex
W4389670555
Cited by
1,203
References
157
Citations by year
2026: 505 · 2025: 451 · 2024: 215 · 2023: 31
FWCI
92.45
Type
journal article
Language
en
Open access
Locations
  • Nature (publishedVersion, cc-by)
  • bioRxiv (Cold Spring Harbor Laboratory) (submittedVersion, cc-by-nc)
  • PubMed Central (publishedVersion, cc-by)
  • PubMed (publishedVersion)
  • PubMed Central (submittedVersion, cc-by-nc)
  • RePEc: Research Papers in Economics (submittedVersion)
Funders
National Institute of Mental Health (U01MH105982, U19MH114830, R24MH117295, U24MH130918)
Topics
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Single-cell and spatial transcriptomics
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Biophysics › Cell Image Analysis Techniques
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Genomics and Chromatin Dynamics
Keywords
mouse brain, cell type atlas, single-cell transcriptomics, spatial transcriptomics, MERFISH, neuronal cell types, neurotransmitter expression

About these checks

Registry, retraction, Cites and statistics results are computed, not judged, with one exception: a printed reference that no index lists is read once by a model before it counts as a citation error. Claims are extracted by a model and kept only when their quote appears verbatim in the text; each cited source is then checked, with a model’s help, for whether it says what the claim relies on. Those verdicts can be wrong — read the claim’s quote and the cited source before relying on one.

A language model drew the claims from the paper on 2026-10-07, under version 3 of the claim definition (model deepseek-v4.1-flash, prompt extract-v3).

A second model labelled each claim on 2026-10-07 (model jev-1.13-20260917, definition revision gold-guideline-v0.2).

Cited sources are checked with a language model’s help (gpt-6-sol).

Statistics and arithmetic are recomputed, not judged.

Confidence bands are shown only for fields calibrated on our annotated test set. None is calibrated yet, so no band is shown.

What counts as a claim ›How we measure confidence ›How we check this ›

Are you an author, or do you see a problem with this page? Request a review