- 1Division of Radiology, Department of Clinical Sciences, Skåne University Hospital, Lund University, Lund, Sweden
- 2Division of Pathology, Department of Clinical Sciences, Skåne University Hospital, Lund University, Lund, Sweden
- 3Department of Medical Imaging and Physiology, Skåne University Hospital, Lund, Sweden
- 4Department of Clinical Sciences, Division of Oncology, Lund University, Lund, Sweden
- 5Division of Neurology, Department of Clinical Sciences, Skåne University Hospital, Lund University, Lund, Sweden
- 6Kamprad Laboratory, Division of Neurosurgery, Department of Clinical Sciences, Skåne University Hospital, Lund University, Lund, Sweden
- 7F. M. Kirby Research Center for Functional Brain Imaging, Kennedy Krieger Institute, Baltimore, MD, United States
- 8Department of Neurology, Johns Hopkins University School of Medicine, Baltimore, MD, United States
- 9Department of Medical Radiation Physics, Lund University, Lund, Sweden
- 10Lund University Bioimaging Center (LBIC), Lund University, Lund, Sweden
Background/objectives: High grade gliomas (HGG) are aggressive brain tumors, most frequently glioblastoma and astrocytoma grade 4. Methylation of O6-methylguanine-DNA methyltransferase (MGMT) promoter in HGG is crucial for temozolomide efficacy. As MGMT promoter methylation (MGMTpm) assessment requires tumor tissue, magnetic resonance imaging (MRI) is of interest for non-invasive prediction. We aimed to analyze volumetric data from edema, contrast-enhancing tumor, necrosis, total-tumor and total-tumor/edema ratio for MGMTpm prediction in HGG. Further we assessed overall survival (OS) and progression free survival (PFS) between groups and volumes.
Methods: Segmentation was performed using deep learning models (DL-models), DeepBraTumIA and Raidionics, on 70 HGG patients (45 males, 32 MGMTpm). Manual segmentation was conducted in 37 for validation of DL-models. Group differences were evaluated using Man-Whitney U tests and receiver operation characteristic (ROC) curves. Multivariate analysis was conducted using logistic regression and bootstrapping. Dice coefficient, intraclass correlation coefficient (ICC) and Kruskal–Wallis test evaluated DL-model performance.
Results: MGMTpm tumors displayed significantly larger edema, segmented by DeepBraTumIA (p = 0.03), and lower total-tumor/edema ratio segmented by both DL-models (p < 0.01). Raidionics segmented total-tumor/edema ratio showed highest univariate predictive ability with area under curve 0.687 (sensitivity 46.2%, specificity 87.5%). Multivariate analysis confirmed this, showing that the ratios from both DL-models were the only ROIs to remain independent, significant predictors (p < 0.05) after controlling for clinical covariates. The overall multivariate models were significant (p = 0.01) and improved prediction over baseline. ICC showed interclass correlation of 0.96 (contrast-enhancing tumor), 0.50 (tumor necrosis) and 0.90 (peritumoral edema). Segmentation methods demonstrated 83–91% median overlap in contrast-enhancing tumor, 67–80% in necrosis and 80–84% in edema regions. Significant OS and PFS differences were observed, notably being longer in MGMTpm tumors and lower tumor volumes.
Conclusion: This study suggests that significant radiological differences in MGMTpm can be found using deep learning models, primarily in tumor edema volume. MGMTpm status and region of interest volumes impact OS and PFS. Future studies should incorporate other molecular imaging sequences for methylation prediction.
1 Introduction
Grade 4 glioblastoma (GBM) and astrocytoma grade 4 are highly aggressive tumors, both being high grade gliomas (HGG) (Ostrom et al., 2018; Louis et al., 2021). The distinction between these two gliomas is based on isocitrate dehydrogenase (IDH) mutation status (Torp et al., 2022). Astrocytoma expresses an IDH mutation, while GBM is IDH-wildtype (Torp et al., 2022). Despite modern treatment, both tumor entities have poor survival rate (Yan et al., 2009).
Standard of care for newly diagnosed GBM and astrocytoma consists of surgical resection followed by radiotherapy (RT) with concomitant temozolomide (TMZ), followed by adjuvant TMZ according to the Stupp protocol (Stupp et al., 2005; Witthayanuwat et al., 2018; Fernandes et al., 2017). However, the combined treatment regimen has shown considerable variances in tumor response, especially regarding TMZ (Hegi et al., 2005; Hegi et al., 2024).
It has been established that this is due to the methylation of the O6-methylguanine-DNA methyltransferase (MGMT) promoter (Hegi et al., 2005; Hegi et al., 2024; Singh et al., 2021). MGMT plays a crucial role in modulating cell survival through repairing DNA damage by removing alkyl groups from the O6 position of guanine (Singh et al., 2021). Since TMZ exerts its cytotoxic effect by adding an alkyl group to this position (Hegi et al., 2024), an active MGMT enzyme can counteract TMZ, rendering it ineffective (Hegi et al., 2024). Methylation of the MGMT-promoter limits MGMT production within the cell (Esteller and Herman, 2004).
Testing for MGMT methylation is done through DNA sequencing methods such as pyrosequencing or methylation specific polymerase chain reaction (PCR) on tumor tissue (Hegi et al., 2005; Wick et al., 2012). Furthermore, MGMT testing is rarely done in HGG recurrence, despite tumor heterogeneity and potential variance in MGMT methylation (Weller et al., 2013; Mansouri et al., 2019). Thus, there is a need for noninvasive, dynamic testing of MGMT promoter methylation in newly diagnosed HGG. Magnetic resonance imaging (MRI) provides a noninvasive approach that accounts for tumor heterogeneity (Suh et al., 2018).
There are some previous studies that, with different methods, have explored the use of MRI for predicting MGMT promoter methylation in HGG. A meta-analysis found that MGMT-methylated tumors typically exhibit less aggressive MRI features on conventional MRI, such as reduced peritumoral edema and tumor burden, compared to non-methylated (Suh et al., 2018). Other studies also found predictive potential in peritumoral edema volume (Ellingson et al., 2012; Kanas et al., 2017). Additionally, other studies found no predictive ability of MGMT promoter methylation in peritumoral edema (Carrillo et al., 2012; Kim et al., 2022), but it is concluded that the edema volume could impact survival if the tumor is MGMT methylated (Carrillo et al., 2012).
It is valuable for the neurosurgeon and the treating physician to have pre-operative knowledge about the MGMT status as it might influence pre-operative decision making. Therefore, using MRI for volumetric tumor characterization might aid in this distinction and could potentially be more effectively done using deep learning (DL) based models (Bonada et al., 2024).
The aim of this study is to evaluate imaging features on conventional MRI for MGMT promoter methylation such as tumor characteristics; tumor volume, tumor necrosis volume, and edema volume by using two different DL models, DeepBraTumIA (DeepBraTumIA, 2025) and Raidionics (Bouget et al., 2023; GitHub, 2023), to assess the DL-models ability to prediction of MGMT methylation using raw segmentation outputs. Thus, the novelty of our work lies in the real-world application of readily available DL-models with assessable user interfaces in predicting MGMTpm methylation. We will also assess overall survival (OS) and progression free survival (PFS) of MGMT promotor methylated vs. non MGMT promoter methylated tumors and higher vs. lower volumes of regions of interest (ROIs).
2 Materials and methods
2.1 Patient recruitment
In this retrospective study 70 patients diagnosed with HGG [GBM grade 4 and astrocytoma grade 4 according to the WHO 2021 classification (Louis et al., 2021)] were included from a larger cohort of an ongoing project at the Department of Clinical Sciences, Division of Diagnostic Radiology, Lund University. The study was conducted in accordance with the Declaration of Helsinki, and approved by the Regional Ethics Board Lund, Sweden and by the Swedish Ethical Review Authority prior to starting any work. Informed consent was obtained from all subjects involved in the study. The inclusion criteria were as follows: >18 years old, undergone biopsy or resection of the brain lesion for pathological analysis, and completion of preoperative MRI examination, including conventional MRI sequences. Exclusion criteria were having a pacemaker or metallic wires in the body incompatible with MRI and being unable to sign informed consent.
2.2 MGMT analyses
MGMT status was routinely assessed using pyrosequencing, performed at the Division of Pathology, Skåne University Hospital, Lund, Sweden (Mikeska et al., 2007). A cut-off ≥10% was defined as MGMT promoter methylation.
2.3 Image acquisition
Thirty-nine patients were examined on a MAGNETOM Prisma 3T MRI scanner (Siemens Healthcare, Erlangen, Germany) equipped with a 20-channel head coil. Whole-brain imaging sequences included: T2 Turbo spin echo with a repetition time (TR) = 6,000 ms, echo time (TE) = 100 ms, in-plane resolution = 1 × 1 mm2, slice thickness = 5 mm; T2-weighted Fluid-Attenuated Inversion Recovery (T2-FLAIR) with a TR = 5,000 ms, inversion time (TI) = 1800 ms, TE = 393 ms, in-plane resolution = 1 × 1 mm2, slice thickness = 1 mm, and pre- and post-gadolinium contrast T1 (T1-GD) magnetization prepared rapid gradient echo (MPRAGE) with a TR = 1900 ms, TI = 900 ms, TE = 2.54 ms, in-plane resolution = 1 × 1 mm2, slice thickness = 1 mm.
Thirty-one patients were examined on a MAGNETOM Skyra 3T system (Siemens Healthcare, Erlangen, Germany). The MRI protocol consisted of morphological imaging; T2 Turbo spin echo (TE/TR = 100 ms/6,870 ms), T2-FLAIR (TE/TR/TI = 81 ms/9,000 ms/2,500 ms) and T1-MPRAGE (TE/TR/TI = 2.54 ms/1,900 ms/900 ms, 1 mm isotropic voxels).
2.4 Measurements
Three separate primary segmentations (manual- and two DL-model segmentations) were performed on the tumors to determine peritumoral edema, tumor necrosis, and contrast-enhancing tumor volumes. As part of the processing pipeline in DeepBraTumIA (DeepBraTumIA, 2025), the images were co-registered to the MNI152 atlas with a rigid transform, 1 mm isotropic voxel size and matrix size 182 × 218 × 182, using the T1-GD scan as reference. Skull-stripping was then performed using the HD Brain Extraction Tool (HD-BET) (Isensee et al., 2019). Manual segmentations were conducted on the registered, skull-stripped images while Raidionics (Bouget et al., 2023) was done on registered images. Voxel-based volume measurements were performed in 3D Slicer (slicer.org, version 5.6.2) using the Quantification of Segmentation module, where one isotropic voxel within the ROI corresponded to one cubic millimeter (Fedorov et al., 2012). Dice-coefficients were calculated using segment comparison module from SlicerRT (version 7168e01) extension in 3D slicer (Pinter et al., 2012). A visual workflow summarization from MRI acquisition to DL-segmentation and final analysis can be seen in Figure 1.
Figure 1. Workflow diagram summarizing the study pipeline, from patient inclusion (n = 70) through image acquisition, pre-processing, segmentation, and final analysis.
2.5 DL-segmentation
Both DeepBraTumIA and Raidionics are comprehensive, ready-to-use software packages that require no additional training or preprocessing on the end-user’s data (DeepBraTumIA, 2025; Bouget et al., 2023). DeepBraTumIA employs a modified U-net architecture, so called no-new-netU-net (nnU-net) for all segmentations tasks (DeepBraTumIA, 2025). NnU-Nets have shown promising results in other segmentation tasks such as post-operative segmentation of GBM (Cepeda et al., 2024). In contrast, Raidionics utilizes an attention-gated U-net (AGU-net) architecture (Bouget et al., 2023). Both models are trained on pre- and post-operative brain tumors (DeepBraTumIA, 2025; Bouget et al., 2023). However, it should be noted that no quality assessment or post processing of individual segmentation from the DL-models was conducted in this study.
2.6 Manual segmentation
For manual segmentation, registered and skull-stripped images were imported into 3D Slicer. Manual segmentation was performed on 37 out of 70 subjects to serve as a validation benchmark for the DL models. Using T1-GD, a ROI for the contrast-enhancing tumor was identified as the hyperintense region relative to the surrounding tissue, excluding any vessels. The tumor necrosis ROI was measured using T1-GD and defined as a hypointense area in T1-GD within the tumor ROI. Lastly, the edema ROI was identified as the hyperintense area on FLAIR images excluding the contrast-enhanced area and necrosis identified using the T1-GD. Manual segmentation was performed by a junior researcher (EZ), under the guidance of a junior consultant (TS) and a board certified neuroradiologist (PS). Individual segmentations were sampled and quality assessed by TS and PCS. During the segmentation process, the author (EZ) was blinded to MGMT-status and to the results of the other segmentation models.
2.7 Statistical analysis
All statistical analyses were performed using IBM SPSS (version 29.0.2.0). The ROIs analyzed were tumor volume, tumor necrosis volume, peritumoral edema volume, total tumor volume consisting of the sum of tumor volume, tumor necrosis volume, and the total-tumor/edema ratio.
Differences in betamethasone (corticosteroid) usage between MGMT groups were assessed using a Mann–Whitney U test while remaining patient characteristics were compared using chi-square tests. To assess differences in ROIs between non-methylated and methylated MGMT promoter tumors the Mann–Whitney U test was applied. ROIs from DeepBraTumIA and Raidionics were separately tested. The ROIs’ volumes predictive ability were assessed using receiver operating characteristic (ROC) curves, also calculating the Youden optimal cutoff. For multivariate testing a logistic regression with bootstrapping (1,000 samples) was performed with an individual ROI, age, sex, tumor type, corticosteroid dose and IDH as covariates. Logistic regression was also performed excluding IDH status due to small sample size causing large confidence intervals, model instability and lack of confidence in result interpretation.
For comparisons between the volumetric data from volumes of manual ROIs and those generated by DL models, a Wilcoxon signed ranks test was used. Statistical significance was defined as p < 0.05. Furthermore, inter-observer reliability among segmentation methods was evaluated using intraclass correlation coefficient (ICC).
Survival analysis of OS and PFS was performed using the Kaplan–Meier method and log-rank test, comparing MGMT-promoter methylated and unmethylated tumors. For each ROIs’ volume the same approach was used to compare OS and PFS between tumors with volumes above the median, defined as higher volume, and those with volumes at or below the median, defined as lower volume.
3 Results
3.1 Patient characteristics
A total of 70 patients (45 males, 25 females; mean age: 60 years) with histologically confirmed HGG (60 GBM and 10 astrocytoma WHO grade 4) and known MGMT promoter methylation status were included. Among these, 32 tumors had a methylated MGMT promoter, while the remaining 38 were unmethylated. Five of the tumors did not have histologically determined IDH status due to lacking sufficient material for the analysis. These tumors were classified as GBM, using the histological morphology and other pathology markers, by a senior board-certified pathologist (XS). For patient demographics and tumor characteristics see Table 1.
Table 1. Patient characteristics, including demographic, histological tumor classification and tumor localization.
Significant differences in tumor type between the MGMT-methylated and non-methylated groups were observed (p = 0.02). The non-methylated group had a higher proportion of GBM [36 (95%) vs. 24 (75%)], while the methylated group had a comparatively higher proportion of grade 4 astrocytomas [8 (25%) vs. 2 (5%)] (Table 1). The remaining characteristics in Table 1 demonstrated no significant differences (data not shown).
3.2 Segmentation assessment
Comparing DeepBraTumIA to manual segmentation, significant differences were found in mean ranks in four out of five ROIs (tumor necrosis, total tumor volume, peritumoral edema and ratio) (Table 2). Compared to volumes from manual segmentation, both mean tumor necrosis and total tumor volume were significantly over segmented by 7,946 mm3 (p < 0.01) and 8,600 mm3 (p < 0.01) respectively (Table 2). Mean peritumoral edema was under segmented by 6,118 mm3 (p = 0.02) and ratio was overestimated (mean +0.144, p < 0.01) (Table 2). Comparing dice-coefficients across all 37 subjects shows a median overlap of 85.0% in contrast-enhancing tumor, 66.7% in tumor necrosis and 83.9% in peritumoral edema (Figure 2).
Table 2. Mean volumes of regions of interest (ROIs) from manually segmented images of 37 patients and deep learning (DL) models segmentation in the same images.
Figure 2. Box plot showing dice coefficient between segmentation methods within region of interest. Dice coefficients to manual from manually segmented images of 37 patients. Dice coefficients between deep learning (DL)-models are based on all 70 patients.
Comparing Raidionics to manual segmentation, significant overestimation of total tumor volume by 1,494 mm3 (p = 0.01) was observed (Table 2). No significant differences were observed for the other ROIs. Comparing dice-coefficients across all 37 subjects shows a median overlap of 83.0% in contrast-enhancing tumor, 73.1% in necrosis and 79.2% in edema (Figure 2).
Comparing segmentations between the two DL-models across all 70 subjects shows a median overlap of 90.7% in contrast-enhancing tumor, 80.3% in necrosis and 82.4% in edema (Figure 2).
Assessing all three segmentation methods using ICC showed interclass correlation of 0.96 (95% CI 0.94–0.98, p < 0.001) for the contrast-enhancing tumor, 0.50 (95% CI 0.14–0.73, p = 0.006) for tumor necrosis, and 0.90 (95% CI, 0.83–0.95, p < 0.001) for peritumoral edema. These results indicate excellent, poor to moderate and good to excellent reliability, respectively. Figure 3 shows a visual assessment of inter-observer agreement among the different segmentation methods. Two assessments are shown, one demonstrates good overlap in all segmented ROIs, and one reveals a poor overlap through mislabeling of the tumor necrosis ROI.
Figure 3. Depiction of the marked region of interest (ROIs) on the same slice in T1-GD and T2 images between manual, DeepBraTumIA and Raidionics segmentation. Pink is contrast enhancing tumor ROI, green is tumor necrosis ROI and light blue is peritumoral edema ROI. The image shows good and poor segmentation overlap. The good segmentation overlap showcases high agreement across all three measured ROIs and segmentation methods. The poor segmentation overlap showcases a case where DeepBraTumIA clearly mislabels an area with edema as necrosis.
3.3 MGMT-prediction
3.3.1 Mixed cohort
Significant differences between non methylated and methylated MGMT tumors were observed in the total tumor/edema ratio ROI in both DL-models (DeepBraTumIA p = 0.02, Raidionics p < 0.01), and in the edema volumes extracted from DeepBraTumIA (p < 0.05) (Table 3). Peritumoral edema was significantly larger in the methylated group (mean 81,934 mm3, Std 8,733 mm3) compared to non-methylated (mean 61,249 mm3, Std 6,471 mm3) (Table 3). The total tumor/edema ratio was significantly lower in MGMT methylated group from both DeepBraTumIA (0.417, Std 0.062) and Raidionics (0.332, Std 0.048) segmentations compared to non MGMT methylated [DeepBraTumIA (0.876, Std 0.213) and Raidionics (0.747, Std 0.120)] (Table 3).
Table 3. Mean volumes from each region of interest (ROI) segmented by both deep learning models separated by O6-methylguanine-DNA methyltransferase (MGMT) methylation status.
ROC analysis indicated potential MGMT promotor methylation prediction in the peritumoral edema ROI from DeepBraTumIA (AUC 0.638, 95% CI 0.504–0.772, sensitivity 68.8%, specificity 61.5%, at Youden index cutoff 62,633 mm3) (Table 4). Visualization of the ROC-curve for DeepBraTumIA segmented peritumoral edema can be seen in Figure 4A. The total tumor/edema ratio was predictive of non MGMT promotor methylation from DeepBraTumIA (AUC 0.667, 95% CI 0.540–0.794, sensitivity 60.5%, specificity 68.7%, at Youden optimal cutoff 0.449) and Raidionics (AUC 0.687, 95% CI 0.562–0.811, sensitivity 44.7%, specificity 87.5%, at Youden optimal cutoff 0.625) (Table 4). Visualization of the ROC-curve for DeepBraTumIA and Raidionics based total-tumor/edema ratio can be seen in Figure 4B.
Table 4. Area under curve (AUC) and Youden’s optimal cut off for all regions of interest (ROI’s) with significant results from receiver operation characteristics (ROC) analysis.
Figure 4. (A) Receiver operation characteristics (ROC) curve for O6-methylguanine-DNA methyltransferase (MGMT)-promoter methylation prediction using peritumoral edema region of interest segmented using DeepBraTumIA. AUC 0.638, 95% CI 0.504–0.772. (B) Receiver operation characteristics (ROC) curve for non O6-methylguanine-DNA methyltransferase (MGMT)-promoter methylation prediction using total-tumor/edema ratio segmented by DeepBraTumIA and Raidionics. DeepBraTumIA AUC 0.667, 95% CI 0.540–0.794, Raidionics AUC 0.687, 95% CI 0.562–0.811.
Multivariate analysis was only performed on ROIs which were found to be significantly different in volume between MGMT methylated and non-methylated tumors. Initial modeling was performed using IDH as a covariate, however this proved to cause an unstable model, with large confidence intervals and skewed significance results as seen in Table 5. However, with and without IDH we observed similar results regarding ROI significance and effects of remaining covariates. When excluding IDH, two out of the three tested ROIs, resulted in overall statistically significant MGMT prediction model, those being ratio segmented by DeepBraTumIA (prediction improved to 67.1% from 54.3% (baseline), Omnibus test of model coefficients p = 0.01) and Raidionics (prediction improved to 65.7% from 54.3% (baseline), Omnibus test of model coefficients p = 0.01). More importantly, the ROIs were the only variables to stand out as independent, statistically significant predictors post bootstrapping with DeepBraTumIA ratio p = 0.03 (B = −2.34, 95% CI −7.63 to −1.03) and Raidionics ratio p = 0.02 (B = −2.09, 95% CI −6.49 to −1.02). Remaining co-variates were found to be non-significant (Table 5).
Table 5. Results of multivariate logistic regression with bootstrapping (1,000 samples) for prediction of O6-methylguanine-DNA methyltransferase (MGMT) methylation in high grade gliomas showing separate models for DeepBraTumIA ratio, DeepBraTum peritumoral edema and Raidionics ratio.
3.3.2 GBM-subgroup
Significant differences between non methylated and methylated MGMT GBMs were observed in the total tumor/edema ratio ROI in both DL-models (DeepBraTumIA p < 0.01, Raidionics p < 0.01). The total tumor/edema ratio was significantly lower in MGMT methylated group from both DeepBraTumIA (0.338, Std 0.045) and Raidionics (0.335, Std 0.052) segmentations compared to non MGMT methylated [DeepBraTumIA (0.891, Std 0.225) and Raidionics (0.770, Std 0.125)].
ROC analysis showed best predictive ability of MGMT promotor methylation in peritumoral edema ROI from DeepBraTumIA (0.616, 95% CI 0.467–0.764, sensitivity 72.0%, specificity 55.6% at Youden optimal cutoff 53,076 mm3). Visualization of the ROC-curve for DeepBraTumIA segmented peritumoral edema can be seen in Figure 5A. Total tumor/edema ratio was predictive of non MGMT promotor methylation from DeepBraTumIA (AUC 0.711, 95% CI 0.582–0.840, sensitivity 44.4%, specificity 92.0%, at Youden optimal cutoff 0.625) and Raidionics (AUC 0.698, 95% CI 0.566–0.829, sensitivity 44.4%, specificity 88.0%, at Youden optimal cutoff 0.638). Visualization of the ROC-curve for DeepBraTumIA and Raidionics based total-tumor/edema ratio can be seen in Figure 5B.
Figure 5. (A) Receiver operation characteristics (ROC) curve for O6-methylguanine-DNA methyltransferase (MGMT)-promoter methylation prediction in GBM subgroup using peritumoral edema region of interest segmented using DeepBraTumIA. 0.616, 95% CI 0.467–0.764. (B) Receiver operation characteristics (ROC) curve for non O6-methylguanine-DNA methyltransferase (MGMT)-promoter methylation prediction in GBM subgroup using total-tumor/edema ratio segmented by DeepBraTumIA and Raidionics. DeepBraTumIA AUC 0.711, 95% CI 0.582–0.840, Raidionics AUC 0.698, 95% CI 0.566–0.829.
3.4 Survival analysis
3.4.1 Mixed cohort
OS and PFS were shown to be significantly higher (p < 0.01 and p = 0.03 respectively) in MGMT-promoter methylated compared to non-methylated tumors. Mean OS in the non-methylated group was 16.9 months (Std 1.89, 95% CI 13.2–20.6) compared to 36.1 (Std 6.57, 95% CI 23.2–49.0) in the methylated group (Table 6). PFS was 6.1 months (Std 0.66, 95% CI 4.8–7.4) and 14.2 months (Std 3.87, 95% CI 6.6–21.7) respectively (Table 7).
Table 6. Overall survival in months between O6-methylguanine-DNA methyltransferase (MGMT) status and between larger/lower than median region of interest (ROI) volume in all ROIs segmented through deep learning methods.
Table 7. Progression free survival in months between O6-methylguanine-DNA methyltransferase (MGMT) status and between larger/lower than median region of interest (ROI) volume in all ROIs segmented through deep learning methods.
Analyzing volumes from DeepBraTumIA we found significantly higher OS in lower volumes of contrast enhancing tumor (p < 0.01) (lower volume mean 37.7 months, Std 6.90, 95% CI 21.2–51.2 vs. higher volume mean 16.0 months, Std 1.59, 95% CI 12.9–19.1) and in lower total-tumor/edema ratio (p = 0.02) (lower ratio mean 31.8 months, Std 5.39, 95% CI 21.2–42.3 vs. higher ratio mean 19.9 months, Std 3.66, 95% CI 12.8–27.1) (Table 6). PFS was significantly higher in lower volumes of contrast enhancing tumor (p = 0.01) (lower volume mean 14.3 months, Std 4.13, 95% CI 6.2–22.4 vs. higher volume mean 6.4 months, Std 0.93, 95% CI 4.6–8.2) (Table 7). The remaining volumes showed no significant differences.
Analyzing volumes from Raidionics we found significantly higher OS in lower contrast enhancing volume (p < 0.01) (lower volume mean 36.2 months, Std 6.74, 95% CI 23.0–49.4 vs. higher volume mean 16.8 months, Std 1.57, 95% CI 13.7–19.9), lower necrosis volume (p = 0.04) (lower volume mean 31.5 months, Std 5.59, 95% CI 20.6–42.5 vs. higher volume mean 20.2 months, Std 3.46, 95% CI 13.5–27.0), lower total tumor volume (p = 0.01) (lower volume mean 35.1 months, Std 6.60, 95% CI 22.2–48.0 vs. higher volume mean 17.3 months, Std 1.69, 95% CI 14.0–20.6) and lower total-tumor/edema ratio (p < 0.01) (lower ratio mean 35.3 months, Std 6.20, 95% CI 23.1–47.4 vs. higher ratio mean 17.7 months, Std 3.09, 95% CI 11.7–23.8) (Table 6). Significantly higher PFS was observed in lower total tumor volume (p = 0.04) (lower volume mean 13.6 months, Std 3.87, 95% CI 6.0–21.2 vs. higher volume mean 6.7 months, Std 1.04, 95% CI 4.7–8.7) and lower total-tumor/edema ratio (p = 0.03) (lower ratio mean 13.8 months, Std 3.92, 95% CI 6.1–21.5 vs. higher ratio mean 6.6 months, Std 1.07, 95% CI 4.5–8.7) (Table 7). The remaining volumes showed no significant differences.
3.4.2 GBM-subgroup
OS and PFS were shown to be significantly higher (p < 0.01 and p = 0.04 respectively) in MGMT-promoter methylated compared to non-methylated tumors. Mean OS in the non-methylated group was 21.9 months (Std 1.75, 95% CI 12.5–19.4) compared to 30.0 (Std 6.29, 95% CI 17.7–42.3) in the methylated group (Table 8). PFS was 5.7 months (Std 0.63, 95% CI 4.5–6.9) and 10.2 months (Std 2.38, 95% CI 5.6–14.9) respectively (Table 9).
Table 8. Overall survival in months in glioblastoma (GBM) subgroup between O6-methylguanine-DNA methyltransferase (MGMT) status and between larger/lower than median region of interest (ROI) volume in all ROIs segmented through deep learning methods.
Table 9. Progression free survival in months in glioblastoma (GBM) subgroup between O6-methylguanine-DNA methyltransferase (MGMT) status and between larger/lower than median region of interest (ROI) volume in all ROIs segmented through deep learning methods.
Analyzing volumes from DeepBraTumIA we found significantly higher OS in lower volumes of contrast enhancing tumor (p < 0.01) (lower volume mean 29.4 months, Std 5.88, 95% CI 17.9–40.9 vs. higher volume mean 15.0 months, Std 1.52, 95% CI 12.0–18.0), in lower tumor necrosis (p = 0.03) (lower volume mean 28.2 months, Std 5.82, 95% CI 16.8–39.6 vs. higher volume mean 15.8 months, Std 1.39, 95% CI 13.1–18.6), in lower total tumor volume (p = 0.02) (lower volume mean 30.0 months, Std 5.96, 95% CI 17.3–40.7 vs. higher volume mean 15.5 months, Std 1.51, 95% CI 12.5–18.5) and in lower total-tumor/edema ratio (p < 0.01) (lower ratio mean 30.9 months, Std 5.73, 95% CI 19.7–42.1 vs. higher ratio mean 13.3 months, Std 1.40, 95% CI 10.6–16.1) (Table 8). PFS was significantly higher in lower volumes of contrast enhancing tumor (p = 0.01) (lower volume mean 9.7 months, Std 2.24, 95% CI 5.3–14.1 vs. higher volume mean 5.8 months, Std 0.61, 95% CI 4.6–7.0), in lower tumor necrosis (p = 0.03) (lower volume mean 9.8 months, Std 2.16, 95% CI 5.6–14.0 vs. higher volume mean 5.5 months, Std 0.54, 95% CI 4.5–6.6), in lower total tumor volume (p = 0.02) (lower volume mean 10.2 months, Std 2.19, 95% CI 5.9–14.5 vs. higher volume mean 5.3 months, Std 0.56, 95% CI 4.2–6.4) and in lower total-tumor/edema ratio (p < 0.01) (lower ratio mean 10.0 months, Std 2.35, 95% CI 5.4–14.6 vs. higher ratio mean 5.7 months, Std 0.69, 95&CI 4.4–7.0) (Table 9). The remaining volumes showed no significant differences.
Analyzing volumes from Raidionics we found significantly higher OS in lower contrast enhancing volume (p < 0.05) (lower volume mean 28.0 months, Std 5.86, 95% CI 16.6–39.5 vs. higher volume mean 16.1 months, Std 1.49, 95% CI 13.2–19.0), lower necrosis volume (p < 0.01) (lower volume mean 28.5 months, Std 5.40, 95% CI 17.9–39.1 vs. higher volume mean 14.7 months, Std 1.32, 95% CI 12.1–17.3), lower total tumor volume (p = 0.02) (lower volume mean 28.4 months, Std 5.80, 95% CI 17.0–39.8 vs. higher volume mean 15.7 months, Std 1.51, 95% CI 12.7–18.6) and lower total-tumor/edema ratio (p < 0.01) (lower ratio mean 31.4 months, Std 5.6, 95% CI 20.4–42.4 vs. higher ratio mean 12.6 months, Std 1.26, 95% CI 10.2–15.1) (Table 8). Significantly higher PFS was observed in lower total tumor volume (p = 0.02) (lower volume mean 9.8 months, Std 2.16, 95% CI 5.6–14.0 vs. higher volume mean 5.5 months, Std 0.55, 95% CI 4.5–6.6) and lower total-tumor/edema ratio (p < 0.01) (lower ratio mean 10.2 months, Std 2.31, 95% CI 5.7–14.7 vs. higher ratio mean 5.5 months, Std 0.70, 95% CI 4.1–6.8) (Table 9). The remaining volumes showed no significant differences.
4 Discussion
The ability to accurately predict MGMT-promoter methylation preoperatively could play a vital role in the treatment of HGG. Using volumes gathered by segmenting ROIs, which is more effectively done using DL-based segmentation models, might aid in this distinction. Therefore, this study examined contrast enhancing tumor, central necrosis and peritumoral edema volumes using both automatic and manual segmentation methods, in HGG patients with and without MGMT promotor methylation.
Our data show that larger peritumoral edema volumes, segmented using DeepBraTumIA, positively predicts MGMT promoter methylation. These findings could be in line with Li et al. (2011) who found that low expression of MGMT was associated with presence of more severe edema. However, they did not explore any volumetric correlations beyond visual inspection. Additionally other studies contradict our findings, with one study not finding any significant differences between methylated and non-methylated groups (Kanas et al., 2017) and others having identified MGMT promotor methylation in tumors with lower edema volume (Suh et al., 2018; Kanas et al., 2017; Koska and Koska, 2025).
Although we did find significant group differences in the peritumoral edema ROI with DeepBraTumIA, the ROC analysis (AUC = 0.638) indicates limited clinical utility. These findings are consistent with previous research which, despite showing lower volume of peritumoral edema in MGMT-methylated tumors, also struggles with achieving high predictive potential on models created through a single parameter (Kanas et al., 2017; Koska and Koska, 2025). These discoveries indicate that the use of peritumoral edema volume alone as a predictive MRI feature is not feasible. This conclusion is further supported by our multivariate analysis, where the peritumoral edema ROI failed to remain a significant predictor when tested against other clinical covariates.
It is important to note that peritumoral edema represents a downstream effect of a tumor, and may not be directly linked to expression of MGMT. It can be speculated that factors such as vascular endothelial growth factor (VEGF) and matrix metalloproteinases (MMPs), expressed by the tumor, may play a more prominent role in mediating peritumoral edema trough their effect on the blood brain barrier (BBB) (Ohmura et al., 2023). Previous studies indicate that expression of VEGF within the brain may increase BBB leakage and the formation of abnormal tumor vessels, causing increased localized edema burden (Ohmura et al., 2023; Zhang et al., 2000). Furthermore, extensive expression of MMPs causes destruction of the basal membrane of the BBB, resulting in edema (Ohmura et al., 2023; Zhang et al., 2019). Consequently, the limited predictive ability of peritumoral edema in our study may reflect that such secondary pathophysiological processes, rather than MGMT promoter methylation itself, exert the dominant influence on edema formation.
Additionally, tumor edema is influenced by symptomatic intake of corticosteroids. Non-MGMT-methylated tumors may have more aggressive clinical features and thus a higher proportion of steroid use pre-operatively. Since most patients in our cohort had either recently initiated corticosteroid therapy or were on a low to no dose at the time of imaging, we expect the impact of such treatment to be minimal on the peritumoral edema. Additionally, we observed no significant difference in mean dosage of betamethasone between patients stratified by MGMT status (Table 1). Our multivariate analysis also showed no association between corticosteroid dosage and MGMT in the tested prediction models. This suggests that corticosteroids did not act as a confounder in our predictive model including the tumor/edema ratio.
Some studies have found a lower survival rate in patients with larger volumes of peritumoral edema (Carrillo et al., 2012; Pope et al., 2005; Hammoud et al., 1996). The lower survival rate is theorized to be related to the increased mass effect of the peritumoral edema on the surrounding tissue (Carlson et al., 2007). As previously mentioned, peritumoral edema could be a consequence of increased expression of VEGF from the HGG. This could potentially indicate a more hypoxic and aggressive tumor (Carlson et al., 2007). Lower survival rates are also reported in patients with VEGF-expressing HGG and high peritumoral edema volumes, whereas VEGF-positive tumors with low or no peritumoral edema had no significant impact on survival (Carlson et al., 2007), supporting the notion that larger peritumoral edema is found within more aggressive tumors. Additionally, pro-tumor associated gene expressions, such as p53 suppression, have been found within peritumoral edema lowering the delay between surgery and tumor recurrence (Luo et al., 2021). However, our findings are contradictory to these findings, suggesting no difference in OS or PFS based on edema burden. Furthermore, we can conclude that MGMT-promoter methylated tumors offer significantly higher OS compared to non-methylated. Considering our findings that methylated tumors express larger volumes of peritumoral edema we theorize that edema burden might not be as impactful on OS compared to the contrast-enhancing tumor volume and total tumor/edema ratio (Carrillo et al., 2012).
As mentioned previously our results are contradictory to some previous studies, primarily those which identified MGMT methylation in tumors with lower amount of peritumoral edema volume (Suh et al., 2018; Kanas et al., 2017; Koska and Koska, 2025). It is possible that patients in our cohort expressing MGMT promoter methylation may show up during later stages of their disease, compared to non-methylated, if the tumor is less aggressive. Furthermore, previous studies have demonstrated that tumors with higher degree of MGMT methylation, measured with pyrosequencing, exhibit a larger volumetric peritumoral edema reduction as response to eventual treatment, and prolonged survival (Hosoya et al., 2022). This agrees with our results, which showed a longer OS in patients with MGMT methylated tumors, despite having larger preoperative peritumoral edema. It is important to note that the MGMT-methylated group included a significantly higher proportion of grade 4 astrocytoma patients compared to the non-methylated group. However, we see similar results in OS and PFS in the GBM subgroup analysis when looking at MGMT methylation. Furthermore, it’s important to note that there is currently no robust imaging-based method to differentiate between GBM and astrocytoma grade 4 pre-operatively when assessing HGG. As a result, it would be of limited clinical value to further assess the effect of the skewed group composition. Additionally, the inclusion of both subtypes in the overall cohort reflects the clinical reality faced at diagnosis while the OS and PFS findings indicate a correlation between the mixed group and GBM subgroup in relation to MGMT-methylation (significantly higher in both) suggesting that the observed increased survival may be driven by MGMT status rather than group composition.
The largest changes to OS and PFS between the GBM subgroup and the mixed cohort were observed in tumor necrosis volume, total tumor volume, and total-tumor/edema ratio segmented by DeepBraTumIA. In the GBM-only analysis, all of these showed significantly longer OS or PFS in patients with lower volumes/ratio, whereas no such associations were seen when astrocytomas were included. However, volumes segmented by Raidionics did not show such a change in results in the subgroup. This pattern suggests that the differences observed are a consequence of DeepBraTumIA ability to segment astrocytomas in the mixed cohort, rather than reflecting true OS and PFS differences between the mixed cohort and the subgroup.
We also found the total-tumor/edema ratio to be higher in non-MGMT promoter methylated tumors. A higher ratio is derived from a larger total tumor or smaller peritumoral edema volume, which aligns with our findings showing larger peritumoral edema volume in MGMT-promoter methylated tumors. Using segmentation from both DL-models significant differences and improved predictive capabilities were found in total-tumor/edema ratio which could indicate that the combination of multiple ROIs could have a larger applicability and be less sensitive to specific segmentation methods, compared to individual ROIs. This is in line with previous studies showing higher AUC values when combining multiple metrics (Li et al., 2024). Furthermore, our multivariate analysis indicates that lower ratio was predictive of MGMTpm methylation and that the predictive findings are truly related to regional variations seen on MRI instead of a consequence of other covariates. Ratios segmented by both DL-models were the only significant factor in the logistic regression models when tested against our clinical covariates. Such results were not present when modeling using peritumoral edema segmented by DeepBraTumIA, strengthening the hypothesis that a singular ROI is not sufficient in capturing the complexity of HGG and predicting MGMT-methylation. Despite all these findings, the ratio’s predictive utility is limited by its low AUC in ROC analysis while improving prediction over baseline in multivariate regression.
We observed improved predictive ability in both total-tumor/edema ratios when analyzing the GBM subgroup compared to the mixed cohort, albite only by a slight margin. Despite the prediction not being clinically applicable, the implication of these findings is important to note. As previously mentioned, no robust pre-operative imaging method currently exists to distinguish HGG gliomas from each other in pre-operative imaging. As such, knowing that we have similar ability to predict MGMT methylation in both a mixed cohort and GBM subgroup indicates that further studies using cohorts that mirror this clinical reality are warranted.
Whereas we report a higher total-tumor/edema ratio to be predictive of non MGMT-methylation, another study reports the opposite (Kanas et al., 2017) while other studies found no correlation between MRI features and MGMT promoter methylation (Kim et al., 2022; Mikkelsen et al., 2020). The discrepancy could indicate other factors in play which could impact peritumoral edema.
We speculate that one such factor could be the immune response surrounding the tumor. Previous studies suggest that MGMT status may modulate and impact the local immune response to the tumor (Kushihara et al., 2024). It is suggested that higher expressions of MGMT may increase CD8 T-cells abundance (Kushihara et al., 2024). This response may in turn change the microenvironment surrounding the tumor in a way which promotes a heightened inflammatory response and subsequently increased peritumoral edema (Engelhorn et al., 2009). This is further complicated by the heterogeneous expression of MGMT within the tumor, as each subpopulation within the tumor may create distinct unique microenvironment due to different levels of immune cell recruitment (Singh et al., 2025). Collectively, these factors highlight the complexity of MGMT prediction in HGG based on a single factor approach. They also raise the possibility that varying degrees of homogeneity and expression of MGMT may impact prediction even within non methylated MGMT tumors, which could be explored in future studies. Some studies have shown that factors such as apparent diffusion coefficient (ADC) (Rundle-Thiele et al., 2015; Kanazawa et al., 2019) and textural differences of the images (Drabycz et al., 2010) could instead be used as a basis for an MGMT-methylation prediction model. Furthermore, amide proton transfer weighted imaging has recently also shown promising results in predicting MGMT in GBM and grade 4 astrocytoma (Durmo et al., 2025). Given the complex nature of both the brain and HGG, we believe that a singular-factor approach may be insufficient for predicting MGMT methylation. Our multivariate analysis confirms this theory. While our models did become statistically significant while accounting for co-variates, the overall predictive ability remained limited with a singular ROI. A more holistic study design, incorporating multiple factors and their relation to each other, could be a better approach to find factors which predict MGMT-promotor methylation. Thus, future studies should take multiple factors into consideration, looking for a linked correlation between the overall tumor image and MGMT status, instead of focusing on a singular ROI.
Due to the limited predictive performance of our models, we chose not to perform PFS or OS analyses stratified by predicted MGMT status. Given the modest AUC values observed for our ROIs, such analyses would not meaningfully reflect the prognostic impact of MGMT methylation but rather the effects of the volumetric features themselves on survival. The effect of higher and lower ROI volumes on PFS and OS is an association we already demonstrated. As a result, we concluded that additional survival analyses based on predicted MGMT status are unlikely to offer added clinical insights. We observed, when comparing DL-model segmentation to manual segmentation, that Raidionics differ in fewer ROIs than DeepBraTumIA. This discrepancy may be due to differences in model architecture or potential biases in Raidionics’ training data that better align with our dataset. However, these findings do not necessarily indicate that Raidionics is a superior segmentation model, only that it corresponds to a higher degree with our specific manual segmentation. Given the well-documented inter-observer variability among experienced radiologists (Erasmus et al., 2003; Zhao et al., 2013; Barboriak et al., 2019), manual segmentation itself shows considerable variation. Rather than comparing DL models solely to human segmentation, their performance should be assessed based on task-specific accuracy. Thus, their agreeability in the ability to predict MGMT-methylation suggests that they are equally suited for that task.
Analysis of the Dice coefficients reveals that tumor necrosis is the most inconsistent ROI across all segmentation methods. No comparison had greater than 73% overlap, in contrast with the remaining ROIs which had a range of 82–91%. These results indicate a large discrepancy in tumor necrosis definition among the segmentation methods. For the manual segmentation the tumor necrosis was defined as a hypointense area, surrounded by contrast enhancing tumor on the T1-GD images. Visual observation reveals that in the DL models’ segmentations, tumor necrosis may be defined as hypointense regions adjacent to contrast enhancing tumor tissue, resulting in the labeling of necrosis outside of the contrast enhancing region. This can be seen when comparing the performed DL segmentations in two different patients (Figure 3). These findings highlight the need for a more robust and standardized segmentation criteria for tumor necrosis. The DL-models rely on their segmentation of training data that incorporates true labels from experienced radiologists. The observed discrepancy between segmentation methods likely reflects an underlying bias and variability within the training data. This is further supported by the results of the ICC, showing excellent inter-observer agreement for contrast-enhancing tumor and edema ROIs, while at best having moderate agreement within the necrosis ROI. Previous studies have also identified similar limitations in using DL-models on segmentation (Bonada et al., 2024). Furthermore, as tumor necrosis volume is the smallest ROI across all segmentation methods, even minor differences in labeling definitions could result in large impact on the dice-coefficient. DeepBraTumIA, judging by volume, seems to be more likely to over segment tumor necrosis, compared to the other segmentation methods.
Such findings are likely to impact on our ability to reliably predict MGMT promoter methylation within the tumor necrosis ROI. Even if significant group differences were observed, clinical applicability would remain limited, as high interobserver variability would translate to high variability in prediction. However, due to the definition used for other ROIs in our study, especially total tumor volume and in turn total-tumor/edema ratio, we expect that the downstream effect on them to be minimal.
Our study has some limitations. Firstly, our study employs a single center study design with a cohort of 70 patients, limiting statistical power and generalizability. For future studies a wider dataset from multiple centers is recommended to enhance the generalizability of the results and improve statistical power. Secondly, the manual segmentations were in large performed by a single assessor. It could have been performed by multiple readers independently, to account for inter-observer variance and improve the comparison between the manual and DL segmentations. However, the manual segmentation performed by EZ was reviewed by both an experienced researcher (TS) and a neuroradiologist with more than 30 years of experience (PS). Additionally, the manual segmentation in our study was only performed in a subgroup of patients as an initial validation of the DL models’ performance, as the main aim was to evaluate imaging features in the predictive ability of MGMT status. Lastly, due to the complexity of segmentation and deep learning, a large dataset, using labeling and segmenting from different radiologists, would be recommended to accurately compare DL-models to each other. Recent studies have highlighted this methodology by using large multi-institutional datasets to develop and validate a post-operative DL-model against other tools, such as Raidionics (Cepeda et al., 2024).
5 Conclusion
This study demonstrated significant differences between MGMT promoter methylated and non-methylated HGG in the peritumoral edema volume and total-tumor/edema ratio. Total-tumor/edema ratio seems to be the single most predictive factor for MGMT-status using multivariate analysis. However, the diagnostic performance of both peritumoral edema and ratio were limited. OS and PFS were longer in the MGMT-promoter methylated group. Volumes obtained from both DL-models were both equally suited for predicting MGMT status, highlighting the potential use of such models clinically as a tool to help radiologists.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
The studies involving humans were approved by the Regional Ethics Board Lund, Sweden and the Swedish Ethical Review Authority prior to starting any work. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.
Author contributions
EZ: Data curation, Formal analysis, Investigation, Methodology, Visualization, Writing – original draft, Writing – review & editing. TS: Data curation, Methodology, Resources, Software, Supervision, Validation, Writing – review & editing. MK: Investigation, Validation, Writing – review & editing. XS: Resources, Writing – review & editing. JL: Data curation, Writing – review & editing. SK: Investigation, Resources, Writing – review & editing. MB: Resources, Writing – review & editing. AR: Resources, Writing – review & editing. JB: Resources, Writing – review & editing. LK: Writing – review & editing, Methodology. PS: Conceptualization, Funding acquisition, Methodology, Project administration, Resources, Supervision, Validation, Writing – review & editing.
Funding
The author(s) declared that financial support was received for this work and/or its publication. The study was supported by the Swedish Cancer Society CAN 2024 24 3568 (PS) and CAN 23 2937 Pj 01 H (JB), Swedish Research Council VR 2023-02412 (PS), Regional Research Funds (ALF) F 2022/0108 (PS), the Fru Berta Kamprad Foundation (JB), and the Skåne University Hospital donation funds (JB).
Acknowledgments
The authors acknowledge Associate Professor Markus Nilsson, Department of Clinical Sciences/Division of Radiology, Lund University, Skåne University Hospital, Lund, Sweden for valuable insight and input. Research nurse Linda Wennberg, MR radiographer, PhD student, Department of Clinical Sciences/Division of Radiology, Lund University, Skåne University Hospital, Lund, Sweden for her support in patient recruitment and MR scanning.
Conflict of interest
The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declared that Generative AI was not used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
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.
References
Barboriak, D. P., Zhang, Z., Desai, P., Snyder, B. S., Safriel, Y., McKinstry, R. C., et al. (2019). Interreader variability of dynamic contrast-enhanced MRI of recurrent glioblastoma: the multicenter ACRIN 6677/RTOG 0625 study. Radiology 290, 467–476. doi: 10.1148/radiol.2019181296,
Bonada, M., Rossi, L. F., Carone, G., Panico, F., Cofano, F., Fiaschi, P., et al. (2024). Deep learning for MRI segmentation and molecular subtyping in glioblastoma: critical aspects from an emerging field. Biomedicine 12:1878. doi: 10.3390/biomedicines12081878,
Bouget, D., Alsinan, D., Gaitan, V., Helland, R. H., Pedersen, A., Solheim, O., et al. (2023). Raidionics: an open software for pre- and postoperative central nervous system tumor segmentation and standardized reporting. Sci. Rep. 13:15570. doi: 10.1038/s41598-023-42048-7,
Carlson, M. R., Pope, W. B., Horvath, S., Braunstein, J. G., Nghiemphu, P., Tso, C. L., et al. (2007). Relationship between survival and edema in malignant gliomas: role of vascular endothelial growth factor and neuronal pentraxin 2. Clin. Cancer Res. 13, 2592–2598. doi: 10.1158/1078-0432.CCR-06-2772,
Carrillo, J. A., Lai, A., Nghiemphu, P. L., Kim, H. J., Phillips, H. S., Kharbanda, S., et al. (2012). Relationship between tumor enhancement, edema, IDH1 mutational status, MGMT promoter methylation, and survival in glioblastoma. AJNR Am. J. Neuroradiol. 33, 1349–1355. doi: 10.3174/ajnr.A2950,
Cepeda, S., Romero, R., Luque, L., Garcia-Perez, D., Blasco, G., Luppino, L. T., et al. (2024). Deep learning-based postoperative glioblastoma segmentation and extent of resection evaluation: development, external validation, and model comparison. Neurooncol Adv. 6:vdae199. doi: 10.1093/noajnl/vdae199,
DeepBraTumIA. (2025). DeepBraTumIA NITRIC 2021. Available online at: https://www.nitrc.org/projects/deepbratumia/. (Accessed June 3, 2025)
Drabycz, S., Roldan, G., de Robles, P., Adler, D., JB, M. I., Magliocco, A. M., et al. (2010). An analysis of image texture, tumor location, and MGMT promoter methylation in glioblastoma using magnetic resonance imaging. NeuroImage 49, 1398–1405. doi: 10.1016/j.neuroimage.2009.09.049,
Durmo, F., Latt, J., Rydelius, A., Englund, E., Salomonsson, T., Liebig, P., et al. (2025). Role of amide proton transfer weighted MRI in predicting MGMTp methylation status, p53-status, Ki-67 index, IDH-status, and ATRX expression in WHO grade 4 high grade glioma. Tomography. 11:64. doi: 10.3390/tomography11060064,
Ellingson, B. M., Cloughesy, T. F., Pope, W. B., Zaw, T. M., Phillips, H., Lalezari, S., et al. (2012). Anatomic localization of O6-methylguanine DNA methyltransferase (MGMT) promoter methylated and unmethylated tumors: a radiographic study in 358 de novo human glioblastomas. NeuroImage 59, 908–916. doi: 10.1016/j.neuroimage.2011.09.076,
Engelhorn, T., Savaskan, N. E., Schwarz, M. A., Kreutzer, J., Meyer, E. P., Hahnen, E., et al. (2009). Cellular characterization of the peritumoral edema zone in malignant brain tumors. Cancer Sci. 100, 1856–1862. doi: 10.1111/j.1349-7006.2009.01259.x,
Erasmus, J. J., Gladish, G. W., Broemeling, L., Sabloff, B. S., Truong, M. T., Herbst, R. S., et al. (2003). Interobserver and intraobserver variability in measurement of non-small-cell carcinoma lung lesions: implications for assessment of tumor response. J. Clin. Oncol. 21, 2574–2582. doi: 10.1200/JCO.2003.01.144,
Esteller, M., and Herman, J. G. (2004). Generating mutations but providing chemosensitivity: the role of O6-methylguanine DNA methyltransferase in human cancer. Oncogene 23, 1–8. doi: 10.1038/sj.onc.1207316,
Fedorov, A., Beichel, R., Kalpathy-Cramer, J., Finet, J., Fillion-Robin, J. C., Pujol, S., et al. (2012). 3D slicer as an image computing platform for the quantitative imaging network. Magn. Reson. Imaging 30, 1323–1341. doi: 10.1016/j.mri.2012.05.001,
Fernandes, C., Costa, A., Osorio, L., Lago, R. C., Linhares, P., Carvalho, B., et al. (2017). “Current standards of care in glioblastoma therapy” in Glioblastoma (Brisbane, QLD: Codon Publications).
GitHub. (2023). Raidionics. Tumor segmentation. Available online at: https://github.com/raidionics/Raidionics/wiki/Tumor-segmentation
Hammoud, M. A., Sawaya, R., Shi, W., Thall, P. F., and Leeds, N. E. (1996). Prognostic significance of preoperative MRI scans in glioblastoma multiforme. J. Neuro-Oncol. 27, 65–73. doi: 10.1007/BF00146086,
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,
Hegi, M. E., Oppong, F. B., Perry, J. R., Wick, W., Henriksson, R., Laperriere, N. J., et al. (2024). No benefit from TMZ treatment in glioblastoma with truly unmethylated MGMT promoter: reanalysis of the CE.6 and the pooled Nordic/NOA-08 trials in elderly glioblastoma patients. Neuro-Oncol. 26, 1867–1875. doi: 10.1093/neuonc/noae108,
Hosoya, T., Takahashi, M., Honda-Kitahara, M., Miyakita, Y., Ohno, M., Yanagisawa, S., et al. (2022). MGMT gene promoter methylation by pyrosequencing method correlates volumetric response and neurological status in IDH wild-type glioblastomas. J. Neuro-Oncol. 157, 561–571. doi: 10.1007/s11060-022-03999-5,
Isensee, F., Schell, M., Pflueger, I., Brugnara, G., Bonekamp, D., Neuberger, U., et al. (2019). Automated brain extraction of multisequence MRI using artificial neural networks. Hum. Brain Mapp. 40, 4952–4964. doi: 10.1002/hbm.24750,
Kanas, V. G., Zacharaki, E. I., Thomas, G. A., Zinn, P. O., Megalooikonomou, V., and Colen, R. R. (2017). Learning MRI-based classification models for MGMT methylation status prediction in glioblastoma. Comput. Methods Prog. Biomed. 140, 249–257. doi: 10.1016/j.cmpb.2016.12.018,
Kanazawa, T., Minami, Y., Jinzaki, M., Toda, M., Yoshida, K., and Sasaki, H. (2019). Predictive markers for MGMT promoter methylation in glioblastomas. Neurosurg. Rev. 42, 867–876. doi: 10.1007/s10143-018-01061-5,
Kim, B. H., Lee, H., Choi, K. S., Nam, J. G., Park, C. K., Park, S. H., et al. (2022). Validation of MRI-based models to predict MGMT promoter methylation in gliomas: BraTS 2021 radiogenomics challenge. Cancers 14:4827. doi: 10.3390/cancers14194827,
Koska, I. O., and Koska, C. (2025). Deep learning classification of MGMT status of glioblastomas using multiparametric MRI with a novel domain knowledge augmented mask fusion approach. Sci. Rep. 15:3273. doi: 10.1038/s41598-025-87803-0,
Kushihara, Y., Tanaka, S., Kobayashi, Y., Nagaoka, K., Kikuchi, M., Nejo, T., et al. (2024). Glioblastoma with high O6-methyl-guanine DNA methyltransferase expression are more immunologically active than tumors with low MGMT expression. Front. Immunol. 15:1328375. doi: 10.3389/fimmu.2024.1328375,
Li, W. B., Tang, K., Zhang, W., Yan, W., You, G., Li, S. W., et al. (2011). Relationship between magnetic resonance imaging and molecular pathology in patients with glioblastoma multiforme. Chin. Med. J. 124, 2589–2592. doi: 10.3760/cma.j.issn.0366-6999.2011.17.005,
Li, L., Xiao, F., Wang, S., Kuang, S., Li, Z., Zhong, Y., et al. (2024). Preoperative prediction of MGMT promoter methylation in glioblastoma based on multiregional and multi-sequence MRI radiomics analysis. Sci. Rep. 14:16031. doi: 10.1038/s41598-024-66653-2,
Louis, D. N., Perry, A., Wesseling, P., Brat, D. J., Cree, I. A., Figarella-Branger, D., et al. (2021). The 2021 WHO classification of tumors of the central nervous system: a summary. Neuro-Oncol. 23, 1231–1251. doi: 10.1093/neuonc/noab106,
Luo, C., Song, K., Wu, S., Hameed, N. U. F., Kudulaiti, N., Xu, H., et al. (2021). The prognosis of glioblastoma: a large, multifactorial study. Br. J. Neurosurg. 35, 555–561. doi: 10.1080/02688697.2021.1907306,
Mansouri, A., Hachem, L. D., Mansouri, S., Nassiri, F., Laperriere, N. J., Xia, D., et al. (2019). MGMT promoter methylation status testing to guide therapy for glioblastoma: refining the approach based on emerging evidence and current challenges. Neuro-Oncol. 21, 167–178. doi: 10.1093/neuonc/noy132,
Mikeska, T., Bock, C., El-Maarri, O., Hubner, A., Ehrentraut, D., Schramm, J., et al. (2007). Optimization of quantitative MGMT promoter methylation analysis using pyrosequencing and combined bisulfite restriction analysis. J. Mol. Diagn. 9, 368–381. doi: 10.2353/jmoldx.2007.060167,
Mikkelsen, V. E., Dai, H. Y., Stensjøen, A. L., Berntsen, E. M., Salvesen, Ø., Solheim, O., et al. (2020). MGMT promoter methylation status is not related to histological or radiological features in IDH wild-type glioblastomas. J. Neuropathol. Exp. Neurol. 79, 855–862. doi: 10.1093/jnen/nlaa060,
Ohmura, K., Tomita, H., and Hara, A. (2023). Peritumoral edema in gliomas: a review of mechanisms and management. Biomedicine 11:2731. doi: 10.3390/biomedicines11102731,
Ostrom, Q. T., Gittleman, H., Truitt, G., Boscia, A., Kruchko, C., and Barnholtz-Sloan, J. S. (2018). CBTRUS statistical report: primary brain and other central nervous system tumors diagnosed in the United States in 2011–2015. Neuro-Oncol. 20, iv1–iv86. doi: 10.1093/neuonc/noy131,
Pinter, C., Lasso, A., Wang, A., Jaffray, D., and Fichtinger, G. (2012). SlicerRT: radiation therapy research toolkit for 3D slicer. Med. Phys. 39, 6332–6338. doi: 10.1118/1.4754659,
Pope, W. B., Sayre, J., Perlina, A., Villablanca, J. P., Mischel, P. S., and Cloughesy, T. F. (2005). MR imaging correlates of survival in patients with high-grade gliomas. AJNR Am. J. Neuroradiol. 26, 2466–2474,
Rundle-Thiele, D., Day, B., Stringer, B., Fay, M., Martin, J., Jeffree, R. L., et al. (2015). Using the apparent diffusion coefficient to identifying MGMT promoter methylation status early in glioblastoma: importance of analytical method. J. Med. Radiat. Sci. 62, 92–98. doi: 10.1002/jmrs.103,
Singh, N., Miner, A., Hennis, L., and Mittal, S. (2021). Mechanisms of temozolomide resistance in glioblastoma—a comprehensive review. Cancer Drug Resist. 4, 17–43. doi: 10.20517/cdr.2020.79,
Singh, S., Nema, R., Banerjee, M., and Kushwah, A. S. (2025). Genetic and epigenetic landscape of O(6)-methylguanine-DNA methyltransferase (MGMT): implications for DNA repair and cancer therapeutics. Explor. Target. Antitumor. Ther. 6:1002335. doi: 10.37349/etat.2025.1002335,
Stupp, R., Mason, W. P., van den Bent, M. J., Weller, M., Fisher, B., Taphoorn, M. J., et al. (2005). Radiotherapy plus concomitant and adjuvant temozolomide for glioblastoma. N. Engl. J. Med. 352, 987–996. doi: 10.1056/NEJMoa043330
Suh, C. H., Kim, H. S., Jung, S. C., Choi, C. G., and Kim, S. J. (2018). Clinically relevant imaging features for MGMT promoter methylation in multiple glioblastoma studies: a systematic review and meta-analysis. AJNR Am. J. Neuroradiol. 39, 1439–1445. doi: 10.3174/ajnr.A5711,
Torp, S. H., Solheim, O., and Skjulsvik, A. J. (2022). The WHO 2021 classification of central nervous system tumours: a practical update on what neurosurgeons need to know—a minireview. Acta Neurochir. 164, 2453–2464. doi: 10.1007/s00701-022-05301-y,
Weller, M., Cloughesy, T., Perry, J. R., and Wick, W. (2013). Standards of care for treatment of recurrent glioblastoma—Are we there yet? Neuro-Oncol. 15, 4–27. doi: 10.1093/neuonc/nos273,
Wick, W., Platten, M., Meisner, C., Felsberg, J., Tabatabai, G., Simon, M., et al. (2012). Temozolomide chemotherapy alone versus radiotherapy alone for malignant astrocytoma in the elderly: the NOA-08 randomised, phase 3 trial. Lancet Oncol. 13, 707–715. doi: 10.1016/S1470-2045(12)70164-X,
Witthayanuwat, S., Pesee, M., Supaadirek, C., Supakalin, N., Thamronganantasakul, K., and Krusun, S. (2018). Survival analysis of glioblastoma multiforme. Asian Pac. J. Cancer Prev. 19, 2613–2617. doi: 10.22034/APJCP.2018.19.9.2613,
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,
Zhang, H., Ma, Y., Wang, H., Xu, L., and Yu, Y. (2019). MMP-2 expression and correlation with pathology and MRI of glioma. Oncol. Lett. 17, 1826–1832. doi: 10.3892/ol.2018.9806,
Zhang, Z. G., Zhang, L., Jiang, Q., Zhang, R., Davies, K., Powers, C., et al. (2000). VEGF enhances angiogenesis and promotes blood-brain barrier leakage in the ischemic brain. J. Clin. Invest. 106, 829–838. doi: 10.1172/JCI9369,
Zhao, B., Tan, Y., Bell, D. J., Marley, S. E., Guo, P., Mann, H., et al. (2013). Exploring intra- and inter-reader variability in uni-dimensional, bi-dimensional, and volumetric measurements of solid tumors on CT scans reconstructed at different slice intervals. Eur. J. Radiol. 82, 959–968. doi: 10.1016/j.ejrad.2013.02.018,
Glossary
ADC - Apparent diffusion coefficient
AUC - Area under curve
BBB - Blood brain barrier
BraTS - Brain tumor segmentation
DL - Deep learning
GBM - Glioblastoma
GTR - Gross total resection
HD-BET - HD Brain Extraction Tool
HGG - High grade glioma
ICC - Intra-class correlation coefficient
IDH - Isocitrate dehydrogenase
MGMT - O6-methylguanine-DNA methyltransferase
MMP - Matrix metalloproteinase
MPRAGE - Magnetization prepared rapid gradient echo
MRI - Magnetic resonance imaging
OS - Overall survival
PCR - Polymerase chain reaction
PFS - Progression free survival
ROC - Receiver operation characteristics
ROI - Region of interest
RT - Radiotherapy
T1-GD - T1 gadolinium contrast
T2-FLAIR - T2-weighted fluid-attenuated inversion recovery
TE - Echo time
TI - Inversion time
TMZ - Temozolomide
TR - Repetition time
TSE - Turbo spin echo
VEGF - Vascular endothelial growth factor
Keywords: MGMT, methylation, MRI, deep learning, high grade glioma (HGG)
Citation: Zahirovic E, Salomonsson T, Knutsson M, Sarda XS, Lätt J, Kinhult S, Belting M, Rydelius A, Bengzon J, Knutsson L and Sundgren PC (2025) Noninvasive MGMT-promotor methylation prediction in high grade gliomas using conventional MRI and deep learning-based segmentations. Front. Neurosci. 19:1689003. doi: 10.3389/fnins.2025.1689003
Edited by:
Dong Zeng, Southern Medical University, ChinaReviewed by:
Miao Chang, The First Affiliated Hospital of China Medical University, ChinaPietro Fiaschi, University of Genoa, Italy
Ping Hu, Second Affiliated Hospital of Nanchang University, China
Copyright © 2025 Zahirovic, Salomonsson, Knutsson, Sarda, Lätt, Kinhult, Belting, Rydelius, Bengzon, Knutsson and Sundgren. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Pia C. Sundgren, cGlhLnN1bmRncmVuQG1lZC5sdS5zZQ==
Xavier Saenz Sarda2