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

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) Ca2+ 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 Ca2+, resulted in a poor fit to the data. However, the HbT time-course was well predicted, pixel-by-pixel, by a weighted sum of Ca2+ and NE time-courses. Consistent with this result, modeling HbT as a weighted sum of stationary Ca2+ - and NE-specific IRFs (IRFCa2+ and IRFNE) convolved with the respective time-courses dramatically improved the fit compared to the global IRF. IRFCa2+ and IRFNE, estimated from the data, were positive and negative, respectively. In contrast to NE, ACh was largely redundant with Ca2+ 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 Ca2+ activity. We conclude that while neurovascular coupling with respect to neuronal Ca2+ 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 IRFCa2+ (dilation) and negative IRFNE (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.