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

Evidence for direct CO2 -mediated alterations in cerebral oxidative metabolism in humans

AIM

How the cerebral metabolic rates of oxygen and glucose utilization (CMRO2 and CMRGlc, respectively) are affected by alterations in arterial PCO2 (PaCO2) is equivocal and therefore was the primary question of this study.

METHODS

This retrospective analysis involved pooled data from four separate studies, involving 41 healthy adults (35 males/6 females). Participants completed stepwise steady-state alterations in PaCO2 ranging between 30 and 60 mmHg. The CMRO2 and CMRGlc were assessed via the Fick approach (CBF × arterial-internal jugular venous difference of oxygen or glucose content, respectively) utilizing duplex ultrasound of the internal carotid artery and vertebral artery to calculate cerebral blood flow (CBF).

RESULTS

The CMRO2 was altered by 0.5 mL × min-1 (95% CI: -0.6 to -0.3) per mmHg change in PaCO2 (p < 0.001) which corresponded to a 9.8% (95% CI: -13.2 to -6.5) change in CMRO2 with a 9 mmHg change in PaCO2 (inclusive of hypo- and hypercapnia). The CMRGlc was reduced by 7.7% (95% CI: -15.4 to -0.08, p = 0.045; i.e., reduction in net glucose uptake) and the oxidative glucose index (ratio of oxygen to glucose uptake) was reduced by 5.6% (95% CI: -11.2 to 0.06, p = 0.049) with a + 9 mmHg increase in PaCO2.

CONCLUSION

Collectively, the CMRO2 is altered by approximately 1% per mmHg change in PaCO2. Further, glucose is incompletely oxidized during hypercapnia, indicating reductions in CMRO2 are either met by compensatory increases in nonoxidative glucose metabolism or explained by a reduction in total energy production.

Anesthesia Depresses Cerebrovascular Reactivity to Acetazolamide in Pediatric Moyamoya Vasculopathy

Measurements of cerebrovascular reactivity (CVR) are essential for treatment decisions in moyamoya vasculopathy (MMV). Since MMV patients are often young or cognitively impaired, anesthesia is commonly used to limit motion artifacts. Our aim was to investigate the effect of anesthesia on the CVR in pediatric MMV. We compared the CVR with multidelay-ASL and BOLD MRI, using acetazolamide as a vascular stimulus, in all awake and anesthesia pediatric MMV scans at our institution. Since a heterogeneity in disease and treatment influences the CVR, we focused on the (unaffected) cerebellum. Ten awake and nine anesthetized patients were included. The post-acetazolamide CBF and ASL-CVR were significantly lower in anesthesia patients (47.1 ± 15.4 vs. 61.4 ± 12.1, p = 0.04; 12.3 ± 8.4 vs. 23.7 ± 12.2 mL/100 g/min, p = 0.03, respectively). The final BOLD-CVR increase (0.39 ± 0.58 vs. 3.6 ± 1.2% BOLD-change (mean/SD), p < 0.0001), maximum slope of increase (0.0050 ± 0.0040%/s vs. 0.017 ± 0.0059%, p < 0.0001), and time to maximum BOLD-increase (~463 ± 136 and ~697 ± 144 s, p = 0.0028) were all significantly lower in the anesthesia group. We conclude that the response to acetazolamide is distinctively different between awake and anesthetized MMV patients, and we hypothesize that these findings can also apply to other diseases and methods of measuring CVR under anesthesia. Considering that treatment decisions heavily depend on CVR status, caution is warranted when assessing CVR under anesthesia.

The many layers of BOLD. The effect of hypercapnic and hyperoxic stimuli on macro- and micro-vascular compartments quantified by CVR, M, and CBV across cortical depth

Ultra-high field functional magnetic resonance imaging (fMRI) offers the spatial resolution to measure neuronal activity at the scale of cortical layers. However, cortical depth dependent vascularization differences, such as a higher prevalence of macro-vascular compartments near the pial surface, have a confounding effect on depth-resolved blood-oxygen-level dependent (BOLD) fMRI signals. In the current study, we use hypercapnic and hyperoxic breathing conditions to quantify the influence of all venous vascular and micro-vascular compartments on laminar BOLD fMRI, as measured with gradient-echo (GE) and spin-echo (SE) scan sequences, respectively. We find that all venous vascular and micro-vascular compartments are capable of comparable theoretical maximum signal intensities, as represented by the M-value parameter. However, the capacity for vessel dilation, as reflected by the cerebrovascular reactivity (CVR), is approximately two and a half times larger for all venous vascular compartments combined compared to the micro-vasculature at superficial layers. Finally, there is roughly a 35% difference in estimates of CBV changes between all venous vascular and micro-vascular compartments, although this relative difference was approximately uniform across cortical depth. Thus, our results suggest that fMRI BOLD signal differences across cortical depth are likely caused by differences in dilation properties between macro- and micro-vascular compartments.

Cerebral oxygen extraction fraction MRI: Techniques and applications

The human brain constitutes 2% of the body's total mass but uses 20% of the oxygen. The rate of the brain's oxygen utilization can be derived from a knowledge of cerebral blood flow and the oxygen extraction fraction (OEF). Therefore, OEF is a key physiological parameter of the brain's function and metabolism. OEF has been suggested to be a useful biomarker in a number of brain diseases. With recent advances in MRI techniques, several MRI-based methods have been developed to measure OEF in the human brain. These MRI OEF techniques are based on the T₂ of blood, the blood signal phase, the magnetic susceptibility of blood-containing voxels, the effect of deoxyhemoglobin on signal behavior in extravascular tissue, and the calibration of the BOLD signal using gas inhalation. Compared to ¹⁵ O PET, which is considered the "gold standard" for OEF measurement, MRI-based techniques are non-invasive, radiation-free, and are more widely available. This article provides a review of these emerging MRI-based OEF techniques. We first briefly introduce the role of OEF in brain oxygen homeostasis. We then review the methodological aspects of different categories of MRI OEF techniques, including their signal mechanisms, acquisition methods, and data analyses. The strengths and limitations of the techniques are discussed. Finally, we review key applications of these techniques in physiological and pathological conditions.

Hemodynamic and metabolic changes during hypercapnia with normoxia and hyperoxia using pCASL and TRUST MRI in healthy adults

Blood oxygenation level-dependent (BOLD) or arterial spin labeling (ASL) MRI with hypercapnic stimuli allow for measuring cerebrovascular reactivity (CVR). Hypercapnic stimuli are also employed in calibrated BOLD functional MRI for quantifying neuronally-evoked changes in cerebral oxygen metabolism (CMRO₂). It is often assumed that hypercapnic stimuli (with or without hyperoxia) are iso-metabolic; increasing arterial CO₂ or O₂ does not affect CMRO₂. We evaluated the null hypothesis that two common hypercapnic stimuli, 'CO₂ in air' and carbogen, are iso-metabolic. TRUST and ASL MRI were used to measure the cerebral venous oxygenation and cerebral blood flow (CBF), from which the oxygen extraction fraction (OEF) and CMRO₂ were calculated for room-air, 'CO₂ in air' and carbogen. As expected, CBF significantly increased (9.9% ± 9.3% and 12.1% ± 8.8% for 'CO₂ in air' and carbogen, respectively). CMRO₂ decreased for 'CO₂ in air' (-13.4% ± 13.0%, p < 0.01) compared to room-air, while the CMRO₂ during carbogen did not significantly change. Our findings indicate that 'CO₂ in air' is not iso-metabolic, while carbogen appears to elicit a mixed effect; the CMRO₂ reduction during hypercapnia is mitigated when including hyperoxia. These findings can be important for interpreting measurements using hypercapnic or hypercapnic-hyperoxic (carbogen) stimuli.

Cardiorespiratory fitness is associated with increased middle cerebral arterial compliance and decreased cerebral blood flow in young healthy adults: A pulsed ASL MRI study

Cardiorespiratory fitness is thought to have beneficial effects on systemic vascular health, in part, by decreasing arterial stiffness. However, in the absence of non-invasive methods, it remains unknown whether this effect extends to the cerebrovasculature. The present study uses a novel pulsed arterial spin labelling (pASL) technique to explore the relationship between cardiorespiratory fitness and arterial compliance of the middle cerebral arteries (MCAC). Other markers of cerebrovascular health, including resting cerebral blood flow (CBF) and cerebrovascular reactivity to CO2 (CVRCO2) were also investigated. Eleven healthy males aged 21 ± 2 years with varying levels of cardiorespiratory fitness (maximal oxygen uptake (V · O2MAX) 38-76 ml/min/kg) underwent MRI scanning at 3 Tesla. Higher V · O2MAX was associated with greater MCAC (R2 = 0.64, p < 0.01) and lower resting grey matter CBF (R2 = 0.75, p < 0.01). However, V · O2MAX was not predictive of global grey matter BOLD-based CVR (R2 = 0.47, p = 0.17) or CBF-based CVR (R2 = 0.19, p = 0.21). The current experiment builds upon the established benefits of exercise on arterial compliance in the systemic vasculature, by showing that increased cardiorespiratory fitness is associated with greater cerebral arterial compliance in early adulthood.

Changes in volumetric and metabolic parameters relate to differences in exposure to sub-concussive head impacts

Structural and calibrated magnetic resonance imaging data were acquired on 44 collegiate football players prior to the season (PRE), following the first four weeks in-season (PTC) and one month after the last game (POST). Exposure data collected from g-Force accelerometers mounted to the helmet of each player were used to split participants into HIGH (N = 22) and LOW (N = 22) exposure groups, based on the frequency of impacts sustained by each athlete. Significant decreases in grey-matter volume specific to the HIGH group were documented at POST (P = 0.009), compared to baseline. Changes in resting cerebral blood flow (CBF0), corrected for partial volume effects, were observed within the HIGH group, throughout the season (P < 0.0001), suggesting that alterations in perfusion may follow exposure to sub-concussive collisions. Co-localized significant increases in cerebral metabolic rate of oxygen consumption (CMRO2|0) mid-season were also documented in the HIGH group, with respect to both PRE- and POST values. No physiological changes were observed in the LOW group. Therefore, cerebral metabolic demand may be elevated in players with greater exposure to head impacts. These results provide novel insight into the effects of sub-concussive collisions on brain structure and cerebrovascular physiology and emphasize the importance of multi-modal imaging for a complete characterization of cerebral health.

A frequency-domain machine learning method for dual-calibrated fMRI mapping of oxygen extraction fraction (OEF) and cerebral metabolic rate of oxygen consumption (CMRO2)

Magnetic resonance imaging (MRI) offers the possibility to non-invasively map the brain's metabolic oxygen consumption (CMRO₂), which is essential for understanding and monitoring neural function in both health and disease. However, in depth study of oxygen metabolism with MRI has so far been hindered by the lack of robust methods. One MRI method of mapping CMRO₂ is based on the simultaneous acquisition of cerebral blood flow (CBF) and blood oxygen level dependent (BOLD) weighted images during respiratory modulation of both oxygen and carbon dioxide. Although this dual-calibrated methodology has shown promise in the research setting, current analysis methods are unstable in the presence of noise and/or are computationally demanding. In this paper, we present a machine learning implementation for the multi-parametric assessment of dual-calibrated fMRI data. The proposed method aims to address the issues of stability, accuracy, and computational overhead, removing significant barriers to the investigation of oxygen metabolism with MRI. The method utilizes a time-frequency transformation of the acquired perfusion and BOLD-weighted data, from which appropriate feature vectors are selected for training of machine learning regressors. The implemented machine learning methods are chosen for their robustness to noise and their ability to map complex non-linear relationships (such as those that exist between BOLD signal weighting and blood oxygenation). An extremely randomized trees (ET) regressor is used to estimate resting blood flow and a multi-layer perceptron (MLP) is used to estimate CMRO₂ and the oxygen extraction fraction (OEF). Synthetic data with additive noise are used to train the regressors, with data simulated to cover a wide range of physiologically plausible parameters. The performance of the implemented analysis method is compared to published methods both in simulation and with in-vivo data (n = 30). The proposed method is demonstrated to significantly reduce computation time, error, and proportional bias in both CMRO₂ and OEF estimates. The introduction of the proposed analysis pipeline has the potential to not only increase the detectability of metabolic difference between groups of subjects, but may also allow for single subject examinations within a clinical context.

The association between BOLD-based cerebrovascular reactivity (CVR) and end-tidal CO2 in healthy subjects

Cerebrovascular reactivity (CVR) mapping using CO₂-inhalation can provide important insight into vascular health. At present, blood-oxygenation-level-dependent (BOLD) MRI acquisition is the most commonly used CVR method due to its high sensitivity, high spatial resolution, and relatively straightforward processing. However, large variations in CVR across subjects and across different sessions of the same subject are often observed, which can cloud the ability of this promising measure in detecting diseases or monitoring treatment responses. The present work aims to identify the physiological components underlying the observed variability in CVR data. When studying the association between CVR value and the subject’s CO₂ levels in a total of N = 253 healthy participants, we found that CVR was lower in individuals with a higher basal end-tidal CO₂, EtCO₂ (slope = −0.0036 ± 0.0008%/mmHg², p < 0.001), or with a greater EtCO₂ change (ΔEtCO₂) with hypercapnic condition (slope = −0.0072 ± 0.0018%/mmHg², p < 0.001). In a within-subject setting, when studying the CVR difference between two repeated scans (with repositioning) in relation to the corresponding differences in basal EtCO₂ and ΔEtCO₂ (n = 11), it was found that CVR values were lower if the basal EtCO₂ or ΔEtCO₂ during that particular scan session was greater. The present work suggests that basal physiological state and the level of hypercapnic stimulus intensity should be considered in application studies of CVR in order to reduce inter-subject and intra-subject variations in the data. Potential approaches to use these findings to reduce noise and augment sensitivity are proposed.

The potential for gas-free measurements of absolute oxygen metabolism during both baseline and activation states in the human brain

Quantitative functional magnetic resonance imaging methods make it possible to measure cerebral oxygen metabolism (CMRO₂) in the human brain. Current methods require the subject to breathe special gas mixtures (hypercapnia and hyperoxia). We tested a noninvasive suite of methods to measure absolute CMRO₂ in both baseline and dynamic activation states without the use of special gases: arterial spin labeling (ASL) to measure baseline and activation cerebral blood flow (CBF), with concurrent measurement of the blood oxygenation level dependent (BOLD) signal as a dynamic change in tissue R₂*; VSEAN to estimate baseline O₂ extraction fraction (OEF) from a measurement of venous blood R₂, which in combination with the baseline CBF measurement yields an estimate of baseline CMRO₂; and FLAIR-GESSE to measure tissue R₂' to estimate the scaling parameter needed for calculating the change in CMRO₂ in response to a stimulus with the calibrated BOLD method. Here we describe results for a study sample of 17 subjects (8 female, mean age = 25.3 years, range 21-31 years). The primary findings were that OEF values measured with the VSEAN method were in good agreement with previous PET findings, while estimates of the dynamic change in CMRO₂ in response to a visual stimulus were in good agreement between the traditional hypercapnia calibration and calibration based on R₂'. These results support the potential of gas-free methods for quantitative physiological measurements.

Neuroimaging of vascular reserve in patients with cerebrovascular diseases

Cerebrovascular reactivity, defined broadly as the ability of brain parenchyma to adjust cerebral blood flow in response to altered metabolic demand or a vasoactive stimulus, is being measured with increasing frequency and may have a use for portending new or recurrent stroke risk in patients with cerebrovascular disease. The purpose of this review is to outline (i) the physiological basis of variations in cerebrovascular reactivity, (ii) available approaches for measuring cerebrovascular reactivity in research and clinical settings, and (iii) clinically-relevant cerebrovascular reactivity findings in the context of patients with cerebrovascular disease, including atherosclerotic arterial steno-occlusion, non-atherosclerotic arterial steno-occlusion, anemia, and aging. Literature references summarizing safety considerations for these procedures and future directions for standardizing protocols and post-processing procedures across centers are presented in the specific context of major unmet needs in the setting of cerebrovascular disease.

Dual-calibrated fMRI measurement of absolute cerebral metabolic rate of oxygen consumption and effective oxygen diffusivity

Dual-calibrated fMRI is a multi-parametric technique that allows for the quantification of the resting oxygen extraction fraction (OEF), the absolute rate of cerebral metabolic oxygen consumption (CMRO2), cerebral vascular reactivity (CVR) and baseline perfusion (CBF). It combines measurements of arterial spin labelling (ASL) and blood oxygenation level dependent (BOLD) signal changes during hypercapnic and hyperoxic gas challenges. Here we propose an extension to this methodology that permits the simultaneous quantification of the effective oxygen diffusivity of the capillary network (DC). The effective oxygen diffusivity has the scope to be an informative biomarker and useful adjunct to CMRO2, potentially providing a non-invasive metric of microvascular health, which is known to be disturbed in a range of neurological diseases. We demonstrate the new method in a cohort of healthy volunteers (n = 19) both at rest and during visual stimulation. The effective oxygen diffusivity was found to be highly correlated with CMRO2 during rest and activation, consistent with previous PET observations of a strong correlation between metabolic oxygen demand and effective diffusivity. The increase in effective diffusivity during functional activation was found to be consistent with previously reported increases in capillary blood volume, supporting the notion that measured oxygen diffusivity is sensitive to microvascular physiology.

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.

Calibrated fMRI for mapping absolute CMRO2: Practicalities and prospects

Functional magnetic resonance imaging (fMRI) is an essential workhorse of modern neuroscience, providing valuable insight into the functional organisation of the brain. The physiological mechanisms underlying the blood oxygenation level dependent (BOLD) effect are complex and preclude a straightforward interpretation of the signal. However, by employing appropriate calibration of the BOLD signal, quantitative measurements can be made of important physiological parameters including the absolute rate of cerebral metabolic oxygen consumption or oxygen metabolism (CMRO₂) and oxygen extraction (OEF). The ability to map such fundamental parameters has the potential to greatly expand the utility of fMRI and to broaden its scope of application in clinical research and clinical practice. In this review article we discuss some of the practical issues related to the calibrated-fMRI approach to the measurement of CMRO₂. We give an overview of the necessary precautions to ensure high quality data acquisition, and explore some of the pitfalls and challenges that must be considered as it is applied and interpreted in a widening array of diseases and research questions.

Cerebrovascular reactivity (CVR) MRI with CO2 challenge: A technical review

Cerebrovascular reactivity (CVR) is an indicator of cerebrovascular reserve and provides important information about vascular health in a range of brain conditions and diseases. Unlike steady-state vascular parameters, such as cerebral blood flow (CBF) and cerebral blood volume (CBV), CVR measures the ability of cerebral vessels to dilate or constrict in response to challenges or maneuvers. Therefore, CVR mapping requires a physiological challenge while monitoring the corresponding hemodynamic changes in the brain. The present review primarily focuses on methods that use CO2 inhalation as a physiological challenge while monitoring changes in hemodynamic MRI signals. CO2 inhalation has been increasingly used in CVR mapping in recent literature due to its potency in causing vasodilation, rapid onset and cessation of the effect, as well as advances in MRI-compatible gas delivery apparatus. In this review, we first discuss the physiological basis of CVR mapping using CO2 inhalation. We then review the methodological aspects of CVR mapping, including gas delivery apparatus, the timing paradigm of the breathing challenge, the MRI imaging sequence, and data analysis. In addition, we review alternative approaches for CVR mapping that do not require CO2 inhalation.

Multiparametric measurement of cerebral physiology using calibrated fMRI

The ultimate goal of calibrated fMRI is the quantitative imaging of oxygen metabolism (CMRO₂), and this has been the focus of numerous methods and approaches. However, one underappreciated aspect of this quest is that in the drive to measure CMRO₂, many other physiological parameters of interest are often acquired along the way. This can significantly increase the value of the dataset, providing greater information that is clinically relevant, or detail that can disambiguate the cause of signal variations. This can also be somewhat of a double-edged sword: calibrated fMRI experiments combine multiple parameters into a physiological model that requires multiple steps, thereby providing more opportunity for error propagation and increasing the noise and error of the final derived values. As with all measurements, there is a trade-off between imaging time, spatial resolution, coverage, and accuracy. In this review, we provide a brief overview of the benefits and pitfalls of extracting multiparametric measurements of cerebral physiology through calibrated fMRI experiments.

Current understanding of chronic total occlusion of the internal carotid artery

At present, there is limited understanding of chronic total occlusion (CTO) of the internal carotid artery (ICA). Therefore, the present report collected related cases from PubMed and reviewed the literature. Cerebral vessels may form collateral circulation immediately or gradually following CTO of the ICA. The natural history of CTO of the ICA includes a variety of outcomes, all of which are biased toward a non-benign progressive process and are characterized by insufficient cerebral perfusion, embolus detachment and cognitive dysfunction. The majority of cases of CTO of the ICA require treatment. In early studies, the results of external-ICA bypass were unsatisfactory, while recanalization is now considered the only viable option. The current treatment indications mainly depend on the degree of injury to the cerebrovascular reserve and the extent to which the oxygen extraction fraction is increased. The length, height and duration of ICA occlusion are also relevant, though more frequently, the condition depends on multiple factors. Endovascular interventional recanalization, carotid endarterectomy (CEA) and hybrid surgery may be conducted in a select group of patients. As novel materials are developed, the success rate of simple recanalization may gradually increase; however, hybrid surgery may be more representative of the current trend, as advanced CEA can remove carotid atherosclerosis plaques, thus reducing the technological demands of the subsequent interventional recanalization. There are many complications that may result from recanalization following CTO of the ICA, including hyperperfusion and technical errors; therefore, the operation must be conducted carefully. If the recanalization is successful, it typically results in a stable improvement of patient condition in the long term. However, despite these conclusions, more studies are required in the future to further improve current understanding of CTO of the ICA.

Validation and optimization of hypercapnic-calibrated fMRI from oxygen-sensitive two-photon microscopy

Hypercapnic-calibrated fMRI allows the estimation of the relative changes in the cerebral metabolic rate of oxygen (rCMRO2) from combined BOLD and arterial spin labelling measurements during a functional task, and promises to permit more quantitative analyses of brain activity patterns. The estimation relies on a macroscopic model of the BOLD effect that balances oxygen delivery and consumption to predict haemoglobin oxygenation and the BOLD signal. The accuracy of calibrated fMRI approaches has not been firmly established, which is limiting their broader adoption. We use our recently developed microscopic vascular anatomical network model in mice as a ground truth simulator to test the accuracy of macroscopic, lumped-parameter BOLD models. In particular, we investigate the original Davis model and a more recent heuristic simplification. We find that these macroscopic models are inaccurate using the originally defined parameters, but that the accuracy can be significantly improved by redefining the model parameters to take on new values. In particular, we find that the parameter α that relates cerebral blood-volume changes to cerebral blood-flow changes is significantly smaller than typically assumed and that the optimal value changes with magnetic field strength. The results are encouraging in that they support the use of simple BOLD models to quantify BOLD signals, but further work is needed to understand the physiological interpretation of the redefined model parameters.This article is part of the themed issue 'Interpreting BOLD: a dialogue between cognitive and cellular neuroscience'.

Mapping the pharmacological modulation of brain oxygen metabolism: The effects of caffeine on absolute CMRO2 measured using dual calibrated fMRI

This study aims to map the acute effects of caffeine ingestion on grey matter oxygen metabolism and haemodynamics with a novel MRI method. Sixteen healthy caffeine consumers (8 males, age=24.7±5.1) were recruited to this randomised, double-blind, placebo-controlled study. Each participant was scanned on two days before and after the delivery of an oral caffeine (250mg) or placebo capsule. Our measurements were obtained with a newly proposed estimation approach applied to data from a dual calibration fMRI experiment that uses hypercapnia and hyperoxia to modulate brain blood flow and oxygenation. Estimates were based on a forward model that describes analytically the contributions of cerebral blood flow (CBF) and of the measured end-tidal partial pressures of CO2 and O2 to the acquired dual-echo GRE signal. The method allows the estimation of grey matter maps of: oxygen extraction fraction (OEF), CBF, CBF-related cerebrovascular reactivity (CVR) and cerebral metabolic rate of oxygen consumption (CMRO2). Other estimates from a multi inversion time ASL acquisition (mTI-ASL), salivary samples of the caffeine concentration and behavioural measurements are also reported. We observed significant differences between caffeine and placebo on average across grey matter, with OEF showing an increase of 15.6% (SEM±4.9%, p<0.05) with caffeine, while CBF and CMRO2 showed differences of -30.4% (SEM±1.6%, p<0.01) and -18.6% (SEM±2.9%, p<0.01) respectively with caffeine administration. The reduction in oxygen metabolism found is somehow unexpected, but consistent with a hypothesis of decreased energetic demand, supported by previous electrophysiological studies reporting reductions in spectral power with EEG. Moreover the maps of the physiological parameters estimated illustrate the spatial distribution of changes across grey matter enabling us to localise the effects of caffeine with voxel-wise resolution. CBF changes were widespread as reported by previous findings, while changes in OEF were found to be more restricted, leading to unprecedented mapping of significant CMRO2 reductions mainly in frontal gyrus, parietal and occipital lobes. In conclusion, we propose the estimation framework based on our novel forward model with a dual calibrated fMRI experiment as a viable MRI method to map the effects of drugs on brain oxygen metabolism and haemodynamics with voxel-wise resolution.

Post-stimulus fMRI and EEG responses: Evidence for a neuronal origin hypothesised to be inhibitory

Post-stimulus undershoots, negative responses following cessation of stimulation, are widely observed in functional magnetic resonance (fMRI) blood oxygenation level dependent (BOLD) data. However, the debate surrounding whether the origin of this response phase is neuronal or vascular, and whether it provides functionally relevant information, that is additional to what is contained in the primary response, means that undershoots are widely overlooked. We simultaneously recorded electroencephalography (EEG), BOLD and cerebral blood-flow (CBF) [obtained from arterial spin labelled (ASL) fMRI] fMRI responses to hemifield checkerboard stimulation to test the potential neural origin of the fMRI post-stimulus undershoot. The post-stimulus BOLD and CBF signal amplitudes in both contralateral and ipsilateral visual cortex depended on the post-stimulus power of the occipital 8-13Hz (alpha) EEG neuronal activity, such that trials with highest EEG power showed largest fMRI undershoots in contralateral visual cortex. This correlation in post-stimulus EEG-fMRI responses was not predicted by the primary response amplitude. In the contralateral visual cortex we observed a decrease in both cerebral rate of oxygen metabolism (CMRO₂) and CBF during the post-stimulus phase. In addition, the coupling ratio (n) between CMRO₂ and CBF was significantly lower during the positive contralateral primary response phase compared with the post-stimulus phase and we propose that this reflects an altered balance of excitatory and inhibitory neuronal activity. Together our data provide strong evidence that the post-stimulus phase of the BOLD response has a neural origin which reflects, at least partially, an uncoupling of the neuronal responses driving the primary and post-stimulus responses, explaining the uncoupling of the signals measured in the two response phases. We suggest our results are consistent with inhibitory processes driving the post-stimulus EEG and fMRI responses. We therefore propose that new methods are required to model the post-stimulus and primary responses independently, enabling separate investigation of response phases in cognitive function and neurological disease.

Graded Hypercapnia-Calibrated BOLD: Beyond the Iso-metabolic Hypercapnic Assumption

Calibrated BOLD is a promising technique that overcomes the sensitivity of conventional fMRI to the cerebrovascular state; measuring either the basal level, or the task-induced response of cerebral metabolic rate of oxygen consumption (CMRO₂). The calibrated BOLD method is susceptible to errors in the measurement of the calibration parameter M, the theoretical BOLD signal change that would occur if all deoxygenated hemoglobin were removed. The original and most popular method for measuring M uses hypercapnia (an increase in arterial CO₂), making the assumption that it does not affect CMRO₂. This assumption has since been challenged and recent studies have used a corrective term, based on literature values of a reduction in basal CMRO₂ with hypercapnia. This is not ideal, as this value may vary across subjects and regions of the brain, and will depend on the level of hypercapnia achieved. Here we propose a new approach, using a graded hypercapnia design and the assumption that CMRO₂ changes linearly with hypercapnia level, such that we can measure M without assuming prior knowledge of the scale of CMRO₂ change. Through use of a graded hypercapnia gas challenge, we are able to remove the bias caused by a reduction in basal CMRO₂ during hypercapnia, whilst simultaneously calculating the dose-wise CMRO₂ change with hypercapnia. When compared with assuming no change in CMRO₂, this approach resulted in significantly lower M-values in both visual and motor cortices, arising from significant dose-dependent hypercapnia reductions in basal CMRO₂ of 1.5 ± 0.6%/mmHg (visual) and 1.8 ± 0.7%/mmHg (motor), where mmHg is the unit change in end-tidal CO₂ level. Variability in the basal CMRO₂ response to hypercapnia, due to experimental differences and inter-subject variability, is accounted for in this approach, unlike previous correction approaches, which use literature values. By incorporating measurement of, and correction for, the reduction in basal CMRO₂ during hypercapnia in the measurement of M-values, application of our approach will correct for an overestimation in both CMRO₂ task-response values and absolute CMRO₂.

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.

Regional Reproducibility of BOLD Calibration Parameter M, OEF and Resting-State CMRO2 Measurements with QUO2 MRI

The current generation of calibrated MRI methods goes beyond simple localization of task-related responses to allow the mapping of resting-state cerebral metabolic rate of oxygen (CMRO₂) in micromolar units and estimation of oxygen extraction fraction (OEF). Prior to the adoption of such techniques in neuroscience research applications, knowledge about the precision and accuracy of absolute estimates of CMRO₂ and OEF is crucial and remains unexplored to this day. In this study, we addressed the question of methodological precision by assessing the regional inter-subject variance and intra-subject reproducibility of the BOLD calibration parameter M, OEF, O₂ delivery and absolute CMRO₂ estimates derived from a state-of-the-art calibrated BOLD technique, the QUantitative O2 (QUO2) approach. We acquired simultaneous measurements of CBF and R2* at rest and during periods of hypercapnia (HC) and hyperoxia (HO) on two separate scan sessions within 24 hours using a clinical 3 T MRI scanner. Maps of M, OEF, oxygen delivery and CMRO₂, were estimated from the measured end-tidal O₂, CBF_0, CBF_HC/HO and R2*_HC/HO. Variability was assessed by computing the between-subject coefficients of variation (bwCV) and within-subject CV (wsCV) in seven ROIs. All tests GM-averaged values of CBF₀, M, OEF, O₂ delivery and CMRO₂ were: 49.5 ± 6.4 mL/100 g/min, 4.69 ± 0.91%, 0.37 ± 0.06, 377 ± 51 μmol/100 g/min and 143 ± 34 μmol/100 g/min respectively. The variability of parameter estimates was found to be the lowest when averaged throughout all GM, with general trends toward higher CVs when averaged over smaller regions. Among the MRI measurements, the most reproducible across scans was R2*₀ (wsCV_GM = 0.33%) along with CBF₀ (wsCV_GM = 3.88%) and R2*_HC (wsCV_GM = 6.7%). CBF_HC and R2*_HO were found to have a higher intra-subject variability (wsCV_GM = 22.4% and wsCV_GM = 16% respectively), which is likely due to propagation of random measurement errors, especially for CBF_HC due to the low contrast-to-noise ratio intrinsic to ASL. Reproducibility of the QUO2 derived estimates were computed, yielding a GM intra-subject reproducibility of 3.87% for O₂ delivery, 16.8% for the M value, 13.6% for OEF and 15.2% for CMRO₂. Although these results focus on the precision of the QUO2 method, rather than the accuracy, the information will be useful for calculation of statistical power in future validation studies and ultimately for research applications of the method. The higher test-retest variability for the more extensively modeled parameters (M, OEF, and CMRO₂) highlights the need for further improvement of acquisition methods to reduce noise levels.