How much luxury is there in 'luxury perfusion'? An analysis of the BOLD response in the visual areas V1 and V2

We re-analyzed the functional magnetic resonance imaging data from a study involving awake, adult, human volunteers in order to examine the influence of vascular density on the blood oxygenation level-dependent (BOLD) response. We employed a flashed and reversing stimulus paradigm where the latter stimulated twice the number of receptive fields and with it doubled the neuronal metabolic load (CMRO2) compared to the former stimulus. The blood flow increase to these stimuli was identical, so that differences in the BOLD response are due to differences in the oxygen extraction fraction. By comparing the BOLD response in human striate cortex (V1) and its neighbor, extra-striate area V2 to the two stimuli, we were able to determine the influence of the higher vascular density of striate cortex on the BOLD response. In striate cortex, the extent of activation, as measured by the number of activated voxels, was larger for the flashed than for the reversing stimulus. In extra-striate area V2, no such difference in the extent of activation was noted. Gauging the local concentration of HbR using deltaR2*, we found it to be significantly lower for the flashed than for the reversing checkerboard. We estimated the HbR concentration in extra-striate area V2 to be double that of striate cortex independent of the stimulus presented. A frequency distribution of the deltaR2* values for the flashed and reversing checkerboard revealed a shift consistent with an increase in the HbR concentration between areas V1 and V2. The metabolically most demanding stimulus, the reversing checkerboard was associated with the highest HbR concentration and with the largest number of voxels with a negative BOLD response.     

When more means less: a paradox BOLD response in human visual cortex

The predictions of the 'Linear Transfer Model' (LTM) have been tested only by modulating the frequency of the action potentials while keeping the size of the activated neuronal population constant. The LTM states that the blood oxygenation level-dependent contrast (BOLD) signal is directly proportional to the neuronal activity averaged over milliseconds or seconds. We examined the influence on the BOLD response, of manipulating the size of the activated neuronal population while maintaining the electrical discharge activity constant. We performed functional MR measurements on 30 awake, healthy adult volunteers (15 male and 15 female) using a flashed and reversing checkerboard. These stimuli induced the same vascular response and the same increase in the electrical discharge activity but varied in the size of the neuronal population being activated. The BOLD response measured by the extent of activation and the BOLD signal amplitude, was larger for the flashed than to the reversing checkerboard. An assessment of the local deoxyhemoglobin (HbR) concentration indicated that the neuronal activity was lower during the flashed checkerboard than the reversing checkerboard. Because the checkerboard associated with the lower neuronal activity yielded the larger number of activated voxels and the larger BOLD signal, our results run contrary to the predictions of the 'Linear Transfer Model' and for this reason we refer to them as paradoxical. Stimuli defined by luminance contrast or a chromatic contrast yielded identical results. We conclude that the 'LTM' may apply to stimuli that modulate the electrical discharge activity but not to stimuli that modulate the size of the activated neuronal population.     

Functional MRI of the human motor cortex using single-shot, multiple gradient-echo spiral imaging

In this study, we combined the advantages of a fast multi-slice spiral imaging approach with a multiple gradient-echo sampling scheme at high magnetic field strength to improve quantification of BOLD and inflow effects and to estimate T2* relaxation times in functional brain imaging. Eight echoes are collected with echo time (TE) ranging from 5 to 180 ms. Acquisition time per slice and echo time is 25 ms for a nominal resolution of 4 x 4 x 4 mm3. Evaluation of parameter images during rest and stimulation yields no significant activation on the inflow sensitive spin-density images (rho or I0-maps) whereas clear activation patterns in primary human motor cortex (M1) and supplementary motor area (SMA) are detected on BOLD sensitive T2*-maps. The calculation of relaxation times and rates of the activated areas over all subjects yields an average T2* +/- standard deviation (SD) of 46.1+/-4.5 ms (R2* of 21.8+/-2.2 s(-1)) and an average increase (deltaT2* +/- SD) of 0.93+/-0.47 ms (deltaR2* of -0.4+/-0.14 s(-1)). Our findings demonstrate the usefulness of a multiple gradient echo data acquisition approach in separating various vascular contributions to brain activation in fMRI.     

Quantification of venous blood signal contribution to BOLD functional activation in the auditory cortex at 3 T

Most modern techniques for functional magnetic resonance imaging (fMRI) rely on blood-oxygen-level-dependent (BOLD) contrast as the basic principle for detecting neuronal activation. However, the measured BOLD effect depends on a transfer function related to neurophysiological changes accompanying electrical neural activation. The spatial accuracy and extension of the region of interest are determined by vascular effect, which introduces incertitude on real neuronal activation maps. Our efforts have been directed towards the development of a new methodology that is capable of combining morphological, vascular and functional information; obtaining new insight regarding foci of activation; and distinguishing the nature of activation on a pixel-by-pixel basis. Six healthy volunteers were studied in a parametric auditory functional experiment at 3 T; activation maps were overlaid on a high-resolution brain venography obtained through a novel technique. The BOLD signal intensities of vascular and nonvascular activated voxels were analyzed and compared: it was shown that nonvascular active voxels have lower values for signal peak (P<10(-7)) and area (P<10(-8)) with respect to vascular voxels. The analysis showed how venous blood influenced the measured BOLD signals, supplying a technique to filter possible venous artifacts that potentially can lead to misinterpretation of fMRI results. This methodology, although validated in the auditory cortex activation, maintains a general applicability to any cortical fMRI study, as the basic concepts on which it relies on are not limited to this cortical region. The results obtained in this study can represent the basis for new methodologies and tools that are capable of adding further characterization to the BOLD signal properties.     

Non-neural BOLD variability in block and event-related paradigms

Block and event-related stimulus designs are typically used in fMRI studies depending on the importance of detection power or estimation efficiency. The extent of vascular contribution to variability in block and event-related fMRI-BOLD response is not known. With scaling, the extent of vascular variability in the fMRI-BOLD response during block and event-related design tasks was investigated. Blood oxygen level-dependent (BOLD) contrast data from healthy volunteers performing a block design motor task and an event-related memory task requiring performance of a motor response were analyzed from the regions of interest (ROIs) surrounding the primary and supplementary motor cortices. Average BOLD signal change was significantly larger during the block design compared to the event-related design. In each subject, BOLD signal change across voxels in the ROIs had higher variation during the block design task compared to the event-related design task. Scaling using the resting state fluctuation of amplitude (RSFA) and breath-hold (BH), which minimizes BOLD variation due to vascular origins, reduced the within-subject BOLD variability in every subject during both tasks but significantly reduced BOLD variability across subjects only during the block design task. The strong non-neural source of intra- and intersubject variability of BOLD response during the block design compared to event-related task indicates that study designs optimizing for statistical power through enhancement of the BOLD contrast (for, e.g., block design) can be affected by enhancement of non-neural sources of BOLD variability.     

Influence of speech stimuli intensity on the activation of auditory cortex investigated with functional magnetic resonance imaging.

Understanding the impact of variations in the acoustic signal is critical for the development of auditory and language fMRI as an experimental tool. We describe the dependence of the BOLD signal and speech intelligibility on the intensity of auditory stimuli. Eighteen subjects were imaged on a 1.5-T MRI scanner. Speech stimuli were English monosyllabic words played at five intensity levels. Intrasubject reproducibility was measured on one subject by presenting the stimulus five times at the same intensity level. Intelligibility was measured during data acquisition as subjects signaled when hearing two targets. Each functional trial consisted of four cycles (30 s off-30 s on). Five oblique slices covering primary and association auditory areas were imaged. Activated voxels were identified by cross-correlation analysis and their percent signal change (delta S) was measured. Intersubject differences in activation extent, asymmetry, and dependence on intensity were striking. Volume of activation was significantly greater in the left than in the right hemisphere. Intrasubject reproducibility for delta S was higher than for volume of activation. delta S and intelligibility showed a similar dependence on intensity suggesting that not only intensity but also intelligibility affect the fMRI signal.     

Component structure of event-related fMRI responses in the different neurovascular compartments

In most functional magnetic resonance imaging (fMRI) studies, brain activity is localized by observing changes in the blood oxygenation level-dependent (BOLD) signal that are believed to arise from capillaries, venules and veins in and around the active neuronal population. However, the contribution from veins can be relatively far downstream from active neurons, thereby limiting the ability of BOLD imaging methods to precisely pinpoint neural generators. Hemodynamic measures based on apparent diffusion coefficients (ADCs) have recently been used to identify more upstream functional blood flow changes in the capillaries, arterioles and arteries. In particular, we recently showed that, due to the complementary vascular sensitivities of ADC and BOLD signals, the voxels conjointly activated by both measures may identify the capillary networks of the active neuronal areas. In this study, we first used simultaneously acquired ADC and BOLD functional imaging signals to identify brain voxels activated by ADC only, by both ADC and BOLD and by BOLD only, thereby delineating voxels relatively dominated by the arterial, capillary, and draining venous neurovascular compartments, respectively. We then examined the event-related fMRI BOLD responses in each of these delineated neurovascular compartments, hypothesizing that their event-related responses would show different temporal componentries. In the regions activated by both the BOLD and ADC contrasts, but not in the BOLD-only areas, we observed an initial transient signal reduction (an initial dip), consistent with the local production of deoxyhemoglobin by the active neuronal population. In addition, the BOLD-ADC overlap areas and the BOLD-only areas showed a clear poststimulus undershoot, whereas the compartment activated by only ADC did not show this component. These results indicate that using ADC contrast in conjunction with BOLD imaging can help delineate the various neurovascular compartments, improve the localization of active neural populations, and provide insight into the physiological mechanisms underlying the hemodynamic signals.     

Effects of mild hypoxic hypoxia on poststimulus undershoot of blood-oxygenation-level-dependent fMRI signal in the human visual cortex

Characteristics of the blood-oxygenation-level-dependent (BOLD) functional magnetic resonance imaging (fMRI) signal poststimulus undershoot in the visual cortex were studied at varying levels of arterial blood oxygen saturation (Ysat). Undershoot with an amplitude of -0.6+/-0.2% appeared after positive BOLD response (+1.7+/-0.5%) under control conditions. Cerebral blood volume (CBV), as determined with vascular-space-occupancy-dependent fMRI, increased by 26-43% during the positive BOLD peak, but the CBV proceeded at baseline level during the BOLD poststimulus undershoot. Mild hypoxic hypoxia (Ysat ranging from 0.82 to 0.89) had no effect on the amplitude or duration of poststimulus undershoot in activated BOLD pixels. Hypoxia did not influence CBV during the BOLD poststimulus undershoot. In contrast, the positive BOLD signal at the level of all activated pixels was smaller in hypoxia than in normoxia. The present results show that the BOLD poststimulus undershoot is not influenced by curtailed oxygen availability and that, during the undershoot, CBV is not different from resting state.     

Quantitative NumART2* mapping in functional MRI studies at 1.5 T

Quantitative mapping of the effective transverse relaxation time, T₂* and proton density was performed in a motor activation functional MRI (fMRI) study using multi-echo, echo planar imaging (EPI) and NumART₂* (Numerical Algorithm for Real time T₂*). Comparisons between NumART₂* and conventional single echo EPI with an echo time of 64 ms were performed for five healthy participants examined twice. Simulations were also performed to address specific issues associated with the two techniques, such as echo time-dependent signal variation. While the single echo contrast varied with the baseline T₂* value, relative changes in T₂* remained unaffected. Statistical analysis of the T₂* maps yielded fMRI activation patterns with an improved statistical detection relative to conventional EPI but with less activated voxels, suggesting that NumART₂* has superior spatial specificity. Two effects, inflow and dephasing, that may explain this finding were investigated. Particularly, a statistically significant increase in proton density was found in a brain area that was detected as activated by conventional EPI but not by NumART₂* while no such changes were observed in brain areas that showed stimulus correlated signal changes on T₂* maps.     

Influence of baseline hematocrit and hemodilution on BOLD fMRI activation.

Current understanding of blood oxygenation level dependent (BOLD) fMRI physiology predicts a close relationship between BOLD signal and blood hematocrit level. However, neither this relationship nor its effect on BOLD percent activation (BPA) has been empirically examined in man. To that end, BPA in primary visual cortex in response to photic stimulation was determined in a group of 24 normal subjects. A positive linear relationship between BPA and hematocrit was seen, particularly in men. To evaluate the effect of change in hematocrit on BPA, 9 men were studied before and following isotonic saline hemodilution, resulting in an average 6% reduction in hematocrit and an 8-31% reduction in BPA. No significant change in the number of activated pixels was seen. A model of predicted BPA as a function of hematocrit and vessel size was developed, and results from this model closely mirrored the empiric data. These results suggest that hematocrit significantly influences the magnitude of BPA and that such baseline factors should be accounted for when comparing BOLD data across groups of subjects, particularly in the many instances in which hematocrit may vary systematically. Such instances include several disease states as well as studies involving sex differences, drug administration, stress and other factors. Finally, the robust agreement between predicted and empiric data serves to validate a semiquantitative approach to the analysis of BOLD fMRI data.     

Simultaneous BOLD/perfusion measurement using dual-echo FAIR and UNFAIR: sequence comparison at 1.5T and 3.0T.

Functional MRI (fMRI) studies designed for simultaneously measuring Blood Oxygenation Level Dependent (BOLD) and Cerebral Blood Flow (CBF) signal often employ the standard Flow Alternating Inversion Recovery (FAIR) technique. However, some sensitivity is lost in the BOLD data due to inherent T1 relaxation. We sought to minimize the preceding problem by employing a modified UN-inverted FAIR (UNFAIR) technique, which (in theory) should provide identical CBF signal as FAIR with minimal degradation of the BOLD signal. UNFAIR BOLD maps acquired from human subjects (n = 8) showed significantly higher mean z-score of approximately 17% (p < 0.001), and number of activated voxels at 1.5T. On the other hand, the corresponding FAIR perfusion maps were superior to the UNFAIR perfusion maps as reflected in a higher mean z-score of approximately 8% (p = 0.013), and number of activated voxels. The reduction in UNFAIR sensitivity for perfusion is attributed to increased motion sensitivity related to its higher background signal, and, T2 related losses from the use of an extra inversion pulse. Data acquired at 3.0T demonstrating similar trends are also presented.     

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.     

Temporal and spatial MRI responses to subsecond visual activation

The temporal and spatial characteristics of oxygenation-sensitive MRI responses to very brief visual stimuli (five Hz reversing black and white checkerboard pattern versus darkness) were investigated (nine subjects) by means of serial single-shot gradient-echo echo-planar imaging (2.0 T, TR = 400 ms, mean TE = 54 ms, flip angle 30°). The use of a 0.2-s stimulus and a 90-s control phase resulted in an initial latency phase (about 2 s, no signal change), a positive MRI response (2.5% signal increase peaking at 5 s after stimulus onset), and a post-stimulus undershoot (1% signal decrease peaking at 15 s after stimulus onset) lasting for about 50–60 s. The finding that a subsecond visual stimulus elicits both a strong positive MRI response and a long-lasting undershoot provides further evidence for the neuronal origin of slow signal fluctuations seen in the absence of functional challenge and their utility for mapping functional connectivity. The additional observation that a reduction of the inter-stimulus control phase from 90 s to 9.8 s does not seem to affect the spatial extent of cortical activation in pertinent maps is of major relevance for the design and analysis of “event-related” MRI studies.     

The quantification of DeltaR2(*) under brain activation: dependence on relaxation rate at rest and significance threshold.

Transverse relaxation rates R2(*) were measured in subjects performing a motor task using a segmented EPI double gradient echo sequence (TE = 23/70 ms) with five different voxel sizes between 1.8 mm(3) and 41.7 mm(3). An analysis of the errors involved in the calculation of the change of the transverse relaxation rate--DeltaR2(*) and of the consequences of defining an arbitrary threshold of statistical significance in the data analysis was performed. Correlations between the magnitude of the BOLD effect and the significance level on one hand and between the transverse relaxation time at rest and its change under activation on the other, both referenced in the literature, can be understood as a consequence of this procedure. Analysing histograms of parameter changes rather than average values alone allows for an estimate of the contribution of false positive voxels. Furthermore, while the averaged signal change increases in proportion to the selection threshold the histograms of activated voxels remain insensitive to the latter.     

Functional magnetic resonance imaging with intermolecular multiple-quantum coherences

For the first time, we demonstrate here functional magnetic resonance imaging (fMRI) using intermolecular multiple-quantum coherences (iMQCs). iMQCs are normally not observed in liquid-state NMR because dipolar interactions between spins average to zero. If the magnetic isotropy of the sample is broken through the use of magnetic field gradients, dipolar couplings can reappear, and hence iMQCs can be observed. Conventional (BOLD) fMRI measures susceptibility variations averaged over each voxel. In the experiment performed here, the sensitivity of iMQCs to frequency variations over mesoscopic and well-defined distances is exploited. We show that iMQC contrast is qualitatively and quantitatively different from BOLD contrast in a visual stimulation task. While the number of activated pixels is smaller in iMQC contrast, the intensity change in some pixels exceeds that of BOLD contrast severalfold.     

Input permutation method to detect active voxels in fMRI study

Correctly identifying voxels or regions of interest (ROI) that actively respond to a given stimulus is often an important objective/step in many functional magnetic resonance imaging (fMRI) studies. In this article, we study a nonparametric method to detect active voxels, which makes minimal assumption about the distribution of blood oxygen level-dependent (BOLD) signals. Our proposal has several interesting features. It uses time lagged correlation to take into account the delay in response to the stimulus, due to hemodynamic variations. We introduce an input permutation method (IPM), a type of block permutation method, to approximate the null distribution of the test statistic. Also, we propose to pool the permutation-derived statistics of preselected voxels for a better approximation to the null distribution. Finally, we control multiple testing error rate using the local false discovery rate (FDR) by Efron [Correlation and large-scale simultaneous hypothesis testing. J Am Stat Assoc 102 (2007) 93–103] and Park et al. [Estimation of empirical null using a mixture of normals and its use in local false discovery rate. Comput Stat Data Anal 55 (2011) 2421–2432] to select the active voxels.     

A comparison of Gamma and Gaussian dynamic convolution models of the fMRI BOLD response

Blood oxygenation level-dependent (BOLD) contrast-based functional magnetic resonance imaging (fMRI) has been widely utilized to detect brain neural activities and great efforts are now stressed on the hemodynamic processes of different brain regions activated by a stimulus. The focus of this paper is the comparison of Gamma and Gaussian dynamic convolution models of the fMRI BOLD response. The convolutions are between the perfusion function of the neural response to a stimulus and a Gaussian or Gamma function. The parameters of the two models are estimated by a nonlinear least-squares optimal algorithm for the fMRI data of eight subjects collected in a visual stimulus experiment. The results show that the Gaussian model is better than the Gamma model in fitting the data. The model parameters are different in the left and right occipital regions, which indicate that the dynamic processes seem different in various cerebral functional regions.     

The effect of task block arrangement on the detectability of activation in fMRI

The effect of task block arrangements on the detection of brain activation was investigated. Sessions of functional magnetic resonance imaging (fMRI) including the same number of two different task conditions but with different arrangements were compared. The two task conditions were, A) Ellipse-shaped black and white checkerboard flicker stimulation at 4.2 Hz covering the bilateral visual field, and B) the same flicker stimuli covering only the left visual field. In the rest blocks (0), the subjects looked at a fixation point. Four different task block arrangements were compared, 1) A0 (0A0A0A0) and B0 (0B0B0B0), 2) A0B0 (0A0B0A0B0A0B0), 3) AB0 (0AB0AB0AB0) and 4) AB (0ABABAB). Bilateral V1, V2, V3 and the left V5 were activated by condition A, and the right V1 and V2 by B. The activation in the left visual field by A0 was larger than in the other three conditions. In a differential analysis between conditions A and B, activation in the left V3 and V5 was declined by AB0 or AB. When rest blocks were located in the post-stimulus undershoot phase, the % signal change of the BOLD signal was emphasized, which caused augmented significance in the detection of the activity. It was indicated that the outcome of the activation map was influenced by the arrangement of task blocks, even though the same number of task blocks were repeated within the sessions. In fMRI studies, task conditions should be carefully compared within or across sessions considering the characteristics of hemodynamic response functions.     

Increased BOLD sensitivity in the orbitofrontal cortex using slice-dependent echo times at 3 T

Functional magnetic resonance imaging (fMRI) exploits the blood oxygenation level dependent (BOLD) effect to detect neuronal activation related to various experimental paradigms. Some of these, such as reversal learning, involve the orbitofrontal cortex and its interaction with other brain regions like the amygdala, striatum or dorsolateral prefrontal cortex. These paradigms are commonly investigated with event-related methods and gradient echo-planar imaging (EPI) with short echo time of 27 ms. However, susceptibility-induced signal losses and image distortions in the orbitofrontal cortex are still a problem for this optimized sequence as this brain region consists of several slices with different optimal echo times. An EPI sequence with slice-dependent echo times is suitable to maximize BOLD sensitivity in all slices and might thus improve signal detection in the orbitofrontal cortex. To test this hypothesis, we first optimized echo times via BOLD sensitivity simulation. Second, we measured 12 healthy volunteers using a standard EPI sequence with an echo time of 27 ms and a modified EPI sequence with echo times ranging from 22 ms to 47 ms. In the orbitofrontal cortex, the number of activated voxels increased from 87±44 to 549±83 and the maximal t-value increased from 4.4±0.3 to 5.4±0.3 when the modified EPI was used. We conclude that an EPI with slice-dependent echo times may be a valuable tool to mitigate susceptibility artifacts in event-related whole-brain fMRI studies with a focus on the orbitofrontal cortex.     

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.