Assessing the repeatability of absolute CMRO2, OEF and haemodynamic measurements from calibrated fMRI

As energy metabolism in the brain is largely oxidative, the measurement of cerebral metabolic rate of oxygen consumption (CMRO2) is a desirable biomarker for quantifying brain activity and tissue viability. Currently, PET techniques based on oxygen isotopes are the gold standard for obtaining whole brain CMRO2 maps. Among MRI techniques that have been developed as an alternative are dual calibrated fMRI (dcFMRI) methods, which exploit simultaneous measurements of BOLD and ASL signals during a hypercapnic-hyperoxic experiment to modulate brain blood flow and oxygenation. In this study we quantified the repeatability of a dcFMRI approach developed in our lab, evaluating its limits and informing its application in studies aimed at characterising the metabolic state of human brain tissue over time. Our analysis focussed on the estimates of oxygen extraction fraction (OEF), cerebral blood flow (CBF), CBF-related cerebrovascular reactivity (CVR) and CMRO2 based on a forward model that describes analytically the acquired dual echo GRE signal. Indices of within- and between-session repeatability are calculated from two different datasets both at a bulk grey matter and at a voxel-wise resolution and finally compared with similar indices obtained from previous MRI and PET measurements. Within- and between-session values of intra-subject coefficient of variation (CVintra) calculated from bulk grey matter estimates 6.7 ± 6.6% (mean ± std.) and 10.5 ± 9.7% for OEF, 6.9 ± 6% and 5.5 ± 4.7% for CBF, 12 ± 9.7% and 12.3 ± 10% for CMRO2. Coefficient of variation (CV) and intraclass correlation coefficient (ICC) maps showed the spatial distribution of the repeatability metrics, informing on the feasibility limits of the method. In conclusion, results show an overall consistency of the estimated physiological parameters with literature reports and a satisfactory level of repeatability considering the higher spatial sensitivity compared to other MRI methods, with varied performance depending on the specific parameter under analysis, on the spatial resolution considered and on the study design.

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 influence of noise on BOLD-mediated vessel size imaging analysis methods

Vessel size imaging (VSI) is a magnetic resonance imaging (MRI) technique that aims to provide quantitative measurements of tissue microvasculature. An emerging variation of this technique uses the blood oxygenation level-dependent (BOLD) effect as the source of the imaging contrast. Gas challenges have the advantage over contrast injection techniques in that they are noninvasive and easily repeatable because of the fast washout of the contrast. However, initial results from BOLD-VSI studies are somewhat contradictory, with substantially different estimates of the mean vessel radius. Owing to BOLD-VSI being an emerging technique, there is not yet a standard processing methodology, and different techniques have been used to calculate the mean vessel radius and reject uncertain estimates. In addition, the acquisition methodology and signal modeling vary from group to group. Owing to these differences, it is difficult to determine the source of this variation. Here we use computer modeling to assess the impact of noise on the accuracy and precision of different BOLD-VSI calculations. Our results show both potential overestimates and underestimates of the mean vessel radius, which is confirmed with a validation study at 3T.

An analysis of the use of hyperoxia for measuring venous cerebral blood volume: comparison of the existing method with a new analysis approach

Hyperoxia is known to cause an increase in the blood oxygenation level dependent (BOLD) signal that is primarily localised to the venous vasculature. This contrast mechanism has been proposed as a way to measure venous cerebral blood volume (CBVv) without the need for more invasive contrast media. In the existing method the analysis modelled the data as a dynamic contrast agent experiment, with the assumption that the BOLD signal of tissue was dominated by intravascular signal. The effects on the accuracy of the method due to extravascular BOLD signal changes, as well as signal modulation by intersubject differences in baseline physiology, such as haematocrit and oxygen extraction fraction, have so far been unexplored. In this study the effect of extravascular signal and intersubject physiological variability was investigated by simulating the hyperoxia CBVv experiment using a detailed BOLD signal model. This analysis revealed substantial uncertainty in the measurement of CBVv using the existing analysis based on dynamic contrast agent experiments. Instead, the modelling showed a simple and direct relationship between the BOLD signal change and CBVv, and an alternative analysis method with much reduced uncertainty was proposed based on this finding. Both methods were tested experimentally, with the new method producing results that are consistent with the limited literature in this area.