Architecture of the Serotonergic Projectome in the Mouse Brain
Jun Ho Song, Drew FriedmanniD, Yunming WuiD, Ashley Moses, Jalal BaruniiD, Yuan YuaniD, Isabella Rana, Kelly Zhao, Tom Hindmarsh SteniD, Qian WangiD, Shuyun Alina XiaoiD, Xiaoke CheniD, Scott W. LindermaniD, Liqun LuoiD
Abstract
The serotonin system innervates nearly the entire brain to support diverse functions, and its dysfunction is implicated in multiple psychiatric disorders. Although recent studies suggested that this anatomically diffuse system supports differentiated rather than uniform modulation, its overall organization is unclear. Here, using systematic whole-brain axon tracing and integrated analyses, we show that the serotonin projectome is organized by functional relatedness rather than physical proximity of its targets. Dorsal and median raphe serotonin neurons partition into five projectomic groups, each preferentially innervating the hippocampus, basal ganglia, cortex, medial interbrain, or brainstem/lateral thalamus, with within-group structure ranging from discrete subgroups to continuous variation. Spatial transcriptomic analyses reveal that each projectomic group corresponds to a distinct combination of transcriptomic clusters, with transcriptome and somatic position jointly predicting projectomic identity. Projection-specific serotonin depletion produces dissociable behavioral phenotypes. Together, projectomic identity emerges as a central axis linking axon collateralization, molecular diversity, and behavioral function.
02 /Specification & Classification
Click any tag to browse matching papers- Citations
- 1
- FWCI
- 2.51
- Field percentile
- top 8%
- Organism
- Assay
- Research Focus
- Approach
In-depth summary
Claude Opus 5.5 · xhigh effort · 8 Nov 2026 · from the full text of version 2 (3 Sept 2026) · CC BY 4.0
Whole-brain axon tracing in 110 mouse brains divides raphe serotonin neurons into five projectomic groups, tied to transcriptome plus position and to dissociable behavioral effects of depletion.
Most experiments on serotonin neurons have treated them (mostly from the dorsal raphe) as one population, and recent work on their diversity covered only small subsets. This paper maps collateralization systematically across the brain and links it to molecular identity and behavior.
Key results
ShownSerotonin projections from dorsal and median raphe (DR, MR) are well described by five basis patterns. Hierarchical clustering of brains in this space gives five projectomic groups: hippocampal, basal ganglia, cortical, medial interbrain, and brainstem/lateral thalamic.
Matrix of 110 brains x 280 regions. Held-out reconstruction error from speckled cross-validation plateaued at k = 5, at the replicate-variability baseline and well below the shuffle baseline. The factorization explained 73.5% of variance. 101 of 110 brains kept their group under 280-D clustering; 79.2% of 250 subsampling trials gave five groups. · Determining whole-brain serotonin axon…; DR and MR serotonin neurons comprise five…; Axon data analysis using unsupervised NMF; Figure legends
SuggestedGroup membership appears to follow functional relatedness rather than physical proximity. Central amygdala injections fell in the basal ganglia group but basolateral amygdala injections in the cortical group. VTA and nucleus accumbens core injections fell in the medial interbrain group, not the basal ganglia group.
Central amygdala injections clustered with basal ganglia, basolateral amygdala with cortical; VTA and accumbens core with medial interbrain, not basal ganglia. The authors call this consistent with functional divergence. Same-group different-site cosine similarity matched same-site similarity (except cortical) and exceeded between-group similarity. Marker regions per group: 16, 8, 18, 28, 81. · DR and MR serotonin neurons comprise five…; Organization of the mouse brain from the…; Marker region analysis
ShownWithin groups, structure runs from discrete subgroups to continuous variation. Basal ganglia innervation splits into subpatterns that pair topographically across caudoputamen and globus pallidus externus, while cortical group brains vary continuously.
NMF gave k = 4 subpatterns in caudoputamen (17 brains) and k = 3 in globus pallidus externus, with dorsal, ventral and caudal subpatterns correlated across the two regions. The 27 cortical-group brains (k = 3) did not form distinct clusters. No significant subpatterns for the medial interbrain or brainstem/lateral thalamic groups. · Within-group organization: from discrete…; Axon data analysis using unsupervised NMF; Figure legends
SuggestedTranscriptomic cluster and somatic position jointly predict group identity. Each group maps to a combination of transcriptomic clusters rather than a single one. Neurons collateralize preferentially within a group.
Spatial transcriptomics with retrograde barcodes: 8,715 Tph2+ neurons, 635 barcode-positive. Multi-barcode rate 15.2% (62/408) in within-group mice vs 4.8% (11/227) in cross-group mice. Random forest on single barcoded cells: 74.0% accuracy combined vs 41.4% transcriptome-only and 58.4% position-only. Supervised NMF on raphe cell-body maps: 83.6% (chance 20%). · Projection-specific cell body distributions in…; Transcriptomic and spatial identity jointly…; STARmap data processing
SuggestedProjection-specific Tph2 depletion gives dissociable behavioral profiles across eight paradigms. The basal ganglia group diverged most. A multivariate discriminant model separated groups and generalized to a second batch using different marker regions.
54 features. Basal ganglia group: fewer marbles buried, more food intake after fasting, less immobility in forced swim. Cross-validated LDA accuracy 79% (females) and 60% (males) vs permutation chance 25.9 ± 7.2% and 26.1 ± 7.0%. Batch 2 within-group similarity exceeded between-group similarity in all groups in both sexes, except controls. · Projection-specific Tph2 depletion produces…; Serotonergic projectomic groups organize the…; Multivariate behavioral classification using…
Methods, data and reproducibility
Projection-defined serotonin neurons were labeled with retrograde AAV in intersectional reporter mice. Cleared brains were imaged by light-sheet microscopy, axons segmented and registered to Allen CCFv3, and NMF with clustering was applied. This was combined with STARmap plus retrograde barcodes, and projection-specific Tph2 deletion with video-based behavior scoring and LDA.
- Light-sheet whole-brain axon maps of projection-defined serotonin neurons, quantified in 280 CCFv3 regions · mouse · n = 110 brains; 35 injection regions · doi 10.35077/g.1213 (Brain Image Library, raw images)
- STARmap with retrograde circular barcodes in DR and MR · mouse · n = 14 mice; 8,715 Tph2+ neurons; 635 barcode-positive
- Behavior after projection-specific Tph2 deletion, Batch 1, sample sizes given as (male, female) per injected region plus controls · mouse · n = dCA3 (5, 4), GPe (6, 4), PIR (4, 6), LPO (5, 5), SC (5, 5), control (5, 5)
- Behavior, Batch 2 with different marker regions, sample sizes given as (male, female) per injected region plus controls · mouse · n = vCA1 (2, 2), SNr (3, 1), AI (3, 1), PVT (1, 3), LG (2, 2), control (3, 3)
- Piezo2-Cre;Sert-Flp intersectional validation of Cluster 11 projections · mouse · n = 3 mice; 102 mCherry+ neurons
- Methods & tools
- TrailMap, Allen CCFv3, Elastix, Ilastik, NMF (MATLAB), STARmap, Cellpose, Baysor, Harmony, MoSeq, SlowFast, YOLOv8
- Data & accessions
- https://doi.org/10.35077/g.1213
- Software & environment
- Analyses in MATLAB R2021b; behavior pipeline in Python with PyTorch, PyTorchVideo, Ultralytics YOLOv8 and OpenCV; Harmony, scanpy and umap-learn for STARmap. No container or requirements file mentioned.
Limitations (5)
- Projections were mapped by bulk tracing of neurons sharing a target, which can miss heterogeneity inside groups. Additional groups or subgroups might appear with injections in more regions. · stated by the authors
- The cytosolic reporter does not separate terminals from axons-in-passage; axons near terminals may still contribute to the projectome. · stated by the authors
- Behavior used only loss of function (Tph2 deletion); no gain-of-function or in vivo recordings of the targeted populations. · stated by the authors
- Batch 2 behavior groups are very small (1 to 3 animals per sex per group), so the cross-batch generalization rests on few animals. · noted in this summary
- The transcriptome-projection analysis relies on 439 consensus-classified barcode-positive cells from 14 mice, and 430 single-group cells for the classifier. · noted in this summary
Who should read:Read this if you study neuromodulator circuit organization, or want a worked example of NMF-based projectome analysis integrated with spatial transcriptomics and multivariate behavioral phenotyping.
Reproducibility ledger
Nobody has reported trying to reproduce this paper yet. Did you replicate an experiment or run the code? Your report saves the next lab a week.
Sign in to file a reproducibility report.
Related papers
- Serotonin Neurons in the Dorsal and Medial Raphe Nuclei: from Single-Cell Transcriptomes to Whole-Brain Projections
Jing Ren, Alina Isakova, … Liqun Luo·eLifePeer-reviewed·Oct 24, 2019
- Anatomical, Physiological, and Functional Heterogeneity of the Dorsal Raphe Serotonin System
Jing Ren, Drew Friedmann, … Liqun Luo·bioRxiv·Jan 31, 2018
- Comprehensive Topographical Map of the Serotonergic Fibers in the Mouse Brain
Janak R. Awasthi, Kota Tamada, … Toru Takumi·The Journal of Comparative NeurologyPeer-reviewed·Sep 6, 2020
- Qualitative Observations of Transgenically Individualized Serotonergic Fibers in the Mouse Telencephalon
Justin H. Haiman, Kasie C. Mays, … Skirmantas Janušonis·bioRxiv·Nov 25, 2023
- High-Resolution Spatiotemporal Analysis of Single Serotonergic Axons in an In Vitro System
Melissa Hingorani, Adele M.L. Viviani, … Skirmantas Janušonis·Frontiers in NeurosciencePeer-reviewed·Oct 24, 2022
- NEUROTRANSMITTER TRANSPORTER/RECEPTOR CO-EXPRESSION SHARES ORGANIZATIONAL TRAITS WITH BRAIN STRUCTURE AND FUNCTION
Benjamin Hänisch, Justine Y. Hansen, … Sofie Louise Valk·eLifePeer-reviewed·Jul 13, 2023