Assessment of cerebrovascular reactivity with functional magnetic resonance imaging: comparison of CO(2) and breath holding

Cerebral blood flow (CBF) and oxygenation changes following both a simple breath holding test (BHT) and a CO(2) challenge can be detected with functional magnetic resonance imaging techniques. The BHT has the advantage of not requiring a source of CO(2) and acetazolamide and therefore it can easily be performed during a routine MR examination. In this study we compared global hemodynamic changes induced by breath holding and CO(2) inhalation with blood oxygenation level dependent (BOLD) and CBF sensitized fMRI techniques. During each vascular challenge BOLD and CBF signals were determined simultaneously with a combined BOLD and flow-sensitive alternating inversion recovery (FAIR) pulse sequence. There was a good correlation between the global BOLD signal intensity changes during breath holding and CO(2) inhalation supporting the notion that the BHT is equivalent to CO(2) inhalation in evaluating the hemodynamic reserve capacity with BOLD fMRI. In contrast, there was no correlation between relative CBF changes during both vascular challenges, which was probably due to the reduced temporal resolution of the combined BOLD and FAIR pulse sequence.     

Nonlinear blood oxygen level-dependent responses for transient activations and deactivations in V1 - insights into the hemodynamic response function with the balloon model

We report studies of the nonlinear nature of blood oxygen level-dependent (BOLD) responses to short transient deactivations in human visual cortex. Both functional magnetic resonance imaging (fMRI) and near-infrared spectroscopy (NIRS) have been used to compare and contrast the hemodynamic response functions (HRFs) associated with transient activation and deactivation in primary visual cortex. We show that signal decreases for short duration deactivations are smaller than corresponding signal increases in activation studies. Moreover, the standard balloon model of BOLD effects may be modified to account for the observed nonlinear nature of deactivations by appropriate changes to simple hemodynamic parameters without recourse to new assumptions about the nature of the coupling between activity and oxygen use.     

Dynamics and nonlinearities of the BOLD response at very short stimulus durations

In designing a functional imaging experiment or analyzing data, it is typically assumed that task duration and hemodynamic response are linearly related to each other. However, numerous human and animal studies have previously reported a deviation from linearity for short stimulus durations (<4 s). Here, we investigated nonlinearities of blood-oxygenation-level-dependent (BOLD) signals following visual stimulation of 5 to 1000 ms duration at two different luminance levels in human subjects. It was found that (a) a BOLD response to stimulus durations as short as 5 ms can be reliably detected; this stimulus duration is shorter than employed in any previous study investigating BOLD signal time courses; (b) the responses are more nonlinear than in any other previous study: the BOLD response to 1000 ms stimulation is only twice as large as the BOLD response to 5 ms stimulation although 200 times more photons were projected onto the retina; (c) the degree of nonlinearity depends on stimulus intensity; that is, nonlinearities have to be characterized not only by stimulus duration but also by stimulus features like luminance. These findings are especially of most practical importance in rapid event-related functional magnetic resonance imaging (fMRI) experimental designs. In addition, an 'initial dip' response--thought to be generated by a rapid increase in cerebral metabolic rate of oxygen metabolism (CMRO2) relative to cerebral blood flow--was observed and shown to colocalize well with the positive BOLD response. Highly intense stimulation, better tolerated by human subjects for short stimulus durations, causes early CMRO2 increase, and thus, the experimental design utilized in this study is better for detecting the initial dip than standard fMRI designs. These results and those from other groups suggest that short stimulation combined with appropriate experimental designs allows neuronal events and interactions to be examined by BOLD signal analysis, despite its slow evolution.     

An investigation of the impulse functions for the nonlinear BOLD response in functional MRI

Functional MRI (fMRI) based on blood oxygenation level dependent (BOLD) contrast can be used to detect hemodynamic responses to a broad range of stimuli. It however remains unclear in what fashion the BOLD response is a linear system, and how the impulse function differs with stimulation of varying duration. To address this question, fMRI using visual stimulation with a wide range of duration (0.5-12 s) was performed in six human volunteers. A strong linear correlation was shown on the full width at half maximum (r = 0.998) of the BOLD response curves and the area under the curves (r = 0.999) to the duration of stimulation. However, comparing the errors of the measured and predicted response curves, our results showed a poorer linearity at stimuli of shorter duration. By examining the impulse functions derived from different stimuli, based on the assumption that a linear convolution relationship existed, a higher differentiation was shown in the experiments with shorter stimuli (<3 s). Compared to the area under the impulse function derived from 12 s stimulation, with that obtained from 0.5, 1, 2, 3, 4, 8 s stimuli resulted in differences of 66.2, 33.5, 15.1, 5.4, 0.9, 7.9%, respectively. This study suggests a higher degree of nonlinearity in the BOLD signal changes due to stimuli of shorter duration, in agreement with earlier work.     

BOLD signal sign and transient vessels volume variation

The purpose of this work was to investigate the relation between BOLD signal sign and transient vessels volume variation induced by apnea. This stimulus consisting of breath holding after inspiration is able to induce a light slowing down in venous blood flow like in a sort of Valsalva maneuver. We observed diffuse negative BOLD responding areas at cortical level and a stronger negative response in correspondence of the main sinuses. These phenomena seem to be unrelated to a specific neural activity, appearing to be expressions of a mechanical variation in the hemodynamics. Our study suggests that particular care must be considered in the interpretation of fMRI findings, especially when patients with vascular-related cerebral diseases are involved.     

Origins of intersubject variability of blood oxygenation level dependent and arterial spin labeling fMRI: implications for quantification of brain activity

Accurate localization of brain activity using blood oxygenation level-dependent (BOLD) functional magnetic resonance imaging (fMRI) has been challenged because of the large BOLD signal within distal veins. Arterial spin labeling (ASL) techniques offer greater sensitivity to the microvasculature but possess low temporal resolution and limited brain coverage. In this study, we show that the physiological origins of BOLD and ASL depend on whether percent change or statistical significance is being considered. For BOLD and ASL fMRI data collected during a simple unilateral hand movement task, we found that in the area of the contralateral motor cortex the centre of gravity (CoG) of the intersubject coefficient of variation (CV) of BOLD fMRI was near the brain surface for percent change in signal, whereas the CoG of the intersubject CV for Z-score was in close proximity of sites of brain activity for both BOLD and ASL. These findings suggest that intersubject variability of BOLD percent change is vascular in origin, whereas the origin of inter-subject variability of Z-score is neuronal for both BOLD and ASL. For longer duration tasks (12 s or greater), however, there was a significant correlation between BOLD and ASL percent change, which was not evident for short duration tasks (6 s). These findings suggest that analyses directly comparing percent change in BOLD signal between pre-defined regions of interest using short duration stimuli, as for example in event-related designs, may be heavily weighted by large-vessel responses rather than neuronal responses.     

Nonlinear responses of cerebral blood volume, blood flow and blood oxygenation signals during visual stimulation

Hemodynamic-based functional magnetic resonance imaging (fMRI) techniques provide a great utility for noninvasive functional brain mapping. However, because the hemodynamic signals reflect underlying neural activity indirectly, characterization of these signals following brain activation is essential for experimental design and data interpretation. In this report, the linear (or nonlinear) responses to neuronal activation of three hemodynamic parameters based primarily on changes of cerebral blood volume, blood flow and blood oxygenation were investigated by testing these hemodynamic responses' additivity property. Using a recently developed fMRI technique that acquires vascular space occupancy (VASO), arterial spin labeling (ASL) perfusion and blood oxygenation level-dependent (BOLD) signals simultaneously, the additivity property of the three hemodynamic responses in human visual cortex was assessed using various visual stimulus durations. Experiments on healthy volunteers showed that all three hemodynamic-weighted signals responded nonlinearly to stimulus durations less than 4 s, with the degree of nonlinearity becoming more severe as the stimulus duration decreased. Vascular space occupancy and ASL perfusion signals showed similar nonlinearity properties, whereas the BOLD signal was the most nonlinear. These data suggest that caution should be taken in the interpretation of hemodynamic-based signals in fMRI.     

Regional homogeneity of resting-state fMRI contributes to both neurovascular and task activation variations

The task induced blood oxygenation level dependent signal changes observed using functional magnetic resonance imaging (fMRI) are critically dependent on the relationship between neuronal activity and hemodynamic response. Therefore, understanding the nature of neurovascular coupling is important when interpreting fMRI signal changes evoked via task. In this study, we used regional homogeneity (ReHo), a measure of local synchronization of the BOLD time series, to investigate whether the similarities of one voxel with the surrounding voxels are a property of neurovascular coupling. FMRI scans were obtained from fourteen subjects during bilateral finger tapping (FTAP), digit–symbol substitution (DSST) and periodic breath holding (BH) paradigm. A resting-state scan was also obtained for each of the subjects for 4 min using identical imaging parameters. Inter-voxel correlation analyses were conducted between the resting-state ReHo, resting-state amplitude of low frequency fluctuations (ALFF), BH responses and task activations within the masks related to task activations. There was a reliable mean voxel-wise spatial correlation between ReHo and other neurovascular variables (BH responses and ALFF). We observed a moderate correlation between ReHo and task activations (FTAP: r = 0.32; DSST: r = 0.22) within the task positive network and a small yet reliable correlation within the default mode network (DSST: r = − 0.08). Subsequently, a linear regression was used to estimate the contribution of ReHo, ALFF and BH responses to the task activated voxels. The unique contribution of ReHo was minimal. The results suggest that regional synchrony of the BOLD activity is a property that can explain the variance of neurovascular coupling and task activations; but its contribution to task activations can be accounted for by other neurovascular factors such as the ALFF.     

Visual cortex reactivity in sedated children examined with perfusion MRI (FAIR)

Sleeping and sedated children can respond to visual stimulation with a decrease in blood oxygenation level dependent (BOLD) functional MRI signal response. The contribution of metabolic and hemodynamic parameters to this inverse signal response is incompletely understood. It has been hypothesized that it is caused by a relatively greater increase of oxygen consumption compared to rCBF (regional cerebral blood flow) increase. We studied the rCBF changes during visual stimulation in four sedated children, aged 4-71 months, and four alert adults, with an arterial water spin labeling technique (FAIR) and BOLD fMRI in a 1.5T MR scanner. In the children, FAIR signal decreased by a mean of 0.96% (range 0.77-1.05) of the baseline periods of the non-selective images, while BOLD signal decreased by 2.03% (range 1.99-2.93). In the adults, FAIR and BOLD signal increased by 0.88% (range 0.8-0.99) and 2.63% (range 1.99-2.93), respectively. Thus, in the children, an rCBF increase could not be detected by perfusion MRI, but indications of a FAIR signal decrease were found. An rCBF decrease in the primary visual cortex during stimulation has not been reported previously, but it is a possible explanation for the negative BOLD response. Future studies will have to address if this response pattern is a consequence of age or sleep/sedation.     

Failure to direct detect magnetic field dephasing corresponding to ERP generation

Functional magnetic resonance imaging (fMRI) has become the method of choice for mapping brain activity in human subjects and detects changes in regional blood oxygenation and volume associated with local changes in neuronal activity. While imaging based on blood oxygenation level dependent (BOLD) contrast has good spatial resolution and sensitivity, the hemodynamic signal develops relatively slowly and is only indirectly related to neuronal activity. An alternative approach termed magnetic source magnetic resonance imaging (msMRI) is based on the premise that neural activity may be mapped by magnetic resonance imaging (MRI) with greater temporal resolution by detecting the local magnetic field perturbations associated with local neuronal electric currents. We used a hybrid ms/BOLD MRI method to investigate whether msMRI could detect signal changes that occur simultaneously at the time of the production of well-defined event-related potentials, the P300 and N170, in regions that previously have been identified as generators of these electrical signals. Robust BOLD activations occurred after some seconds, but we were unable to detect any significant changes in the T2*-weighted signal in these locations that correlated temporally with the timings of the evoked response potentials (ERPs).     

Deconvolution filtering: Temporal smoothing revisited

Inferences made from analysis of BOLD data regarding neural processes are potentially confounded by multiple competing sources: cardiac and respiratory signals, thermal effects, scanner drift, and motion-induced signal intensity changes. To address this problem, we propose deconvolution filtering, a process of systematically deconvolving and reconvolving the BOLD signal via the hemodynamic response function such that the resultant signal is composed of maximally likely neural and neurovascular signals. To test the validity of this approach, we compared the accuracy of BOLD signal variants (i.e., unfiltered, deconvolution filtered, band-pass filtered, and optimized band-pass filtered BOLD signals) in identifying useful properties of highly confounded, simulated BOLD data: (1) reconstructing the true, unconfounded BOLD signal, (2) correlation with the true, unconfounded BOLD signal, and (3) reconstructing the true functional connectivity of a three-node neural system. We also tested this approach by detecting task activation in BOLD data recorded from healthy adolescent girls (control) during an emotion processing task.     

Neural and vascular variability and the fMRI-BOLD response in normal aging

Neural, vascular and structural variables contributing to the blood oxygen level-dependent (BOLD) signal response variability were investigated in younger and older humans. Twelve younger healthy human subjects (six male and six female; mean age: 24 years; range: 19–27 years) and 12 older healthy subjects (five male and seven female; mean age: 58 years; range: 55–71 years) with no history of head trauma and neurological disease were scanned. Functional magnetic resonance imaging measurements using the BOLD contrast were made when participants performed a motor, cognitive or a breath hold (BH) task. Activation volume and the BOLD response amplitude were estimated for the younger and older at both group and subject levels. Mean activation volume was reduced by 45%, 40% and 38% in the elderly group during the motor, cognitive and BH tasks, respectively, compared to the younger. Reduction in activation volume was substantially higher compared to the reduction in the gray matter volume of 14% in the older compared to the younger. A significantly larger variability in the intersubject BOLD signal change occurred during the motor task, compared to the cognitive task. BH-induced BOLD signal change between subjects was significantly less-variable in the motor task-activated areas in the younger compared to older whereas such a difference between age groups was not observed during the cognitive task. Hemodynamic scaling using the BH signal substantially reduced the BOLD signal variability during the motor task compared to the cognitive task. The results indicate that the origin of the BOLD signal variability between subjects was predominantly vascular during the motor task while being principally a consequence of neural variability during the cognitive task. Thus, in addition to gray matter differences, the type of task performed can have different vascular variability weighting that can influence age-related differences in brain functional response.     

Direct measurement of oxygen extraction with fMRI using 6% CO2 inhalation

The blood-oxygenation-level-dependent (BOLD) signal is an indirect hemodynamic signal that is sensitive to cerebral blood flow (CBF), cerebral blood volume (CBV) and cerebral metabolic rate of oxygen. Therefore, the BOLD signal amplitude and dynamics cannot be interpreted unambiguously without additional physiological measurements, and thus, there remains a need for a functional magnetic resonance imaging (fMRI) signal, which is more closely related to the underlying neuronal activity. In this study, we measured CBF with continuous arterial spin labeling, CBV with an exogenous contrast agent and BOLD combined with intracortical electrophysiological recording in the primary visual cortex of the anesthetized monkey. During inhalation of 6% CO2, it was observed that CBF and CBV are not further increased by a visual stimulus, although baseline CBF for 6% CO2 is below the maximal value of CBF. In contrast, the electrophysiological response to the stimulation was found to be preserved during hypercapnia. As a consequence, the simultaneously measured BOLD signal responds negatively to a visual stimulation for 6% CO2 inhalation in the same voxels responding positively during normocapnia. These observations suggest that the fMRI response to a sensory stimulus for 6% CO2 inhalation occurs in the absence of a hemodynamic response, and it therefore directly reflects oxygen extraction into the tissue.     

Illuminating the BOLD signal: combined fMRI-fNIRS studies

Functional magnetic resonance imaging (fMRI) is currently combined with electrophysiological methods to identify the relationship between neuronal activity and the blood oxygenation level-dependent (BOLD) signal. Several processes like neuronal activity, synaptic activity, vascular dilation, blood volume and oxygenation changes underlie both response modalities, that is, the electrophysiological signal and the vascular response. However, accessing single process relationships is absolutely mandatory when aiming at a deeper understanding of neurovascular coupling and necessitates studies on the individual building blocks of the vascular response. Combined fMRI and functional near-infrared spectroscopy studies have been performed to validate the correlation of the BOLD signal to the hemodynamic changes in the brain. Here we review the current status of the integration of both technologies and judge these studies in the light of recent findings on neurovascular coupling.     

Neural activity-induced modulation of BOLD poststimulus undershoot independent of the positive signal

Despite intense research on the blood oxygenation level-dependent (BOLD) signal underlying functional magnetic resonance imaging, our understanding of its physiological basis is far from complete. In this study, it was investigated whether the so-called poststimulus BOLD signal undershoot is solely a passive vascular effect or actively induced by neural responses. Prolonged static and flickering black-white checkerboard stimulation with isoluminant grey screen as baseline condition were employed on eight human subjects. Within the same region of interest, the positive BOLD time courses for static and flickering stimuli were identical over the entire stimulus duration. In contrast, the static stimuli exhibited no poststimulus BOLD signal undershoot, whereas the flickering stimuli caused a strong BOLD poststimulus undershoot. To ease the interpretation, we performed an additional study measuring both BOLD signal and cerebral blood flow (CBF) using arterial spin labeling. Also for CBF, a difference in the poststimulus period was found for the two stimuli. Thus, a passive blood volume effect as the only contributor to the poststimulus undershoot comes short in explaining the BOLD poststimulus undershoot phenomenon for this particular experiment. Rather, an additional active neuronal activation or deactivation can strongly modulate the BOLD poststimulus behavior. In summary, the poststimulus time course of BOLD signal could potentially be used to differentiate neuronal activity patterns that are otherwise indistinguishable using the positive evoked response.     

Hemodynamic scaling of fMRI-BOLD signal: validation of low-frequency spectral amplitude as a scalability factor

Functional magnetic resonance imaging blood-oxygenation-level-dependent (fMRI-BOLD) signal representing neural activity may be optimized by discriminating MR signal components related to neural activity and those related to intrinsic properties of the cortical vasculature. The objective of this study was to reduce the hemodynamic change independent of neural activity to obtain a scaled fMRI-BOLD response using two factors, namely, low-frequency spectral amplitude (LFSA) and breath-hold amplitude (BHA). Ten subjects (age range, 22–38 years) were scanned during four task conditions: (a) rest while breathing room air, (b) bilateral finger tapping while breathing room air, (c) rest during a partial inspirational breath-hold, and (d) rest during moderate hypercapnia (breathing 5% CO₂, 20% O₂ and 75% N₂). In all subjects who breathed 5% CO₂, regions with significant BOLD response during breath-hold correlated significantly with the percent signal increase during 5% CO₂ inhalation. Finger-tapping-induced responses in the motor cortex were diminished to a similar extent after scaling using either LFSA or BHA. Inter- and intrasubject variation in the amplitude of the BOLD signal response reduced after hemodynamic scaling using LFSA or BHA. The results validated the hemodynamic amplitude scaling using LFSA with the earlier established BHA. LFSA free from motor-task contamination can be used to calibrate the fMRI-BOLD response in lieu of BHA or hypercapnia to minimize intra- and intersubject variation arising from vascular anatomy and vasodilative capacity.     

Post-stimulus response in hemodynamics observed by functional magnetic resonance imaging--difference between the primary sensorimotor area and the supplementary motor area

The blood oxygen level dependency (BOLD) contrast is a useful tool for functional neuroimaging based on the hemodynamic response to neuronal activation. We observed different hemodynamic responses in the BOLD signal between the primary sensorimotor area (SM1) and the supplementary motor area (SMA) in the sequential finger movement task. In the SMA, a stronger initial overshoot and a post-stimulus overshoot were observed. It was hypothesized from the time course analysis that the stronger initial overshoot reflected the activation of the SMA for motor control programming in the initial phase. Although the post-stimulus overshoot may be partially explained by cerebral blood flow (CBF) cerebral blood volume (CBV) uncoupling, its mechanism remained unknown. In the SM1, only the initial overshoot was observed and the level of BOLD signal was almost constant after the initial overshoot during the task period. These observations suggested that the BOLD signal is characterized by both CBF-CBV uncoupling and the neuronal activation characteristics in each region.     

Decoding neural events from fMRI BOLD signal: A comparison of existing approaches and development of a new algorithm

Neuroimaging methodology predominantly relies on the blood oxygenation level dependent (BOLD) signal. While the BOLD signal is a valid measure of neuronal activity, variances in fluctuations of the BOLD signal are not only due to fluctuations in neural activity. Thus, a remaining problem in neuroimaging analyses is developing methods that ensure specific inferences about neural activity that are not confounded by unrelated sources of noise in the BOLD signal. Here, we develop and test a new algorithm for performing semiblind (i.e., no knowledge of stimulus timings) deconvolution of the BOLD signal that treats the neural event as an observable, but intermediate, probabilistic representation of the system's state. We test and compare this new algorithm against three other recent deconvolution algorithms under varied levels of autocorrelated and Gaussian noise, hemodynamic response function (HRF) misspecification and observation sampling rate. Further, we compare the algorithms' performance using two models to simulate BOLD data: a convolution of neural events with a known (or misspecified) HRF versus a biophysically accurate balloon model of hemodynamics. We also examine the algorithms' performance on real task data. The results demonstrated good performance of all algorithms, though the new algorithm generally outperformed the others (3.0% improvement) under simulated resting-state experimental conditions exhibiting multiple, realistic confounding factors (as well as 10.3% improvement on a real Stroop task). The simulations also demonstrate that the greatest negative influence on deconvolution accuracy is observation sampling rate. Practical and theoretical implications of these results for improving inferences about neural activity from fMRI BOLD signal are discussed.     

Changes in MRI signal intensity during hypercapnic challenge under conscious and anesthetized conditions

Most functional magnetic resonance imaging (fMRI) studies in animals are conducted under anesthesia to minimize motion artifacts. However, methods and techniques have been developed recently for imaging fully conscious rats. Functional MRI studies on conscious animals report enhanced BOLD signal changes as compared to the anesthetized condition. In this study, rats were exposed to different concentrations of carbon dioxide (CO(2)) while conscious and anesthetized to test whether cerebrovascular reactivity may be contributing to these enhanced BOLD signal changes. Hypercapnia produced significantly greater increases in MRI signal intensity in fully conscious animals (6.7-13.3% changes) as when anesthetized with 1% isoflurane (3.2-4.9% changes). In addition, the response to hypercapnia was more immediate in the conscious condition (< 30s) with signal risetimes twice as fast as in the anesthetized state (60s). Both cortical and subcortical brain regions showed a robust, dose- dependent increase in MRI signal intensity with hypercapnic challenge while the animals were conscious but little or no change when anesthetized. Baseline variations in MRI signal were higher while animals were conscious but this was off set by greater signal intensity changes leading to a greater contrast-to-noise ratio, 13.1 in conscious animals, as compared to 8.0 in the anesthetized condition. In summary, cerebral vasculature appears to be more sensitive to hypercapnic challenge in the conscious condition resulting in enhanced T2* MRI signal intensity and the potential for better BOLD signal changes during functional imaging.     

Characterization of BOLD-fMRI signal during a verbal fluency paradigm in patients with intracerebral tumors affecting the frontal lobe

Previous studies have indicated that the BOLD-fMRI signal can be modified by tumor processes in close vicinity to functional brain areas. This effect has been investigated primarily for the perirolandic area but there is only a limited number of studies concerning frontal cortical regions. Therefore, the aim of the current study was to characterize BOLD-fMRI signal and activation patterns in patients with frontal brain tumors while performing a verbal fluency task. Six patients (ages 31-56 years) suffering from frontal (5 left sided and 1 right sided) intracerebral tumors were examined with fMRI while performing a verbal fluency task in a blocked paradigm design. Eight healthy volunteers served as the control group. The patients (5 right and 1 left handed) demonstrated left frontal activation which could be clearly located outside the tumor area and adjacent edema with varying degrees of additional right frontal activation. In the predominant left frontal activation cluster, the mean voxel based z-score and cluster size were not statistically different between patients and controls. The present fMRI study is indicating that language related BOLD signal changes in the frontal cortex of patients with tumors close to functional areas were comparable to the signal in normal controls. Additionally, the temporal hemodynamic response characteristic was comparable in both groups. This is an important finding consistent with PET results and corroborates the feasibility of functional mapping approaches in patients with tumors affecting the frontal lobe. Additional studies investigating alterations of the hemodynamic response depending on tumor location and histology are required in order to further elucidate the association between pathophysiology and BOLD fMRI signal.