Modeling the role of osmotic forces in the cerebrovascular response to CO2.

Increases in blood osmolarity have been shown to exert a vasodilatory effect on cerebral and other vasculature, with accompanying increases in blood flow. It has also been shown that, through an influence on blood concentration of the bicarbonate ion and pH, changes in blood levels of CO2 can alter blood osmolarity sufficiently to have an impact on vessel diameter. We propose here that this phenomenon plays a previously unappreciated role in CO2-mediated vasodilation, and present a biophysical model of osmotically driven vasodilation. Our model, which is based on literature data describing CO2-dependent changes in blood osmolarity and hydraulic conductivity (Lp) of the blood-brain barrier, is used to predict the change in cerebral blood flow (CBF) associated with osmotic forces arising from a specific hypercapnic challenge. Modeled changes were then compared with actual CBF changes determined using arterial spin-labeling (ASL) MRI. For changes in the arterial partial pressure of CO2 (PaCO2) of 20 mmHg, our model predicted increases of 80% from baseline CBF with a temporal evolution that was comparable to the measured hemodynamic responses. Our modeling results suggest that osmotic forces could play a significant role in the cerebrovascular response to CO2.

Absolute quantification of resting oxygen metabolism and metabolic reactivity during functional activation using QUO2 MRI.

We have recently described an extension of calibrated MRI, which we term QUO2 (for QUantitative O(2) imaging), providing absolute quantification of resting oxidative metabolism (CMRO(2)) and oxygen extraction fraction (OEF(0)). By combining BOLD, arterial spin labeling (ASL) and end-tidal O(2) measurements in response to hypercapnia, hyperoxia and combined hyperoxia/hypercapnia manipulations, and the same MRI measurements during a task, a comprehensive set of vascular and metabolic measurements can be obtained using a generalized calibration model (GCM). These include the baseline absolute CBF in units of ml/100g/min, cerebrovascular reactivity (CVR) in units of %Δ CBF/mm Hg, M in units of percent, OEF(0) and CMRO(2) at rest in units of µmol/100g/min, percent evoked CMRO(2) during the task and n, the value for flow-metabolic coupling associated with the task. The M parameter is a calibration constant corresponding to the maximal BOLD signal that would occur upon removal of all deoxyhemoglobin. We have previously shown that the GCM provides estimates of the above resting parameters in grey matter that are in excellent agreement with literature. Here we demonstrate the method using functionally-defined regions-of-interest in the context of an activation study. We applied the method under high and low signal-to-noise conditions, corresponding respectively to a robust visual stimulus and a modified Stroop task. The estimates fall within the physiological range of literature values, showing the general validity of the GCM approach to yield non-invasively an extensive array of relevant vascular and metabolic parameters.

Reduction of physiological noise with independent component analysis improves the detection of nociceptive responses with fMRI of the human spinal cord.

The evaluation of spinal cord neuronal activity in humans with functional magnetic resonance imaging (fMRI) is technically challenging. Major difficulties arise from cardiac and respiratory movement artifacts that constitute significant sources of noise. In this paper we assessed the Correction of Structured noise using spatial Independent Component Analysis (CORSICA). FMRI data of the cervical spinal cord were acquired in 14 healthy subjects using gradient-echo EPI. Nociceptive electrical stimuli were applied to the thumb. Additional data with short TR (250 ms, to prevent aliasing) were acquired to generate a spatial map of physiological noise derived from Independent Component Analysis (ICA). Physiological noise was subsequently removed from the long-TR data after selecting independent components based on the generated noise map. Stimulus-evoked responses were analyzed using the general linear model, with and without CORSICA and with a regressor generated from the cerebrospinal fluid region. Results showed higher sensitivity to detect stimulus-related activation in the targeted dorsal segment of the cord after CORSICA. Furthermore, fewer voxels showed stimulus-related signal changes in the CSF and outside the spinal region, suggesting an increase in specificity. ICA can be used to effectively reduce physiological noise in spinal cord fMRI time series.

Task-related BOLD responses and resting-state functional connectivity during physiological clamping of end-tidal CO(2).

Carbon dioxide (CO(2)), a potent vasodilator, is known to have a significant impact on the blood-oxygen level dependent (BOLD) signal. With the growing interest in studying synchronized BOLD fluctuations during the resting state, the extent to which the apparent synchrony is due to variations in the end-tidal pressure of CO(2) (PETCO(2)) is an important consideration. CO(2)-related fluctuations in BOLD signal may also represent a potential confound when studying task-related responses, especially if breathing depth and rate are affected by the task. While previous studies of the above issues have explored retrospective correction of BOLD fluctuations related to arterial PCO(2), here we demonstrate an alternative approach based on physiological clamping of the arterial CO(2) level to a near-constant value. We present data comparing resting-state functional connectivity within the default-mode-network (DMN), as well as task-related BOLD responses, acquired in two conditions in each subject: 1) while subject's PETCO(2) was allowed to vary spontaneously; and 2) while controlling subject's PETCO(2) within a narrow range. Strong task-related responses and areas of maximal signal correlation in the DMN were not significantly altered by suppressing fluctuations in PETCO(2). Controlling PETCO(2) did, however, improve the performance of retrospective physiological noise correction techniques, allowing detection of additional regions of task-related response and resting-state connectivity in highly vascularized regions such as occipital cortex. While these results serve to further rule out systemic physiological fluctuations as a significant source of apparent resting-state network connectivity, they also demonstrate that fluctuations in arterial CO(2) are one of the factors limiting sensitivity in task-based and resting-state fMRI, particularly in regions of high vascular density. This must be considered when comparing subject groups who might exhibit differences in respiratory physiology or breathing patterns.

Magnetic resonance imaging of resting OEF and CMRO2 using a generalized calibration model for hypercapnia and hyperoxia.

We present a method allowing determination of resting cerebral oxygen metabolism (CMRO2) from MRI and end-tidal O2 measurements acquired during a pair of respiratory manipulations producing different combinations of hypercapnia and hyperoxia. The approach is based on a recently introduced generalization of calibrated MRI signal models that is valid for arbitrary combinations of blood flow and oxygenation change. Application of this model to MRI and respiratory data during a predominantly hyperoxic gas manipulation yields a specific functional relationship between the resting BOLD signal M and the resting oxygen extraction fraction OEF0. Repeating the procedure using a second, primarily hypercapnic, manipulation provides a different functional form of M vs. OEF0. These two equations can be readily solved for the two unknowns M and OEF0. The procedure also yields the resting arterial O2 content, which when multiplied by resting cerebral blood flow provides the total oxygen delivery in absolute physical units. The resultant map of oxygen delivery can be multiplied by the map of OEF0 to obtain a map of the resting cerebral metabolic rate of oxygen consumption (CMRO2) in absolute physical units. Application of this procedure in a group of seven human subjects provided average values of 0.35 ± 0.04 and 6.0 ± 0.7% for OEF0 and M, respectively in gray-matter (M valid for 30 ms echo-time at 3T). Multiplying OEF0 estimates by the individual values of resting gray-matter CBF (mean 52 ± 5 ml/100 g/min) and the measured arterial O2 content gave a group average resting CMRO2 value of 145 ± 30 µmol/100 g/min. The method also allowed the generation of maps depicting resting OEF, BOLD signal, and CMRO2.

Elimination of visually evoked BOLD responses during carbogen inhalation: implications for calibrated MRI.

Breathing a mixture of 10% CO(2) with 90% O(2) (referred to here as carbogen-10) increases blood flow due to the vasodilatory effect of CO(2), and raises blood O(2) saturation due to the enriched oxygen level. These effects both tend to reduce the level of deoxygenated hemoglobin in brain tissues, thereby reducing the potential for further increases in BOLD contrast. In the present study, blocks of intense visual stimulation (60s) were presented amid longer blocks (180s) during which subjects breathed various fractional concentrations (0-100%) of carbogen-10 diluted with medical air. When breathing undiluted carbogen-10, the BOLD response to visual stimulation was reduced below the level of noise against the background of the carbogen-10 response. At these concentrations, the total (visual+carbogen) BOLD response amplitude (7.5±1.0%, n=6) converged toward that seen with carbogen alone (7.5±1.0%, n=6). In spite of the almost complete elimination of the visual BOLD response, pseudo-continuous arterial spin-labeling on a separate cohort indicated a largely preserved perfusion response (89±34%, n=5) to the visual stimulus during inhalation of carbogen-10. The previously discussed observations suggest that venous saturation can be driven to very high levels during carbogen inhalation, a finding which has significant implications for calibrated MRI techniques. The latter methods involve estimation of the relative change in venous O(2) saturation by expressing activation-induced BOLD signal increases as a fraction of the maximal BOLD signal M that would be observed as venous saturation approaches 100%. While the value of M has generally been extrapolated from much smaller BOLD responses induced using hypercapnia or hyperoxia, our results suggest that these effects could be combined through carbogen inhalation to obtain estimates of M based on larger BOLD increases. Using a hybrid BOLD calibration model taking into account changes in both blood flow and arterial oxygenation, we estimated that inhalation of carbogen-10 led to an average venous saturation of 91%, allowing us to compute an estimated M value of 9.5%.

BOLD signal responses to controlled hypercapnia in human spinal cord.

Functional MRI of the spinal cord is challenging due to the small cross section of the cord and high level of physiological noise. Though blood oxygenation level-dependent (BOLD) contrast has been used to study specific responses of the spinal cord to various stimuli, it has not been demonstrated using a controlled stimulus. In this paper, we use hypercapnic manipulation to study the sensitivity and specificity of functional MRI in the human cervical spinal cord. Simultaneous MR imaging in the brain and spinal cord was performed for direct comparison with the brain, in which responses to hypercapnia have been more extensively characterized. Original contributions include: (i) prospectively controlled hypercapnic changes in end-tidal PCO(2), (ii) simultaneous recording of BOLD responses in the brain and spinal cord, and (iii) generation of statistical maps of BOLD responses throughout the brain and spinal cord, taking into account physiological noise sources. Results showed significant responses in all subjects both in the brain and the spinal cord. In anatomically-defined regions of interest, mean percent changes were 0.6% in the spinal cord and 1% in the brain. Analysis of residual variance demonstrated significantly larger contribution of physiological noise in the spinal cord (P<0.005). To obtain more reliable results from fMRI in the spinal cord, it will be necessary to improve sensitivity through the use of highly parallelized coil arrays and better modeling of physiological noise. Finely, we believe that the use of controlled global stimuli, such as hypercapnia, will help assess the effectiveness of new acquisition techniques.

BDNF, relative preference, and reward circuitry responses to emotional communication.

Brain derived neurotrophic factor (BDNF) regulates neural development and synaptic transmission. We have tested the hypothesis that functional variation in the BDNF gene (Val66Met polymorphism, rs6265) affects brain reward circuitry encoding human judgment and decision-making regarding relative preference. We quantified relative preference among faces with emotional expressions (angry, fearful, sad, neutral, and happy) by a keypress procedure performed offline to measure effort traded for viewing time. Keypress-based relative preferences across the ensemble of faces were mirrored significantly by fMRI signal in the orbitofrontal cortex, amygdala, and hippocampus when passively viewing these faces. For these three brain regions, there was also a statistically significant group difference by BDNF genotype in the fMRI responses to the emotional expressions. In comparison with Val/Met heterozygotes, Val/Val individuals preferentially sought exposure to positive emotions (e.g., happy faces) and had stronger regional fMRI activation to aversive stimuli (e.g., angry, fearful, and sad faces). BDNF genotype accounted for approximately 30% of the variance in fMRI signal that mirrors keypress responses to these stimuli. This study demonstrates that functional allelic variation in BDNF modulates human brain circuits processing reward/aversion information and relative preference transactions.

Investigations on spinal cord fMRI of cats under ketamine.

Functional magnetic resonance imaging (fMRI) of the spinal cord has been the subject of intense research for the last ten years. An important motivation for this technique is its ability to detect non-invasively neuronal activity in the spinal cord related to sensorimotor functions in various conditions, such as after spinal cord lesions. Although promising results of spinal cord fMRI have arisen from previous studies, the poor reproducibility of BOLD activations and their characteristics remain a major drawback. In the present study we investigated the reproducibility of BOLD fMRI in the spinal cord of cats (N=9) by repeating the same stimulation protocol over a long period (approximately 2 h). Cats were anaesthetized with ketamine, and spinal cord activity was induced by electrical stimulation of cutaneous nerves of the hind limbs. As a result, task-related signals were detected in most cats with relatively good spatial specificity. However, BOLD response significantly varied within and between cats. This variability was notably attributed to the moderate intensity of the stimulus producing a low amplitude haemodynamic response, variation in end-tidal CO(2) during the session, low signal-to-noise ratio (SNR) in spinal fMRI time series and animal-specific vascular anatomy. Original contributions of the present study are: (i) first spinal fMRI experiment in ketamine-anaesthetized animals, (ii) extensive study of intra- and inter-subject variability of activation, (iii) characterisation of static and temporal SNR in the spinal cord and (iv) investigation on the impact of CO(2) end-tidal level on the amplitude of BOLD response.

Detection of multiple pathways in the spinal cord using q-ball imaging.

Magnetic resonance diffusion tensor imaging (DTI) has been extensively applied to the spinal cord for depicting its architecture and for assessing its integrity following spinal lesions. However, DTI is limited in representing complex white matter architecture, notably in the presence of crossing fibres. Recently, q-ball imaging (QBI) has been proposed as a new method for recovering complex white matter architecture. We applied this technique to both ex vivo and in vivo spinal cords of cats using a 3T scanner. For the purpose of comparison, gradients have been applied in 55 and 100 encoding directions and b-values varied from 800 to 3000 s/mm(2). As a result, QBI was able to retrieve crossing fibre information, where the DTI approach was constrained in a unique diffusion direction. To our knowledge, this is the first study demonstrating the benefits of QBI for detecting the presence of longitudinal, commissural and dorso-ventral fibres in the spinal cord. It is a first step towards in vivo characterization of the healthy and injured human spinal cord using high angular resolution diffusion imaging and QBI.

In vivo DTI of the healthy and injured cat spinal cord at high spatial and angular resolution.

Spinal cord diffusion tensor imaging (DTI) is challenging in many ways: the small size of the cord, physiological motion and susceptibility artifacts pose daunting obstacles to the acquisition of high-quality data. Here, we present DTI results computed from in vivo studies of the healthy and injured spinal cord of five cats. Both high spatial (1.1 mm3) and angular (55 directions) resolutions were used to optimise modelling of the diffusion process. Also, particular effort was directed towards a strategy that limits susceptibility artifacts. For validation purposes, acquisitions were repeated in two cats before and after making a spinal lesion. As a result, various axonal trajectories were identified by tractography including dorsal and ventral columns as well as lateral tracts. Also, fibre bundles showed robust disruption at the site of spinal cord injuries (partial and complete) via tractography, accompanied with significantly lower fractional anisotropy values at the site of lesions. Important outcomes of this work are (i) tractography-based localisation of anatomical tracts in the thoraco-lumbar spinal cord and (ii) in vivo assessment of axonal integrity following experimental spinal cord injury.

A temporal comparison of BOLD, ASL, and NIRS hemodynamic responses to motor stimuli in adult humans.

In this study, we have preformed simultaneous near-infrared spectroscopy (NIRS) along with BOLD (blood oxygen level dependent) and ASL (arterial spin labeling)-based fMRI during an event-related motor activity in human subjects in order to compare the temporal dynamics of the hemodynamic responses recorded in each method. These measurements have allowed us to examine the validity of the biophysical models underlying each modality and, as a result, gain greater insight into the hemodynamic responses to neuronal activation. Although prior studies have examined the relationships between these two methodologies through similar experiments, they have produced conflicting results in the literature for a variety of reasons. Here, by employing a short-duration, event-related motor task, we have been able to emphasize the subtle temporal differences between the hemodynamic parameters with a high contrast-to-noise ratio. As a result of this improved experimental design, we are able to report that the fMRI measured BOLD response is more correlated with the NIRS measure of deoxy-hemoglobin (R = 0.98; P < 10(-20)) than with oxy-hemoglobin (R = 0.71), or total hemoglobin (R = 0.53). This result was predicted from the theoretical grounds of the BOLD response and is in agreement with several previous works [Toronov, V.A.W., Choi, J.H., Wolf, M., Michalos, A., Gratton, E., Hueber, D., 2001. "Investigation of human brain hemodynamics by simultaneous near-infrared spectroscopy and functional magnetic resonance imaging." Med. Phys. 28 (4) 521-527.; MacIntosh, B.J., Klassen, L.M., Menon, R.S., 2003. "Transient hemodynamics during a breath hold challenge in a two part functional imaging study with simultaneous near-infrared spectroscopy in adult humans". NeuroImage 20 1246-1252.; Toronov, V.A.W., Walker, S., Gupta, R., Choi, J.H., Gratton, E., Hueber, D., Webb, A., 2003. "The roles of changes in deoxyhemoglobin concentration and regional cerebral blood volume in the fMRI BOLD signal" Neuroimage 19 (4) 1521-1531]. These data have also allowed us to examine more detailed measurement models of the fMRI signal and comment on the roles of the oxygen saturation and blood volume contributions to the BOLD response. In addition, we found high correlation between the NIRS measured total hemoglobin and ASL measured cerebral blood flow (R = 0.91; P < 10(-10)) and oxy-hemoglobin with flow (R = 0.83; P < 10(-05)) as predicted by the biophysical models. Finally, we note a significant amount of cross-modality, correlated, inter-subject variability in amplitude change and time-to-peak of the hemodynamic response. The observed co-variance in these parameters between subjects is in agreement with hemodynamic models and provides further support that fMRI and NIRS have similar vascular sensitivity.

Comparison of physiological noise at 1.5 T, 3 T and 7 T and optimization of fMRI acquisition parameters.

Previous studies have shown that under some conditions, noise fluctuations in an fMRI time-course are dominated by physiological modulations of the image intensity with secondary contributions from thermal image noise and that these two sources scale differently with signal intensity, susceptibility weighting (TE) and field strength. The SNR of the fMRI time-course was found to be near its asymptotic limit for moderate spatial resolution measurements at 3 T with only marginal gains expected from acquisition at higher field strengths. In this study, we investigate the amplitude of image intensity fluctuations in the fMRI time-course at magnetic field strengths of 1.5 T, 3 T, and 7 T as a function of image resolution, flip angle and TE. The time-course SNR was a similar function of the image SNR regardless of whether the image SNR was modulated by flip angle, image resolution, or field strength. For spatial resolutions typical of those currently used in fMRI (e.g., 3 x 3 x 3 mm(3)), increases in image SNR obtained from 7 T acquisition produced only modest increases in time-course SNR. At this spatial resolution, the ratio of physiological noise to thermal image noise was 0.61, 0.89, and 2.23 for 1.5 T, 3 T, and 7 T. At a resolution of 1 x 1 x 3 mm(3), however, the physiological to thermal noise ratio was 0.34, 0.57, and 0.91 for 1.5 T, 3 T and 7 T for TE near T2*. Thus, by reducing the signal strength using higher image resolution, the ratio of physiologic to image noise could be reduced to a regime where increased sensitivity afforded by higher field strength still translated to improved SNR in the fMRI time-series.

Simultaneous recording of task-induced changes in blood oxygenation, volume, and flow using diffuse optical imaging and arterial spin-labeling MRI.

Increased neural activity in brain tissue is accompanied by an array of supporting physiological processes, including increases in blood flow and the rates at which glucose and oxygen are consumed. These responses lead to secondary effects such as alterations in blood oxygenation and blood volume, and are ultimately the primary determinants of the amplitude and temporal signature of the blood oxygenation level-dependent (BOLD) signal used prevalently to map brain function. We have performed experiments using a combination of optical and MRI-based imaging methods to develop a more comprehensive picture of the physiological events accompanying activation of primary motor cortex during a finger apposition task. Temporal profiles for changes in tissue hemoglobin concentrations were qualitatively similar to those observed for MRI-based flow and oxygenation signals. Quantitative analysis of these signals revealed peak changes of +16 +/- 2% for HbO, -13 +/- 2% for HbR, +8 +/- 3% for total Hb, +83 +/- 9% for cerebral blood flow, and +1.4 +/- 0.1% for the BOLD MRI signal. A mass balance model was used to estimate the change in rate of oxidative metabolism implied by the optical and flow measurements, leading to a computed value of +47 +/- 5%. It should be noted that the optical and MRI observations may in general reflect changes over different volumes of tissue. The ratio of fractional changes in oxidative metabolism to fractional change in blood flow was found to be 0.56 +/- 0.08, in general agreement with previous studies of flow-metabolism coupling.

Can the cerebral metabolic rate of oxygen be estimated with near-infrared spectroscopy?

We have measured the changes in oxy-haemoglobin and deoxy-haemoglobin in the adult human brain during a brief finger tapping exercise using near-infrared spectroscopy (NIRS). The cerebral metabolic rate of oxygen (CMRO2) can be estimated from these NIRS data provided certain model assumptions. The change in CMRO2 is related to changes in the total haemoglobin concentration, deoxy-haemoglobin concentration and blood flow. As NIRS does not provide a measure of dynamic changes in blood flow during brain activation, we relied on a Windkessel model that relates dynamic blood volume and flow changes, which has been used previously for estimating CMRO2 from functional magnetic resonance imaging (fMRI) data. Because of the partial volume effect we are unable to quantify the absolute changes in the local brain haemoglobin concentrations with NIRS and thus are unable to obtain an estimate of the absolute CMRO2 change. An absolute estimate is also confounded by uncertainty in the flow-volume relationship. However, the ratio of the flow change to the CMRO2 change is relatively insensitive to these uncertainties. For the linger tapping task, we estimate a most probable flow-consumption ratio ranging from 1.5 to 3 in agreement with previous findings presented in the literature, although we cannot exclude the possibility that there is no CMRO2 change. The large range in the ratio arises from the large number of model parameters that must be estimated from the data. A more precise estimate of the flow-consumption ratio will require better estimates of the model parameters or flow information, as can be provided by combining NIRS with fMRI.

Oxidative metabolism and the detection of neuronal activation via imaging.

Recent years have witnessed a great growth of interest in non-invasive imaging methods, such as functional magnetic resonance imaging (fMRI) and positron emission tomography (PET), permitting identification of brain structures that mediate specific cognitive and behavioural tasks in humans. Because these techniques use physiological responses such as increased perfusion or metabolism as surrogate indicators of evoked neuronal electrical activity, understanding the role of these processes in sustaining the information processing function of neurons is vital to the proper interpretation of functional neuroimaging data. An ultimate goal of these non-invasive techniques is to approach the sensitivity and spatial resolution of earlier autoradiographic methods, which have repeatedly demonstrated exquisitely detailed delineations of neuronal response patterns using metabolic glucose uptake as a physiological tag. Although glucose is generally metabolized in conjunction with oxygen, technical challenges in imaging tissue oxygen consumption in vivo have limited the use of this complementary process in the detection of neuronal activation. In this article we review concepts linking cerebral blood flow and metabolism to neuronal activation, and compare functional imaging techniques that exploit these relationships. We also describe recently introduced MRI based methods for measurement of oxygen consumption and assess the relative contributions of different metabolic pathways during neuronal activation. Our calculations suggest that the bulk of the energy demand evoked during stimulation of neurons in visual cortex is met through oxidative metabolism of glucose, supporting the use of oxygen uptake as a marker for increased neuronal electrical activity.

Relating axonal injury to functional recovery in MS.

A patient was followed after the new onset of hemiparesis from relapse of MS with serial MR spectroscopy and functional MRI. The association of clinical improvement with recovery of N-acetylaspartate, a marker of neuronal integrity, and progressive reduction of abnormally large functional MRI cortical activation with movement demonstrates that dynamic reorganization of the motor cortex accompanies remission of MS.

Investigation of BOLD signal dependence on cerebral blood flow and oxygen consumption: the deoxyhemoglobin dilution model.

The relationship between blood oxygenation level-dependent (BOLD) MRI signals, cerebral blood flow (CBF), and oxygen consumption (CMR(O2)) in the physiological steady state was investigated. A quantitative model, based on flow-dependent dilution of metabolically generated deoxyhemoglobin, was validated by measuring BOLD signals and relative CBF simultaneously in the primary visual cortex (V1) of human subjects (N = 12) during graded hypercapnia at different levels of visual stimulation. BOLD and CBF responses to specific conditions were averaged across subjects and plotted as points in the BOLD-CBF plane, tracing out lines of constant CMR(O2). The quantitative deoxyhemoglobin dilution model could be fit to these measured iso-CMR(O2) contours without significant (P

Linear coupling between cerebral blood flow and oxygen consumption in activated human cortex.

The aim of this study was to test the hypothesis that, within a specific cortical unit, fractional changes in cerebral blood flow (CBF) and cerebral metabolic rate of oxygen consumption (CMR(O(2))) are coupled through an invariant relationship during physiological stimulation. This aim was achieved by simultaneously measuring relative changes in these quantities in human primary visual cortex (V1) during graded stimulation with patterns designed to selectively activate different populations of V1 neurons. Primary visual cortex was delineated individually in each subject by using phase-encoded retinotopic mapping. Flow-sensitive alternating inversion recovery MRI, in conjunction with blood oxygenation-sensitive MRI and hypercapnic calibration, was used to monitor CBF and CMR(O(2)). The stimuli used included (i) diffuse isoluminant chromatic displays; (ii) high spatial-frequency achromatic luminance gratings; and (iii) radial checkerboard patterns containing both color and luminance contrast modulated at different temporal rates. Perfusion responses to each pattern were graded by varying luminance and/or color modulation amplitudes. For all stimulus types, fractional changes in blood flow and oxygen uptake were found to be linearly coupled in a consistent ratio of approximately 2:1. The most potent stimulus produced CBF and CMR(O(2)) increases of 48 +/- 5% and 25 +/- 4%, respectively, with no evidence of a plateau for oxygen consumption. Estimation of aerobic ATP yields from the observed CMR(O(2)) increases and comparison with the maximum possible anaerobic ATP contribution indicate that elevated energy demands during brain activation are met largely through oxidative metabolism.

Stimulus-dependent BOLD and perfusion dynamics in human V1.

Blood oxygenation level-dependent (BOLD) fMRI signals often exhibit pronounced over- or undershoot upon changes in stimulation state. Current models postulate that this is due to the delayed onset or decay of perfusion-dependent attenuating responses such as increased cerebral blood volume or oxygen consumption, which are presumed to lag behind the rapid adjustment of blood flow rate to a new steady-state level. If this view is correct, then BOLD overshoot amplitudes in a specific tissue volume should be correlated with steady-state increases in perfusion, independent of stimulus type. To test this prediction, we simultaneously recorded BOLD and relative perfusion signals in primary visual cortex while inducing graded perfusion increases with three types of visual stimulus. Two of these, a diffuse chromatic stimulus with no luminance variation and a very high spatial frequency luminance grating, did not produce detectable BOLD overshoot (or undershoot) when an equal mean luminance baseline was used. Radial checkerboard stimuli, however, caused pronounced over/undershoot of both BOLD and perfusion signals even when temporal mean luminance was held constant and stimulus contrast was adjusted to produce the same steady-state blood flow increases evoked by the other stimuli. Transient amplitudes were relatively invariant in spite of large changes in steady-state response, demonstrating nonlinear BOLD and perfusion step responses in human V1. These findings suggest that, rather than a purely tissue-specific biomechanical or metabolic phenomenon, BOLD overshoot and undershoot represent transient features in the perfusion signal whose effects may be amplified by slowly evolving blood volume changes.