BOLD and Perfusion Response to Finger-Thumb Apposition after Acetazolamide Administration: Differential Relationship to Global Perfusion

Abstract

The authors studied the effects of altering global cerebral blood flow on both blood oxygen level–dependent (BOLD) response and perfusion response to finger-thumb apposition. A PICORE/QUIPSS II protocol was used to collect interleaved BOLD-weighted and perfusion-weighted images on eight finger-thumb apposition trials. Subjects were studied on a drug-free day and on a day when acetazolamide was administered between the second and third trials. After acetazolamide administration, resting cortical perfusion increased an average of 20% from preadministration levels, whereas the BOLD response to finger-thumb apposition decreased by an average of 35% in the S1M1 hand area. Contrary to predictions from the exhausted cerebrovascular reserve hypothesis and the oxygen limitation model, an effect of acetazolamide on cerebral blood flow response in the S1M1 hand area was not observed. Across the acetazolamide trials, BOLD response was inversely correlated with resting cortical perfusion for individual subject data. These results suggest that resting perfusion affects the magnitude of the BOLD response and is thus an important confounding factor in fMRI studies, and that the physiologic systems that increase cerebral blood flow in response to acetazolamide administration and systems that increase cerebral blood flow in response to altered neural activity appear to have additive effects.

Keywords

Cerebral blood flow BOLD contrast Acetazolamide Arterial spin labeling Functional magnetic resonance imaging

Most functional magnetic resonance imaging (fMRI) studies use blood oxygen level–dependent (BOLD) contrast to measure changes in neuronal functioning associated with alterations of behavioral state (Kwong, 1995; Ogawa et al., 1992). BOLD contrast reflects variation in the deoxyhemoglobin content of cerebral vessels induced by changes in the level of neuronal activity. Ultimately, deoxyhemoglobin levels depend on CMRo₂, CBF, and cerebral blood volume (CBV) (Buxton et al., 1998). Thus, changes in neural activity are transduced through a hemodynamic response system to produce the observed BOLD effect. Most analyses of BOLD contrast do not correct for global levels of blood flow and blood volume. However, several recent papers have discussed the use of the global T₂*-weighted signal as a covariate in BOLD contrast studies (Aquirre et al., 1998; Desjardins et al., 2001). The issue of whether BOLD response is independent of the global signal is central to the question of how to account for the effects of global signal variation on local BOLD response to varying experimental conditions (Aquirre et al., 1998; Cohen et al., 2002). In the current study, we examined the more circumscribed issue of whether the BOLD and perfusion responses to behavioral challenge are independent of manipulations of global, resting, cortical perfusion.

Results of the few fMRI studies that have manipulated baseline brain perfusion are contradictory. One set of findings indicated a reduction of BOLD response to behavioral challenge with an increase in baseline signal (Bandettini et al., 1995; Bruhn et al., 1994; Cohen et al., 2002), whereas a study by Corfield et al. (2001) found no interaction between BOLD response to visual stimulation and the partial pressure of carbon dioxide (Pco₂) when comparing normocapnia and hypercapnia. In another study, BOLD contrast to visual stimulation peaked during mild hypercapnia (Posse et al., 2001).

Like the results of Corfield and et al.'s BOLD study, several human studies have found no effect of the manipulation of global flow on the size of the CBF response associated with behavioral challenges. The independence of CBF response and global flow has been found using ¹³³Xe-inhalation CBF measurements (Maximilian et al., 1980), positron emission tomography (Ramsay et al., 1993), and perfusion-weighted magnetic resonance imaging (Li et al., 2000). That manipulations of global flow do not alter the magnitude of CBF response to behavioral activation was also observed between normocapnic and hypercapnic conditions (Shimosegawa et al., 1995) and among 40, 50, and 60 mm Hg CO₂ conditions (Kemna et al., 2001). A coupling of global perfusion and CBF response to behavioral activation may exist at extreme values of Pco₂, adding to the complexity of results (Kemna et al., 2001; Shimosegawa et al., 1995). It appears that the independence of experimentally manipulated global flow and CBF response to behavioral activation occurs within a window ranging from normocapnia to moderate hypercapnia.

Two issues underlie much of the contradiction described above: (1) the untested assumption that alterations of basal CBF ought to have similar effects on perfusion response and BOLD response to behavior challenge, and (2) the comparison of simple change scores (ΔCBF) with percent signal change (%BOLD) data. In the present study, we measured cortical cerebral perfusion and BOLD response to behavioral challenge under drug-free and acetazolamide conditions. We used multislice PICORE/QUIPSS II to obtain interleaved perfusion-weighted and BOLD-weighted images (Wong et al., 1998). Acetazolamide has cerebrovascular effects similar to CO₂. It is a selective inhibitor of erythrocyte carbonic anhydrase (EC 4.2.1.1) that inhibits the removal of CO₂ in the brain (Vorstrup et al., 1984), decreases tissue pH (Vorstrup et al., 1984), and increases both CBF (Ehrenreich et al., 1961; Vorstrup et al., 1984) and venous oxygen saturation (Vorstrup et al., 1984). Yet acetazolamide appears to have little effect on CMRo₂ and little effect on neuronal function (Vorstrup et al., 1984). Acetazolamide has been shown to reliably increase cerebral perfusion in studies using magnetic resonance methods (Yen et al., 2002). The effects of acetazolamide on BOLD and perfusion response to finger apposition was measured using both simple and percent-change scores.

MATERIALS AND METHODS

Subjects

Five healthy volunteers (two women and three men) with a mean (SD) age of 32.2 (8.3) years participated in the study. All were right handed and had earned graduate and/or undergraduate degrees. To minimize risks associated with acetazolamide exposure, subjects were excluded if they had a history of abnormal potassium or calcium blood levels, kidney or liver disease, glaucoma, heart disease, allergic reaction to sulfonamides, migraine, or if they were currently pregnant. Subjects gave written informed consent to participate. The University of California San Diego Institutional Review Board approved the study.

Procedure

Study protocol. Subjects were scanned on two different occasions, separated by at least 24 hours. Both scanning procedures were identical except that on one occasion, subjects received intravenous acetazolamide (Diamox; Lederle Laboratories, Div. American Cyanamid Co., Wayne, NJ, U.S.A.). For two subjects, the drug-free procedure occurred before the acetazolamide procedure.

On each study day, an axial localizer scan of thirteen 5-mm slices was acquired to assure that the subject's head was aligned in the transverse plane. Sagittal localizer images (thirteen 5-mm slices) were then acquired to identify the locations of the functional slices to be acquired. Five contiguous axial functional slices (8-mm thick) were chosen to cover the motor strip for each subject and generally spanned from the most superior point of the arch of the corpus callosum to the most superior point of the brain. A neuroradiologist then obtained and reviewed a magnetic resonance angiogram of the Circle of Willis to screen for vascular abnormalities that would contra-indicate acetazolamide injection.

Eight functional runs, spaced about 5 minutes apart, were acquired on each occasion. During each run, 110 interleaved control and inversion-tagged images were acquired for each of five slices over a 3-minute 40-second time period using the multislice PICORE/QUIPSS II protocol (Wong et al., 1998). The image-acquisition protocol weighted mean images (control + tag) towards BOLD contrast and difference images (control – tag) towards perfusion contrast. Subjects alternated between 24 seconds of rest and 24 seconds of right-handed, sequential, finger-to-thumb apposition performed at a brisk but comfortable pace. Auditory cues signaled the start of each rest or tapping period. On the acetazolamide day, two runs (pre-ACZ) were acquired before acetazolamide administration. The remaining six runs (post-ACZ) were acquired starting approximately 5 minutes after acetazolamide injection and were spaced 5 minutes apart.

Acetazolamide administration. Before scanning, a neuroradiologist placed an intravenous line in the subject's left arm while the subject was lying in the magnet. The line was periodically flushed with saline. After the second pre-ACZ run, a 1-g bolus of acetazolamide in a saline solution was delivered over a period of 2 to 3 minutes. One gram of acetazolamide produced an effective dose of 15 mg/kg for subject 1, 12 mg/kg for subject 2, and 13.5 mg/kg for all other subjects. Each subject received a dose sufficient to inhibit completely carbonic anhydrase activity in most tissues (Vorstrup et al., 1984). All subjects reported tingling in the perioral region soon after injection, and subject 1 reported auditory distortions similar to those occasionally noted in migraine auras.

Scanning parameters. Images were collected on a GE Signa 1.5T system using an inserted local head gradient coil of our own design (Wong et al., 1992). High-resolution anatomical images were acquired using a spoiled gradient recalled echo (SPGR) protocol (128 sagittal slices, 1.4-mm thickness, 0.938 × 0.938-mm in-plane resolution, 11.4-millisecond repetition time, 4.4-millisecond echo time, 10° flip angle, 240-mm field of view).

We acquired interleaved control and inversion-tagged images using PICORE tagging with the QUIPSS II modification and single-shot blipped echo planar imaging (25-millisecond echo time, 2,000-millisecond repetition time, TI₁/TI₂ = 600/1,200 milliseconds) (Wong et al., 1998). The in-plane resolution was 3.75 × 3.75 mm with a 240-mm field of view. The inversion pulse was a 15-millisecond adiabatic hyperbolic secant (μ = 10, β = 800 s⁻¹). A 10-cm tag inversion band was used, with a 1-cm gap to the edge of the imaging region.

Image construction

Gray matter segmentation. Gray matter voxels were segmented from voxels containing white matter and CSF using a three-step process. First, MPRAGE images were passed through a three-dimensional, high-pass filter in two iterations to reduce equipment-related signal inhomogeneities across the field of view (Jernigan et al., 2001). The skull was then stripped from each subject's filtered high-resolution anatomical images using the Analysis of Functional NeuroImage (AFNI) 3dIntracranial routine (Cox, 1996). Finally, a three-compartment Gaussian mixture model representing gray matter, white matter, and CSF was fit to the signal intensity histogram produced by the filtered SPGR images. A mean, standard deviation, and weight for each compartment were estimated. The downhill simplex method of Nelder and Mead (1965) was used to estimate the nine parameters that minimized the least-squares fit of the three-compartment model to the histogram. Cutting scores to segment gray matter, white matter, and CSF were derived from the two cusp points defined by the overlap of the three Gaussian compartments (one cusp was between gray and white matter compartments; the other was between gray matter and CSF compartments). A composite image of the three compartments was color-coded and the resulting three-dimensional image was compared visually with the filtered SPGR image to verify the segmentation. Occasionally this comparison indicated that the threshold between gray matter and CSF had to be altered to match the anatomical image. A gray matter mask was derived from the composite segmented image.

Perfusion and BOLD images. Subtracting from each image the average of the previous and subsequent images (and alternating the sign so that each difference corresponded to control – tag) generated raw perfusion images. Temporally local subtraction rather than subtraction of the mean of all tagged images from the mean of all control images was used in order to reduce the confounding of the flow signal by signal drift across the image-acquisition period. As with other arterial spin labeling methods, the magnitude of the PICORE/QUIPSS II difference image (control – tag) is directly proportional to local CBF, lacking only a global calibration factor to convert imager units into the physiologic units of CBF. No attempt was made to calculate the calibration factor, and the data were analyzed as perfusion-weighted imager units. Averaging across temporally adjacent tag and control images generated BOLD images. Cortical perfusion and BOLD images were coregistered with the anatomical images using the visualization tool 3dNudge from the AFNI library (Cox, 1996).

Dependent variables

BOLD response (ΔS and %BOLD). We removed abnormal signal spikes in the BOLD response in the left sensorimotor hand area (see definition below) by calculating on a case-by-case basis the standard deviation of the BOLD response across time separately for each finger-thumb apposition and rest block. For each rest and tapping block, we used box and whisker plots of the standard deviations across subjects to identify outliers. Blocks with extremely high standard deviations (i.e., more than three interquartile ranges outside the distribution of values across subjects for a given block) were considered to be too variable for analysis and were not processed further. Three blocks from one subject were excluded. Additionally, the BOLD and perfusion response volumes were separately motion-corrected using a three-dimensional iterated least-squares minimization process (AFNI 3dvolreg; Cox, 1996) with Fourier interpolation. Each time point was registered to the central time point in the run.

A reference vector representing the experimental design (zero coded rest, one coded figure-thumb apposition) was convolved with a two-step impulse-response function, yielding a simple trapezoidal model of the hemodynamic response with a graded lag of two repetition time periods. For each subject, this hemodynamic model and a linear term to adjust for linear drifts across the time series were fit to the BOLD time series in each voxel using least-squares estimation. The fit coefficient associated with the hemodynamic model measured the BOLD response to finger-thumb apposition in imager signal units. With the previously described coding scheme, the hemodynamic fit coefficient for the BOLD time series is equal to ΔS, the difference of the mean image signal averaged over rest trials subtracted from the mean image signal averaged over finger-thumb apposition trials, expressed in local magnetic resonance intensity units. The fit coefficient (ΔS) maps were normalized to Talairach coordinates and resampled into 1-mm³ voxels before further analyses (Talairach and Tournoux, 1988). Because the large majority of voxels in the functionally defined S1M1 hand area were in the gray matter, we calculated mean ΔS for positive fit coefficients across all voxels in the S1M1 hand area rather than cortically segment this volume as described previously. Percent BOLD response (%BOLD) was calculated by dividing mean ΔS by the average magnetic resonance signal in the S1M1 hand area obtained during the rest trials.

Perfusion response (ΔCBF and %CBF). We removed abnormal signals producing negative flow values from the perfusion time series. These abnormal signals represented less than 1% of the total data. To derive a perfusion fit coefficient, the regression model applied to the BOLD time series was applied to the perfusion time series. The resulting fit coefficient was a measure of perfusion response to finger-thumb apposition ΔCBF) expressed in local imager units. The method for calculating %BOLD was used to calculate %CBF for voxels with positive fit coefficients. The ΔCBF and %CBF response maps were normalized to Talairach coordinates before further analyses. Perfusion response measures were averaged over all voxels in the S1M1 hand area.

Resting global cortical perfusion. The study protocol involved comparing BOLD signals across alternating blocks of rest and finger-apposition trials. To account for delayed hemodynamic response to motor activation, we defined a rest period as the epoch beginning with the third trial after the previous motor activation block and extending into the first two trials of the subsequent block. To avoid noise from the large volume of low flow white matter voxels from attenuating the global CBF measure, we examined resting global perfusion only in cortically segmented voxels. The resting cortical perfusion measure for each subject was derived in two steps. First, perfusion values for each voxel were averaged over all trials within rest periods and the subsequent mean slices were combined into a three-dimensional volume. As with the ΔCBF response measure, we removed abnormal signal spikes producing negative flow signals from the time series. Second, we multiplied the resulting three-dimensional volume by the segmented gray matter mask to derive a resting, global cortical-perfusion image.

Resting cortical magnetic resonance signal S. Magnetic resonance signal values (S) from the resting periods were identified and processed as described for resting global cortical perfusion. Using methods described for the BOLD response measures, we removed abnormal signals from the resting magnetic resonance signal time series. Because acetazolamide predominantly affects the gray matter T₂* signal (Corfield et al., 2001), the segmented gray matter mask was used to derive resting, global cortical-signal values.

S1M1 region of interest. Four of the study's dependent measures, ΔS, %BOLD, ΔCBF, and %CBF, were obtained for the region of interest in the left sensorimotor area (S1M1). In order to identify the outer bounds of intersubject variability in the location of the hand S1M1 area in the left hemisphere, we constructed a mask that reflected the normal range of anatomical variation of the S1M1 region (mask volume = 8,428 mm³). The mask, based on Talairach landmarks derived from Carey et al.'s study (Carey et al., 2000), was superimposed on each subject's three-dimensional BOLD and perfusion response maps. Study measures were averaged over all voxels in the mask. These analyses were performed using routines from the AFNI library (Cox, 1996).

Statistical analyses

We analyzed the mean %BOLD, %CBF, ΔS, and ΔCBF in the S1M1 area of the left hemisphere, the global resting cortical perfusion values, and global resting cortical magnetic resonance signal values. Because the distribution of these variables tended to be skewed, we analyzed the natural logarithm of these values (Miller, 1986). Effects of acetazolamide and trial on study measures were analyzed in a two (condition) by two (mean of the pre-ACZ trials versus mean of the post-ACZ administration trials) within-subjects repeated-measures analysis of variance, treating subjects as a random factor. Partial η² was used to measure effect sizes as the amount of variance accounted for by the acetazolamide condition and by its interaction with the trial factor (Cohen, 1973).

RESULTS

Five minutes after acetazolamide administration in trial 2, global resting cortical perfusion increased 19% from preadministration values and averaged 20% above preadministration resting values throughout the entire postadministration period (F_1,4 = 20.68, P = 0.01, partial η² = 0.84; Fig. 1B). Similarly, global resting cortical-image values, S, were significantly greater after acetazolamide administration than in the drug-free condition (F_1,4 = 9.44, P = 0.037, partial η² = 0.70; Fig. 1A).

FIG. 1.

(A) Effect of acetazolamide on MR signal intensity. High resolution, axial MPRAGE image most proximal to the tagging slice for a representative subject (i). Mean perfusion image averaged across the two baseline runs before acetazolamide administration for the same subject and slice level (ii). Mean perfusion image averaged across the six runs after acetazolamide injection for the same subject and slice level (iii). (B) Resting cortical perfusion. Effects of acetazolamide and trial global resting cortical perfusion (mean ± SD, N = 5).

Although acetazolamide increased resting cortical perfusion, it had no significant effect on the magnitude of perfusion response (ΔCBF) to finger-thumb apposition in the S1M1 hand area of the left hemisphere (F_1,4 = 0.044, P = 0.845, partial η² = 0.01; Fig. 2A). Furthermore, acetazolamide administration did not significantly influence relative perfusion (%CBF: F_1,4 = 1.078, P = 0.358, partial η² = 0.21; Fig. 2B). To determine whether the PICORE/QUIPSS II measure of perfusion was insensitive to physiologic manipulations, we tested whether perfusion differed from zero during finger-thumb apposition. Both the ΔCBF and %CBF measures showed robust finger-thumb apposition effects averaged over trials and drug condition (ΔCBF: F_1,4 = 686.86, P<0.001, partial η² = 0.99; %CBF: F_1,4 = 1,272.96, P<0.001, partial η² = 0.99).

FIG. 2.

Effects of acetazolamide and trial on the perfusion and BOLD response to finger-thumb apposition compared with rest in the primary sensorimotor hand area of the left hemisphere. (A) Perfusion response (ΔCBF). (B) Percent perfusion response (%CBF). (C) BOLD response (ΔS). (D) Relative BOLD response (%BOLD). Data are mean ± SD; N = 5. Acetazolamide was administered after trial 2.

Although acetazolamide did not significantly affect CBF response, it did decrease the BOLD response expressed either as raw magnetic resonance signal change (ΔS) or as percentage signal change (ΔS/S_resting) in the S1M1 hand area (drug by trial interaction, ΔS: F_1,4 = 9.05, P = 0.04, partial η² = 0.69; %BOLD: F_1,4 = 11.34, P = 0.028, partial η² = 0.74; Figs. 2C and 2D). Five minutes after acetazolamide administration, ΔS decreased 28% from preadministration values and %BOLD decreased 29%. Averaged across all acetazolamide trials, both ΔS and %BOLD declined 35% from preadministration values.

To observe how global resting cortical perfusion and %BOLD response were related across trials in individual participants, we regressed individual %BOLD values onto resting cortical perfusion. All five regression weights were negative; greater resting global cortical perfusion was associated with smaller %BOLD values in the acetazolamide condition. Regression of ΔS onto resting global cortical perfusion produced the same pattern of association. That all five regression weights would be less than zero is unlikely to occur by chance (sign test: P = 0.031, one-tailed). Median correlation between %BOLD and resting cortical perfusion among the five subjects during the acetazolamide run was −0.73.

DISCUSSION

An intravenous bolus of 1 g acetazolamide reduced the BOLD response by about 35% and increased global resting cortical perfusion by 20%, though perfusion response was unaltered by acetazolamide administration. These findings confirm Bruhn et al.'s (1994) prediction that the reduction of BOLD response after acetazolamide injection would be related to increased cerebral perfusion. An inverse correlation between global cortical perfusion and BOLD response is compatible with the increased BOLD response observed in the somatosensory region of Sprague-Dawley rats when hyperventilation-induced hypocapnia was compared with normocapnia (Hsu et al., 1998) and with evidence that caffeine ingestion reduces baseline CBF and increases BOLD response (Mulderink et al., 2002). The finding that ΔCBF and %CBF values were unaltered by manipulations of global resting cortical perfusion adds to the growing literature showing the independence between pharmacologically or physiologically induced changes in global brain perfusion and local CBF response to behavioral challenge (Li et al., 2000; Maximilian et al., 1980; Ramsay et al., 1993; Shimosegawa et al., 1995). Given that the effects of 1 g acetazolamide on global CBF are similar to the effects of mild to moderate hypercapnia (Cotev et al., 1968; Severinghaus and Cotev, 1968), the independence of global CBF and perfusion response to behavioral challenge observed in our study is compatible with the Pco₂ window found in the study conducted by Kemna et al. (2001). Acetazolamide administration also increased the resting T₂*-weighted signal compared with the drug-free condition, replicating findings from other laboratories of widespread T₂*-weighted signal activation associated with hypercapnia (Corfield et al., 2001; Posse et al., 2001).

The exhaustion of cerebrovascular reserve hypothesis (Bruhn et al., 1994) and the oxygen limitation model (Buxton and Frank, 1997) have been advanced to explain the dependence of BOLD response on global perfusion. A nearly exhausted cerebrovascular reserve would produce a ceiling effect on the CBF response to behavioral activation, and the constrained CBF response would in turn constrain the BOLD response. The oxygen limitation model is based on the observation that there is little oxygen reserve for neurons to use when increased neural activity increases oxygen demand (Buxton and Frank, 1997; Fennema et al., 1989). Consequently, increased cerebral oxygen consumption during behavioral activation would require an increased flux of oxygen from the capillaries to neural tissue (Raichle, 1998). Furthermore, the rate of oxygen delivery to mitochondria in neurons and glia depends on partial pressure of oxygen (Po₂) gradients. If CBF in the resting baseline were increased by pharmacologic means without altering CMRo₂, the Po₂ gradient would increase and smaller changes in CBF would be required to deliver a sufficient amount of oxygen to meet the increased metabolic demand produced by increased behavioral activity. The diminished change in CBF would attenuate the change in oxygenated hemoglobin delivery and lead to smaller BOLD responses. Like the exhaustion of cerebrovascular reserve hypothesis, the oxygen limitation model predicts a joint alteration of BOLD and perfusion response with acetazolamide-induced increases in resting cortical perfusion. However, the prediction of a coupled BOLD and perfusion response to acetazolamide administration was not supported by the current study.

If neither the oxygen limitation model nor the exhausted cerebrovascular reserve hypothesis can account for the inverse relationship between BOLD response and global cortical perfusion, how can an inverse coupling between BOLD response and perfusion be explained? An understanding of this coupling involves an analysis of the transverse relaxation parameter,

R_{2}^{*} = \frac{1}{T_{2}^{*}} = R_{2 o}^{*} + R_{2 d H b *}^{*}

In this model, the total relaxation rate constant is taken to be the sum of an intrinsic relaxation rate (R*_2o) and a second term, R*_2dHb, due to the local concentration of deoxyhemoglobin. Specifically, in the absence of deoxyhemoglobin the relaxation rate would be R*_2o. For simplicity, we will refer to the deoxyhemoglobin dependent term as R*₂ with the understanding that the full relaxation rate constant includes an additive term. Monte Carlo simulations (Boxerman et al., 1995a), analytical computations (Yablonsky and Haacke, 1994), and model-system experiments (Weiskoff et al., 1994) confirm a simple proportional relationship between R*₂ on the one hand, and deoxyhemoglobin concentration in blood (dHb) and blood volume (V) on the other:

R_{2}^{*} \infty V {[dHb]}^{β}

(1)

When β>1, the superlinear effects of deoxyhemoglobin concentration on the transverse relaxation rate account for the diffusion of water molecules through magnetic field gradients around blood vessels, especially capillaries, and for intravascular magnetic field offsets (Boxerman et al., 1995b; Davis et al., 1998; Grubb et al., 1974). The change in the transverse relaxation rate, ΔR₂*, due to behavioral, pharmacologic, or physiologic challenge is proportional to the difference between the activated relaxation rate, _aR₂*, and the baseline resting rate, _bR₂*. Reexpressing ΔR in terms of Eq. 1 and assuming signal relaxation is small leads to the following approximation (Davis et al., 1998):

\frac{Δ S}{S} = M {1 - (F^{(α - β)} m^{β})}

(2)

where m and f are the ratio of activated to resting baseline CMRo₂ and CBF, respectively. We assume that CBV and CBF are related by the power law CBV = CBF^α.

The BOLD response in Eq. 2 is determined by a scaling parameter, M, which depends on the baseline resting state, and by a dynamic term representing activation changes. Because changes in CMRo₂ due to neural activation tend to be small compared with CBF changes, the dynamic term in Eq. 2 is dominated by changes in flow (Davis et al., 1998). Given the usual values of β = 1.5 and α ≂ 0.38 at 1.5 T for normal physiologic conditions (Boxerman et al., 1995b; Grubb et al., 1974), f ^(α–β) is roughly the reciprocal of the flow ratio, making %BOLD approximately proportional to the fractional change in flow from resting to activated levels. For a fixed flow response to behavioral challenge, the ratio of baseline resting flow to activated flow, CBF_rest/(CBF_rest + ΔCBF), approaches unity as resting CBF levels increase. As Eq. 2 shows, the BOLD response decreases as the ratio of resting to activation flow approaches 1. The preceding discussion formalizes the intuitive argument that if the hemodynamic response to behavioral activation, Δresponse, remains constant after acetazolamide administration while the resting flow signal increases, then the ratio of Δresponse to the resting value will decrease. Notice that this account of the uncoupling of BOLD and perfusion response explains how %BOLD and ΔCBF can become uncoupled. Yet the account, which depends on the dynamic term in Eq. 2, fails to explain how %BOLD and %CBF might become uncoupled, as observed in the current study, because the argument predicts smaller changes in %CBF as resting perfusion increases.

A more satisfactory explanation of the relationship between baseline resting perfusion levels and BOLD response operates through parameter M in Eq. 2. M, the fractional reduction of the local signal at baseline from what it would be if no deoxyhemoglobin was present, modulates the dynamic range of BOLD response to behavioral challenge. Increases in resting CBF decrease deoxyhemoglobin content and M, whereas increases in resting CBV increase deoxyhemoglobin content and M. Thus whether alteration of M drives

\frac{Δ S}{S}

up or down depends on relative changes in resting CBF and CBV and not on changes in %CBF, which would operate through the dynamic term in Eq. 2. Because the exponential term in Grubb's power law relating flow to blood volume is less than 1, the finding that acetazolamide alters CBF more than CBV (Okazawa et al., 2001) implies that acetazolamide administration would reduce M, causing %BOLD to decline. The restricted dynamic range account of the relationship between %BOLD and resting perfusion explains the diminished BOLD response as a result of a floor effect for deoxyhemoglobin concentration under conditions of high resting CBF. Our data suggest that this floor effect is reached before perfusion response to behavioral challenge becomes constrained by an exhausted cerebrovascular reserve.

Because ΔS = S·M{1–(f^{−((β–α)} m^β)}, the effects of acetazolamide on the simple magnetic resonance signal difference between baseline resting and activation conditions, ΔS, is determined by resting baseline levels of S, in addition to baseline CBF, and CBV, and by behaviorally driven changes in CBF, CBV, and CMRo₂. Because acetazolamide increases S, the impact of acetazolamide on ΔS is determined by a balance of its effect on S and its effect on the model's physiological parameters. Given the multiplicity of determinants of ΔS, it is not surprising that studies of the effects of altered perfusion levels on behaviorally driven BOLD response have produced somewhat contradictory results. Whereas some studies have found no effect of perfusion level on ΔS (Corfield et al., 2001), some find increases in ΔS with presumed increased perfusion (Posse et al., 2001) and others find an inverse relationship of ΔS with measured or presumed perfusion levels (Bruhn et al., 1994; Hsu et al., 1998; Mulderink et al., 2002), as in the current study.

The gradient-recalled arterial spin labeling protocols typically used to obtain interleaved BOLD/perfusion images do not always separate BOLD effects from flow effects (Buxton, 2002). When averaging tagged and control images to obtain the BOLD signal, flow weighting of BOLD contrast can occur as inverted spins in the tagged image are averaged with fully relaxed spins from the control images. As spins flow into a slice, longitudinal magnetization is reduced in perfused vessels. Less longitudinal magnetization leads to smaller transverse magnetization after the longitudinal component has been tipped towards the transverse plane. With basic EPISTAR and PICORE protocols, higher rates of flow into an imaging slice would reduce BOLD contrast, producing a method-induced dependence of BOLD response on global perfusion (Buxton, 2002). To minimize flow weighting of BOLD images, the net magnetization delivered by flow must match the magnetization of the tissue across the image slice throughout the image-acquisition period. Fortunately, the temporal dynamics of the net magnetization flowing into the slice is known; it follows a saturation recovery curve (Buxton, 2002). The PICORE/QUIPSS II protocol minimizes flow weighting of the average signal by delivering a 90° saturation pulse to the image plane before applying the tag (Buxton, 2002; Wong et al., 1998). The saturation pulse matches the recovery of longitudinal magnetization to the flow weighting of BOLD images.

The validity of flow estimates for multislice arterial spin labeling methods can be corrupted by poorly chosen control images and by spatially varying delays in the transit of blood from the tagging region into the imaging slice (Buxton, 2002; Wong et al., 1998). The PICORE/QUIPSS II method used in the current study applies saturation pulses to the image slice in both the tagged and control states to minimize differences in slice profile and tissue magnetization that might be caused by differences in the pulse protocols generating the tagged and control states (Buxton, 2002; Wong et al., 1998). The dual saturation of the image slice minimizes errors in the subtraction of the slice images in tagged and control states. To minimize errors related to spatially varying transit delays, PICORE/QUIPSS II applies a saturation pulse to the tagged region after a specified inversion interval (TI1) has elapsed, effectively trimming the tail of the inversion bolus flowing out of the tagged region (Buxton, 2002; Wong et al., 1998). The sharp temporal demarcation of the tagged input bolus simplifies the timing of successive slice acquisitions so that tagged spins are delivered equally to each of the successively acquired slices.

Because factors that affect brain CBF separately from neural activity can alter the magnitude of regional BOLD response to a behavioral challenge, the interpretation of differences in BOLD response across experimental conditions needs to account for differences in baseline CBF (Cohen et al., 2002). This caveat applies most obviously to studies comparing healthy individuals and patients who are taking medications that could alter cerebrovascular tone or to those studied on and off a vasoactive medication (Bertel et al., 1987; Goldman et al., 1996; Miller et al., 1997; Shimamoto and Shimamoto, 1995). Although the vascular and neural effects of psychoactive drugs are hard to disentangle, our results suggest that the potential vascular effects of pharmacologic interventions need to be considered in drug-treatment studies of psychotropic medications using fMRI methods. Furthermore, given the complexity of factors that determine the relationship between BOLD response and perfusion levels, no one approach to adjusting BOLD response by perfusion level might work in all cases. However, by combining CBF data with a model of the BOLD signal, the effects of various adjustments can be evaluated quantitatively. This approach will be explored in future work.

Footnotes

Acknowledgments:

The authors thank Elizabeth Yoder for her critical review of an earlier version of this paper.

References

Aquirre

Zarahn

D'Esposito

(1998) The inferential impact of global signal covariates in functional neuroimaging analyses. Neuroimage 8:302–306

Bandettini

Wong

Binder

Rao

Jesmanowicz

Aaron

Lowry

Forster

Hinks

Hyde

(1995) Functional MR imaging using the BOLD approach: Dynamic characteristics and data analysis methods. In: Diffusion and perfusion magnetic resonance imaging ( Le Bihan

, ed), New York: Raven Press, Ltd., pp 335–349

Bertel

Marx

Conen

(1987) Effects of antihypertensive treatment on cerebral perfusion. Am J Med 82(Suppl 3B):29–36

Boxerman

Bandettini

Kwong

Baker

Davis

Rosen

Weisskoff

(1995a) The intravascular contribution to fMRI signal change: Monte Carlo modeling and diffusion-weighted studies in vivo. Magn Reson Med 34:4–10

Boxerman

Hamberg

Rosen

Weisskoff

(1995b) MR contrast due to intravascular magnetic susceptibility perturbations. Magn Reson Med 34:555–566

Bruhn

Kleinschmidt

Boecker

Merboldt

K-D

Hanicke

Frahm

(1994) The effect of acetazolamide on regional cerebral blood oxygenation at rest and under stimulation as assessed by MRI. J Cereb Blood Flow Metab 14:742–748

Buxton

(2002) Introduction to functional magnetic resonance imaging: principles and techniques. Cambridge UK: Cambridge University Press

Buxton

Frank

(1997) A model for the coupling between cerebral blood flow and oxygen metabolism during neural stimulation. J Cereb Blood Flow Metab 17:64–72

Buxton

Wong

Frank

(1998) Dynamics of blood flow and oxygenation changes during brain activation: The balloon model. Magn Reson Med 39:855–864

10.

Carey

Abbott

Egan

Tochon-Danguy

Bonnan

(2000) The functional neuroanatomy and long-term reproducibility of brain activation associated with as simple finger tapping task in older healthy volunteers: A serial PET study. Neuroimage 11:124–144

11.

Cohen

Ugurbil

Kim

S-G

(2002). Effect of basal conditions on the magnitude and dynamics of the blood oxygenation level-dependent fMRI response. J Cereb Blood Flow Metab 22:1042–1053

12.

Cohen

(1973) Eta-squared and partial eta-squared in fixed factor ANOVA designs. Educ Psychol Meas 33:107–112

13.

Corfield

Josephs

Adams

Turner

(2001) Does hypercapnia-induced cerebral vasodilation modulate the hemodynamic response to neural activation? Neuroimage 12:1207–1211

14.

Cotev

Lee

Severinghaus

(1968) The effects of acetazolamide on cerebral blood flow and cerebral tissue PO2. Anesthesiology 29:471–477

15.

Cox

(1996) AFNI: Software for analysis and visualization of functional magnetic resonance neuroimages. Comput Biomed Res 29:162–173

16.

Davis

Kwong

Weisskoff

Rosen

(1998) Calibrated functional MRI: Mapping the dynamics of oxidative metabolism. Proc Natl Acad Sci U S A 95:1834–1839

17.

Desjardins

Kiehl

Liddle

(2001) Removal of confounding effects of global signal in functional MRI analyses. Neuroimage 13:751–758

18.

Ehrenreich

Burns

Alman

Frazekas

(1961) Influence of acetazolamide on cerebral blood flow. Arch Neurol 5:227–232

19.

Fennema

Wessel

Faithful

Erdmann

(1989) Tissue oxygen tension in the cerebral cortex of the rabbit. Adv Exp Med Biol 248:451–460

20.

Goldman

Alexander

Zemishlany

Mukherjee

Sackeim

Prohovnik

(1996) Acute effects of haloperidol on cerebral cortex blood flow in normal and schizophrenic subjects. Biol Psychiatry 40:604–608

21.

Grubb

Raichle

Eichling

Ter-Pogossian

(1974) The effects of changes in PaCO₂ on cerebral blood volume, blood flow, and vascular mean transit time. Stroke 5:630–639

22.

Hsu

Hedlund

MacFall

(1998) Functional MRI of the rat somatosensory cortex: Effects of hyperventilation. Magn Reson Med 40:421–426

23.

Jernigan

Archibald

Fennema-Notestine

Gamst

Stout

Bonner

Hesselink

(2001) Effects of age on tissues and regions of the cerebrum and cerebellum. Neurobiol Aging 22:581–594

24.

Kemna

Posse

Tellmann

Schmitz

Herzog

(2001) Interdependence of regional and global cerebral blood flow during visual stimulation: An O-15-butanol positron emission tomography study. J Cereb Blood Flow Metab 21:664–670

25.

Kwong

(1995) Functional magnetic resonance imaging with echo planar imaging. Magn Reson Q 11:1–20

26.

T-Q

Kastrup

Moseley

Glover

(2000) Changes in baseline cerebral blood flow in humans do not influence regional cerebral blood flow response to photic stimulation. J Magn Reson Imaging 12:757–762

27.

Maximilian

Prohovnik

Risberg

(1980) Cerebral hemodynamic response to mental activation in normo- and hypercapnia. Stroke 11:342–347

28.

Miller

Rezai

Alliger

Andreasen

(1997) The effect of antipsychotic medication on relative cerebral blood perfusion in schizophrenia: Assessment with technitium-99m hexamethylpropyleneamine oxime single photon emission computed tomography. Biol Psychiatry 41:550–559

29.

Miller

(1986) Beyond ANOVA: basics of applied statistics. New York: John Wiley and Sons

30.

Mulderink

Gitelman

Mesulam

Parrish

(2002) On the use of caffeine as a contrast booster for BOLD fMRI studies. Neuroimage 15:37–44

31.

Nelder

Mead

(1965) Simplex method for function minimization. Comput J 7:308–313

32.

Ogawa

Tank

Menon

Ellermann

Kim

Merkle

Ugurbil

(1992). Intrinsic signal changes accompanying sensory stimulation: Functional brain mapping with magnetic resonance imaging. Proc Natl Acad Sci U S A 89:5951–5955

33.

Okazawa

Yamauchi

Sugimoto

Toyoda

Kishbe

Takahasi

(2001) Effects of acetazolamide on cerebral blood flow, blood volume, and oxygen metabolism: A positron emission tomography study with healthy volunteers. J Cereb Blood Flow Metab 21:1472–1479

34.

Posse

Kemna

Elghahwagi

Wiese

Kiselev

(2001) Effect of graded hypo- and hypercapnia on fMRI contrast in visual cortex: Quantification of T₂* changes by multiecho EPI. Magn Reson Med 46:264–271

35.

Raichle

(1998) Behind the scenes of functional brain imaging. Proc Natl Acad Sci U S A 95:765–772

36.

Ramsay

Murphy

Shea

Friston

Lammertsma

Clark

Adams

Guz

Frackowiak

(1993) Changes in global cerebral blood flow in humans: Effect on regional cerebral blood flow during a neural activation task. J Physiol 471:521–534

37.

Severinghaus

Cotev

(1968) Carbonic acidosis and cerebral vasodilation after Diamox. Scand J Clin Lab Invest Suppl 102:I:E

38.

Shimamoto

(1995) Nilvadipine increases cerebral blood flow in elderly hypertensives: Comparison with nifedipine. J Hum Hypertens 9:271–279

39.

Shimosegawa

Kanno

Hatazawa

Fujita

Iida

Miura

Murakami

Inugami

Ogawa

Itoh

(1995) Photic stimulation study of changing the arterial partial pressure level of carbon dioxide. J Cereb Blood Flow Metab 15:111–114

40.

Talairach

Tournoux

(1988) Co-planar stereotaxic atlas of the human brain, New York: Thieme

41.

Vorstrup

Henriksen

Paulson

(1984) Effect of acetazolamide on cerebral blood flow and cerebral metabolic rate for oxygen. J Clin Invest 74:1634–1639

42.

Weiskoff

Zuo

Boxerman

Rosen

(1994) Microscopic susceptibility variation and transverse relaxation: Theory and experiment. Magn Reson Med 31:601–610

43.

Wong

Bandettini

Hyde

(1992) Echoplanar imaging of the human brain using a three axis local gradient coil [abstract]. In: Proceedings of the 11th Annual Meeting of the Society of Magnetic Resonance Medicine, 1:105

44.

Wong

Buxton

Frank

(1998) Quantitative imaging of perfusion using a single subtraction (QUIPSS and QUIPSS II). Magn Reson Med 39:702–708

45.

Yablonsky

Haacke

(1994) Theory of NMR signal behavior in magnetically inhomogenous tissues: The static dephasing regime. Magn Reson Med 32:749–763

46.

Yen

Y-F

Field

Martin

Ari

Burdette

Moody

Takahasi

(2002) Test-retest reproducibility of quantitative CBF measurements using FAIR perfusion MRI and acetazolamide challenge. Magn Reson Med 40:421–426