EASI-PASS: An accessible pipeline for linking functional imaging and mRNA profiling
Jonnathan Singh AlvaradoiD, Crystian I MassengilliD, Josh Stern, Oren AmsalemiD, Bettina Francheska Ventura, Ahram JangiD, Sarah CookiD, Ana Veliche, Praneel S. Sunkavalli, Diti PateliD, Juliana Colaccino, Kathryn E. EvansiD, Yuhan WangiD, Mark L. AndermanniD
Abstract
We developed EASI-PASS, a reliable, high-throughput method for estimating the molecular identity of functionally characterized cells by merging live imaging with subsequent fixed-tissue imaging using conventional microscopes. Our method matches the shapes and locations of thousands of densely imaged cells between large (>1 mm2) functional images and a thick, expanded, and cleared EASI-FISH tissue volume to assess gene expression. This approach is more efficient than alignment to thin sections and recovers the molecular identity of ~78% of cells. In acute brain slice imaging from the mouse parabrachial nucleus during optogenetic stimulation of long-range spinal inputs, we observed fine-scale specificity in the molecular identity of spinorecipient neurons. In the awake mouse visual cortex, we observed distinct arousal modulation and spatial falloff in correlations within and across interneuron classes. Thus, EASI-PASS provides reliable and efficient alignment of cellular activity with molecular identity.
02 /Specification & Classification
Click any tag to browse matching papers- Citations
- 0
- Organism
- Assay
- Research Focus
- Approach
In-depth summary
Claude Opus 5.5 · xhigh effort · 8 Nov 2026 · from the full text of version 2 (28 Aug 2026) · CC BY 4.0
Matching cell masks between two-photon calcium images and thick, expanded EASI-FISH volumes recovers molecular identity for 72% (slice) to 83% (in vivo) of imaged neurons on conventional microscopes.
Linking calcium imaging to post-hoc mRNA labelling in the same cells is hampered by tissue deformation and cross-modal differences in cell appearance, and existing strategies typically need specialized hardware and yield few matched cells from small (<0.5 mm) images. EASI-PASS instead matches cells by segmented-mask location and shape, not pixel-wise alignment, against one intact expanded slab (up to 500 µm thick) rather than many thin sections, using standard two-photon and confocal microscopes.
Key results
ShownEASI-PASS found a matching EASI-FISH cell for most functionally imaged neurons in both preparations, giving thousands of matched cells per experiment.
In vivo visual cortex: 22,421 of 27,152 two-photon cells matched (82.6%; n = 5 mice). Acute parabrachial slices: 14,728 of 20,397 cells matched (72.2%; 7 slabs from 7 mice). Each in vivo experiment yielded ∼150-250 simultaneously recorded interneurons per class plus ∼4,000 putative excitatory neurons. · Abstract; Cell matching across 2P and EASI-FISH; Distinct arousal tuning and spatial correlations…; Figure legends
ShownAfter manual landmarking, stepwise refinement by mask intersection over union (IoU) improved alignment in every sample. Soma-print, which matches cells by the geometry of their neighbours, recovered more cells than a fixed IoU threshold and mostly picked the same partners.
Compared with the landmark-only baseline, the final alignment raised median IoU (paired t test, slice p = 1×10⁻⁵, in vivo p = 0.002), the Soma-print matched fraction (p = 0.001; p = 0.002) and the AUC separating best from second-best scores (p = 0.03; p = 0.008). Soma-print and IoU > 0.3 chose the same partner in ∼94% of cells both matched; Soma-print matched ∼9% more cells. · Cross-modal registration approach; Cell matching across 2P and EASI-FISH; Figure legends
ShownIn mouse parabrachial nucleus (PBN) slices, neurons excited by optogenetic stimulation of lumbar spinoparabrachial axons were enriched for Syt10. Syt10 stayed the best predictor of the evoked response after accounting for each neuron's position and local axon density.
1,242 of 9,024 matched PBN neurons (13.7%; 6 mice) were activated. Syt10 (p = 0.006) and Th (p = 0.007) were higher in activated neurons; % activated was higher in the top vs bottom decile (Syt10 p = 0.007, Th p = 0.012). Syt10 led both conditional probabilities and a mixed model controlling for position and ChrimsonR signal; Syt10+ cells out-responded low-expressing neighbours (median 45.3 µm away). · Spinoparabrachial input recruits a molecularly…; Regression against neuron location and ChrimsonR…; Comparison with neighboring neurons; Figure legends
ShownIn awake mouse visual cortex, coupling to pupil-linked arousal differed across Pvalb, Sst and Vip interneurons recorded at the same time. Calb2+ VIP neurons were more positively coupled than Calb2- VIP neurons.
Pvalb+ cells split roughly evenly between positive and negative pupil correlation, Vip+ cells were mostly positive and Sst+ cells were intermediate (within-mouse permutation Kruskal-Wallis H = 190.79, p < 0.001, n = 2,385 cells; Friedman χ² = 11.10, p = 0.011, 4 mice). Calb2+ minus Calb2- mean r: VIP 0.092 (p < 0.001, 680 neurons), SST 0.062 (p = 0.162, 600 neurons), 4 mice. · Distinct arousal tuning and spatial correlations…; Statistical testing; Figure legends
SuggestedCorrelations between neurons of the same class fell off with distance most slowly for VIP, at an intermediate rate for SST and steeply for PV and putative excitatory neurons. Calb2+ SST and VIP subtypes had longer decay lengths, and longer decay went with stronger arousal coupling mainly in VIP cells.
The decay constant λ differed across classes (H = 107.10, p < 0.001, n = 1,864 neurons; Friedman χ² = 12.00, p = 0.007, 4 mice), but the fitted amplitude did not (H = 5.16, p = 0.082). Calb2+ vs Calb2- Δmean log λ: SST 0.386, VIP 0.638 (both p < 0.001; 456 and 360 neurons, 3 mice). Within each class, the λ-arousal correlation was weaker in PV and SST and reversed sign in unlabeled neurons. · Distinct arousal tuning and spatial correlations…; Spatial correlation decay; Statistical testing; Figure legends
Methods, data and reproducibility
Two-photon planes (soma-restricted Ribo-GCaMP8s) are warped into a thick, cleared, expanded EASI-FISH volume from manual BigWarp landmarks on a sparse tdTomato marker, then a 3D rigid search, global affine and progressively local affine fits maximizing Cellpose4 mask IoU. Cells are matched by IoU or Soma-print, later FISH rounds are registered to the first via cytoDAPI, and per-cell mRNA intensity is tabulated.
- In vivo three-plane two-photon calcium imaging of layer 2/3 visual cortex (lateral V1, AL, LM) during spontaneous behaviour with pupil tracking, then EASI-FISH for Pvalb, Sst, Vip and Calb2 · mouse (Sst-Cre;Ai14, both sexes) · n = 27,152 imaged cells; 5 mice
- Acute coronal slice two-photon imaging of the parabrachial nucleus during optogenetic stimulation of lumbar spinoparabrachial axons, then up to seven EASI-FISH rounds (∼20 marker genes) · mouse (wild-type males, 9 to 15 weeks) · n = 20,397 imaged cells from 7 mice; 9,024 PBN neurons from 6 mice analysed
- Methods & tools
- EASI-PASS, EASI-FISH, HCR v3.0 (Molecular Instruments), BigWarp (ImageJ), Cellpose4, Soma-print, Suite2p, OASIS, DeepLabCut, Ribo-GCaMP8s, ChrimsonR-tdTomato, Leica SP5 confocal microscope
- Experimental conditions
- Global round-to-round FISH registration takes <20 seconds per registration on a standard computer; confocal acquisition takes about 1.5 h for 5 channels and 4 tiles; landmark placement usually takes under 30 minutes when tissue orientation records are kept.
- Reproducibility gaps
- The text has no data availability statement or accessions. Three steps need manual input: BigWarp landmarks, curation of the 2P masks, and checking the round-to-round FISH alignment. Whether the custom Cellpose4 models are available is not stated.
Limitations (6)
- Placing the initial landmarks by hand can be time-consuming unless accurate records of tissue orientation are kept through post-hoc processing (with such records, usually under 30 minutes). · stated by the authors
- Registration quality depends on uniform expansion and complete clearing during EASI-FISH; the protocol may need adapting for other brain tissues. · stated by the authors
- mRNA is summarized as mean fluorescence per cell rather than counted puncta. This is adequate for high-abundance cell-type markers, but lower-abundance genes such as G-protein coupled receptors may need single-molecule counting. · stated by the authors
- The authors cannot confirm that optogenetically evoked PBN activation comes only from monosynaptic spinoparabrachial connections; they suggest future experiments could test this by combining EASI-FISH with anterograde transsynaptic tracing. · stated by the authors
- Per-animal cortical tests rest on 4 mice (3 for the Calb2 decay comparison). The SST Calb2+ arousal difference was not significant (p = 0.162), although the text says Calb2+ cells were more arousal-coupled within each class. · noted in this summary
- Match correctness rests on internal metrics (IoU, Soma-print scores, mutual nearest neighbours); the tdTomato overlap (over 75% of Sst+ vs 6.8% of Pvalb+ cells) used all EASI-FISH cells, not only matched ones, so it does not test the matches. No comparison with other published matching methods is reported. · noted in this summary
Who should read:Read this if you want to tie calcium activity to marker-gene expression for thousands of neurons per experiment with standard two-photon and confocal microscopes, or study parabrachial spinal inputs or cortical interneuron correlations.
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
- A Scalable Cross-Modal Workflow for Linking In Vivo Neural Activity to Neuronal Cell Types via Spatial Gene Profiling
Risa Yamazaki, Mark Eddison, … Yuhan Wang·bioRxiv·Aug 27, 2026
- Discovery of neuronal cell types by pairing whole cell reconstructions with RNA expression profiles
Tiago Ferreira, Mark Eddison, … MouseLight Project Team·bioRxiv·Dec 30, 2024
- Expansion-Assisted Iterative-FISH defines lateral hypothalamus spatio-molecular organization
Yuhan Wang, Mark Eddison, … Paul W. Tillberg·bioRxiv·Mar 8, 2021
- Correlating cell function and morphology by performing fluorescent immunocytochemical staining on the light-microscope stage
Hiroyuki Kawano, Yasuhiro Kakazu, … Nobutoshi C. Harata·bioRxiv·Jun 30, 2020
- 2P-NucTag: on-demand phototagging for molecular analysis of functionally identified cortical neurons
Jingcheng Shi, Boaz Nutkovich, … Ivo Spiegel·NeuronPeer-reviewed·Apr 2, 2026
- A transcriptomic axis aligns with in vivo functional dynamics in hippocampal inhibitory circuits
Hyun Choong Yong, Stephanie A. Herrlinger, … Attila Losonczy·bioRxiv·Apr 8, 2026