Abstract
Human neuroimaging studies report that psychedelics induce serotonin-2A receptor-dependent changes in functional brain reorganization, presumably reflecting neuromodulation. However, these studies often overlook the potent vasoactive effects of serotonin. Here we identified psilocybin-induced alterations in hemodynamic response functions during human functional magnetic resonance imaging, suggesting potential disruptions in neurovascular coupling. We then used wide-field optical imaging in awake Thy1-jRGECO1a mice to determine whether psychedelic-induced changes in hemodynamics arise from neuronal, vascular or neurovascular effects. Exposure to the psychedelic 2,5-dimethoxy-4-iodoamphetamine (DOI) differentially altered coupling between cortical excitatory neuronal versus hemodynamic activity, both during whisker stimulation and in the resting state. Furthermore, DOI resulted in discordant changes between neuronal-based versus hemodynamic-based assessments of functional connectivity. A selective serotonin-2A receptor antagonist (MDL100907) reversed many of the effects of DOI. Our results demonstrate a dissociation between DOI-induced neuronal and hemodynamic signals, indicating a need to consider neurovascular effects of psychedelics when interpreting blood-based measures of brain function.
This is a preview of subscription content, access via your institution
Access options
Access Nature and 54 other Nature Portfolio journals
Get Nature+, our best-value online-access subscription
$32.99 / 30 days
cancel any time
Subscribe to this journal
Receive 12 print issues and online access
$259.00 per year
only $21.58 per issue
Buy this article
- Purchase on SpringerLink
- Instant access to full article PDF
Prices may be subject to local taxes which are calculated during checkout
Similar content being viewed by others
Data availability
The Allen Mouse Brain Atlas was downloaded from https://alleninstitute.github.io/AllenSDK/reference_space.html. Processed WFOI data are available via Zenodo at https://doi.org/10.5281/zenodo.15857641 (ref. 155). Raw HTR data are available via Zenodo at https://doi.org/10.5281/zenodo.15857233 (ref. 156). All other data supporting the findings of this study will be made available upon request.
Code Availability
MATLAB processing code is available via GitHub at https://github.com/BauerLabCodebase/WFOI-Textbook-Chapter.
References
Thomas, K., Malcolm, B. & Lastra, D. Psilocybin-assisted therapy: a review of a novel treatment for psychiatric disorders. J. Psychoact. Drugs 49, 446–455 (2017).
Madras, B. K. Psilocybin in treatment-resistant depression. N. Engl. J. Med. 387, 1708–1709 (2022).
Zafar, R. et al. Psychedelic therapy in the treatment of addiction: the past, present and future. Front. Psychiatry 14, 1183740 (2023).
Watts, R. et al. Patients’ accounts of increased ‘connectedness’ and ‘acceptance’ after psilocybin for treatment-resistant depression. J. Humanist Psychol. 57, 520–564 (2017).
Tagliazucchi, E. et al. Increased global functional connectivity correlates with LSD-induced ego dissolution. Curr. Biol. 26, 1043–1050 (2016).
Viol, A. et al. Shannon entropy of brain functional complex networks under the influence of the psychedelic Ayahuasca. Sci. Rep. 7, 7388 (2017).
Iadecola, C. The neurovascular unit coming of age: a journey through neurovascular coupling in health and disease. Neuron 96, 17–42 (2017).
Rapport, M. M., Green, A. A. & Page, I. H. Serum vasoconstrictor, serotonin; isolation and characterization. J. Biol. Chem. 176, 1243–1251 (1948).
Carhart-Harris, R. L. & Nutt, D. J. Serotonin and brain function: a tale of two receptors. J. Psychopharmacol. 31, 1091–1120 (2017).
Vollenweider, F. X. & Kometer, M. The neurobiology of psychedelic drugs: implications for the treatment of mood disorders. Nat. Rev. Neurosci. 11, 642–651 (2010).
Dyer, D. C. & Gant, D. W. Vasoconstriction produced by hallucinogens on isolated human and sheep umbilical vasculature. J. Pharm. Exp. Ther. 184, 366–375 (1973).
Glennon, R. A., Titeler, M. & McKenney, J. D. Evidence for 5-HT2 involvement in the mechanism of action of hallucinogenic agents. Life Sci. 35, 2505–2511 (1984).
Xu, T. & Pandey, S. C. Cellular localization of serotonin(2A) (5HT(2A)) receptors in the rat brain. Brain Res. Bull. 51, 499–505 (2000).
Cohen, Z. et al. Serotonin in the regulation of brain microcirculation. Prog. Neurobiol. 50, 335–362 (1996).
Subramanian, S. et al. Psilocybin’s acute and persistent brain effects: a precision imaging drug trial. Sci. Data 12, 941 (2025).
Siegel, J. S. et al. Psilocybin desynchronizes the human brain. Nature 632, 131–138 (2024).
Padawer-Curry, J. A. et al. Wide-field optical imaging in mouse models of ischemic stroke. Methods Mol. Biol. 2616, 113–151 (2023).
Dearnley, B. et al. Stretching and squeezing of neuronal log firing rate distribution by psychedelic and intrinsic brain state transitions. Preprint at bioRxiv https://doi.org/10.1101/2021.08.22.457198 (2021).
Hires, S. A., Tian, L. & Looger, L. L. Reporting neural activity with genetically encoded calcium indicators. Brain Cell Biol. 36, 69–86 (2008).
Warwick, C. et al. Cell type-specific calcium imaging of central sensitization in mouse dorsal horn. Nat. Commun. 13, 5199 (2022).
Celada, P. et al. The hallucinogen DOI reduces low-frequency oscillations in rat prefrontal cortex: reversal by antipsychotic drugs. Biol. Psychiatry 64, 392–400 (2008).
Tagliazucchi, E. et al. Baseline power of theta oscillations predicts mystical-type experiences induced by DMT in a natural setting. Front. Psychiatry 12, 720066 (2021).
Timmermann, C. et al. Human brain effects of DMT assessed via EEG-fMRI. Proc. Natl Acad. Sci. USA 120, e2218949120 (2023).
Wang, X. et al. Spatiotemporal relationships between neuronal, metabolic, and hemodynamic signals in the awake and anesthetized mouse brain. Cell Rep. 43, 114723 (2024).
Ma, Y. et al. Wide-field optical mapping of neural activity and brain haemodynamics: considerations and novel approaches. Philos. Trans. R. Soc. Lond. B 371, 20150360 (2016).
Sunil, S. et al. Neurovascular coupling is preserved in chronic stroke recovery after targeted photothrombosis. Neuroimage Clin. 38, 103377 (2023).
Yuan, R. et al. Regional homogeneity of resting-state fMRI contributes to both neurovascular and task activation variations. Magn. Reson. Imaging 31, 1492–1500 (2013).
Pogorelov, V. M. et al. The G protein biased serotonin 5-HT2A receptor agonist lisuride exerts anti-depressant drug-like activities in mice. Front. Mol. Biosci. 10, 1233743 (2023).
Palhano-Fontes, F. et al. The psychedelic state induced by Ayahuasca modulates the activity and connectivity of the default mode network. PLoS ONE 10, e0118143 (2015).
Lebedev, A. V. et al. Finding the self by losing the self: neural correlates of ego-dissolution under psilocybin. Hum. Brain Mapp. 36, 3137–3153 (2015).
Carhart-Harris, R. L. et al. Neural correlates of the LSD experience revealed by multimodal neuroimaging. Proc. Natl Acad. Sci. USA 113, 4853–4858 (2016).
Carhart-Harris, R. L. et al. The entropic brain: a theory of conscious states informed by neuroimaging research with psychedelic drugs. Front. Hum. Neurosci. 8, 20 (2014).
Timmermann, C. et al. Neural correlates of the DMT experience assessed with multivariate EEG. Sci. Rep. 9, 16324 (2019).
Carhart-Harris, R. L. et al. Neural correlates of the psychedelic state as determined by fMRI studies with psilocybin. Proc. Natl Acad. Sci. USA 109, 2138–2143 (2012).
Whitesell, J. D. et al. Regional, layer, and cell-type-specific connectivity of the mouse default mode network. Neuron 109, 545–559 (2021).
Zirkel, R. T. et al. Psilocybin prolongs the neurovascular coupling response in mouse visual cortex. Preprint at bioRxiv https://doi.org/10.1101/2025.07.25.666803 (2025).
Spain, A. et al. Neurovascular and neuroimaging effects of the hallucinogenic serotonin receptor agonist psilocin in the rat brain. Neuropharmacology 99, 210–220 (2015).
Zhang, G. & Stackman, R. W. Jr. The role of serotonin 5-HT2A receptors in memory and cognition. Front. Pharm. 6, 225 (2015).
Willins, D. L., Deutch, A. Y. & Roth, B. L. Serotonin 5-HT2A receptors are expressed on pyramidal cells and interneurons in the rat cortex. Synapse 27, 79–82 (1997).
Gattuso, J. J. et al. Default mode network modulation by psychedelics: a systematic review. Int J. Neuropsychopharmacol. 26, 155–188 (2023).
Kometer, M. et al. Activation of serotonin 2A receptors underlies the psilocybin-induced effects on alpha oscillations, N170 visual-evoked potentials, and visual hallucinations. J. Neurosci. 33, 10544–10551 (2013).
Muthukumaraswamy, S. D. et al. Broadband cortical desynchronization underlies the human psychedelic state. J. Neurosci. 33, 15171–15183 (2013).
Lambe, E. K. & Aghajanian, G. K. Prefrontal cortical network activity: opposite effects of psychedelic hallucinogens and D1/D5 dopamine receptor activation. Neuroscience 145, 900–910 (2007).
McCulloch, D. E. et al. Psychedelic resting-state neuroimaging: a review and perspective on balancing replication and novel analyses. Neurosci. Biobehav Rev. 138, 104689 (2022).
Muller, F. et al. Altered network hub connectivity after acute LSD administration. Neuroimage Clin. 18, 694–701 (2018).
Daws, R. E. et al. Increased global integration in the brain after psilocybin therapy for depression. Nat. Med. 28, 844–851 (2022).
Dai, R. et al. Classical and non-classical psychedelic drugs induce common network changes in human cortex. Neuroimage 273, 120097 (2023).
Jakab, R. L. & Goldman-Rakic, P. S. Segregation of serotonin 5-HT2A and 5-HT3 receptors in inhibitory circuits of the primate cerebral cortex. J. Comp. Neurol. 417, 337–348 (2000).
Zonta, M. et al. Neuron-to-astrocyte signaling is central to the dynamic control of brain microcirculation. Nat. Neurosci. 6, 43–50 (2003).
Takano, T. et al. Astrocyte-mediated control of cerebral blood flow. Nat. Neurosci. 9, 260–267 (2006).
Attwell, D. et al. Glial and neuronal control of brain blood flow. Nature 468, 232–243 (2010).
Lee, L. et al. Key aspects of neurovascular control mediated by specific populations of inhibitory cortical interneurons. Cereb. Cortex 30, 2452–2464 (2020).
Krawchuk, M. B. et al. Optogenetic assessment of VIP, PV, SOM and NOS inhibitory neuron activity and cerebral blood flow regulation in mouse somato-sensory cortex. J. Cereb. Blood Flow. Metab. 40, 1427–1440 (2020).
Vo, T. T. et al. Parvalbumin interneuron activity drives fast inhibition-induced vasoconstriction followed by slow substance P-mediated vasodilation. Proc. Natl Acad. Sci. USA 120, e2220777120 (2023).
Uhlirova, H. et al. Cell type specificity of neurovascular coupling in cerebral cortex. eLife 5, e14315 (2016).
Lee, J. et al. Opposed hemodynamic responses following increased excitation and parvalbumin-based inhibition. J. Cereb. Blood Flow. Metab. 41, 841–856 (2021).
Moreau, A. W. et al. Serotoninergic fine-tuning of the excitation-inhibition balance in rat visual cortical networks. Cereb. Cortex 20, 456–467 (2010).
Pi, H. J. et al. Cortical interneurons that specialize in disinhibitory control. Nature 503, 521–524 (2013).
Markram, H. et al. Interneurons of the neocortical inhibitory system. Nat. Rev. Neurosci. 5, 793–807 (2004).
Han, K. et al. Excitation-inhibition imbalance leads to alteration of neuronal coherence and neurovascular coupling under acute stress. J. Neurosci. 40, 9148–9162 (2020).
Ruff, C. F. et al. Long-range inhibitory neurons mediate cortical neurovascular coupling. Cell Rep. 43, 113970 (2024).
Pokorny, T. et al. Modulatory effect of the 5-HT1A agonist buspirone and the mixed non-hallucinogenic 5-HT1A/2A agonist ergotamine on psilocybin-induced psychedelic experience. Eur. Neuropsychopharmacol. 26, 756–766 (2016).
Inserra, A. Hypothesis: the psychedelic ayahuasca heals traumatic memories via a sigma 1 receptor-mediated epigenetic-mnemonic process. Front. Pharm. 9, 330 (2018).
Jinks, T. (ed.) Psychological Perspectives on Reality, Consciousness and Paranormal Experience (Palgrave Macmillan, 2019).
Hoyer, D., Hannon, J. P. & Martin, G. R. Molecular, pharmacological and functional diversity of 5-HT receptors. Pharm. Biochem Behav. 71, 533–554 (2002).
Butler, J. J. et al. Classical psychedelics’ action on brain monoaminergic systems. Int. J. Biochem. Cell Biol. 176, 106669 (2024).
Lawn, T. et al. Differential contributions of serotonergic and dopaminergic functional connectivity to the phenomenology of LSD. Psychopharmacology (Berl.) 239, 1797–1808 (2022).
Fontanilla, D. et al. The hallucinogen N,N-dimethyltryptamine (DMT) is an endogenous sigma-1 receptor regulator. Science 323, 934–937 (2009).
Glatfelter, G. C. et al. Comparative pharmacological effects of lisuride and lysergic acid diethylamide revisited. ACS Pharmacol. Transl. Sci. 7, 641–653 (2024).
Muñoz, M. F., Puebla, M. & Figueroa, X. F. Control of the neurovascular coupling by nitric oxide-dependent regulation of astrocytic Ca2+ signaling. Front. Cell. Neurosci. 9, 59 (2015).
Banerjee, A. A. & Vaidya, V. A. Differential signaling signatures evoked by DOI versus lisuride stimulation of the 5-HT2A receptor. Biochem. Biophys. Res. Commun. 531, 609–614 (2020).
Kim, K. J. et al. Vasculo-neuronal coupling: retrograde vascular communication to brain neurons. J. Neurosci. 36, 12624–12639 (2016).
Moore, C. I. & Cao, R. The hemo-neural hypothesis: on the role of blood flow in information processing. J. Neurophysiol. 99, 2035–2047 (2008).
Drew, P. J. Neurovascular coupling: motive unknown. Trends Neurosci. 45, 809–819 (2022).
Phillips, A. A. et al. Neurovascular coupling in humans: physiology, methodological advances and clinical implications. J. Cereb. Blood Flow. Metab. 36, 647–664 (2016).
Palmer, R. M., Ferrige, A. G. & Moncada, S. Nitric oxide release accounts for the biological activity of endothelium-derived relaxing factor. Nature 327, 524–526 (1987).
Pape, H. C. & Mager, R. Nitric oxide controls oscillatory activity in thalamocortical neurons. Neuron 9, 441–448 (1992).
Lein, E. S. et al. Genome-wide atlas of gene expression in the adult mouse brain. Nature 445, 168–176 (2007).
Zhu, M. et al. How the body controls brain temperature: the temperature shielding effect of cerebral blood flow. J. Appl. Physiol. (1985) 101, 1481–1488 (2006).
Kacem, K. et al. Structural organization of the perivascular astrocyte endfeet and their relationship with the endothelial glucose transporter: a confocal microscopy study. Glia 23, 1–10 (1998).
Garthwaite, G. et al. Signaling from blood vessels to CNS axons through nitric oxide. J. Neurosci. 26, 7730–7740 (2006).
Islas, L., Pasantes-Morales, H. & Sanchez, J. A. Characterization of stretch-activated ion channels in cultured astrocytes. Glia 8, 87–96 (1993).
Rungta, R. L. et al. Diversity of neurovascular coupling dynamics along vascular arbors in layer II/III somatosensory cortex. Commun. Biol. 4, 855 (2021).
Bauer, A. Q. et al. Effective connectivity measured using optogenetically evoked hemodynamic signals exhibits topography distinct from resting state functional connectivity in the mouse. Cereb. Cortex 28, 370–386 (2018).
Lecrux, C. & Hamel, E. Neuronal networks and mediators of cortical neurovascular coupling responses in normal and altered brain states. Philos. Trans. R. Soc. Lond. B 371, 20150350 (2016).
Cauli, B. & Hamel, E. Revisiting the role of neurons in neurovascular coupling. Front. Neuroenergetics 2, 9 (2010).
Krogsgaard, A. et al. PV interneurons evoke astrocytic Ca(2+) responses in awake mice, which contributes to neurovascular coupling. Glia 71, 1830–1846 (2023).
Liu, T. T. Noise contributions to the fMRI signal: an overview. Neuroimage 143, 141–151 (2016).
Jones, T. B., Bandettini, P. A. & Birn, R. M. Integration of motion correction and physiological noise regression in fMRI. Neuroimage 42, 582–590 (2008).
Birn, R. M. The role of physiological noise in resting-state functional connectivity. Neuroimage 62, 864–870 (2012).
Cardenas, D. P. et al. Functional MRI during hyperbaric oxygen: effects of oxygen on neurovascular coupling and BOLD fMRI signals. Neuroimage 119, 382–389 (2015).
Weckesser, M. et al. Functional imaging of the visual cortex with bold-contrast MRI: hyperventilation decreases signal response. Magn. Reson. Med. 41, 213–216 (1999).
Rostrup, E. et al. Regional differences in the CBF and BOLD responses to hypercapnia: a combined PET and fMRI study. Neuroimage 11, 87–97 (2000).
Birn, R. M. et al. Separating respiratory-variation-related fluctuations from neuronal-activity-related fluctuations in fMRI. Neuroimage 31, 1536–1548 (2006).
Shmueli, K. et al. Low-frequency fluctuations in the cardiac rate as a source of variance in the resting-state fMRI BOLD signal. Neuroimage 38, 306–320 (2007).
Raut, R. V. et al. Global waves synchronize the brain’s functional systems with fluctuating arousal. Sci. Adv. 7, eabf2709 (2021).
Hasler, F. et al. Acute psychological and physiological effects of psilocybin in healthy humans: a double-blind, placebo-controlled dose-effect study. Psychopharmacology (Berl.) 172, 145–156 (2004).
Enager, P. et al. Pathway-specific variations in neurovascular and neurometabolic coupling in rat primary somatosensory cortex. J. Cereb. Blood Flow Metab. 29, 976–986 (2009).
Li, B. & Freeman, R. D. Neurometabolic coupling between neural activity, glucose, and lactate in activated visual cortex. J. Neurochem. 135, 742–754 (2015).
Gorrostieta, C. et al. Investigating brain connectivity using mixed effects vector autoregressive models. Neuroimage 59, 3347–3355 (2012).
Casey, A. B. et al. ‘Selective’ serotonin 5-HT2A receptor antagonists. Biochem. Pharmacol. 200, 115028 (2022).
Jaster, A. M. et al. Differences across sexes on head-twitch behavior and 5-HT(2A) receptor signaling in C57BL/6J mice. Neurosci. Lett. 788, 136836 (2022).
Turner, K. L., Gheres, K. W. & Drew, P. J. Relating pupil diameter and blinking to cortical activity and hemodynamics across arousal states. J. Neurosci. 43, 949–964 (2023).
Xu, Y. X. et al. Restraint stress induced anxiety and sleep in mice. Front. Psychiatry 14, 1090420 (2023).
Eggebrecht, A. T. et al. A quantitative spatial comparison of high-density diffuse optical tomography and fMRI cortical mapping. Neuroimage 61, 1120–1128 (2012).
Gordon, E. M. et al. Generation and evaluation of a cortical area parcellation from resting-state correlations. Cereb. Cortex 26, 288–303 (2016).
Henson, R. & Friston, K. in Statistical Parametric Mapping: the Analysis of Functional Brain Images (eds Friston, K. et al.) 178–192 (Academic, 2007).
Dana, H. et al. Sensitive red protein calcium indicators for imaging neural activity. eLife 5, e12727 (2016).
Park, K. et al. Optical imaging of stimulation-evoked cortical activity using GCaMP6f and jRGECO1a. Quant. Imaging Med. Surg. 11, 998–1009 (2021).
Dana, H. et al. Thy1 transgenic mice expressing the red fluorescent calcium indicator jRGECO1a for neuronal population imaging in vivo. PLoS ONE 13, e0205444 (2018).
Chen, Q. et al. Imaging neural activity using Thy1-GCaMP transgenic mice. Neuron 76, 297–308 (2012).
Padawer-Curry, J. A. et al. in Neural Repair: Methods and Protocols (eds Karamyan, V. T. & Stowe, A. M.) 113–151 (Springer, 2023).
Bice, A. R. et al. Homotopic contralesional excitation suppresses spontaneous circuit repair and global network reconnections following ischemic stroke. eLife 11, e68852 (2022).
Cao, D. et al. Structure-based discovery of nonhallucinogenic psychedelic analogs. Science 375, 403–411 (2022).
Qu, Y. et al. Rapid antidepressant-like effect of non-hallucinogenic psychedelic analog lisuride, but not hallucinogenic psychedelic DOI, in lipopolysaccharide-treated mice. Pharm. Biochem. Behav. 222, 173500 (2023).
de la Fuente Revenga, M. et al. Fully automated head-twitch detection system for the study of 5-HT2A receptor pharmacology in vivo. Sci. Rep. 9, 14247 (2019).
Halberstadt, A. L. et al. Use of the head-twitch response to investigate the structure-activity relationships of 4-thio-substituted 2,5-dimethoxyphenylalkylamines. Psychopharmacology (Berl.) 240, 115–126 (2023).
Halberstadt, A. L. & Geyer, M. A. Characterization of the head-twitch response induced by hallucinogens in mice: detection of the behavior based on the dynamics of head movement. Psychopharmacology (Berl.) 227, 727–739 (2013).
Moliner, R. et al. Psychedelics promote plasticity by directly binding to BDNF receptor TrkB. Nat. Neurosci. 26, 1032–1041 (2023).
Shahar, O. et al. Role of 5-HT2A, 5-HT2C, 5-HT1A and TAAR1 receptors in the head twitch response induced by 5-hydroxytryptophan and psilocybin: translational implications. Int. J. Mol. Sci. 23, 14148 (2022).
Halberstadt, A. L. et al. Correlation between the potency of hallucinogens in the mouse head-twitch response assay and their behavioral and subjective effects in other species. Neuropharmacology 167, 107933 (2020).
Canal, C. E. & Morgan, D. Head-twitch response in rodents induced by the hallucinogen 2,5-dimethoxy-4-iodoamphetamine: a comprehensive history, a re-evaluation of mechanisms, and its utility as a model. Drug Test. Anal. 4, 556–576 (2012).
Sun, Y. et al. An ontogenic study of receptor mechanisms by which acute administration of low-doses of methamphetamine suppresses DOI-induced 5-HT(2A)-receptor mediated head-twitch response in mice. BMC Neurosci. 23, 2 (2022).
Kuo, C. C. et al. Inhibitory interneurons regulate phasic activity of noradrenergic neurons in the mouse locus coeruleus and functional implications. J. Physiol. 598, 4003–4029 (2020).
Kuo, C. C. & McCall, J. G. Neural circuit-selective, multiplexed pharmacological targeting of prefrontal cortex-projecting locus coeruleus neurons drives antinociception. Preprint at bioRxiv https://doi.org/10.1101/2024.06.08.598059 (2024).
Seitzman, B. A. et al. Functional network disorganization and cognitive decline following fractionated whole-brain radiation in mice. Geroscience 46, 543–562 (2024).
Bianciardi, M. et al. Sources of functional magnetic resonance imaging signal fluctuations in the human brain at rest: a 7 T study. Magn. Reson. Imaging 27, 1019–1029 (2009).
Snyder, A. Z. in Neuroscience in the 21st Century: From Basic to Clinical 2nd ed. (eds Pfaff, D. W. & Volkow, N. D.) 1939–1990 (Springer, 2022).
Chen, S. et al. Visual deprivation during mouse critical period reorganizes network-level functional connectivity. J. Neurosci. 44, e1019232024 (2024).
Barron, J. L. et al. Performance of optical flow techniques. In Proc. 1992 IEEE Computer Society Conference on Computer Vision and Pattern Recognition. 236–242 (IEEE Computer Society, 1992).
Barron, J. L., Fleet, D. J. & Beauchemin, S. S. Performance of optical flow techniques. Int. J. Computer Vis. 12, 43–77 (1994).
Reimer, J. et al. Pupil fluctuations track fast switching of cortical states during quiet wakefulness. Neuron 84, 355–362 (2014).
Raut, R. V. et al. Arousal as a universal embedding for spatiotemporal brain dynamics. Nature https://doi.org/10.1038/s41586-025-09544-4 (2025).
Syeda, A. et al. Facemap: a framework for modeling neural activity based on orofacial tracking. Nat. Neurosci. 27, 187–195 (2024).
Kirillov, A. et al. Segment anything. Preprint at https://doi.org/10.48550/arXiv.2304.02643 (2023).
Dosovitskiy, A. et al. An image is worth 16x16 words: transformers for image recognition at scale. Preprint at https://doi.org/10.48550/arXiv.2010.11929 (2020).
White, B. R. et al. Imaging of functional connectivity in the mouse brain. PLoS ONE 6, e16322 (2011).
Wright, P. W. et al. Functional connectivity structure of cortical calcium dynamics in anesthetized and awake mice. PLoS ONE 12, e0185759 (2017).
Prahl, S. A. Tabulated molar extinction coefficient for hemoglobin in water. http://omlc.org/spectra/hemoglobin/summary.html (1999).
Arridge, S. R., Cope, M. & Delpy, D. T. The theoretical basis for the determination of optical pathlengths in tissue: temporal and frequency analysis. Phys. Med. Biol. 37, 1531–1560 (1992).
Ma, Y. et al. Resting-state hemodynamics are spatiotemporally coupled to synchronized and symmetric neural activity in excitatory neurons. Proc. Natl Acad. Sci. USA 113, E8463–E8471 (2016).
Fox, M. D. et al. The global signal and observed anticorrelated resting state brain networks. J. Neurophysiol. 101, 3270–3283 (2009).
Paxinos, G. & Franklin, K. B. Paxinos and Franklin’s the Mouse Brain in Stereotaxic Coordinates (Academic, 2016).
Kraft, A. W. et al. Sensory deprivation after focal ischemia in mice accelerates brain remapping and improves functional recovery through Arc-dependent synaptic plasticity. Sci. Transl. Med. 10, eaag1328 (2018).
Murphy, K. & Fox, M. D. Towards a consensus regarding global signal regression for resting state functional connectivity MRI. Neuroimage 154, 169–173 (2017).
Murphy, K. et al. The impact of global signal regression on resting state correlations: are anti-correlated networks introduced? Neuroimage 44, 893–905 (2009).
Thirion, B. et al. Which fMRI clustering gives good brain parcellations?. Front. Neurosci. 8, 167 (2014).
Ward, J. H. Hierarchical grouping to optimize an objective function. J. Am. Stat. Assoc. 58, 236–244 (1963).
Park, H. et al. Agreement between functional connectivity and cortical thickness-driven correlation maps of the medial frontal cortex. PLoS ONE 12, e0171803 (2017).
Bowman, F. D., Patel, R. & Lu, C. Methods for detecting functional classifications in neuroimaging data. Hum. Brain Mapp. 23, 109–119 (2004).
Lake, E. M. R. et al. Simultaneous cortex-wide fluorescence Ca2+ imaging and whole-brain fMRI. Nat. Methods 17, 1262–1271 (2020).
Petri, G. et al. Homological scaffolds of brain functional networks. J. R. Soc. Interface 11, 20140873 (2014).
Brier, L. M. & Culver, J. P. Open-source statistical and data processing tools for wide-field optical imaging data in mice. Neurophotonics 10, 016601 (2023).
White, B. R. et al. Controlling the familywise error rate in widefield optical neuroimaging of functional connectivity in mice. Neurophotonics 10, 015004 (2023).
Bauer, A. WFOI data—calcium and hemodynamic image sequences, pre/post DOI. Zenodo https://doi.org/10.5281/zenodo.15857640 (2025).
Bauer, A. Head twitch responses, cohort 1. Zenodo https://doi.org/10.5281/zenodo.15857232 (2025).
Acknowledgements
We thank R. Raut, T. Laumann and R. Reneau (Washington University in St. Louis) for helpful discussions and advice. We would also like to express our sincere gratitude to the mice for their vital contributions to this study, which were essential for the findings presented in this paper. We also invite readers to view the bioRxiv version of this manuscript, where the absence of citation limits allowed us to more fully acknowledge the many additional authors, laboratories and scientific contributions that helped shape this work. Finally, we would like to acknowledge the Osage Nation, Missouria, Illinois Confederacy and many other tribes as the ancestral, traditional and contemporary custodians of the land where this work was conducted. This work was supported by National Institutes of Health grants R01NS126326 (A.Q.B.), R01NS102870 (A.Q.B.), RF1AG07950301 (A.Q.B.), R01NS117899 (J.G.M.), R01NS135401 (J.G.M.), F99NS139512 (J.A.P.-C.), T32EB014855 (J.A.P.-C.), T32NS121881 (O.J.K.). This work was also supported by the Center for Holistic Interdisciplinary Research in Psychedelics (CHIRP), a Here and Next Transcend Initiative funded by Washington University in St. Louis.
Author information
Authors and Affiliations
Contributions
J.A.P.-C., J.S.S., J.G.M. and A.Q.B. developed project conceptualization and experimental design. J.A.P.-C., O.J.K., C.-C.K., X.W., J.G.M. and A.Q.B. provided technical resources. J.A.P.-C., O.J.K, C.-C.K. and A.R.B. managed animal preparation and surgeries. J.A.P.-C., O.J.K. and C.-C.K. acquired data. J.A.P.-C., O.J.K., C.-C.K., J.S.S., J.G.M. and A.Q.B. analyzed the data. J.S.S., J.G.M. and A.Q.B. supervised the project. J.A.P.-C, A.Z.S., J.S.S., J.G.M. and A.Q.B. wrote the original manuscript. J.A.P.-C, O.J.K, C.-C.K., X.W., A.R.B., G.E.N., A.Z.S., J.S.S., J.G.M. and A.Q.B. reviewed and approved the final manuscript. J.A.P.-C., G.E.N., J.S.S., J.G.M. and A.Q.B. provided financial support.
Corresponding authors
Ethics declarations
Competing interests
J.S.S. has received consulting fees from Forbes Manhattan. G.E.N. has received research support from Usona Institute (drug only). In the past 36 months, G.E.N. has received salary support from institutional grants supported by the National Center for Translational Sciences (NCATS) and the National Institute of Digestive and Diabetes and Kidney Diseases (NIDDK); research support from the National Institute of Mental Health (NIMH), the Health Resources and Services Administration (HRSA), the Taylor Family Institute for Innovative Psychiatric Research and the Center for Holistic Interdisciplinary Research on Psychedelics (CHIRP) at Washington University. She has served as a Coinvestigator or Principal Investigator for studies funded by COMPASS Pathways, LB Pharmaceuticals, Inc. and Usona Institute. These potential conflicts of interest have been reviewed and are managed by Washington University School of Medicine. The other authors declare no competing interests.
Peer review
Peer review information
Nature Neuroscience thanks Boris Heifets and the other, anonymous, reviewer(s) for their contribution to the peer review of this work.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Extended data
Extended Data Fig. 1 Psychedelics alter hemodynamic response functions (HRFs) in humans as reported by BOLD-fMRI.
a) Previously published functional magnetic resonance imaging (fMRI) data were analyzed from humans undergoing a simple auditory-visual matching task during acute exposure to psilocybin, methylphenidate, or no compound. Generalized linear models were computed for task-evoked activity in left and right visual cortex, left hand, left and right auditory cortex, and left language. HRFs were fit using double Gamma basis functions and characterized by three-parameters: Peak Value (P), Dispersion (D), and Time to Peak (T). Data are presented as mean ± SEM. b) Quantified changes in human hemodynamic response functions. Statistical differences for each parameter across conditions were assessed via 2-way ANOVA; significant effects were evaluated post-hoc via two-sided t-tests. p-values are indicated such that #<0.05, ##<0.01, ###<0.001. Data are presented as quantiles and all outliers (>99.3rd percentile) were excluded for visualization purposes only. Top panel: Peak value decreased in left and right visual regions under psilocybin conditions (ANOVA, p=0.030 and p=0.020, respectively). Middle panel: Time-to-peak decreased in all regions except the right visual region (ANOVA, left visual, p=0.03 left hand, p=0.006; left auditory, p=0.002; right auditory, p=0.044; left language, p=0.007). Bottom panel: Dispersion decreased in the left and right visual regions, and left and right auditory regions (ANOVA, p=0.007, p=0.004, p=0.021, p=0.022, respectively).
Extended Data Fig. 2 Mouse body movements and pupil dynamics do not differ across experimental conditions.
Two cameras were placed in front of the mouse to monitor body movement and pupil dynamics. Image sequences of mouse movement and pupil dynamics were time locked to WFOI. a) Motion monitoring using optical flow. (i) Optical flow (OF) estimates were generated using the Lucas-Kanade method (ii); OF yields a spatiotemporal series of vectors, with each pixel in each frame assigned a vector that reflects local spatiotemporal gradients. At each pixel, the vector magnitude provides a scalar measure of motion, with higher values corresponding to larger movements (iii). The temporal variance of vector magnitudes is shown, with yellow colors indicating larger variance, and demonstrates OF sensitivity to the motion of the felt pouch supporting the mouse (which moves when the mouse moves). (iv) The vector magnitude was computed at each pixel; higher values correspond to larger movement. (v) The root-sum square across all frames was computed to yield a single time trace of estimated movement. b) Distributions of OF-derived movement measures. No differences were observed between conditions (\(n=8\) mice per condition, as assessed by differences between post- minus pre- injection; Kruskal-Wallis: p=0.262). c) Simplified schematic of the automated pupil segmentation algorithm. d) Example of a user-annotated pupil segmentation and the automated pupil segmentation using our algorithm. As demonstrated, user-guided segmentation matches well with automated segmentation. e) Automated outlier detection. Large spatial deviations in the likelihood distributions of the pupil mask’s center-of-mass are attributed to artifact due to the slow rate of pupil diameter fluctuations compared to the imaging frame rate. Outlier frames, along with the frames immediately before and after, were manually labeled and appended to the original training dataset to increase robustness of the algorithm. f) Example frames and time course of pupil area change during imaging. g) Changes in pupil area before and after compound injection. Distributions of pupil area changes did not differ before and after injection of any compound (saline, DOI (4mg/kg), MDL (0.1mg/kg), DOI+MDL, \(n=8\) mice per condition).
Extended Data Fig. 3 5-HT2AR-mediated action potential-independent increase in global calcium signal.
Ex vivo recordings of global calcium signals in acute brain slices upon DOI (10μM) and TTX (1μM) application (\(n=6\) mice). a) Experimental arrangement. Global imaging and cell-attached recordings were conducted on prefrontal cortical region (for example, cingulate, prelimbic and infralimbic cortex) and pyramidal cells within these regions, respectively. b) Subthreshold effects led by direct activation of 5-HT2AR. Representative cell-attached recordings demonstrate lack of spontaneous activity of cortical pyramidal cells in acute brain slice preparation. Application of NMDA, but not DOI or DOI+TTX, resulted in burst firing in recorded cells (aCSF only: 3 cells, NMDA: 5 cells, DOI: 6 cells, DOI upon TTX: 4 cells). c) Tonic, TTX-resistant global calcium increases upon DOI application. Representative traces corrected for photobleaching show a small but steady increase of global calcium signal led by DOI application, which was not affected by TTX pretreatment. (1.55±0.96 vs. 0.87±0.78%, two-way, t-test: p=0.314). d) Elimination of 5-HT2AR-mediated action potential-independent tonic calcium signal by temporal filtering. Representative filtered traces from the same recordings shown in panel c (0.01-5Hz) demonstrate complete removal of action potential-independent DOI-mediated calcium increase (aCSF: 0.00±0.01 vs. 0.01±0.03%, two-sided, Wilcoxon rank-sum test: p=0.343). Abbreviations: aca, anterior commissure (anterior part); AI, anterior insular cortex; Cg1, cingulate cortex (area 1); IL, infralimbic cortex; M1, primary motor cortex; M2, secondary motor cortex; PrL, prelimbic cortex.
Extended Data Fig. 4 5-HT2AR activation does not contribute to calcium fluctuations in electrically-evoked neural activity.
Ex vivo recordings of electrical stimulation-evoked calcium signal in acute brain slice upon DOI (10μM) application (\(n=4\) mice). a-b) Experimental arrangement. Simultaneous cell-attached recordings and calcium imaging was conducted on layer V pyramidal cells in prefrontal cortical region (for example, cingulate, prelimbic and infralimbic cortices) with a stimulation electrode placed in layer II/III. c) Direct activation of 5-HT2AR does not affect the number of evoked action potentials. (top) Representative cell-attached recordings demonstrate evoked action potentials in both baseline and DOI administration. (bottom) Summarized plot showing the input-output relationship between stimulation intensity and the number of evoked action potentials. Note that all representative traces in c-e are from recordings with 50% of stimulation intensity. d) Activation of 5-HT2AR does not affect the evoked calcium fluctuation in recorded cell. (top) Representative stimulus-evoked traces (solid: mean; pale: individual) demonstrate evoked calcium fluctuation in recorded cells shown in b. (bottom) Summarized plot showing the input-output relationship between stimulation intensity and evoked calcium fluctuation. e) Activation of 5-HT2AR does not affect the evoked calcium fluctuation in recorded slices. (top) Representative peri-stimulus traces (solid: mean; pale: individual) demonstrate evoked calcium fluctuation in the whole FOV. (bottom) Summarized plot showing the input-output relationship between stimulation intensity and evoked calcium fluctuation. f) Activation of 5-HT2AR does not affect the evoked calcium fluctuation under blockade of synaptic transmission. (upper left) Global calcium imaging recording was conducted on layer V pyramidal cells in prefrontal cortices with stimulation electrode locally placed in layer V under kynurenic acid (5mM) and picrotoxin (100μM). (lower left) Representative peri-stimulus traces (solid: mean; pale: individual) demonstrate evoked global calcium fluctuation. (right) No significant difference in evoked calcium fluctuation. Abbreviations: aca, anterior commissure (anterior part); AI, anterior insular cortex; AP, action potential; Cg1, cingulate cortex (area 1); IL, infralimbic cortex; M1, primary motor cortex; M2, secondary motor cortex; PrL, prelimbic cortex.
Extended Data Fig. 5 Hallucinogenic 5-HT2A receptor agonism alters resting-state activity.
Global (that is, regionally averaged) PSDEs were integrated over double-octave bins and compared between compounds (saline, DOI (4mg/kg), MDL (0.1mg/kg), DOI+MDL, \(n=8\)) for both calcium (a) and hemodynamic (b) activity. The lowest frequency bin was limited by the lower limit of 0.02Hz due to the window length used to segment resting-state from stimulus-evoked data. All measures were averaged over the cortex and comparisons were made between pre- and post-injection (post- minus pre-) and across compounds using two-sided t-tests (*<0.05, **<0.01, and ***<0.001) and one-way ANOVAs with post-hoc, two-sided, t-tests (#<0.05, ##<0.01, ###<0.001), respectively. Multiple comparisons were corrected for using the Bonferroni method. c) Fractional changes in broad-band activity attributed to hallucinogenic 5-HT2A receptor agonism. Fractional changes in regional power spectra are displayed as the ratio of post-injection values divided by pre-injection values and regions are organized from anterior-to-posterior and from left-to-right. DOI differentially affected the spectral content of calcium (left) and hemodynamic (right) activity in a region-dependent manner. ISA: DOI increased calcium ISA power in frontal, cingulate, and motor cortices while hemodynamic ISA activity increased primarily in somatosensory regions. Intermediate: Calcium activity decreased in all cortical regions excluding frontal and cingulate. In contrast, hemodynamic activity decreased in frontal and cingulate cortex. Delta: Both calcium and hemodynamic activity exhibited increased delta band activity in nearly every cortical region examined with the largest increases occurring at ~0.8Hz hemodynamic activity. Across all frequencies examined, MDL largely reversed the effects of DOI. Saline and MDL minimally altered hemodynamics and calcium spectral content. Average fractional changes over the cortex are visualized on top of each column and plotted as mean ± std across mice.
Extended Data Fig. 6 Hallucinogenic 5-HT2A receptor agonism alters resting-state neurovascular coupling.
a) Neurovascular coupling parameterization and quantification. Regional hemodynamic response functions (HRFs) were characterized by their peak value (peak), time to peak amplitude (TTP), and full width at half maximum (FWHM). These parameters were averaged across regions and compared across compounds (saline, DOI (4mg/kg), MDL (0.1mg/kg), DOI+MDL, \(n=8\)). Predicted hemodynamic activity was calculated by convolving each region’s solved HRF with its corresponding measured calcium signal. Model fit was assessed via Pearson’s r between predicted and measured hemodynamic activity. Band-limited frequency transduction was calculated over three frequency bands: infraslow activity band (ISA, 0.03-0.08Hz), intermediate activity band (0.08-0.50Hz), and delta activity band (0.5-4.0Hz). The lowest frequency bin of the TF (that is, 0.03Hz) was determined by the 30s long window used for estimating HRFs. All global measures were averaged over the cortex and comparisons were made between pre- and post-injection and across compounds (post- minus pre-) using two-tailed t-tests (*<0.05, **<0.01, and ***<0.001) and one-way ANOVAs evaluated post-hoc via two-tailed t-test (#<0.05, ##<0.01, ###<0.001), respectively. Multiple comparisons were corrected for using the Bonferroni method. b) Hallucinogenic 5-HT2A receptor agonism alters region-specific NVC. Estimated HRFs were averaged across mice for each region and compound. Regions are organized from anterior-to-posterior and from left-to-right. Prior to compound injection, NVC differed (parameterized in Fig. 5c) across the cortex. For instance, somatosensory and parietal regions exhibited the strongest coupling (for example, large peak value), while frontal and cingulate regions exhibited more modest coupling. DOI caused a zero-lag notch across the entire cortex (that is, the acausal feature seen in Fig. 5a, black arrow), strongest motor, retrosplenial, and visual regions. Additionally, DOI increased coupling in the somatosensory and parietal regions. These effects were dampened when DOI was administered with MDL. Cortical-averaged HRFs are visualized above each column and plotted as mean +/- std across mice. The color axis has units of \(\Delta \mu\)Mol/\(\Delta\)F/F(%).
Extended Data Fig. 7 Lagged cross-covariance and coherence between calcium and total hemoglobin recapitulate DOI-induced changes in neurovascular coupling.
a) Global hemodynamic response functions: The global HRF estimated using deconvolution following DOI injection (4mg/kg, \(n=8\), top; Fig. 4a) and the global lagged cross-covariance function (CCF; bottom) between resting-state calcium and hemodynamics. Pre-injection CCFs display a simple, quasi-causal relation (calcium leads hemodynamics). DOI decreased coupling (peak value, -51% [-66%, -13%], pre- vs. post-injection: p=0.039; Kruskal-Wallis: p=0.004; saline vs. DOI: p=0.020; DOI vs. DOI+MDL: p=0.004; DOI vs. MDL (0.1mg/kg): p<0.001). DOI induced a zero-lag notch in the CCF (black arrow). Covariance has units of \(\Delta \mu\)Mol\(\Delta\)F/F(%) and the HRF has units of \(\Delta \mu\)Mol/\(\Delta\)F/F(%). b) DOI-induced alterations in CCF over the cortex. Regional distribution of cross covariance after the injection of DOI reveals a negative-lag peak across the cortex. CCFs for all regions and compounds are displayed in Supplementary Fig. S6a. c) Magnitude Squared Coherence (MSC) between calcium and hemodynamic activity. Before compound injection, MSC exhibited a large 0.2Hz peak, demonstrating the band-limited nature of NVC. After injection of DOI, a coherence peak ~0.8Hz (~0.5Hz half-bandwidth) emerged (0.5-2.0Hz: +170% [110% 230%]; pre- vs. post-injection: p=0.008; Kruskal-Wallis: p=0.004; saline vs. DOI: p=0.009; DOI vs. DOI+MDL: p=0.030; DOI vs. MDL: p=0.025). Moreover, the 0.2Hz peak present before injection largely diminished (<0.5Hz: -47% [-53%, 49%], pre- vs. post-injection: p=0.008). This phenomenon signifies a DOI-induced shift in the coherent frequencies contained in both neuronal and hemodynamic activity. DOI+MDL largely reversed the effects of DOI alone (Supplementary Fig. S6). Coherence for all regions and compounds are displayed in Supplementary Fig. S6b. d) Phase relation between calcium and hemodynamics. Prior to injection, global calcium activity led hemodynamic activity (positive slope at frequencies below ~1Hz; left). After DOI injection, phase relations between calcium and hemodynamic activity flipped sign (for example, from positive to negative slope in primary motor and primary somatosensory hindpaw; right), indicating that hemodynamics precede calcium over these frequencies. Supplementary Fig. S6c reports phase relations for all regions and compounds. All data are presented mean ± std across mice. Significance was determined via Kruskal-Wallis tests and post-hoc, two-sided Wilcoxon’s sign-rank corrected post-hoc using Bonferroni correction.
Extended Data Fig. 8 Non-hallucinogenic doses of DOI and the non-hallucinogenic, 5-HT2AR agonist, Lisuride do not alter NVC.
a) Head twitch responses (HTRs) recorded during sub-hallucinogenic DOI (low-dose, 0.04mg/kg) and the non-hallucinogenic, psychedelic ligand, Lisuride (0.1mg/kg). No differences in HTRs were observed between groups (Kruskal-Wallis, p>0.99). HTRs were recorded for 30 minutes (\(n=4\)). Significance was determined via Kruskal-Wallis tests and evaluated post-hoc via two-sided Wilcoxon’s sign-rank test corrected for multiple comparisons using Bonferroni correction. a-f) WFOI imaging of non-hallucinogenic ligands. Two groups of Thy1-jRGECO1a mice were imaged under resting-state conditions for 30min before and 30min after injection. Cohort 1 (4M, 5F) received saline and lisuride (0.1mg/kg); cohort 2 (2M, 2F) received saline and a sub-hallucinogenic dose of DOI (low-dose DOI, 0.04mg/kg. b) Hemodynamic response functions before and after injection of low-dose DOI. Low-dose DOI did not alter global estimates of NVC. c) Regional changes NVC following low-dose DOI. Full-width-at-half-maximum (FWHM), time-to-peak (TTP), and peak value are visualized across the cortex. Sub-hallucinogenic doses of DOI did not alter any regional HRF parameters compared to saline. d) Hemodynamic response functions before and after Lisuride injection. No significant changes to global hemodynamic response functions were observed after the injection of lisuride. e) Regional changes NVC following Lisuride. Parameterized hemodynamic response functions reveal no differences in regional NVC after injection of low-dose DOI. f) Global power spectral density estimates before and after injection of lisuride and d) Integrated, band-limited power over double octave frequency bins. Spectra are reported as median and shaded 25th and 75th percentiles. Band-limited power is consistently decreased after injection of both saline and lisuride, like results reported in Extended Data Fig. 5a,b. Lisuride modestly affected power over the 0.32-1.28Hz bin (p=0.009). Band-limited power was tested for significance via Kruskal-Wallis tests and evaluated post-hoc via two-tailed Wilcoxon’s sign-rank corrected for multiple comparisons using Bonferroni correction.
Extended Data Fig. 9 Hallucinogenic 5-HT2A receptor agonism alters inter- vs. intra-network connectivity.
a) Whole-cortex RSFC matrices. Average, pre-injection, whole-cortex RSFC matrices across mice (\(n=8\)) for infralow (ISA; 0.02-0.08Hz) calcium and HbT activity and delta (0.5-4Hz) calcium activity. Matrices represent all pixel-pair comparisons and are organized by network assignment. b) Hierarchical clustering of RSFC. RSFC matrices were hierarchically clustered using Ward’s method and visualized as a dendrogram. c) Resultant clusters at a specific hierarchical level. A distance threshold was selected to create a parcellation containing 13 clusters. These parcels were used to assess inter- and intra-regional relationships. d) Pre-injection Silhouette scores. Silhouette scores were computed at each pixel to assess the ratio between cohesion (intra-cluster distance metric) and separation (inter-cluster metric). Silhouette scores range from -1 to +1, such that +1 indicates a point is well-clustered and far from neighboring clusters (strongly intra-connected, weakly inter-connected), 0 indicates that the point is on, or near, the boundary between two clusters (equally intra- and inter-connected/equally likely to belong to any neighboring cluster), and -1 indicates a point is misclassified and likely should be assigned to a different cluster than its own (weakly intra-connected, strongly inter-connected). i) A graphical demonstration of cohesion (black cluster) and separation to the nearest cluster (red cluster). ii) Pre-injection silhouette scores for each species and frequency band. e) Changes in Silhouette scores after compound injection. DOI (4mg/kg) predominantly decreases ISA calcium silhouette scores everywhere (regions are less intra-connected and more inter-connected), an effect not fully reversed in under DOI+MDL. Importantly, ISA calcium scores in the cingulate and retrosplenial cortices—constellations of the mouse default mode network—showed opposite effects: cingulate scores marginally increased, while retrosplenial scores markedly decreased. This divergence was substantially reduced when scores were computed using infralow hemodynamic activity. Delta calcium activity under DOI suggests that the compound primarily alters network boundaries rather than modulating inter- or intra-network interactions.
Supplementary information
Supplementary Information
Supplementary Results, Discussion, Table 1, Figs. 1–6 and References.
Rights and permissions
Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.
About this article
Cite this article
Padawer-Curry, J.A., Krentzman, O.J., Kuo, CC. et al. Psychedelic 5-HT2A receptor agonism alters neurovascular coupling and differentially affects neuronal and hemodynamic measures of brain function. Nat Neurosci 28, 2330–2343 (2025). https://doi.org/10.1038/s41593-025-02069-z
Received:
Accepted:
Published:
Version of record:
Issue date:
DOI: https://doi.org/10.1038/s41593-025-02069-z