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

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.

Laminar differences in functional oxygen metabolism in monkey visual cortex measured with calibrated fMRI

Blood-oxygenation-level-dependent functional magnetic resonance imaging (BOLD fMRI) of cortical layers relies on the hemodynamic response and is biased toward large veins on the cortical surface. Functional changes in the cerebral metabolic rate of oxygen (ΔCMRO2) may reflect neural cortical function better than BOLD fMRI, but it is unknown whether the calibrated BOLD model for functional CMRO2 measurement remains valid at high resolution. Here, we measure laminar ΔCMRO2 elicited by visual stimulation in macaque primary visual cortex (V1) and find that ΔCMRO2 peaks in the middle of the cortex, in agreement with autoradiographic measures of metabolism. ΔCMRO2 values in gray matter are similar as found previously. Reductions in CMRO2 are associated with veins at the cortical surface, suggesting that techniques for vein removal may improve the accuracy of the model at very high resolution. However, our results show feasibility of laminar ΔCMRO2 measurement, providing a physiologically meaningful metric of laminar functional metabolism.

Spatiotemporal features of neurovascular (un)coupling with stimulus-induced activity and hypercapnia challenge in cerebral cortex and olfactory bulb

Carbon dioxide (CO2) is traditionally considered as metabolic waste, yet its regulation is critical for brain function. It is well accepted that hypercapnia initiates vasodilation, but its effect on neuronal activity is less clear. Distinguishing how stimulus- and CO2-induced vasodilatory responses are (dis)associated with neuronal activity has profound clinical and experimental relevance. We used an optical method in mice to simultaneously image fluorescent calcium (Ca2+) transients from neurons and reflectometric hemodynamic signals during brief sensory stimuli (i.e., hindpaw, odor) and CO2 exposure (i.e., 5%). Stimuli-induced neuronal and hemodynamic responses swiftly increased within locally activated regions exhibiting robust neurovascular coupling. However, hypercapnia produced slower global vasodilation which was temporally uncoupled to neuronal deactivation. With trends consistent across cerebral cortex and olfactory bulb as well as data from GCaMP6f/jRGECO1a mice (i.e., green/red Ca2+ fluorescence), these results unequivocally reveal that stimuli and CO2 generate comparable vasodilatory responses but contrasting neuronal responses. In summary, observations of stimuli-induced regional neurovascular coupling and CO2-induced global neurovascular uncoupling call for careful appraisal when using CO2 in gas mixtures to affect vascular tone and/or neuronal excitability, because CO2 is both a potent vasomodulator and a neuromodulator.

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.

Modeling the dynamics of cerebrovascular reactivity to carbon dioxide in fMRI under task and resting-state conditions

Conventionally, cerebrovascular reactivity (CVR) is estimated as the amplitude of the hemodynamic response to vascular stimuli, most commonly carbon dioxide (CO₂). While the CVR amplitude has established clinical utility, the temporal characteristics of CVR (dCVR) have been increasingly explored and may yield even more pathology-sensitive parameters. This work is motivated by the current need to evaluate the feasibility of dCVR modeling in various experimental conditions. In this work, we present a comparison of several recently published/utilized model-based deconvolution (response estimation) approaches for estimating the CO₂ response function h(t), including maximum a posteriori likelihood (MAP), inverse logit (IL), canonical correlation analysis (CCA), and basis expansion (using Gamma and Laguerre basis sets). To aid the comparison, we devised a novel simulation framework that incorporates a wide range of SNRs, ranging from 10 to -7 dB, representative of both task and resting-state CO₂ changes. In addition, we built ground-truth h(t) into our simulation framework, overcoming the conventional limitation that the true h(t) is unknown. Moreover, to best represent realistic noise found in fMRI scans, we extracted noise from in-vivo resting-state scans. Furthermore, we introduce a simple optimization of the CCA method (CCA_opt) and compare its performance to these existing methods. Our findings suggest that model-based methods can accurately estimate dCVR even amidst high noise (i.e. resting-state), and in a manner that is largely independent of the underlying model assumptions for each method. We also provide a quantitative basis for making methodological choices, based on the desired dCVR parameters, the estimation accuracy and computation time. The BEL method provided the highest accuracy and robustness, followed by the CCA_opt and IL methods. Of the three, the CCA_opt method has the lowest computational requirements. These findings lay the foundation for wider adoption of dCVR estimation in CVR mapping.

Mapping oxidative metabolism in the human brain with calibrated fMRI in health and disease

Conventional functional MRI (fMRI) with blood-oxygenation level dependent (BOLD) contrast is an important tool for mapping human brain activity non-invasively. Recent interest in quantitative fMRI has renewed the importance of oxidative neuroenergetics as reflected by cerebral metabolic rate of oxygen consumption (CMR_O2) to support brain function. Dynamic CMR_O2 mapping by calibrated fMRI require multi-modal measurements of BOLD signal along with cerebral blood flow (CBF) and/or volume (CBV). In human subjects this "calibration" is typically performed using a gas mixture containing small amounts of carbon dioxide and/or oxygen-enriched medical air, which are thought to produce changes in CBF (and CBV) and BOLD signal with minimal or no CMR_O2 changes. However non-human studies have demonstrated that the "calibration" can also be achieved without gases, revealing good agreement between CMR_O2 changes and underlying neuronal activity (e.g., multi-unit activity and local field potential). Given the simpler set-up of gas-free calibrated fMRI, there is evidence of recent clinical applications for this less intrusive direction. This up-to-date review emphasizes technological advances for such translational gas-free calibrated fMRI experiments, also covering historical progression of the calibrated fMRI field that is impacting neurological and neurodegenerative investigations of the human brain.

Lag-Optimized Blood Oxygenation Level Dependent Cerebrovascular Reactivity Estimates Derived From Breathing Task Data Have a Stronger Relationship With Baseline Cerebral Blood Flow

Cerebrovascular reactivity (CVR), an important indicator of cerebrovascular health, is commonly studied with the Blood Oxygenation Level Dependent functional MRI (BOLD-fMRI) response to a vasoactive stimulus. Theoretical and empirical evidence suggests that baseline cerebral blood flow (CBF) modulates BOLD signal amplitude and may influence BOLD-CVR estimates. We address how acquisition and modeling choices affect the relationship between baseline cerebral blood flow (bCBF) and BOLD-CVR: whether BOLD-CVR is modeled with the inclusion of a breathing task, and whether BOLD-CVR amplitudes are optimized for hemodynamic lag effects. We assessed between-subject correlations of average GM values and within-subject spatial correlations across cortical regions. Our results suggest that a breathing task addition to a resting-state acquisition, alongside lag-optimization within BOLD-CVR modeling, can improve BOLD-CVR correlations with bCBF, both between- and within-subjects, likely because these CVR estimates are more physiologically accurate. We report positive correlations between bCBF and BOLD-CVR, both between- and within-subjects. The physiological explanation of this positive correlation is unclear; research with larger samples and tightly controlled vasoactive stimuli is needed. Insights into what drives variability in BOLD-CVR measurements and related measurements of cerebrovascular function are particularly relevant when interpreting results in populations with altered vascular and/or metabolic baselines or impaired cerebrovascular reserve.

Alterations of Cerebral Perfusion and Functional Connectivity in Children With Idiopathic Generalized Epilepsy

Background

Studies have demonstrated that adults with idiopathic generalized epilepsy (IGE) have functional abnormalities; however, the neuropathological pathogenesis differs between adults and children. This study aimed to explore alterations in the cerebral blood flow (CBF) and functional connectivity (FC) to comprehensively elucidate the neuropathological mechanisms of IGE in children.

Methods

We obtained arterial spin labeling (ASL) and resting state functional magnetic resonance imaging data of 28 children with IGE and 35 matched controls. We used ASL to determine differential CBF regions in children with IGE. A seed-based whole-brain FC analysis was performed for regions with significant CBF changes. The mean CBF and FC of brain areas with significant group differences was extracted, then its correlation with clinical variables in IGE group was analyzed by using Pearson correlation analysis.

Results

Compared to controls, children with IGE had CBF abnormalities that were mainly observed in the right middle temporal gyrus, right middle occipital gyrus (MOG), right superior frontal gyrus (SFG), left inferior frontal gyrus (IFG), and triangular part of the left IFG (IFGtriang). We observed that the FC between the left IFGtriang and calcarine fissure (CAL) and that between the right MOG and bilateral CAL were decreased in children with IGE. The CBF in the right SFG was correlated with the age at IGE onset. FC in the left IFGtriang and left CAL was correlated with the IGE duration.

Conclusion

This study found that CBF and FC were altered simultaneously in the left IFGtriang and right MOG of children with IGE. The combination of CBF and FC may provide additional information and insight regarding the pathophysiology of IGE from neuronal and vascular integration perspectives.

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.

The neural-vascular basis of age-related processing speed decline

Most studies examining neurocognitive aging are based on the blood-oxygen level-dependent signal obtained during functional magnetic resonance imaging (fMRI). The physiological basis of this signal is neural-vascular coupling, the process by which neurons signal cerebrovasculature to dilate in response to an increase in active neural metabolism due to stimulation. These fMRI studies of aging rely on the hemodynamic equivalence assumption that this process is not disrupted by physiologic deterioration associated with aging. Studies of neural-vascular coupling challenge this assumption and show that neural-vascular coupling is closely related to cognition. In this review, we put forward a theory of processing speed decline in aging and how it is related to age-related neural-vascular coupling changes based on the results of studies elucidating the relationships between cognition, cerebrovascular dynamics, and aging.

The Role of Cerebrovascular-Reactivity Mapping in Functional MRI: Calibrated fMRI and Resting-State fMRI

Task and resting-state functional MRI (fMRI) is primarily based on the same blood-oxygenation level-dependent (BOLD) phenomenon that MRI-based cerebrovascular reactivity (CVR) mapping has most commonly relied upon. This technique is finding an ever-increasing role in neuroscience and clinical research as well as treatment planning. The estimation of CVR has unique applications in and associations with fMRI. In particular, CVR estimation is part of a family of techniques called calibrated BOLD fMRI, the purpose of which is to allow the mapping of cerebral oxidative metabolism (CMRO2) using a combination of BOLD and cerebral-blood flow (CBF) measurements. Moreover, CVR has recently been shown to be a major source of vascular bias in computing resting-state functional connectivity, in much the same way that it is used to neutralize the vascular contribution in calibrated fMRI. Furthermore, due to the obvious challenges in estimating CVR using gas challenges, a rapidly growing field of study is the estimation of CVR without any form of challenge, including the use of resting-state fMRI for that purpose. This review addresses all of these aspects in which CVR interacts with fMRI and the role of CVR in calibrated fMRI, provides an overview of the physiological biases and assumptions underlying hypercapnia-based CVR and calibrated fMRI, and provides a view into the future of non-invasive CVR measurement.

Cerebral oxygen extraction fraction: Comparison of dual-gas challenge calibrated BOLD with CBF and challenge-free gradient echo QSM+qBOLD

PURPOSE

To compare cortical gray matter oxygen extraction fraction (OEF) estimated from 2 MRI methods: (1) the quantitative susceptibility mapping (QSM) plus quantitative blood oxygen level dependent imaging (qBOLD) (QSM+qBOLD or QQ), and (2) the dual-gas calibrated-BOLD (DGCB) in healthy subjects; and to investigate the validity of iso-cerebral metabolic rate of oxygen consumption assumption during hypercapnia using QQ.

METHODS

In 10 healthy subjects, 3 tesla MRI including a multi-echo gradient echo sequence at baseline and hypercapnia for QQ, as well as an EPI dual-echo pseudo-continuous arterial spin labeling for DGCB, were performed under a hypercapnic and a hyperoxic condition. OEFs from QQ and DGCB were compared using region of interest analysis and paired t test. For QQ, cerebral metabolic rate of oxygen consumption = cerebral blood flow*OEF*arterial oxygen content was generated for both baseline and hypercapnia, which were compared.

RESULTS

Average OEF in cortical gray matter across 10 subjects from QQ versus DGCB was 35.5 ± 6.7% versus 38.0 ± 9.1% (P = .49) at baseline and 20.7 ± 4.4% versus 28.4 ± 7.6% (P = .02) in hypercapnia: OEF in cortical gray matter was significantly reduced as measured in QQ (P < .01) and in DGCB (P < .01). Cerebral metabolic rate of oxygen consumption (in μmol O₂ /min/100 g) was 168.2 ± 54.1 at baseline from DGCB and was 153.1 ± 33.8 at baseline and 126.4 ± 34.2 (P < .01) in hypercapnia from QQ.

CONCLUSION

The differences in OEF obtained from QQ and DGCB are small and nonsignificant at baseline but are statistically significant during hypercapnia. In addition, QQ shows a cerebral metabolic rate of oxygen consumption decrease (17.4%) during hypercapnia.

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.

The impact of hyperoxia on brain activity: A resting-state and task-evoked electroencephalography (EEG) study

A better understanding of the effect of oxygen on brain electrophysiological activity may provide a more mechanistic insight into clinical studies that use oxygen treatment in pathological conditions, as well as in studies that use oxygen to calibrate functional magnetic resonance imaging (fMRI) signals. This study applied electroencephalography (EEG) in healthy subjects and investigated how high a concentration of oxygen in inhaled air (i.e., normobaric hyperoxia) alters brain activity under resting-state and task-evoked conditions. Study 1 investigated its impact on resting EEG and revealed that hyperoxia suppressed α (8-13Hz) and β (14-35Hz) band power (by 15.6±2.3% and 14.1±3.1%, respectively), but did not change the δ (1-3Hz), θ (4-7Hz), and γ (36-75Hz) bands. Sham control experiments did not result in such changes. Study 2 reproduced these findings, and, furthermore, examined the effect of hyperoxia on visual stimulation event-related potentials (ERP). It was found that the main peaks of visual ERP, specifically N1 and P2, were both delayed during hyperoxia compared to normoxia (P = 0.04 and 0.02, respectively). In contrast, the amplitude of the peaks did not show a change. Our results suggest that hyperoxia has a pronounced effect on brain neural activity, for both resting-state and task-evoked potentials.