Retrospective correction for induced magnetic field inhomogeneity in measurements of large-vessel hemoglobin oxygen saturation by MR susceptometry

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.

Sleep-stage-dependent alterations in cerebral oxygen metabolism quantified by magnetic resonance

A key function of sleep is to provide a regular period of reduced brain metabolism, which is critical for maintenance of healthy brain function. The purpose of this work was to quantify the sleep-stage-dependent changes in brain energetics in terms of cerebral metabolic rate of oxygen (CMRO2 ) as a function of sleep stage using quantitative magnetic resonance imaging (MRI) with concurrent electroencephalography (EEG) during sleep in the scanner. Twenty-two young and older subjects with regular sleep hygiene and Pittsburgh Sleep Quality Index (PSQI) in the normal range were recruited for the study. Cerebral blood flow (CBF) and venous oxygen saturation (SvO2 ) were obtained simultaneously at 3 Tesla field strength and 2.7-s temporal resolution during an 80-min time series using OxFlow, an in-house developed imaging sequence. The method yields whole-brain CMRO2 in absolute physiologic units via Fick's Principle. Nineteen subjects yielded evaluable data free of subject motion artifacts. Among these subjects, 10 achieved slow-wave (N3) sleep, 16 achieved N2 sleep, and 19 achieved N1 sleep while undergoing the MRI protocol during scanning. Mean CMRO2 was 98 ± 7(μmol min-1 )/100 g awake, declining progressively toward deepest sleep stage: 94 ± 10.8 (N1), 91 ± 11.4 (N2), and 76 ± 9.0 μmol min-1 /100 g (N3), with each level differing significantly from the wake state. The technology described is able to quantify cerebral oxygen metabolism in absolute physiologic units along with non-REM sleep stage, indicating brain oxygen consumption to be closely associated with depth of sleep, with deeper sleep stages exhibiting progressively lower CMRO2 levels.

MRI-based quantification of whole-organ renal metabolic rate of oxygen

During the early stages of diabetes, kidney oxygen utilization increases. The mismatch between oxygen demand and supply contributes to tissue hypoxia, a key driver of chronic kidney disease. Thus, whole-organ renal metabolic rate of oxygen (rMRO2 ) is a potentially valuable biomarker of kidney function. The key parameters required to determine rMRO2 include the renal blood flow rate (RBF) in the feeding artery and oxygen saturation in the draining renal vein (SvO2 ). However, there is currently no noninvasive method to quantify rMRO2 in absolute physiologic units. Here, a new MRI pulse sequence, Kidney Metabolism of Oxygen via T2 and Interleaved Velocity Encoding (K-MOTIVE), is described, along with evaluation of its performance in the human kidney in vivo. K-MOTIVE interleaves a phase-contrast module before a background-suppressed T2 -prepared balanced steady-state-free-precession (bSSFP) readout to measure RBF and SvO2 in a single breath-hold period of 22 s, yielding rMRO2 via Fick's principle. Variants of K-MOTIVE to evaluate alternative bSSFP readout strategies were studied. Kidney mass was manually determined from multislice gradient recalled echo images. Healthy subjects were recruited to quantify rMRO2 of the left kidney at 3-T field strength (N = 15). Assessments of repeat reproducibility and comparisons with individual measurements of RBF and SvO2 were performed, and the method's sensitivity was evaluated with a high-protein meal challenge (N = 8). K-MOTIVE yielded the following metabolic parameters: T2  = 157 ± 19 ms; SvO2  = 92% ± 6%; RBF = 400 ± 110 mL/min; and rMRO2  = 114 ± 117(μmol O2 /min)/100 g tissue. Reproducibility studies of T2 and RBF (parameters directly measured by K-MOTIVE) resulted in coefficients of variation less than 10% and intraclass correlation coefficients more than 0.75. The high-protein meal elicited an increase in rMRO2 , which was corroborated by serum biomarkers. The K-MOTIVE sequence measures SvO2 and RBF, the parameters necessary to quantify whole-organ rMRO2 , in a single breath-hold. The present work demonstrates that rMRO2 quantification is feasible with good reproducibility. rMRO2 is a potentially valuable physiological biomarker.

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.

Repeatability assessment for simultaneous measurement of arterial blood flow, venous oxygen saturation, and muscle perfusion following dynamic exercise

The purpose of the present study was to demonstrate a new sequence and determine the repeatability of simultaneous dynamic measurements of blood flow, venous oxygen saturation (SvO₂ ), and relative perfusion (change from resting perfusion) in calf muscle during recovery from plantar flexion exercise. The feasibility of near simultaneous measurement of bio-energetic parameters was also demonstrated. A sequence was developed to simultaneously measure arterial blood flow using flow-encoded projection, SvO₂ using susceptibility-based oximetry, and relative perfusion using arterial spin labeling in combination with dynamic plantar flexion exercise. The parameters were determined at rest and during recovery from single leg plantar flexion exercise. Test-retest repeatability was analyzed using Bland-Altman analysis and intraclass correlation coefficients (ICC). The mitochondrial capacity of skeletal muscle was also measured immediately afterwards with dynamic phosphorus magnetic resonance spectroscopy. Eight healthy subjects participated in the study for test-retest repeatability. Popliteal artery blood flow at rest was 1.79 ± 0.58 ml/s and increased to 11.18 ± 3.02 ml/s immediately after exercise. Popliteal vein SvO₂ decreased to 45.93% ± 6.5% from a resting value of 70.46% ± 4.76% following exercise. Relative perfusion (change from rest value) was 51.83 ± 15.00 ml/100 g/min at the cessation of exercise. The recovery of blood flow and SvO₂ was modeled as a single exponential with time constants of 38.03 ± 6.91 and 71.19 ± 14.53 s, respectively. All the measured parameters exhibited good repeatability with ICC ranging from 0.8 to 0.95. Bioenergetics measurements were within normal range, demonstrating the feasibility of near simultaneous measurement of hemodynamic and energetic parameters. Clinical feasibility was assessed with Barth syndrome patients, demonstrating reduced oxygen extraction from the blood and reduced mitochondrial oxidative capacity compared with healthy controls. The proposed protocol allows rapid imaging of multiple parameters in skeletal muscle that might be affected in disease.

Cerebral hemodynamic and metabolic dysregulation in the postradiation brain

Technological advances in the delivery of radiation and other novel cancer therapies have significantly improved the 5-year survival rates over the last few decades. Although recent developments have helped to better manage the acute effects of radiation, the late effects such as impairment in cognition continue to remain of concern. Accruing data in the literature have implicated derangements in hemodynamic parameters and metabolic activity of the irradiated normal brain as predictive of cognitive impairment. Multiparametric imaging modalities have allowed us to precisely quantify functional and metabolic information, enhancing the anatomic and morphologic data provided by conventional MRI sequences, thereby contributing as noninvasive imaging-based biomarkers of radiation-induced brain injury. In this review, we have elaborated on the mechanisms of radiation-induced brain injury and discussed several novel imaging modalities, including MR spectroscopy, MR perfusion imaging, functional MR, SPECT, and PET that provide pathophysiological and functional insights into the postradiation brain, and its correlation with radiation dose as well as clinical neurocognitive outcomes. Additionally, we explored some innovative imaging modalities, such as quantitative blood oxygenation level-dependent imaging, susceptibility-based oxygenation measurement, and T2-based oxygenation measurement, that hold promise in delineating the potential mechanisms underlying deleterious neurocognitive changes seen in the postradiation setting. We aim that this comprehensive review of a range of imaging modalities will help elucidate the hemodynamic and metabolic injury mechanisms underlying cognitive impairment in the irradiated normal brain in order to optimize treatment regimens and improve the quality of life for these patients.

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.

MRI evaluation of cerebral metabolic rate of oxygen (CMRO2) in obstructive sleep apnea

Patients with obstructive sleep apnea (OSA) are at elevated risk of developing systemic vascular disease and cognitive dysfunction. Here, cerebral oxygen metabolism was assessed in patients with OSA by means of a magnetic resonance-based method involving simultaneous measurements of cerebral blood flow rate and venous oxygen saturation in the superior sagittal sinus for a period of 10 minutes at an effective temporal resolution of 1.3 seconds before, during, and after repeated 24-second breath-holds mimicking spontaneous apneas, yielding, along with pulse oximetry-derived arterial saturation, whole-brain CMRO2 via Fick's Principle. Enrolled subjects were classified based on their apnea-hypopnea indices into OSA (N = 31) and non-sleep apnea reference subjects (NSA = 21), and further compared with young healthy subjects (YH, N = 10). OSA and NSA subjects were matched for age and body mass index. CMRO2 was lower in OSA than in the YH group during normal breathing (105.6 ± 14.1 versus 123.7 ± 22.8 μmol O2/min/100g, P = 0.01). Further, the fractional change in CMRO2 in response to a breath-hold challenge was larger in OSA than in the YH group (15.2 ± 9.2 versus 8.5 ± 3.4%, P = 0.04). However, there was no significant difference in CMRO2 between OSA and NSA subjects. The data suggest altered brain oxygen metabolism in OSA and possibly in NSA as well.

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.

Comparison of MRI methods for measuring whole‐brain oxygen extraction fraction under different geometric conditions at 7T

Background and Purpose

Methods

Results

Conclusion

Evaluation of Vascular Reactivity of Maternal Vascular Adaptations of Pregnancy With Quantitative MRI: Pilot Study

BACKGROUND

Abnormal maternal vascular function during pregnancy stemming from systemic endothelial dysfunction (EDF) has a central role in the pathophysiology of preeclampsia (PE).

PURPOSE

To utilize quantitative MRI to investigate changes in physiological measures of vascular reactivity during normal pregnancy, and to explore EDF associated with preeclampsia.

STUDY TYPE

Prospective.

POPULATION

Healthy pregnant (HP) (n = 14, mean GA = 26 ± 7 weeks) and nonpregnant women (NP; n = 14); newly postpartum (PP <48 hours) women with severe PE (PP-PE; n = 4) and normotensive pregnancy (PP-HP; n = 5).

FIELD STRENGTH/SEQUENCE

1.5T/3T. RF spoiled multiecho gradient-recalled echo, 1D phase-contrast MRI, time-of-flight.

ASSESSMENT

The micro- and macrovascular function (vasodilatory capacity of arterioles and conduit arteries, respectively) of the femoral vascular bed was evaluated with MRI-based venous oximetry, arterial velocimetry, and luminal flow-mediated dilation quantification, during cuff-induced reactive hyperemia. Aortic arch pulse-wave velocity (aPWV) was quantified to assess arterial stiffness using an ungated 1D technique.

STATISTICAL TESTS

Two-tailed unpaired t-tests were performed to address our two, primary a priori comparisons, HP vs. NP, and PP-PE vs. PP-HP. Given the pilot nature of this study, adjustments for multiple comparisons were not performed.

RESULTS

In HP, microvascular function was attenuated compared to NP by a significant increase in the washout time (10 ± 2 vs. 8 ± 2 sec; P < 0.05) and reduced upslope (2.1 ± 0.5 vs. 3.2 ± 0.8%HbO₂ /s; P < 0.05), time of forward flow (28 ± 5 vs. 33 ± 6 sec, P < 0.05), and hyperemic index (11 ± 3 vs. 16 ± 4 cm/s² ; P < 0.05), but luminal flow-mediated dilatation (FMD_L )was comparable between HP and NP. PP-PE exhibited significant vascular dysfunction compared to PP-HP, as evidenced by differences in upslope (2.2 ± 0.6 vs. 1.3 ± 0.2%HbO₂ /s, P < 0.05), overshoot (16 ± 5 vs. 7 ± 3%HbO₂ , P < 0.05), time of forward flow (28 ± 6 vs. 15 ± 7 s, P < 0.05), and aPWV (7 ± 1 vs. 8 ± 1 m/s, P < 0.05).

DATA CONCLUSION

Attenuated vascular reactivity during pregnancy suggests that the systemic vasodilatory state partially depletes nitric oxide bioavailability. Preliminary data support the potential for MRI to identify vascular dysfunction in vivo that underlies PE. Level of Evidence 2 Technical Efficacy Stage 1 J. MAGN. RESON. IMAGING 2021;53:447-455.

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.

Calibrated fMRI for dynamic mapping of CMRO2 responses using MR-based measurements of whole-brain venous oxygen saturation

Functional MRI (fMRI) can identify active foci in response to stimuli through BOLD signal fluctuations, which represent a complex interplay between blood flow and cerebral metabolic rate of oxygen (CMRO₂) changes. Calibrated fMRI can disentangle the underlying contributions, allowing quantification of the CMRO₂ response. Here, whole-brain venous oxygen saturation (Y_v) was computed alongside ASL-measured CBF and BOLD-weighted data to derive the calibration constant, M, using the proposed Y_v-based calibration. Data were collected from 10 subjects at 3T with a three-part interleaved sequence comprising background-suppressed 3D-pCASL, 2D BOLD-weighted, and single-slice dual-echo GRE (to measure Y_v via susceptometry-based oximetry) acquisitions while subjects breathed normocapnic/normoxic, hyperoxic, and hypercapnic gases, and during a motor task. M was computed via Y_v-based calibration from both hypercapnia and hyperoxia stimulus data, and results were compared to conventional hypercapnia or hyperoxia calibration methods. Mean M in gray matter did not significantly differ between calibration methods, ranging from 8.5 ± 2.8% (conventional hyperoxia calibration) to 11.7 ± 4.5% (Y_v-based calibration in response to hyperoxia), with hypercapnia-based M values between (p = 0.56). Relative CMRO₂ changes from finger tapping were computed from each M map. CMRO₂ increased by ∼20% in the motor cortex, and good agreement was observed between the conventional and proposed calibration methods.

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.

Imaging the physiological evolution of the ischemic penumbra in acute ischemic stroke

We review the hemodynamic, metabolic and cellular parameters affected during early ischemia and their changes as a function of approximate cerebral blood flow ( CBF) thresholds. These parameters underlie the current practical definition of an ischemic penumbra, namely metabolically affected but still viable brain tissue. Such tissue is at risk of infarction under continuing conditions of reduced CBF, but can be rescued through timely intervention. This definition will be useful in clinical diagnosis only if imaging techniques exist that can rapidly, and with sufficient accuracy, visualize the existence of a mismatch between such a metabolically affected area and regions that have suffered cell depolarization. Unfortunately, clinical data show that defining the outer boundary of the penumbra based solely on perfusion-related thresholds may not be sufficiently accurate. Also, thresholds for CBF and cerebral blood volume ( CBV) differ for white and gray matter and evolve with time for both inner and outer penumbral boundaries. As such, practical penumbral imaging would involve parameters in which the physiology is immediately displayed in a manner independent of baseline CBF or CBF threshold, namely pH, oxygen extraction fraction ( OEF), diffusion constant and mean transit time ( MTT). Suitable imaging technologies will need to meet this requirement in a 10-20 min exam.

Quantitative Susceptibility Mapping (QSM) Algorithms: Mathematical Rationale and Computational Implementations

Quantitative susceptibility mapping (QSM) solves the magnetic field-to-magnetization (tissue susceptibility) inverse problem under conditions of noisy and incomplete field data acquired using magnetic resonance imaging. Therefore, sophisticated algorithms are necessary to treat the ill-posed nature of the problem and are reviewed here. The forward problem is typically presented as an integral form, where the field is the convolution of the dipole kernel and tissue susceptibility distribution. This integral form can be equivalently written as a partial differential equation (PDE). Algorithmic challenges are to reduce streaking and shadow artifacts characterized by the fundamental solution of the PDE. Bayesian maximum a posteriori estimation can be employed to solve the inverse problem, where morphological and relevant biomedical knowledge (specific to the imaging situation) are used as priors. As the cost functions in Bayesian QSM framework are typically convex, solutions can be robustly computed using a gradient-based optimization algorithm. Moreover, one can not only accelerate Bayesian QSM, but also increase its effectiveness at reducing shadows using prior knowledge based preconditioners. Improving the efficiency of QSM is under active development, and a rigorous analysis of preconditioning needs to be carried out for further investigation.Quantitative susceptibility mapping (QSM) solves the magnetic field-to-magnetization (tissue susceptibility) inverse problem under conditions of noisy and incomplete field data acquired using magnetic resonance imaging. Therefore, sophisticated algorithms are necessary to treat the ill-posed nature of the problem and are reviewed here. The forward problem is typically presented as an integral form, where the field is the convolution of the dipole kernel and tissue susceptibility distribution. This integral form can be equivalently written as a partial differential equation (PDE). Algorithmic challenges are to reduce streaking and shadow artifacts characterized by the fundamental solution of the PDE. Bayesian maximum a posteriori estimation can be employed to solve the inverse problem, where morphological and relevant biomedical knowledge (specific to the imaging situation) are used as priors. As the cost functions in Bayesian QSM framework are typically convex, solutions can be robustly computed using a gradient-based optimization algorithm. Moreover, one can not only accelerate Bayesian QSM, but also increase its effectiveness at reducing shadows using prior knowledge based preconditioners. Improving the efficiency of QSM is under active development, and a rigorous analysis of preconditioning needs to be carried out for further investigation.

In vivo whole-blood T2 versus HbO2 calibration by modulating blood oxygenation level in the femoral vein through intermittent cuff occlusion

PURPOSE

To investigate the feasibility of estimating calibration constants (K and T_2o ) in vivo for converting whole-blood T₂ to blood hemoglobin oxygen saturation (HbO₂ ) according to the Luz-Meiboom model, 1/T2=1/T2o+K(1-HbO2)2, where K and T_2o are relaxivity and transverse relaxation time of fully saturated blood, respectively.

METHODS

A range of HbO₂ values was achieved in the superficial femoral vein with intermittent cuff occlusion in seven healthy adults (four males) to establish a calibration curve between blood T₂ and HbO₂ at 1.5T. HbO₂ was derived via MR susceptometry, a technique previously validated, and the transverse relaxation time was quantified with an optimized T₂ -prepared balanced steady-state free precession pulse sequence. To evaluate the accuracy of the in vivo calibration method, T₂ and HbO₂ were quantified in the superior sagittal sinus in six additional subjects and compared with susceptometry.

RESULTS

Two sets of gender-specific calibration constants were derived, one for each gender corresponding to hematocrits of 0.47 ± 0.02 for males and 0.38 ± 0.01 for females, yielding K/T_2o  = 41 Hz/260 ms and 26 Hz/280 ms, respectively. The in vivo calibration returned physiologically plausible superior sagittal sinus SvO₂ values (65 ± 5% HbO₂ ), and there was no significant difference between the results from the two methods (average difference -0.3% HbO₂ ).

CONCLUSION

The results show feasibility of performing in vivo calibration for converting whole-blood T₂ to HbO₂ . The proposed approach bypasses the involved and cumbersome processes associated with in vitro calibration. Magn Reson Med 79:2290-2296, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

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.

High-speed whole-brain oximetry by golden-angle radial MRI

PURPOSE

To determine whole-brain cerebral metabolic rate of oxygen (CMRO₂ ), an improved imaging approach, based on radial encoding, termed radial OxFlow (rOxFlow), was developed to simultaneously quantify draining vein venous oxygen saturation (SvO₂ ) and total cerebral blood flow (tCBF).

METHODS

To evaluate the efficiency and precision of the rOxFlow sequence, 10 subjects were studied during a paradigm of repeated breath-holds with both rOxFlow and Cartesian OxFlow (cOxFlow) sequences. CMRO₂ was calculated at baseline from OxFlow-measured data assuming an arterial O₂ saturation of 97%, and the SvO₂ and tCBF breath-hold responses were quantified.

RESULTS

Average neurometabolic-vascular parameters across the 10 subjects for cOxFlow and rOxFlow were, respectively: SvO₂ (%) baseline: 64.6 ± 8.0 versus 64.2 ± 6.6; SvO₂ peak: 70.5 ± 8.5 versus 72.6 ± 5.4; tCBF (mL/min/100 g) baseline: 39.2 ± 3.8 versus 40.6 ± 8.0; tCBF peak: 53.2 ± 5.1 versus 56.1 ± 11.7; CMRO₂ (µmol O₂ /min/100 g) baseline: 111.5 ± 26.8 versus 120.1 ± 19.6. The above measures were not significantly different between sequences (P > 0.05).

CONCLUSION

There was good agreement between the two methods in terms of the physiological responses measured. Comparing the two, rOxFlow provided higher temporal resolution and greater flexibility for reconstruction while maintaining high SNR. Magn Reson Med 79:217-223, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

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.

Evaluating the Role of Reduced Oxygen Saturation and Vascular Damage in Traumatic Brain Injury Using Magnetic Resonance Perfusion-Weighted Imaging and Susceptibility-Weighted Imaging and Mapping

The cerebral vasculature, along with neurons and axons, is vulnerable to biomechanical insult during traumatic brain injury (TBI). Trauma-induced vascular injury is still an underinvestigated area in TBI research. Cerebral blood flow and metabolism could be important future treatment targets in neural critical care. Magnetic resonance imaging offers a number of key methods to probe vascular injury and its relationship with traumatic hemorrhage, perfusion deficits, venous blood oxygen saturation changes, and resultant tissue damage. They make it possible to image the hemodynamics of the brain, monitor regional damage, and potentially show changes induced in the brain's function not only acutely but also longitudinally following treatment. These methods have recently been used to show that even mild TBI (mTBI) subjects can have vascular abnormalities, and thus they provide a major step forward in better diagnosing mTBI patients.

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.

Acute hypoxia increases the cerebral metabolic rate - a magnetic resonance imaging study

The aim of the present study was to examine changes in cerebral metabolism by magnetic resonance imaging of healthy subjects during inhalation of 10% O2 hypoxic air. Hypoxic exposure elevates cerebral perfusion, but its effect on energy metabolism has been less investigated. Magnetic resonance imaging techniques were used to measure global cerebral blood flow and the venous oxygen saturation in the sagittal sinus. Global cerebral metabolic rate of oxygen was quantified from cerebral blood flow and arteriovenous oxygen saturation difference. Concentrations of lactate, glutamate, N-acetylaspartate, creatine and phosphocreatine were measured in the visual cortex by magnetic resonance spectroscopy. Twenty-three young healthy males were scanned for 60 min during normoxia, followed by 40 min of breathing hypoxic air. Inhalation of hypoxic air resulted in an increase in cerebral blood flow of 15.5% (p = 0.058), and an increase in cerebral metabolic rate of oxygen of 8.5% (p = 0.035). Cerebral lactate concentration increased by 180.3% ([Formula: see text]), glutamate increased by 4.7% ([Formula: see text]) and creatine and phosphocreatine decreased by 15.2% (p[Formula: see text]). The N-acetylaspartate concentration was unchanged (p = 0.36). In conclusion, acute hypoxia in healthy subjects increased perfusion and metabolic rate, which could represent an increase in neuronal activity. We conclude that marked changes in brain homeostasis occur in the healthy human brain during exposure to acute hypoxia.

Fat-water MRI of a diet-induced obesity mouse model at 15.2T

Quantitative fat-water MRI (FWMRI) methods provide valuable information about the distribution, volume, and composition of adipose tissue (AT). Ultra high field FWMRI of animal models may have the potential to provide insights into the progression of obesity and its comorbidities. Here, we present quantitative FWMRI with all known confounder corrections on a 15.2T preclinical scanner for noninvasive in vivo monitoring of an established diet-induced obesity mouse model. Male C57BL/6J mice were placed on a low-fat (LFD) or a high-fat diet (HFD). Three-dimensional (3-D) multiple gradient echo MRI at 15.2T was performed at baseline, 4, 8, 12, and 16 weeks after diet onset. A 3-D fat-water separation algorithm and additional processing were used to generate proton-density fat fraction (PDFF), local magnetic field offset, and [Formula: see text] maps. We examined these parameters in perirenal AT ROIs from LFD and HFD mice. The data suggest that PDFF, local field offset, and [Formula: see text] have different time course behaviors between LFD and HFD mice over 16 weeks. This work suggests FWMRI at 15.2T may be a useful tool for longitudinal studies of adiposity due to the advantages of ultra high field although further investigation is needed to understand the observed time course behavior.

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

A recently reported quantitative magnetic resonance imaging (MRI) method denoted OxFlow has been shown to be able to quantify whole-brain cerebral metabolic rate of oxygen (CMRO2) by simultaneously measuring oxygen saturation (SvO2) in the superior sagittal sinus and cerebral blood flow (CBF) in the arteries feeding the brain in 30 seconds, which is adequate for measurement at baseline but not necessarily in response to neuronal activation. Here, we present an accelerated version of the method (referred to as F-OxFlow) that quantifies CMRO2 in 8 seconds scan time under full retention of the parent method's capabilities and compared it with its predecessor at baseline in 10 healthy subjects. Results indicate excellent agreement between both sequences, with mean bias of 2.2% (P=0.18, two-tailed t-test), 3.4% (P=0.08, two-tailed t-test), and 2.0% (P=0.56, two-tailed t-test) for SvO2, CBF, and CMRO2, respectively. F-OxFlow's potential to monitor dynamic changes in SvO2, CBF, and CMRO2 is illustrated in a paradigm of volitional apnea applied to five of the study subjects. The sequence captured an average increase in SvO2, CBF, and CMRO2 of 10.1±2.5%, 43.2±9.2%, and 7.1±2.2%, respectively, in good agreement with literature values. The method may therefore be suited for monitoring alterations in CBF and SvO2 in response to neurovascular stimuli.

Correction and optimization of a T2-based approach to map blood oxygenation in small cerebral veins

PURPOSE

Cerebral venous blood oxygenation (Yv ) is an important biomarker in brain physiology and function. The present study proposes a procedure to provide a quantitative map of the brain's intravascular Yv.

THEORY AND METHODS

The method is based on a pulse sequence, T2 -Relaxation-Under-Phase-Contrast (TRU-PC) MRI, with postprocessing approaches to correct eddy-current effects. A complete scan protocol consists of four TRU-PC scans sensitized to large and small vessels with anterior-posterior and foot-head flow-encoding directions, and the data are analyzed conjunctively. Eddy-current correction was performed by fitting the tissue phase to a hyperplane, and then subtracting the eddy-current phase from the measured vessel phase. The reproducibility of the Yv-maps was examined in five participants. Sensitivity of the Yv map to a caffeine challenge was studied in another five participants.

RESULTS

Removal of eddy-current induced artifact allowed for the correction of T2 measurements, as demonstrated in vivo and with simulation. A Yv-map depicting all vessels in the slice can be obtained with the proposed protocol. Test-retest variability of the Yv -map was 3.7 ± 1.2%. Yv reduction can be reliably detected (P < 0.001) following the caffeine ingestion.

CONCLUSION

With the proposed TRU-PC protocol and eddy-current correction procedure, an accurate, vessel-specific Yv map of the human brain can be obtained.

The effect of dissolved oxygen on the susceptibility of blood

PURPOSE

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

THEORY AND METHODS

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

RESULTS

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

CONCLUSION

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

Effects of age and smoking on endothelial function assessed by quantitative cardiovascular magnetic resonance in the peripheral and central vasculature

BACKGROUND

Both age and smoking promote endothelial dysfunction and impair vascular reactivity. Here, we tested this hypothesis by quantifying new cardiovascular magnetic resonance (CMR)-based biomarkers in smokers and nonsmokers.

METHODS

Study population: young non-smokers (YNS: N = 45, mean age = 30.2 ± 0.7 years), young smokers (YS: N = 39 mean age 32.1 ± 0.7 years), older non-smokers (ONS: N = 45, mean age = 57.8 ± 0.6 years), and older smokers (OS: N = 40, mean age = 56.3 ± 0.6 years), all without overt cardiovascular disease. Vascular reactivity was evaluated following cuff-induced hyperemia via time-resolved blood flow velocity and oxygenation (SvO2) in the femoral artery and vein, respectively. SvO2 dynamics yielded washout time (time to minimum SvO2), resaturation rate (upslope) and maximum change from baseline (overshoot). Arterial parameters included pulse ratio (PR), hyperemic index (HI) and duration of hyperemia (TFF). Pulse-wave velocity (PWV) was assessed in aortic arch, thoracoabdominal aorta and iliofemoral arteries. Ultrasound-based carotid intimal-medial thickness (IMT) and brachial flow-mediated dilation were measured for comparison.

RESULTS

Age and smoking status were independent for all parameters. Smokers had reduced upslope (-28.4%, P < 0.001), increased washout time (+15.3%, P < 0.01), and reduced HI (-19.5%, P < 0.01). Among non-smokers, older subjects had lower upslope (-22.7%, P < 0.01) and overshoot (-29.4%, P < 0.01), elevated baseline pulse ratio (+14.9%, P < 0.01), central and peripheral PWV (all P < 0.05). Relative to YNS, YS had lower upslope (-23.6%, P < 0.01) and longer washout time (13.5%, P < 0.05). Relative to ONS, OS had lower upslope (-33.0%, P < 0.01). IMT was greater in ONS than in YNS (+45.6%, P < 0.001), and also in YS compared to YNS (+14.7%, P < 0.05).

CONCLUSIONS

Results suggest CMR biomarkers of endothelial function to be sensitive to age and smoking independent of each other.

Cerebral hemodynamic changes of mild traumatic brain injury at the acute stage

Mild traumatic brain injury (mTBI) is a significant public health care burden in the United States. However, we lack a detailed understanding of the pathophysiology following mTBI and its relation to symptoms and recovery. With advanced magnetic resonance imaging (MRI), we can investigate brain perfusion and oxygenation in regions known to be implicated in symptoms, including cortical gray matter and subcortical structures. In this study, we assessed 14 mTBI patients and 18 controls with susceptibility weighted imaging and mapping (SWIM) for blood oxygenation quantification. In addition to SWIM, 7 patients and 12 controls had cerebral perfusion measured with arterial spin labeling (ASL). We found increases in regional cerebral blood flow (CBF) in the left striatum, and in frontal and occipital lobes in patients as compared to controls (p = 0.01, 0.03, 0.03 respectively). We also found decreases in venous susceptibility, indicating increases in venous oxygenation, in the left thalamostriate vein and right basal vein of Rosenthal (p = 0.04 in both). mTBI patients had significantly lower delayed recall scores on the standardized assessment of concussion, but neither susceptibility nor CBF measures were found to correlate with symptoms as assessed by neuropsychological testing. The increased CBF combined with increased venous oxygenation suggests an increase in cerebral blood flow that exceeds the oxygen demand of the tissue, in contrast to the regional hypoxia seen in more severe TBI. This may represent a neuroprotective response following mTBI, which warrants further investigation.

Quantitative susceptibility mapping (QSM): Decoding MRI data for a tissue magnetic biomarker

In MRI, the main magnetic field polarizes the electron cloud of a molecule, generating a chemical shift for observer protons within the molecule and a magnetic susceptibility inhomogeneity field for observer protons outside the molecule. The number of water protons surrounding a molecule for detecting its magnetic susceptibility is vastly greater than the number of protons within the molecule for detecting its chemical shift. However, the study of tissue magnetic susceptibility has been hindered by poor molecular specificities of hitherto used methods based on MRI signal phase and T2* contrast, which depend convolutedly on surrounding susceptibility sources. Deconvolution of the MRI signal phase can determine tissue susceptibility but is challenged by the lack of MRI signal in the background and by the zeroes in the dipole kernel. Recently, physically meaningful regularizations, including the Bayesian approach, have been developed to enable accurate quantitative susceptibility mapping (QSM) for studying iron distribution, metabolic oxygen consumption, blood degradation, calcification, demyelination, and other pathophysiological susceptibility changes, as well as contrast agent biodistribution in MRI. This paper attempts to summarize the basic physical concepts and essential algorithmic steps in QSM, to describe clinical and technical issues under active development, and to provide references, codes, and testing data for readers interested in QSM.

Feasibility and reproducibility of measurement of whole muscle blood flow, oxygen extraction, and VO2 with dynamic exercise using MRI

PURPOSE

Develop an MRI method to estimate skeletal muscle oxygen consumption (VO2 ) with dynamic exercise using simultaneous measurement of venous blood flow (VBF) and venous oxygen saturation (SvO2 ).

METHODS

Real-time imaging of femoral VBF using a complex-difference method was interleaved with imaging of venous hemoglobin oxygen saturation (SvO2 ) using magnetic susceptometry to estimate muscle VO2 (Fick principle). Nine healthy subjects performed repeated 5-watt knee-extension (quadriceps) exercise within the bore of a 1.5 Tesla MRI scanner, for test/re-test comparison. VBF, SvO2 , and derived VO2 were estimated at baseline and immediately (<1 s) postexercise and every 2.4 s for 4 min.

RESULTS

Quadriceps muscle mass was 2.43 ± 0.31 kg. Mean baseline values were VBF = 0.13 ± 0.06 L/min/kg, SvO2 = 69.4 ± 10.1%, and VO2 = 6.8 ± 4.1 mL/min/kg. VBF, SvO2 , and VO2 values from peak exercise had good agreement between trials (VBF = 0.9 ± 0.1 versus 1.0 ± 0.1 L/min/kg, R(2) = 0.83, CV = 7.6%; SvO2 = 43.2 ± 13.5 versus 40.9 ± 13.1%, R(2) = 0.88, CV = 15.6%; VO2 = 95.7 ± 18.0 versus 108.9 ± 17.3 mL/min/kg, R(2) = 0.88, CV = 12.3%), as did the VO2 recovery time constant (26.1 ± 3.5 versus 26.0 ± 4.0 s, R(2) = 0.85, CV = 6.0%). CV = coefficient of variation.

CONCLUSION

Rapid imaging of VBF and SvO2 for the estimation of whole muscle VO2 is compatible with dynamic exercise for the estimation of peak values and recovery dynamics following exercise with good reproducibility.

Rapid T2- and susceptometry-based CMRO2 quantification with interleaved TRUST (iTRUST)

Susceptometry-based oximetry (SBO) and T2-relaxation-under-spin-tagging (TRUST) are two promising methods for quantifying the cerebral metabolic rate of oxygen (CMRO2), a critical parameter of brain function. We present a combined method, interleaved TRUST (iTRUST), which achieves rapid, simultaneous quantification of both susceptometry- and T2-based CMRO2 via insertion of a flow-encoded, dual-echo gradient-recalled echo (OxFlow) module within the T1 recovery portion of the TRUST sequence. In addition to allowing direct comparison between SBO- and TRUST-derived venous oxygen saturation (Yv) values, iTRUST substantially improves TRUST temporal resolution for CMRO2 quantification and obviates the need for a separate blood flow measurement following TRUST acquisition. iTRUST was compared directly to TRUST and OxFlow alone in three resting subjects at baseline, exhibiting close agreement with the separate techniques and comparable precision. These baseline data as well as simulation results support the use of two instead of the traditional four T2 preparation times for T2 fitting, allowing simultaneous quantification of susceptometry- and T2-based Yv (and CMRO2) with three- and six-second temporal resolution, respectively. In 10 young healthy subjects, iTRUST was applied during a 5% CO2 gas mixture-breathing paradigm. T2-based Yv values were lower at baseline relative to susceptometry (62.3 ± 3.1 vs. 66.7 ± 5.1 %HbO2, P<0.05), but increased more in response to hypercapnia. As a result, T2-based CMRO2 decreased from 140.4 ± 9.7 to 120.0 ± 9.5 μMol/100g/min, a significant -14.6 ± 3.6% response (P < 0.0001), whereas susceptometry-based CMRO2 changed insignificantly from 123.4 ± 18.7 to 127.9 ± 25.7, a 3.3 ± 9.7% response (P = 0.31). These differing results are in accord with previous studies applying the parent OxFlow or TRUST sequences individually, thus supporting the reliability of iTRUST but also strongly suggesting that a systematic bias exists between the susceptometry- and T2-based Yv quantification techniques.

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

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

Comparison of MRI methods for measuring whole-brain venous oxygen saturation

PURPOSE

In this work, we compare susceptometry-based oximetry (SBO) and two T2 -based methods for estimating resting baseline SvO2 in the superior sagittal sinus (SSS).

METHODS

SBO is a field-mapping technique whereas in T2 -based methods the intravascular blood signal is isolated either with velocity-encoded projections [projection-based T2 (PT2 )] or a tag-control scheme [T2 -relaxation under spin tagging (TRUST)] after T2 -preparation. The measurements were performed on twelve healthy subjects (mean age = 33 ± 6 years) at 3 Tesla field strength. The reliability, precision, and reproducibility were examined for the three techniques.

RESULTS

The mean (± standard deviation) SvO2 quantified by SBO, PT2 , and TRUST were found to be 65.9 ± 3.3, 65.6 ± 3.5, and 63.2 ± 4.1%. The standard deviation (SD) for 10 consecutive measurements in the quantified SvO2 was less than 2.7%, 4.7%, and 5.0% for SBO, PT2 , and TRUST across all subjects. In testing reproducibility across different days, the resulting SDs were 2.6, 3.5, and 2.0% for SBO, PT2 , and TRUST.

CONCLUSION

The results indicate that all three SvO2 quantification techniques to be reliable with good agreement between PT2 and SBO while TRUST yielded slightly lower values compared with the other two techniques.

Measuring venous blood oxygenation in fetal brain using susceptibility-weighted imaging

PURPOSE

To evaluate fetal cerebral venous blood oxygenation, Yv, using principles of MR susceptometry.

MATERIALS AND METHODS

A cohort of 19 pregnant subjects, with a mean gestational age of 31.6 ± 4.7 weeks were imaged using a modified susceptibility-weighted imaging (SWI) sequence. Data quality was first assessed for feasibility of oxygen saturation measurement, and data from five subjects (mean ± std gestational age of 33.7 ± 3.6 weeks) were then chosen for further quantitative analysis. SWI phase in the superior sagittal sinus was used to evaluate oxygen saturation using the principles of MR susceptometry. Systematic error in the measured Y(v) values was studied through simulations.

RESULTS

Simulations showed that the systematic error in Yv depended upon the assumed angle of the vessel, θ, relative to the main magnetic field and the error in that vessel angle δθ. For the typical vessel angle of θ = 30° encountered in the fetal data analyzed, a δθ as large as ±20° led to an absolute error, δYv, of less than 11%. The measured mean oxygen saturation across the five fetuses was 66% ± 9.4%. This average cerebral venous blood oxygenation value is in close agreement with values in the published literature.

CONCLUSION

We have reported the first in vivo measurement of human fetal cerebral venous oxygen saturation using MRI.

Background field removal by solving the Laplacian boundary value problem

The removal of the background magnetic field is a critical step in generating phase images and quantitative susceptibility maps, which have recently been receiving increasing attention. Although it is known that the background field satisfies Laplace's equation, the boundary values of the background field for the region of interest have not been explicitly addressed in the existing methods, and they are not directly available from MRI measurements. In this paper, we assume simple boundary conditions and remove the background field by explicitly solving the boundary value problems of Laplace's or Poisson's equation. The proposed Laplacian boundary value (LBV) method for background field removal retains data near the boundary and is computationally efficient. Tests on a numerical phantom and an experimental phantom showed that LBV was more accurate than two existing methods.

Susceptibility-based analysis of dynamic gadolinium bolus perfusion MRI

PURPOSE

An algorithm is developed for the reconstruction of dynamic, gadolinium (Gd) bolus MR perfusion images of the human brain, based on quantitative susceptibility mapping (QSM).

METHODS

The method is evaluated in five perfusion scans obtained from four different patients scanned at 3 Tesla, and compared with the conventional analysis based on changes in the transverse relaxation rate ΔR2 * and to theoretical predictions. QSM images were referenced to ventricular cerebrospinal fluid (CSF) for each dynamic of the perfusion sequence.

RESULTS

Images of cerebral blood flow and blood volume were successfully reconstructed from the QSM-analysis, and were comparable to those reconstructed using ΔR2 *. The magnitudes of the Gd-associated susceptibility effects in gray and white matter were consistent with theoretical predictions.

CONCLUSION

QSM-based analysis may have some theoretical advantages compared with ΔR2 *, including a simpler relationship between signal change and Gd concentration. However, disadvantages are its much lower contrast-to-noise ratio, artifacts due to respiration and other effects, and more complicated reconstruction methods. More work is required to optimize data acquisition protocols for QSM-based perfusion imaging.

Time-resolved MRI oximetry for quantifying CMRO(2) and vascular reactivity

This brief review of magnetic resonance susceptometry summarizes the methods conceived in the authors' laboratory during the past several years. This article shows how venous oxygen saturation is quantified in large draining veins by field mapping and how this information, in concert with simultaneous measurement of cerebral blood flow, yields cerebral metabolic rate of oxygen, the brain's rate of oxygen consumption. The accuracy of this model-based approach in which the blood vessel is approximated as a long, straight cylinder, for which an analytical solution for the induced field exists, is discussed. It is shown that the approach is remarkably robust, allowing for time-resolved quantification of whole-brain metabolism at rest and in response to stimuli, thereby providing detailed information on cerebral physiology in health and disease not previously amenable by noninvasive methods.

High temporal resolution MRI quantification of global cerebral metabolic rate of oxygen consumption in response to apneic challenge

We present a technique for quantifying global cerebral metabolic rate of oxygen consumption (CMRO2) in absolute physiologic units at 3-second temporal resolution and apply the technique to quantify the dynamic CMRO2 response to volitional apnea. Temporal resolution of 3 seconds was achieved via a combination of view sharing and superior sagittal sinus-based estimation of total cerebral blood flow (tCBF) rather than tCBF measurement in the neck arteries. These modifications were first validated in three healthy adults and demonstrated to produce minimal errors in image-derived blood flow and venous oxygen saturation (SvO2) values. The technique was then applied in 10 healthy adults during an apnea paradigm of three repeated 30-second breath-holds. Subject-averaged baseline tCBF, arteriovenous oxygen difference (AVO2D), and CMRO2 were 48.6 ± 7.0 mL/100 g per minute, 29.4 ± 3.4 %HbO2, and 125.1 ± 11.4 μmol/100 g per minute, respectively. Subject-averaged maximum changes in tCBF and AVO2D were 43.5 ± 9.4% and -32.1 ± 5.7%, respectively, resulting in a small (6.0 ± 3.5%) but statistically significant (P=0.00044, two-tailed t-test) increase in average end-apneic CMRO2. This method could be used to investigate neurometabolic-hemodynamic relationships in normal physiology, to better define the biophysical origins of the BOLD signal, and to quantify neurometabolic responsiveness in diseases of altered neurovascular reactivity.

Quantitative oxygenation venography from MRI phase

PURPOSE

To demonstrate acquisition and processing methods for quantitative oxygenation venograms that map in vivo oxygen saturation (SvO2 ) along cerebral venous vasculature.

METHODS

Regularized quantitative susceptibility mapping (QSM) is used to reconstruct susceptibility values and estimate SvO2 in veins. QSM with ℓ1 and ℓ2 regularization are compared in numerical simulations of vessel structures with known magnetic susceptibility. Dual-echo, flow-compensated phase images are collected in three healthy volunteers to create QSM images. Bright veins in the susceptibility maps are vectorized and used to form a three-dimensional vascular mesh, or venogram, along which to display SvO2 values from QSM.

RESULTS

Quantitative oxygenation venograms that map SvO2 along brain vessels of arbitrary orientation and geometry are shown in vivo. SvO2 values in major cerebral veins lie within the normal physiological range reported by (15) O positron emission tomography. SvO2 from QSM is consistent with previous MR susceptometry methods for vessel segments oriented parallel to the main magnetic field. In vessel simulations, ℓ1 regularization results in less than 10% SvO2 absolute error across all vessel tilt orientations and provides more accurate SvO2 estimation than ℓ2 regularization.

CONCLUSION

The proposed analysis of susceptibility images enables reliable mapping of quantitative SvO2 along venograms and may facilitate clinical use of venous oxygenation imaging.

Combined measurement of perfusion, venous oxygen saturation, and skeletal muscle T2* during reactive hyperemia in the leg

BACKGROUND

The function of the peripheral microvascular may be interrogated by measuring perfusion, tissue oxygen concentration, or venous oxygen saturation (SvO2) recovery dynamics following induced ischemia. The purpose of this work is to develop and evaluate a magnetic resonance (MR) technique for simultaneous measurement of perfusion, SvO2, and skeletal muscle T2*.

METHODS

Perfusion, Intravascular Venous Oxygen saturation, and T2* (PIVOT) is comprised of interleaved pulsed arterial spin labeling (PASL) and multi-echo gradient-recalled echo (GRE) sequences. During the PASL post-labeling delay, images are acquired with a multi-echo GRE to quantify SvO2 and T2* at a downstream slice location. Thus time-courses of perfusion, SvO2, and T2* are quantified simultaneously within a single scan. The new sequence was compared to separately measured PASL or multi-echo GRE data during reactive hyperemia in five young healthy subjects. To explore the impairment present in peripheral artery disease patients, five patients were evaluated with PIVOT.

RESULTS

Comparison of PIVOT-derived data to the standard techniques shows that there was no significant bias in any of the time-course-derived metrics. Preliminary data show that PAD patients exhibited alterations in perfusion, SvO2, and T2* time-courses compared to young healthy subjects.

CONCLUSION

Simultaneous quantification of perfusion, SvO2, and T2* is possible with PIVOT. Kinetics of perfusion, SvO2, and T2* during reactive hyperemia may help to provide insight into the function of the peripheral microvasculature in patients with PAD.

Blood oxygenation level-dependent (BOLD)-based techniques for the quantification of brain hemodynamic and metabolic properties - theoretical models and experimental approaches

The quantitative evaluation of brain hemodynamics and metabolism, particularly the relationship between brain function and oxygen utilization, is important for the understanding of normal human brain operation, as well as the pathophysiology of neurological disorders. It can also be of great importance for the evaluation of hypoxia within tumors of the brain and other organs. A fundamental discovery by Ogawa and coworkers of the blood oxygenation level-dependent (BOLD) contrast opened up the possibility to use this effect to study brain hemodynamic and metabolic properties by means of MRI measurements. Such measurements require the development of theoretical models connecting the MRI signal to brain structure and function, and the design of experimental techniques allowing MR measurements to be made of the salient features of theoretical models. In this review, we discuss several such theoretical models and experimental methods for the quantification of brain hemodynamic and metabolic properties. The review's main focus is on methods for the evaluation of the oxygen extraction fraction (OEF) based on the measurement of the blood oxygenation level. A combination of the measurement of OEF and the cerebral blood flow (CBF) allows an evaluation to be made of the cerebral metabolic rate of oxygen consumption (CMRO2 ). We first consider in detail the magnetic properties of blood - magnetic susceptibility, MR relaxation and theoretical models of the intravascular contribution to the MR signal under different experimental conditions. We then describe a 'through-space' effect - the influence of inhomogeneous magnetic fields, created in the extravascular space by intravascular deoxygenated blood, on the formation of the MR signal. Further, we describe several experimental techniques taking advantage of these theoretical models. Some of these techniques - MR susceptometry and T2 -based quantification of OEF - utilize the intravascular MR signal. Another technique - quantitative BOLD - evaluates OEF by making use of through-space effects. In this review, we target both scientists just entering the MR field and more experienced MR researchers interested in the application of advanced BOLD-based techniques to the study of the brain in health and disease.

Susceptibility weighted imaging of cartilage canals in porcine epiphyseal growth cartilage ex vivo and in vivo

PURPOSE

High-resolution visualization of cartilage canals has been restricted to histological methods and contrast-enhanced imaging. In this study, the feasibility of non-contrast-enhanced susceptibility weighted imaging (SWI) for visualization of the cartilage canals was investigated ex vivo at 9.4 T, further explored at 7 and 3 T and demonstrated in vivo at 7 T, using a porcine animal model.

METHODS

SWI scans of specimens of distal femur and humerus from 1 to 8 week-old piglets were conducted at 9.4 T using 3D-GRE sequence and SWI post-processing. The stifle joints of a 2-week old piglet were scanned ex vivo at 7 and 3 T. Finally, the same sites of a 3-week-old piglet were scanned, in vivo, at 7 T under general anesthesia using the vendor-provided sequences.

RESULTS

High-contrast visualization of the cartilage canals was obtained ex vivo, especially at higher field strengths; the results were confirmed histologically. In vivo feasibility was demonstrated at 7 T and comparison of ex vivo scans at 3 and 7 T indicated feasibility of using SWI at 3 T.

CONCLUSIONS

High-resolution 3D visualization of cartilage canals was demonstrated using SWI. This demonstration of fully noninvasive visualization opens new avenues to explore skeletal maturation and the role of vascular supply for diseases such as osteochondrosis.

Quantitative CMR markers of impaired vascular reactivity associated with age and peripheral artery disease

BACKGROUND

The aim of this study was to develop and evaluate an integrated CMR method incorporating dynamic oximetry, blood flow and pulse-wave velocimetry to assess vascular reactivity in patients with peripheral artery disease (PAD) and healthy controls.

METHODS AND RESULTS

The study population consisted of young healthy subjects (YH, 30 ± 7 yrs, N = 19),PAD (71 ± 9 yrs, N = 38), and older healthy controls (OHC, 68 ± 9 yrs, N = 43). Peripheral vascular reactivity was evaluated with two methods, time-resolved quantification of blood flow velocity and oxygenation level in the femoral artery and vein, respectively, performed simultaneously both at rest and hyperemia. Aortic stiffness was assessed via pulse-wave velocity. Oximetric data showed that compared to OHC, the time-course of the hemoglobin oxygen saturation in PAD patients had longer washout time (28.6 ± 1.2 vs 16.9 ± 1.1 s, p < 0.0001), reduced upslope (0.60 ± 0.1 vs 1.3 ± 0.08 HbO2/sec, p < 0.0001) and lower overshoot (8 ± 1.4 vs 14 ± 1.2 HbO2, p = 0.0064). PAD patients also had longer-lasting antegrade femoral artery flow during hyperemia (51 ± 2.1 vs 24 ± 1.8 s, p < 0.0001), and reduced peak-to-baseline flow rate (3.1 ± 0.5 vs 7.4 ± 0.4, p < 0.0001). Further, the pulsatility at rest was reduced (0.75 ± 0.32 vs 5.2 ± 0.3, p < 0.0001), and aortic PWV was elevated (10.2 ± 0.4 vs 8.1 ± 0.4 m/s, p = 0.0048).

CONCLUSION

The proposed CMR protocol quantifies multiple aspects of vascular reactivity and represents an initial step toward development of a potential tool for evaluating interventions, extrapolating clinical outcomes and predicting functional endpoints based on quantitative parameters.

Imaging brain oxygenation with MRI using blood oxygenation approaches: methods, validation, and clinical applications

SUMMARY

In many pathophysiologic situations, including brain neoplasms, neurodegenerative disease, and chronic and acute ischemia, an imbalance exists between oxygen tissue consumption and delivery. Furthermore, oxygenation changes following a stress challenge, such as with carbogen gas or acetazolamide, can yield information about cerebrovascular reactivity. The unique sensitivity of the BOLD effect to the presence of deoxyhemoglobin has led to its widespread use in the field of cognitive neurosciences. However, the high spatial and temporal resolution afforded by BOLD imaging does not need to be limited to the study of healthy brains. While the complex relationship between the MR imaging signal and tissue oxygenation hinders a direct approach, many different methods have been developed during the past decade to obtain specific oxygenation measurements. These include qBOLD, phase- and susceptibility-based imaging, and intravascular T2-based approaches. The aim of this review is to give an overview of the theoretic basis of these methods as well as their application to measure oxygenation in both healthy subjects and those with disease.

Investigating the magnetic susceptibility properties of fresh human blood for noninvasive oxygen saturation quantification

Quantification of blood oxygen saturation on the basis of a measurement of its magnetic susceptibility demands knowledge of the difference in volume susceptibility between fully oxygenated and fully deoxygenated blood (Δχ(do) ). However, two very different values of Δχ(do) are currently in use. In this work we measured Δχ(do) as well as the susceptibility of oxygenated blood relative to water, Δχ(oxy) , by MR susceptometry in samples of freshly drawn human blood oxygenated to various levels, from 6 to 98% as determined by blood gas analysis. Regression analysis yielded 0.273 ± 0.006 and -0.008 ± 0.003 ppm (cgs) respectively, for Δχ(do) and Δχ(oxy) , in excellent agreement with previous work by Spees et al. (Magn Reson Med 2001;45:533-542).

A novel background field removal method for MRI using projection onto dipole fields (PDF)

For optimal image quality in susceptibility-weighted imaging and accurate quantification of susceptibility, it is necessary to isolate the local field generated by local magnetic sources (such as iron) from the background field that arises from imperfect shimming and variations in magnetic susceptibility of surrounding tissues (including air). Previous background removal techniques have limited effectiveness depending on the accuracy of model assumptions or information input. In this article, we report an observation that the magnetic field for a dipole outside a given region of interest (ROI) is approximately orthogonal to the magnetic field of a dipole inside the ROI. Accordingly, we propose a nonparametric background field removal technique based on projection onto dipole fields (PDF). In this PDF technique, the background field inside an ROI is decomposed into a field originating from dipoles outside the ROI using the projection theorem in Hilbert space. This novel PDF background removal technique was validated on a numerical simulation and a phantom experiment and was applied in human brain imaging, demonstrating substantial improvement in background field removal compared with the commonly used high-pass filtering method.