Method for rapid MRI quantification of global cerebral metabolic rate of oxygen

Recent Advances in MR Imaging-based Quantification of Brain Oxygen Metabolism

The metabolic rate of oxygen (MRO2) is fundamental to tissue metabolism. Determination of MRO2 demands knowledge of the arterio-venous difference in hemoglobin-bound oxygen concentration, typically expressed as oxygen extraction fraction (OEF), and blood flow rate (BFR). MRI is uniquely suited for measurement of both these quantities, yielding MRO2 in absolute physiologic units of µmol O2 min-1/100 g tissue. Two approaches are discussed, both relying on hemoglobin magnetism. Emphasis will be on cerebral oxygen metabolism expressed in terms of the cerebral MRO2 (CMRO2), but translation of the relevant technologies to other organs, including kidney and placenta will be touched upon as well. The first class of methods exploits the blood's bulk magnetic susceptibility, which can be derived from field maps. The second is based on measurement of blood water T2, which is modulated by diffusion and exchange in the local-induced fields within and surrounding erythrocytes. Some whole-organ methods achieve temporal resolution adequate to permit time-series studies of brain energetics, for instance, during sleep in the scanner with concurrent electroencephalogram (EEG) sleep stage monitoring. Conversely, trading temporal for spatial resolution has led to techniques for spatially resolved approaches based on quantitative blood oxygen level dependent (BOLD) or calibrated BOLD models, allowing regional assessment of vascular-metabolic parameters, both also exploiting deoxyhemoglobin paramagnetism like their whole-organ counterparts.

Superior sagittal sinus flow as a proxy for tracking global cerebral blood flow dynamics during wakefulness and sleep

Sleep, a state of reduced consciousness, affects brain oxygen metabolism and lowers cerebral metabolic rate of oxygen (CMRO2). Previously, we quantified CMRO2 during sleep via Fick's Principle, with a single-band MRI sequence measuring both hemoglobin O2 saturation (SvO2) and superior sagittal sinus (SSS) blood flow, which was upscaled to obtain total cerebral blood flow (tCBF). The procedure involves a brief initial calibration scan to determine the upscaling factor (fc), assumed state-invariant. Here, we used a dual-band sequence to simultaneously provide SvO2 in SSS and tCBF in the neck every 16 seconds, allowing quantification of fc dynamically. Ten healthy subjects were scanned by MRI with simultaneous EEG for 80 minutes, yielding 300 temporal image frames per subject. Four volunteers achieved slow-wave sleep (SWS), as evidenced by increased δ-wave activity (per American Academy of Sleep Medicine criteria). SWS was maintained for 13.5 ± 7.0 minutes, with CMRO2 28.6 ± 5.5% lower than pre-sleep wakefulness. Importantly, there was negligible bias between tCBF obtained by upscaling SSS-blood flow, and tCBF measured directly in the inflowing arteries of the neck (intra-class correlation 0.95 ± 0.04, averaged across all subjects), showing that the single-band approach is a valid substitute for quantifying tCBF, simplifying image data collection and analysis without sacrificing accuracy.

Cerebral oxygen metabolism from MRI susceptibility

This article provides an overview of MRI methods exploiting magnetic susceptibility properties of blood to assess cerebral oxygen metabolism, including the tissue oxygen extraction fraction (OEF) and the cerebral metabolic rate of oxygen (CMRO2). The first section is devoted to describing blood magnetic susceptibility and its effect on the MRI signal. Blood circulating in the vasculature can have diamagnetic (oxyhemoglobin) or paramagnetic properties (deoxyhemoglobin). The overall balance between oxygenated and deoxygenated hemoglobin determines the induced magnetic field which, in turn, modulates the transverse relaxation decay of the MRI signal via additional phase accumulation. The following sections of this review then illustrate the principles underpinning susceptibility-based techniques for quantifying OEF and CMRO2. Here, it is detailed whether these techniques provide global (OxFlow) or local (Quantitative Susceptibility Mapping - QSM, calibrated BOLD - cBOLD, quantitative BOLD - qBOLD, QSM+qBOLD) measurements of OEF or CMRO2, and what signal components (magnitude or phase) and tissue pools they consider (intravascular or extravascular). Validations studies and potential limitations of each method are also described. The latter include (but are not limited to) challenges in the experimental setup, the accuracy of signal modeling, and assumptions on the measured signal. The last section outlines the clinical uses of these techniques in healthy aging and neurodegenerative diseases and contextualizes these reports relative to results from gold-standard PET.

A flow-diffusion model of oxygen transport for quantitative mapping of cerebral metabolic rate of oxygen (CMRO2) with single gas calibrated fMRI

One promising approach for mapping CMRO₂ is dual-calibrated functional MRI (dc-fMRI). This method exploits the Fick Principle to combine estimates of CBF from ASL, and OEF derived from BOLD-ASL measurements during arterial O₂ and CO₂ modulations. Multiple gas modulations are required to decouple OEF and deoxyhemoglobin-sensitive blood volume. We propose an alternative single gas calibrated fMRI framework, integrating a model of oxygen transport, that links blood volume and CBF to OEF and creates a mapping between the maximum BOLD signal, CBF and OEF (and CMRO₂). Simulations demonstrated the method's viability within physiological ranges of mitochondrial oxygen pressure, P_mO₂, and mean capillary transit time. A dc-fMRI experiment, performed on 20 healthy subjects using O₂ and CO₂ challenges, was used to validate the approach. The validation conveyed expected estimates of model parameters (e.g., low P_mO₂), with spatially uniform OEF maps (grey matter, GM, OEF spatial standard deviation ≈ 0.13). GM OEF estimates obtained with hypercapnia calibrated fMRI correlated with dc-fMRI (r = 0.65, p = 2·10^-3). For 12 subjects, OEF measured with dc-fMRI and the single gas calibration method were correlated with whole-brain OEF derived from phase measures in the superior sagittal sinus (r = 0.58, p = 0.048; r = 0.64, p = 0.025 respectively). Simplified calibrated fMRI using hypercapnia holds promise for clinical application.

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.

A Noninvasive Method for Quantifying Cerebral Metabolic Rate of Oxygen by Hybrid PET/MRI: Validation in a Porcine Model

The gold standard for imaging the cerebral metabolic rate of oxygen (CMRO2) is positron emission tomography (PET); however, it is an invasive and complex procedure that also requires correction for recirculating 15O-H2O and the blood-borne activity. We propose a noninvasive reference-based hybrid PET/magnetic resonance imaging (MRI) method that uses functional MRI techniques to calibrate 15O-O2-PET data. Here, PET/MR imaging of oxidative metabolism (PMROx) was validated in an animal model by comparison to PET-alone measurements. Additionally, we investigated if the MRI-perfusion technique arterial spin labelling (ASL) could be used to further simplify PMROx by replacing 15O-H2O-PET, and if the PMROx was sensitive to anesthetics-induced changes in metabolism. Methods: 15O-H2O and 15O-O2 PET data were acquired in a hybrid PET/MR scanner (3 T Siemens Biograph mMR), together with simultaneous functional MRI (OxFlow and ASL), from juvenile pigs (n = 9). Animals were anesthetized with 3% isoflurane and 6 mL/kg/h propofol for the validation experiments and arterial sampling was performed for PET-alone measurements. PMROx estimates were obtained using whole-brain (WB) CMRO2 from OxFlow and local cerebral blood flow (CBF) from either noninvasive 15O-H2O-PET or ASL (PMROxASL). Changes in metabolism were investigated by increasing the propofol infusion to 20 mL/kg/h. Results: Good agreement and correlation were observed between regional CMRO2 measurements from PMROx and PET-alone. No significant differences were found between OxFlow and PET-only measurements of WB oxygen extraction fraction (0.30 ± 0.09 and 0.31 ± 0.09) and CBF (54.1 ± 16.7 and 56.6 ± 21.0 mL/100 g/min), or between PMROx and PET-only CMRO2 estimates (1.89 ± 0.16 and 1.81 ± 0.10 mLO2/100 g/min). Moreover, PMROx and PMROxASL were sensitive to propofol-induced reduction in CMRO2 Conclusion: This study provides initial validation of a noninvasive PET/MRI technique that circumvents many of the complexities of PET CMRO2 imaging. PMROx does not require arterial sampling and has the potential to reduce PET imaging to 15O-O2 only; however, future validation involving human participants are required.

A non-invasive reference-based method for imaging the cerebral metabolic rate of oxygen by PET/MR: theory and error analysis

Positron emission tomography (PET) remains the gold standard for quantitative imaging of the cerebral metabolic rate of oxygen (CMRO₂); however, it is an invasive and complex procedure that requires accounting for recirculating [¹⁵O]H₂O (RW) and the cerebral blood volume (CBV). This study presents a non-invasive reference-based technique for imaging CMRO₂ that was developed for PET/magnetic resonance imaging (MRI) with the goal of simplifying the PET procedure while maintaining its ability to quantify metabolism. The approach is to use whole-brain (WB) measurements of oxygen extraction fraction (OEF) and cerebral blood flow (CBF) to calibrate [¹⁵O]O₂-PET data, thereby avoiding the need for invasive arterial sampling. Here we present the theoretical framework, along with error analyses, sensitivity to PET noise and inaccuracies in input parameters, and initial assessment on PET data acquired from healthy participants. Simulations showed that neglecting RW and CBV corrections caused errors in CMRO₂ of less than ±10% for changes in regional OEF of ±25%. These predictions were supported by applying the reference-based approach to PET data, which resulted in remarkably similar CMRO₂ images to those generated by analyzing the same data using a modeling approach that incorporated the arterial input functions and corrected for CBV contributions. Significant correlations were observed between regional CMRO₂ values from the two techniques (slope = 1.00 ± 0.04, R ² > 0.98) with no significant differences found for integration times of 3 and 5 min. In summary, results demonstrate the feasibility of producing quantitative CMRO₂ images by PET/MRI without the need for invasive blood sampling.

Measurements of venous oxygen saturation in the superior sagittal sinus using conventional 3D multiple gradient-echo MRI: Effects of flow velocity and acceleration

PURPOSE

This work investigates the effects of flow acceleration in the superior sagittal sinus on slice-dependent variations in venous oxygen saturation (SvO₂ ) estimations using susceptibility-based MR oximetry.

METHODS

Three-dimensional multiple gradient-echo images, with first-order flow compensation along the anterior-posterior readout direction for the first echo, were acquired twice from 15 healthy volunteers. For all slices, phases within the superior sagittal sinus were fitted using linear regression across four TEs to obtain the Pearson's correlation coefficients (PCCs), the largest of which corresponded to minimum acceleration influence. SvO₂ derived from odd echoes on this slice was used to assess interscan difference, and compared with the central 15th slice for slice-dependent difference, both using Bland-Altman analysis. Within-scan interslice SvO₂ consistency was examined versus PCC. Multislice-averaged SvO₂ values were then computed from slices with PCCs above a certain threshold.

RESULTS

Slice-dependent difference in SvO₂ varied from -16.2% to 21.5% at two SDs, in agreement with a recent report, and about twice larger than interscan differences for the automatically selected slice (-7.5% to 10.3%) and for the central 15th slice (-8.0% to 8.8%). For slices with PCCs higher than -0.98, interslice SvO₂ deviations were all found to be less than 5.0%. Multislice-averaged SvO₂ with PCCs higher than -0.98 further reduced interscan difference to -4.7% to 8.2%.

CONCLUSION

Slice-dependent variations in SvO₂ may partly be explained by the effects of flow acceleration. Our method may enable conventional 3D multiple gradient echo to be used for SvO₂ estimations in the presence of pulsatile flow.

In vivo validation of T2- and susceptibility-based Sv O2 measurements with jugular vein catheterization under hypoxia and hypercapnia

PURPOSE

To investigate the mutual agreement of T2-based and susceptibility-based methods as well as their agreement with jugular catheterization, for quantifying venous oxygen saturation (S_v O₂ ) at a broad range of brain oxygenation levels.

METHODS

S_v O₂ measurements using T2-relaxation-under-spin-tagging (TRUST) and susceptibility-based oximetry (SBO) were performed in 13 healthy subjects under room air, hypoxia, and hypercapnia conditions. Agreement between TRUST and SBO was quantitatively evaluated. In two of the subjects, TRUST and SBO were compared against the clinical gold standard, co-oximeter measurement via internal jugular vein catheterization.

RESULTS

Absolute S_v O₂ measurements using TRUST and SBO were highly correlated across a range of saturations from 45% to 84% (Pearson r = 0.91, P < .0001). S_v O₂ -TRUST was significantly lower than S_v O₂ -SBO under hypoxia and room air conditions, but the two were comparable under hypercapnia. TRUST demonstrated a larger S_v O₂ increase under hypercapnia than SBO and had good agreement with jugular catheterization under hypercapnia but significantly underestimated S_v O₂ under room air and hypoxia. The agreement between S_v O₂ -SBO and the reference did not depend on the physiological state.

CONCLUSION

A systematic bias was observed between T2-based and susceptibility-based methods that depended on the oxygenation state. In vivo validation with jugular catheterization indicated potential underestimation of TRUST under room air and hypoxia conditions. Our findings suggested that caution should be employed in comparison of absolute S_v O₂ measurements using either TRUST or SBO.

T2 -prepared balanced steady-state free precession (bSSFP) for quantifying whole-blood oxygen saturation at 1.5T

PURPOSE

To establish a calibration equation to convert human blood T2 to the full range of oxygen saturation levels (HbO2 ) and physiologic hematocrit (Hct) values using a T2 -prepared balanced steady-state free precession sequence (T2 -SSFP) at 1.5T.

METHODS

Blood drawn from 10 healthy donors (29.1 ± 3.9 years old) was prepared into samples of varying HbO2 and Hct (n = 79), and imaged using T2 -SSFP sequence at 37°C and interrefocusing interval τ180  = 12 ms. The relationship between blood T2 , HbO2 , and Hct was established based on the model R2=R2,plasma+Hct (R2,RBC-R2,plasma)+k·Hct·(1-Hct)·(1-HbO2)2. Measured R2 and HbO2 levels were fit by the model yielding values of R2,plasma, R2,RBC, and k. T2 -SSFP and the established calibration equation were applied to extract HbO2 at the superior sagittal sinus (SSS) in vivo and were compared with susceptometry-based oximetry.

RESULTS

Constants derived from the fit were: k = 74.2 [s-1 ], R2,plasma  = 1.5 [s-1 ], R2,RBC  = 11.6 [s-1 ], the R2 of the fit was 0.95. Average HbO2 at the SSS in seven healthy volunteers was 65% ± 7% and 66% ± 7% via T2 - and susceptometry-based oximetry, respectively. Bland-Altman analysis indicated agreement between the two oximetric methods with no significant bias.

CONCLUSION

The calibration constants presented here should ensure improved accuracy for whole-blood oximetry based on T2 -SSFP at 1.5T. Magn Reson Med 79:1893-1900, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Cerebral venous collaterals: A new fort for fighting ischemic stroke?

Stroke therapy has entered a new era highlighted by the use of endovascular therapy in addition to intravenous thrombolysis. However, the efficacy of current therapeutic regimens might be reduced by their associated adverse events. For example, over-reperfusion and futile recanalization may lead to large infarct, brain swelling, hemorrhagic complication and neurological deterioration. The traditional pathophysiological understanding on ischemic stroke can hardly address these occurrences. Accumulating evidence suggests that a functional cerebral venous drainage, the major blood reservoir and drainage system in brain, may be as critical as arterial infusion for stroke evolution and clinical sequelae. Further exploration of the multi-faceted function of cerebral venous system may add new implications for stroke outcome prediction and future therapeutic decision-making. In this review, we emphasize the anatomical and functional characteristics of the cerebral venous system and illustrate its necessity in facilitating the arterial infusion and maintaining the cerebral perfusion in the pathological stroke content. We then summarize the recent critical clinical studies that underscore the associations between cerebral venous collateral and outcome of ischemic stroke with advanced imaging techniques. A novel three-level venous system classification is proposed to demonstrate the distinct characteristics of venous collaterals in the setting of ischemic stroke. Finally, we discuss the current directions for assessment of cerebral venous collaterals and provide future challenges and opportunities for therapeutic strategies in the light of these new concepts.

Susceptibility-based time-resolved whole-organ and regional tissue oximetry

The magnetism of hemoglobin - being paramagnetic in its deoxy and diamagnetic in its oxy state - offers unique opportunities to probe oxygen metabolism in blood and tissues. The magnetic susceptibility χ of blood scales linearly with blood oxygen saturation, which can be obtained by measuring the magnetic field ΔB of the intravascular MR signal relative to tissue. In contrast to χ, the induced field ΔB is non-local. Therefore, to obtain the intravascular susceptibility Δχ relative to adjoining tissue from the measured ΔB demands solution of an inverse problem. Fortunately, for ellipsoidal structures, to which a straight, cylindrically shaped blood vessel segment conforms, the solution is trivial. The article reviews the principle of MR susceptometry-based blood oximetry. It then discusses applications for quantification of whole-brain oxygen extraction - typically on the basis of a measurement in the superior sagittal sinus - and, in conjunction with total cerebral blood flow, the cerebral metabolic rate of oxygen (CMRO₂ ). By simultaneously measuring flow and venous oxygen saturation (SvO₂ ) a temporal resolution of a few seconds can be achieved, allowing the study of the response to non-steady-state challenges such as volitional apnea. Extensions to regional measurements in smaller cerebral veins are also possible, as well as voxelwise quantification of venous blood saturation in cerebral veins accomplished by quantitative susceptibility mapping (QSM) techniques. Applications of susceptometry-based oximetry to studies of metabolic and degenerative disorders of the brain are reviewed. Lastly, the technique is shown to be applicable to other organ systems such as the extremities using SvO₂ as a dynamic tracer to monitor the kinetics of the microvascular response to induced ischemia. Copyright © 2016 John Wiley & Sons, Ltd.

Multiplexed MRI methods for rapid estimation of global cerebral metabolic rate of oxygen consumption

The global cerebral metabolic rate of oxygen (CMRO₂), which reflects metabolic activity of the brain under various physiologic conditions, can be quantified using a method, referred to as 'OxFlow', which simultaneously measures hemoglobin oxygen saturation in a draining vein (Y_v) and total cerebral blood flow (tCBF). Conventional OxFlow (Conv-OxFlow) entails four interleaves incorporated in a single pulse sequence - two for phase-contrast based measurement of tCBF in the supplying arteries of the neck, and two to measure the intra- to extravascular phase difference in the superior sagittal sinus to derive Y_v [Jain et al., JCBFM 2010]. However, this approach limits achievable temporal resolution thus precluding capture of rapid changes of brain metabolic states such as the response to apneic stimuli. Here, we developed a time-efficient, multiplexed OxFlow method and evaluated its potential for measuring dynamic alterations in global CMRO₂ during a breath-hold challenge. Two different implementations of multiplexed OxFlow were investigated: 1) simultaneous-echo-refocusing based OxFlow (SER-OxFlow) and 2) simultaneous-multi-slice imaging-based dual-band OxFlow (DB-OxFlow). The two sequences were implemented on 3T scanners (Siemens TIM Trio and Prisma) and their performance was evaluated in comparison to Conv-OxFlow in ten healthy subjects for baseline CMRO₂ quantification. Comparison of measured parameters (Y_v, tCBF, CMRO₂) revealed no significant bias of SER-OxFlow and DB-OxFlow, with respect to the reference Conv-OxFlow while improving temporal resolution two-fold (12.5 versus 25s). Further acceleration shortened scan time to 8 and 6s for SER and DB-OxFlow, respectively, for time-resolved CMRO₂ measurement. The two sequences were able of capturing smooth transitions of Y_v, tCBF, and CMRO₂ over the time course consisting of 30s of normal breathing, 30s of volitional apnea, and 90s of recovery. While both SER- and DB-OxFlow techniques provide significantly improved temporal resolution (by a factor of 3 - 4 relative to Conv-OxFlow), DB-OxFlow was found to be superior for the study of short physiologic stimuli.

Cerebral Hemodynamics in Patients with Hemolytic Uremic Syndrome Assessed by Susceptibility Weighted Imaging and Four-Dimensional Non-Contrast MR Angiography

Background and Purpose

Conventional magnetic resonance imaging (MRI) of patients with hemolytic uremic syndrome (HUS) and neurological symptoms performed during an epidemic outbreak of Escherichia coli O104:H4 in Northern Europe has previously shown pathological changes in only approximately 50% of patients. In contrast, susceptibility-weighted imaging (SWI) revealed a loss of venous contrast in a large number of patients. We hypothesized that this observation may be due to an increase in cerebral blood flow (CBF) and aimed to identify a plausible cause.

Materials and Methods

Baseline 1.5T MRI scans of 36 patients (female, 26; male, 10; mean age, 38.2±19.3 years) were evaluated. Venous contrast was rated on standard SWI minimum intensity projections. A prototype four-dimensional (time resolved) magnetic resonance angiography (4D MRA) assessed cerebral hemodynamics by global time-to-peak (TTP), as a surrogate marker for CBF. Clinical parameters studied were hemoglobin, hematocrit, creatinine, urea levels, blood pressure, heart rate, and end-tidal CO₂.

Results

SWI venous contrast was abnormally low in 33 of 36 patients. TTP ranged from 3.7 to 10.2 frames (mean, 7.9 ± 1.4). Hemoglobin at the time of MRI (n = 35) was decreased in all patients (range, 5.0 to 12.6 g/dL; mean, 8.2 ± 1.4); hematocrit (n = 33) was abnormally low in all but a single patient (range, 14.3 to 37.2%; mean, 23.7 ± 4.2). Creatinine was abnormally high in 30 of 36 patients (83%) (range, 0.8 to 9.7; mean, 3.7 ± 2.2). SWI venous contrast correlated significantly with hemoglobin (r = 0.52, P = 0.0015), hematocrit (r = 0.65, P < 0.001), and TTP (r = 0.35, P = 0.036). No correlation of SWI with blood pressure, heart rate, end-tidal CO₂, creatinine, and urea level was observed. Findings suggest that the loss of venous contrast is related to an increase in CBF secondary to severe anemia related to HUS. SWI contrast of patients with pathological conventional MRI findings was significantly lower compared to patients with normal MRI (mean SWI score, 1.41 and 2.05, respectively; P = 0.04). In patients with abnormal conventional MRI, mean TTP (7.45), mean hemoglobin (7.65), and mean hematocrit (22.0) were lower compared to patients with normal conventional MRI scans (mean TTP = 8.28, mean hemoglobin = 8.63, mean hematocrit = 25.23).

Conclusion

In contrast to conventional MRI, almost all patients showed pathological changes in cerebral hemodynamics assessed by SWI and 4D MRA. Loss of venous contrast on SWI is most likely the result of an increase in CBF and may be related to the acute onset of anemia. Future studies will be needed to assess a possible therapeutic effect of blood transfusions in patients with HUS and neurological symptoms.