Spatial relationship between neuronal activity and BOLD functional MRI
Introduction
Since its introduction in 1992 Bandettini et al., 1992, Kwong et al., 1992, Ogawa et al., 1992 blood oxygenation level-dependent (BOLD) functional magnetic resonance imaging (fMRI) has revolutionized cognitive neurosciences by allowing the foci of cortical “activity” to be visualized in vivo in a noninvasive manner. The BOLD contrast originates from the intravoxel magnetic field inhomogeneity induced by paramagnetic deoxyhemoglobin (deoxyHb) sequestered in red blood cells that are compartmentalized within blood vessels. The magnetic susceptibility differences between the deoxyHb-containing compartments and the surrounding space generate magnetic field gradients around the boundaries of these compartments. Therefore, perturbation of regional deoxyHb content will alter the signal intensities in MR images sensitized to BOLD contrast. Such regional perturbation occurs as the result of enhanced neuronal activity and metabolism during sensory (Engel et al., 1994), motor (Kim et al., 1993), or cognitive (Wagner et al., 1998) functions.
While BOLD-based neuroimaging studies have provided unprecedented amount of insights into the workings of the human brain in vivo, the explanatory power of BOLD fMRI is currently limited since there is a fundamental gap in our understanding of the linkage between the observed BOLD contrast and the underlying neuronal physiology. In particular, the extent to which the magnitude and spatial scale of the BOLD signal correlates with neuronal physiology remains elusive Attwell and Iadecola, 2002, Ugurbil et al., 2003 To this end, a growing body of results suggests a predominantly linear coupling between BOLD and neuronal activity. For example, a recent study by Ogawa et al. (2000) demonstrated a linear relationship between somatosensory-evoked potentials and BOLD signals for brief stimulation durations. Rees et al. (2000) and Heeger et al. (2000) demonstrated a linear correlation between BOLD contrast in humans and averaged spike rate in a monkey cortical area during stimulation with nearly identical stimuli. A similarly linear relationship was observed also in anesthetized monkeys by Logothetis et al. (2001), a study in which single-unit responses were acquired simultaneously with BOLD signals inside the MRI scanner.
While the above results suggest that the fundamental coupling between BOLD and the underlying neuronal activity is approximately linear, important questions remain about the spatial scale over which the linear coupling remains valid. Is the hypothesized linear coupling between BOLD and neuronal activity invariant across the different spatial scales of the cortical architecture? Can we assume a universal linearity from the spatial scale of entire cortical areas (several millimeters to centimeters) to individual cortical columns (sub-millimeter)? Such universal linearity would be surprising, since the BOLD signal is not a direct measure of neuronal activity per se. Rather, it is a complex convolution of changes in cerebral metabolic rate of oxygen (CMRO2), cerebral blood flow (CBF), and cerebral blood volume (CBV) following focal neuronal activity. Therefore, spatial specificity of BOLD response becomes a critical issue in examining its relationship to the underlying neuronal activity. It is in principle possible to record from a site where neuronal activity is absent but there exists a large BOLD response as the latter effect “spills over” from an adjacent active region due to lack of spatial specificity. Consequently, it cannot be assumed a priori that the correlation between the BOLD signals and the underlying neural activity is linear irrespective of spatial considerations, especially at the sub-millimeter scale at which the finest details of spatial organization in cortex have been observed.
In the present study, we utilize a novel recording scheme to characterize the relationship between BOLD fMRI signals and the underlying neuronal modulation within the same animals from multiple electrode recording sites over a large cortical area. Our study extends recent results (Logothetis et al., 2001) in which correlations between BOLD and neuronal activity were characterized from a single electrode penetration site. With multiple recording sites spaced over the entire cat area 18 on the lateral gyrus, the results of our study provide a detailed look into the spatial relationship between neuronal activity and BOLD functional signals.
Section snippets
Animal preparation
Juvenile cats between postnatal weeks 5 and 12 were used for the present studies. All procedures conformed to federal and state guidelines for treatment of animals and were approved by the University of Minnesota IACUC. Animals were initially anesthetized with ketamine (10–25 mg/kg, i.v.) and xylazine (2.5 mg/kg). The animals were intubated and artificially ventilated under isoflurane anesthesia (0.8–1.5%) in a N2O/O2 mixture of 70:30 throughout the experiment. Contact lenses were used to focus
Results and discussion
It is known that across the surface of visual cortex, orientation and spatial frequency preferences change on a scale of hundreds of micrometers Hubel and Wiesel, 1962, Hubener et al., 1997, Issa et al., 2000, Maffei and Fiorentini, 1977; therefore, it is necessary to co-register single-unit recordings and fMRI data within that level of precision. For such accurate registration that would allow the evaluation of both the relationship between magnitudes of the neuronal and fMRI signals and the
Conclusions
In this study, we have developed a method to directly compare the BOLD and responses from multiple single-unit sites from the same cortical area. Identical stimuli were used for BOLD and single-unit studies to directly correlate the spatiotemporal coupling between the BOLD responses and the underlying neuronal activity. By taking advantage of the consecutive recording scheme, we were able to obtain correlative data from a large number of recording sites across the cat area 18.
The results of our
Acknowledgements
We thank Drs. T. Duong and N. Harel for their assistance during some of the initial experiments. N. Port from NIH provided the initial recording chambers for testing. This work was supported by the NIH (MH57180, MH67530, NS38295, RR08079), the W.M. Keck Foundation, and the Whitaker Foundation, NARSAD.
References (44)
- et al.
The neural basis of functional brain imaging signals
Trends Neurosci.
(2002) The elusive initial dip
NeuroImage
(2001)Linearity of synaptic interactions in the assembly of receptive fields in cat visual cortex
Curr. Opin. Neurobiol.
(1994)- et al.
Spatial frequency rows in the striate visual cortex
Vision Res.
(1977) - et al.
Noise, neural codes and cortical organization
Curr. Opin. Neurobiol.
(1994) Time as coding space?
Curr. Opin. Neurobiol.
(1999)- et al.
How accurate is magnetic resonance imaging of brain function?
Trends Neurosci.
(2003) - et al.
GABA-mediated representation of temporal information in rat barrel cortex
NeuroReport
(1999) - et al.
Time course EPI of human brain function during task activation
Magn. Reson. Med.
(1992) - et al.
Localized cerebral blood flow response at submillimeter columnar resolution
Proc. Natl. Acad. Sci. U. S. A.
(2001)
fMRI of human visual cortex
Nature
Focal physiological uncoupling of cerebral blood flow and oxidative metabolism during somatosensory stimulation in human subjects
Proc. Natl. Acad. Sci. U. S. A.
Cortical point-spread function and long-range lateral interactions revealed by real-time optical imaging of macaque monkey primary visual cortex
J. Neurosci.
The suppressive influence of moving textured backgrounds on responses of cat striate neurons to moving bars
J. Neurophysiol.
Origin of negative blood oxygenation level-dependent fMRI signals
J. Cereb. Blood Flow Metab.
Spikes versus BOLD: what does neuroimaging tell us about neuronal activity?
Nat. Neurosci.
Evaluation of the early response in fMRI in individual subjects using short stimulus duration
Magn. Reson. Med.
Receptive field, binocular interaction and functional architecture in the cat's visual cortex
J. Physiol. (London)
Spatial relationships among three columnar systems in cat area 17
J. Neurosci.
Spatial frequency maps in cat visual cortex
J. Neurosci.
Linearity of summation of synaptic potentials underlying direction selectivity in simple cells of the cat visual cortex
Science
Functional imaging of human motor cortex at high magnetic field
J. Neurophysiol.
Cited by (0)
- 1
Present address: University of Pittsburgh, Pittsburgh, PA, USA.