Sources of systematic bias in hypercapnia-calibrated functional MRI estimation of oxygen metabolism

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.

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.

Quantitative relations between BOLD responses, cortical energetics and impulse firing across cortical depth

The blood oxygen level-dependent (BOLD) functional magnetic resonance imaging (fMRI) signal arises as a consequence of changes in cerebral blood flow (CBF) and cerebral metabolic rate of oxygen ( CMR O 2 ) that in turn are modulated by changes in neural activity. Recent advances in imaging have achieved sub-millimetre resolution and allowed investigation of the BOLD response as a function of cortical depth. Here, we adapt our previous theory relating the BOLD signal to neural activity to produce a quantitative model that incorporates venous blood draining between cortical layers. The adjustable inputs to the model are the neural activity and a parameter governing this blood draining. A three-layer version for transient neural inputs and a multi-layer version for constant or tonic neural inputs are able to account for a variety of experimental results, including negative BOLD signals.

Resting CMRO2 fluctuations show persistent network hyper-connectivity following exposure to sub-concussive collisions

Exposure to head impacts may alter brain connectivity within cortical hubs such as the default-mode network (DMN). However, studies have yet to consider the confounding effects of altered resting cerebral blood flow (CBF0) and cerebrovascular reactivity (CVR) on changes in connectivity following sub-concussive impacts. Here, 23 Canadian collegiate football players were followed during a season using calibrated resting-state MRI and helmet accelerometers to examine the interplay between the neural and vascular factors that determine functional connectivity (FC). Connectivity-based analyses using blood oxygen level dependent (BOLD) and cerebral metabolic rate of oxygen consumption (CMRO2) mapping were used to study the DMN longitudinally. Network-specific decreases in CBF0 were observed one month following the season, while impaired CVR was documented at both mid-season and one month following the season, compared to pre-season baseline. Alterations in CBF0 and BOLD-based CVR throughout the season suggest that neurophysiological markers may show different susceptibility timelines following head impacts. DMN connectivity was increased throughout the season, independent of changes in cerebrovascular physiology, suggesting that alterations in FC following sub-concussive impacts are robust and independent of changes in brain hemodynamics. No significant correlations between impact kinematics and DMN connectivity changes were documented in this study. Altogether, these findings create a strong paradigm for future studies to examine the underlying neural and vascular mechanisms associated with increases in network connectivity following repeated exposure to sub-concussive collisions, in an effort to improve management of head impacts in contact sports.

Quantitative relations between BOLD responses, cortical energetics, and impulse firing

The blood oxygen level-dependent (BOLD) functional magnetic resonance imaging signal arises as a consequence of changes in blood flow and oxygen usage that in turn are modulated by changes in neural activity. Much attention has been given to both theoretical and experimental aspects of the energetics but not to the neural activity. Here we identify the best energetic theory for the steady-state BOLD signal on the basis of correct predictions of experimental observations. This theory is then used, together with the recently determined relationship between energetics and neural activity, to predict how the BOLD signal changes with activity. Unlike existing treatments, this new theory incorporates a nonzero baseline activity in a completely consistent way and is thus able to account for both sustained positive and negative BOLD signals. We also show that the increase in BOLD signal for a given increase in activity is significantly smaller the larger the baseline activity, as is experimentally observed. Furthermore, the decline of the positive BOLD signal arising from deeper cortical laminae in response to an increase in neural firing is shown to arise as a consequence of the larger baseline activity in deeper laminae. Finally, we provide quantitative relations integrating BOLD responses, energetics, and impulse firing, which among other predictions give the same results as existing theories when the baseline activity is zero. NEW & NOTEWORTHY We use a recently established relation between energetics and neural activity to give a quantitative account of BOLD dependence on neural activity. The incorporation of a nonzero baseline neural activity accounts for positive and negative BOLD signals, shows that changes in neural activity give BOLD changes that are smaller the larger the baseline, and provides a basis for the observed inverse relation between BOLD responses and the depth of cortical laminae giving rise to them.

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

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.

Age-related changes in brain hemodynamics; A calibrated MRI study

INTRODUCTION

Blood oxygenation-level dependent (BOLD) magnetic resonance imaging signal changes in response to stimuli have been used to evaluate age-related changes in neuronal activity. Contradictory results from these types of experiments have been attributed to differences in cerebral blood flow (CBF) and cerebral metabolic rate of oxygen (CMRO2 ). To clarify the effects of these physiological parameters, we investigated the effect of age on baseline CBF and CMRO2 .

MATERIALS AND METHODS

Twenty young (mean ± sd age, 28 ± 3 years), and 45 older subjects (66 ± 4 years) were investigated. A dual-echo pseudocontinuous arterial spin labeling (ASL) sequence was performed during normocapnic, hypercapnic, and hyperoxic breathing challenges. Whole brain and regional gray matter values of CBF, ASL cerebrovascular reactivity (CVR), BOLD CVR, oxygen extraction fraction (OEF), and CMRO2 were calculated.

RESULTS

Whole brain CBF was 49 ± 14 and 40 ± 9 ml/100 g/min in young and older subjects respectively (P < 0.05). Age-related differences in CBF decreased to the point of nonsignificance (B=-4.1, SE=3.8) when EtCO2 was added as a confounder. BOLD CVR was lower in the whole brain, in the frontal, in the temporal, and in the occipital of the older subjects (P<0.05). Whole brain OEF was 43 ± 8% in the young and 39 ± 6% in the older subjects (P = 0.066). Whole brain CMRO2 was 181 ± 60 and 133 ± 43 µmol/100 g/min in young and older subjects, respectively (P<0.01).

DISCUSSION

Age-related differences in CBF could potentially be explained by differences in EtCO2 . Regional CMRO2 was lower in older subjects. BOLD studies should take this into account when investigating age-related changes in neuronal activity.

Regional variations in vascular density correlate with resting-state and task-evoked blood oxygen level-dependent signal amplitude

Functional magnetic resonance imaging (fMRI) has become one of the primary tools used for noninvasively measuring brain activity in humans. For the most part, the blood oxygen level-dependent (BOLD) contrast is used, which reflects the changes in hemodynamics associated with active brain tissue. The main advantage of the BOLD signal is that it is relatively easy to measure and thus is often used as a proxy for comparing brain function across population groups (i.e., control vs. patient). However, it is particularly weighted toward veins whose structural architecture is known to vary considerably across the brain. This makes it difficult to interpret whether differences in BOLD between cortical areas reflect true differences in neural activity or vascular structure. We therefore investigated how regional variations of vascular density (VAD) relate to the amplitude of resting-state and task-evoked BOLD signals. To address this issue, we first developed an automated method for segmenting veins in images acquired with susceptibility-weighted imaging, allowing us to visualize the venous vascular tree across the brain. In 19 healthy subjects, we then applied voxel-based morphometry (VBM) to T1-weighted images and computed regional measures of gray matter density (GMD). We found that, independent of spatial scale, regional variations in resting-state and task-evoked fMRI amplitudes were better correlated to VAD compared to GMD. Using a general linear model (GLM), it was observed that the bulk of regional variance in resting-state activity could be modeled by VAD. Cortical areas whose resting-state activity was most suppressed by VAD correction included Cuneus, Precuneus, Culmen, and BA 9, 10, and 47. Taken together, our results suggest that resting-state BOLD signals are significantly related to the underlying structure of the brain vascular system. Calibrating resting BOLD activity by venous structure may result in a more accurate interpretation of differences observed between cortical areas and/or individuals.

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.

Neural mechanisms of age-related slowing: the ΔCBF/ΔCMRO2 ratio mediates age-differences in BOLD signal and human performance

The precise mechanisms that give rise to the blood-oxygen-level-dependent (BOLD) activation differences that accompany age-related cognitive slowing remain fundamentally unknown. We sought to isolate the origin of age-related BOLD changes by comparing blood-flow and oxygen-metabolic constituents of the BOLD response using dual-echo arterial spin labeling during visual stimulation and CO2 ingestion. We hypothesized, and our results confirmed, that age-related changes in the ratio of fractional cerebral blood flow to fractional cerebral metabolic rate of oxygen consumption (ΔCBF/ΔCMRO2) lead to the BOLD changes that are observed in older adults. ΔCBF/ΔCMRO2 was also significantly related to performance, suggesting that age-related cognitive slowing results from neural cell assemblies that operate less efficiently, requiring greater oxygen metabolism that is not matched by blood-flow changes relative to younger adults. Age-related changes in ΔCBF/ΔCMRO2 are sufficient to explain variations in BOLD responding and performance cited throughout the literature, assuming no bias based on physiological baseline CMRO2.

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.

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.

Poststimulus undershoots in cerebral blood flow and BOLD fMRI responses are modulated by poststimulus neuronal activity

fMRI is the foremost technique for noninvasive measurement of human brain function. However, its utility is limited by an incomplete understanding of the relationship between neuronal activity and the hemodynamic response. Though the primary peak of the hemodynamic response is modulated by neuronal activity, the origin of the typically negative poststimulus signal is poorly understood and its amplitude assumed to covary with the primary response. We use simultaneous recordings of EEG with blood oxygenation level-dependent (BOLD) and cerebral blood flow (CBF) fMRI during unilateral median nerve stimulation to show that the poststimulus fMRI signal is neuronally modulated. We observe high spatial agreement between concurrent BOLD and CBF responses to median nerve stimulation, with primary signal increases in contralateral sensorimotor cortex and primary signal decreases in ipsilateral sensorimotor cortex. During the poststimulus period, the amplitude and directionality (positive/negative) of the BOLD signal in both contralateral and ipsilateral sensorimotor cortex depends on the poststimulus synchrony of 8-13 Hz EEG neuronal activity, which is often considered to reflect cortical inhibition, along with concordant changes in CBF and metabolism. Therefore we present conclusive evidence that the fMRI time course represents a hemodynamic signature of at least two distinct temporal phases of neuronal activity, substantially improving understanding of the origin of the BOLD response and increasing the potential measurements of brain function provided by fMRI. We suggest that the poststimulus EEG and fMRI responses may be required for the resetting of the entire sensory network to enable a return to resting-state activity levels.

A BOLD Perspective on Age-Related Neurometabolic-Flow Coupling and Neural Efficiency Changes in Human Visual Cortex

Age-related performance declines in visual tasks have been attributed to reductions in processing efficiency. The neural basis of these declines has been explored by comparing the blood-oxygen-level-dependent (BOLD) index of neural activity in older and younger adults during visual task performance. However, neural activity is one of many factors that change with age and lead to BOLD signal differences. We investigated the origin of age-related BOLD changes by comparing blood flow and oxygen metabolic constituents of BOLD signal. Subjects periodically viewed flickering annuli and pressed a button when detecting luminance changes in a central fixation cross. Using magnetic resonance dual-echo arterial spin labeling and CO2 ingestion, we observed age-equivalent (i.e., similar in older and younger groups) fractional cerebral blood flow (ΔCBF) in the presence of age-related increases in fractional cerebral metabolic rate of oxygen (ΔCMRO2). Reductions in ΔCBF responsiveness to increased ΔCMRO2 in elderly led to paradoxical age-related BOLD decreases. Age-related ΔCBF/ΔCMRO2 ratio decreases were associated with reaction times, suggesting that age-related slowing resulted from less efficient neural activity. We hypothesized that reduced vascular responsiveness to neural metabolic demand would lead to a reduction in ΔCBF/ΔCMRO2. A simulation of BOLD relative to ΔCMRO2 for lower and higher neurometabolic-flow coupling ratios (approximating those for old and young, respectively) indicated less BOLD signal change in old than young in relatively lower CMRO2 ranges, as well as greater BOLD signal change in young compared to old in relatively higher CMRO2 ranges. These results suggest that age-comparative studies relying on BOLD signal might be misinterpreted, as age-related BOLD changes do not merely reflect neural activity changes. Age-related declines in neurometabolic-flow coupling might lead to neural efficiency reductions that can adversely affect visual task performance.

Relationships between hypercarbic reactivity, cerebral blood flow, and arterial circulation times in patients with moyamoya disease

PURPOSE

To evaluate the correlation between angiographic measures of Moyamoya disease and tissue-level impairment from measurements of tissue perfusion and cerebrovascular reactivity (CVR).

MATERIALS AND METHODS

The relationship between perfusion-weighted arterial spin labeling (ASL) and hypercarbic blood oxygenation-level dependent (BOLD) CVR and time-to-peak (TTP) were compared with angiographically measured risk factors, including arterial circulation time (ACT) and modified Suzuki Score (mSS) in patients (n = 15) with Moyamoya disease.

RESULTS

Hemodynamic contrasts provided information not apparent from structural or angiographic imaging. Mean z-statistics demonstrate that BOLD is significantly (P = 0.017) higher in low mSS hemispheres (z-statistic = 5.0 ± 2.5) compared with high mSS hemispheres (z-statistic = 3.7 ± 1.7), suggesting that regions with less advanced stages of Moyamoya disease have higher reactivity. After correcting for multiple comparisons, a strong trend for a direct relationship (R = 0.38; P = 0.03) between BOLD TTP and ACT was observed, and a significant inverse relationship between CBF and ACT (R = -0.47; P = 0.01) was found, demonstrating that BOLD and ASL contrasts reflect DSA measures of vascular compromise in Moyamoya disease, albeit with different sensitivity.

CONCLUSION

Correlative measures between angiography and hemodynamic methods suggest that BOLD and ASL could be used for expanding the diagnostic imaging infrastructure in Moyamoya patients and potentially tracking tissue response to revascularization.

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.

Calibrating the BOLD signal during a motor task using an extended fusion model incorporating DOT, BOLD and ASL data

Multimodal imaging improves the accuracy of the localization and the quantification of brain activation when measuring different manifestations of the hemodynamic response associated with cerebral activity. In this study, we incorporated cerebral blood flow (CBF) changes measured with arterial spin labeling (ASL), Diffuse Optical Tomography (DOT) and blood oxygen level-dependent (BOLD) recordings to reconstruct changes in oxy- (ΔHbO(2)) and deoxyhemoglobin (ΔHbR). Using the Grubb relation between relative changes in CBF and cerebral blood volume (CBV), we incorporated the ASL measurement as a prior to the total hemoglobin concentration change (ΔHbT). We applied this ASL fusion model to both synthetic data and experimental multimodal recordings during a 2-s finger-tapping task. Our results show that the new approach is very powerful in estimating ΔHbO(2) and ΔHbR with high spatial and quantitative accuracy. Moreover, our approach allows the computation of baseline total hemoglobin concentration (HbT(0)) as well as of the BOLD calibration factor M on a single subject basis. We obtained an average HbT(0) of 71 μM, an average M value of 0.18 and an average increase of 13% in cerebral metabolic rate of oxygen (CMRO(2)), all of which are in agreement with values previously reported in the literature. Our method yields an independent measurement of M, which provides an alternative measurement to validate the hypercapnic calibration of the BOLD signal.

Evolution of the dynamic changes in functional cerebral oxidative metabolism from tissue mitochondria to blood oxygen

The dynamic properties of the cerebral metabolic rate of oxygen consumption (CMR(O2)) during changes in brain activity remain unclear. Therefore, the spatial and temporal evolution of functional increases in CMR(O2) was investigated in the rat somato-sensory cortex during forelimb stimulation under a suppressed blood flow response condition. Temporally, stimulation elicited a fast increase in tissue mitochondria CMR(O2) described by a time constant of ~1 second measured using flavoprotein autofluorescence imaging. CMR(O2)-driven changes in the tissue oxygen tension measured using an oxygen electrode and blood oxygenation measured using optical imaging of intrinsic signal followed; however, these changes were slow with time constants of ~5 and ~10 seconds, respectively. This slow change in CMR(O2)-driven blood oxygenation partly explains the commonly observed post-stimulus blood oxygen level-dependent (BOLD) undershoot. Spatially, the changes in mitochondria CMR(O2) were similar to the changes in blood oxygenation. Finally, the increases in CMR(O2) were well correlated with the evoked multi-unit spiking activity. These findings show that dynamic CMR(O2) calculations made using only blood oxygenation data (e.g., BOLD functional magnetic resonance imaging (fMRI)) do not directly reflect the temporal changes in the tissue's mitochondria metabolic rate; however, the findings presented can bridge the gap between the changes in cellular oxidative rate and blood oxygenation.

Indication of BOLD-specific venous flow-volume changes from precisely controlled hyperoxic vs. hypercapnic calibration

Deriving cerebral metabolic rate of oxygen consumption (CMRO(2)) from blood oxygenation level-dependent (BOLD) signals involves a flow-volume parameter (α), reflecting total cerebral blood volume changes, and a calibration constant (M). Traditionally, the former is assumed a fixed value and the latter is measured under alterations in fixed inspired fractional concentrations of carbon dioxide. We recently reported on reductions in M-variability via precise control of end-tidal pressures of both hypercapnic (HC) and hyperoxic (HO) gases. In light of these findings, our aim was to apply the improved calibration alternatives to neuronal activation, making use of their distinct vasoactive natures to evaluate the α-value. Nine healthy volunteers were imaged at 3 T while simultaneously measuring BOLD and arterial spin-labeling signals during controlled, graded, HC, and HO, followed by visual (VC) and sensorimotor cortices (SMC) activation. On the basis of low M- and CMRO(2)-variability, the comparison of these calibration alternatives accurately highlighted a reduced venous flow-volume relationship (α=0.16±0.02, with α(VC)=0.12±0.04, and α(SMC)=0.20±0.02), as appropriate for BOLD modeling.

A comparison of physiologic modulators of fMRI signals

One of the main obstacles in quantitative interpretation of functional magnetic resonance imaging (fMRI) signal is that this signal is influenced by non-neural factors such as vascular properties of the brain, which effectively increases signal variability. One approach to account for non-neural components is to identify and measure these confounding factors and to include them as covariates in data analysis or interpretation. Previously, several research groups have independently identified four potential physiologic modulators of fMRI signals, including baseline venous oxygenation (Yv ), cerebrovascular reactivity (CVR), resting state BOLD fluctuation amplitude (RSFA), and baseline cerebral blood flow (CBF). This study sought to directly compare the modulation effects of these indices in the same fMRI session. The physiologic parameters were measured with techniques comparable with those used in the previous studies except for CBF, which was determined globally with a velocity-based phase-contrast MRI (instead of arterial-spin-labeling MRI). Using an event-related, scene-categorization fMRI task, we showed that the fMRI signal amplitude was positively correlated with CVR (P < 0.0001) and RSFA (P = 0.002), while negatively correlated with baseline Yv (P < 0.0001). The fMRI-CBF correlation did not reach significance, although the (negative) sign of the correlation was consistent with the earlier study. Furthermore, among the physiologic modulators themselves, significant correlations were observed between baseline Yv and baseline CBF (P = 0.01), and between CVR and RSFA (P = 0.05), suggesting that some of the modulators may partly be of similar physiologic origins. These observations as well as findings in recent literature suggest that additional measurement of physiologic modulator(s) in an fMRI session may provide a practical approach to control for inter-subject variations and to improve the ability of fMRI in detecting disease or medication related differences. Hum Brain Mapp 34:2078-2088, 2013. © 2011 Wiley Periodicals, Inc.

Characterizing brain oxygen metabolism in patients with multiple sclerosis with T2-relaxation-under-spin-tagging MRI

In this study, venous oxygen saturation and oxygen metabolic changes in multiple sclerosis (MS) patients were assessed using a recently developed T2-relaxation-under-spin-tagging (TRUST) magnetic resonance imaging (MRI), which measures the superior sagittal venous sinus blood oxygenation (Yv) and cerebral metabolic rate of oxygen (CMRO(2)), an index of global oxygen consumption. Thirty patients with relapsing-remitting MS and 30 age-matched healthy controls were studied using TRUST at 3 T MR. The mean expanded disability status scale (EDSS) of the patients was 2.3 (range, 0 to 5.5). We found significantly increased Yv (P<0.0001) and decreased CMRO(2) (P=0.003) in MS patients (mean±s.d.: 65.9%±5.1% and 138.8±35.4 μmol per 100 g per minute) as compared with healthy control subjects (60.2%±4.0% and 180.2±24.8 μmol per 100 g per minute, respectively), implying decrease of oxygen consumption in MS. There was a significant positive correlation between Yv and EDSS and between Yv and lesion load in MS patients (n=30); on the contrary, there was a significant negative correlation between CMRO(2) and EDSS and between CMRO(2) and lesion load (n=12). There was no correlation between Yv and brain atrophy measures. This study showed preliminary evidence of the potential utility of TRUST in global oxygen metabolism. Our results of significant underutilization of oxygen in MS raise important questions regarding mitochondrial respiratory dysfunction and neurodegeneration of the disease.

Calibrated FMRI

Functional magnetic resonance imaging with blood oxygenation level-dependent (BOLD) contrast has had a tremendous influence on human neuroscience in the last twenty years, providing a non-invasive means of mapping human brain function with often exquisite sensitivity and detail. However the BOLD method remains a largely qualitative approach. While the same can be said of anatomic MRI techniques, whose clinical and research impact has not been diminished in the slightest by the lack of a quantitative interpretation of their image intensity, the quantitative expression of BOLD responses as a percent of the baseline T2*- weighted signal has been viewed as necessary since the earliest days of fMRI. Calibrated MRI attempts to dissociate changes in oxygen metabolism from changes in blood flow and volume, the latter three quantities contributing jointly to determine the physiologically ambiguous percent BOLD change. This dissociation is typically performed using a "calibration" procedure in which subjects inhale a gas mixture containing small amounts of carbon dioxide or enriched oxygen to produce changes in blood flow and BOLD signal which can be measured under well-defined hemodynamic conditions. The outcome is a calibration parameter M which can then be substituted into an expression providing the fractional change in oxygen metabolism given changes in blood flow and BOLD signal during a task. The latest generation of calibrated MRI methods goes beyond fractional changes to provide absolute quantification of resting-state oxygen consumption in micromolar units, in addition to absolute measures of evoked metabolic response. This review discusses the history, challenges, and advances in calibrated MRI, from the personal perspective of the author.

A generalized procedure for calibrated MRI incorporating hyperoxia and hypercapnia

Calibrated MRI techniques use the changes in cerebral blood flow (CBF) and blood oxygenation level-dependent (BOLD) signal evoked by a respiratory manipulation to extrapolate the total BOLD signal attributable to deoxyhemoglobin at rest (M). This parameter can then be used to estimate changes in the cerebral metabolic rate of oxygen consumption (CMRO(2)) based on task-induced BOLD and CBF signals. Different approaches have been described previously, including addition of inspired CO(2) (hypercapnia) or supplemental O(2) (hyperoxia). We present here a generalized BOLD signal model that reduces under appropriate conditions to previous models derived for hypercapnia or hyperoxia alone, and is suitable for use during hybrid breathing manipulations including simultaneous hypercapnia and hyperoxia. This new approach yields robust and accurate M maps, in turn allowing more reliable estimation of CMRO(2) changes evoked during a visual task. The generalized model is valid for arbitrary flow changes during hyperoxia, thus benefiting from the larger total oxygenation changes produced by increased blood O(2) content from hyperoxia combined with increases in flow from hypercapnia. This in turn reduces the degree of extrapolation required to estimate M. The new procedure yielded M estimates that were generally higher (7.6 ± 2.6) than those obtained through hypercapnia (5.6 ± 1.8) or hyperoxia alone (4.5 ± 1.5) in visual areas. These M values and their spatial distribution represent a more accurate and robust depiction of the underlying distribution of tissue deoxyhemoglobin at rest, resulting in more accurate estimates of evoked CMRO(2) changes.

Quantitative measurement of cerebral physiology using respiratory-calibrated MRI

Functional magnetic resonance imaging typically measures signal increases arising from changes in the transverse relaxation rate over small regions of the brain and associates these with local changes in cerebral blood flow, blood volume and oxygen metabolism. Recent developments in pulse sequences and image analysis methods have improved the specificity of the measurements by focussing on changes in blood flow or changes in blood volume alone. However, FMRI is still unable to match the physiological information obtainable from positron emission tomography (PET), which is capable of quantitative measurements of blood flow and volume, and can indirectly measure resting metabolism. The disadvantages of PET are its cost, its availability, its poor spatial resolution and its use of ionising radiation. The MRI techniques introduced here address some of these limitations and provide physiological data comparable with PET measurements. We present an 18-minute MRI protocol that produces multi-slice whole-brain coverage and yields quantitative images of resting cerebral blood flow, cerebral blood volume, oxygen extraction fraction, CMRO(2), arterial arrival time and cerebrovascular reactivity of the human brain in the absence of any specific functional task. The technique uses a combined hyperoxia and hypercapnia paradigm with a modified arterial spin labelling sequence.

The Relationship Between M in "Calibrated fMRI" and the Physiologic Modulators of fMRI

The “calibrated fMRI” technique requires a hypercapnia calibration experiment in order to estimate the factor “M”. It is desirable to be able to obtain the M value without the need of a gas challenge calibration. According to the analytical expression of M, it is a function of several baseline physiologic parameters, such as baseline venous oxygenation and CBF, both of which have recently been shown to be significant modulators of fMRI signal. Here we studied the relationship among hypercapnia-calibrated M, baseline venous oxygenation and CBF, and assessed the possibility of estimating M from the baseline physiologic parameters. It was found that baseline venous oxygenation and CBF are highly correlated (R²=0.77, P<0.0001) across subjects. However, the hypercapnia-calibrated M was not correlated with baseline venous oxygenation or CBF. The hypercapnia-calibrated M was not correlated with an estimation of M based on analytical expression either. The lack of correlation may be explained by the counteracting effect of venous oxygenation and CBF on the M factor, such that the actual M value of an individual may be mostly dependent on other parameters such as hematocrit. Potential biases in hypercapnia-based M estimation were also discussed in the context of possible reduction of CMRO₂ during hypercapnia.

Vascular component analysis of hyperoxic and hypercapnic BOLD contrast

Hyperoxia or hypercapnia provides a useful experimental tool to systematically alter the blood oxygenation level dependent (BOLD) contrast. Typical applications include calibrated functional magnetic resonance imaging (fMRI), BOLD sensitivity mapping, vessel size imaging or cerebrovascular reactivity mapping. This article describes a novel biophysical model of hyperoxic and hypercapnic BOLD contrast, which accounts for the magnetic susceptibility effects of molecular oxygen that is dissolved in blood and tissue, in addition to the well-established effects caused by the paramagnetic properties of deoxyhaemoglobin. Furthermore, the concept of vascular component analysis (VCA) is introduced and is shown to provide a computationally efficient tool for investigating the vascular specificity of hyperoxic and hypercapnic BOLD contrast. A theoretical investigation of gradient and spin echo BOLD contrast based on computer simulations was performed to compare three different conditions (hypercapnia induced by breathing 6% CO2, hyperoxia induced by breathing 100% O2, and simultaneous hypercapnia and hyperoxia induced by breathing carbogen, i.e. 5% CO2 in 95% CO2) with baseline (breathing air). Simulations were carried out for different levels of metabolic oxygen extraction fraction (OEF) ranging from 0 to 0.5. The key findings can be summarised as follows: (i) for hyperoxia the susceptibility of dissolved O2 may lead to a significant arterial BOLD contrast; (ii) under normoxic conditions the susceptibility of dissolved O2 is negligible; (iii) an almost complete loss of BOLD sensitivity may occur at lower OEF values in all parts of the vascular tree, whereas hyperoxic BOLD sensitivity is largely maintained; (iv) under hyperoxic conditions, a transition from positive to negative BOLD contrast occurs with decreasing OEF values. These findings have important implications for experimental applications of hyperoxic and hypercapnic BOLD contrast and may enable new clinical applications in ischemic stroke and other forms of acquired brain injury.

A pathophysiological framework of hippocampal dysfunction in ageing and disease

The hippocampal formation has been implicated in a growing number of disorders, from Alzheimer's disease and cognitive ageing to schizophrenia and depression. How can the hippocampal formation, a complex circuit that spans the temporal lobes, be involved in a range of such phenotypically diverse and mechanistically distinct disorders? Recent neuroimaging findings indicate that these disorders differentially target distinct subregions of the hippocampal circuit. In addition, some disorders are associated with hippocampal hypometabolism, whereas others show evidence of hypermetabolism. Interpreted in the context of the functional and molecular organization of the hippocampal circuit, these observations give rise to a unified pathophysiological framework of hippocampal dysfunction.

HIV and chronic methamphetamine dependence affect cerebral blood flow

Human immunodeficiency virus (HIV) and methamphetamine (METH) dependence are independently associated with neuronal dysfunction. The coupling between cerebral blood flow (CBF) and neuronal activity is the basis of many task-based functional neuroimaging techniques. We examined the interaction between HIV infection and a previous history of METH dependence on CBF within the lenticular nuclei (LN). Twenty-four HIV-/METH-, eight HIV-/METH+, 24 HIV+/METH-, and 15 HIV+/METH+ participants performed a finger tapping paradigm. A multiple regression analysis of covariance assessed associations and two-way interactions between CBF and HIV serostatus and/or previous history of METH dependence. HIV+ individuals had a trend towards a lower baseline CBF (-10%, p = 0.07) and greater CBF changes for the functional task (+32%, p = 0.01) than HIV- subjects. Individuals with a previous history of METH dependence had a lower baseline CBF (-16%, p = 0.007) and greater CBF changes for a functional task (+33%, p = 0.02). However, no interaction existed between HIV serostatus and previous history of METH dependence for either baseline CBF (p = 0.53) or CBF changes for a functional task (p = 0.10). In addition, CBF and volume in the LN were not correlated. A possible additive relationship could exist between HIV infection and a history of METH dependence on CBF with a previous history of METH dependence having a larger contribution. Abnormalities in CBF could serve as a surrogate measure for assessing the chronic effects of HIV and previous METH dependence on brain function.

Depression of cortical activity in humans by mild hypercapnia

The effects of neural activity on cerebral hemodynamics underlie human brain imaging with functional magnetic resonance imaging and positron emission tomography. However, the threshold and characteristics of the converse effects, wherein the cerebral hemodynamic and metabolic milieu influence neural activity, remain unclear. We tested whether mild hypercapnia (5% CO2 ) decreases the magnetoencephalogram response to auditory pattern recognition and visual semantic tasks. Hypercapnia induced statistically significant decreases in event-related fields without affecting behavioral performance. Decreases were observed in early sensory components in both auditory and visual modalities as well as later cognitive components related to memory and language. Effects were distributed across cortical regions. Decreases were comparable in evoked versus spontaneous spectral power. Hypercapnia is commonly used with hemodynamic models to calibrate the blood oxygenation level-dependent response. Modifying model assumptions to incorporate the current findings produce a modest but measurable decrease in the estimated cerebral metabolic rate for oxygen change with activation. Because under normal conditions, low cerebral pH would arise when bloodflow is unable to keep pace with neuronal activity, the cortical depression observed here may reflect a homeostatic mechanism by which neuronal activity is adjusted to a level that can be sustained by available bloodflow. Animal studies suggest that these effects may be mediated by pH-modulating presynaptic adenosine receptors. Although the data is not clear, comparable changes in cortical pH to those induced here may occur during sleep apnea, sleep, and exercise. If so, these results suggest that such activities may in turn have generalized depressive effects on cortical activity.

The influence of carbon dioxide on brain activity and metabolism in conscious humans

A better understanding of carbon dioxide (CO(2)) effect on brain activity may have a profound impact on clinical studies using CO(2) manipulation to assess cerebrovascular reserve and on the use of hypercapnia as a means to calibrate functional magnetic resonance imaging (fMRI) signal. This study investigates how an increase in blood CO(2), via inhalation of 5% CO(2), may alter brain activity in humans. Dynamic measurement of brain metabolism revealed that mild hypercapnia resulted in a suppression of cerebral metabolic rate of oxygen (CMRO(2)) by 13.4% ± 2.3% (N=14) and, furthermore, the CMRO(2) change was proportional to the subject's end-tidal CO(2) (Et-CO(2)) change. When using functional connectivity MRI (fcMRI) to assess the changes in resting-state neural activity, it was found that hypercapnia resulted in a reduction in all fcMRI indices assessed including cluster volume, cross-correlation coefficient, and amplitude of the fcMRI signal in the default-mode network (DMN). The extent of the reduction was more pronounced than similar indices obtained in visual-evoked fMRI, suggesting a selective suppression effect on resting-state neural activity. Scalp electroencephalogram (EEG) studies comparing hypercapnia with normocapnia conditions showed a relative increase in low frequency power in the EEG spectra, suggesting that the brain is entering a low arousal state on CO(2) inhalation.

Simulation study of brain blood flow regulation by intra-cortical arterioles in an anatomically accurate large human vascular network. Part II: flow variations induced by global or localized modifications of arteriolar diameters

In a companion paper (Lorthois et al., Neuroimage, in press), we perform the first simulations of blood flow in an anatomically accurate large human intra-cortical vascular network (~10000 segments), using a 1D non-linear model taking into account the complex rheological properties of blood flow in microcirculation. This model predicts blood pressure, blood flow and hematocrit distributions, volumes of functional vascular territories, regional flow at voxel and network scales, etc. Using the same approach, we study flow reorganizations induced by global arteriolar vasodilations (an isometabolic global increase in cerebral blood flow). For small to moderate global vasodilations, the relationship between changes in volume and changes in flow is in close agreement with Grubb's law, providing a quantitative tool for studying the variations of its exponent with underlying vascular architecture. A significant correlation between blood flow and vascular structure at the voxel scale, practically unchanged with respect to baseline, is demonstrated. Furthermore, the effects of localized arteriolar vasodilations, representative of a local increase in metabolic demand, are analyzed. In particular, localized vasodilations induce flow changes, including vascular steal, in the neighboring arteriolar trunks at small distances (<300 μm), while their influence in the neighboring veins is much larger (about 1 mm), which provides an estimate of the vascular point spread function. More generally, for the first time, the hemodynamic component of various functional neuroimaging techniques has been isolated from metabolic and neuronal components, and a direct relationship with several known characteristics of the BOLD signal has been demonstrated.

Test-retest stability of calibrated BOLD-fMRI in HIV- and HIV+ subjects

Subject performance, scanner hardware, or biological factors can affect single session neuroimaging measures. Stability studies using calibrated blood oxygenation level-dependent functional magnetic resonance imaging (BOLD-fMRI) have been performed in health but not disease. We utilized calibrated BOLD-fMRI to determine the effects of HIV on neurovascular coupling. Six clinically stable HIV-infected patients (HIV+) and 10 seronegative controls (HIV-) were scanned at two separate sessions approximately 3 months apart. Both mild hypercapnia (5% CO(2)) exposure and a visual functional activation task were performed. Intraclass correlation coefficients (ICC) and inter-subject variance were determined for calibrated BOLD-fMRI measures (baseline cerebral blood flow (CBF), functional CBF, BOLD, and cerebral metabolic rate of oxygen consumption (CMRO(2)) changes) for HIV+ and HIV- subjects. The two groups did not differ in age, sex, or education. HIV+ subjects had lower mean baseline CBF (p<0.04, Cohen's d=-1.07) and functional BOLD responses (p<0.001, Cohen's d=-2.47) and a trend towards a decrease in mean functional CBF responses (p=0.07, Cohen's d=-0.92) despite similar mean functional CMRO(2) changes (p=0.71, Cohen's d=0.19). The stability of each calibrated BOLD-fMRI measure, as assessed by ICC, was significantly lower for HIV+ subjects. In addition, HIV+ participants had greater inter-subject variability for baseline CBF (p<0.02), functional BOLD (p<0.001), CBF (p<0.001), and CMRO(2) (p<0.002) responses. Our results demonstrate that calibrated BOLD-fMRI measures have excellent stability within healthy controls. In contrast, these values have greater variability in clinically stable HIV+ subjects and may reflect alterations in coupling between CBF and CMRO(2) with disease.

Improved fMRI calibration: precisely controlled hyperoxic versus hypercapnic stimuli

The calibration of functional magnetic resonance imaging (fMRI) for the estimation of neuronal activation-induced changes in cerebral metabolic rate of oxygen (CMRO(2)) has been achieved through hypercapnic-induced iso-metabolic increases in cerebral blood flow (CBF). Hypercapnia (HC) has been traditionally implemented through alterations in the fixed inspired fractional concentrations of carbon dioxide (F(I)CO(2)) without otherwise controlling end-tidal partial pressures of carbon dioxide (P(ET)CO(2)) or oxygen (P(ET)O(2)). There are several shortcomings to the use of this manual HC method that may be improved by using precise targeting of P(ET)CO(2) while maintaining iso-oxia. Similarly, precise control of blood gases can be used to induce isocapnic hyperoxia (HO) to reduce venous deoxyhaemoglobin (dHb) and thus increase BOLD signals, without appreciably altering CMRO(2) or CBF. The aim of our study was to use precise end-tidal targeting to compare the calibration of BOLD signals under an isocapnic hyperoxic protocol (HOP) (rises in P(ET)O(2) to 140, 240 and 340 mm Hg from baseline) to that of an iso-oxic hypercapnic protocol (HCP) (rises in P(ET)CO(2) of 3, 5, 7 and 9 mm Hg from baseline). Nine healthy volunteers were imaged at 3T while monitoring end-tidal gas concentrations and simultaneously measuring BOLD and CBF signals, via arterial spin labeling (ASL), during graded HCP and HOP, alternating with normocapnic states in a blocked experimental design. The variability of the calibration constant obtained under HOP (M(HOP)) was 0.3-0.5 that of the HCP one (M(HCP)). In addition, M-variances with precise gas targeting (M(HCP) and M(HOP)) were less than those reported in studies using traditional F(I)CO(2) and F(I)O(2) methods (M(HC) and M(HO), respectively). We conclude that precise controlled gas delivery markedly improves BOLD-calibration for fMRI studies of oxygen metabolism with both the HCP and the more precise HOP-alternative.

Interpreting oxygenation-based neuroimaging signals: the importance and the challenge of understanding brain oxygen metabolism

Functional magnetic resonance imaging is widely used to map patterns of brain activation based on blood oxygenation level dependent (BOLD) signal changes associated with changes in neural activity. However, because oxygenation changes depend on the relative changes in cerebral blood flow (CBF) and cerebral metabolic rate of oxygen (CMRO(2)), a quantitative interpretation of BOLD signals, and also other functional neuroimaging signals related to blood or tissue oxygenation, is fundamentally limited until we better understand brain oxygen metabolism and how it is related to blood flow. However, the positive side of the complexity of oxygenation signals is that when combined with dynamic CBF measurements they potentially provide the best tool currently available for investigating the dynamics of CMRO(2). This review focuses on the problem of interpreting oxygenation-based signals, the challenges involved in measuring CMRO(2) in general, and what is needed to put oxygenation-based estimates of CMRO(2) on a firm foundation. The importance of developing a solid theoretical framework is emphasized, both as an essential tool for analyzing oxygenation-based multimodal measurements, and also potentially as a way to better understand the physiological phenomena themselves. The existing data, integrated within a simple theoretical framework of O(2) transport, suggests the hypothesis that an important functional role of the mismatch of CBF and CMRO(2) changes with neural activation is to prevent a fall of tissue pO(2). Future directions for better understanding brain oxygen metabolism are discussed.

Quantification of CMRO(2) without hypercapnia using simultaneous near-infrared spectroscopy and fMRI measurements

Estimation of the cerebral metabolic rate of oxygen (CMRO(2)) and cerebral blood flow (CBF) is important to investigate the neurovascular coupling and physiological components in blood oxygenation level-dependent (BOLD) signals quantitatively. Although there are methods that can determine CMRO(2) changes using functional MRI (fMRI) or near-infrared spectroscopy (NIRS), current approaches require a separate hypercapnia calibration process and have the potential to incur bias in many assumed model parameters. In this paper, a novel method to estimate CMRO(2) without hypercapnia is described using simultaneous measurements of NIRS and fMRI. Specifically, an optimization framework is proposed that minimizes the differences between the two forms of the relative CMRO(2)-CBF coupling ratio from BOLD and NIRS biophysical models, from which hypercapnia calibration and model parameters are readily estimated. Based on the new methods, we found that group average CBF, CMRO(2), cerebral blood volume (CBV), and BOLD changes within activation of the primary motor cortex during a finger tapping task increased by 39.5 +/- 21.4%, 18.4 +/- 8.7%, 12.9 +/- 6.7%, and 0.5 +/- 0.2%, respectively. The group average estimated flow-metabolism coupling ratio was 2.38 +/- 0.65 and the hypercapnia parameter was 7.7 +/- 1.7%. These values are within the range of values reported from other literatures. Furthermore, the activation maps from CBF and CMRO(2) were well localized on the primary motor cortex, which is the main target region of the finger tapping task.

Improving fMRI sensitivity by normalization of basal physiologic state

The power of fMRI in assessing neural activities is hampered by inter-subject variations in basal physiologic parameters, which may not be related to neural activation but has a modulatory effect on fMRI signals. Therefore, normalization of fMRI signals with these parameters is useful in reducing variations and improving sensitivity of this important technique. Recently, we have shown that basal venous oxygenation is a significant modulator of fMRI signals and individuals with higher venous oxygenation tend to have lower fMRI signals. In this study, we aim to test the utility of venous oxygenation normalization in distinguishing subject groups. A "model" condition was used in which two visual stimuli with different flashing frequencies were used to stimulate two subject groups, respectively, thereby simulating the situation of control and patient groups. It was found that visual-evoked BOLD signal is significantly correlated with baseline venous T2 (P = 0.0003) and inclusion of physiologic modulator in the regression analysis can substantially reduce P values of group-level statistical tests. When applied to voxel-wise analysis, the normalization process can allow the detection of more significant voxels. The utility of other basal parameters, including blood pressure, heart rate, arterial oxygenation, and end-tidal CO(2), in BOLD normalization was also assessed and it was found that the improvement was less significant. Time-to-peak of the BOLD responses was also studied and it was found that subjects with higher basal venous oxygenation tend to slower BOLD responses.

Noninvasive quantification of whole-brain cerebral metabolic rate of oxygen (CMRO2) by MRI

Cerebral metabolic rate of oxygen (CMRO(2)) is an important marker for brain function and brain health. Existing techniques for quantification of CMRO(2) with positron emission tomography (PET) or MRI involve special equipment and/or exogenous agents, and may not be suitable for routine clinical studies. In the present study, a noninvasive method is developed to estimate whole-brain CMRO(2) in humans. This method applies phase-contrast MRI for quantitative blood flow measurement and T(2)-relaxation-under-spin-tagging (TRUST) MRI for venous oxygenation estimation, and uses the Fick principle of arteriovenous difference for the calculation of CMRO(2). Whole-brain averaged CMRO(2) values in young, healthy subjects were 132.1 +/- 20.0 micromol/100 g/min, in good agreement with literature reports using PET. Various acquisition strategies for phase-contrast and TRUST MRI were compared, and it was found that nongated phase-contrast and sagittal sinus (SS) TRUST MRI were able to provide the most efficient and accurate estimation of CMRO(2). In addition, blood flow and venous oxygenation were found to be positively correlated across subjects. Owing to the noninvasive nature of this method, it may be a convenient and useful approach for assessment of brain metabolism in brain disorders as well as under various physiologic conditions.

On the assessment of cerebrovascular reactivity using hypercapnia BOLD MRI

Cerebrovascular reactivity (CVR) reflects the capacity of blood vessels to dilate and is an important marker for brain vascular reserve. It may provide a useful addition to the traditional baseline blood flow measurement when assessing vascular factors in brain disorders. Blood-oxygenation-level-dependent MRI under CO(2) inhalation offers a non-invasive and quantitative means to estimate CVR in humans. In this study, we investigated several important methodological aspects of this technique with the goal of optimizing the experimental and data processing strategies for clinical use. Comparing 4 min of 5% CO(2) inhalation (less comfortable) to a 1 min inhalation (more comfortable) duration, it was found that the CVR values were 0.31 +/- 0.05%/mmHg (N = 11) and 0.31 +/- 0.08%/mmHg (N = 9), respectively, showing no significant differences between the two breathing paradigms. Therefore, the 1 min paradigm is recommended for future application studies for patient comfort and tolerability. Furthermore, we have found that end-tidal CO(2) recording was useful for accurate quantification of CVR because it provided both timing and amplitude information regarding the input function to the brain vascular system, which can be subject-dependent. Finally, we show that inter-subject variations in CVR are of physiologic origin and affect the whole brain in a similar fashion. Based on this, it is proposed that relative CVR (normalized against the CVR of the whole brain or a reference tissue) may be a more sensitive biomarker than absolute CVR in clinical applications as it minimizes inter-subject variations. With these technological optimizations, CVR mapping may become a useful method for studies of neurological and psychiatric diseases.

Sensitivity of neural-hemodynamic coupling to alterations in cerebral blood flow during hypercapnia

The relationship between measurements of cerebral blood oxygenation and neuronal activity is highly complex and depends on both neurovascular and neurometabolic biological coupling. While measurements of blood oxygenation changes via optical and MRI techniques have been developed to map functional brain activity, there is evidence that the specific characteristics of these signals are sensitive to the underlying vascular physiology and structure of the brain. Since baseline blood flow and oxygen saturation may vary between sessions and across subjects, functional blood oxygenation changes may be a less reliable indicator of brain activity in comparison to blood flow and metabolic changes. In this work, we use a biomechanical model to examine the relationships between neural, vascular, metabolic, and hemodynamic responses to parametric whisker stimulation under both normal and hypercapnic conditions in a rat model. We find that the relationship between neural activity and oxy- and deoxyhemoglobin changes is sensitive to hypercapnia-induced changes in baseline cerebral blood flow. In contrast, the underlying relationships between evoked neural activity, blood flow, and model-estimated oxygen metabolism changes are unchanged by the hypercapnic challenge. We conclude that evoked changes in blood flow and cerebral oxygen metabolism are more closely associated with underlying evoked neuronal responses.

Effects of aging on cerebral blood flow, oxygen metabolism, and blood oxygenation level dependent responses to visual stimulation

Calibrated functional magnetic resonance imaging (fMRI) provides a noninvasive technique to assess functional metabolic changes associated with normal aging. We simultaneously measured both the magnitude of the blood oxygenation level dependent (BOLD) and cerebral blood flow (CBF) responses in the visual cortex for separate conditions of mild hypercapnia (5% CO(2)) and a simple checkerboard stimulus in healthy younger (n = 10, mean: 28-years-old) and older (n = 10, mean: 53-years-old) adults. From these data we derived baseline CBF, the BOLD scaling parameter M, the fractional change in the cerebral metabolic rate of oxygen consumption (CMRO(2)) with activation, and the coupling ratio n of the fractional changes in CBF and CMRO(2). For the functional activation paradigm, the magnitude of the BOLD response was significantly lower for the older group (0.57 +/- 0.07%) compared to the younger group (0.95 +/- 0.14%), despite the finding that the fractional CBF and CMRO(2) changes were similar for both groups. The weaker BOLD response for the older group was due to a reduction in the parameter M, which was significantly lower for older (4.6 +/- 0.4%) than younger subjects (6.5 +/- 0.8%), most likely reflecting a reduction in baseline CBF for older (41.7 +/- 4.8 mL/100 mL/min) compared to younger (59.6 +/- 9.1 mL/100 mL/min) subjects. In addition to these primary responses, for both groups the BOLD response exhibited a post-stimulus undershoot with no significant difference in this magnitude. However, the post-undershoot period of the CBF response was significantly greater for older compared to younger subjects. We conclude that when comparing two populations, the BOLD response can provide misleading reflections of underlying physiological changes. A calibrated approach provides a more quantitative reflection of underlying metabolic changes than the BOLD response alone.

Mapping functional connectivity based on synchronized CMRO2 fluctuations during the resting state

Synchronized low-frequency fluctuations in the resting state functional MRI (fMRI) signal have been suggested to be associated with functional connectivity in brain networks. However, the underlying mechanism of this connectivity is still poorly understood, with the synchronized fluctuations could either originate from hemodynamic oscillations or represent true neuronal signaling. To better interpret the resting signal, in the current work, we examined spontaneous fluctuations at the level of cerebral metabolic rate of oxygenation (CMRO2), an index reflecting regional oxygen consumption and metabolism, and thus less sensitive to vascular dynamics. The CMRO2 signal was obtained based on a biophysical model with data acquired from simultaneous blood oxygenation level dependent (BOLD) and perfusion signals. CMRO2-based functional connectivity maps were generated in three brain networks: visual, default-mode, and hippocampus. Experiments were performed on twelve healthy participants during 'resting state' and as a comparison, with a visual task. CMRO2 signals in each of the above mentioned brain networks showed significant correlations. Functional connectivity maps from the CMRO2 signal are, in general, similar to those from BOLD and perfusion. In addition, we demonstrated that the three parameters (M, alpha and beta) in the biophysical model for calculating CMRO2 have negligible effects on the determination of the CMRO2-based connectivity strength. This study provides evidence that the spontaneous fluctuations in fMRI at rest likely originate from dynamic changes of cerebral metabolism reflecting neuronal activity.