Neurovascular Impulse Response Function (IRF) during spontaneous activity differentially reflects intrinsic neuromodulation across cortical regions
Bradley C. Rauscher,Natalie Fomin-Thunemann,Sreekanth Kura,Patrick R. Doran,Pablo D. Perez,Kıvılcım Kılıç,Emily A. Martin,Dora Balog,Nathan X. Chai,Francesca A. Froio,Patrick F. Bloniasz,Kate E. Herrema,Rockwell Tang,Scott G. Knudstrup,Andrew Garcia,John X. Jiang,Jeffrey P. Gavornik,David Kleinfeld,Michael E. Hasselmo,Laura D. Lewis,Sava Sakadzic,Lei Tian,Gal Mishne,Emily P. Stephen,Martin Thunemann,David A. Boas,Anna Devor
DOI: https://doi.org/10.1101/2024.09.14.612514
2024-09-15
Abstract:Ascending neuromodulatory projections from deep brain nuclei generate internal brain states that differentially engage specific neuronal cell types. Because neurovascular coupling is cell-type specific and neuromodulatory transmitters have vasoactive properties, we hypothesized that the impulse response function (IRF) linking spontaneous neuronal activity with hemodynamics would depend on neuromodulation. To test this hypothesis, we used optical imaging to measure (1) release of neuromodulatory transmitters norepinephrine (NE) or acetylcholine (ACh), (2) Ca activity of local cortical neurons, and (3) changes in hemoglobin concentration and oxygenation across the dorsal surface of cerebral cortex during spontaneous neuronal activity in awake mice. A canonical convolution model with a global, stationary IRF (i.e., the convolution kernel) describing evolution of total hemoglobin (HbT, reflective of dilation dynamics) with respect to Ca , resulted in a poor fit to the data. However, the HbT time-course was well predicted, pixel-by-pixel, by a weighted sum of Ca and NE time-courses. Consistent with this result, modeling HbT as a weighted sum of stationary Ca - and NE-specific IRFs (IRF and IRF ) convolved with the respective time-courses dramatically improved the fit compared to the global IRF. IRF and IRF , estimated from the data, were positive and negative, respectively. In contrast to NE, ACh was largely redundant with Ca and therefore did not improve HbT estimation. Because NE covaried with arousal, we observed instances of the diminished hemodynamic coherence between cortical regions during high arousal despite coherent behavior of the underlying neuronal Ca activity. We conclude that while neurovascular coupling with respect to neuronal Ca is a dynamic and seemingly complex phenomenon, hemodynamic fluctuations can be captured by a simple linear model with stationary IRFs with respect to the underlying dilatory and constrictive forces. In the current study, these forces were captured by the positive IRF (dilation) and negative IRF (constriction). Without accounting for NE neuromodulation and the associated vasoconstriction, diminished hemodynamic coherence, commonly referred to as “functional (dys)connectivity” in BOLD fMRI studies, can be falsely interpreted as neuronal desynchronizations.
Neuroscience