Coupling of the cortical hemodynamic response to cortical and thalamic neuronal activity

Accurate interpretation of functional MRI (fMRI) signals requires knowledge of the relationship between the hemodynamic response and the neuronal activity that underlies it. Here we address the question of coupling between pre- and postsynaptic neuronal activity and the hemodynamic response in rodent somatosensory (Barrel) cortex in response to single-whisker deflection. Using full-field multiwavelength optical imaging of hemoglobin oxygenation and electrophysiological recordings of spiking activity and local field potentials, we demonstrate that a point hemodynamic measure is influenced by neuronal activity across multiple cortical columns. We demonstrate that the hemodynamic response is a spatiotemporal convolution of the neuronal activation. Therefore, positive hemodynamic response in one cortical column might be explained by neuronal activity not only in that column but also in the neighboring columns. Thus, attempts at characterizing the neurovascular relationship based on point measurements of electrophysiology and hemodynamics may yield inconsistent results, depending on the spatial extent of neuronal activation. The finding that the hemodynamic signal observed at a given location is a function of electrophysiological activity over a broad spatial region helps explain a previously observed increase of local vascular response beyond the saturation of local neuronal activity. We also demonstrate that the oxy- and total-hemoglobin hemodynamic responses can be well approximated by space-time separable functions with an antagonistic center-surround spatial pattern extending over several millimeters. The surround "negative" hemodynamic activity did not correspond to observable changes in neuronal activity. The complex spatial integration of the hemodynamic response should be considered when interpreting fMRI data.