Validation of a new 3D quantitative BOLD based cerebral oxygen extraction mapping

Quantitative BOLD (qBOLD) MRI allows evaluation of oxidative metabolism of the brain based purely on an endogenous contrast mechanism. The method quantifies deoxygenated blood volume (DBV) and hemoglobin oxygen saturation level of venous blood (Yv), yielding oxygen extraction fraction (OEF), and along with a separate measurement of cerebral blood flow, cerebral metabolic rate of oxygen (CMRO2) maps. Here, we evaluated our recently reported 3D qBOLD method that rectifies a number of deficiencies in prior qBOLD approaches in terms of repeat reproducibility and sensitivity to hypercapnia on the metabolic parameters, and in comparison to dual-gas calibrated BOLD (cBOLD) MRI for determining resting-state oxygen metabolism. Results suggested no significant difference between test-retest qBOLD scans in either DBV and OEF. Exposure to hypercapnia yielded group averages of 38 and 28% for OEF and 151 and 146 µmol/min/100 g for CMRO2 in gray matter at baseline and hypercapnia, respectively. The decrease of OEF during hypercapnia was significant (p ≪ 0.01), whereas CMRO2 did not change significantly (p = 0.25). Finally, baseline OEF (37 vs. 39%) and CMRO2 (153 vs. 145 µmol/min/100 g) in gray matter using qBOLD and dual-gas cBOLD were found to be in good agreement with literature values, and were not significantly different from each other (p > 0.1).

Field strength dependence of grey matter R2* on venous oxygenation

The relationship between venous blood oxygenation and change in transverse relaxation rate (ΔR2*) plays a key role in calibrated BOLD fMRI. This relationship, defined by the parameter β, has previously been determined using theoretical simulations and experimental measures. However, these earlier studies have been confounded by the change in venous cerebral blood volume (CBV) in response to functional tasks. This study used a double-echo gradient echo EPI scheme in conjunction with a graded isocapnic hyperoxic sequence to assess quantitatively the relationship between the fractional venous blood oxygenation (1−Y_v) and transverse relaxation rate of grey matter (ΔR2⁢GM*), without inducing a change in vCBV.

The results demonstrate that the relationship between ΔR2* and fractional venous oxygenation at all magnet field strengths studied was adequately described by a linear relationship. The gradient of this relationship did not increase monotonically with field strength, which may be attributed to the relative contributions of intravascular and extravascular signals which will vary with both field strength and blood oxygenation.

•

We assess the relationship between grey matter R2* and venous oxygenation.

•

Isocapnic hyperoxia prevented confounding changes in cerebral blood volume.

•

A linear dependency is an appropriate assumption at 1.5, 3 and 7 T.

•

Intravascular/extravascular signal ratios will vary with both B0 and oxygenation.

Calibrated bold fMRI with an optimized ASL-BOLD dual-acquisition sequence

Calibrated fMRI techniques estimate task-induced changes in the cerebral metabolic rate of oxygen (CMRO₂) based on simultaneous measurements of cerebral blood flow (CBF) and blood-oxygen-level-dependent (BOLD) signal changes evoked by stimulation. To determine the calibration factor M (corresponding to the maximum possible BOLD signal increase), BOLD signal and CBF are measured in response to a gas breathing challenge (usually CO₂ or O₂). Here we describe an ASL dual-acquisition sequence that combines a background-suppressed 3D-GRASE readout with 2D multi-slice EPI. The concatenation of these two imaging sequences allowed separate optimization of the acquisition for CBF and BOLD data. The dual-acquisition sequence was validated by comparison to an ASL sequence with a dual-echo EPI readout, using a visual fMRI paradigm. Results showed a 3-fold increase in temporal signal-to-noise ratio (tSNR) of the ASL time-series data while BOLD tSNR was similar to that obtained with the dual-echo sequence. The longer TR of the proposed dual-acquisition sequence, however, resulted in slightly lower T-scores (by 30%) in the BOLD activation maps. Further, the potential of the dual-acquisition sequence for M-mapping on the basis of a hypercapnia gas breathing challenge and for quantification of CMRO₂ changes in response to a motor activation task was assessed. In five subjects, an average gray matter M-value of 8.71±1.03 and fractional changes of CMRO₂ of 12.5±5% were found. The new sequence remedies the deficiencies of prior combined BOLD-ASL acquisition strategies by substantially enhancing perfusion tSNR, which is essential for accurate BOLD calibration.

Measurement of oxygen extraction fraction (OEF): An optimized BOLD signal model for use with hypercapnic and hyperoxic calibration

Several techniques have been proposed to estimate relative changes in cerebral metabolic rate of oxygen consumption (CMRO2) by exploiting combined BOLD fMRI and cerebral blood flow data in conjunction with hypercapnic or hyperoxic respiratory challenges. More recently, methods based on respiratory challenges that include both hypercapnia and hyperoxia have been developed to assess absolute CMRO2, an important parameter for understanding brain energetics. In this paper, we empirically optimize a previously presented "original calibration model" relating BOLD and blood flow signals specifically for the estimation of oxygen extraction fraction (OEF) and absolute CMRO2. To do so, we have created a set of synthetic BOLD signals using a detailed BOLD signal model to reproduce experiments incorporating hypercapnic and hyperoxic respiratory challenges at 3T. A wide range of physiological conditions was simulated by varying input parameter values (baseline cerebral blood volume (CBV0), baseline cerebral blood flow (CBF0), baseline oxygen extraction fraction (OEF0) and hematocrit (Hct)). From the optimization of the calibration model for estimation of OEF and practical considerations of hypercapnic and hyperoxic respiratory challenges, a new "simplified calibration model" is established which reduces the complexity of the original calibration model by substituting the standard parameters α and β with a single parameter θ. The optimal value of θ is determined (θ=0.06) across a range of experimental respiratory challenges. The simplified calibration model gives estimates of OEF0 and absolute CMRO2 closer to the true values used to simulate the experimental data compared to those estimated using the original model incorporating literature values of α and β. Finally, an error propagation analysis demonstrates the susceptibility of the original and simplified calibration models to measurement errors and potential violations in the underlying assumptions of isometabolism. We conclude that using the simplified calibration model results in a reduced bias in OEF0 estimates across a wide range of potential respiratory challenge experimental designs.

The absolute CBF response to activation is preserved during elevated perfusion: Implications for neurovascular coupling measures

Functional magnetic resonance imaging (fMRI) techniques in which the blood oxygenation level dependent (BOLD) and cerebral blood flow (CBF) response to a neural stimulus are measured, can be used to estimate the fractional increase in the cerebral metabolic rate of oxygen consumption (CMRO2) that accompanies evoked neural activity. A measure of neurovascular coupling is obtained from the ratio of fractional CBF and CMRO2 responses, defined as n, with the implicit assumption that relative rather than absolute changes in CBF and CMRO2 adequately characterise the flow-metabolism response to neural activity. The coupling parameter n is important in terms of its effect on the BOLD response, and as potential insight into the flow-metabolism relationship in both normal and pathological brain function. In 10 healthy human subjects, BOLD and CBF responses were measured to test the effect of baseline perfusion (modulated by a hypercapnia challenge) on the coupling parameter n during graded visual stimulation. A dual-echo pulsed arterial spin labelling (PASL) sequence provided absolute quantification of CBF in baseline and active states as well as relative BOLD signal changes, which were used to estimate CMRO2 responses to the graded visual stimulus. The absolute CBF response to the visual stimuli were constant across different baseline CBF levels, meaning the fractional CBF responses were reduced at the hyperperfused baseline state. For the graded visual stimuli, values of n were significantly reduced during hypercapnia induced hyperperfusion. Assuming the evoked neural responses to the visual stimuli are the same for both baseline CBF states, this result has implications for fMRI studies that aim to measure neurovascular coupling using relative changes in CBF. The coupling parameter n is sensitive to baseline CBF, which would confound its interpretation in fMRI studies where there may be significant differences in baseline perfusion between groups. The absolute change in CBF, as opposed to the change relative to baseline, may more closely match the underlying increase in neural activity in response to a stimulus.

A novel Bayesian approach to accounting for uncertainty in fMRI-derived estimates of cerebral oxygen metabolism fluctuations

Calibrated blood oxygenation level dependent (BOLD) imaging is a multimodal functional MRI technique designed to estimate changes in cerebral oxygen metabolism from measured changes in cerebral blood flow and the BOLD signal. This technique addresses fundamental ambiguities associated with quantitative BOLD signal analysis; however, its dependence on biophysical modeling creates uncertainty in the resulting oxygen metabolism estimates. In this work, we developed a Bayesian approach to estimating the oxygen metabolism response to a neural stimulus and used it to examine the uncertainty that arises in calibrated BOLD estimation due to the presence of unmeasured model parameters. We applied our approach to estimate the CMRO2 response to a visual task using the traditional hypercapnia calibration experiment as well as to estimate the metabolic response to both a visual task and hypercapnia using the measurement of baseline apparent R2' as a calibration technique. Further, in order to examine the effects of cerebral spinal fluid (CSF) signal contamination on the measurement of apparent R2', we examined the effects of measuring this parameter with and without CSF-nulling. We found that the two calibration techniques provided consistent estimates of the metabolic response on average, with a median R2'-based estimate of the metabolic response to CO2 of 1.4%, and R2'- and hypercapnia-calibrated estimates of the visual response of 27% and 24%, respectively. However, these estimates were sensitive to different sources of estimation uncertainty. The R2'-calibrated estimate was highly sensitive to CSF contamination and to uncertainty in unmeasured model parameters describing flow-volume coupling, capillary bed characteristics, and the iso-susceptibility saturation of blood. The hypercapnia-calibrated estimate was relatively insensitive to these parameters but highly sensitive to the assumed metabolic response to CO2.

Sources of systematic error in calibrated BOLD based mapping of baseline oxygen extraction fraction

Recently a new class of calibrated blood oxygen level dependent (BOLD) functional magnetic resonance imaging (MRI) methods were introduced to quantitatively measure the baseline oxygen extraction fraction (OEF). These methods rely on two respiratory challenges and a mathematical model of the resultant changes in the BOLD functional MRI signal to estimate the OEF. However, this mathematical model does not include all of the effects that contribute to the BOLD signal, it relies on several physiological assumptions and it may be affected by intersubject physiological variability. The aim of this study was to investigate these sources of systematic error and their effect on estimating the OEF. This was achieved through simulation using a detailed model of the BOLD signal. Large ranges for intersubject variability in baseline physiological parameters such as haematocrit and cerebral blood volume were considered. Despite this the uncertainty in the relationship between the measured BOLD signals and the OEF was relatively low. Investigations of the physiological assumptions that underlie the mathematical model revealed that OEF measurements are likely to be overestimated if oxygen metabolism changes during hypercapnia or cerebral blood flow changes under hyperoxia. Hypoxic hypoxia was predicted to result in an underestimation of the OEF, whilst anaemic hypoxia was found to have only a minimal effect.

Calibrating the BOLD response without administering gases: comparison of hypercapnia calibration with calibration using an asymmetric spin echo

The calibrated BOLD (blood oxygen level dependent) technique was developed to quantify the BOLD signal in terms of changes in oxygen metabolism. In order to achieve this a calibration experiment must be performed, which typically requires a hypercapnic gas mixture to be administered to the participant. However, an emerging technique seeks to perform this calibration without administering gases using a refocussing based calibration. Whilst hypercapnia calibration seeks to emulate the physical removal of deoxyhaemoglobin from the blood, the aim of refocussing based calibration is to refocus the dephasing effect of deoxyhaemoglobin on the MR signal using a spin echo. However, it is not possible to refocus all of the effects that contribute to the BOLD signal and a scale factor is required to estimate the BOLD scaling parameter M. In this study the feasibility of a refocussing based calibration was investigated. The scale factor relating the refocussing calibration to M was predicted by simulations to be approximately linear and empirically measured to be 0.88±0.36 for the visual cortex and 0.93±0.32 for a grey matter region of interest (mean±standard deviation). Refocussing based calibration is a promising approach for greatly simplifying the calibrated BOLD methodology by eliminating the need for the subject to breathe special gas mixtures, and potentially provides the basis for a wider implementation of quantitative functional MRI.

The coupling of cerebral blood flow and oxygen metabolism with brain activation is similar for simple and complex stimuli in human primary visual cortex

Quantitative functional MRI (fMRI) experiments to measure blood flow and oxygen metabolism coupling in the brain typically rely on simple repetitive stimuli. Here we compared such stimuli with a more naturalistic stimulus. Previous work on the primary visual cortex showed that direct attentional modulation evokes a blood flow (CBF) response with a relatively large oxygen metabolism (CMRO2) response in comparison to an unattended stimulus, which evokes a much smaller metabolic response relative to the flow response. We hypothesized that a similar effect would be associated with a more engaging stimulus, and tested this by measuring the primary human visual cortex response to two contrast levels of a radial flickering checkerboard in comparison to the response to free viewing of brief movie clips. We did not find a significant difference in the blood flow-metabolism coupling (n=%ΔCBF/%ΔCMRO2) between the movie stimulus and the flickering checkerboards employing two different analysis methods: a standard analysis using the Davis model and a new analysis using a heuristic model dependent only on measured quantities. This finding suggests that in the primary visual cortex a naturalistic stimulus (in comparison to a simple repetitive stimulus) is either not sufficient to provoke a change in flow-metabolism coupling by attentional modulation as hypothesized, that the experimental design disrupted the cognitive processes underlying the response to a more natural stimulus, or that the technique used is not sensitive enough to detect a small difference.

Effects of background suppression on the sensitivity of dual-echo arterial spin labeling MRI for BOLD and CBF signal changes

Dual-echo arterial spin labeling (DE-ASL) enables the simultaneous acquisition of BOLD and CBF fMRI data and is often used for calibrated BOLD and cerebrovascular CO2 reactivity measurements. DE-ASL, like all ASL techniques, suffers from a low intrinsic CBF SNR, which can be improved by suppressing the background signal via the inclusion of additional inversion pulses. However, until now this approach has been considered to be undesirable for DE-ASL, because the BOLD signal is extracted from the background signal and attenuating the background signal could decrease the sensitivity of DE-ASL scans for BOLD changes. In this study, the effect of background suppression on the sensitivity of DE-ASL MRI for BOLD and CBF signal changes with a visual stimulation paradigm was studied. Results showed that with an average background suppression level of 70% the BOLD sensitivity of DE-ASL MRI decreases slightly (15%), while the CBF sensitivity of the scans increased by almost a factor-of-two (81%). These findings support the conclusion that the gains in CBF sensitivity of DE-ASL MRI due to background suppression outweigh the slight decrease in sensitivity of these scans for BOLD changes, and thus that background suppression is highly recommended for DE-ASL.