Measuring brain oxygenation in humans using a multiparametric quantitative blood oxygenation level dependent MRI approach

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.

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

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.

Advances in PET and MRI imaging of tumor hypoxia

Tumor hypoxia is a complex and evolving phenomenon both in time and space. Molecular imaging allows to approach these variations, but the tracers used have their own limitations. PET imaging has the disadvantage of low resolution and must take into account molecular biodistribution, but has the advantage of high targeting accuracy. The relationship between the signal in MRI imaging and oxygen is complex but hopefully it would lead to the detection of truly oxygen-depleted tissue. Different ways of imaging hypoxia are discussed in this review, with nuclear medicine tracers such as [¹⁸F]-FMISO, [¹⁸F]-FAZA, or [⁶⁴Cu]-ATSM but also with MRI techniques such as perfusion imaging, diffusion MRI or oxygen-enhanced MRI. Hypoxia is a pejorative factor regarding aggressiveness, tumor dissemination and resistance to treatments. Therefore, having accurate tools is particularly important.

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.

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.

Radiophysiomics: Brain Tumors Classification by Machine Learning and Physiological MRI Data

Simple Summary

The pretreatment diagnosis of contrast-enhancing brain tumors is still challenging in clinical neuro-oncology due to their very similar appearance on conventional MRI. A precise initial characterization, however, is essential to initiate appropriate treatment management, which can substantially differ between brain tumor entities. To overcome the disadvantage of the low specificity of conventional MRI, several new neuroimaging methods have been developed and validated over the past decades. This increasing amount of diagnostic information makes a timely evaluation without computational support impossible in a clinical setting. Artificial intelligence methods such as machine learning offer new options to support clinicians. In this study, we combined nine common machine learning algorithms with a physiological MRI technique (we named this approach “radiophysiomics”) to investigate the effectiveness of the multiclass classification of contrast-enhancing brain tumors in a clinical setting. We were able to demonstrate that radiophysiomics could be helpful in the routine diagnostics of contrast-enhancing brain tumors, but further automation using deep neural networks is required.

Abstract

The precise initial characterization of contrast-enhancing brain tumors has significant consequences for clinical outcomes. Various novel neuroimaging methods have been developed to increase the specificity of conventional magnetic resonance imaging (cMRI) but also the increased complexity of data analysis. Artificial intelligence offers new options to manage this challenge in clinical settings. Here, we investigated whether multiclass machine learning (ML) algorithms applied to a high-dimensional panel of radiomic features from advanced MRI (advMRI) and physiological MRI (phyMRI; thus, radiophysiomics) could reliably classify contrast-enhancing brain tumors. The recently developed phyMRI technique enables the quantitative assessment of microvascular architecture, neovascularization, oxygen metabolism, and tissue hypoxia. A training cohort of 167 patients suffering from one of the five most common brain tumor entities (glioblastoma, anaplastic glioma, meningioma, primary CNS lymphoma, or brain metastasis), combined with nine common ML algorithms, was used to develop overall 135 classifiers. Multiclass classification performance was investigated using tenfold cross-validation and an independent test cohort. Adaptive boosting and random forest in combination with advMRI and phyMRI data were superior to human reading in accuracy (0.875 vs. 0.850), precision (0.862 vs. 0.798), F-score (0.774 vs. 0.740), AUROC (0.886 vs. 0.813), and classification error (5 vs. 6). The radiologists, however, showed a higher sensitivity (0.767 vs. 0.750) and specificity (0.925 vs. 0.902). We demonstrated that ML-based radiophysiomics could be helpful in the clinical routine diagnosis of contrast-enhancing brain tumors; however, a high expenditure of time and work for data preprocessing requires the inclusion of deep neural networks.

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.

Association between tissue hypoxia, perfusion restrictions, and microvascular architecture alterations with lesion-induced impairment of neurovascular coupling

Functional magnetic resonance imaging (fMRI) has been mainly utilized for the preoperative localization of eloquent cortical areas. However, lesion-induced impairment of neurovascular coupling (NVC) in the lesion border zone may lead to false-negative fMRI results. The purpose of this study was to determine physiological factors impacting the NVC. Twenty patients suffering from brain lesions were preoperatively examined using multimodal neuroimaging including fMRI, magnetoencephalography (MEG) during language or sensorimotor tasks (depending on lesion location), and a novel physiologic MRI approach for the combined quantification of oxygen metabolism, perfusion state, and microvascular architecture. Congruence of brain activity patterns between fMRI and MEG were found in 13 patients. In contrast, we observed missing fMRI activity in perilesional cortex that demonstrated MEG activity in seven patients, which was interpreted as lesion-induced impairment of NVC. In these brain regions with impaired NVC, physiologic MRI revealed significant brain tissue hypoxia, as well as significantly decreased macro- and microvascular perfusion and microvascular architecture. We demonstrated that perilesional hypoxia with reduced vascular perfusion and architecture is associated with lesion-induced impairment of NVC. Our physiologic MRI approach is a clinically applicable method for preoperative risk assessment for the presence of false-negative fMRI results and may prevent severe postoperative functional deficits.

Imaging effective oxygen diffusivity in the human brain with multiparametric magnetic resonance imaging

Cerebrovascular diseases can impair blood circulation and oxygen extraction from the blood. The effective oxygen diffusivity (EOD) of the capillary bed is a potential biomarker of microvascular function that has gained increasing interest, both for clinical diagnosis and for elucidating oxygen transport mechanisms. Models of capillary oxygen transport link EOD to measurable oxygen extraction fraction (OEF) and cerebral blood flow (CBF). In this work, we confirm that two well established mathematical models of oxygen transport yield nearly equivalent EOD maps. Furthermore, we propose an easy-to-implement and clinically applicable multiparametric magnetic resonance imaging (MRI) protocol for quantitative EOD mapping. Our approach is based on imaging OEF and CBF with multiparametric quantitative blood oxygenation level dependent (mq-BOLD) MRI and pseudo-continuous arterial spin labeling (pCASL), respectively. We evaluated the imaging protocol by comparing MRI-EOD maps of 12 young healthy volunteers to PET data from a published study in different individuals. Our results show comparably good correlation between MRI- and PET-derived cortical EOD, OEF and CBF. Importantly, absolute values of MRI and PET showed high accordance for all three parameters. In conclusion, our data indicates feasibility of the proposed MRI protocol for EOD mapping, rendering the method promising for future clinical evaluation of patients with cerebrovascular diseases.

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.

Metabolic Tumor Microenvironment Characterization of Contrast Enhancing Brain Tumors Using Physiologic MRI

The tumor microenvironment is a critical regulator of cancer development and progression as well as treatment response and resistance in brain neoplasms. The available techniques for investigation, however, are not well suited for noninvasive in vivo characterization in humans. A total of 120 patients (59 females; 61 males) with newly diagnosed contrast-enhancing brain tumors (64 glioblastoma, 20 brain metastases, 15 primary central nervous system (CNS) lymphomas (PCNSLs), and 21 meningiomas) were examined with a previously established physiological MRI protocol including quantitative blood-oxygen-level-dependent imaging and vascular architecture mapping. Six MRI biomarker maps for oxygen metabolism and neovascularization were fused for classification of five different tumor microenvironments: glycolysis, oxidative phosphorylation (OxPhos), hypoxia with/without neovascularization, and necrosis. Glioblastoma showed the highest metabolic heterogeneity followed by brain metastasis with a glycolysis-to-OxPhos ratio of approximately 2:1 in both tumor entities. In addition, glioblastoma revealed a significant higher percentage of hypoxia (24%) compared to all three other brain tumor entities: brain metastasis (7%; p < 0.001), PCNSL (8%; p = 0.001), and meningioma (8%; p = 0.003). A more aggressive biological brain tumor behavior was associated with a higher percentage of hypoxia and necrosis and a lower percentage of remaining vital tumor tissue and aerobic glycolysis. The proportion of oxidative phosphorylation, however, was rather similar (17–26%) for all four brain tumor entities. Tumor microenvironment (TME) mapping provides insights into neurobiological differences of contrast-enhancing brain tumors and deserves further clinical cancer research attention. Although there is a long roadmap ahead, TME mapping may become useful in order to develop new diagnostic and therapeutic approaches.

Physiological MRI of microvascular architecture, neovascularization activity, and oxygen metabolism facilitate early recurrence detection in patients with IDH-mutant WHO grade 3 glioma

PURPOSE

This study aimed to determine the diagnostic performance of physiological MRI biomarkers including microvascular perfusion and architecture, neovascularization activity, tissue oxygen metabolism, and tension for recurrence detection of IDH-mutant WHO grade 3 glioma.

METHODS

Sixty patients with IDH-mutant WHO grade 3 glioma who received overall 288 follow-up MRI examinations at 3 Tesla after standard treatment were retrospectively evaluated. A conventional MRI protocol was extended with a physiological MRI approach including vascular architecture mapping and quantitative blood-oxygen-level-dependent imaging which required 7 min extra data acquisition time. Custom-made MATLAB software was used for the calculation of MRI biomarker maps of microvascular perfusion and architecture, neovascularization activity, tissue oxygen metabolism, and tension. Statistical procedures included receiver operating characteristic analysis.

RESULTS

Overall, 34 patients showed recurrence of the WHO grade 3 glioma; of these, in 15 patients, recurrence was detected one follow-up examination (averaged 160 days) earlier by physiological MRI data than by conventional MRI. During this time period, the tumor volume increased significantly (P = 0.001) on average 7.4-fold from 1.5 to 11.1 cm³. Quantitative analysis of MRI biomarkers demonstrated microvascular but no macrovascular hyperperfusion in early recurrence. Neovascularization activity (AUC = 0.833), microvascular perfusion (0.682), and oxygen metabolism (0.661) showed higher diagnostic performance for early recurrence detection of WHO grade 3 glioma compared to conventional MRI including cerebral blood volume (0.649).

CONCLUSION

This study demonstrated that the targeted assessment of microvascular features and tissue oxygen tension as an early sign of neovascularization activity provided valuable information for recurrence diagnostic of WHO grade 3 glioma.

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.

Hypoxia and Microvascular Alterations Are Early Predictors of IDH-Mutated Anaplastic Glioma Recurrence

Simple Summary

Anaplastic gliomas (AGs) are considered the most common and aggressive primary brain tumors of young adults with inevitable recurrence and treatment failure. The aim of this study was to investigate whether the imaging biomarkers of hypoxia, microvascular architecture and neovascularization activity can be of assistance to detect pathophysiological changes in the early developmental stages of isocitrate-dehydrogenase (IDH) mutated AG recurrence. We evaluated 142 physiological magnetic resonance imaging follow-up examinations as a part of the conventional magnetic resonance imaging (MRI) protocol in 60 AG patients after standard therapy. Physiological MRI biomarkers showed intensifying local tissue hypoxia 250 days prior to radiological recurrence with following upregulation of neovascularization activity 50 to 70 days later. Integration of physiological MRI in the monitoring of AG patients may be of clinical significance to make personalized decision of early tumor recurrence without an additional delay for multimodal therapy.

Abstract

Anaplastic gliomas (AG) represents aggressive brain tumors that often affect young adults. Although isocitrate-dehydrogenase (IDH) gene mutation has been identified as a more favorable prognostic factor, most IDH-mutated AG patients are confronted with tumor recurrence. Hence, increased knowledge about pathophysiological precursors of AG recurrence is urgently needed in order to develop precise diagnostic monitoring and tailored therapeutic approaches. In this study, 142 physiological magnetic resonance imaging (phyMRI) follow-up examinations in 60 AG patients after standard therapy were evaluated and magnetic resonance imaging (MRI) biomarker maps for microvascular architecture and perfusion, neovascularization activity, oxygen metabolism, and hypoxia calculated. From these 60 patients, 34 patients developed recurrence of the AG, and 26 patients showed no signs for AG recurrence during the study period. The time courses of MRI biomarker changes were analyzed regarding early pathophysiological alterations over a one-year period before radiological AG recurrence or a one-year period of stable disease for patients without recurrence, respectively. We detected intensifying local tissue hypoxia 250 days prior to radiological recurrence which initiated upregulation of neovascularization activity 50 to 70 days later. These changes were associated with a switch from an avascular infiltrative to a vascularized proliferative phenotype of the tumor cells another 30 days later. The dynamic changes of blood perfusion, microvessel density, neovascularization activity, and oxygen metabolism showed a close physiological interplay in the one-year period prior to radiological recurrence of IDH-mutated AG. These findings may path the wave for implementing both new MR-based imaging modalities for routine follow-up monitoring of AG patients after standard therapy and furthermore may support the development of novel, tailored therapy options in recurrent AG.

Advanced Imaging Techniques for Radiotherapy Planning of Gliomas

The accuracy of target delineation in radiation treatment (RT) planning of cerebral gliomas is crucial to achieve high tumor control, while minimizing treatment-related toxicity. Conventional magnetic resonance imaging (MRI), including contrast-enhanced T1-weighted and fluid-attenuated inversion recovery (FLAIR) sequences, represents the current standard imaging modality for target volume delineation of gliomas. However, conventional sequences have limited capability to discriminate treatment-related changes from viable tumors, owing to the low specificity of increased blood-brain barrier permeability and peritumoral edema. Advanced physiology-based MRI techniques, such as MR spectroscopy, diffusion MRI and perfusion MRI, have been developed for the biological characterization of gliomas and may circumvent these limitations, providing additional metabolic, structural, and hemodynamic information for treatment planning and monitoring. Radionuclide imaging techniques, such as positron emission tomography (PET) with amino acid radiopharmaceuticals, are also increasingly used in the workup of primary brain tumors, and their integration in RT planning is being evaluated in specialized centers. This review focuses on the basic principles and clinical results of advanced MRI and PET imaging techniques that have promise as a complement to RT planning of gliomas.

Alternating unbalanced SSFP for 3D R 2 ' mapping of the human brain

PURPOSE

Measuring the transverse-relaxation rate R 2 ' provides valuable information in quantitative evaluation of tissue microstructure, for example, in terms of oxygenation levels. Here, we propose an alternating unbalanced SSFP pulse sequence for rapid whole-brain 3D R 2 ' mapping.

METHODS

Unlike currently practiced, spin echo-based R 2 ' measurement techniques, the proposed method alternates between SSFP-FID and SSFP-ECHO modes for rapid 3D encoding of transverse relaxation rates expressed as R₂ + R 2 ' and R₂ - R 2 ' . Z-shimming gradients embedded into multi-echo trains of each SSFP module are designed to achieve relative immunity to large-scale magnetic-field variations (ΔB₀ ). Appropriate models for the temporal evolution of the two groups of SSFP signals were constructed with ΔB₀ -induced modulations accounted for, leading to ΔB₀ -corrected estimation of R₂ , R 2 ' , and R 2 ∗ (= R₂ + R 2 ' ). Additionally, relative magnetic susceptibility (Δχ) maps were obtained by quantitative susceptibility mapping of the phase data. Numerical simulations were performed to optimize scan parameters, followed by in vivo studies at 3 T in 7 healthy subjects. Measured parameters were evaluated in six brain regions, and subjected to interparameter correlation analysis.

RESULTS

The resultant maps of R 2 ' and additionally derived R₂ , R 2 ∗ , and Δχ all demonstrated the expected contrast across brain territories (eg, deep brain structures versus cortex), with the measured values in good agreement with previous reports. Furthermore, regression analyses yielded strong linear relationships for the transverse relaxation parameters ( R 2 ' , R₂ , and R 2 ∗ ) against Δχ.

CONCLUSION

Results suggest feasibility of the proposed method as a practical and reliable means for measuring R 2 ' , R₂ , R 2 ∗ , and Δχ across the entire brain.

Tissue Hypoxia and Alterations in Microvascular Architecture Predict Glioblastoma Recurrence in Humans

PURPOSE

Insufficient control of infiltrative glioblastoma (GBM) cells is a major cause of treatment failure and tumor recurrence. Hence, detailed insights into pathophysiologic changes that precede GBM recurrence are needed to develop more precise neuroimaging modalities for tailored diagnostic monitoring and therapeutic approaches.

EXPERIMENTAL DESIGN

Overall, 168 physiologic MRI follow-up examinations of 56 patients with GBM who developed recurrence after standard therapy were retrospectively evaluated, that is, two post-standard-therapeutic follow-ups before and one at radiological recurrence. MRI biomarkers for microvascular architecture and perfusion, neovascularization activity, oxygen metabolism, and hypoxia were determined for brain areas that developed in the further course into recurrence and for the recurrent GBM itself. The temporal pattern of biomarker changes was fitted with locally estimated scatterplot smoothing functions and analyzed for pathophysiologic changes preceding radiological GBM recurrence.

RESULTS

Our MRI approach demonstrated early pathophysiologic changes prior to radiological GBM recurrence in all patients. Analysis of the time courses revealed a model for the pathophysiology of GBM recurrence: 190 days prior to radiological recurrence, vascular cooption by GBM cells induced vessel regression, detected as decreasing vessel density/perfusion and increasing hypoxia. Seventy days later, neovascularization activity was upregulated, which reincreased vessel density and perfusion. Hypoxia, however, continued to intensify for 30 days and peaked 90 days before radiological recurrence.

CONCLUSIONS

Hypoxia may represent an early sign for GBM recurrence. This might become useful in the development of new combined diagnostic-therapeutic approaches for tailored clinical management of recurrent GBM. Further preclinical and in-human studies are required for validation and evaluation.

Relative oxygen extraction fraction (rOEF) MR imaging reveals higher hypoxia in human epidermal growth factor receptor (EGFR) amplified compared with non-amplified gliomas

PURPOSE

Epidermal growth factor receptor (EGFR) amplification promotes gliomagenesis and is linked to lack of oxygen within the tumor microenvironment. Using hypoxia-sensitive spin-and-gradient echo echo-planar imaging and perfusion MRI, we investigated the influence of EGFR amplification on tissue oxygen availability and utilization in human gliomas.

METHODS

This study included 72 histologically confirmed EGFR-amplified and non-amplified glioma patients. Reversible transverse relaxation rate (R₂'), relative cerebral blood volume (rCBV), and relative oxygen extraction fraction (rOEF) were calculated for the contrast-enhancing and non-enhancing tumor regions. Using Student t test or Wilcoxon rank-sum test, median R₂', rCBV, and rOEF were compared between EGFR-amplified and non-amplified gliomas. ROC analysis was performed to assess the ability of imaging characteristics to discriminate EGFR amplification status. Overall survival (OS) was determined using univariate and multivariate cox models. Kaplan-Meier survival curves were plotted and compared using the log-rank test.

RESULTS

EGFR amplified gliomas exhibited significantly higher median R₂' and rOEF than non-amplified gliomas. ROC analysis suggested that R₂' (AUC = 0.7190; P = 0.0048) and rOEF (AUC = 0.6959; P = 0.0156) could separate EGFR status. Patients with EGFR-amplified gliomas had a significantly shorter OS than non-amplified patients. Univariate cox regression analysis determined both R₂' and rOEF significantly influence OS. No significant difference was observed in rCBV between patient cohorts nor was rCBV found to be an effective differentiator of EGFR status.

CONCLUSION

Imaging of tumor oxygen characteristics revealed EGFR-amplified gliomas to be more hypoxic and contribute to shorter patient survival than EGFR non-amplified gliomas.

Oxygen extraction fraction mapping with multi-parametric quantitative BOLD MRI: Reduced transverse relaxation bias using 3D-GraSE imaging

Magnetic resonance imaging (MRI)-based quantification of the blood-oxygenation-level-dependent (BOLD) effect allows oxygen extraction fraction (OEF) mapping. The multi-parametric quantitative BOLD (mq-BOLD) technique facilitates relative OEF (rOEF) measurements with whole brain coverage in clinically applicable scan times. Mq-BOLD requires three separate scans of cerebral blood volume and transverse relaxation rates measured by gradient-echo (1/T₂*) and spin-echo (1/T₂). Although the current method is of clinical merit in patients with stroke, glioma and internal carotid artery stenosis (ICAS), there are relaxation measurement artefacts that impede the sensitivity of mq-BOLD and artificially elevate reported rOEF values.

We posited that T₂-related biases caused by slice refocusing imperfections during rapid 2D-GraSE (Gradient and Spin Echo) imaging can be reduced by applying 3D-GraSE imaging sequences, because the latter requires no slice selective pulses. The removal of T₂-related biases would decrease overestimated rOEF values measured by mq-BOLD. We characterized effects of T₂-related bias in mq-BOLD by comparing the initially employed 2D-GraSE and two proposed 3D-GraSE sequences to multiple single spin-echo reference measurements, both in vitro and in vivo. A phantom and 25 participants, including young and elderly healthy controls as well as ICAS-patients, were scanned. We additionally proposed a procedure to reliably identify and exclude artefact affected voxels. In the phantom, 3D-GraSE derived T₂ values had 57% lower deviation from the reference. For in vivo scans, the formerly overestimated rOEF was reduced by −27% (p < 0.001). We obtained rOEF = 0.51, which is much closer to literature values from positron emission tomography (PET) measurements. Furthermore, increased sensitivity to a focal rOEF elevation in an ICAS-patient was demonstrated.

In summary, the application of 3D-GraSE improves the mq-BOLD-based rOEF quantification while maintaining clinically feasible scan times. Thus, mq-BOLD with non-slice selective T₂ imaging is highly promising to improve clinical diagnostics of cerebrovascular diseases such as ICAS.

No Added Neuroprotective Effect of Remote Ischemic Postconditioning and Therapeutic Hypothermia After Mild Hypoxia-Ischemia in a Piglet Model

Introduction: Hypoxic ischemic encephalopathy (HIE) is a major cause of death and disability in children worldwide. Apart from supportive care, the only established treatment for HIE is therapeutic hypothermia (TH). As TH is only partly neuroprotective, there is a need for additional therapies. Intermittent periods of limb ischemia, called remote ischemic postconditioning (RIPC), have been shown to be neuroprotective after HIE in rats and piglets. However, it is unknown whether RIPC adds to the effect of TH. We tested the neuroprotective effect of RIPC with TH compared to TH alone using magnetic resonance imaging and spectroscopy (MRI/MRS) in a piglet HIE model. Methods: Thirty-two male and female piglets were subjected to 45-min global hypoxia-ischemia (HI). Twenty-six animals were randomized to TH or RIPC plus TH; six animals received supportive care only. TH was induced through whole-body cooling. RIPC was induced 1 h after HI by four cycles of 5 min of ischemia and 5 min of reperfusion in both hind limbs. Primary outcome was Lac/NAA ratio at 24 h measured by MRS. Secondary outcomes were NAA/Cr, diffusion-weighted imaging (DWI), arterial spin labeling, aEGG score, and blood oxygen dependent (BOLD) signal measured by MRI/MRS at 6, 12, and 24 h after the hypoxic-ischemic insult. Results: All groups were subjected to a comparable but mild insult. No difference was found between the two intervention groups in Lac/NAA ratio, NAA/Cr ratio, DWI, arterial spin labeling, or BOLD signal. NAA/Cr ratio at 24 h was higher in the two intervention groups compared to supportive care only. There was no difference in aEEG score between the three groups. Conclusion: Treatment with RIPC resulted in no additional neuroprotection when combined with TH. However, insult severity was mild and only evaluated at 24 h after HI with a short MRS echo time. In future studies more subtle neurological effects may be detected with increased MRS echo time and post mortem investigations, such as brain histology. Thus, the possible neuroprotective effect of RIPC needs further evaluation.

Development of a Non-invasive Assessment of Hypoxia and Neovascularization with Magnetic Resonance Imaging in Benign and Malignant Breast Tumors: Initial Results

PURPOSE

To develop a novel magnetic resonance imaging (MRI) approach for the noninvasive assessment of hypoxia and neovascularization in breast tumors.

PROCEDURES

In this IRB-approved prospective study, 20 patients with suspicious breast lesions (BI-RADS 4/5) underwent multiparametric breast MRI including quantitative BOLD (qBOLD) and vascular architecture mapping (VAM). Custom-made in-house MatLab software was used for qBOLD and VAM data postprocessing and calculation of quantitative MRI biomarker maps of oxygen extraction fraction (OEF), metabolic rate of oxygen (MRO₂), and mitochondrial oxygen tension (mitoPO₂) to measure tissue hypoxia and neovascularization including vascular architecture including microvessel radius (VSI), density (MVD), and type (MTI). Histopathology was used as standard of reference. Appropriate statistics were performed to assess and compare correlations between MRI biomarkers for hypoxia and neovascularization.

RESULTS

qBOLD and VAM data with good quality were obtained from all patients with 13 invasive ductal carcinoma (IDC) and 7 benign breast tumors with a lesion diameter of at least 10 mm in all spatial directions. MRI biomarker maps of oxygen metabolism and neovascularization demonstrated intratumoral spatial heterogeneity with a broad range of biomarker values. Bulk tumor neovasculature consisted of draining venous microvasculature with slow flowing blood. High OEF and low mitoPO₂ were associated with low MVD and vice versa. The heterogeneous pattern of MRO₂ values showed spatial congruence with VSI. IDCs showed significantly higher MRO₂ (P = 0.007), lower mitoPO₂ (P = 0.021), higher MVD (P = 0.005), and lower (i.e., more pathologic) MTI (P = 0.001) compared with benign breast tumors. These results indicate that IDCs consume more oxygen and are more hypoxic and neovascularized than benign tumors.

CONCLUSIONS

We developed a novel MRI approach for the noninvasive assessment of hypoxia and neovascularization in benign and malignant breast tumors that can be easily integrated in a diagnostic MRI protocol and provides insight into intratumoral heterogeneity.

Predicting Glioblastoma Response to Bevacizumab Through MRI Biomarkers of the Tumor Microenvironment

PURPOSE

Glioblastoma (GB) is one of the most vascularized of all solid tumors and, therefore, represents an attractive target for antiangiogenic therapies. Many lesions, however, quickly develop escape mechanisms associated with changes in the tumor microenvironment (TME) resulting in rapid treatment failure. To prevent patients from adverse effects of ineffective therapy, there is a strong need to better predict and monitor antiangiogenic treatment response.

PROCEDURES

We utilized a novel physiological magnetic resonance imaging (MRI) method combining the visualization of oxygen metabolism and neovascularization for classification of five different TME compartments: necrosis, hypoxia with/without neovascularization, oxidative phosphorylation, and aerobic glycolysis. This approach, termed TME mapping, was used to monitor changes in tumor biology and pathophysiology within the TME in response to bevacizumab treatment in 18 patients with recurrent GB.

RESULTS

We detected dramatic changes in the TME by rearrangement of its compartments after the onset of bevacizumab treatment. All patients showed a decrease in active tumor volume and neovascularization as well as an increase in hypoxia and necrosis in the first follow-up after 3 months. We found that recurrent GB with a high percentage of neovascularization and active tumor before bevacizumab onset showed a poor or no treatment response.

CONCLUSIONS

TME mapping might be useful to develop strategies for patient stratification and response prediction before bevacizumab onset.

Characterizing white matter fiber orientation effects on multi-parametric quantitative BOLD assessment of oxygen extraction fraction

Relative oxygen extraction fraction (rOEF) is a fundamental indicator of cerebral metabolic function. An easily applicable method for magnetic resonance imaging (MRI) based rOEF mapping is the multi-parametric quantitative blood oxygenation level dependent (mq-BOLD) approach with separate acquisitions of transverse relaxation times   T 2 * and T2 and dynamic susceptibility contrast (DSC) based relative cerebral blood volume (rCBV). Given that transverse relaxation and rCBV in white matter (WM) strongly depend on nerve fiber orientation, mq-BOLD derived rOEF is expected to be affected as well. To investigate fiber orientation related rOEF artefacts, we present a methodological study characterizing anisotropy effects of WM as measured by diffusion tensor imaging (DTI) on mq-BOLD in 30 healthy volunteers. Using a 3T clinical MRI-scanner, we performed a comprehensive correlation of all parameters (  T 2 * , T2, R 2 ' , rCBV, rOEF, where R 2 ' =1/  T 2 * -1/T2) with DTI-derived fiber orientation towards the main magnetic field (B0). Our results confirm strong dependencies of transverse relaxation and rCBV on the nerve fiber orientation towards B0, with anisotropy-driven variations up to 37%. Comparably weak orientation-dependent variations of mq-BOLD derived rOEF (3.8%) demonstrate partially counteracting influences of R 2 ' and rCBV effects, possibly suggesting applicability of rOEF as an oxygenation sensitive biomarker. However, unresolved issues warrant caution when applying mq-BOLD to WM.

Physiologic MR imaging of the tumor microenvironment revealed switching of metabolic phenotype upon recurrence of glioblastoma in humans

Treating recurrent glioblastoma (GB) is one of the challenges in modern neurooncology. Hypoxia, neovascularization, and energy metabolism are of crucial importance for therapy failure and recurrence. Twenty-one patients with initially untreated GB who developed recurrence were examined with a novel MRI approach for noninvasive visualization of the tumor microenvironment (TME). Imaging biomarker information about oxygen metabolism (mitochondrial oxygen tension) and neovascularization (microvascular density and type) were fused for classification of five different TME compartments: necrosis, hypoxia with/without neovascularization, oxidative phosphorylation, and glycolysis. Volume percentages of these TME compartments were compared between untreated and recurrent GB. At initial diagnosis, all 21 GB showed either the features of a glycolytic dominant phenotype with a high percentage of functional neovasculature (N = 12) or those of a necrotic/hypoxic dominant phenotype with a high percentage of defective tumor neovasculature (N = 9). At recurrence, all 21 GB revealed switching of the initial metabolic phenotype: either from the glycolytic to the necrotic/hypoxic dominant phenotype or vice-versa. A necrotic/hypoxic phenotype at recurrence was associated with a higher rate of multifocality of the recurrent lesions. Our MRI approach may be helpful for a better understanding of treatment-induced metabolic phenotype switching and for future studies developing targeted therapeutic strategies for recurrent GB.

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.

Quantifying Reoxygenation in Pancreatic Cancer During Stereotactic Body Radiotherapy

Hypoxia, the state of low oxygenation that often arises in solid tumours due to their high metabolism and irregular vasculature, is a major contributor to the resistance of tumours to radiation therapy (RT) and other treatments. Conventional RT extends treatment over several weeks or more, and nominally allows time for oxygen levels to increase ("reoxygenation") as cancer cells are killed by RT, mitigating the impact of hypoxia. Recent advances in RT have led to an increase in the use stereotactic body radiotherapy (SBRT), which delivers high doses in five or fewer fractions. For cancers such as pancreatic adenocarcinoma for which hypoxia varies significantly between patients, SBRT might not be optimal, depending on the extent to which reoxygenation occurs during its short duration. We used fluoro-5-deoxy-α-D-arabinofuranosyl)-2-nitroimidazole positron-emission tomography (FAZA-PET) imaging to quantify hypoxia before and after 5-fraction SBRT delivered to patient-derived pancreatic cancer xenografts orthotopically implanted in mice. An imaging technique using only the pre-treatment FAZA-PET scan and repeat dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) scans throughout treatment was able to predict the change in hypoxia. Our results support the further testing of this technique for imaging of reoxygenation in the clinic.

Flow-metabolism uncoupling in patients with asymptomatic unilateral carotid artery stenosis assessed by multi-modal magnetic resonance imaging

Oxygen extraction (OEF), oxidative metabolism (CMRO₂), and blood flow (CBF) in the brain, as well as the coupling between CMRO₂ and CBF due to cerebral autoregulation are fundamental to brain's health. We used a clinically feasible MRI protocol to assess impairments of these parameters in the perfusion territories of stenosed carotid arteries. Twenty-nine patients with unilateral high-grade carotid stenosis and thirty age-matched healthy controls underwent multi-modal MRI scans. Pseudo-continuous arterial spin labeling (pCASL) yielded absolute CBF, whereas multi-parametric quantitative blood oxygenation level dependent (mqBOLD) modeling allowed imaging of relative OEF and CMRO₂. Both CBF and CMRO₂ were significantly reduced in the stenosed territory compared to the contralateral side, while OEF was evenly distributed across both hemispheres similarly in patients and controls. The CMRO₂-CBF coupling was significantly different between both hemispheres in patients, i.e. significant interhemispheric flow-metabolism uncoupling was observed in patients compared to controls. Given that CBF and CMRO₂ are intimately linked to brain function in health and disease, the proposed easily applicable MRI protocol of pCASL and mqBOLD imaging might serve as a valuable tool for early diagnosis of potentially harmful cerebral hemodynamic and metabolic states with the final aim to select clinically asymptomatic patients who would benefit from carotid revascularization therapy.

Intratumoral heterogeneity of oxygen metabolism and neovascularization uncovers 2 survival-relevant subgroups of IDH1 wild-type glioblastoma

Background

The intratumoral heterogeneity of oxygen metabolism in combination with variable patterns of neovascularization (NV) as well as reprogramming of energy metabolism affects the landscape of tumor microenvironments (TMEs) in glioblastoma. Knowledge of the hypoxic and perivascular niches within the TME is essential for understanding treatment failure.

Methods

Fifty-two patients with untreated glioblastoma (isocitrate dehydrogenase 1 wild type [IDH1wt]) were examined with a physiological MRI protocol including a multiparametric quantitative blood oxygen level dependent (qBOLD) approach and vascular architecture mapping (VAM). Imaging biomarker information about oxygen metabolism (mitochondrial oxygen tension) and neovascularization (microvascular density and type) were fused for classification of 6 different TMEs: necrosis, hypoxia with/without neovascularization, oxidative phosphorylation (OxPhos), and glycolysis with/without neovascularization. Association of the different TME volume fractions with progression-free survival (PFS) was assessed using Kaplan-Meier analysis and Cox proportional hazards models.

Results

A common spatial structure of TMEs was detected: central necrosis surrounded by tumor hypoxia (with defective and functional neovasculature) and different TMEs with a predominance of OxPhos and glycolysis for energy production, respectively. The percentage of the different TMEs on the total tumor volume uncovered 2 clearly different subtypes of glioblastoma IDH1wt: a glycolytic dominated phenotype with predominantly functional neovasculature and a necrotic/hypoxic dominated phenotype with approximately 50% of defective neovasculature. Patients with a necrotic/hypoxic dominated phenotype showed significantly shorter PFS (P = 0.035).

Conclusions

Our non-invasive mapping approach allows for classification of the TME and detection of tumor-supportive niches in glioblastoma which may be helpful for both clinical patient management and research.

Using an artificial neural network for fast mapping of the oxygen extraction fraction with combined QSM and quantitative BOLD

PURPOSE

To apply an artificial neural network (ANN) for fast and robust quantification of the oxygen extraction fraction (OEF) from a combined QSM and quantitative BOLD analysis of gradient echo data and to compare the ANN to a traditional quasi-Newton (QN) method for numerical optimization.

METHODS

Random combinations of OEF, deoxygenated blood volume ( ν ), R₂ , and nonblood magnetic susceptibility ( χ nb ) with each parameter following a Gaussian distribution that represented physiological gray matter and white matter values were used to simulate quantitative BOLD signals and QSM values. An ANN was trained with the simulated data with added Gaussian noise. The ANN was applied to multigradient echo brain data of 7 healthy subjects, and the reconstructed parameters and maps were compared to QN results using Student t test and Bland-Altman analysis.

RESULTS

Intersubject means and SDs of gray matter were OEF = 43.5 ± 0.8 %, R 2 = 13.5 ± 0.3 Hz, ν = 3.4 ± 0.1 %, χ nb = - 25 ± 5 ppb for ANN; and OEF = 43.8 ± 5.2 %, R 2 = 12.2 ± 0.8 Hz, ν = 4.2 ± 0.6 %, χ nb = - 39 ± 7 ppb for QN, with a significant difference ( P < 0.05 ) for R 2 , ν , and χ nb . For white matter, they were OEF = 47.5 ± 1.1 %, R 2 = 17.1 ± 0.4 Hz, ν = 2.5 ± 0.2 %, χ nb = - 38 ± 5 ppb for ANN; and OEF = 42.3 ± 5.6 %, R 2 = 16.7 ± 0.7 Hz, ν = 2.9 ± 0.3 %, χ nb = - 45 ± 9 ppb for QN, with a significant difference ( P < 0.05 ) for OEF and ν . ANN revealed more gray-white matter contrast but less intersubject variation in OEF than QN. In contrast to QN, the ANN reconstruction did not need an additional sequence for parameter initialization and took approximately 1 s rather than roughly 1 h.

CONCLUSION

ANNs allow faster and, with regard to initialization, more robust reconstruction of OEF maps with lower intersubject variation than QN approaches.

Prospects for investigating brain oxygenation in acute stroke: Experience with a non-contrast quantitative BOLD based approach

Metabolic markers of baseline brain oxygenation and tissue perfusion have an important role to play in the early identification of ischaemic tissue in acute stroke. Although well established MRI techniques exist for mapping brain perfusion, quantitative imaging of brain oxygenation is poorly served. Streamlined-qBOLD (sqBOLD) is a recently developed technique for mapping oxygenation that is well suited to the challenge of investigating acute stroke. In this study a noninvasive serial imaging protocol was implemented, incorporating sqBOLD and arterial spin labelling to map blood oxygenation and perfusion, respectively. The utility of these parameters was investigated using imaging based definitions of tissue outcome (ischaemic core, infarct growth and contralateral tissue). Voxel wise analysis revealed significant differences between all tissue outcomes using pairwise comparisons for the transverse reversible relaxation rate (R ₂ '), deoxygenated blood volume (DBV) and deoxyghaemoglobin concentration ([dHb]; p < 0.01 in all cases). At the patient level (n = 9), a significant difference was observed for [dHb] between ischaemic core and contralateral tissue. Furthermore, serial analysis at the patient level (n = 6) revealed significant changes in R ₂ ' between the presentation and 1 week scans for both ischaemic core (p < 0.01) and infarct growth (p < 0.05). In conclusion, this study presents evidence supporting the potential of sqBOLD for imaging oxygenation in stroke.

The relationship between blood flow impairment and oxygen depletion in acute ischemic stroke imaged with magnetic resonance imaging

Oxygenation-sensitive spin relaxation time T2' and relaxation rate R2' (1/T2') are presumed to be markers of the cerebral oxygen extraction fraction (OEF) in acute ischemic stroke. In this study, we investigate the relationship of T2'/R2' with dynamic susceptibility contrast-based relative cerebral blood flow (rCBF) in acute ischemic stroke to assess their plausibility as surrogate markers of the ischemic penumbra. Twenty-one consecutive patients with internal carotid artery and/or middle cerebral artery occlusion were studied at 3.0 T. A physiological model of the cerebral vasculature (VM) was used to process PWI raw data in addition to a conventional deconvolution technique. T2', R2', and rCBF values were extracted from the ischemic core and hypoperfused areas. Within hypoperfused tissue, no correlation was found between deconvolved rCBF and T2' ( r = -0.05, p = 0.788), or R2' ( r = 0.039, p = 0.836). In contrast, we found a strong positive correlation with T2' ( r = 0.444, p = 0.006) and negative correlation with R2' ( r = -0.494, p = 0.0025) for rCBF_VM, indicating increasing OEF with decreasing CBF and that rCBF based on the vascular model may be more closely related to metabolic disturbances. Further research to refine and validate these techniques may enable their use as MRI-based surrogate markers of the ischemic penumbra for selecting stroke patients for interventional treatment strategies.

Oxygen metabolism in acute ischemic stroke

Gaining insights into brain oxygen metabolism has been one of the key areas of research in neurosciences. Extensive efforts have been devoted to developing approaches capable of providing measures of brain oxygen metabolism not only under normal physiological conditions but, more importantly, in various pathophysiological conditions such as cerebral ischemia. In particular, quantitative measures of cerebral metabolic rate of oxygen using positron emission tomography (PET) have been shown to be capable of discerning brain tissue viability during ischemic insults. However, the complex logistics associated with oxygen-15 PET have substantially hampered its wide clinical applicability. In contrast, magnetic resonance imaging (MRI)-based approaches have provided quantitative measures of cerebral oxygen metabolism similar to that obtained using PET. Given the wide availability, MRI-based approaches may have broader clinical impacts, particularly in cerebral ischemia, when time is a critical factor in deciding treatment selection. In this article, we review the pathophysiological basis of altered cerebral hemodynamics and oxygen metabolism in cerebral ischemia, how quantitative measures of cerebral metabolism were obtained using the Kety-Schmidt approach, the physical concepts of non-invasive oxygen metabolism imaging approaches, and, finally, clinical applications of the discussed imaging approaches.

Interleaved quantitative BOLD: Combining extravascular R2' - and intravascular R2-measurements for estimation of deoxygenated blood volume and hemoglobin oxygen saturation

Quantitative BOLD (qBOLD), a non-invasive MRI method for assessment of hemodynamic and metabolic properties of the brain in the baseline state, provides spatial maps of deoxygenated blood volume fraction (DBV) and hemoglobin oxygen saturation (HbO₂) by means of an analytical model for the temporal evolution of free-induction-decay signals in the extravascular compartment. However, mutual coupling between DBV and HbO₂ in the signal model results in considerable estimation uncertainty precluding achievement of a unique set of solutions. To address this problem, we developed an interleaved qBOLD method (iqBOLD) that combines extravascular R₂' and intravascular R₂ mapping techniques so as to obtain prior knowledge for the two unknown parameters. To achieve these goals, asymmetric spin echo and velocity-selective spin-labeling (VSSL) modules were interleaved in a single pulse sequence. Prior to VSSL, arterial blood and CSF signals were suppressed to produce reliable estimates for cerebral venous blood volume fraction (CBV_v) as well as venous blood R₂ (to yield HbO₂). Parameter maps derived from the VSSL module were employed to initialize DBV and HbO₂ in the qBOLD processing. Numerical simulations and in vivo experiments at 3 T were performed to evaluate the performance of iqBOLD in comparison to the parent qBOLD method. Data obtained in eight healthy subjects yielded plausible values averaging 60.1 ± 3.3% for HbO₂ and 3.1 ± 0.5 and 2.0 ± 0.4% for DBV in gray and white matter, respectively. Furthermore, the results show that prior estimates of CBV_v and HbO₂ from the VSSL component enhance the solution stability in the qBOLD processing, and thus suggest the feasibility of iqBOLD as a promising alternative to the conventional technique for quantifying neurometabolic parameters.

Recurrence of glioblastoma is associated with elevated microvascular transit time heterogeneity and increased hypoxia

Dynamic susceptibility contrast (DSC) perfusion MRI provide information about differences in macro- and microvasculature when executed with gradient-echo (GE; sensitive to macrovasculature) and spin-echo (SE; sensitive to microvasculature) contrast. This study investigated whether there are differences between macro- and microvascular transit time heterogeneity (MVTH and µVTH) and tissue oxygen tension (PO2_mit) in newly-diagnosed and recurrent glioblastoma. Fifty-seven patients with glioblastoma (25 newly-diagnosed/32 recurrent) were examined with GE- and SE-DSC perfusion sequences, and a quantitative blood-oxygen-level-dependent (qBOLD) approach. Maps of MVTH, µVTH and coefficient of variation (MCOV and µCOV) were calculated from GE- and SE-DSC data, respectively, using an extended flow-diffusion equation. PO2_mit maps were calculated from qBOLD data. Newly-diagnosed and recurrent glioblastoma showed significantly lower ( P ≤ 0.001) µCOV values compared to both normal brain and macrovasculature (MCOV) of the lesions. Recurrent glioblastoma had significantly higher µVTH ( P = 0.014) and µCOV ( P = 0.039) as well as significantly lower PO2_mit values ( P = 0.008) compared to newly-diagnosed glioblastoma. The macrovasculature, however, showed no significant differences. Our findings provide evidence of microvascular adaption in the disorganized tumor vasculature for retaining the metabolic demands in stress response of therapeutically-uncontrolled glioblastomas. Thus, µVTH and PO2_mit mapping gives insight into the tumor microenvironment (vascular and hypoxic niches) responsible for therapy resistance.

Characterizing hypoxia in human glioma: A simultaneous multimodal MRI and PET study

Hypoxia plays an important role for the prognosis and therapy response of cancer. Thus, hypoxia imaging would be a valuable tool for pre-therapeutic assessment of tumor malignancy. However, there is no standard validated technique for clinical application available yet. Therefore, we performed a study in 12 patients with high-grade glioma, where we directly compared the two currently most promising techniques, namely the MR-based relative oxygen extraction fraction (MR-rOEF) and the PET hypoxia marker H-1-(3-[¹⁸ F]-fluoro-2-hydroxypropyl)-2-nitroimidazole ([¹⁸ F]-FMISO). MR-rOEF was determined from separate measurements of T₂ , T₂ * and relative cerebral blood volume (rCBV) employing a multi-parametric approach for quantification of the blood-oxygenation-level-dependent (BOLD) effect. With respect to [¹⁸ F]-FMISO-PET, besides the commonly used late uptake between 120 and 130 min ([¹⁸ F]-FMISO_{120-130 min} ), we also analyzed the hypoxia specific uptake rate [¹⁸ F]-FMISO-k₃ , as obtained by pharmacokinetic modeling of dynamic uptake data. Since pharmacokinetic modeling of partially acquired dynamic [¹⁸ F]-FMISO data was sensitive to a low signal-to-noise-ratio, analysis was restricted to high-uptake tumor regions. Individual spatial analyses of deoxygenation and hypoxia-related parameter maps revealed that high MR-rOEF values clustered in (edematous) peritumoral tissue, while areas with high [¹⁸ F]-FMISO_{120-130 min} concentrated in and around active tumor with disrupted blood-brain barrier, i.e. contrast enhancement in T₁ -weighted MRI. Volume-of-interest-based correlations between MR-rOEF and [¹⁸ F]-FMISO_{120-130 min} as well as [¹⁸ F]-FMISO-k₃ , and voxel-wise analyses in individual patients, yielded limited correlations, supporting the notion that [¹⁸ F]-FMISO uptake, even after 2 h, might still be influenced by perfusion while [¹⁸ F]-FMISO-k₃ was severely hampered by noise. According to these results, vascular deoxygenation, as measured by MR-rOEF, and severe tissue hypoxia, as measured by [¹⁸ F]-FMISO, show a poor spatial correspondence. Overall, the two methods appear to rather provide complementary than redundant information about high-grade glioma biology.

Delineation of Tumor Habitats based on Dynamic Contrast Enhanced MRI

Tumor heterogeneity can be elucidated by mapping subregions of the lesion with differential imaging characteristics, called habitats. Dynamic Contrast Enhanced (DCE-)MRI can depict the tumor microenvironments by identifying areas with variable perfusion and vascular permeability, since individual tumor habitats vary in the rate and magnitude of the contrast uptake and washout. Of particular interest is identifying areas of hypoxia, characterized by inadequate perfusion and hyper-permeable vasculature. An automatic procedure for delineation of tumor habitats from DCE-MRI was developed as a two-part process involving: (1) statistical testing in order to determine the number of the underlying habitats; and (2) an unsupervised pattern recognition technique to recover the temporal contrast patterns and locations of the associated habitats. The technique is examined on simulated data and DCE-MRI, obtained from prostate and brain pre-clinical cancer models, as well as clinical data from sarcoma and prostate cancer patients. The procedure successfully identified habitats previously associated with well-perfused, hypoxic and/or necrotic tumor compartments. Given the association of tumor hypoxia with more aggressive tumor phenotypes, the obtained in vivo information could impact management of cancer patients considerably.

Contributions to magnetic susceptibility of brain tissue

This review discusses the major contributors to the subtle magnetic properties of brain tissue and how they affect MRI contrast. With the increased availability of high-field scanners, the use of magnetic susceptibility contrast for the study of human brain anatomy and function has increased dramatically. This has not only led to novel applications, but has also improved our understanding of the complex relationship between MRI contrast and magnetic susceptibility. Chief contributors to the magnetic susceptibility of brain tissue have been found to include myelin as well as iron. In the brain, iron exists in various forms with diverse biological roles, many of which are now only starting to be uncovered. An interesting aspect of magnetic susceptibility contrast is its sensitivity to the microscopic distribution of iron and myelin, which provides opportunities to extract information at spatial scales well below MRI resolution. For example, in white matter, the myelin sheath that surrounds the axons can provide tissue contrast that is dependent on the axonal orientation and reflects the relative size of intra- and extra-axonal water compartments. The extraction of such ultrastructural information, together with quantitative information about iron and myelin concentrations, is an active area of research geared towards the characterization of brain structure and function, and their alteration in disease. Copyright © 2016 John Wiley & Sons, Ltd.

Contributions to magnetic susceptibility of brain tissue

This review discusses the major contributors to the subtle magnetic properties of brain tissue and how they affect MRI contrast. With the increased availability of high-field scanners, the use of magnetic susceptibility contrast for the study of human brain anatomy and function has increased dramatically. This has not only led to novel applications, but has also improved our understanding of the complex relationship between MRI contrast and magnetic susceptibility. Chief contributors to the magnetic susceptibility of brain tissue have been found to include myelin as well as iron. In the brain, iron exists in various forms with diverse biological roles, many of which are now only starting to be uncovered. An interesting aspect of magnetic susceptibility contrast is its sensitivity to the microscopic distribution of iron and myelin, which provides opportunities to extract information at spatial scales well below MRI resolution. For example, in white matter, the myelin sheath that surrounds the axons can provide tissue contrast that is dependent on the axonal orientation and reflects the relative size of intra- and extra-axonal water compartments. The extraction of such ultrastructural information, together with quantitative information about iron and myelin concentrations, is an active area of research geared towards the characterization of brain structure and function, and their alteration in disease. Copyright © 2016 John Wiley & Sons, Ltd.

Intraoperative Magnetic Resonance Imaging of Cerebral Oxygen Metabolism During Resection of Brain Lesions

OBJECTIVE

Tissue oxygen tension is an important parameter for brain tissue viability and its noninvasive intraoperative monitoring in the whole brain is of highly clinical relevance. The purpose of this study was the introduction of a multiparametric quantitative blood oxygenation dependent magnetic resonance imaging (MRI) approach for intraoperative examination of oxygen metabolism during the resection of brain lesions.

METHODS

Sixteen patients suffering from brain lesions were examined intraoperatively twice (before craniotomy and after gross-total resection) via the quantitative blood oxygenation dependent technique and a 1.5-Tesla MRI scanner, which is installed in an operating room. The MRI protocol included T2*- and T2 mapping and dynamic susceptibility weighted perfusion. Data analysis was performed with a custom-made, in-house MatLab software for calculation of maps of oxygen extraction fraction (OEF) and cerebral metabolic rate of oxygen (CMRO₂) as well as of cerebral blood volume and cerebral blood flow.

RESULTS

Perilesional edema showed a significant increase in both perfusion (cerebral blood volume +21%, cerebral blood flow +13%) and oxygen metabolism (OEF +32%, CMRO₂ +16%) after resection of the lesions. In perilesional nonedematous tissue only, however, oxygen metabolism (OEF +19%, CMRO₂ +11%) was significantly increased, but not perfusion. No changes were found in normal brain. Fortunately, no neurovascular adverse events were observed.

CONCLUSIONS

This approach for intraoperative examination of oxygen metabolism in the whole brain is a new application of intraoperative MRI additionally to resection control (residual tumor detection) and updating of neuronavigation (brain shift detection). It may help to detect neurovascular adverse events early during surgery.

MR Imaging-derived Oxygen Metabolism and Neovascularization Characterization for Grading and IDH Gene Mutation Detection of Gliomas

Purpose To explore the diagnostic performance of physiological magnetic resonance (MR) imaging of oxygen metabolism and neovascularization activity for grading and characterization of isocitrate dehydrogenase (IDH) gene mutation status of gliomas. Materials and Methods This retrospective study had institutional review board approval; written informed consent was obtained from all patients. Eighty-three patients with histopathologically proven glioma (World Health Organization [WHO] grade II-IV) were examined with quantitative blood oxygen level-dependent imaging and vascular architecture mapping. Biomarker maps of neovascularization activity (microvessel radius, microvessel density, and microvessel type indicator [MTI]) and oxygen metabolism (oxygen extraction fraction [OEF] and cerebral metabolic rate of oxygen [CMRO₂]) were calculated. Receiver operating characteristic analysis was used to determine diagnostic performance for grading and detection of IDH gene mutation status. Results Low-grade (WHO grade II) glioma showed areas with increased OEF (+18%, P < .001, n = 20), whereas anaplastic glioma (WHO grade III) and glioblastoma (WHO grade IV) showed decreased OEF when compared with normal brain tissue (-54% [P < .001, n = 21] and -49% [P < .001, n = 41], respectively). This allowed clear differentiation between low- and high-grade glioma (area under the receiver operating characteristic curve [AUC], 1) for the patient cohort. MTI had the highest diagnostic performance (AUC, 0.782) for differentiation between gliomas of grades III and IV among all biomarkers. CMRO₂ was decreased (P = .037) in low-grade glioma with a mutated IDH gene, and MTI was significantly increased in glioma grade III with IDH mutation (P = .013) when compared with the IDH wild-type counterparts. CMRO₂ showed the highest diagnostic performance for IDH gene mutation detection in low-grade glioma (AUC, 0.818) and MTI in high-grade glioma (AUC, 0.854) and for all WHO grades (AUC, 0.899) among all biomarkers. Conclusion MR imaging-derived oxygen metabolism and neovascularization characterization may be useful for grading and IDH mutation detection of gliomas and requires only 7 minutes of extra imaging time. ^© RSNA, 2016 Online supplemental material is available for this article.

A streamlined acquisition for mapping baseline brain oxygenation using quantitative BOLD

Quantitative BOLD (qBOLD) is a non-invasive MR technique capable of producing quantitative measurements of the haemodynamic and metabolic properties of the brain. Here we propose a refinement of the qBOLD methodology, dubbed streamlined-qBOLD, in order to provide a clinically feasible method for mapping baseline brain oxygenation. In streamlined-qBOLD confounding signal contributions are minimised during data acquisition through the application of (i) a Fluid Attenuated Inversion Recovery (FLAIR) preparation to remove cerebral spinal fluid (CSF) signal contamination, (ii) a Gradient Echo Slice Excitation Profile Imaging (GESEPI) acquisition to reduce the effect of macroscopic magnetic field gradients and (iii) an Asymmetric Spin Echo (ASE) pulse sequence to directly measure the reversible transverse relaxation rate, R₂'. Together these features simplify the application of the qBOLD model, improving the robustness of the resultant parametric maps. A theoretical optimisation framework was used to optimise acquisition parameters in relation to signal to noise ratio. In a healthy subject group (n = 7) apparent elevations in R₂' caused by partial volumes of CSF were shown to be reduced with the application of CSF nulling. Significant decreases in R₂' (p < 0.001) and deoxygenated blood volume (p < 0.01) were seen in cortical grey matter, across the group, with the application of CSF suppression. Quantitative baseline brain oxygenation parameter maps were calculated using qBOLD modelling and compared with literature values.

Developmental trajectories of cerebral blood flow and oxidative metabolism at baseline and during working memory tasks

The neurobiological interpretation of developmental BOLD fMRI findings remains difficult due to the confounding issues of potentially varied baseline of brain function and varied strength of neurovascular coupling across age groups. The central theme of the present research is to study the development of brain function and neuronal activity through in vivo assessments of cerebral blood flow (CBF), oxygen extraction fraction (OEF) and cerebral metabolic rate of oxygen (CMRO2) both at baseline and during the performance of a working memory task in a cohort of typically developing children aged 7 to 18years. Using a suite of 4 emerging MRI technologies including MR blood oximetry, phase-contrast MRI, pseudo-continuous arterial spin labeling (pCASL) perfusion MRI and concurrent CBF/BOLD fMRI, we found: 1) At baseline, both global CBF and CMRO2 showed an age related decline while global OEF was stable across the age group; 2) During the working memory task, neither BOLD nor CBF responses showed significant variations with age in the activated fronto-parietal brain regions. Nevertheless, detailed voxel-wise analyses revealed sub-regions within the activated fronto-parietal regions that show significant decline of fractional CMRO2 responses with age. These findings suggest that the brain may become more "energy efficient" with age during development.

Arterial Spin Labeling Techniques 2009-2014

PURPOSE

Arterial spin labeling (ASL) techniques have been implemented across a diverse range of clinical and experimental applications. This review aims to evaluate the current feasibility of ASL in clinical neuroradiology based on recent improvements to ASL sequences and highlight areas for potential clinical applications.

METHODS AND MATERIALS

In December 2014, a literature search was conducted on PubMed Central, EMBASE, and Scopus using the search terms: "arterial spin labeling, neuroradiology," for studies published between 2009 and 2014 (inclusive). Of 483 studies matching the inclusion criteria, the number of studies using continuous, pseudocontinuous, pulsed, and velocity-selective ASL sequences was 42, 209, 226, and 3, respectively. Studies were classified based on several common clinical applications according to the type of ASL sequence used. Studies using pulsed ASL and pseudo-continuous ASL were grouped based on common sequences.

RESULTS

The number of clinical studies was 264. Numerous studies applied ASL to stroke management (43 studies), drug testing (21 studies), neurodegenerative diseases (40 studies), and psychiatric disorders (26 studies).

CONCLUSIONS

This review discusses several factors hindering the implementation of clinical ASL and ASL-related radiofrequency safety issues encountered in clinical practice. However, a limited number of search terms were used. Further development of robust sequences with multislice imaging capabilities and reduced radiofrequency energy deposition will hopefully improve the clinical acceptance of ASL.

Reduced CMRO₂ and cerebrovascular reserve in patients with severe intracranial arterial stenosis: a combined multiparametric qBOLD oxygenation and BOLD fMRI study

Multiparametric quantitative blood oxygenation level dependent (mqBOLD) magnetic resonance Imaging (MRI) approach allows mapping tissular oxygen saturation (StO2 ) and cerebral metabolic rate of oxygen (CMRO2 ). To identify hemodynamic alteration related to severe intracranial arterial stenosis (SIAS), functional MRI of cerebrovascular reserve (CVR BOLD fMRI) to hypercapnia has been proposed. Diffusion imaging suggests chronic low grade ischemia in patients with impaired CVR. The aim of the present study was to evaluate how oxygen parameters (StO2 and CMRO2 ), assessed with mqBOLD approach, correlate with CVR in patients (n = 12) with SIAS and without arterial occlusion. The perfusion (dynamic susceptibility contrast), oxygenation, and CVR were compared. The MRI protocol conducted at 3T lasted approximately 1 h. Regions of interest measures on maps were delineated on segmented gray matter (GM) of middle cerebral artery territories. We have shown that decreased CVR is spatially associated with decreased CMRO2 in GM of patients with SIAS. Further, the degree of ipsilateral CVR reduction was well-correlated with the amplitude of the CMRO2 deficit. The altered CMRO2 suggests the presence of a moderate ischemia explained by both a decrease in perfusion and in CVR. CVR and mqBOLD method may be helpful in the selection of patients with SIAS to advocate for medical therapy or percutaneous transluminal angioplasty-stenting.

Technical considerations on the validity of blood oxygenation level-dependent-based MR assessment of vascular deoxygenation

A blood oxygenation level-dependent (BOLD)-based apparent relative oxygen extraction fraction (rOEF) as a semi-quantitative marker of vascular deoxygenation has recently been introduced in clinical studies of patients with glioma and stroke, yielding promising results. These rOEF measurements are based on independent quantification of the transverse relaxation times T2 and T2* and relative cerebral blood volume (rCBV). Simulations demonstrate that small errors in any of the underlying measures may result in a large deviation of the calculated rOEF. Therefore, we investigated the validity of such measurements. For this, we evaluated the quantitative measurements of T2 and T2* at 3 T in a gel phantom, in healthy subjects and in healthy tissue of patients with brain tumors. We calculated rOEF maps covering large portions of the brain from T2, T2* and rCBV [routinely measured in patients using dynamic susceptibility contrast (DSC)], and obtained rOEF values of 0.63 ± 0.16 and 0.90 ± 0.21 in healthy-appearing gray matter (GM) and white matter (WM), respectively; values of about 0.4 are usually reported. Quantitative T2 mapping using the fast, clinically feasible, multi-echo gradient spin echo (GRASE) approach yields significantly higher values than much slower multiple single spin echo (SE) experiments. Although T2* mapping is reliable in magnetically homogeneous tissues, uncorrectable macroscopic background gradients and other effects (e.g. iron deposition) shorten T2*. Cerebral blood volume (CBV) measurement using DSC and normalization to WM yields robust estimates of rCBV in healthy-appearing brain tissue; absolute quantification of the venous fraction of CBV, however, is difficult to achieve. Our study demonstrates that quantitative measurements of rOEF are currently biased by inherent difficulties in T2 and CBV quantification, but also by inadequacies of the underlying model. We argue, however, that standardized, reproducible measurements of apparent T2, T2* and rCBV may still allow the estimation of a meaningful apparent rOEF, which requires further validation in clinical studies.

Non-invasive multiparametric qBOLD approach for robust mapping of the oxygen extraction fraction

INTRODUCTION

The quantitative blood oxygenation level-dependent (qBOLD) method has not become clinically established yet because long acquisition times are necessary to achieve an acceptable certainty of the parameter estimates. In this work, a non-invasive multiparametric (nimp) qBOLD approach based on a simple analytical model is proposed to facilitate robust oxygen extraction fraction (OEF) mapping within clinically acceptable acquisition times by using separate measurements.

METHODS

The protocol consisted of a gradient-echo sampled spin-echo sequence (GESSE), a T2-weighted Carr-Purcell-Meiboom-Gill (CPMG) sequence, and a T2(*)-weighted multi-slice multi-echo gradient echo (MMGE) sequence. The GESSE acquisition time was less than 5 minutes and the extra measurement time for CPMG/MMGE was below 2 minutes each. The proposed nimp-qBOLD approach was validated in healthy subjects (N = 5) and one patient.

RESULTS

The proposed nimp-qBOLD approach facilitated more robust OEF mapping with significantly reduced inter- and intra-subject variability compared to the standard qBOLD method. Thereby, an average OEF in all subjects of 27±2% in white matter (WM) and 29±2% in gray matter (GM) using the nimp-qBOLD method was more stable compared to 41±10% (WM) and 46±10% (GM) with standard qBOLD. Moreover, the spatial variance in the image slice (i.e. standard deviation divided by mean) was on average reduced from 35% to 25%. In addition, the preliminary results of the patient are encouraging.

CONCLUSION

The proposed nimp-qBOLD technique provides a promising tool for robust OEF mapping within clinically acceptable acquisition times and could therefore provide an important contribution for analyzing tumors or monitoring the success of radio and chemo therapies.

Numerical modeling of susceptibility-related MR signal dephasing with vessel size measurement: phantom validation at 3T

PURPOSE

MRI is used to obtain quantitative oxygenation and blood volume information from the susceptibility-related MR signal dephasing induced by blood vessels. However, analytical models that fit the MR signal are usually not accurate over the range of small blood vessels. Moreover, recent studies have demonstrated limitations in the simultaneous assessment of oxygenation and blood volume. In this study, a multiparametric MRI framework that aims to measure vessel radii in addition to magnetic susceptibility and volume fraction was introduced.

METHODS

The protocol consisted of gradient-echo sampling of the spin-echo, diffusion, T2, and B0 acquisitions. After correction steps, the data were postprocessed with a versatile numerical model of the MR signal. An important analytical model was implemented for comparison. The approach was validated in phantoms with coiling strings as proxy for blood vessels.

RESULTS

The feasibility of the vessel radius measurement is demonstrated. The numerical model shows an improved accuracy compared with the analytical approach. However, both methods overestimate the radius. The simultaneous measurement of the magnetic susceptibility and the volume fraction remains challenging.

CONCLUSION

The results suggest that this approach could be interesting in vivo to better characterize the microvasculature without contrast agent.