Improved fMRI calibration: precisely controlled hyperoxic versus hypercapnic stimuli

Contrast-free dynamic susceptibility contrast using sinusoidal and bolus oxygenation challenges

Deoxygenation-based dynamic susceptibility contrast (dDSC) MRI uses respiratory challenges as a source of endogenous contrast as an alternative to gadolinium injection. These gas challenges induce T2*-weighted MRI signal losses, after which tracer kinetics modeling was applied to calculate cerebral perfusion. This work compares three gas challenges, desaturation (transient hypoxia), resaturation (transient normoxia), and SineO2 (sinusoidal modulation of end-tidal oxygen pressures) in a cohort of 10 healthy volunteers (age 37 ± 11 years; 60% female). Perfusion estimates consisted of cerebral blood flow (CBF), cerebral blood volume (CBV), and mean transit time (MTT). Calculations were computed using a traditional tracer kinetics model in the time domain for desaturation and resaturation and in the frequency domain for SineO2. High correlations and limits of agreement were observed among the three deoxygenation-based paradigms for CBV, although MTT and CBF estimates varied with the hypoxic stimulus. Cross-modality correlation with gadolinium DSC was lower, particularly for MTT, but on a par with agreement between the other perfusion references. Overall, this work demonstrated the feasibility and reliability of oxygen respiratory challenges to measure brain perfusion. Additional work is needed to assess the utility of dDSC in the diagnostic evaluation of various pathologies such as ischemic strokes, brain tumors, and neurodegenerative diseases.

Multi-Echo Investigations of Positive and Negative CBF and Concomitant BOLD Changes: Positive and negative CBF and BOLD changes

Unlike the positive blood oxygenation level-dependent (BOLD) response (PBR), commonly taken as an indication of an 'activated' brain region, the physiological origin of negative BOLD signal changes (i.e. a negative BOLD response, NBR), also referred to as 'deactivation' is still being debated. In this work, an attempt was made to gain a better understanding of the underlying mechanism by obtaining a comprehensive measure of the contributing cerebral blood flow (CBF) and its relationship to the NBR in the human visual cortex, in comparison to a simultaneously induced PBR in surrounding visual regions. To overcome the low signal-to-noise ratio (SNR) of CBF measurements, a newly developed multi-echo version of a center-out echo planar-imaging (EPI) readout was employed with pseudo-continuous arterial spin labeling (pCASL). It achieved very short echo and inter-echo times and facilitated a simultaneous detection of functional CBF and BOLD changes at 3 T with improved sensitivity. Evaluations of the absolute and relative changes of CBF and the effective transverse relaxation rate,R₂^* the coupling ratios, and their dependence on CBF at rest, CBF_rest indicated differences between activated and deactivated regions. Analysis of the shape of the respective functional responses also revealed faster negative responses with more pronounced post-stimulus transients. Resulting differences in the flow-metabolism coupling ratios were further examined for potential distinctions in the underlying neuronal contributions.

The Brain in Motion II Study: study protocol for a randomized controlled trial of an aerobic exercise intervention for older adults at increased risk of dementia

BACKGROUND

There remains no effective intervention capable of reversing most cases of dementia. Current research is focused on prevention by addressing risk factors that are shared between cardiovascular disease and dementia (e.g., hypertension) before the cognitive, functional, and behavioural symptoms of dementia manifest. A promising preventive treatment is exercise. This study describes the methods of a randomized controlled trial (RCT) that assesses the effects of aerobic exercise and behavioural support interventions in older adults at increased risk of dementia due to genetic and/or cardiovascular risk factors. The specific aims are to determine the effect of aerobic exercise on cognitive performance, explore the biological mechanisms that influence cognitive performance after exercise training, and determine if changes in cerebrovascular physiology and function persist 1 year after a 6-month aerobic exercise intervention followed by a 1-year behavioural support programme (at 18 months).

METHODS

We will recruit 264 participants (aged 50-80 years) at elevated risk of dementia. Participants will be randomly allocated into one of four treatment arms: (1) aerobic exercise and health behaviour support, (2) aerobic exercise and no health behaviour support, (3) stretching-toning and health behaviour support, and (4) stretching-toning and no health behaviour support. The aerobic exercise intervention will consist of three supervised walking/jogging sessions per week for 6 months, whereas the stretching-toning control intervention will consist of three supervised stretching-toning sessions per week also for 6 months. Following the exercise interventions, participants will receive either 1 year of ongoing telephone behavioural support or no telephone support. The primary aim is to determine the independent effect of aerobic exercise on a cognitive composite score in participants allocated to this intervention compared to participants allocated to the stretching-toning group. The secondary aims are to examine the effects of aerobic exercise on a number of secondary outcomes and determine whether aerobic exercise-related changes persist after a 1-year behavioural support programme (at 18 months).

DISCUSSION

This study will address knowledge gaps regarding the underlying mechanisms of the pro-cognitive effects of exercise by examining the potential mediating factors, including cerebrovascular/physiological, neuroimaging, sleep, and genetic factors that will provide novel biologic evidence on how aerobic exercise can prevent declines in cognition with ageing.

TRIAL REGISTRATION

ClinicalTrials.gov NCT03035851 . Registered on 30 January 2017.

Cerebrovascular Reactivity Mapping Without Gas Challenges: A Methodological Guide

Cerebrovascular reactivity (CVR) is defined as the ability of vessels to alter their caliber in response to vasoactive factors, by means of dilating or constricting, in order to increase or decrease regional cerebral blood flow (CBF). Importantly, CVR may provide a sensitive biomarker for pathologies where vasculature is compromised. Furthermore, the spatiotemporal dynamics of CVR observed in healthy subjects, reflecting regional differences in cerebral vascular tone and response, may also be important in functional MRI studies based on neurovascular coupling mechanisms. Assessment of CVR is usually based on the use of a vasoactive stimulus combined with a CBF measurement technique. Although transcranial Doppler ultrasound has been frequently used to obtain global flow velocity measurements, MRI techniques are being increasingly employed for obtaining CBF maps. For the vasoactive stimulus, vasodilatory hypercapnia is usually induced through the manipulation of respiratory gases, including the inhalation of increased concentrations of carbon dioxide. However, most of these methods require an additional apparatus and complex setups, which not only may not be well-tolerated by some populations but are also not widely available. For these reasons, strategies based on voluntary breathing fluctuations without the need for external gas challenges have been proposed. These include the task-based methodologies of breath holding and paced deep breathing, as well as a new generation of methods based on spontaneous breathing fluctuations during resting-state. Despite the multitude of alternatives to gas challenges, existing literature lacks definitive conclusions regarding the best practices for the vasoactive modulation and associated analysis protocols. In this work, we perform an extensive review of CVR mapping techniques based on MRI and CO2 variations without gas challenges, focusing on the methodological aspects of the breathing protocols and corresponding data analysis. Finally, we outline a set of practical guidelines based on generally accepted practices and available data, extending previous reports and encouraging the wider application of CVR mapping methodologies in both clinical and academic MRI settings.

Quantitative perfusion mapping with induced transient hypoxia using BOLD MRI

PURPOSE

Gadolinium-based dynamic susceptibility contrast (DSC) is commonly used to characterize blood flow in patients with stroke and brain tumors. Unfortunately, gadolinium contrast administration has been associated with adverse reactions and long-term accumulation in tissues. In this work, we propose an alternative deoxygenation-based DSC (dDSC) method that uses a transient hypoxia gas paradigm to deliver a bolus of paramagnetic deoxygenated hemoglobin to the cerebral vasculature for perfusion imaging.

METHODS

Through traditional DSC tracer kinetic modeling, the MR signal change induced by this hypoxic bolus can be used to generate regional perfusion maps of cerebral blood flow, cerebral blood volume, and mean transit time. This gas paradigm and blood-oxygen-level-dependent (BOLD)-MRI were performed concurrently on a cohort of 66 healthy and chronically anemic subjects (age 23.5 ± 9.7, female 64%).

RESULTS

Our results showed reasonable global and regional agreement between dDSC and other flow techniques, such as phase contrast and arterial spin labeling.

CONCLUSION

In this proof-of-concept study, we demonstrated the feasibility of using transient hypoxia to generate a contrast bolus that mimics the effect of gadolinium and yields reasonable perfusion estimates. Looking forward, optimization of the hypoxia boluses and measurement of the arterial-input function is necessary to improve the accuracy of dDSC. Additionally, a cross-validation study of dDSC and DSC in brain tumor and ischemic stroke subjects is warranted to evaluate the clinical diagnostic utility of this approach.

Cerebral Metabolic Changes During Visuomotor Adaptation Assessed Using Quantitative fMRI

The brain retains a lifelong ability to adapt through learning and in response to injury or disease-related damage, a process known as functional neuroplasticity. The neural energetics underlying functional brain plasticity have not been thoroughly investigated experimentally in the healthy human brain. A better understanding of the blood flow and metabolic changes that accompany motor skill acquisition, and which facilitate plasticity, is needed before subsequent translation to treatment interventions for recovery of function in disease. The aim of the current study was to characterize cerebral blood flow (CBF) and oxygen consumption (relative CMRO2) responses, using calibrated fMRI conducted in 20 healthy participants, during performance of a serial reaction time task which induces rapid motor adaptation. Regions of interest (ROIs) were defined from areas showing task-induced BOLD and CBF responses that decreased over time. BOLD, CBF and relative CMRO2 responses were calculated for each block of the task. Motor and somatosensory cortices and the cerebellum showed statistically significant positive responses to the task compared to baseline, but with decreasing amplitudes of BOLD, CBF, and CMRO2 response as the task progressed. In the cerebellum, there was a sustained positive BOLD response in the absence of a significant CMRO2 increase from baseline, for all but the first task blocks. This suggests that the brain may continue to elevate the supply energy even after CMRO2 has returned to near baseline levels. Relying on BOLD fMRI data alone in studies of plasticity may not reveal the nature of underlying metabolic responses and their changes over time. Calibrated fMRI approaches may offer a more complete picture of the energetic changes supporting plasticity and learning.

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.

Improving the Assessment of Breath-Holding Induced Cerebral Vascular Reactivity Using a Multiband Multi-echo ASL/BOLD Sequence

Breath holding (BH) is a viable vasodilatory stimulus for calculating functional MRI-derived cerebral vascular reactivity (CVR). The BH technique suffers from reduced repeatability compared with gas inhalation techniques; however, extra equipment is needed to perform gas inhalation techniques, and this equipment is not available at all institutions. This study aimed to determine the sensitivity and repeatability of BH activation and CVR using a multiband multi-echo simultaneous arterial spin labelling/blood oxygenation level dependent (ASL/BOLD) sequence. Whole-brain images were acquired in 14 volunteers. Ten subjects returned for repeat imaging. Each subject performed four cycles of 16 s BH on expiration interleaved with paced breathing. Following standard preprocessing, the echoes were combined using a T2*-weighted approach. BOLD and ASL BH activation was computed, and CVR was then determined as the percent signal change related to the activation. The "M" parameter from the Davis Model was also computed by incorporating the ASL signal. Our results showed higher BH activation strength, volume, and repeatability for the combined multi-echo (MEC) data compared with the single-echo data. MEC CVR also had higher repeatability, sensitivity, specificity, and reliability compared with the single-echo BOLD data. These data support the usefulness of an MBME ASL/BOLD acquisition for BH CVR and M measurements.

Cerebral blood volume changes during the BOLD post-stimulus undershoot measured with a combined normoxia/hyperoxia method

Cerebral blood flow (CBF) and blood oxygenation level dependent (BOLD) signal measurements make it possible to estimate steady-state changes in the cerebral metabolic rate of oxygen (CMRO2) with a calibrated BOLD method. However, extending this approach to measure the dynamics of CMRO2 requires an additional assumption: that deoxygenated cerebral blood volume (CBVdHb) follows CBF in a predictable way. A test-case for this assumption is the BOLD post-stimulus undershoot, for which one proposed explanation is a strong uncoupling of flow and blood volume with an elevated level of CBVdHb during the post-stimulus period compared to baseline due to slow blood volume recovery (Balloon Model). A challenge in testing this model is that CBVdHb differs from total blood volume, which can be measured with other techniques. In this study, the basic hypothesis of elevated CBVdHb during the undershoot was tested, based on the idea that the BOLD signal change when a subject switches from breathing a normoxic gas to breathing a hyperoxic gas is proportional to the absolute CBVdHb. In 19 subjects (8F), dual-echo BOLD responses were measured in primary visual cortex during a flickering radial checkerboard stimulus in normoxia, and the identical experiment was repeated in hyperoxia (50% O2/balance N2). The BOLD signal differences between normoxia and hyperoxia for the pre-stimulus baseline, stimulus, and post-stimulus periods were compared using an equivalent BOLD signal calculated from measured R2* changes to eliminate signal drifts. Relative to the pre-stimulus baseline, the average BOLD signal change from normoxia to hyperoxia was negative during the undershoot period (p = 0.0251), consistent with a reduction of CBVdHb and contrary to the prediction of the Balloon Model. Based on these results, the BOLD post-stimulus undershoot does not represent a case of strong uncoupling of CBVdHb and CBF, supporting the extension of current calibrated BOLD methods to estimate the dynamics of CMRO2.

Modeling of dynamic cerebrovascular reactivity to spontaneous and externally induced CO2 fluctuations in the human brain using BOLD-fMRI

In this work, we investigate the regional characteristics of the dynamic interactions between arterial CO₂ and BOLD (dynamic cerebrovascular reactivity - dCVR) during normal breathing and hypercapnic, externally induced step CO₂ challenges. To obtain dCVR curves at each voxel, we use a custom set of basis functions based on the Laguerre and gamma basis sets. This allows us to obtain robust dCVR estimates both in larger regions of interest (ROIs), as well as in individual voxels. We also implement classification schemes to identify brain regions with similar dCVR characteristics. Our results reveal considerable variability of dCVR across different brain regions, as well as during different experimental conditions (normal breathing and hypercapnic challenges), suggesting a differential response of cerebral vasculature to spontaneous CO₂ fluctuations and larger, externally induced CO₂ changes that are possibly associated with the underlying differences in mean arterial CO₂ levels. The clustering results suggest that anatomically distinct brain regions are characterized by different dCVR curves that in some cases do not exhibit the standard, positive valued curves that have been previously reported. They also reveal a consistent set of dCVR cluster shapes for resting and forcing conditions, which exhibit different distribution patterns across brain voxels.

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₂.

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.

MRI-based methods for quantification of the cerebral metabolic rate of oxygen

The brain depends almost entirely on oxidative metabolism to meet its significant energy requirements. As such, the cerebral metabolic rate of oxygen (CMRO2) represents a key measure of brain function. Quantification of CMRO2 has helped elucidate brain functional physiology and holds potential as a clinical tool for evaluating neurological disorders including stroke, brain tumors, Alzheimer's disease, and obstructive sleep apnea. In recent years, a variety of magnetic resonance imaging (MRI)-based CMRO2 quantification methods have emerged. Unlike positron emission tomography - the current "gold standard" for measurement and mapping of CMRO2 - MRI is non-invasive, relatively inexpensive, and ubiquitously available in modern medical centers. All MRI-based CMRO2 methods are based on modeling the effect of paramagnetic deoxyhemoglobin on the magnetic resonance signal. The various methods can be classified in terms of the MRI contrast mechanism used to quantify CMRO2: T2*, T2', T2, or magnetic susceptibility. This review article provides an overview of MRI-based CMRO2 quantification techniques. After a brief historical discussion motivating the need for improved CMRO2 methodology, current state-of-the-art MRI-based methods are critically appraised in terms of their respective tradeoffs between spatial resolution, temporal resolution, and robustness, all of critical importance given the spatially heterogeneous and temporally dynamic nature of brain energy requirements.

Reproducibility of Brachial Vascular Changes with Alterations in End-Tidal Carbon Dioxide

The purpose of this study was to examine the reproducibility of the peripheral vascular response to hypercapnia. Healthy college-aged men (n = 7) and women (n = 10) underwent an iso-oxic 10-mm Hg increase in PetCO2 for 12 min. Brachial artery diameter changes were measured using ultrasound imaging. Two tests were completed on day 1 with 15 min of rest between tests. Tests were repeated on day 2. Paired t-tests, Bland-Altman plots and intra-class correlations (ICCs) determined reproducibility. There were no significant differences in peak dilation within day (5.33 ± 3.73% vs. 4.52 ± 2.49%, p = 0.378). The within-day ICC was poor (0.213). Within-day time-to-peak dilation did not significantly differ (660.0 ± 231.8 s vs. 602.7 ± 259.9 s, p = 0.379), and the ICC was fair (0.416, p = 0.113). Between-day peak dilation did not significantly differ (5.24 ± 3.84% vs. 4.71 ± 3.17%, p = 0.123), and the ICC was fair (0.419). Hypercapnia-induced brachial artery dilation is similar within day and between days. The ICC for peak dilation suggests the methodology is not reproducible.

New horizons in neurometabolic and neurovascular coupling from calibrated fMRI

Neurovascular coupling relates changes in neuronal activity to constriction/dilation of microvessels. However neurometabolic coupling, which is less well known, relates alterations in neuronal activity with metabolic demands. The link between the blood oxygenation level dependent (BOLD) signal and neural activity opened doors for functional MRI (fMRI) to be a powerful neuroimaging tool in the neurosciences. But due to the complex makeup of BOLD contrast, researchers began to investigate the relationship between BOLD signal and blood flow and/or volume changes during functional brain activation, which together provided the tools to measure oxygen consumption on the basis of the biophysical model of BOLD. This field is called calibrated fMRI, thereby allowed probing of both neurometabolic and neurovascular couplings for a variety of health conditions in animals and humans. Calibrated fMRI may provide brain disorder biomarkers that could be used for monitoring effective therapies.

Comparison of brachial dilatory responses to hypercapnia and reactive hyperemia

Flow-mediated dilation (FMD) relies on reactive hyperemia to stimulate the endothelium to release nitric oxide, causing smooth muscle relaxation. Hypercapnia also produces vasodilation, which is thought to be nitric oxide-independent. The purpose of this study was to compare and contrast the effects of hypercapnia and reactive hyperemia as stimuli for brachial artery dilation. On separate days, twenty-five participants underwent vasodilation studies via reactive hyperemia or hypercapnia (i.e. 10 mmHg increase in end-tidal carbon dioxide [PetCO2)]). During both studies changes in brachial artery diameter were recorded using continuous ultrasound imaging. Heart rate (HR) was measured throughout both tests. Resting HR (63  ±  11 versus 68  ±  14 beats min(-1), p  =  0.0027) and baseline brachial artery diameter measurements (4.57  ±  1.51 versus 5.28  ±  1.86 mm, p  =  0.022) were significantly different between reactive hyperemia and hypercapnia, respectively. HR at peak dilation (65  ±  11 versus 76  ±  14 beats min(-1), p  <  0.0001), peak vessel dilation (8.68  ±  4.50 versus 5.28  ±  1.86%, p  =  0.002), and time to peak dilation (90.8  ±  120.1 versus 658.3  ±  226.6 s, p  <  0.0001) were also significantly different between reactive hyperemia and hypercapnia. The dynamics by which reactive hyperemia and hypercapnia stimulate vasodilation appear to differ. Hypercapnia produces a smaller and slower vasodilatory effect than reactive hyperemia. Further research is necessary to better understand the mechanisms of vasodilation under hypercapnic conditions.

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.

Quantitative β mapping for calibrated fMRI

The metabolic and hemodynamic dependencies of the blood oxygenation level-dependent (BOLD) signal form the basis for calibrated fMRI, where the focus is on oxidative energy demanded by neural activity. An important part of calibrated fMRI is the power-law relationship between the BOLD signal and the deoxyhemoglobin concentration, which in turn is related to the ratio between oxidative demand (CMRO2) and blood flow (CBF). The power-law dependence between BOLD signal and deoxyhemoglobin concentration is signified by a scaling exponent β. Until recently most studies assumed a β value of 1.5, which is based on numerical simulations of the extravascular BOLD component. Since the basal value of CMRO2 and CBF can vary from subject-to-subject and/or region-to-region, a method to independently measure β in vivo should improve the accuracy of calibrated fMRI results. We describe a new method for β mapping through characterizing R2' - the most sensitive relaxation component of BOLD signal (i.e., the reversible magnetic susceptibility component that is predominantly of extravascular origin at high magnetic field) - as a function of intravascular magnetic susceptibility induced by an FDA-approved superparamagnetic contrast agent. In α-chloralose anesthetized rat brain, at 9.4 T, we measured β values of ~0.8 uniformly across large neocortical swathes, with lower magnitude and more heterogeneity in subcortical areas. Comparison of β maps in rats anesthetized with medetomidine and α-chloralose revealed that β is independent of neural activity levels at these resting states. We anticipate that this method for β mapping can help facilitate calibrated fMRI for clinical studies.

Brain region and activity-dependent properties of M for calibrated fMRI

Calibrated fMRI extracts changes in oxidative energy demanded by neural activity based on hemodynamic and metabolic dependencies of the blood oxygenation level-dependent (BOLD) response. This procedure requires the parameter M, which is determined from the dynamic range of the BOLD signal between deoxyhemoglobin (paramagnetic) and oxyhemoglobin (diamagnetic). Since it is unclear if the range of M-values in human calibrated fMRI is due to regional/state differences, we conducted a 9.4T study to measure M-values across brain regions in deep (α-chloralose) and light (medetomidine) anesthetized rats, as verified by electrophysiology. Because BOLD signal is captured differentially by gradient-echo (R2*) and spin-echo (R2) relaxation rates, we measured M-values by the product of the fMRI echo time and R2' (i.e., the reversible magnetic susceptibility component), which is given by the absolute difference between R2* and R2. While R2' mapping was shown to be dependent on the k-space sampling method used, at nominal spatial resolutions achieved at high magnetic field of 9.4T the M-values were quite homogenous across cortical gray matter. However cortical M-values varied in relation to neural activity between brain states. The findings from this study could improve precision of future calibrated fMRI studies by focusing on the global uniformity of M-values in gray matter across different resting activity levels.

The influence of mild carbon dioxide on brain functional homotopy using resting-state fMRI

Homotopy reflects the intrinsic functional architecture of the brain through synchronized spontaneous activity between corresponding bilateral regions, measured as voxel mirrored homotopic connectivity (VMHC). Hypercapnia is known to have clear impact on brain hemodynamics through vasodilation, but have unclear effect on neuronal activity. This study investigates the effect of hypercapnia on brain homotopy, achieved by breathing 5% carbon dioxide (CO2 ) gas mixture. A total of 14 healthy volunteers completed three resting state functional MRI (RS-fMRI) scans, the first and third under normocapnia and the second under hypercapnia. VMHC measures were calculated as the correlation between the BOLD signal of each voxel and its counterpart in the opposite hemisphere. Group analysis was performed between the hypercapnic and normocapnic VMHC maps. VMHC showed a diffused decrease in response to hypercapnia. Significant regional decreases in VMHC were observed in all anatomical lobes, except for the occipital lobe, in the following functional hierarchical subdivisions: the primary sensory-motor, unimodal, heteromodal, paralimbic, as well as in the following functional networks: ventral attention, somatomotor, default frontoparietal, and dorsal attention. Our observation that brain homotopy in RS-fMRI is affected by arterial CO2 levels suggests that caution should be used when comparing RS-fMRI data between healthy controls and patients with pulmonary diseases and unusual respiratory patterns such as sleep apnea or chronic obstructive pulmonary disease.

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.

Investigating the field-dependence of the Davis model: Calibrated fMRI at 1.5, 3 and 7T

Gas calibrated fMRI in its most common form uses hypercapnia in conjunction with the Davis model to quantify relative changes in the cerebral rate of oxygen consumption (CMRO2) in response to a functional stimulus. It is most commonly carried out at 3T but, as 7T research scanners are becoming more widespread and the majority of clinical scanners are still 1.5T systems, it is important to investigate whether the model used remains accurate across this range of field strengths. Ten subjects were scanned at 1.5, 3 and 7T whilst performing a bilateral finger-tapping task as part of a calibrated fMRI protocol, and the results were compared to a detailed signal model. Simulations predicted an increase in value and variation in the calibration parameter M with field strength. Two methods of defining experimental regions of interest (ROIs) were investigated, based on (a) BOLD signal and (b) BOLD responses within grey matter only. M values from the latter ROI were in closer agreement with theoretical predictions; however, reassuringly, ROI choice had less impact on CMRO2 than on M estimates. Relative changes in CMRO2 during motor tasks at 3 and 7T were in good agreement but were over-estimated at 1.5T as a result of the lower signal to noise ratio. This result is encouraging for future studies at 7T, but also highlights the impact of imaging and analysis choices (such as ASL sequence and ROI definition) on the calibration parameter M and on the calculation of CMRO2.

The effect of dissolved oxygen on the susceptibility of blood

PURPOSE

It has been predicted that, during hyperoxia, excess O2 dissolved in arterial blood will significantly alter the blood's magnetic susceptibility. This would confound the interpretation of the hyperoxia-induced blood oxygenation level-dependent signal as arising solely from changes in deoxyhemoglobin. This study, therefore, aimed to determine how dissolved O2 affects the susceptibility of blood.

THEORY AND METHODS

We present a comprehensive model for the effect of dissolved O2 on the susceptibility of blood and compare it with another recently published model, referred to here as the ideal gas model (IGM). For validation, distilled water and samples of bovine plasma were oxygenated over a range of hyperoxic O2 concentrations and their susceptibilities were determined using multiecho gradient echo phase imaging.

RESULTS

In distilled water and plasma, the measured changes in susceptibility were very linear, with identical slopes of 0.062 ppb/mm Hg of O2. This change was dramatically less than previously predicted using the IGM and was close to that predicted by our model. The primary source of error in the IGM is the overestimation of the volume fraction occupied by dissolved O2.

CONCLUSION

Under most physiological conditions, the susceptibility of dissolved O2 can be disregarded in MRI studies employing hyperoxia.

Assessing cerebrovascular reactivity abnormality by comparison to a reference atlas

Attribution of vascular pathophysiology to reductions in cerebrovascular reactivity (CVR) is confounded by subjective assessment and the normal variation between anatomic regions. This study aimed to develop an objective scoring assessment of abnormality. CVR was measured as the ratio of the blood-oxygen-level-dependent magnetic resonance signal response divided by an increase in CO2, standardized to eliminate variability. A reference normal atlas was generated by coregistering the CVR maps from 46 healthy subjects into a standard space and calculating the mean and standard deviation (s.d.) of CVR for each voxel. Example CVR studies from 10 patients with cerebral vasculopathy were assessed for abnormality, by normalizing each patient's CVR to the same standard space as the atlas, and assigning a z-score to each voxel relative to the mean and s.d. of the corresponding atlas voxel. Z-scores were color coded and superimposed on their anatomic scans to form CVR z-maps. We found the CVR z-maps provided an objective evaluation of abnormality, enhancing our appreciation of the extent and distribution of pathophysiology compared with CVR maps alone. We concluded that CVR z-maps provide an objective, improved form of evaluation for comparisons of voxel-specific CVR between subjects, and across tests sites.

Global intravascular and local hyperoxia contrast phase-based blood oxygenation measurements

The measurement of venous cerebral blood oxygenation (Yv) has potential applications in the study of patient groups where oxygen extraction and/or metabolism are compromised. It is also useful for fMRI studies to assess the stimulus-induced changes in Yv, particularly since basal Yv partially accounts for inter-subject variation in the haemodynamic response to a stimulus. A range of MRI-based methods of measuring Yv have been developed recently. Here, we use a method based on the change in phase in the MR image arising from the field perturbation caused by deoxygenated haemoglobin in veins. We build on the existing phase based approach (Method I), where Yv is measured in a large vein (such as the superior sagittal sinus) based on the field shift inside the vein with assumptions as to the vein's shape and orientation. We demonstrate two novel modifications which address limitations of this method. The first modification (Method II), maps the actual form of the vein, rather than assume a given shape and orientation. The second modification (Method III) uses the intra and perivascular phase change in response to a known change in Yv on hyperoxia to measure normoxic Yv in smaller veins. Method III can be applied to veins whose shape, size and orientation are not accurately known, thus allowing more localised measures of venous oxygenation. Results demonstrate that the use of an overly fine spatial filter caused an overestimation in Yv for Method I, whilst the measurement of Yv using Method II was less sensitive to this bias, giving Yv = 0.62 ± 0.03. Method III was applied to mapping of Yv in local veins across the brain, yielding a distribution of values with a mode of Yv = 0.661 ± 0.008.

Evidence that the negative BOLD response is neuronal in origin: a simultaneous EEG-BOLD-CBF study in humans

Unambiguous interpretation of changes in the BOLD signal is challenging because of the complex neurovascular coupling that translates changes in neuronal activity into the subsequent haemodynamic response. In particular, the neurophysiological origin of the negative BOLD response (NBR) remains incompletely understood. Here, we simultaneously recorded BOLD, EEG and cerebral blood flow (CBF) responses to 10 s blocks of unilateral median nerve stimulation (MNS) in order to interrogate the NBR. Both negative BOLD and negative CBF responses to MNS were observed in the same region of the ipsilateral primary sensorimotor cortex (S1/M1) and calculations showed that MNS induced a decrease in the cerebral metabolic rate of oxygen consumption (CMRO2) in this NBR region. The ∆CMRO2/∆CBF coupling ratio (n) was found to be significantly larger in this ipsilateral S1/M1 region (n=0.91±0.04, M=10.45%) than in the contralateral S1/M1 (n=0.65±0.03, M=10.45%) region that exhibited a positive BOLD response (PBR) and positive CBF response, and a consequent increase in CMRO2 during MNS. The fMRI response amplitude in ipsilateral S1/M1 was negatively correlated with both the power of the 8-13 Hz EEG mu oscillation and somatosensory evoked potential amplitude. Blocks in which the largest magnitude of negative BOLD and CBF responses occurred therefore showed greatest mu power, an electrophysiological index of cortical inhibition, and largest somatosensory evoked potentials. Taken together, our results suggest that a neuronal mechanism underlies the NBR, but that the NBR may originate from a different neurovascular coupling mechanism to the PBR, suggesting that caution should be taken in assuming the NBR simply represents the neurophysiological inverse of the PBR.

Measurement of OEF and absolute CMRO2: MRI-based methods using interleaved and combined hypercapnia and hyperoxia

Blood oxygenation level dependent (BOLD) functional magnetic resonance imaging (fMRI) is most commonly used in a semi-quantitative manner to infer changes in brain activity. Despite the basis of the image contrast lying in the cerebral venous blood oxygenation level, quantification of absolute cerebral metabolic rate of oxygen consumption (CMRO2) has only recently been demonstrated. Here we examine two approaches to the calibration of fMRI signal to measure absolute CMRO2 using hypercapnic and hyperoxic respiratory challenges. The first approach is to apply hypercapnia and hyperoxia separately but interleaved in time and the second is a combined approach in which we apply hyperoxic challenges simultaneously with different levels of hypercapnia. Eleven healthy volunteers were studied at 3T using a dual gradient-echo spiral readout pulsed arterial spin labelling (ASL) imaging sequence. Respiratory challenges were conducted using an automated system of dynamic end-tidal forcing. A generalised BOLD signal model was applied, within a Bayesian estimation framework, that aims to explain the effects of modulation of CBF and arterial oxygen content to estimate venous deoxyhaemoglobin concentration ([dHb]0). Using CBF measurements combined with the estimated oxygen extraction fraction (OEF), absolute CMRO2 was calculated. The interleaved approach to hypercapnia and hyperoxia, as well as yielding estimates of CMRO2 and OEF demonstrated a significant increase in regional CBF, venous oxygen saturation (SvO2) (a decrease in OEF) and absolute CMRO2 in visual cortex in response to a continuous (20 min) visual task, demonstrating the potential for the method in measuring long term changes in CMRO2. The combined approach to oxygen and carbon dioxide modulation, as well as taking less time to acquire data, yielded whole brain grey matter estimates of CMRO2 and OEF of 184±45 μmol/100 g/min and 0.42±0.12 respectively, along with additional estimates of the vascular parameters α=0.33±0.06, the exponent relating relative increases in CBF to CBV, and β=1.35±0.13, the exponent relating deoxyhaemoglobin concentration to the relaxation rate R2*. Maps of cerebrovascular and cerebral metabolic parameters were also calculated. We show that combined modulation of oxygen and carbon dioxide can offer an experimentally more efficient approach to estimating OEF and absolute CMRO2 along with the additional vascular parameters that form an important part of the commonly used calibrated fMRI signal model.

Measuring cerebrovascular reactivity: what stimulus to use?

Cerebrovascular reactivity is the change in cerebral blood flow in response to a vasodilatory or vasoconstrictive stimulus. Measuring variations of cerebrovascular reactivity between different regions of the brain has the potential to not only advance understanding of how the cerebral vasculature controls the distribution of blood flow but also to detect cerebrovascular pathophysiology. While there are standardized and repeatable methods for estimating the changes in cerebral blood flow in response to a vasoactive stimulus, the same cannot be said for the stimulus itself. Indeed, the wide variety of vasoactive challenges currently employed in these studies impedes comparisons between them. This review therefore critically examines the vasoactive stimuli in current use for their ability to provide a standard repeatable challenge and for the practicality of their implementation. Such challenges include induced reductions in systemic blood pressure, and the administration of vasoactive substances such as acetazolamide and carbon dioxide. We conclude that many of the stimuli in current use do not provide a standard stimulus comparable between individuals and in the same individual over time. We suggest that carbon dioxide is the most suitable vasoactive stimulus. We describe recently developed computer-controlled MRI compatible gas delivery systems which are capable of administering reliable and repeatable vasoactive CO2 stimuli.

Using carbogen for calibrated fMRI at 7Tesla: comparison of direct and modelled estimation of the M parameter

Task-evoked changes in cerebral oxygen metabolism can be measured using calibrated functional Magnetic Resonance Imaging (fMRI). This technique requires the use of breathing manipulations such as hypercapnia, hyperoxia or a combination of both to determine a calibration factor M. The M-value is usually obtained by extrapolating the BOLD signal measured during the gas manipulation to its upper theoretical physiological limit using a biophysical model. However, a recently introduced technique uses a combination of increased inspired concentrations of O2 and CO2 to saturate the BOLD signal completely. In this study, we used this BOLD saturation technique to measure M directly at 7Tesla (T). Simultaneous carbogen-7 (7% CO2 in 93% O2) inhalation and visuo-motor task performance were used to elevate venous oxygen saturation in visual and motor areas close to their maximum, and the BOLD signal measured during this manipulation was used as an estimate of M. As accurate estimation of M is crucial for estimation of valid oxidative metabolism values, these directly estimated M-values were assessed and compared with M-values obtained via extrapolation modelling using the generalized calibration model (GCM) on the same dataset. Average M-values measured using both methods were 10.4±3.9% (modelled) and 7.5±2.2% (direct) for a visual-related ROI, and 11.3±5.2% (modelled) and 8.1±2.6% (direct) for a motor-related ROI. Results from this study suggest that, for the CO2 concentration used here, modelling is necessary for the accurate estimation of the M parameter. Neither gas inhalation alone, nor gas inhalation combined with a visuo-motor task, was sufficient to completely saturate venous blood in most subjects. Calibrated fMRI studies should therefore rely on existing models for gas inhalation-based calibration of the BOLD signal.

The physics of functional magnetic resonance imaging (fMRI)

Functional magnetic resonance imaging (fMRI) is a methodology for detecting dynamic patterns of activity in the working human brain. Although the initial discoveries that led to fMRI are only about 20 years old, this new field has revolutionized the study of brain function. The ability to detect changes in brain activity has a biophysical basis in the magnetic properties of deoxyhemoglobin, and a physiological basis in the way blood flow increases more than oxygen metabolism when local neural activity increases. These effects translate to a subtle increase in the local magnetic resonance signal, the blood oxygenation level dependent (BOLD) effect, when neural activity increases. With current techniques, this pattern of activation can be measured with resolution approaching 1 mm(3) spatially and 1 s temporally. This review focuses on the physical basis of the BOLD effect, the imaging methods used to measure it, the possible origins of the physiological effects that produce a mismatch of blood flow and oxygen metabolism during neural activation, and the mathematical models that have been developed to understand the measured signals. An overarching theme is the growing field of quantitative fMRI, in which other MRI methods are combined with BOLD methods and analyzed within a theoretical modeling framework to derive quantitative estimates of oxygen metabolism and other physiological variables. That goal is the current challenge for fMRI: to move fMRI from a mapping tool to a quantitative probe of brain physiology.

A review of calibrated blood oxygenation level-dependent (BOLD) methods for the measurement of task-induced changes in brain oxygen metabolism

The dynamics of the blood oxygenation level-dependent (BOLD) response are dependent on changes in cerebral blood flow, cerebral blood volume and the cerebral metabolic rate of oxygen consumption. Furthermore, the amplitude of the response is dependent on the baseline physiological state, defined by the haematocrit, oxygen extraction fraction and cerebral blood volume. As a result of this complex dependence, the accurate interpretation of BOLD data and robust intersubject comparisons when the baseline physiology is varied are difficult. The calibrated BOLD technique was developed to address these issues. However, the methodology is complex and its full promise has not yet been realised. In this review, the theoretical underpinnings of calibrated BOLD, and issues regarding this theory that are still to be resolved, are discussed. Important aspects of practical implementation are reviewed and reported applications of this methodology are presented.

A New Functional MRI Approach for Investigating Modulations of Brain Oxygen Metabolism

Functional MRI (fMRI) using the blood oxygenation level dependent (BOLD) signal is a common technique in the study of brain function. The BOLD signal is sensitive to the complex interaction of physiological changes including cerebral blood flow (CBF), cerebral blood volume (CBV), and cerebral oxygen metabolism (CMRO₂). A primary goal of quantitative fMRI methods is to combine BOLD imaging with other measurements (such as CBF measured with arterial spin labeling) to derive information about CMRO₂. This requires an accurate mathematical model to relate the BOLD signal to the physiological and hemodynamic changes; the most commonly used of these is the Davis model. Here, we propose a new nonlinear model that is straightforward and shows heuristic value in clearly relating the BOLD signal to blood flow, blood volume and the blood flow-oxygen metabolism coupling ratio. The model was tested for accuracy against a more detailed model adapted for magnetic fields of 1.5, 3 and 7T. The mathematical form of the heuristic model suggests a new ratio method for comparing combined BOLD and CBF data from two different stimulus responses to determine whether CBF and CMRO₂ coupling differs. The method does not require a calibration experiment or knowledge of parameter values as long as the exponential parameter describing the CBF-CBV relationship remains constant between stimuli. The method was found to work well for 1.5 and 3T but is prone to systematic error at 7T. If more specific information regarding changes in CMRO₂ is required, then with accuracy similar to that of the Davis model, the heuristic model can be applied to calibrated BOLD data at 1.5T, 3T and 7T. Both models work well over a reasonable range of blood flow and oxygen metabolism changes but are less accurate when applied to a simulated caffeine experiment in which CBF decreases and CMRO₂ increases.

Coupling of cerebral blood flow and oxygen metabolism is conserved for chromatic and luminance stimuli in human visual cortex

The ratio of the changes in cerebral blood flow (CBF) and cerebral metabolic rate of oxygen (CMRO(2)) during brain activation is a critical determinant of the magnitude of the blood oxygenation level dependent (BOLD) response measured with functional magnetic resonance imaging (fMRI). Cytochrome oxidase (CO), a key component of oxidative metabolism in the mitochondria, is non-uniformly distributed in visual area V1 in distinct blob and interblob regions, suggesting significant spatial variation in the capacity for oxygen metabolism. The goal of this study was to test whether CBF/CMRO(2) coupling differed when these subpopulations of neurons were preferentially stimulated, using chromatic and luminance stimuli to preferentially stimulate either the blob or interblob regions. A dual-echo spiral arterial spin labeling (ASL) technique was used to measure CBF and BOLD responses simultaneously in 7 healthy human subjects. When the stimulus contrast levels were adjusted to evoke similar CBF responses (mean 65.4% ± 19.0% and 64.6% ± 19.9%, respectively for chromatic and luminance contrast), the BOLD responses were remarkably similar (1.57% ± 0.39% and 1.59% ± 0.35%) for both types of stimuli. We conclude that CBF-CMRO(2) coupling is conserved for the chromatic and luminance stimuli used, suggesting a consistent coupling for blob and inter-blob neuronal populations despite the difference in CO concentration.

Luminance contrast of a visual stimulus modulates the BOLD response more than the cerebral blood flow response in the human brain

The blood oxygenation level dependent (BOLD) response measured with functional magnetic resonance imaging (fMRI) depends on the evoked changes in cerebral blood flow (CBF) and cerebral metabolic rate of oxygen (CMRO(2)) in response to changes in neural activity. This response is strongly modulated by the CBF/CMRO(2) coupling relationship with activation, defined as n, the ratio of the fractional changes. The reliability of the BOLD signal as a quantitative reflection of underlying physiological changes depends on the stability of n in response to different stimuli. The effect of visual stimulus contrast on this coupling ratio was tested in 9 healthy human subjects, measuring CBF and BOLD responses to a flickering checkerboard at four visual contrast levels. The theory of the BOLD effect makes a robust prediction-independent of details of the model-that if the CBF/CMRO(2) coupling ratio n remains constant, then the response ratio between the lowest and highest contrast levels should be higher for the BOLD response than the CBF response because of the ceiling effect on the BOLD response. Instead, this response ratio was significantly lower for the BOLD response (BOLD response: 0.23 ± 0.13, mean ± SD; CBF response: 0.42 ± 0.18; p=0.0054). This data is consistent with a reduced dynamic range (strongest/weakest response ratio) of the CMRO(2) response (~1.7-fold) compared to that of the CBF response (~2.4-fold) as luminance contrast increases, corresponding to an increase of n from 1.7 at the lowest contrast level to 2.3 at the highest contrast level. The implication of these results for fMRI studies is that the magnitude of the BOLD response does not accurately reflect the magnitude of underlying physiological processes.

Calibrated BOLD using direct measurement of changes in venous oxygenation

Calibration of the BOLD signal is potentially of great value in providing a closer measure of the underlying changes in brain function related to neuronal activity than the BOLD signal alone, but current approaches rely on an assumed relationship between cerebral blood volume (CBV) and cerebral blood flow (CBF). This is poorly characterised in humans and does not reflect the predominantly venous nature of BOLD contrast, whilst this relationship may vary across brain regions and depend on the structure of the local vascular bed. This work demonstrates a new approach to BOLD calibration which does not require an assumption about the relationship between cerebral blood volume and cerebral blood flow. This method involves repeating the same stimulus both at normoxia and hyperoxia, using hyperoxic BOLD contrast to estimate the relative changes in venous blood oxygenation and venous CBV. To do this the effect of hyperoxia on venous blood oxygenation has to be calculated, which requires an estimate of basal oxygen extraction fraction, and this can be estimated from the phase as an alternative to using a literature estimate. Additional measurement of the relative change in CBF, combined with the blood oxygenation change can be used to calculate the relative change in CMRO₂ due to the stimulus. CMRO₂ changes of 18 ± 8% in response to a motor task were measured without requiring the assumption of a CBV/CBF coupling relationship, and are in agreement with previous approaches.