Sign in

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

In-depth summary

✦ Multi-agent verified full-text synthesis

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

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

  2. 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

  3. 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

  4. 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

  5. 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
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.

Similarity Map
SelectedSame fieldCross-field
THISRen '19Ren '18Awasthi '20Haiman '23Hingorani '22Hänisch '23
70% semantic matchNeuroscience
Jing Ren, Alina Isakova, Drew Friedmann+14 more·eLife·Oct 24, 2019· 304 citations