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.

Ultra-slow Oscillations in fMRI and Resting-State Connectivity: Neuronal and Vascular Contributions and Technical Confounds

Global Fluctuations of Cerebral Blood Flow Indicate a Global Brain Network Independent of Systemic Factors.

Associations of Resting-State fMRI Functional Connectivity with Flow-BOLD Coupling and Regional Vasculature

Low Frequency Systemic Hemodynamic “Noise” in Resting State BOLD fMRI: Characteristics, Causes, Implications, Mitigation Strategies, and Applications

Resting-state "physiological networks".

Arousal Contributions to Resting-State Fmri Connectivity and Dynamics

An Orderly Sequence of Autonomic and Neural Events at Transient Arousal Changes

Linking brain vascular physiology to hemodynamic response in ultra-high field MRI

Low-frequency calcium oscillations accompany deoxyhemoglobin oscillations in rat somatosensory cortex

Resting-state hemodynamics are spatiotemporally coupled to synchronized and symmetric neural activity in excitatory neurons

Ongoing physiological processes in the cerebral cortex

Brain-wide functional connectivity artifactually inflates throughout functional magnetic resonance imaging scans

Ultra-slow fMRI fluctuations in the fourth ventricle as a marker of drowsiness

Functional connectivity is dominated by aperiodic, rather than oscillatory, coupling

On the Stability of BOLD fMRI Correlations

Neurovascular Impulse Response Function (IRF) during spontaneous activity differentially reflects intrinsic neuromodulation across cortical regions

Brain-wide Neural Co-Activations in Resting Human

The Oscillating Brain: Complex and Reliable.

Detecting Perfusion Pattern Based on the Background Low-Frequency Fluctuation in Resting-State Functional Magnetic Resonance Imaging Data and Its Influence on Resting-State Networks: An Iterative Postprocessing Approach

Coherent Slow Cortical Potentials Reveal a Superior Localization of Resting-State Functional Connectivity Using Voltage-Sensitive Dye Imaging

Synchronous Multiscale Neuroimaging Environment for Critically Sampled Physiological Analysis of Brain Function: Hepta-Scan Concept.