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.

Neurovascular coupling: insights from multi-modal dynamic causal modelling of fMRI and MEG

Bayesian fusion and multimodal DCM for EEG and fMRI

Physiologically informed dynamic causal modeling of fMRI data

Investigating Human Neurovascular Coupling Using Functional Neuroimaging: A Critical Review of Dynamic Models

Determining Excitatory and Inhibitory Neuronal Activity from Multimodal fMRI Data Using a Generative Hemodynamic Model

Neurovascular Coupling Analysis Based on Multivariate Variational Gaussian Process Convergent Cross-Mapping

Nonlinear Local Electrovascular Coupling. I: A Theoretical Model.

Quantifying neurovascular coupling through a concurrent assessment of arterial, capillary, and neuronal activation in humans: A multimodal EEG-fNIRS-TCD investigation

The Neural Basis of the Hemodynamic Response Nonlinearity in Human Primary Visual Cortex: Implications for Neurovascular Coupling Mechanism.

[Neurovascular Coupling Analysis of Working Memory Based on Electroencephalography and Functional Near-Infrared Spectroscopy].

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

Network response of brain microvasculature to neuronal stimulation

Neural Coupling Mechanism in Fmri Hemodynamics

Bayesian Modelling Approaches for Breath-Hold Induced Cerebrovascular Reactivity

A model of neurovascular coupling and its application to cortical spreading depolarization.

Fusing Eeg And Fmri Based On A Bottom-Up Model: Inferring Activation And Effective Connectivity In Neural Masses

Quantitatively Interpreting fMRI signal

Neurophysiological investigation of the basis of the fMRI signal

Transmodal Analysis of Neural Signals

Human electromagnetic and haemodynamic networks systematically converge in unimodal cortex and diverge in transmodal cortex

Two distinct modes of hemodynamic responses across the human cortex