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.

Introduction

Reprogramming oxygen metabolism and inducing angiogenesis are counted among the hallmarks of cancer (Hanahan and Weinberg, 2011). High-grade gliomas, such as anaplastic gliomas and glioblastomas (GBMs), account for the majority of malignant brain tumors in adults and have the characteristics of extensive hypoxia, high vascularization, and heterogeneity (Kaur et al., 2005; Hardee and Zagzag, 2012). Importantly, early identification of distinct genetic profiles such as IDH and O⁶-methylguanine-DNA-methyltransferase (MGMT) in combination with relevant receptor tyrosine kinases (RTKs) highlights the mutational profile of oxygen metabolism or angiogenesis, thus enabling early therapeutic intervention in patients (Chahal et al., 2010; Gluck et al., 2015; Kickingereder et al., 2015). Gliomas with IDH1 mutation were found to have significantly reduced hypoxia-inducible factor 1α (HIF-1α) and decreased neovascularization (Semenza, 2003). Methylation status of the MGMT promoter is also involved in carcinogenesis, as is highlighted by the association with angiogenic profile in GBM (Chahal et al., 2010). Additionally, accumulating evidence suggests that dynamic changes in coactivated RTK pathway involved in GBM may account for the oncogenic processes where tumor hypoxia and angiogenesis frequently coexist (Gluck et al., 2015).

Although molecular imaging probes under development for positron emission tomography and other imaging modalities aim to more directly evaluate tumor oxygenation, these methods are limited by the need of (i) specialized imaging agents with high costs and (ii) complex physiologic assumptions in data interpretation (la Fougere et al., 2011; Bulte et al., 2012; Paech et al., 2020). Moreover, this information currently confronts tissue sampling errors, as a single specimen/sample might not reflect the presence of the mutation in such a heterogeneous tumor. Consequently, alternative non-invasive magnetic resonance imaging (MRI) sequences to predict genetic profiles via its reliable oxygen information are urgently needed.

Nowadays, a novel developed approach, cluster analysis of time evolution (CAT), for quantitative susceptibility mapping and quantitative blood oxygen level–dependent magnitude (QSM + qBOLD or QQ)–based oxygen extraction fraction (OEF) mapping, might have the potential to cope with this problem (Cho et al., 2018, 2020b). The robustness of QQ-based OEF mapping has been substantially improved by introduction of an unsupervised machine learning method, CAT, which may enable clinically practical evaluation of oxygen information without vascular challenges. Because of the relative simplicity of the susceptibility measurement and an acceptable time requirement, CAT for QQ represents a non-invasive measure of quantitative OEF that meets the immediate applicability in clinical practice (Cho et al., 2020a,b). Furthermore, dynamic contrast-enhanced (DCE) MRI has long been clinically used to investigate tumor angiogenesis and plays a pivotal role for characterization of tumor microvessel proliferation and permeability, cerebral blood flow (CBF), or cerebral blood volume (CBV) (Jensen et al., 2014; Arevalo-Perez et al., 2015).

The purpose of this study was to investigate the diagnostic performance of CAT for QQ-based OEF mapping and apply DCE-MRI to predict not only the molecular parameters such as IDH1 mutation and MGMT promoter methylation status in glioma but also differentiation of World Health Organization (WHO) grade II/III lower-grade glioma (LGG) vs. WHO grade IV GBM. Impressively, we also tried to highlight the association of physiological MRI and molecular stratification via the combination of RTK genetic aberrations rather than assessing the status of individual mutations, potentially providing new therapeutic opportunities against these deadly brain tumors. Altogether, we hypothesized that CAT for QQ-based OEF mapping might be a potential method to specifically quantify oxygen metabolism in glioma, this method in combination with characterization of tumor angiogenesis assessed by DCE-MRI can potentially be used as an imaging biomarker to non-invasively identify genetic profiles in the preoperative workup of glioma patients. The established imaging signatures are derived from routine clinically acquired MRI, and therefore, it is easily translatable to the clinic.

Materials and Methods

Patient Selection

This retrospective study was performed with approval of the local institutional review board, and written informed consent was waived. Between July 2016 and December 2019, 138 consecutive patients with suspected primary gliomas were enrolled in this study. Forty-seven patients (34%) were excluded because of the following criteria: (a) diagnosis other than glioma (21 patients), (b) recurrent glioma undergoing therapy (12 patients), and (c) insufficient data quality, in the form of patient motion during MRI (three patients with severe neurological deficits) or poor contrast material injection (11 patients). A total of 91 pathologically confirmed gliomas were identified in our analysis. Detailed patient characteristics are further profiled in Table 1.

TABLE 1

Table 1. Patient demographics and genetic information.

Magnetic Resonance Imaging Acquisition and Data Analysis

Images including conventional MRI scans, OEF images, and DCE perfusion images were acquired during clinical workup on a 3.0-T MR scanner (Discovery 750, GE Healthcare, United States). All images were coregistered to the three-dimensional (3D) T1-BRAVO images, using FMRIB’s Linear Image Registration Tool algorithm (FSL).¹ Details on MRI acquisition parameters and the complete analysis workflow (performed by Drs. Shen and Xie, with 8 and 5 years of experience in brain tumor image processing and interpretation of glioma imaging data, respectively) are outlined in Supplementary Methods 1.

Magnetic Resonance Imaging-Based Oxygen Extraction Fraction Mapping Post-processing

To achieve OEF = 1 − Y/Y_a, the following steps were used for MRI-based OEF postprocessing. (a) QSM images were reconstructed from multi-echo gradient echo (mGRE) data using a fully automated zero-referenced morphology-enabled dipole inversion method as described previously (Liu et al., 2018). QSM model was then generated from the decomposed susceptibility sources by using the following equation:

$${\mathrm{\chi }}_{\text{QSM}}(\mathrm{Y},\mathrm{\nu },{\mathrm{\chi }}_{\text{nb}})=\left[\frac{{\mathrm{\chi }}_{\text{ba}}}{\mathrm{\alpha }}+{\mathrm{\Psi }}_{\text{Hb}}\cdot \mathrm{△}{\mathrm{\chi }}_{\text{Hb}}\cdot \left(-\text{Y}+\frac{1-\left(1-\mathrm{\alpha }\right)\cdot {\text{Y}}_{\text{a}}}{\mathrm{\alpha }}\right)\right]$$

$$\cdot \mathrm{\nu }+\left(1-\frac{\mathrm{\nu }}{\mathrm{\alpha }}\right)\cdot {\mathrm{\chi }}_{\text{nb}}(1)$$

where α is the vein volume fraction assumed to be constant (0.77), Ψ_Hb is the hemoglobin volume fraction (0.0909 for tissue and 0.1197 for vein), and △χ_Hb is the susceptibility difference between deoxyhemoglobin and oxyhemoglobin (12,522 ppb). (b) The qBOLD model of |s_j| based on a voxel spread function method (Cho et al., 2018) was fitted into the following equation:

$$\left|s_j\right|=F_{\text{qBOLD}}(Y,\mathrm{\nu },{\mathrm{\chi }}_{\text{nb}},{\text{s}}^0,R_2,\mathrm{j}\mathrm{△}\text{TE})={\text{s}}^0{\text{e}}^{-R_2\cdot \text{j}\mathrm{△}\text{TE}}$$

$${\text{e}}^{-\mathrm{\nu }\text{f}\left({\mathrm{\delta }}_{\mathrm{\omega }}(\mathrm{Y},{\mathrm{\chi }}_{\text{nb}})\cdot \text{j}\mathrm{△}\text{TE}\right)}g\left(\text{j}\mathrm{△}\text{TE}\right)(2)$$

$${\mathrm{\delta }}_{\mathrm{\omega }}(Y,{\mathrm{\chi }}_{\text{nb}})=\frac{1}{3}\cdot \mathrm{\gamma }\cdot {\text{B}}_0\cdot \left[\text{Hct}\cdot \mathrm{△}{\mathrm{\chi }}_0\cdot \left(1-Y\right)+{\mathrm{\chi }}_{\text{ba}}-{\mathrm{\chi }}_{\text{nb}}\right](3)$$

where γ is the gyromagnetic ratio (267.513 MHz/T), B_0 is the main magnetic field (3.0T), Hct is the hematocrit (0.357), △χ₀is susceptibility difference between fully oxygenated and fully deoxygenated blood [4π × 0.27 ppm (Yablonskiy and Haacke, 1994)], Y is the oxygenation, χ_ba is the purely oxygenated blood susceptibility (−108.3 ppb), χ_nb is the non-blood susceptibility, ν is the vein blood volume fraction, and g accounts for the macroscopic contributions due to voxel sensitivity function. (c) The QSM model of the phase analysis and the qBOLD model of the magnitude analysis were combined by using a denoising regularization R (Cho et al., 2018):

$${\text{argmin}}_{Y,\mathrm{\nu },{\mathrm{\chi }}_{\text{nb}},{\text{s}}^0,{\text{R}}_2}{\sum }_{\text{j}}{|||s_j|-F_{\text{qBOLD}}(Y,\mathrm{\nu },{\mathrm{\chi }}_{\text{nb}},{\text{s}}^0,R_2,\mathrm{j}\mathrm{△}\text{TE})||}_2^2+\text{w}{||\mathrm{\chi }-{\mathrm{\chi }}_{\text{QSM}}(Y,\mathrm{\nu },{\mathrm{\chi }}_{\text{nb}})||}_2^2+\mathrm{\lambda }\text{R}(Y,\mathrm{\nu },{\mathrm{\chi }}_{\text{nb}},{\text{s}}^0,R_2)(4)$$

(d) A CAT method, where voxels with similar mGRE magnitude time evolutions were grouped into a cluster, was applied to improve the effective signal-to-noise ratio (SNR) for QSM + qBOLD model (Cho et al., 2020b).

Histological Examination

All diagnoses were histopathologically proven after surgical resection or tumor biopsy, according to the 2016 WHO classification of CNS tumors, by neuropathologists who were blinded to the MRI data. Genomic sequencing analysis was performed based on tissue availability. IDH1 mutation status was determined by using next-generation sequencing and/or immunohistochemistry (IHC). MGMT promoter methylation status (methylated vs. unmethylated) was detected by using MGMT pyrosequencing from the PyroMark Q24 sequencer. In order to further characterize genetic heterogeneity potentially driving oxygen metabolism or angiogenesis in glioma, we examined the association of established preoperative MRI parameters with coactivated RTKs based on somatic genomic alterations and hallmark copy number variation of glioma as described previously (Supplementary Methods 2; Stommel et al., 2007; Nakada et al., 2020). We also performed IHC analysis of Ki67 proteins routinely analyzed in current clinical practice, using a standard immunohistochemical staining procedure (n = 81).

Statistical Analysis

Normal distribution of the parameters was tested by Shapiro–Wilk test (Razali and Wah, 2011; Das and Rahmatullah Imon, 2016). Data normality was also checked and reported by using Kolmogorov–Smirnov test. Detailed results for testing normality are shown in Supplementary Table 1 and Supplementary Figure 1. In general, the majority of the results of the two methods were consistent for testing the normality. The Shapiro–Wilk test provides greater power than the Kolmogorov–Smirnov test (even with its Lilliefors correction). Additionally, the Shapiro–Wilk test is more appropriate for small sample sizes (N ≤ 50), but it can also be validly applied with large sample sizes. For these reasons, the Shapiro–Wilk test was used for assessing data normality in our study. Then, the subgroups of participants with LGG and GBM, molecular parameters [IDH1 (mutation vs. wild-type), MGMT (methylated vs. unmethylated), activation of RTKs subgroup (RTK vs. non-RTK)] were compared by using Student t-test (for normally distributed data) or Wilcoxon rank-sum test (for non-normally distributed data). There is evidence that the most investigated parameters are non-normally distributed based on the Shapiro–Wilk test (Supplementary Table 1 and Supplementary Figure 1). The DeLong et al. (1988) procedure is commonly used to non-parametrically test the hypothesis of the equality of the areas under the curve (AUCs) and performs well regardless of the variable distribution. Thus, non-parametric receiver operating characteristic (ROC) AUC testing based on the Delong approach was performed to assess the predictive performance for grading and molecular detection. The cutoff values were selected by using the maximized values of the Youden indexes. Then, the sensitivity and specificity at the threshold values for each parameter were determined to evaluate the diagnostic performance. The statistical significance of the single AUC was also calculated to test the null hypothesis that the AUC really equals 0.50. The correlations between MRI biomarkers and the Ki67 expression were evaluated with Pearson correlation. For each outcome, we calculated the false discovery rate using the Benjamini–Hochberg method to account for multiple hypothesis tests (Benjamini and Hochberg, 1995). p < 0.05 was considered statistically significant. Analyses were performed using R software v3.6.1.

Results

Predictability of IDH1 Mutation Status

The representative MRI-based OEF and DCE perfusion images are shown in Figures 1, 2. For oxygen metabolism of all patients with LGG, OEF was found to be remarkably decreased (p = 0.001; Figure 3A) in IDH1 mutation than in IDH1 wild-type lesions. OEF did not enable IDH1 status detection in patients with GBM (p = 0.645; Table 2).

FIGURE 1

Figure 1. MRI-based OEF and DCE perfusion images in representative glioma patients with and without IDH mutation. Two exemplary patients with IDH1-wt (astrocytoma II; A₁–F₁) and IDH1-mut (astrocytoma II; A₂–F₂) lower-grade glioma shown: (A₁) gadolinium contrast-enhanced (GdCE) 3DT1BRAVO, (B₁) T2-FLAIR, (C₁) QQ-based OEF maps by CAT, DCE-MRI with separated K^trans (D₁), CBV (E₁), and CBF (F₁) maps. OEF map (C₁) allows to delineate significantly inhomogeneous hyperintensity in the contrast-enhanced region of the IDH1 wild-type astrocytoma; DCE-derived maps (D₁–F₁) clearly show the area of hyperperfusion and high permeability in the IDH1 wild-type astrocytoma.

FIGURE 2

Figure 2. Representative MRI-based OEF and DCE perfusion images of glioma patients with and without methylated MGMT promoter. Two GBM patients with unmethylated (51-year-old woman with MGMT^–: A₁–F₁) and methylated (58-year-old man MGMT⁺: A₂–F₂) promoter shown: (A₁) GDCE 3D T1BRAVO, (B₁) T2-FLAIR, (C₁) QQ-based OEF maps with CAT, DCE-MRI with separated K^trans (D₁), CBV (E₁), and CBF (F₁) maps. OEF map (C₁) allows to identify a focal area of higher values within the lesion for patients with unmethylated MGMT promoter; DCE-derived maps (D₁–F₁) show a tendency toward higher signal intensities.

FIGURE 3

Figure 3. Quantitative comparison of oxygen metabolism and vascular metrics of genetic profiles. (A–D) Comparison of OEF, K^trans, CBV, and CBF measurements between IDH1 mutant and wild type in lower-grade glioma (left panel), as well as between MGMT methylated and unmethylated glioblastoma (right panel). Comparison of OEF (E) and CBV (F) values between RTK and non-RTK subgroup. The violin plot marked with asterisk (*) indicates that the values were significantly different between the groups (*P < 0.05; **P < 0.01).

TABLE 2

Table 2. Comparison of biomarkers for subgroups of IDH mutation status.

For angiogenesis of all patients with LGG, markedly decreased volume transfer constant (K^trans), CBV, and CBF in IDH1 mutation were found compared with IDH1 wild-type tumors (all p < 0.01; Figures 3B–D). None of the investigated perfusion metrics allowed for differentiation of IDH1 gene mutation status in patients with GBM (p > 0.05).

Predictability of O⁶-Methylguanine-DNA Methyltransferase Promoter Methylation Status

GBM with MGMT promoter methylation showed predominantly decreased OEF compared with unmethylated tumors (p = 0.01; Figure 3A). LGG demonstrated no significant differences in OEF regarding MGMT status (p = 0.408; Table 3).

TABLE 3

Table 3. Comparison of biomarkers for subgroups of MGMT status.

Perfusion analysis showed a tendency toward higher values in patients with unmethylated MGMT promoter but did not allow for significant differentiation of methylated vs. unmethylated gliomas (p > 0.05; Figures 3B–D).

Predictability of Receptor Tyrosine Kinase Aberrations in a Subset of Patients With Next-Generation Sequencing

The association of each imaging parameter and the corresponding RTK aberrations demonstrated a significant positive relationship between RTK subgroup and increased OEF (p = 0.003, Figure 3E and Table 4). Perfusion metrics (CBV and CBF) showed a tendency toward higher values in the RTK subgroup (p = 0.004 and p = 0.011, respectively, Figure 3F).

TABLE 4

Table 4. Comparison of biomarkers for RTK subgroup.

Predictability of World Health Organization Tumor Grade

When comparing imaging biomarkers in tumor cores between WHO grades, OEF demonstrated significant difference between LGG and GBM (p < 0.0001; Table 5). Perfusion analysis in patients with GBM showed areas with increased K^trans, CBV, and CBF when compared with LGG (p < 0.0001 for each).

TABLE 5

Table 5. Comparison of biomarkers for WHO tumor grades.

Diagnostic Performance of Imaging Biomarkers

Regarding classification of molecular profiling for IDH1 mutation status, best test performance was achieved using OEF values with sensitivity and specificity of 74.20 and 86.67%, respectively, at a threshold of 0.196 of all imaging biomarkers in LGG (AUC = 0.828; p < 0.001; Figure 4A). Results of perfusion ROC analysis were good to excellent, albeit with slightly inferior performance (the AUCs of these parameters ranged from 0.781 to 0.815; p < 0.01).

FIGURE 4

Figure 4. Predictability of genetic profiles and WHO tumor grade in glioma patients. ROC analyses performed for OEF and vascular metrics in order to assess the contrast ability to predict (A) IDH mutation status in LGG (IDH₁–R132H mutation vs. wild type), (B) MGMT promoter methylation status in GBM (methylated vs. unmethylated), (C) RTK subgroup in a subset of patients with NGS, and (D) WHO tumor grade (GBM vs. LGG).

Regarding MGMT promoter methylation status classification in patients with GBM, only OEF allowed for significant prediction of MGMT promoter methylation yielding an AUC of 0.784 (p = 0.005; Figure 4B).

In comparison of diagnostic abilities of imaging biomarkers, OEF, and perfusion metrics, best classification of the RTK subgroup was achieved using the OEF values with sensitivity and specificity of 83.33 and 72.41%, respectively, at a cutoff value of 19.28 (AUC = 0.764; p = 0.001). CBV also exhibited moderate performance in identification of RTKs, yielding an AUC of 0.754 (p = 0.002; Figure 4C).

The differentiation ability of LGG and GBM had a range of AUCs from 0.794 to 0.819 (p < 0.0001); reliable prediction of WHO tumor grade was possible using the OEF metric, yielding areas under the ROC curve of 0.810 (p < 0.0001; Figure 4D), with a sensitivity and specificity of 76.74 and 79.17%, respectively. The results of all ROC curve analyses are provided in Supplementary Table 2.

The Mutational Landscape in the Subset of Patients With Next-Generation Sequencing

A more detailed view of molecular biology revealed that the most significantly prominent genes of epidermal growth factor receptor (EGFR) (p = 0.001), PTEN (p = 0.002), and TP53 (p = 0.009) between two groups (RTK vs. non-RTK) separated by the optimized thresholds of imaging parameters were identified (Figures 5A–C). These genes are closely linked to angiogenesis or tumor hypoxia (Amelio and Melino, 2015; Gluck et al., 2015; An et al., 2018). Of note is that the Kyoto Encyclopedia of Genes and Genomes pathway analysis revealed biologically interesting insights mainly enriched in cancer- and angiogenesis-related pathways including glioma, vascular endothelial growth factor (VEGF) signaling, and EGFR tyrosine kinase inhibitor resistance in the RTK subgroup (p < 0.05; Figure 5D).

FIGURE 5

Figure 5. The mutational landscape and biological signatures for 50 tumor samples with NGS, 49 of which harbored genetic alterations. (A–C) Samples from patients have been stratified according to the cutoff value of imaging parameters (CBV, CBF, and OEF) for RTK aberrations. The heat map below displays the somatic mutated genes ordered by their respective mutation frequencies. The significant genes with a p-value < 0.05 were labeled on the left. (D) Kyoto Encyclopedia of Genes and Genomes pathway analysis enriched in the RTK subgroup.

Correlation Between Magnetic Resonance Imaging Measurements With Immunohistochemistry Results

Specifically, histological results in our study also indicated that OEF values (r = 0.42, p = 0.002) or less relevant K^trans values (r = 0.39, p = 0.005) were positively correlated with the cell cycle marker Ki67, but only in MGMT unmethylated gliomas (Figure 6). Apart from that, Ki67 in IDH mutant and/or wild-type gliomas did not correlate with OEF or any perfusion metrics (p > 0.05).

FIGURE 6

Figure 6. Correlation between MRI measures and Ki67 expression in glioma patients. Scatterplot with linear regression line shows a weak to moderate association between MRI measures of OEF (A) and K^trans (B) with Ki67 expression in MGMT promoter methylated and unmethylated gliomas followed by Ki-67 immunohistochemistry on the lower panel (C,D).

Discussion

In this study, we investigated CAT for QQ-based oxygenation extraction fraction mapping; meanwhile, DCE-MRI was explored to demonstrate the non-invasive predictability of genetic profiles and WHO tumor grade in glioma patients. CAT for QQ-based OEF mapping that is respiratory challenge-free is a clinically compatible and potentially valuable imaging technique for oxygenation mapping. The current study revealed the feasibility of CAT for QQ-based OEF mapping, along with characterization of tumor angiogenesis assessed by using DCE-MRI, can potentially be used as non-invasive imaging biomarkers for genetic profile prediction and glioma grading in the clinical routine.

Regional OEF is an essential biomarker for investigating tissue metabolism and function in various diseases such as stroke (Derdeyn et al., 1998, 2002), cerebral tumors (Ito et al., 1982), and Alzheimer’s disease (Ishii et al., 1996). The OEF mapping is achieved via application of CAT to denoise QSM + qBOLD (QQ) estimation of OEF based on 3D multiecho gradient echo acquisition. The CAT algorithm based on unsupervised machine learning substantially improves the robustness against noise of QQ-based mapping of OEF within a few minutes. This CAT for QQ-based OEF mapping eliminates assumptions and is theoretically more accurate than other OEF mapping methods (Kudo et al., 2016; Zhang et al., 2018). Therefore, rapid OEF mapping can be included in routine clinical MRI protocols to accurately evaluate oxygen metabolism in glioma patients and, upon further validation, would help patient stratification for targeted therapy.

QQ-based OEF mapping with CAT performed markedly better in prediction of IDH1 mutation status than WHO tumor grade in this study cohort. Compared with perfusion ROC analyses, the OEF mapping yielded good prediction of IDH mutation status. Therefore, we infer that genetic status, by means of IDH1 mutation status, is more reflective than the histologic class, regarding oxygen metabolism or tumor vascularization. The latest WHO classification of CNS tumors termed gliomas of WHO grades I–III together as “lower-grade gliomas” with a great majority of IDH1 mutations and a wide range of overall survival within this group (Cancer Genome Atlas Research Network et al., 2015; Louis et al., 2016). Mutations in IDH1 gene–encoded enzymes are expected to cause widespread disturbances of cellular metabolism, including DNA hypermethylation with subsequently increased tumor metabolism and degradation of HIF, leading to downstream inhibition of vasculogenesis- and angiogenesis-related signaling (Ye et al., 2013). Hence, the investigated OEF values and vascular metrics are expected to be significantly lower in IDH1-mutated tumors, in agreement with our study findings. It is conceivable that lower OEF and perfusion values in IDH1-mutated gliomas could be due to lower proliferation rates (Bralten et al., 2011) and less angiogenesis (Kickingereder et al., 2015) compared with their wild-type counterparts.

OEF mapping allowed for prediction of MGMT promoter methylation status that is difficult to detect in perfusion MRI. Tumors with unmethylated MGMT promoter showed a trend toward hyperperfusion in our data, in agreement with previous results obtained from perfusion imaging (Choi et al., 2017; Han et al., 2018). MGMT, the gene encoding the DNA repair enzyme O⁶-methylguanine-DNA methyltransferase, is methylated in 30–40% of GBM and 80% of IDH1-mutated LGGs (Hegi et al., 2005). Unmethylated MGMT is associated with high levels of MGMT expression and consequently induces upregulation of HIF-1α, elevated hypoxia, and increased angiogenic potential (Cao et al., 2009; Pistollato et al., 2010; Jue et al., 2017). Transferred to our findings, this may explain higher OEF and perfusion in the GBMs with unmethylated MGMT promoter compared with methylated phenotype. This observation appears consistent with results from Chahal et al. (2010) showing that a direct link between MGMT expression and decreased angiogenicity and tumorigenicity of GBM cells. Most importantly, these distinguishable molecular signatures such as IDH and MGMT can be translated into distinct phenotypes and are non-invasively predictable with physiologic MRI in a clinical setting.

Interestingly, we found moderate associations of increased OEF and perfusion metrics in the designated RTK subgroup, showing a performance of imaging biomarkers to predict relevant RTKs implicated in glioma. Our finding of elevated imaging parameters in tumors may correspond to the promotion of angiogenesis and hypoxia in a subset of patients via amplification and mutational activation of RTK-encoding genes (Stommel et al., 2007). Aberrant activation of RTKs promotes increased angiogenesis through multiple downstream effectors, which mediate a variety of vascular activities, including endothelial cell proliferation, migration, and new vessel formation (Gluck et al., 2015). As to signaling interplays between RTK activity and tumor hypoxia, the activation of VEGF receptor (VEGFR) via HIF-mediated VEGF transcriptional activation and consequent accelerated tumor angiogenesis represents one of the commonly described mechanisms that confer aggressive manifestation of hypoxic tumors (Forsythe et al., 1996). Moreover, activation of the RTK EGFR is a critical pathogenetic event involved in GBMs, with amplification and mutation, which promotes angiogenesis via distinct signaling pathways. In line with the finding of elevated relative CBV in patients with EGFR amplification by Kickingereder et al. (2016), such research efforts are essential to develop targeted therapies for more personalized cancer management.

In the current study, quantitative analysis regarding WHO grade showed increased OEF, vascular permeability, and perfusion metrics in GBM compared to LGG. OEF mapping has a fair prediction ability of WHO tumor grade. The preliminary correlation between HIF-1α expression (as per IHC) and increased OEF areas (Tóth et al., 2013) might well reflect hypoxic areas, highly characterized GBM, and thus aid in the differentiation between LGG and GBM. Our findings are in accordance with generally accepted knowledge about the mechanisms of glioma-associated neovascularization and the role of oxygen metabolism in the following processes (Hardee and Zagzag, 2012). Initially, tumor cells infiltrate through the central nervous system as well as increase oxygen extraction and nutrient supplies through vascular co-option of intact native blood vessels. Involution of co-opted vessels resulted in tumor hypoxia, upregulation of proangiogenic factors, and shift toward an angiogenic phenotype (Holash et al., 1999). Interestingly, higher OEF values were most prevalent in the peripheral area, which also fits with recent observations of Pistollato et al. (2010), who demonstrated expression of hypoxia markers (HIF-1α and VEGF) in this region. Therefore, multiparametric MRI-based physiological imaging may provide a better characterization of tumor microenvironment and yield grading of gliomas in clinical routine.

There are several limitations in the current study. The numbers of patients for the GBM with mutant IDH1 gene were relatively small. Up to now, this study recruited the largest patient cohort with glioma to obtain QQ-based OEF mapping with CAT. IDH1, however, is mutated in only ∼10% of GBMs (Yan et al., 2009). A further limitation of our study is that we did not include a validation of our approach. Biological validation of the MRI-based imaging for OEF and vascular parameters is required by correlation with findings from invasive methods or other imaging modalities in independent study cohorts.

Conclusion

CAT for QQ-based OEF mapping, along with DCE perfusion MRI, enables superior predictability of the molecular characteristics including IDH mutation, MGMT promoter methylation status and RTK subgroup, and fair performance in differentiation of GBM vs. LGG. Impressively, the imaging biomarkers provided macroscopic insights into the heterogeneity of the biological changes in oxygen metabolism or angiogenesis that is potentially driven by genetic aberrations and shed a new light onto potential patient subsets for targeted therapy. Therefore, QQ−based OEF mapping by CAT and tumor-associated angiogenesis assessed by DCE-MRI enables prediction of genetic features and glioma grading, significantly extending the existing repertoire of non-invasive imaging biomarkers used in the preoperative workup and treatment guidance for patients with glioma.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.

Ethics Statement

The studies involving human participants were reviewed and approved by the Tongji Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China. The ethics committee waived the requirement of written informed consent for participation.

Author Contributions

WZ and YW contributed to the conception and design of the study. NS, SZ, and JC carried out the data collection and evaluation. SL and JZ performed the statistical analyses and visualization. NS and YX wrote the manuscript. All authors critically reviewed and approved the manuscript.

Funding

This work was supported in part by the National Natural Science Foundation (No. 81730049).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnins.2021.736891/full#supplementary-material

Abbreviations

OEF, oxygen extraction fraction; CAT, cluster analysis of time evolution; K^trans, volume transfer constant; CBV, cerebral blood volume; CBF, cerebral blood flow; NGS, next-generation sequencing; IDH, isocitrate dehydrogenase; MGMT, O⁶-methylguanine-DNA-methyltransferase; RTK, receptor tyrosine kinase; GBM, glioblastoma; LGG, lower-grade glioma; HIF-1 α, hypoxia-inducible factor 1 α.

Footnotes

1. ^ http://fsl.fmrib.ox.ac.uk/fsl/fslwiki/FSL

References

Amelio, I., and Melino, G. (2015). The p53 family and the hypoxia-inducible factors (HIFs): determinants of cancer progression. Trends Biochem. Sci. 40, 425–434. doi: 10.1016/j.tibs.2015.04.007

PubMed Abstract | CrossRef Full Text | Google Scholar

An, Z., Aksoy, O., Zheng, T., Fan, Q. W., and Weiss, W. A. (2018). Epidermal growth factor receptor and EGFRvIII in glioblastoma: signaling pathways and targeted therapies. Oncogene 37, 1561–1575. doi: 10.1038/s41388-017-0045-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Arevalo-Perez, J., Thomas, A. A., Kaley, T., Lyo, J., Peck, K. K., Holodny, A. I., et al. (2015). T1-weighted dynamic contrast-enhanced MRI as a noninvasive biomarker of epidermal growth factor receptor vIII status. AJNR Am. J. Neuroradiol. 36, 2256–2261. doi: 10.3174/ajnr.A4484

PubMed Abstract | CrossRef Full Text | Google Scholar

Benjamini, Y., and Hochberg, Y. (1995). Controlling the false discovery rate: a practical and powerful approach to multiple testing. J. R. Stat. Soc. Ser. B Methodol. 57, 289–300. doi: 10.1111/j.2517-6161.1995.tb02031.x

CrossRef Full Text | Google Scholar

Bralten, L. B., Kloosterhof, N. K., Balvers, R., Sacchetti, A., Lapre, L., Lamfers, M., et al. (2011). IDH1 R132H decreases proliferation of glioma cell lines in vitro and in vivo. Ann. Neurol. 69, 455–463. doi: 10.1002/ana.22390

PubMed Abstract | CrossRef Full Text | Google Scholar

Bulte, D. P., Kelly, M., Germuska, M., Xie, J., Chappell, M. A., Okell, T. W., et al. (2012). Quantitative measurement of cerebral physiology using respiratory-calibrated MRI. Neuroimage 60, 582–591. doi: 10.1016/j.neuroimage.2011.12.017

PubMed Abstract | CrossRef Full Text | Google Scholar

Cancer Genome Atlas Research Network, Brat, D. J., Verhaak, R. G., Aldape, K. D., Yung, W. K. A., Salama, S. R., et al. (2015). Comprehensive, integrative genomic analysis of diffuse lower-grade gliomas. N. Engl. J. Med. 372, 2481–2498. doi: 10.1056/NEJMoa1402121

PubMed Abstract | CrossRef Full Text | Google Scholar

Cao, V. T., Jung, T. Y., Jung, S., Jin, S. G., Moon, K. S., Kim, I. Y., et al. (2009). The correlation and prognostic significance of MGMT promoter methylation and MGMT protein in glioblastomas. Neurosurgery 65, 866–875; discussion 875. doi: 10.1227/01.NEU.0000357325.90347.A1

CrossRef Full Text | Google Scholar

Chahal, M., Xu, Y., Lesniak, D., Graham, K., Famulski, K., Christensen, J. G., et al. (2010). MGMT modulates glioblastoma angiogenesis and response to the tyrosine kinase inhibitor sunitinib. Neuro Oncol. 12, 822–833. doi: 10.1093/neuonc/noq017

PubMed Abstract | CrossRef Full Text | Google Scholar

Cho, J., Kee, Y., Spincemaille, P., Nguyen, T. D., Zhang, J., Gupta, A., et al. (2018). Cerebral metabolic rate of oxygen (CMRO2) mapping by combining quantitative susceptibility mapping (QSM) and quantitative blood oxygenation level-dependent imaging (qBOLD). Magn. Reson. Med. 80, 1595–1604. doi: 10.1002/mrm.27135

PubMed Abstract | CrossRef Full Text | Google Scholar

Cho, J., Zhang, S., Kee, Y., Spincemaille, P., Nguyen, T. D., Hubertus, S., et al. (2020b). 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. Magn. Reson. Med. 83, 844–857. doi: 10.1002/mrm.27967

PubMed Abstract | CrossRef Full Text | Google Scholar

Cho, J., Lee, J., An, H., Goyal, M. S., Su, Y., and Wang, Y. (2020a). Cerebral oxygen extraction fraction (OEF): comparison of challenge-free gradient echo QSM+qBOLD (QQ) with (15)O PET in healthy adults. J. Cereb. Blood Flow Metab. 41, 1658–1668. doi: 10.1177/0271678X20973951

PubMed Abstract | CrossRef Full Text | Google Scholar

Choi, Y. S., Ahn, S. S., Lee, H. J., Chang, J. H., Kang, S. G., Kim, E. H., et al. (2017). The initial area under the curve derived from dynamic contrast-enhanced MRI improves prognosis prediction in glioblastoma with unmethylated MGMT promoter. AJNR Am. J. Neuroradiol. 38, 1528–1535. doi: 10.3174/ajnr.A5265

PubMed Abstract | CrossRef Full Text | Google Scholar

Das, K. R., and Rahmatullah Imon, A. H. M. (2016). A brief review of tests for normality. Am. J. Theor. Appl. Stat. 5, 5–12. doi: 10.11648/j.ajtas.20160501.12

CrossRef Full Text | Google Scholar

DeLong, E. R., DeLong, D. M., and Clarke-Pearson, D. L. (1988). Comparing the areas under two or more correlated receiver operating characteristic curves: a nonparametric approach. Biometrics 44, 837–845. doi: 10.2307/2531595

CrossRef Full Text | Google Scholar

Derdeyn, C. P., Videen, T. O., Yundt, K. D., Fritsch, S. M., Carpenter, D. A., Grubb, R. L., et al. (2002). Variability of cerebral blood volume and oxygen extraction: stages of cerebral haemodynamic impairment revisited. Brain 125(Pt 3), 595–607. doi: 10.1093/brain/awf047

PubMed Abstract | CrossRef Full Text | Google Scholar

Derdeyn, C. P., Yundt, K. D., Videen, T. O., Carpenter, D. A., Grubb, R. L. Jr., and Powers, W. J. (1998). Increased oxygen extraction fraction is associated with prior ischemic events in patients with carotid occlusion. Stroke 29, 754–758. doi: 10.1161/01.STR.29.4.754

CrossRef Full Text | Google Scholar

Forsythe, J. A., Jiang, B. H., Iyer, N. V., Agani, F., Leung, S. W., Koos, R. D., et al. (1996). Activation of vascular endothelial growth factor gene transcription by hypoxia-inducible factor 1. Mol. Cell. Biol. 16, 4604–4613. doi: 10.1128/MCB.16.9.4604

PubMed Abstract | CrossRef Full Text | Google Scholar

Gluck, A. A., Aebersold, D. M., Zimmer, Y., and Medova, M. (2015). Interplay between receptor tyrosine kinases and hypoxia signaling in cancer. Int. J. Biochem. Cell Biol. 62, 101–114. doi: 10.1016/j.biocel.2015.02.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Han, Y., Yan, L. F., Wang, X. B., Sun, Y. Z., Zhang, X., Liu, Z. C., et al. (2018). Structural and advanced imaging in predicting MGMT promoter methylation of primary glioblastoma: a region of interest based analysis. BMC Cancer 18:215. doi: 10.1186/s12885-018-4114-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Hanahan, D., and Weinberg, R. A. (2011). Hallmarks of cancer: the next generation. Cell 144, 646–674. doi: 10.1016/j.cell.2011.02.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Hardee, M. E., and Zagzag, D. (2012). Mechanisms of glioma-associated neovascularization. Am. J. Pathol. 181, 1126–1141. doi: 10.1016/j.ajpath.2012.06.030

PubMed Abstract | CrossRef Full Text | Google Scholar

Hegi, M. E., Diserens, A. C., Gorlia, T., Hamou, M. F., de Tribolet, N., Weller, M., et al. (2005). MGMT gene silencing and benefit from temozolomide in glioblastoma. N. Engl. J. Med. 352, 997–1003. doi: 10.1056/NEJMoa043331

PubMed Abstract | CrossRef Full Text | Google Scholar

Holash, J., Maisonpierre, P. C., Compton, D., Boland, P., Alexander, C. R., Zagzag, D., et al. (1999). Vessel cooption, regression, and growth in tumors mediated by angiopoietins and VEGF. Science 284, 1994–1998. doi: 10.1126/science.284.5422.1994

PubMed Abstract | CrossRef Full Text | Google Scholar

Ishii, K., Kitagaki, H., Kono, M., and Mori, E. (1996). Decreased medial temporal oxygen metabolism in Alzheimer’s disease shown by PET. J. Nucl. Med. 37, 1159–1165.

Google Scholar

Ito, M., Lammertsma, A. A., Wise, R. J., Bernardi, S., Frackowiak, R. S., Heather, J. D., et al. (1982). Measurement of regional cerebral blood flow and oxygen utilisation in patients with cerebral tumours using 15O and positron emission tomography: analytical techniques and preliminary results. Neuroradiology 23, 63–74. doi: 10.1007/BF00367239

PubMed Abstract | CrossRef Full Text | Google Scholar

Jensen, R. L., Mumert, M. L., Gillespie, D. L., Kinney, A. Y., Schabel, M. C., and Salzman, K. L. (2014). Preoperative dynamic contrast-enhanced MRI correlates with molecular markers of hypoxia and vascularity in specific areas of intratumoral microenvironment and is predictive of patient outcome. Neuro Oncol. 16, 280–291. doi: 10.1093/neuonc/not148

PubMed Abstract | CrossRef Full Text | Google Scholar

Jue, T. R., Nozue, K., Lester, A. J., Joshi, S., Schroder, L. B., Whittaker, S. P., et al. (2017). Veliparib in combination with radiotherapy for the treatment of MGMT unmethylated glioblastoma. J. Transl. Med. 15:61. doi: 10.1186/s12967-017-1164-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Kaur, B., Khwaja, F. W., Severson, E. A., Matheny, S. L., Brat, D. J., and Van Meir, E. G. (2005). Hypoxia and the hypoxia-inducible-factor pathway in glioma growth and angiogenesis. Neuro Oncol. 7, 134–153. doi: 10.1215/S1152851704001115

PubMed Abstract | CrossRef Full Text | Google Scholar

Kickingereder, P., Bonekamp, D., Nowosielski, M., Kratz, A., Sill, M., Burth, S., et al. (2016). Radiogenomics of glioblastoma: machine learning-based classification of molecular characteristics by using multiparametric and multiregional MR imaging features. Radiology 281, 907–918. doi: 10.1148/radiol.2016161382

PubMed Abstract | CrossRef Full Text | Google Scholar

Kickingereder, P., Sahm, F., Radbruch, A., Wick, W., Heiland, S., van Deimling, A., et al. (2015). IDH mutation status is associated with a distinct hypoxia/angiogenesis transcriptome signature which is non-invasively predictable with rCBV imaging in human glioma. Sci. Rep. 5:16238. doi: 10.1038/srep16238

PubMed Abstract | CrossRef Full Text | Google Scholar

Kudo, K., Liu, T., Murakami, T., Goodwin, J., Uwano, I., Yamashita, F., et al. (2016). Oxygen extraction fraction measurement using quantitative susceptibility mapping: comparison with positron emission tomography. J. Cereb. Blood Flow Metab. 36, 1424–1433. doi: 10.1177/0271678X15606713

PubMed Abstract | CrossRef Full Text | Google Scholar

la Fougere, C., Suchorska, B., Bartenstein, P., Kreth, F. W., and Tonn, J. C. (2011). Molecular imaging of gliomas with PET: opportunities and limitations. Neuro Oncol. 13, 806–819. doi: 10.1093/neuonc/nor054

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, Z., Spincemaille, P., Yao, Y., Zhang, Y., and Wang, Y. (2018). MEDI+0: morphology enabled dipole inversion with automatic uniform cerebrospinal fluid zero reference for quantitative susceptibility mapping. Magn. Reson. Med. 79, 2795–2803. doi: 10.1002/mrm.26946

PubMed Abstract | CrossRef Full Text | Google Scholar

Louis, D. N., Perry, A., Reifenberger, G., von Deimling, A., Figarella-Branger, D., Cavenee, W. K., et al. (2016). The 2016 World Health Organization classification of tumors of the central nervous system: a summary. Acta Neuropathol. 131, 803–820. doi: 10.1007/s00401-016-1545-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Nakada, M., Kita, D., Teng, L., Pyko, I. V., Watanabe, T., Hayashi, Y., et al. (2020). Receptor tyrosine kinases: principles and functions in glioma invasion. Adv. Exp. Med. Biol. 1202, 151–178. doi: 10.1007/978-3-030-30651-9_8

CrossRef Full Text | Google Scholar

Paech, D., Nagel, A. M., Schultheiss, M. N., Umathum, R., Regnery, S., Scherer, M., et al. (2020). Quantitative dynamic oxygen 17 MRI at 7.0 T for the cerebral oxygen metabolism in glioma. Radiology 295, 181–189. doi: 10.1148/radiol.2020191711

PubMed Abstract | CrossRef Full Text | Google Scholar

Pistollato, F., Abbadi, S., Rampazzo, E., Persano, L., Della Puppa, A., Frasson, C., et al. (2010). Intratumoral hypoxic gradient drives stem cells distribution and MGMT expression in glioblastoma. Stem Cells 28, 851–862. doi: 10.1002/stem.415

PubMed Abstract | CrossRef Full Text | Google Scholar

Razali, N. M., and Wah, Y. B. (2011). Power comparisons of Shapiro-Wilk, Kolmogorov-Smirnov, Lilliefors and Anderson-Darling tests. J. Stat. Model. Anal. 2, 21–33.

Google Scholar

Semenza, G. L. (2003). Targeting HIF-1 for cancer therapy. Nat. Rev. Cancer 3, 721–732. doi: 10.1038/nrc1187

PubMed Abstract | CrossRef Full Text | Google Scholar

Stommel, J. M., Kimmelman, A. C., Ying, H., Nabioullin, R., Ponugoti, A. H., Wiedemeyer, R., et al. (2007). Coactivation of receptor tyrosine kinases affects the response of tumor cells to targeted therapies. Science 318, 287–290. doi: 10.1126/science.1142946

PubMed Abstract | CrossRef Full Text | Google Scholar

Tóth, V., Förschler, A., Hirsch, N. M., den Hollander, J., Kooijman, H., Gempt, J., et al. (2013). MR-based hypoxia measures in human glioma. J. Neurooncol. 115, 197–207. doi: 10.1007/s11060-013-1210-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Yablonskiy, D. A., and Haacke, E. M. (1994). Theory of NMR signal behavior in magnetically inhomogeneous tissues: the static dephasing regime. Magn. Reson. Med. 32, 749–763. doi: 10.1002/mrm.1910320610

PubMed Abstract | CrossRef Full Text | Google Scholar

Yan, H., Parsons, D. W., Jin, G., McLendon, R., Rasheed, B. A., Yuan, W., et al. (2009). IDH1 and IDH2 mutations in gliomas. N. Engl. J. Med. 360, 765–773. doi: 10.1056/NEJMoa0808710

PubMed Abstract | CrossRef Full Text | Google Scholar

Ye, D., Ma, S. H., Xiong, Y., and Guan, K. L. (2013). R-2-hydroxyglutarate as the key effector of IDH mutations promoting oncogenesis. Cancer Cell 23, 274–276. doi: 10.1016/j.ccr.2013.03.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, J., Cho, J., Zhou, D., Nguyen, T. D., Spincemaille, P., Gupta, A., et al. (2018). Quantitative susceptibility mapping-based cerebral metabolic rate of oxygen mapping with minimum local variance. Magn. Reson. Med. 79, 172–179. doi: 10.1002/mrm.26657

PubMed Abstract | CrossRef Full Text | Google Scholar