Cerebral metabolic rate of oxygen (CMRO2 ) mapping by combining quantitative susceptibility mapping (QSM) and quantitative blood oxygenation level-dependent imaging (qBOLD)

Quantitative BOLD (qBOLD) imaging of oxygen metabolism and blood oxygenation in the human body: A scoping review

PURPOSE

There are many approaches to the quantitative BOLD (qBOLD) technique described in the literature, differing in pulse sequences, MRI parameters and data processing. Thus, in this review, we summarized the acquisition methods, approaches used for oxygenation quantification and clinical populations investigated.

METHODS

Three databases were systematically searched (Medline, Embase, and Web of Science) for published research that used qBOLD methods for quantification of oxygen metabolism. Data extraction and synthesis were performed by one author and reviewed by a second author.

RESULTS

A total of 93 relevant papers were identified. Acquisition strategies were summarized, and oxygenation parameters were found to have been investigated in many pathologies such as steno-occlusive diseases, stroke, glioma, and multiple sclerosis disease.

CONCLUSION

A summary of qBOLD approaches for oxygenation measurements and applications could help researchers to identify good practice and provide objective information to inform the development of future consensus recommendations.

Vessel-specific quantification of cerebral venous oxygenation with velocity-encoding preparation and rapid acquisition

PURPOSE

Non-invasive measurement of cerebral venous oxygenation (Yv) is of critical importance in brain diseases. The present work proposed a fast method to quantify regional Yv map for both large and small veins.

METHODS

A new sequence was developed, referred to as TRU-VERA (T2 relaxation under velocity encoding and rapid acquisition, which isolates blood spins from static tissue with velocity-encoding preparation, modulates the T2 weighting of venous signal with T2-preparation and utilizes a bSSFP readout to achieve fast acquisition with high resolution. The sequence was first optimized to achieve best sensitivity for both large and small veins, and then validated with TRUST (T2 relaxation under spin tagging), TRUPC (T2 relaxation under phase contrast), and accelerated TRUPC MRI. Regional difference of Yv was evaluated, and test-retest reproducibility was examined.

RESULTS

Optimal Venc was determined to be 3 cm/s, while recovery time and balanced SSFP flip angle within reasonable range had minimal effect on SNR efficiency. Venous T2 measured with TRU-VERA was highly correlated with T2 from TRUST (R2 = 0.90), and a conversion equation was established for further calibration to Yv. TRU-VERA sequences showed consistent Yv estimation with TRUPC (R2 = 0.64) and accelerated TRUPC (R2 = 0.79). Coefficient of variation was 0.84% for large veins and 2.49% for small veins, suggesting an excellent test-retest reproducibility.

CONCLUSION

The proposed TRU-VERA sequence is a promising method for vessel-specific oxygenation assessment.

Reduced oxygen extraction fraction in deep cerebral veins associated with cognitive impairment in multiple sclerosis

Studying the relationship between cerebral oxygen utilization and cognitive impairment is essential to understanding neuronal functional changes in the disease progression of multiple sclerosis (MS). This study explores the potential of using venous susceptibility in internal cerebral veins (ICVs) as an imaging biomarker for cognitive impairment in relapsing-remitting MS (RRMS) patients. Quantitative susceptibility mapping derived from fully flow-compensated MRI phase data was employed to directly measure venous blood oxygen saturation levels (SvO2) in the ICVs. Results revealed a significant reduction in the susceptibility of ICVs (212.4 ± 30.8 ppb vs 239.4 ± 25.9 ppb) and a significant increase of SvO2 (74.5 ± 1.89% vs 72.4 ± 2.23%) in patients with RRMS compared with age- and sex-matched healthy controls. Both the susceptibility of ICVs (r = 0.508, p = 0.031) and the SvO2 (r = -0.498, p = 0.036) exhibited a moderate correlation with cognitive decline in these patients assessed by the Paced Auditory Serial Addition Test, while no significant correlation was observed with clinical disability measured by the Expanded Disability Status Scale. The findings suggest that venous susceptibility in ICVs has the potential to serve as a specific indicator of oxygen metabolism and cognitive function in RRMS. .

Ischemic perfusion radiomics: assessing neurological impairment in acute ischemic stroke

Introduction

Accurate neurological impairment assessment is crucial for the clinical treatment and prognosis of patients with acute ischemic stroke (AIS). However, the original perfusion parameters lack the deep information for characterizing neurological impairment, leading to difficulty in accurate assessment. Given the advantages of radiomics technology in feature representation, this technology should provide more information for characterizing neurological impairment. Therefore, with its rigorous methodology, this study offers practical implications for clinical diagnosis by exploring the role of ischemic perfusion radiomics features in assessing the degree of neurological impairment.

Methods

This study employs a meticulous methodology, starting with generating perfusion parameter maps through Dynamic Susceptibility Contrast-Perfusion Weighted Imaging (DSC-PWI) and determining ischemic regions based on these maps and a set threshold. Radiomics features are then extracted from the ischemic regions, and the t-test and least absolute shrinkage and selection operator (Lasso) algorithms are used to select the relevant features. Finally, the selected radiomics features and machine learning techniques are used to assess the degree of neurological impairment in AIS patients.

Results

The results show that the proposed method outperforms the original perfusion parameters, radiomics features of the infarct and hypoxic regions, and their combinations, achieving an accuracy of 0.926, sensitivity of 0.923, specificity of 0.929, PPV of 0.923, NPV of 0.929, and AUC of 0.923, respectively.

Conclusion

The proposed method effectively assesses the degree of neurological impairment in AIS patients, providing an objective auxiliary assessment tool for clinical diagnosis.

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.

Validation of a new 3D quantitative BOLD based cerebral oxygen extraction mapping

Quantitative BOLD (qBOLD) MRI allows evaluation of oxidative metabolism of the brain based purely on an endogenous contrast mechanism. The method quantifies deoxygenated blood volume (DBV) and hemoglobin oxygen saturation level of venous blood (Yv), yielding oxygen extraction fraction (OEF), and along with a separate measurement of cerebral blood flow, cerebral metabolic rate of oxygen (CMRO2) maps. Here, we evaluated our recently reported 3D qBOLD method that rectifies a number of deficiencies in prior qBOLD approaches in terms of repeat reproducibility and sensitivity to hypercapnia on the metabolic parameters, and in comparison to dual-gas calibrated BOLD (cBOLD) MRI for determining resting-state oxygen metabolism. Results suggested no significant difference between test-retest qBOLD scans in either DBV and OEF. Exposure to hypercapnia yielded group averages of 38 and 28% for OEF and 151 and 146 µmol/min/100 g for CMRO2 in gray matter at baseline and hypercapnia, respectively. The decrease of OEF during hypercapnia was significant (p ≪ 0.01), whereas CMRO2 did not change significantly (p = 0.25). Finally, baseline OEF (37 vs. 39%) and CMRO2 (153 vs. 145 µmol/min/100 g) in gray matter using qBOLD and dual-gas cBOLD were found to be in good agreement with literature values, and were not significantly different from each other (p > 0.1).

Unveiling metabolic patterns in dementia: Insights from high-resolution quantitative blood-oxygenation-level-dependent MRI

BACKGROUND

Oxygen extraction fraction (OEF) and deoxyhemoglobin (DoHb) levels reflect variations in cerebral oxygen metabolism in demented patients.

PURPOSE

Delineating the metabolic profiles evident throughout different phases of dementia necessitates an integrated analysis of OEF and DoHb levels. This is enabled by leveraging high-resolution quantitative blood oxygenation level dependent (qBOLD) analysis of magnitude images obtained from a multi-echo gradient-echo MRI (mGRE) scan performed on a 3.0 Tesla scanner.

METHODS

Achieving superior spatial resolution in qBOLD necessitates the utilization of an mGRE scan with only four echoes, which in turn limits the number of measurements compared to the parameters within the qBOLD model. Consequently, it becomes imperative to discard non-essential parameters to facilitate further analysis. This process entails transforming the qBOLD model into a format suitable for fitting the log-magnitude difference (L-MDif) profiles of the four echo magnitudes present in each brain voxel. In order to bolster spatial specificity, the log-difference qBOLD model undergoes refinement into a representative form, termed as r-qBOLD, particularly when applied to class-averaged L-MDif signals derived through k-means clustering of L-MDif signals from all brain voxels into a predetermined number of clusters. The agreement between parameters estimated using r-qBOLD for different cluster sizes is validated using Bland-Altman analysis, and the model's goodness-of-fit is evaluated using aχ 2 ${\chi ^2}$ -test. Retrospective MRI data of Alzheimer's disease (AD), mild cognitive impairment (MCI), and non-demented patients without neuropathological disorders, pacemakers, other implants, or psychiatric disorders, who completed a minimum of three visits prior to MRI enrolment, are utilized for the study.

RESULTS

Utilizing a cohort comprising 30 demented patients aged 65-83 years in stages 4-6 representing mild, moderate, and severe stages according to the clinical dementia rating (CDR), matched with an age-matched non-demented control group of 18 individuals, we conducted joint observations of OEF and DoHb levels estimated using r-qBOLD. The observations elucidate metabolic signatures in dementia based on OEF and DoHb levels in each voxel. Our principal findings highlight the significance of spatial patterns of metabolic profiles (metabolic patterns) within two distinct regimes: OEF levels exceeding the normal range (S1-regime), and OEF levels below the normal range (S2-regime). The S1-regime, accompanied by low DoHb levels, predominantly manifests in fronto-parietal and perivascular regions with increase in dementia severity. Conversely, the S2-regime, accompanied by low DoHb levels, is observed in medial temporal (MTL) regions. Other regions with abnormal metabolic patterns included the orbitofrontal cortex (OFC), medial-orbital prefrontal cortex (MOPFC), hypothalamus, ventro-medial prefrontal cortex (VMPFC), and retrosplenial cortex (RSP). Dysfunction in the OFC and MOPFC indicated cognitive and emotional impairment, while hypothalamic involvement potentially indicated preclinical dementia. Reduced metabolic activity in the RSP suggested early-stage AD related functional abnormalities.

CONCLUSIONS

Integrated analysis of OEF and DoHb levels using r-qBOLD reveals distinct metabolic signatures across dementia phases, highlighting regions susceptible to neuronal loss, vascular involvement, and preclinical indicators.

Spatiotemporal patterns of brain iron-oxygen metabolism in patients with Parkinson's disease

OBJECTIVES

Iron deposition and mitochondrial dysfunction are closely associated with the genesis and progression of Parkinson's disease (PD). This study aims to extract susceptibility and oxygen extraction fraction (OEF) values of deep grey matter (DGM) to explore spatiotemporal progression patterns of brain iron-oxygen metabolism in PD.

METHODS

Ninety-five PD patients and forty healthy controls (HCs) were included. Quantitative susceptibility mapping (QSM) and OEF maps were computed from MRI multi-echo gradient echo data. Analysis of covariance (ANCOVA) was used to compare mean susceptibility and OEF values in DGM between early-stage PD (ESP), advanced-stage PD (ASP) patients and HCs. Then Granger causality analysis on the pseudo-time-series of MRI data was applied to assess the causal effect of early altered nuclei on iron content and oxygen extraction in other DGM nuclei.

RESULTS

The susceptibility values in substantia nigra (SN), red nucleus, and globus pallidus (GP) significantly increased in PD patients compared with HCs, while the iron content in GP did not elevate obviously until the late stage. The mean OEF values for the caudate nucleus, putamen, and dentate nucleus were higher in ESP patients than in ASP patients or/and HCs. We also found that iron accumulation progressively expands from the midbrain to the striatum. These alterations were correlated with clinical features and improved AUC for early PD diagnosis to 0.824.

CONCLUSIONS

Abnormal cerebral iron deposition and tissue oxygen utilization in PD measured by QSM and OEF maps could reflect pathological alterations in neurodegenerative processes and provide valuable indicators for disease identification and management.

CLINICAL RELEVANCE STATEMENT

Noninvasive assessment of cerebral iron-oxygen metabolism may serve as clinical evidence of pathological changes in PD and improve the validity of diagnosis and disease monitoring.

KEY POINTS

• Quantitative susceptibility mapping and oxygen extraction fraction maps indicated the cerebral pathology of abnormal iron accumulation and oxygen metabolism in Parkinson's disease. • Iron deposition is mainly in the midbrain, while altered oxygen metabolism is concentrated in the striatum and cerebellum. • The susceptibility and oxygen extraction fraction values in subcortical nuclei were associated with clinical severity.

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.

Multi-Echo Complex Quantitative Susceptibility Mapping and Quantitative Blood Oxygen Level-Dependent Magnitude (mcQSM + qBOLD or mcQQ) for Oxygen Extraction Fraction (OEF) Mapping

Oxygen extraction fraction (OEF), the fraction of oxygen that tissue extracts from blood, is an essential biomarker used to directly assess tissue viability and function in neurologic disorders. In ischemic stroke, for example, increased OEF can indicate the presence of penumbra-tissue with low perfusion yet intact cellular integrity-making it a primary therapeutic target. However, practical OEF mapping methods are not currently available in clinical settings, owing to the impractical data acquisitions in positron emission tomography (PET) and the limitations of existing MRI techniques. Recently, a novel MRI-based OEF mapping technique, termed QQ, was proposed. It shows high potential for clinical use by utilizing a routine sequence and removing the need for impractical multiple gas inhalations. However, QQ relies on the assumptions of Gaussian noise in susceptibility and multi-echo gradient echo (mGRE) magnitude signals for OEF estimation. This assumption is unreliable in low signal-to-noise ratio (SNR) regions like disease-related lesions, risking inaccurate OEF estimation and potentially impacting clinical decisions. Addressing this, our study presents a novel multi-echo complex QQ (mcQQ) that models realistic Gaussian noise in mGRE complex signals. We implemented mcQQ using a deep learning framework (mcQQ-NET) and compared it with the existing QQ-NET in simulations, ischemic stroke patients, and healthy subjects, using identical training and testing datasets and schemes. In simulations, mcQQ-NET provided more accurate OEF than QQ-NET. In the subacute stroke patients, mcQQ-NET showed a lower average OEF ratio in lesions relative to unaffected contralateral normal tissue than QQ-NET. In the healthy subjects, mcQQ-NET provided uniform OEF maps, similar to QQ-NET, but without unrealistically high OEF outliers in areas of low SNR, such as SNR ≤ 15 (dB). Therefore, mcQQ-NET improves OEF accuracy by more accurately reflecting realistic Gaussian noise in complex mGRE signals. Its enhanced sensitivity to OEF abnormalities, based on more realistic biophysics modeling, suggests that mcQQ-NET has potential for investigating tissue variability in neurologic disorders.

Reduced Oxygen Extraction Fraction in Deep Cerebral Veins Associated with Cognitive Impairment in Multiple Sclerosis

Studying the relationship between cerebral oxygen utilization and cognitive impairment is essential to understanding neuronal functional changes in the disease progression of multiple sclerosis (MS). This study explores the potential of using venous susceptibility in internal cerebral veins (ICVs) as an imaging biomarker for cognitive impairment in relapsing-remitting MS (RRMS) patients. Quantitative susceptibility mapping derived from fully flow-compensated MRI phase data was employed to directly measure venous blood oxygen saturation levels (SvO2) in the ICVs. Results revealed a significant reduction in the susceptibility of ICVs (212.4 ± 30.8 ppb vs 239.4 ± 25.9 ppb) and a significant increase of SvO2 (74.5 ± 1.89 % vs 72.4 ± 2.23 %) in patients with RRMS compared with age- and sex-matched healthy controls. Both the susceptibility of ICVs (r = 0.646, p = 0.004) and the SvO2 (r = -0.603, p = 0.008) exhibited a strong correlation with cognitive decline in these patients assessed by the Paced Auditory Serial Addition Test, while no significant correlation was observed with clinical disability measured by the Expanded Disability Status Scale. The findings suggest that venous susceptibility in ICVs has the potential to serve as a specific indicator of oxygen metabolism and cognitive function in RRMS.

Evaluation of whole-brain oxygen metabolism in Alzheimer's disease using QSM and quantitative BOLD

OBJECTIVE

The objective of this study was to evaluate the whole-brain pattern of oxygen extraction fraction (OEF), cerebral blood flow (CBF), and cerebral metabolic rate of oxygen consumption (CMRO2) perturbation in Alzheimer's disease (AD) and investigate the relationship between regional cerebral oxygen metabolism and global cognition.

METHODS

Twenty-six AD patients and 25 age-matched healthy controls (HC) were prospectively recruited in this study. Mini-Mental State Examination (MMSE) was used to evaluate cognitive status. We applied the QQ-CCTV algorithm which combines quantitative susceptibility mapping and quantitative blood oxygen level-dependent models (QQ) for OEF calculation. CBF map was computed from arterial spin labeling and CMRO2 was generated based on Fick's principle. Whole-brain and regional OEF, CBF, and CMRO2 analyses were performed. The associations between these measures in substructures of deep brain gray matter and MMSE scores were assessed.

RESULTS

Whole brain voxel-wise analysis showed that CBF and CMRO2 values significantly decreased in AD predominantly in the bilateral angular gyrus, precuneus gyrus and parieto-temporal regions. Regional analysis showed that CBF value decreased in the bilateral caudal hippocampus and left rostral hippocampus and CMRO2 value decreased in left caudal and rostral hippocampus in AD patients. Considering all subjects in the AD and HC groups combined, the mean CBF and CMRO2 values in the bilateral hippocampus positively correlated with the MMSE score.

CONCLUSION

CMRO2 mapping with the QQ-CCTV method - which is readily available in MR systems for clinical practice - can be a potential biomarker for AD. In addition, CMRO2 in the hippocampus may be a useful tool for monitoring cognitive impairment.

Clinical targeting of the cerebral oxygen cascade to improve brain oxygenation in patients with hypoxic-ischaemic brain injury after cardiac arrest

The cerebral oxygen cascade includes three key stages: (a) convective oxygen delivery representing the bulk flow of oxygen to the cerebral vascular bed; (b) diffusion of oxygen from the blood into brain tissue; and (c) cellular utilisation of oxygen for aerobic metabolism. All three stages may become dysfunctional after resuscitation from cardiac arrest and contribute to hypoxic-ischaemic brain injury (HIBI). Improving convective cerebral oxygen delivery by optimising cerebral blood flow has been widely investigated as a strategy to mitigate HIBI. However, clinical trials aimed at optimising convective oxygen delivery have yielded neutral results. Advances in the understanding of HIBI pathophysiology suggest that impairments in the stages of the oxygen cascade pertaining to oxygen diffusion and cellular utilisation of oxygen should also be considered in identifying therapeutic strategies for the clinical management of HIBI patients. Culprit mechanisms for these impairments may include a widening of the diffusion barrier due to peri-vascular oedema and mitochondrial dysfunction. An integrated approach encompassing both intra-parenchymal and non-invasive neuromonitoring techniques may aid in detecting pathophysiologic changes in the oxygen cascade and enable patient-specific management aimed at reducing the severity of HIBI.

Evaluating Physiological MRI Parameters in Patients with Brain Metastases Undergoing Stereotactic Radiosurgery-A Preliminary Analysis and Case Report

Brain metastases occur in ten to thirty percent of the adult cancer population. Treatment consists of different (palliative) options, including stereotactic radiosurgery (SRS). Sensitive MRI biomarkers are needed to better understand radiotherapy-related effects on cerebral physiology and the subsequent effects on neurocognitive functioning. In the current study, we used physiological imaging techniques to assess cerebral blood flow (CBF), oxygen extraction fraction (OEF), cerebral metabolic rate of oxygen (CMRO2) and cerebrovascular reactivity (CVR) before and three months after SRS in nine patients with brain metastases. The results showed improvement in OEF, CBF and CMRO2 within brain tissue that recovered from edema (all p ≤ 0.04), while CVR remained impacted. We observed a global post-radiotherapy increase in CBF in healthy-appearing brain tissue (p = 0.02). A repeated measures correlation analysis showed larger reductions within regions exposed to higher radiotherapy doses in CBF (rrm = -0.286, p < 0.001), CMRO2 (rrm = -0.254, p < 0.001), and CVR (rrm = -0.346, p < 0.001), but not in OEF (rrm = -0.004, p = 0.954). Case analyses illustrated the impact of brain metastases progression on the post-radiotherapy changes in both physiological MRI measures and cognitive performance. Our preliminary findings suggest no radiotherapy effects on physiological parameters occurred in healthy-appearing brain tissue within 3-months post-radiotherapy. Nevertheless, as radiotherapy can have late side effects, larger patient samples allowing meaningful grouping of patients and longer follow-ups are needed.

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.

Multimodal comparisons of QSM and PET in neurodegeneration and aging

Quantitative susceptibility mapping (QSM) has been used to study susceptibility changes that may occur based on tissue composition and mineral deposition. Iron is a primary contributor to changes in magnetic susceptibility and of particular interest in applications of QSM to neurodegeneration and aging. Iron can contribute to neurodegeneration through inflammatory processes and via interaction with aggregation of disease-related proteins. To better understand the local susceptibility changes observed on QSM, its signal has been studied in association with other imaging metrics such as positron emission tomography (PET). The associations of QSM and PET may provide insight into the pathophysiology of disease processes, such as the role of iron in aging and neurodegeneration, and help to determine the diagnostic utility of QSM as an indirect indicator of disease processes typically evaluated with PET. In this review we discuss the proposed mechanisms and summarize prior studies of the associations of QSM and amyloid PET, tau PET, TSPO PET, FDG-PET, 15O-PET, and F-DOPA PET in evaluation of neurologic diseases with a focus on aging and neurodegeneration.

Cerebral oxygen extraction fraction declines with ventricular enlargement in patients with normal pressure hydrocephalus

OBJECTIVE

Normal pressure hydrocephalus (NPH) is a neurodegenerative disease that is potentially reversible by shunt surgery in approximately 60% of patients. Imaging may provide a means to investigate brain tissue viability and oxygen metabolism in NPH patients.

METHODS

Oxygen extraction fraction (OEF) mapping was generated from 3D multi-echo gradient echo MRI (mGRE) data using QQ-CCTV algorithm and cerebral blood flow (CBF) using 3D arterial spin labeling (ASL) MRI data, thereby calculating the cerebral metabolic rate of oxygen (CMRO2 = CBF × OEF × [H]a) in 16 NPH patients. Regression analyses using cortical gray matter and deep gray matter regions were conducted with age, gender, CSF stroke volume and normalized ventricular volume as independent variables.

RESULTS

OEF showed significant negative correlations with normalized brain ventricular volumes in the whole brain (p = 0.004, q = 0.01), cortical gray matter (p = 0.004, q = 0.01), caudate (p = 0.02, q = 0.04), and pallidum (p = 0.03, q = 0.04), but no significant correlation with CSF stroke volume (q > 0.05). There was no significant finding with CBF or CMRO2.

CONCLUSION

In NPH patients, low OEF in several regions was significantly correlated with large ventricular volumes, indicating decreased tissue oxygen metabolism with increased NPH severity. OEF mapping may provide a functional understanding of neurodegeneration in NPH and may improve monitoring of disease course and treatment outcomes.

Assessing Cerebral Oxygen Metabolism Changes in Patients With Preeclampsia Using Voxel-Based Morphometry of Oxygen Extraction Fraction Maps in Magnetic Resonance Imaging

OBJECTIVE

The objective of this study was to analyze the different brain oxygen metabolism statuses in preeclampsia using magnetic resonance imaging and investigate the factors that affect cerebral oxygen metabolism in preeclampsia.

MATERIALS AND METHODS

Forty-nine women with preeclampsia (mean age 32.4 years; range, 18-44 years), 22 pregnant healthy controls (PHCs) (mean age 30.7 years; range, 23-40 years), and 40 non-pregnant healthy controls (NPHCs) (mean age 32.5 years; range, 20-42 years) were included in this study. Brain oxygen extraction fraction (OEF) values were computed using quantitative susceptibility mapping (QSM) plus quantitative blood oxygen level-dependent magnitude-based OEF mapping (QSM + quantitative blood oxygen level-dependent imaging or QQ) obtained with a 1.5-T scanner. Voxel-based morphometry (VBM) was used to investigate the differences in OEF values in the brain regions among the groups.

RESULTS

Among the three groups, the average OEF values were significantly different in multiple brain areas, including the parahippocampus, multiple gyri of the frontal lobe, calcarine, cuneus, and precuneus (all P-values were less than 0.05, after correcting for multiple comparisons). The average OEF values of the preeclampsia group were higher than those of the PHC and NPHC groups. The bilateral superior frontal gyrus/bilateral medial superior frontal gyrus had the largest size of the aforementioned brain regions, and the OEF values in this area were 24.2 ± 4.6, 21.3 ± 2.4, and 20.6 ± 2.8 in the preeclampsia, PHC, and NPHC groups, respectively. In addition, the OEF values showed no significant differences between NPHC and PHC. Correlation analysis revealed that the OEF values of some brain regions (mainly involving the frontal, occipital, and temporal gyrus) were positively correlated with age, gestational week, body mass index, and mean blood pressure in the preeclampsia group (r = 0.361-0.812).

CONCLUSION

Using whole-brain VBM analysis, we found that patients with preeclampsia had higher OEF values than controls.

MRI assessment of cerebral oxygen extraction fraction in the medial temporal lobe

The medial temporal lobe (MTL) is a key area implicated in many brain diseases, such as Alzheimer's disease. As a functional biomarker, the oxygen extraction fraction (OEF) of MTL may be more sensitive than structural atrophy of MTL, especially at the early stages of diseases. However, there is a lack of non-invasive techniques to measure MTL-OEF in humans. The goal of this work is to develop an MRI technique to assess MTL-OEF in a clinically practical time without using contrast agents. The proposed method measures venous oxygenation (Yv) in the basal veins of Rosenthal (BVs), which are the major draining veins of the MTL. MTL-OEF can then be estimated as the arterio-venous difference in oxygenation. We developed an MRI sequence, dubbed arterial-suppressed accelerated T2-relaxation-under-phase-contrast (AS-aTRUPC), to quantify the blood T2 of the BVs, which was then converted to Yv through a well-established calibration model. MTL-OEF was calculated as (Ya-Yv)/Ya × 100%, where Ya was the arterial oxygenation. The feasibility of AS-aTRUPC to quantify MTL-OEF was evaluated in 16 healthy adults. The sensitivity of AS-aTRUPC in detecting OEF changes was assessed by a caffeine ingestion (200 mg) challenge. For comparison, T2-relaxation-under-spin-tagging (TRUST) MRI, which is a widely used global OEF technique, was also acquired. The dependence of MTL-OEF on age was examined by including another seven healthy elderly subjects. The results showed that in healthy adults, MTL-OEF of the left and right hemispheres were correlated (P=0.005). MTL-OEF was measured to be 23.9±3.6% (mean±standard deviation) and was significantly lower (P<0.0001) than the OEF of 33.3±2.9% measured in superior sagittal sinus (SSS). After caffeine ingestion, there was an absolute percentage increase of 9.1±4.0% in MTL-OEF. Additionally, OEF in SSS measured with AS-aTRUPC showed a strong correlation with TRUST OEF (intra-class correlation coefficient=0.94 with 95% confidence interval [0.91, 0.96]), with no significant bias (P=0.12). MTL-OEF was found to increase with age (MTL-OEF=20.997+0.100 × age; P=0.02). In conclusion, AS-aTRUPC MRI provides non-invasive assessments of MTL-OEF and may facilitate future clinical applications of MTL-OEF as a disease biomarker.

Measurement of Cerebral Oxygen Extraction Fraction Using Quantitative BOLD Approach: A Review

Quantification of brain oxygenation and metabolism, both of which are indicators of the level of brain activity, plays a vital role in understanding the cerebral perfusion and the pathophysiology of brain disorders. Magnetic resonance imaging (MRI), a widely used clinical imaging technique, which is very sensitive to magnetic susceptibility, has the possibility of substituting positron emission tomography (PET) in measuring oxygen metabolism. This review mainly focuses on the quantitative blood oxygenation level-dependent (qBOLD) method for the evaluation of oxygen extraction fraction (OEF) in the brain. Here, we review the theoretic basis of qBOLD, as well as existing acquisition and quantification methods. Some published clinical studies are also presented, and the pros and cons of qBOLD method are discussed as well.

Oxygen extraction fraction (OEF) assesses cerebral oxygen metabolism of deep gray matter in patients with pre-eclampsia

OBJECTIVES

The objective of this study was to compare oxygen extraction fraction (OEF) values in the deep gray matter (GM) of pre-eclampsia (PE) patients, pregnant healthy controls (PHCs), and non-pregnant healthy controls (NPHCs) to explore their brain oxygen metabolism differences in GM.

METHODS

Forty-seven PE patients, forty NPHCs, and twenty-one PHCs were included. Brain OEF values were computed from quantitative susceptibility mapping (QSM) plus quantitative blood oxygen level-dependent magnitude (QSM + qBOLD = QQ)-based mapping. One-way ANOVA was used to compare mean OEF values in the three groups. The area under the curve of the mean OEF value in each region of interest was estimated using a receiver operating characteristic curve analysis.

RESULTS

We found that the mean OEF values in the thalamus, putamen, caudate nucleus, pallidum, and substantia nigra were significantly different in these three groups (F = 5.867, p = 0.004; F = 5.142, p = 0007; F = 6.158, p = 0.003; F = 6.319, p = 0.003; F = 5.491, p = 0.005). The mean OEF values for these 5 regions were higher in PE patients than in NPHCs and in PHCs (p < 0.05). The AUC of these ROIs ranged from 0.673 to 0.692 (p < 0.01) and cutoff values varied from 35.1 to 36.6%, indicating that the OEF values could discriminate patients with and without PE. Stepwise multivariate analysis revealed that the OEF values correlated with hematocrit in pregnant women (r = 0.353, p = 0.003).

CONCLUSION

OEF values in the brains of pregnant women can be measured in clinical practice using QQ-based OEF mapping for noninvasive assessment of hypertensive disorders.

KEY POINTS

• Pre-eclampsia is a hypertensive disorder associated with abnormalities in brain oxygen extraction. • Oxygen extraction fraction (OEF) is an indicator of brain tissue viability and function. QQ-based mapping of OEF is a new MRI technique that can noninvasively quantify brain oxygen metabolism. • OEF values in the brains of pregnant women can be measured for noninvasive assessment of hypertensive disorders in clinical practice.

Cerebral Venous Oxygen Saturation in Hypoperfusion Regions May Become a New Imaging Indicator to Predict the Clinical Outcome of Stroke

To automatically and quantitatively evaluate the venous oxygen saturation (SvO2) in cerebral ischemic tissues and explore its value in predicting prognosis. A retrospective study was conducted on 48 AIS patients hospitalized in our hospital from 2015−2018. Based on quantitative susceptibility mapping and perfusion-weighted imaging, this paper measured the cerebral SvO2 in hypoperfusion tissues and its change after intraarterial rt-PA treatment. The cerebral SvO2 in different hypoperfusion regions between the favorable and unfavorable clinical outcome groups was analyzed using an independent t-test. Relationships between cerebral SvO2 and clinical scores were determined using the Pearson correlation coefficient. The receiver operating characteristic process was conducted to evaluate the accuracy of cerebral SvO2 in predicting unfavorable clinical outcomes. Cerebral SvO2 in hypoperfusion (Tmax > 4 and 6 s) was significantly different between the two groups at follow-up (p < 0.05). Cerebral SvO2 and its changes before and after treatment were negatively correlated with clinical scores. The positive predictive value, negative predictive value, accuracy, and area under the curve of the cerebral SvO2 were higher than those predicted by the ischemic core. Therefore, the cerebral SvO2 of hypoperfusion regions was a stronger imaging predictor of unfavorable clinical outcomes after stroke.

Non-Invasive measurement of the cerebral metabolic rate of oxygen using MRI in rodents

Malfunctions of oxygen metabolism are suspected to play a key role in a number of neurological and psychiatric disorders, but this hypothesis cannot be properly investigated without an in-vivo non-invasive measurement of brain oxygen consumption. We present a new way to measure the Cerebral Metabolic Rate of Oxygen (CMRO2) by combining two existing magnetic resonance imaging techniques, namely arterial spin-labelling and oxygen extraction fraction mapping. This method was validated by imaging rats under different anaesthetic regimes and was strongly correlated to glucose consumption measured by autoradiography.

Non-Invasive measurement of the cerebral metabolic rate of oxygen using MRI in rodents

Malfunctions of oxygen metabolism are suspected to play a key role in a number of neurological and psychiatric disorders, but this hypothesis cannot be properly investigated without an in-vivo non-invasive measurement of brain oxygen consumption. We present a new way to measure the Cerebral Metabolic Rate of Oxygen (CMRO 2) by combining two existing magnetic resonance imaging techniques, namely arterial spin-labelling and oxygen extraction fraction mapping. This method was validated by imaging rats under different anaesthetic regimes and was strongly correlated to glucose consumption measured by autoradiography.

Non-Invasive measurement of the cerebral metabolic rate of oxygen using MRI in rodents

Malfunctions of oxygen metabolism are suspected to play a key role in a number of neurological and psychiatric disorders, but this hypothesis cannot be properly investigated without an in-vivo non-invasive measurement of brain oxygen consumption. We present a new way to measure the Cerebral Metabolic Rate of Oxygen (CMRO2) by combining two existing magnetic resonance imaging techniques, namely arterial spin-labelling and oxygen extraction fraction mapping. This method was validated by imaging rats under different anaesthetic regimes and was strongly correlated to glucose consumption measured by autoradiography.

T2 -oximetry-based cerebral venous oxygenation mapping using Fourier-transform-based velocity-selective pulse trains

Purpose

To develop a T₂‐oximetry method for quantitative mapping of cerebral venous oxygenation fraction (Y_v) using Fourier‐transform–based velocity‐selective (FT‐VS) pulse trains.

Methods

The venous isolation preparation was achieved by using an FT‐VS inversion plus a nonselective inversion (NSI) pulse to null the arterial blood signal while minimally affected capillary blood flows out into the venular vasculature during the outflow time (TO), and then applying an Fourier transform based velocity selective saturation (FT‐VSS) pulse to suppress the tissue signal. A multi‐echo readout was employed to obtain venous T₂ (T_2,v) efficiently with the last echo used to detect the residual CSF signal and correct its contamination in the fitting. Here we compared the performance of this FT‐VS–based venous isolation preparations with a traditional velocity‐selective saturation (VSS)–based approach (quantitative imaging of extraction of oxygen and tissue consumption [QUIXOTIC]) with different cutoff velocities for Y_v mapping on 6 healthy volunteers at 3 Tesla.

Results

The FT‐VS–based methods yielded higher venous blood signal and temporal SNR with less CSF contamination than the velocity‐selective saturation–based results. The averaged Y_v values across the whole slice measured in different experiments were close to the global Y_v measured from the individual internal jugular vein.

Conclusion

The feasibility of the FT‐VS–based Y_v estimation was demonstrated on healthy volunteers. The obtained high venous signal as well as the mitigation of CSF contamination led to a good agreement between the T_2,v and Y_v measured in the proposed method with the values in the literature.

Click here for author‐reader discussions

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.

Brain oxygen extraction and neural tissue susceptibility are associated with cognitive impairment in older individuals

BACKGROUND AND PURPOSE

We investigated the effects of aging, white matter hyperintensities (WMH), and cognitive impairment on brain iron levels and cerebral oxygen metabolism, known to be altered in Alzheimer's disease (AD), using quantitative susceptibility mapping and MR-based cerebral oxygen extraction fraction (OEF).

METHODS

In 100 individuals over the age of 50 (68/32 cognitively impaired/intact), OEF and neural tissue susceptibility (χ_n ) were computed retrospectively from MRI multi-echo gradient echo data, obtained on a 3 Tesla MRI scanner. The effects of age and WMH on OEF and χ_n were assessed within groups, and OEF and χ_n were assessed between groups, using multivariate regression analyses.

RESULTS

Cognitively impaired subjects were found to have 19% higher OEF and 34% higher χ_n than cognitively intact subjects in the cortical gray matter and several frontal, temporal, and parietal regions (p < .05). Increased WMH burden was significantly associated with decreased OEF in the cognitively impaired, but not in the cognitively intact. Older age had a stronger association with decreased OEF in the cognitively intact group. Both older age and increased WMH burden were significantly associated with increased χ_n in temporoparietal regions in the cognitively impaired.

CONCLUSIONS

Higher brain OEF and χ_n in cognitively impaired older individuals may reflect altered oxygen metabolism and iron in areas with underlying AD pathology. Both age and WMH have associations with OEF and χ_n but are modified by the presence of cognitive impairment.

Learning-based motion artifact removal networks for quantitative R 2 ∗ mapping

PURPOSE

To introduce two novel learning-based motion artifact removal networks (LEARN) for the estimation of quantitative motion- and B 0 -inhomogeneity-corrected R 2 ∗ maps from motion-corrupted multi-Gradient-Recalled Echo (mGRE) MRI data.

METHODS

We train two convolutional neural networks (CNNs) to correct motion artifacts for high-quality estimation of quantitative B 0 -inhomogeneity-corrected R 2 ∗ maps from mGRE sequences. The first CNN, LEARN-IMG, performs motion correction on complex mGRE images, to enable the subsequent computation of high-quality motion-free quantitative R 2 ∗ (and any other mGRE-enabled) maps using the standard voxel-wise analysis or machine learning-based analysis. The second CNN, LEARN-BIO, is trained to directly generate motion- and B 0 -inhomogeneity-corrected quantitative R 2 ∗ maps from motion-corrupted magnitude-only mGRE images by taking advantage of the biophysical model describing the mGRE signal decay.

RESULTS

We show that both CNNs trained on synthetic MR images are capable of suppressing motion artifacts while preserving details in the predicted quantitative R 2 ∗ maps. Significant reduction of motion artifacts on experimental in vivo motion-corrupted data has also been achieved by using our trained models.

CONCLUSION

Both LEARN-IMG and LEARN-BIO can enable the computation of high-quality motion- and B 0 -inhomogeneity-corrected R 2 ∗ maps. LEARN-IMG performs motion correction on mGRE images and relies on the subsequent analysis for the estimation of R 2 ∗ maps, while LEARN-BIO directly performs motion- and B 0 -inhomogeneity-corrected R 2 ∗ estimation. Both LEARN-IMG and LEARN-BIO jointly process all the available gradient echoes, which enables them to exploit spatial patterns available in the data. The high computational speed of LEARN-BIO is an advantage that can lead to a broader clinical application.

Penumbra Detection With Oxygen Extraction Fraction Using Magnetic Susceptibility in Patients With Acute Ischemic Stroke

Background

The oxygen extraction fraction (OEF) has been applied to identify ischemic penumbral tissue, but is difficult to use in an urgent care setting. This study aimed to investigate whether an OEF map generated via magnetic resonance quantitative susceptibility mapping (QSM) could help identify the ischemic penumbra in patients with acute ischemic stroke.

Materials and Methods

This prospective imaging study included 21 patients with large anterior circulation vessel occlusion who were admitted <24 h after stroke onset and 21 age-matched healthy controls. We identified the ischemic penumbra as the region with a Tmax of >6 s during dynamic susceptibility contrast-magnetic resonance imaging (DSC-MRI) and calculated the perfusion-core mismatch ratio between the ischemic penumbra and infarct core volumes. The OEF values were measured based on magnetic susceptibility differences between the venous structures and brain tissues using rapid QSM acquisition. Volumes with increased OEF values were compared to the ischemic penumbra volumes using an anatomical template.

Results

Eleven patients had a perfusion-core mismatch ratio of ≥1.8, and reperfusion therapy was recommended. In these patients, the volumes with increased OEF values of >51.5%, which was defined using the anterior circulation territory OEF values from the 21 healthy controls, were positively correlated with the ischemic penumbra volumes (r = 0.636, 95% CI: 0.059 to 0.895, P = 0.035) and inversely correlated with the 30-day change in the National Institutes of Health Stroke Scale scores (r = −0.624, 95% CI: −0.891 to −0.039, P = 0.041).

Conclusion

Tissue volumes with increased OEF values could predict ischemic penumbra volumes based on DSC-MRI, highlighting the potential of the QSM-derived OEF map as a penumbra biomarker to guide treatment selection in patients with acute ischemic stroke.

Brain oxygen extraction fraction mapping in patients with multiple sclerosis

We aimed to demonstrate the feasibility of whole brain oxygen extraction fraction (OEF) mapping for measuring lesion specific and regional OEF abnormalities in multiple sclerosis (MS) patients. In 22 MS patients and 11 healthy controls (HC), OEF and neural tissue susceptibility (χn) maps were computed from MRI multi-echo gradient echo data. In MS patients, 80 chronic active lesions with hyperintense rim on quantitative susceptibility mapping were identified, and the mean OEF and χn within the rim and core were compared using linear mixed-effect model analysis. The rim showed higher OEF and χn than the core: relative to their adjacent normal appearing white matter, OEF contrast = -6.6 ± 7.0% vs. -9.8 ± 7.8% (p < 0.001) and χn contrast = 33.9 ± 20.3 ppb vs. 25.7 ± 20.5 ppb (p = 0.017). Between MS and HC, OEF and χn were compared using a linear regression model in subject-based regions of interest. In the whole brain, compared to HC, MS had lower OEF, 30.4 ± 3.3% vs. 21.4 ± 4.4% (p < 0.001), and higher χn, -23.7 ± 7.0 ppb vs. -11.3 ± 7.7 ppb (p = 0.018). Our feasibility study suggests that OEF may serve as a useful quantitative marker of tissue oxygen utilization in MS.

Whole-brain 3D mapping of oxygen metabolism using constrained quantitative BOLD

Quantitative BOLD (qBOLD) MRI permits noninvasive evaluation of hemodynamic and metabolic states of the brain by quantifying parametric maps of deoxygenated blood volume (DBV) and hemoglobin oxygen saturation level of venous blood (Y_v), and along with a measurement of cerebral blood flow (CBF), the cerebral metabolic rate of oxygen (CMRO₂). The method, thus should have potential to provide important information on many neurological disorders as well as normal cerebral physiology. One major challenge in qBOLD is to separate de-oxyhemoglobin’s contribution to R₂′ from other sources modulating the voxel signal, for instance, R₂, R₂′ from non-heme iron (R′_2,nh), and macroscopic magnetic field variations. Further, even with successful separation of the several confounders, it is still challenging to extract DBV and Y_v from the heme-originated R₂′ because of limited sensitivity of the qBOLD model. These issues, which have not been fully addressed in currently practiced qBOLD methods, have so far precluded 3D whole-brain implementation of qBOLD. Thus, the purpose of this work was to develop a new 3D MRI oximetry technique that enables robust qBOLD parameter mapping across the entire brain. To achieve this goal, we employed a rapid, R₂′-sensitive, steady-state 3D pulse sequence (termed ‘AUSFIDE’) for data acquisition, and implemented a prior-constrained qBOLD processing pipeline that exploits a plurality of preliminary parameters obtained via AUSFIDE, along with additionally measured cerebral venous blood volume. Numerical simulations and in vivo studies at 3 T were performed to evaluate the performance of the proposed, constrained qBOLD mapping in comparison to the parent qBOLD method. Measured parameters (Y_v, DBV, R′_2,nh, nonblood magnetic susceptibility) in ten healthy subjects demonstrate the expected contrast across brain territories, while yielding group-averages of 64.0 ± 2.3 % and 62.2 ± 3.1 % for Y_v and 2.8 ± 0.5 % and 1.8 ± 0.4 % for DBV in cortical gray and white matter, respectively. Given the Y_v measurements, additionally quantified CBF in seven of the ten study subjects enabled whole-brain 3D CMRO₂ mapping, yielding group averages of 134.2 ± 21.1 and 79.4 ± 12.6 µmol/100 g/min for cortical gray and white matter, in good agreement with literature values. The results suggest feasibility of the proposed method as a practical and reliable means for measuring neurometabolic parameters over an extended brain coverage.

QQ-NET - using deep learning to solve quantitative susceptibility mapping and quantitative blood oxygen level dependent magnitude (QSM+qBOLD or QQ) based oxygen extraction fraction (OEF) mapping

PURPOSE

To improve accuracy and speed of quantitative susceptibility mapping plus quantitative blood oxygen level-dependent magnitude (QSM+qBOLD or QQ) -based oxygen extraction fraction (OEF) mapping using a deep neural network (QQ-NET).

METHODS

The 3D multi-echo gradient echo images were acquired in 34 ischemic stroke patients and 4 healthy subjects. Arterial spin labeling and diffusion weighted imaging (DWI) were also performed in the patients. NET was developed to solve the QQ model inversion problem based on Unet. QQ-based OEF maps were reconstructed with previously introduced temporal clustering, tissue composition, and total variation (CCTV) and NET. The results were compared in simulation, ischemic stroke patients, and healthy subjects using a two-sample Kolmogorov-Smirnov test.

RESULTS

In the simulation, QQ-NET provided more accurate and precise OEF maps than QQ-CCTV with 150 times faster reconstruction speed. In the subacute stroke patients, OEF from QQ-NET had greater contrast-to-noise ratio (CNR) between DWI-defined lesions and their unaffected contralateral normal tissue than with QQ-CCTV: 1.9 ± 1.3 vs 6.6 ± 10.7 (p = 0.03). In healthy subjects, both QQ-CCTV and QQ-NET provided uniform OEF maps.

CONCLUSION

QQ-NET improves the accuracy of QQ-based OEF with faster reconstruction.

Application of Cluster Analysis of Time Evolution for Magnetic Resonance Imaging -Derived Oxygen Extraction Fraction Mapping: A Promising Strategy for the Genetic Profile Prediction and Grading of Glioma

Background: The intratumoral heterogeneity of oxygen metabolism and angiogenesis are core hallmarks of glioma, unveiling that genetic aberrations associated with magnetic resonance imaging (MRI) phenotypes may aid in the diagnosis and treatment of glioma.

Objective: To explore the predictability of MRI-based oxygen extraction fraction (OEF) mapping using cluster analysis of time evolution (CAT) for genetic profiling and glioma grading.

Methods: Ninety-one patients with histopathologically confirmed glioma were examined with CAT for quantitative susceptibility mapping and quantitative blood oxygen level–dependent magnitude-based OEF mapping and dynamic contrast-enhanced (DCE) MRI. Imaging biomarkers, including oxygen metabolism (OEF) and angiogenesis [volume transfer constant, cerebral blood volume (CBV), and cerebral blood flow], were investigated to predict IDH mutation, O⁶-methylguanine-DNA-methyltransferase (MGMT) promoter methylation status, receptor tyrosine kinase (RTK) subgroup, and differentiation of glioblastoma (GBM) vs. lower-grade glioma (LGG). The corresponding DNA sequencing was also obtained. Results were compared with DCE-MRI using receiver operating characteristic (ROC) analysis.

Results: IDH1-mutated LGGs exhibited significantly lower OEF and hypoperfusion than IDH wild-type tumors (all p < 0.01). OEF and perfusion metrics showed a tendency toward higher values in MGMT unmethylated GBM, but only OEF retained significance (p = 0.01). Relative prevalence of RTK alterations was associated with increased OEF (p = 0.003) and perfusion values (p < 0.05). ROC analysis suggested OEF achieved best performance for IDH mutation detection [area under the curve (AUC) = 0.828]. None of the investigated parameters enabled prediction of MGMT status except OEF with a moderate AUC of 0.784. Predictive value for RTK subgroup was acceptable by using OEF (AUC = 0.764) and CBV (AUC = 0.754). OEF and perfusion metrics demonstrated excellent performance in glioma grading. Moreover, mutational landscape revealed hypoxia or angiogenesis-relevant gene signatures were associated with specific imaging phenotypes.

Conclusion: CAT for MRI-based OEF mapping is a promising technology for oxygen measurement and along with perfusion MRI can predict genetic profiles and tumor grade in a non-invasive and clinically relevant manner.

Clinical Impact: Physiological imaging provides an in vivo portrait of genetic alterations in glioma and offers a potential strategy for non-invasively selecting patients for individualized therapies.

Impaired cerebral vascular and metabolic responses to parametric N-back tasks in subjective cognitive decline

Previous studies reported abnormally increased and/or decreased blood oxygen level-dependent (BOLD) activations during functional tasks in subjective cognitive decline (SCD). The neurophysiological basis underlying these functional aberrations remains debated. This study aims to investigate vascular and metabolic responses and their dependence on cognitive processing loads during functional tasks in SCD. Twenty-one SCD and 18 control subjects performed parametric N-back working-memory tasks during MRI scans. Task-evoked percentage changes (denoted as δ) in cerebral blood volume (δCBV), cerebral blood flow (δCBF), BOLD signal (δBOLD) and cerebral metabolic rate of oxygen (δCMRO₂) were evaluated. In the frontal lobe, trends of decreased δCBV, δCBF and δCMRO₂ and increased δBOLD were observed in SCD compared with control subjects under lower loads, and these trends increased to significant differences under the 3-back load. δCBF was significantly correlated with δCMRO₂ in controls, but not in SCD subjects. As N-back loads increased, the differences between SCD and control subjects in δCBF and δCMRO₂ tended to enlarge. In the parietal lobe, no significant between-group difference was observed. Our findings suggested that impaired vascular and metabolic responses to functional tasks occurred in the frontal lobe of SCD, which contributed to unusual BOLD hyperactivation and was modulated by cognitive processing loads.

The Spatiotemporal Evolution of MRI-Derived Oxygen Extraction Fraction and Perfusion in Ischemic Stroke

Purpose

This study aimed to assess the spatiotemporal evolution of oxygen extraction fraction (OEF) in ischemic stroke with a newly developed cluster analysis of time evolution (CAT) for a combined quantitative susceptibility mapping and quantitative blood oxygen level-dependent model (QSM + qBOLD, QQ).

Method

One hundred and fifteen patients in different ischemic stroke phases were retrospectively collected for measurement of OEF of the infarcted area defined on diffusion-weighted imaging (DWI). Clinical severity was assessed using the National Institutes of Health Stroke Scale (NIHSS). Of the 115 patients, 11 underwent two longitudinal MRI scans, namely, three-dimensional (3D) multi-echo gradient recalled echo (mGRE) and 3D pseudo-continuous arterial spin labeling (pCASL), to evaluate the reversal region (RR) of the initial diffusion lesion (IDL) that did not overlap with the final infarct (FI). The temporal evolution of OEF and the cerebral blood flow (CBF) in the IDL, the RR, and the FI were assessed.

Results

Compared to the contralateral mirror area, the OEF of the infarcted region was decreased regardless of stroke phases (p < 0.05) and showed a declining tendency from the acute to the chronic phase (p = 0.022). Five of the 11 patients with longitudinal scans showed reversal of the IDL. Relative oxygen extraction fraction (rOEF, compared to the contralateral mirror area) of the RR increased from the first to the second MRI (p = 0.044). CBF was about 1.5-fold higher in the IDL than in the contralateral mirror area in the first MRI. Two patients showed penumbra according to the enlarged FI volume. The rOEF of the penumbra fluctuated around 1.0 at earlier scan times and then decreased, while the CBF decreased continuously.

Conclusion

The spatiotemporal evolution of OEF and perfusion in ischemic lesions is heterogeneous, and the CAT-based QQ method is feasible to capture cerebral oxygen metabolic information.

Brain metabolism in tau and amyloid mouse models of Alzheimer's disease: An MRI study

Alzheimer's disease (AD) is the leading cause of cognitive impairment and dementia in elderly individuals. According to the current biomarker framework for "unbiased descriptive classification", biomarkers of neurodegeneration, "N", constitute a critical component in the tri-category "A/T/N" system. Current biomarkers of neurodegeneration suffer from potential drawbacks such as requiring invasive lumbar puncture, involving ionizing radiation, or representing a late, irreversible marker. Recent human studies have suggested that reduced brain oxygen metabolism may be a new functional marker of neurodegeneration in AD, but the heterogeneity and the presence of mixed pathology in human patients did not allow a full understanding of the role of oxygen extraction and metabolism in AD. In this report, global brain oxygen metabolism and related physiological parameters were studied in two AD mouse models with relatively pure pathology, using advanced MRI techniques including T2 -relaxation-under-spin-tagging (TRUST) and phase contrast (PC) MRI. Additionally, regional cerebral blood flow (CBF) was determined with pseudocontinuous arterial spin labeling. Reduced global oxygen extraction fraction (by -18.7%, p = 0.008), unit-mass cerebral metabolic rate of oxygen (CMRO2 ) (by -17.4%, p = 0.04) and total CMRO2 (by -30.8%, p < 0.001) were observed in Tau4RΔK mice-referred to as the tau AD model-which manifested pronounced neurodegeneration, as measured by diminished brain volume (by -15.2%, p < 0.001). Global and regional CBF in these mice were not different from those of wild-type mice (p > 0.05), suggesting normal vascular function. By contrast, in B6;SJL-Tg [APPSWE]2576Kha (APP) mice-referred to as the amyloid AD model-no brain volume reduction, as well as relatively intact brain oxygen extraction and metabolism, were found (p > 0.05). Consistent with the imaging data, behavioral measures of walking distance were impaired in Tau4RΔK mice (p = 0.004), but not in APP mice (p = 0.88). Collectively, these findings support the hypothesis that noninvasive MRI measurement of brain oxygen metabolism may be a promising biomarker of neurodegeneration in AD.

Cerebral oxygen extraction fraction (OEF): Comparison of challenge-free gradient echo QSM+qBOLD (QQ) with 15O PET in healthy adults

We aimed to validate oxygen extraction fraction (OEF) estimations by quantitative susceptibility mapping plus quantitative blood oxygen-level dependence (QSM+qBOLD, or QQ) using ¹⁵O-PET. In ten healthy adult brains, PET and MRI were acquired simultaneously on a PET/MR scanner. PET was acquired using C[¹⁵O], O[¹⁵O], and H₂[¹⁵O]. Image-derived arterial input functions and standard models of oxygen metabolism provided quantification of PET. MRI included T1-weighted imaging, time-of-flight angiography, and multi-echo gradient-echo imaging that was processed for QQ. Region of interest (ROI) analyses compared PET OEF and QQ OEF. In ROI analyses, the averaged OEF differences between PET and QQ were generally small and statistically insignificant. For whole brains, the average and standard deviation of OEF was 32.8 ± 6.7% for PET; OEF was 34.2 ± 2.6% for QQ. Bland-Altman plots quantified agreement between PET OEF and QQ OEF. The interval between the 95% limits of agreement was 16.9 ± 4.0% for whole brains. Our validation study suggests that respiratory challenge-free QQ-OEF mapping may be useful for non-invasive clinical assessment of regional OEF impairment.

Temporal clustering, tissue composition, and total variation for mapping oxygen extraction fraction using QSM and quantitative BOLD

PURPOSE

To improve the accuracy of quantitative susceptibility mapping plus quantitative blood oxygen level-dependent magnitude (QSM+qBOLD or QQ) based mapping of oxygen extraction fraction (OEF) and cerebral metabolic rate of oxygen (CMRO₂ ) using temporal clustering, tissue composition, and total variation (CCTV).

METHODS

Three-dimensional multi-echo gradient echo and arterial spin labeling images were acquired from 11 healthy subjects and 33 ischemic stroke patients. Diffusion-weighted imaging (DWI) was also obtained from patients. The CCTV mapping was developed for incorporating tissue-type information into clustering of the previous cluster analysis of time evolution (CAT) and applying total variation (TV). The QQ-based OEF and CMRO₂ were reconstructed with CAT, CAT+TV (CATV), and the proposed CCTV, and results were compared using region-of-interest analysis, Kruskal-Wallis test, and post hoc Wilcoxson rank sum test.

RESULTS

In simulation, CCTV provided more accurate and precise OEF than CAT or CATV. In healthy subjects, QQ-based OEF was less noisy and more uniform with CCTV than CAT. In subacute stroke patients, OEF with CCTV had a greater contrast-to-noise ratio between DWI-defined lesions and the unaffected contralateral side than with CAT or CATV: 1.9 ± 1.3 versus 1.1 ± 0.7 (P = .01) versus 0.7 ± 0.5 (P < .001).

CONCLUSION

The CCTV mapping significantly improves the robustness of QQ-based OEF against noise.

Multiecho complex total field inversion method (mcTFI) for improved signal modeling in quantitative susceptibility mapping

PURPOSE

Typical quantitative susceptibility mapping (QSM) reconstruction steps consist of first estimating the magnetization field from the gradient-echo images, and then reconstructing the susceptibility map from the estimated field. The errors from the field-estimation steps may propagate into the final QSM map, and the noise in the estimated field map may no longer be zero-mean Gaussian noise, thus, causing streaking artifacts in the resulting QSM. A multiecho complex total field inversion (mcTFI) method was developed to compute the susceptibility map directly from the multiecho gradient echo images using an improved signal model that retains the Gaussian noise property in the complex domain. It showed improvements in QSM reconstruction over the conventional field-to-source inversion.

METHODS

The proposed mcTFI method was compared with the nonlinear total field inversion (nTFI) method in a numerical brain with hemorrhage and calcification, the numerical brains provided by the QSM Challenge 2.0, 18 brains with intracerebral hemorrhage scanned at 3T, and 6 healthy brains scanned at 7T.

RESULTS

Compared with nTFI, the proposed mcTFI showed more accurate QSM reconstruction around the lesions in the numerical simulations. The mcTFI reconstructed QSM also showed the best image quality with the least artifacts in the brains with intracerebral hemorrhage scanned at 3T and healthy brains scanned at 7T.

CONCLUSION

The proposed multiecho complex total field inversion improved QSM reconstruction over traditional field-to-source inversion through better signal modeling.

Non-Invasive Evaluation of Cerebral Microvasculature Using Pre-Clinical MRI: Principles, Advantages and Limitations

Alterations to the cerebral microcirculation have been recognized to play a crucial role in the development of neurodegenerative disorders. However, the exact role of the microvascular alterations in the pathophysiological mechanisms often remains poorly understood. The early detection of changes in microcirculation and cerebral blood flow (CBF) can be used to get a better understanding of underlying disease mechanisms. This could be an important step towards the development of new treatment approaches. Animal models allow for the study of the disease mechanism at several stages of development, before the onset of clinical symptoms, and the verification with invasive imaging techniques. Specifically, pre-clinical magnetic resonance imaging (MRI) is an important tool for the development and validation of MRI sequences under clinically relevant conditions. This article reviews MRI strategies providing indirect non-invasive measurements of microvascular changes in the rodent brain that can be used for early detection and characterization of neurodegenerative disorders. The perfusion MRI techniques: Dynamic Contrast Enhanced (DCE), Dynamic Susceptibility Contrast Enhanced (DSC) and Arterial Spin Labeling (ASL), will be discussed, followed by less established imaging strategies used to analyze the cerebral microcirculation: Intravoxel Incoherent Motion (IVIM), Vascular Space Occupancy (VASO), Steady-State Susceptibility Contrast (SSC), Vessel size imaging, SAGE-based DSC, Phase Contrast Flow (PC) Quantitative Susceptibility Mapping (QSM) and quantitative Blood-Oxygenation-Level-Dependent (qBOLD). We will emphasize the advantages and limitations of each strategy, in particular on applications for high-field MRI in the rodent's brain.

Non-Invasive measurement of the cerebral metabolic rate of oxygen using MRI in rodents

Malfunctions of oxygen metabolism are suspected to play a key role in a number of neurological and psychiatric disorders, but this hypothesis cannot be properly investigated without an in-vivo non-invasive measurement of brain oxygen consumption. We present a new way to measure the Cerebral Metabolic Rate of Oxygen (CMRO2) by combining two existing magnetic resonance imaging techniques, namely arterial spin-labelling and oxygen extraction fraction mapping. This method was validated by imaging rats under different anaesthetic regimes and was strongly correlated to glucose consumption measured by autoradiography.

Prediction of an oxygen extraction fraction map by convolutional neural network: validation of input data among MR and PET images

Purpose

Oxygen extraction fraction (OEF) is a biomarker for the viability of brain tissue in ischemic stroke. However, acquisition of the OEF map using positron emission tomography (PET) with oxygen-15 gas is uncomfortable for patients because of the long fixation time, invasive arterial sampling, and radiation exposure. We aimed to predict the OEF map from magnetic resonance (MR) and PET images using a deep convolutional neural network (CNN) and to demonstrate which PET and MR images are optimal as inputs for the prediction of OEF maps.

Methods

Cerebral blood flow at rest (CBF) and during stress (sCBF), cerebral blood volume (CBV) maps acquired from oxygen-15 PET, and routine MR images (T1-, T2-, and T2*-weighted images) for 113 patients with steno-occlusive disease were learned with U-Net. MR and PET images acquired from the other 25 patients were used as test data. We compared the predicted OEF maps and intraclass correlation (ICC) with the real OEF values among combinations of MRI, CBF, CBV, and sCBF.

Results

Among the combinations of input images, OEF maps predicted by the model learned with MRI, CBF, CBV, and sCBF maps were the most similar to the real OEF maps (ICC: 0.597 ± 0.082). However, the contrast of predicted OEF maps was lower than that of real OEF maps.

Conclusion

These results suggest that the deep CNN learned useful features from CBF, sCBF, CBV, and MR images and predict qualitatively realistic OEF maps. These findings suggest that the deep CNN model can shorten the fixation time for ¹⁵O PET by skipping ¹⁵O₂ scans. Further training with a larger data set is required to predict accurate OEF maps quantitatively.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11548-021-02356-7.

Venous cerebral blood volume mapping in the whole brain using venous-spin-labeled 3D turbo spin echo

PURPOSE

Venous cerebral blood volume (CBVv ) is a major contributor to BOLD contrast, and therefore is an important parameter for understanding the underlying mechanism. Here, we propose a velocity-selective venous spin labeling (VS-VSL)-prepared 3D turbo spin echo pulse sequence for whole-brain baseline CBVv mapping.

METHODS

Unlike previous CBVv measurement techniques that exploit the interrelationship between BOLD signals and CBVv , in the proposed VS-VSL technique both arterial blood and cerebrospinal fluid (CSF) signals were suppressed before the VS pulse train for exclusive labeling of venous blood, while a single-slab 3D turbo spin echo readout was used because of its relative immunity to magnetic field variations. Furthermore, two approximations were made to the VS-VSL signal model for simplified derivation of CBVv . In vivo studies were performed at 3T field strength in 8 healthy subjects. The performance of the proposed VS-VSL method in baseline CBVv estimation was first evaluated in comparison to the existing, hyperoxia-based method. Then, data were also acquired using VS-VSL under hypercapnic and hyperoxic gas breathing challenges for further validation of the technique.

RESULTS

The proposed technique yielded physiologically plausible baseline CBVv values, and when compared with the hyperoxia-based method, showed no statistical difference. Furthermore, data acquired using VS-VSL yielded average CBVv of 2.89%/1.78%, 3.71%/2.29%, and 2.88%/1.76% for baseline, hypercapnia, and hyperoxia, respectively, in gray/white matter regions. As expected, hyperoxia had negligible effect (P > .8), whereas hypercapnia demonstrated vasodilation (P << .01).

CONCLUSION

Upon further validation of the quantification model, the method is expected to have merit for 3D CBVv measurements across the entire brain.

Initial Experience of Challenge-Free MRI-Based Oxygen Extraction Fraction Mapping of Ischemic Stroke at Various Stages: Comparison With Perfusion and Diffusion Mapping

MRI-based oxygen extraction fraction imaging has a great potential benefit in the selection of clinical strategies for ischemic stroke patients. This study aimed to evaluate the performance of a challenge-free oxygen extraction fraction (OEF) mapping in a cohort of acute and subacute ischemic stroke patients. Consecutive ischemic stroke patients (a total of 30 with 5 in the acute stage, 19 in the early subacute stage, and 6 in the late subacute stage) were recruited. All subjects underwent MRI including multi-echo gradient echo (mGRE), diffusion weighted imaging (DWI), and 3D-arterial spin labeling (ASL). OEF maps were generated from mGRE phase + magnitude data, which were processed using quantitative susceptibility mapping (QSM) + quantitative blood oxygen level-dependent (qBOLD) imaging with cluster analysis of time evolution. Cerebral blood flow (CBF) and apparent diffusion coefficient (ADC) maps were reconstructed from 3D-ASL and DWI, respectively. Further, cerebral metabolic rate of oxygen (CMRO2) was calculated as the product of CBF and OEF. OEF, CMRO2, CBF, and ADC values in the ischemic cores (absolute values) and their contrasts to the contralateral regions (relative values) were evaluated. One-way analysis of variance (ANOVA) was used to compare OEF, CMRO2, CBF, and ADC values and their relative values among different stroke stages. The OEF value of infarct core showed a trend of decrease from acute, to early subacute, and to late subacute stages of ischemic stroke. Significant differences among the three stroke stages were only observed in the absolute OEF (F = 6.046, p = 0.005) and relative OEF (F = 5.699, p = 0.009) values of the ischemic core, but not in other measurements (absolute and relative CMRO2, CBF, ADC values, all values of p > 0.05). In conclusion, the challenge-free QSM + qBOLD-generated OEF mapping can be performed on stroke patients. It can provide more information on tissue viability that was not available with CBF and ADC and, thus, may help to better manage ischemic stroke patients.

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSION

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

Elevated brain oxygen extraction fraction measured by MRI susceptibility relates to perfusion status in acute ischemic stroke

Recent clinical trials of new revascularization therapies in acute ischemic stroke have highlighted the importance of physiological imaging to identify optimal treatments for patients. Oxygen extraction fraction (OEF) is a hallmark of at-risk tissue in stroke, and can be quantified from the susceptibility effect of deoxyhemoglobin molecules in venous blood on MRI phase scans. We measured OEF within cerebral veins using advanced quantitative susceptibility mapping (QSM) MRI reconstructions in 20 acute stroke patients. Absolute OEF was elevated in the affected (29.3 ± 3.4%) versus the contralateral hemisphere (25.5 ± 3.1%) of patients with large diffusion-perfusion lesion mismatch (P = 0.032). In these patients, OEF negatively correlated with relative CBF measured by dynamic susceptibility contrast MRI (P = 0.004), suggesting compensation for reduced flow. Patients with perfusion-diffusion match or no hypo-perfusion showed less OEF difference between hemispheres. Nine patients received longitudinal assessment and showed OEF ratio (affected to contralateral) of 1.2 ± 0.1 at baseline that normalized (decreased) to 1.0 ± 0.1 at follow-up three days later (P = 0.03). Our feasibility study demonstrates that QSM MRI can non-invasively quantify OEF in stroke patients, relates to perfusion status, and is sensitive to OEF changes over time. Clinical trial registration: Longitudinal MRI examinations of patients with brain ischemia and blood brain barrier permeability; clinicaltrials.org :NCT02077582.

Cluster analysis of time evolution (CAT) for quantitative susceptibility mapping (QSM) and quantitative blood oxygen level-dependent magnitude (qBOLD)-based oxygen extraction fraction (OEF) and cerebral metabolic rate of oxygen (CMRO2 ) mapping

PURPOSE

To improve the accuracy of QSM plus quantitative blood oxygen level-dependent magnitude (QSM + qBOLD or QQ)-based mapping of the oxygen extraction fraction (OEF) and cerebral metabolic rate of oxygen (CMRO2 ) using cluster analysis of time evolution (CAT).

METHODS

3D multi-echo gradient echo and arterial spin labeling images were acquired in 11 healthy subjects and 5 ischemic stroke patients. DWI was also carried out on patients. CAT was developed for analyzing signal evolution over TE. QQ-based OEF and CMRO2 were reconstructed with and without CAT, and results were compared using region of interest analysis and a paired t-test.

RESULTS

Simulations demonstrated that CAT substantially reduced noise error in QQ-based OEF. In healthy subjects, QQ-based OEF appeared less noisy and more uniform with CAT than without CAT; average OEF with and without CAT in cortical gray matter was 32.7 ± 4.0% and 37.9 ± 4.5%, with corresponding CMRO2 of 148.4 ± 23.8 and 171.4 ± 22.4 μmol/100 g/min, respectively. In patients, regions of low OEF were confined within the ischemic lesions defined on DWI when using CAT, which was not observed without CAT.

CONCLUSION

The cluster analysis of time evolution (CAT) significantly improves the robustness of QQ-based OEF against noise.

Quantitative Dynamic Oxygen 17 MRI at 7.0 T for the Cerebral Oxygen Metabolism in Glioma

Background Altered metabolism is a characteristic of cancer. Because of a shift in glucose metabolism from oxidative phosphorylation to lactate production for energy generation, malignant tumors are characterized by increased glycolysis followed by lactic acid fermentation, even in the presence of abundant oxygen (the Warburg effect). Purpose To quantitatively investigate dynamic oxygen 17 (¹⁷O) MRI in healthy participants and participants with untreated glioma to understand altered cerebral oxygen metabolism in glioma. Materials and Methods In this prospective study conducted from September 2016 to June 2018, individuals with newly diagnosed previously untreated glioma (World Health Organization grade II-IV) and healthy volunteers were included. Dynamic ¹⁷O MRI was performed with a 7.0-T whole-body system. ¹⁷O₂ gas inhalation enabled dynamic measurement of the cerebral metabolic rate of oxygen (CMRO₂) consumption. In healthy volunteers and participants with glioma, CMRO₂ values in gray matter and white matter volumes were compared by using Wilcoxon signed rank tests. In participants with glioma, the tumor volume and tumor subcompartments were compared with normal-appearing gray matter and white matter by using Friedman test followed by Holm-Sidak post hoc tests. Results Ten participants (mean age, 42 years ± 18 [standard deviation]; nine men) with glioma and three healthy volunteers (mean age, 44 years ± 21; all men) were evaluated. CMRO₂ was higher in normal-appearing gray matter compared with white matter in both participants with glioma (2.36 μmol/g/min ± 0.22 vs 0.75 μmol/g/min ± 0.10, respectively) and healthy volunteers (2.38 μmol/g/min ± 0.15 vs 0.63 μmol/g/min ± 0.05, respectively) (P < .001 and P = .03, respectively). In the tumor region, CMRO₂ was reduced (high-grade tumor CMRO₂, 0.23 μmol/g/min ± 0.07; low-grade tumor CMRO₂, 0.39 μmol/g/min ± 0.16; overall CMRO₂, 0.34 μmol/g/min ± 0.16) compared with normal-appearing gray matter (P < .001) and normal-appearing white matter (P < .001) in accordance with the Warburg theorem. Conclusion Dynamic oxygen 17 MRI method at 7.0 T as a direct metabolic imaging technique in glioma enabled quantitative visualization of the Warburg effect. A general reduction in oxidative glycolysis was observed in accordance with the Warburg theorem. © RSNA, 2020 Online supplemental material is available for this article. See also the editorial by Rapalino in this issue.