The effect of dissolved oxygen on the relaxation rates of blood plasma: Implications for hyperoxia calibrated BOLD

MR-oximetry with fat DESPOT

PURPOSE

The R₁ relaxation rate of fat is a promising marker of tissue oxygenation. Existing techniques to map fat R₁ in MR-oximetry offer limited spatial coverage, require long scan times, or pulse sequences that are not readily available on clinical scanners. This work addresses these limitations with a 3D voxel-wise fat R₁ mapping technique for MR-oximetry based on a variable flip angle (VFA) approach at 3 T.

METHODS

Varying levels of dissolved oxygen (O₂) were generated in a phantom consisting of vials of safflower oil emulsion, used to approximate human fat. Joint voxel-wise mapping of fat and water R₁ was performed with a two-compartment VFA model fitted to multi-echo gradient-echo magnitude data acquired at four flip angles, referred to as Fat DESPOT. Global R₁ was also calculated. Variations of fat, water, and global R₁ were investigated as a function of the partial pressure of O₂ (pO₂). Inversion-prepared stimulated echo magnetic resonance spectroscopy was used as the reference technique for R₁ measurements.

RESULTS

Fat R₁ from Fat DESPOT was more sensitive than water R₁ and global R₁ to variations in pO₂, consistent with previous studies performed with different R₁ mapping techniques. Fat R₁ sensitivity to pO₂ variations with Fat DESPOT (median O₂ relaxivity r_{1, O2} = 1.57× 10^-3 s^-1 mmHg^-1) was comparable to spectroscopy-based measurements for methylene, the main fat resonance (median r_{1, O2}= 1.80 × 10^-3 s^-1 mmHg^-1).

CONCLUSION

Fat and water R₁ can be measured on a voxel-wise basis using a two-component fit to multi-echo 3D VFA magnitude data in a clinically acceptable scan time. Fat and water R₁ measured with Fat DESPOT were sensitive to variations in pO₂. These observations suggest an approach to 3D in vivo MR oximetry.

The many layers of BOLD. The effect of hypercapnic and hyperoxic stimuli on macro- and micro-vascular compartments quantified by CVR, M, and CBV across cortical depth

Ultra-high field functional magnetic resonance imaging (fMRI) offers the spatial resolution to measure neuronal activity at the scale of cortical layers. However, cortical depth dependent vascularization differences, such as a higher prevalence of macro-vascular compartments near the pial surface, have a confounding effect on depth-resolved blood-oxygen-level dependent (BOLD) fMRI signals. In the current study, we use hypercapnic and hyperoxic breathing conditions to quantify the influence of all venous vascular and micro-vascular compartments on laminar BOLD fMRI, as measured with gradient-echo (GE) and spin-echo (SE) scan sequences, respectively. We find that all venous vascular and micro-vascular compartments are capable of comparable theoretical maximum signal intensities, as represented by the M-value parameter. However, the capacity for vessel dilation, as reflected by the cerebrovascular reactivity (CVR), is approximately two and a half times larger for all venous vascular compartments combined compared to the micro-vasculature at superficial layers. Finally, there is roughly a 35% difference in estimates of CBV changes between all venous vascular and micro-vascular compartments, although this relative difference was approximately uniform across cortical depth. Thus, our results suggest that fMRI BOLD signal differences across cortical depth are likely caused by differences in dilation properties between macro- and micro-vascular compartments.

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

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

Hemodynamic and metabolic changes during hypercapnia with normoxia and hyperoxia using pCASL and TRUST MRI in healthy adults

Blood oxygenation level-dependent (BOLD) or arterial spin labeling (ASL) MRI with hypercapnic stimuli allow for measuring cerebrovascular reactivity (CVR). Hypercapnic stimuli are also employed in calibrated BOLD functional MRI for quantifying neuronally-evoked changes in cerebral oxygen metabolism (CMRO₂). It is often assumed that hypercapnic stimuli (with or without hyperoxia) are iso-metabolic; increasing arterial CO₂ or O₂ does not affect CMRO₂. We evaluated the null hypothesis that two common hypercapnic stimuli, 'CO₂ in air' and carbogen, are iso-metabolic. TRUST and ASL MRI were used to measure the cerebral venous oxygenation and cerebral blood flow (CBF), from which the oxygen extraction fraction (OEF) and CMRO₂ were calculated for room-air, 'CO₂ in air' and carbogen. As expected, CBF significantly increased (9.9% ± 9.3% and 12.1% ± 8.8% for 'CO₂ in air' and carbogen, respectively). CMRO₂ decreased for 'CO₂ in air' (-13.4% ± 13.0%, p < 0.01) compared to room-air, while the CMRO₂ during carbogen did not significantly change. Our findings indicate that 'CO₂ in air' is not iso-metabolic, while carbogen appears to elicit a mixed effect; the CMRO₂ reduction during hypercapnia is mitigated when including hyperoxia. These findings can be important for interpreting measurements using hypercapnic or hypercapnic-hyperoxic (carbogen) stimuli.

Feasibility of glioblastoma tissue response mapping with physiologic BOLD imaging using precise oxygen and carbon dioxide challenge

OBJECTIVES

Innovative physiologic MRI development focuses on depiction of heterogenous vascular and metabolic features in glioblastoma. For this feasibility study, we employed blood oxygenation level-dependent (BOLD) MRI with standardized and precise carbon dioxide (CO₂) and oxygen (O₂) modulation to investigate specific tumor tissue response patterns in patients with newly diagnosed glioblastoma.

MATERIALS AND METHODS

Seven newly diagnosed untreated patients with suspected glioblastoma were prospectively included to undergo a BOLD study with combined CO₂ and O₂ standardized protocol. %BOLD signal change/mmHg during hypercapnic, hypoxic, and hyperoxic stimulus was calculated in the whole brain, tumor lesion and segmented volumes of interest (VOI) [contrast-enhancing (CE) - tumor, necrosis and edema] to analyze their tissue response patterns.

RESULTS

Quantification of BOLD signal change after gas challenges can be used to identify specific responses to standardized stimuli in glioblastoma patients. Integration of this approach with automatic VOI segmentation grants improved characterization of tumor subzones and edema. Magnitude of BOLD signal change during the 3 stimuli can be visualized at voxel precision through color-coded maps overlayed onto whole brain and identified VOIs.

CONCLUSIONS

Our preliminary investigation shows good feasibility of BOLD with standardized and precise CO₂ and O₂ modulation as an emerging physiologic imaging technique to detail specific glioblastoma characteristics. The unique tissue response patterns generated can be further investigated to better detail glioblastoma lesions and gauge treatment response.

Quantitative T1 and T2 mapping by magnetic resonance fingerprinting (MRF) of the placenta before and after maternal hyperoxia

INTRODUCTION

MR relaxometry has been used to assess placental exchange function, but methods to date are not sufficiently fast to be robust to placental motion. Magnetic resonance fingerprinting (MRF) permits rapid, voxel-wise, intrinsically co-registered T₁ and T₂ mapping. After characterizing measurement error, we scanned pregnant women during air and oxygen breathing to demonstrate MRF's ability to detect placental oxygenation changes.

METHODS

The accuracy of FISP-based, sliding-window reconstructed MRF was tested on phantoms. MRF scans in 9-s breath holds were acquired at 3T in 31 pregnant women during air and oxygen breathing. A mixed effects model was used to test for changes in placenta relaxation times between physiological states, to assess the dependency on gestational age (GA), and the impact of placental motion.

RESULTS

MRF estimates of known phantom relaxation times resulted in mean absolute errors for T₁ of 92 ms (4.8%), but T₂ was less accurate at 16 ms (13.6%). During normoxia, placental T₁ = 1825 ± 141 ms (avg ± standard deviation) and T₂ = 60 ± 16 ms (gestational age range 24.3-36.7, median 32.6 weeks). In the statistical model, placental T₂ rose and T₁ remained contant after hyperoxia, and no GA dependency was observed for T₁ or T₂.

DISCUSSION

Well-characterized, motion-robust MRF was used to acquire T₁ and T₂ maps of the placenta. Changes with hyperoxia are consistent with a net increase in oxygen saturation. Toward the goal of whole-placenta quantitative oxygenation imaging over time, we aim to implement 3D MRF with integrated motion correction to improve T₂ accuracy.

Postmortem neuroimaging: Temporal and spatial sensitivity of manganese-enhanced magnetic resonance imaging (MEMRI) and impact of Mn2+ uptake

Magnetic resonance imaging data collection and analysis have been challenges in the field of auditory neuroscience. Recent studies have addressed these concerns by using manganese-enhanced magnetic resonance imaging (MEMRI). Basic challenges for in vivo application of MEMRI in rodents includes how to set inclusion criteria for adequate Mn²⁺ uptake and whether valid data can be collected from brains postmortem. Since brain Mn²⁺ uptake is complete within 2-4 h and clearance can take 2-4 weeks, one assumption has been that Mn²⁺-enhanced R1 values continue to reliably reflect the degree of Mn²⁺-uptake for some indeterminate time after death. To address these issues, the impact of death on R1 values was determined in rats administered Mn²⁺ and rats that were not. Images of auditory nuclei were collected at fixed intervals from rats before and after death for up to 10 h postmortem. By taking a ratio of pituitary and muscle T1-W intensities (P/M), a reliable quantitative method for assessing adequate brain Mn²⁺ uptake was created and suggest that P/M ratios should be adopted to objectively measure the quality of the Mn²⁺ injection. Postmortem R1 values decreased in all brain regions in both the After Mn²⁺ and No Mn²⁺ groups. However, the time-course of postmortem changes in R1 was dependent on brain region and degree of Mn²⁺ uptake. Thus, postmortem R1 values not only differ after death, but vary with time and across brain regions. Postmortem R1 values in unfixed brain tissue, including the auditory nuclei, should be interpreted with caution.

Compromised resting cerebral metabolism after sport-related concussion: A calibrated MRI study

Altered resting cerebral blood flow (CBF₀) in the acute phase post-concussion may contribute to neurobehavioral deficiencies, often reported weeks after the injury. However, in addition to changes in CBF₀, little is known about other physiological mechanisms that may be disturbed within the cerebrovasculature. The aim of this study was to assess whether changes in baseline perfusion following sport-related concussion (SRC) were co-localized with changes in cerebral metabolic demand. Forty-two subjects (15 SRC patients 8.0 ± 4.6 days post-injury and 27 age-matched healthy control athletes) were studied cross-sectionally. CBF₀, cerebrovascular reactivity (CVR), resting oxygen extraction (OEF₀) and cerebral metabolic rate of oxygen consumption (CMRO_2|0) were measured using a combination of hypercapnic and hyperoxic breathing protocols, and the biophysical model developed in calibrated MRI. Blood oxygenation level dependent and perfusion data were acquired simultaneously using a dual-echo arterial spin labelling sequence. SRC patients showed significant decreases in CBF₀ spread across the grey-matter (P < 0.05, corrected), and these differences were also confounded by the effects of baseline end-tidal CO₂ (P < 0.0001). Lower perfusion was co-localized with reductions in regional CMRO_2|0 (P = 0.006) post-SRC, despite finding no group-differences in OEF₀ (P = 0.800). Higher CVR within voxels showing differences in CBF was also observed in the SRC group (P = 0.001), compared to controls. Reductions in metabolic demand despite no significant changes in OEF₀ suggests that hypoperfusion post-SRC may reflect compromised metabolic function after the injury. These results provide novel insight about the possible pathophysiological mechanisms underlying concussion that may affect the clinical recovery of athletes after sport-related head injuries.

Changes in volumetric and metabolic parameters relate to differences in exposure to sub-concussive head impacts

Structural and calibrated magnetic resonance imaging data were acquired on 44 collegiate football players prior to the season (PRE), following the first four weeks in-season (PTC) and one month after the last game (POST). Exposure data collected from g-Force accelerometers mounted to the helmet of each player were used to split participants into HIGH (N = 22) and LOW (N = 22) exposure groups, based on the frequency of impacts sustained by each athlete. Significant decreases in grey-matter volume specific to the HIGH group were documented at POST (P = 0.009), compared to baseline. Changes in resting cerebral blood flow (CBF0), corrected for partial volume effects, were observed within the HIGH group, throughout the season (P < 0.0001), suggesting that alterations in perfusion may follow exposure to sub-concussive collisions. Co-localized significant increases in cerebral metabolic rate of oxygen consumption (CMRO2|0) mid-season were also documented in the HIGH group, with respect to both PRE- and POST values. No physiological changes were observed in the LOW group. Therefore, cerebral metabolic demand may be elevated in players with greater exposure to head impacts. These results provide novel insight into the effects of sub-concussive collisions on brain structure and cerebrovascular physiology and emphasize the importance of multi-modal imaging for a complete characterization of cerebral health.

Morphology of perivascular spaces and enclosed blood vessels in young to middle-aged healthy adults at 7T: Dependences on age, brain region, and breathing gas

Perivascular spaces (PVSs) are fluid-filled spaces surrounding penetrating blood vessels in the brain and are an integral pathway of the glymphatic system. A PVS and the enclosed blood vessel are commonly visualized as a single vessel-like complex (denoted as PVSV) in high-resolution MRI images. Quantitative characterization of the PVSV morphology in MRI images in healthy subjects may serve as a reference for detecting disease related PVS and/or blood vessel alterations in patients with brain diseases.

To this end, we evaluated the age dependences, spatial heterogeneities, and dynamic properties of PVSV morphological features in 45 healthy subjects (21–55 years old), using an ultra-high-resolution three-dimensional transverse relaxation time weighted MRI sequence (0.41 × 0.41 × 0.4 mm³) at 7T. Quantitative PVSV parameters, including apparent diameter, count, volume fraction (VF), and relative contrast to noise ratio (rCNR) were calculated in the white matter and subcortical structures. Dynamic changes were induced by carbogen breathing which are known to induce vasodilation and increase the blood oxygenation level in the brain.

PVSV count and VF significantly increased with age in basal ganglia (BG), so did rCNR in BG, midbrain, and white matter (WM). Apparent PVSV diameter also showed a positive association with age in the three brain regions, although it did not reach statistical significance. The PVSV VF and count showed large inter-subject variations, with coefficients of variation ranging from 0.17 to 0.74 after regressing out age and gender effects. Both apparent diameter and VF exhibited significant spatial heterogeneity, which cannot be explained solely by radio-frequency field inhomogeneities. Carbogen breathing significantly increased VF in BG and WM, and rCNR in thalamus, BG, and WM compared to air breathing.

Our results are consistent with gradual dilation of PVSs with age in healthy adults. The PVSV morphology exhibited spatial heterogeneity and large inter-subject variations and changed during carbogen breathing compared to air breathing.

Determination of oxygen r1 at 3 Tesla using samples with a concentration range of dissolved oxygen

OBJECTIVE

To investigate the sensitivity of modified Look-Locker inversion recovery (MOLLI) to measure changes in dissolved oxygen (DO) concentrations in water samples and to calculate sequence-specific relaxivity (r_1m) and limit of detection (LOD).

MATERIALS AND METHODS

Ten water samples with a range of DO concentrations were scanned at 3 T using two variations of MOLLI schemes. Using linear regression the r₁ of DO was estimated from the measured DO concentrations and T₁ relaxation rates (R₁). The results were combined with previously reported values on in vivo stability measures of the MOLLI sequences and used to estimate a LOD.

RESULTS

DO concentrations ranged from 0.5 to 21.6 mg L^-1. A linear correlation between DO and R₁ was obtained with both MOLLI sequences, with an average correlation coefficient (R²) 0.9 and an average estimated r₁ ([Formula: see text]) of 4.45 × 10^-3 s^-1 mg^-1 L. Estimated LOD was ≈ 10 mg L^-1.

CONCLUSION

MOLLI T₁-mapping sequences may be used for detecting dissolved oxygen in vivo at 3 T with an [Formula: see text] in the range 4.18-4.8 × 10^-3 s^-1 mg^-1 L and a corresponding LOD for dissolved oxygen of approximately 10 mg L^-1. MOLLI-based T₁ mapping may be a useful non-invasive tool for quantification of in vivo changes of DO concentration during oxygen challenges.

Placental MRI: Developing Accurate Quantitative Measures of Oxygenation

The Human Placenta Project has focused attention on the need for noninvasive magnetic resonance imaging (MRI)-based techniques to diagnose and monitor placental function throughout pregnancy. The hope is that the management of placenta-related pathologies would be improved if physicians had more direct, real-time measures of placental health to guide clinical decision making. As oxygen alters signal intensity on MRI and oxygen transport is a key function of the placenta, many of the MRI methods under development are focused on quantifying oxygen transport or oxygen content of the placenta. For example, measurements from blood oxygen level-dependent imaging of the placenta during maternal hyperoxia correspond to outcomes in twin pregnancies, suggesting that some aspects of placental oxygen transport can be monitored by MRI. Additional methods are being developed to accurately quantify baseline placental oxygenation by MRI relaxometry. However, direct validation of placental MRI methods is challenging and therefore animal studies and ex vivo studies of human placentas are needed. Here we provide an overview of the current state of the art of oxygen transport and quantification with MRI. We suggest that as these techniques are being developed, increased focus be placed on ensuring they are robust and reliable across individuals and standardized to enable predictive diagnostic models to be generated from the data. The field is still several years away from establishing the clinical benefit of monitoring placental function in real time with MRI, but the promise of individual personalized diagnosis and monitoring of placental disease in real time continues to motivate this effort.

CMR-based blood oximetry via multi-parametric estimation using multiple T2 measurements

BACKGROUND

Measurement of blood oxygen saturation (O2 saturation) is of great importance for evaluation of patients with many cardiovascular diseases, but currently there are no established non-invasive methods to measure blood O2 saturation in the heart. While T2-based CMR oximetry methods have been previously described, these approaches rely on technique-specific calibration factors that may not generalize across patient populations and are impractical to obtain in individual patients. We present a solution that utilizes multiple T2 measurements made using different inter-echo pulse spacings. These data are jointly processed to estimate all unknown parameters, including O2 saturation, in the Luz-Meiboom (L-M) model. We evaluated the accuracy of the proposed method against invasive catheterization in a porcine hypoxemia model.

METHODS

Sufficient data diversity to estimate the various unknown parameters of the L-M model, including O2 saturation, was achieved by acquiring four T2 maps, each at a different τ 180 (12, 15, 20, and 25 ms). Venous and arterial blood T2 values from these maps, together with hematocrit and arterial O2 saturation, were jointly processed to derive estimates for venous O2 saturation and other nuisance parameters in the L-M model. The technique was validated by a progressive graded hypoxemia experiment in seven pigs. CMR estimates of O2 saturation in the right ventricle were compared against a reference O2 saturation obtained by invasive catheterization from the right atrium in each pig, at each hypoxemia stage. O2 saturation derived from the proposed technique was also compared against the previously described method of applying a global calibration factor (K) to the simplified L-M model.

RESULTS

Venous O2 saturation results obtained using the proposed CMR oximetry method exhibited better agreement (y = 0.84× + 12.29, R2 = 0.89) with invasive blood gas analysis when compared to O2 saturation estimated by a global calibration method (y = 0.69× + 27.52, R2 = 0.73).

CONCLUSIONS

We have demonstrated a novel, non-invasive method to estimate O2 saturation using quantitative T2 mapping. This technique may provide a valuable addition to the diagnostic utility of CMR in patients with congenital heart disease, heart failure, and pulmonary hypertension.

The z-spectrum from human blood at 7T

Chemical Exchange Saturation Transfer (CEST) has been used to assess healthy and pathological tissue in both animals and humans. However, the CEST signal from blood has not been fully assessed. This paper presents the CEST and nuclear Overhauser enhancement (NOE) signals detected in human blood measured via z-spectrum analysis. We assessed the effects of blood oxygenation levels, haematocrit, cell structure and pH upon the z-spectrum in ex vivo human blood for different saturation powers at 7T. The data were analysed using Lorentzian difference (LD) model fitting and AREX (to compensate for changes in T₁), which have been successfully used to study CEST effects in vivo. Full Bloch-McConnell fitting was also performed to provide an initial estimate of exchange rates and transverse relaxation rates of the various pools. CEST and NOE signals were observed at 3.5 ppm, −1.7 ppm and −3.5 ppm and were found to originate primarily from the red blood cells (RBCs), although the amide proton transfer (APT) CEST effect, and NOEs showed no dependence upon oxygenation levels. Upon lysing, the APT and NOE signals fell significantly. Different pH levels in blood resulted in changes in both the APT and NOE (at −3.5 ppm), which suggests that this NOE signal is in part an exchange relayed process. These results will be important for assessing in vivo z-spectra.

The impact of inspired oxygen levels on calibrated fMRI measurements of M, OEF and resting CMRO2 using combined hypercapnia and hyperoxia

Recent calibrated fMRI techniques using combined hypercapnia and hyperoxia allow the mapping of resting cerebral metabolic rate of oxygen (CMRO₂) in absolute units, oxygen extraction fraction (OEF) and calibration parameter M (maximum BOLD). The adoption of such technique necessitates knowledge about the precision and accuracy of the model-derived parameters. One of the factors that may impact the precision and accuracy is the level of oxygen provided during periods of hyperoxia (HO). A high level of oxygen may bring the BOLD responses closer to the maximum M value, and hence reduce the error associated with the M interpolation. However, an increased concentration of paramagnetic oxygen in the inhaled air may result in a larger susceptibility area around the frontal sinuses and nasal cavity. Additionally, a higher O₂ level may generate a larger arterial blood T₁ shortening, which require a bigger cerebral blood flow (CBF) T₁ correction. To evaluate the impact of inspired oxygen levels on M, OEF and CMRO₂ estimates, a cohort of six healthy adults underwent two different protocols: one where 60% of O₂ was administered during HO (low HO or LHO) and one where 100% O₂ was administered (high HO or HHO). The QUantitative O2 (QUO2) MRI approach was employed, where CBF and R2* are simultaneously acquired during periods of hypercapnia (HC) and hyperoxia, using a clinical 3 T scanner. Scan sessions were repeated to assess repeatability of results at the different O₂ levels. Our T₁ values during periods of hyperoxia were estimated based on an empirical ex-vivo relationship between T₁ and the arterial partial pressure of O₂. As expected, our T₁ estimates revealed a larger T₁ shortening in arterial blood when administering 100% O₂ relative to 60% O₂ (T_1LHO = 1.56±0.01 sec vs. T_1HHO = 1.47±0.01 sec, P < 4*10⁻¹³). In regard to the susceptibility artifacts, the patterns and number of affected voxels were comparable irrespective of the O₂ concentration. Finally, the model-derived estimates were consistent regardless of the HO levels, indicating that the different effects are adequately accounted for within the model.