Advertisement
Research ArticleAdult Brain
Open Access

MRI-Based Deep-Learning Method for Determining Glioma MGMT Promoter Methylation Status

C.G.B. Yogananda, B.R. Shah, S.S. Nalawade, G.K. Murugesan, F.F. Yu, M.C. Pinho, B.C. Wagner, B. Mickey, T.R. Patel, B. Fei, A.J. Madhuranthakam and J.A. Maldjian
American Journal of Neuroradiology May 2021, 42 (5) 845-852; DOI: https://doi.org/10.3174/ajnr.A7029

This article has a correction. Please see:

Abstract

Editorial expression of concern:

In the May 2021 edition, the American Journal of Neuroradilogy published the article “MRI-Based Deep-Learning Method for Determining Glioma MGMT Promoter Methylation Status” by Yogananda CGB, et al.1 On August 22, 2022, the authors self-reported data errors related to the computer code and the training and testing data sets. The authors are now in the process of re-evaluating the accuracies using the correct test set. This notice of concern is to inform readers about these possible issues related to this articles results. After additional tests from the authors on the correct dataset are available, we will determine what additional action is warranted, such as an erratum.

1. Yogananda CGB, Shah BR, Nalawade SS, et al. MRI-based deep-learning method for determining glioma MGMT promoter methylation status. American Journal of Neuroradilogy. 2021;42(5):845-852. doi:10.3174/AJNR.A7029

BACKGROUND AND PURPOSE: O6-Methylguanine-DNA methyltransferase (MGMT) promoter methylation confers an improved prognosis and treatment response in gliomas. We developed a deep learning network for determining MGMT promoter methylation status using T2 weighted Images (T2WI) only.

MATERIALS AND METHODS: Brain MR imaging and corresponding genomic information were obtained for 247 subjects from The Cancer Imaging Archive and The Cancer Genome Atlas. One hundred sixty-three subjects had a methylated MGMT promoter. A T2WI-only network (MGMT-net) was developed to determine MGMT promoter methylation status and simultaneous single-label tumor segmentation. The network was trained using 3D-dense-UNets. Three-fold cross-validation was performed to generalize the performance of the networks. Dice scores were computed to determine tumor-segmentation accuracy.

RESULTS: The MGMT-net demonstrated a mean cross-validation accuracy of 94.73% across the 3 folds (95.12%, 93.98%, and 95.12%, [SD, 0.66%]) in predicting MGMT methylation status with a sensitivity and specificity of 96.31% [SD, 0.04%] and 91.66% [SD, 2.06%], respectively, and a mean area under the curve of 0.93 [SD, 0.01]. The whole tumor-segmentation mean Dice score was 0.82 [SD, 0.008].

CONCLUSIONS: We demonstrate high classification accuracy in predicting MGMT promoter methylation status using only T2WI. Our network surpasses the sensitivity, specificity, and accuracy of histologic and molecular methods. This result represents an important milestone toward using MR imaging to predict prognosis and treatment response.

ABBREVIATIONS:

IDH
isocitrate dehydrogenase
MGMT
O6-methylguanine-DNA methyltransferase
PCR
polymerase chain reaction
T2WI
T2 weighted Images
TCGA
The Cancer Genome Atlas
TCIA
The Cancer Imaging Archive

O6-methylguanine-DNA methyltransferase (MGMT) promoter methylation is a molecular biomarker of gliomas that has prognostic and therapeutic implications. Unlike isocitrate dehydrogenase (IDH) mutations and 1p/19q co-deletions, MGMT promoter methylation is an epigenetic event. Epigenetic events are functionally relevant but do not involve a change in the nucleotide sequence. Therefore, while MGMT promoter methylation is an important prognostic marker, it does not define a distinct subset of gliomas. MGMT is a DNA repair enzyme that protects normal and glioma cells from alkylating chemotherapeutic agents. The methylation of the MGMT promoter is an example of epigenetic silencing, which results in a loss of function of the MGMT enzyme and its protective effect on glioma cells. The survival benefit incurred by MGMT promoter methylation in patients treated with temozolomide (TMZ) was determined in 2005.1 Subsequent work by Stupp et al2 has shown that in patients who received both radiation and temozolomide, MGMT promoter methylation improved median survival compared with patients with unmethylated gliomas (21.7 versus 12.7 months).2 Long-term follow-up from that initial study has further substantiated the survival benefit.2,3 As a result, determining MGMT promoter methylation status is an important step in predicting survival and determining treatment.

Currently, the only reliable way to determine MGMT promoter methylation status requires analysis of glioma tissue obtained either via an invasive brain biopsy or following open surgical resection. Surgical procedures carry the risk of neurologic injury and complications. Therefore, considerable attention has been dedicated to developing noninvasive, image-based diagnostic methods to determine MGMT promoter methylation status. A meta-analysis of MR imaging features demonstrated that glioblastomas with methylated MGMT promoters were associated with less edema, high ADC, and low perfusion. However, the summary sensitivity and specificity of these clinical features was only 79% and 78%, respectively.4 Although multiple radiomic approaches have also been attempted for MGMT prediction, none, to date, have achieved accuracies sufficient for clinical viability.5⇓⇓⇓-9 Sasaki et al10 attempted to establish an MR imaging–based radiomic model for predicting MGMT promoter status of the tumor, but it reached a predictive accuracy of only 67%. Wei et al11 extracted radiomic features from the tumor and peritumoral edema using multisequence, postcontrast MR imaging but only achieved an accuracy of 51%–74% in predicting MGMT promoter methylation status in astrocytomas. Drabycz et al5 performed texture analysis on MR images to predict MGMT promoter methylation status in glioblastoma, but it reached an accuracy of only 71%. Korfiatis et al9 combined texture features with supervised classification schemes as potential imaging biomarkers for predicting the MGMT methylation status of glioblastoma multiforme, but they achieved a sensitivity and specificity of only 0.803 and 0.813, respectively. Ahn et al7 used dynamic contrast-enhanced MR imaging and diffusion tensor imaging to predict MGMT promoter methylation in glioblastoma, but this method achieved a sensitivity and specificity of only 56.3% and 85.2%, respectively.

Recent advances in deep learning methods have also been used for noninvasive, image-based molecular profiling. Our group has previously demonstrated highly accurate, MR imaging–based, voxelwise, deep learning networks for determining IDH classification and 1p/19q co-deletion status using only T2WI.12,13 The benefits of using T2WIs are that they are routinely acquired, they can be obtained quickly, and high quality T2WI can even be obtained in the setting of motion degradation. Because MGMT promoter methylation in gliomas is such an important biomarker, we sought to develop a highly accurate, fully automated deep learning 3D network for MGMT promoter determination of methylation status using only T2WI.

MATERIALS AND METHODS

Data and Preprocessing

Multiparametric MR images of patients with brain gliomas were obtained from The Cancer Imaging Archive (TCIA) data base.14,15 The genomic information was obtained from both The Cancer Genome Atlas (TCGA) and TCIA data bases.14,16,17 Subject datasets were screened for the availability of preoperative MR images, T2WI, and known MGMT promoter status. The final dataset of 247 subjects included 163 methylated cases and 84 unmethylated cases. TCGA subject identification, MGMT status, and tumor grade are listed in Table 1 of the Online Supplemental Data.

View this table:

Cross-validation results

Tumor masks for 179 subjects were available through previous expert segmentation.18⇓-20 Tumor masks for the remaining 68 subjects were generated by our previously trained 3D-IDH network and were reviewed by 2 neuroradiologists for accuracy.20 These tumor masks were used as ground truth for tumor segmentation in the training step. Ground truth whole-tumor masks for methylated and unmethylated MGMT promoter type were labeled with 1's and 2's, respectively (Fig 1). Data preprocessing steps included the following: 1) the Advanced Normalization Tools software package (http://stnava.github.io/ANTs/) affine coregistration21 to the SRI24 T2 template,22 2) skull stripping using the Brain Extraction Tool (BET; http://fsl.fmrib.ox.ac.uk/fsl/fslwiki/BET)23 from FSL,23⇓⇓-26 3) removing radiofrequency inhomogeneity using N4 Bias Field Correction (https://simpleitk.readthedocs.io/en/master/link_N4BiasFieldCorrection_docs.html),27 and 4) normalizing intensity to zero-mean and unit variance. The preprocessing took <5 minutes per dataset.

FIG 1.
FIG 1.

Ground truth whole-tumor masks. Red voxels represent methylated MGMT promoter status (values of 1) and blue voxels represent unmethylated MGMT promoter status (values of 2). The ground truth labels have the same MGMT promoter status for all voxels in each tumor.

Network Details

Transfer learning for determination of MGMT promoter status was implemented using our previously trained 3D-IDH network.20 The decoder part of the network was fine-tuned for a voxelwise dual-class segmentation of the whole tumor, with 1 and 2 representing methylated and unmethylated MGMT promoter types, respectively. The network architecture is shown in Fig 2B. A detailed schematic of the network is provided in the Online Supplemental Data.

FIG 2.
FIG 2.

A, MGMT-net overview. Voxelwise classification of MGMT promoter status is performed to create 2 volumes (methylated and unmethylated MGMT promoter). Volumes are combined using Dual Volume Fusion to eliminate false-positives and generate a tumor-segmentation volume. Majority voting across voxels is used to determine the overall MGMT promoter status. B, Network architecture for MGMT-net. 3D-dense-UNets are used with 7 dense blocks, 3 transition-down (TD) blocks, and 3 transition-up (TU) blocks. Conv indicates convolution layer.

Network Implementation and Cross-Validation

To generalize the network's performance, we performed a 3-fold cross-validation. The dataset of 247 subjects was randomly shuffled and distributed into 3 groups (approximately 82 subjects for each group). Group 1 had 82 subjects (54 methylated, 28 unmethylated), group 2 had 83 subjects (55 methylated, 28 unmethylated), and group 3 had 82 subjects (54 methylated, 28 unmethylated). The 3 groups alternated among training, in-training validation, and held-out testing groups so that each fold of the cross-validation was a new training phase based on a unique combination of the 3 groups. The network uses the in-training validation dataset to evaluate its learning after each training round and updates model parameters to improve performance. However, the network performance is reported only on the held-out testing group for each fold because it is never seen by the network. The group membership for each cross-validation fold is listed in the Online Supplemental Data.

Seventy-five percent overlapping 3D patches (size: 32 × 32 × 32 voxels) were extracted from the training and in-training validation dataset. The patch extraction was performed as a translation in the x-y-z-plane. During training, only patches with at least 1 tumor voxel were included; thus, the number of patches included per training cases varied depending on the size of the tumor. For testing however, the entire image was sampled, including background masked voxels (of value zero). No patch from the same subject was mixed with the training, in-training validation, or testing datasets to prevent the problem of data leakage.28,29 Data augmentation steps included horizontal and vertical flipping, random and translational rotation, the addition of salt and pepper noise, the addition of Gaussian noise, and projective transformation. Additional data augmentation steps included down-sampling images by 50% and 25% (reducing the voxel resolution to 2 and 4 mm3). The data augmentation provided a total of approximately 300,000 patches for training and 300,000 patches for in-training validation for each fold. The networks were implemented using the Tensorflow30 backend engine, the Keras31 Python package, and an Adaptive Moment Estimation optimizer (Adam).32 The initial learning rate was set to 10−5 with a batch size of 15 and maximal epochs of 100 for each fold.

MGMT-net outputs 2 segmentation volumes (V1 and V2), which are combined to generate the voxelwise prediction of methylated and unmethylated MGMT promoter tumor voxels, respectively. The 2 volumes are fused, and the largest connected component (the 3D-connected component algorithm in Matlab [MathWorks]) is obtained as the single tumor-segmentation map. Majority voting over the voxelwise classes of methylated or unmethylated type provided a single MGMT promoter classification for each subject. Tesla V100s, P100, P40, and K80 NVIDIA-GPUs were used to implement the networks. This MGMT promoter determination process is fully automated, and a tumor segmentation map is a natural output of the voxelwise classification approach.

Statistical Analysis

Statistical analysis of the network's performance was performed in Matlab and R statistical and computing software (http://www.r-project.org/). Network accuracies were evaluated using majority voting (ie, a voxelwise cutoff of 50%). The accuracy, sensitivity, specificity, positive predictive value, and negative predictive value of the model for each fold of the cross-validation procedure were calculated using this threshold. Receiver operating characteristic curves for each fold were generated separately. Dice scores were calculated to evaluate the tumor-segmentation performance of the networks. The Dice score calculates the spatial overlap between the ground truth segmentation and the network segmentation.

RESULTS

The network achieved a mean cross-validation testing accuracy of 94.73% across the 3 folds (95.12%, 93.98%, and 95.12% [SD, 0.66%]). Mean cross-validation sensitivity, specificity, positive predictive value, negative predictive value, and area under the curve for the MGMT-net was 96.31% [SD, 0.04%], 91.66% [SD, 2.06%], 95.74% [SD, 0.95%], 92.76% [SD, 0.15%], and 0.93 [SD, 0.03], respectively. The mean cross-validation Dice score for tumor segmentation was 0.82 [SD, 0.008] (Table). The network misclassified 4 cases for fold one, 5 cases for fold 2, and 4 cases for fold three (13 total of 247 subjects). Six subjects were misclassified as unmethylated, and 7, as methylated.

Receiver Operating Characteristic Analysis

The receiver operating characteristic curves for each cross-validation fold for the network are provided in Fig 3. The network demonstrated very good performance with high sensitivities and specificities.

FIG 3.
FIG 3.

Receiver operating characteristic (ROC) analysis for MGMT-net. Separate curves are plotted for each cross-validation fold along with corresponding area under the curve (AUC) value.

Voxelwise Classification

The network is a voxelwise classifier with the tumor segmentation map being a natural output. Figure 4 shows examples of the voxelwise classification for methylated and unmethylated MGMT promoter types, respectively. The volume-fusion procedure was effective in removing false-positives and improving the Dice scores by approximately 6%. We also computed the voxelwise accuracy for the network. The mean voxelwise accuracies were 81.68% [SD, 0.02%] for methylated type and 70.83% [SD, 0.04%] for unmethylated type.

FIG 4.
FIG 4.

A, An example of voxelwise segmentation for a tumor with a methylated MGMT promoter: native T2WI (a), ground truth segmentation (b), and network output after Dual Volume Fusion (c). Red voxels correspond to MGMT methylated class, and blue voxels correspond to MGMT unmethylated class. B, An example of voxelwise segmentation for a tumor with an unmethylated MGMT promoter. The sharp borders visible between methylated and unmethylated types result from the patch-wise classification approach.

Training and Segmentation Times

Fine-tuning the network took approximately 1 week. The trained network took approximately 3 minutes to segment the whole tumor and determine the MGMT status for each subject.

DISCUSSION

We developed a fully-automated, highly accurate, deep learning network for determining the methylation status of the MGMT promoter that outperforms previously reported algorithms.33⇓-35 Our network is able to determine MGMT promoter methylation status from T2WI alone. This eliminates potential failures from image-acquisition artifacts and makes clinical translation straightforward because T2WI is routinely obtained as part of standard clinical brain MR imaging. Previous approaches have required multicontrast input, which can be compromised due to patient motion from lengthier examination times and the need for gadolinium contrast. Obviating the need for intravenous contrast makes our algorithm applicable to patients with contrast allergies and renal failure. Compared with previously published algorithms, our methodology is fully automated and uses minimal preprocessing. The time required for the MGMT-net to segment the whole tumor and predict the MGMT promoter methylation status for 1 subject is approximately 3 minutes on a K80 or P40 NVIDA-GPU.

Other groups have also proposed deep learning methods for noninvasive, image-based MGMT molecular profiling, but each of these has several limitations. Korfiatis et al9 implemented a 2D-based slice-wise network, pre-selecting only cases of glioblastoma multiforme for training and prediction. While they achieved a high slice-wise accuracy, their average subject-wise MGMT prediction accuracy was only 90%. Most important, in clinical practice, the tumor grade is unknown a priori. Thus, the approach of Korfiatis et al is a nonviable clinical method from the outset. Our approach of using a mix of low-grade and high-grade gliomas is a better approximation of the real-world clinical workflow in which tissue is not yet available.

Similar to the work of Korfiatis et al, Chang et al35 also implemented a 2D-network, but instead used a case mix like ours (low-grade and high-grade gliomas from the TCIA/TCGA). However, they were only able to achieve an MGMT prediction accuracy of 83% (range, 76%–88%), and their network required tumor presegmentation. Our algorithm far outperformed the approach of Chang et al on a similar dataset without the need for presegmentation. Additionally, it is unclear whether 2D algorithms of either Korfiatis et al9 or Chang et al35 addressed the issue of “data leakage.”28,29 This is a potentially significant limitation for 2D networks that can occur during the slice-wise randomization process if different slices of the same tumor from the same subject are mixed among training, validation, and testing datasets. Unless this is explicitly addressed during the slice-randomization procedure, the reported accuracies can be upwardly biased. Our approach outperforms all prior reports on noninvasive determination of MGMT status and is the first to achieve tissue-level performance, representing a milestone in the clinical viability of MR imaging–based MGMT promoter status prediction.

The higher performance achieved by our network compared with previous image-based classification studies can be explained by several factors. The dense connections in our 3D network architecture are easier to train, carry information from the previous layers to the following layers, and can reduce over-fitting.36,37 3D networks also interpolate between slices to maintain interslice information more accurately. The Dual Volume Fusion postprocessing step improved the Dice scores by approximately 6% by eliminating extraneous voxels not connected to the tumor. Our approach also uses voxelwise classifiers and provides a classification for each voxel in the image. These steps provide simultaneous single-label tumor segmentation. The cross-validation single-label whole-tumor segmentation performance for the MGMT network provided excellent Dice scores of 0.82 [SD, 0.008].

The ability to determine MGMT promoter methylation status on the basis of MR images alone is clinically significant because it helps determine whether the glioma will be susceptible to temozolomide (TMZ). Alkylating agents such as temozolomide damage DNA by methylating the oxygen at position 6 of the guanine nucleotide (O6-methylguanine). The process by which many DNA repair enzymes remove O6-methylguanine, results in DNA breaks, culminating in cell death. However, MGMT works differently by restoring the normal guanine residue and rescuing the glioma cell. Therefore, MGMT activity leads to resistance to therapy. Methylation of the MGMT promoter leads to inactivation of MGMT and loss of resistance of glioma cells to alkylating agents. The MGMT protein is encoded on the long arm of chromosome 10 at position 26 (10q26). Transcription of the MGMT gene is regulated by several promoters.29

Although incompletely understood, at least 9 specific regions within the promoter's gene determine whether a cell will express or not express MGMT.29 However, some regions have been shown to be more important for loss of MGMT expression.38 In the clinical setting, methods for determining MGMT methylation focus on these regions in the promoter gene. The 4 most prevalent methods to detect MGMT methylation are the following: immunohistochemistry, pyrosequencing, quantitative methylation-specific polymerase chain reaction (PCR), and methylation-specific PCR. Pyrosequencing is considered the theoretic criterion standard but is not readily available, and although it is quantitative, there is no agreement on what cutoff values to use when determining MGMT promoter methylation status.30 Therefore, although it is not quantitative, methylation-specific PCR is the most widely used method.39 Additionally, most centers perform MGMT methylation detection on formalin-fixed or paraffin-embedded tissue specimens. These methods have several limitations. Evaluating multiple different methylation sites is technically challenging on a single tissue specimen.39 Tumor heterogeneity poses a substantial limitation of these methods because sampling bias can lead to inaccurate determinations. The presence of hemorrhage, necrosis, or nonmalignant cells contaminates the specimen.39 Therefore, some institutions mandate that at least 50% of the sample to be analyzed contains tumor cells. Prior to PCR, several tissue-processing steps are required. Bisulfite treatment is the most critical step because it will produce the modified DNA that will be used for PCR; however, it also degrades the amount of DNA available, and incomplete treatment can lead to false-positive results.39 The reported sensitivity and specificity of methylation-specific PCR is 91% and 75%, respectively, while the reported sensitivity and specificity of pyrosequencing is 78% and 90%.32

Our noninvasive, MR imaging–based deep learning algorithm outperformed these methods with a sensitivity and specificity of 96.3% and 91.6%, respectively. The overall determination of MGMT promoter methylation status is based on the majority voxels in the tumor. Given the variability in the cutoff values for pyrosequencing-based detection, we performed a Youden statistical index analysis to determine whether the optimal cutoff for our deep learning algorithm was different from majority voting (>50%). The analysis demonstrated that maximum accuracy, sensitivity, specificity, positive predictive value, and negative predictive value were obtained at an optimal cutoff of 50%, the same as majority voting.

Our algorithm was trained on ground truth obtained from the TCGA data base. TCGA uses the Infinium Methylation Assay (https://www.illumina.com/science/technology/microarray/infinium-methylation-assay.html) to determine MGMT promoter methylation status.40⇓-42 Infinium Methylation Assays are an immunofluorescence method that uses next-generation high-throughput microchip arrays and probes. While these methods have been reported to be more sensitive and specific than the most widely available clinical assays, they require pre-existing probes to detect specific methylation sites.42 The sensitivity and specificity values change depending on the probe and analytic model used to interpret the results.42 The sensitivities for the best probes range from 87.5% to 90.6%, while the specificity is 94.4%.42 The overall accuracy of these probes with an optimized analytic model ranges from 91.24% to 93.6%.34 The accuracy of the commercially available Infinium Methylation Assay with the best analytic model is 92%.34 Our algorithm outperforms this assay with a mean cross-validation testing accuracy of 94.73%. While the algorithm appears to outperform the ground truth, there are additional factors that need to be considered for this dataset. The TCGA data base used very stringent tissue screening before molecular testing, including review of tissue to ensure a minimum of 80% tumor nuclei and a maximum of 50% necrosis with additional quality-control measurements of the extracted DNA and RNA before analyses. Additionally, the MGMT determinations made in the TCGA data base were verified by a secondary test.43 Thus, the reported accuracy of the Infinium Methylation Assay is not necessarily comparable with the accuracy in TCIA/TCGA datasets. It is also possible that the algorithm learns features that allow it to perform better than the single-site tissue-biopsy sample ground truth performance because the algorithm “samples” the entire tumor and learns imaging features that are specific to MGMT mutation.

Tissue-based methods for determining MGMT promoter methylation status remain a complex, multistep process that is susceptible to failure and inaccuracy even after an adequate tissue sample has been obtained. Thus, the ability to determine MGMT promoter methylation status on the basis of routine T2WI alone is highly desirable. Additionally, because our algorithm was trained and evaluated on the multi-institutional TCIA database, it is a better representative of algorithm robustness, real-world performance, and potential clinical use than the previously reported methods.25

The algorithm misclassified 13 cases: Six subjects were misclassified as unmethylated, and 7, as methylated. Despite these misclassifications, our network achieved a mean cross-validation testing accuracy of 94.73%, which is higher than that for the methylation-specific PCR, pyrosequencing (PYR), and Infinium Methylation Assays.42 While these tissue-based methods require an invasive procedure and subsequent tissue processing for at least 48 hours, our deep learning algorithm can segment the entire glioma and determine MGMT promoter methylation status in 3 minutes. The deep learning algorithm can also be fine-tuned to variations in institutional MR imaging scanners, while other tissue-based methods currently lack standardization as mentioned above.

The limitations of our study are that deep learning studies require large amounts of data and the relative number of subjects with MGMT promoter methylation is small in the TCGA database. While the number of subjects may seem small, we used a patch-based algorithm with data augmentation, which provided well over 300,000 samples (patches) for training and validation. Additionally, acquisition parameters and imaging vendor platforms vary across imaging centers that contribute data, though this may also be a regarded as a desirable aspect for the generalizability of the approach. Our current classification approach uses a largest connected component step to limit false-positives. As a consequence, multifocal tumors represent a potential limitation. Despite these caveats, our algorithm demonstrated high accuracy in determining MGMT promoter methylation status approaching tissue-level performance.

CONCLUSIONS

We demonstrate high accuracy in determining MGMT promoter methylation status using only T2WI. This represents an important milestone toward using MR imaging to predict glioma histology, prognosis, and appropriate treatment.

ACKNOWLEDGMENTS

We thank Yin Xi, PhD, a statistician, for help with the receiver operating characteristics and areas under the curve.

Footnotes

  • This work was supported by National Institutes of Health/National Cancer Institute U01CA207091 (A.J.M., J.A.M.).

  • Disclosures: Chandan Ganesh Bangalore Yogananda—UNRELATED: Employment: University of Texas Southwestern Medical Center. Baowei Fei—RELATED: Grant: National Institutes of Health, Comments: This research was supported, in part, by the US National Institutes of Health grants (R01CA156775, R01CA204254, R01HL140325, and R21CA231911) and by the Cancer Prevention and Research Institute of Texas grant RP190588.* Ananth J. Madhuranthakam—RELATED: Grant: National Institutes of Health/National Cancer Institute, Comments: U01CA207091.* Joseph A. Maldjian—RELATED: Grant: National Institutes of Health/National Cancer Institute grant*; UNRELATED: Consultancy: BioClinica, Comments: blinded clinical trial reader. *Money paid to the institution.

Indicates open access to non-subscribers at www.ajnr.org

References

  1. 1.↵
    2. Hegi ME,
    3. Diserens AC,
    4. Gorlia T, et al
    . MGMT gene silencing and benefit from temozolomide in glioblastoma. N Engl J Med 2005;352:997–1003 doi:10.1056/NEJMoa043331 pmid:15758010
    CrossRefPubMed
  2. 2.↵
    2. Stupp R,
    3. Hegi ME,
    4. Mason WP, et al
    . European Organisation for Research and Treatment of Cancer Brain Tumour and Radiation Oncology Groups. Effects of radiotherapy with concomitant and adjuvant temozolomide versus radiotherapy alone on survival in glioblastoma in a randomised phase III study: 5-year analysis of the EORTC-NCIC trial. Lancet Oncol 2009;10:459–66 doi:10.1016/S1470-2045(09)70025-7 pmid:19269895
    CrossRefPubMed
  3. 3.↵
    2. Chen R,
    3. Smith-Cohn M,
    4. Cohen AL, et al
    . Glioma subclassifications and their clinical significance. Neurotherapeutics 2017;14:284–97 doi:10.1007/s13311-017-0519-x pmid:28281173
    CrossRefPubMed
  4. 4.↵
    2. Suh CH,
    3. Kim HS,
    4. Jung SC, et al
    . Clinically relevant imaging features for MGMT promoter methylation in multiple glioblastoma studies: a systematic review and meta-analysis. AJNR Am J Neuroradiol 2018;39:1439–45 doi:10.3174/ajnr.A5711 pmid:30002055
    Abstract/FREE Full Text
  5. 5.↵
    2. Drabycz S,
    3. Roldan G,
    4. de Robles P, et al
    . An analysis of image texture, tumor location, and MGMT promoter methylation in glioblastoma using magnetic resonance imaging. Neuroimage 2010;49:1398–1405 doi:10.1016/j.neuroimage.2009.09.049 pmid:19796694
    CrossRefPubMed
  6. 6.↵
    2. Moon WJ,
    3. Choi JW,
    4. Roh HG, et al
    . Imaging parameters of high grade gliomas in relation to the MGMT promoter methylation status: the CT, diffusion tensor imaging, and perfusion MR imaging. Neuroradiology 2012;54:555–63 doi:10.1007/s00234-011-0947-y pmid:21833736
    CrossRefPubMed
  7. 7.↵
    2. Ahn SS,
    3. Shin NY,
    4. Chang JH, et al
    . Prediction of methylguanine methyltransferase promoter methylation in glioblastoma using dynamic contrast-enhanced magnetic resonance and diffusion tensor imaging. J Neurosurg 2014;121:367–73 doi:10.3171/2014.5.JNS132279 pmid:24949678
    CrossRefPubMed
  8. 8.↵
    2. Kanas VG,
    3. Zacharaki EI,
    4. Thomas GA, et al
    . Learning MRI-based classification models for MGMT methylation status prediction in glioblastoma. Comput Methods Programs Biomed 2017;140:249–57 doi:10.1016/j.cmpb.2016.12.018 pmid:28254081
    CrossRefPubMed
  9. 9.↵
    2. Korfiatis P,
    3. Kline TL,
    4. Coufalova L, et al
    . MRI texture features as biomarkers to predict MGMT methylation status in glioblastomas. Med Phys 2016;43:2835–44 doi:10.1118/1.4948668 pmid:27277032
    CrossRefPubMed
  10. 10.↵
    2. Sasaki T,
    3. Kinoshita M,
    4. Fujita K, et al
    . Radiomics and MGMT promoter methylation for prognostication of newly diagnosed glioblastoma. Sci Rep 2019;9:14435 doi:10.1038/s41598-019-50849-y pmid:31594994
    CrossRefPubMed
  11. 11.↵
    2. Wei J,
    3. Yang G,
    4. Hao X, et al
    . A multi-sequence and habitat-based MRI radiomics signature for preoperative prediction of MGMT promoter methylation in astrocytomas with prognostic implication. Eur Radiol 2019;29:877–88 doi:10.1007/s00330-018-5575-z pmid:30039219
    CrossRefPubMed
  12. 12.↵
    2. Yogananda CG,
    3. Shah BR,
    4. Yu FF, et al
    . A novel fully automated MRI-based deep-learning method for classification of 1p/19q co-deletion status in brain gliomas. Neurooncol Adv 2020;2:vdaa066 doi:10.1093/noajnl/vdaa066 pmid:32705083
    CrossRefPubMed
  13. 13.↵
    2. Yogananda CG,
    3. Shah BR,
    4. Vejdani-Jahromi M, et al
    . A novel fully automated MRI-based deep-learning method for classification of IDH mutation status in brain gliomas. Neuro Oncol 2019;22:402–11 doi:10.1093/neuonc/noz199 pmid:31637430
    CrossRefPubMed
  14. 14.↵
    2. Clark K,
    3. Vendt B,
    4. Smith K, et al
    . The Cancer Imaging Archive (TCIA): maintaining and operating a public information repository. J Digit Imaging 2013;26:1045–57 doi:10.1007/s10278-013-9622-7 pmid:23884657
    CrossRefPubMed
  15. 15.↵
    2. Puchalski RB,
    3. Shah N,
    4. Miller J, et al
    . An anatomic transcriptional atlas of human glioblastoma. Science 2018;360:660–63 doi:10.1126/science.aaf2666 pmid:29748285
    Abstract/FREE Full Text
  16. 16.↵
    2. Erickson B,
    3. Akkus Z,
    4. Sedlar J, et al
    . Data From LGG-1p19qDeletion. The Cancer Imaging Archive. 2017;76. https://doi.org/10.7937/K9/TCIA.2017.dwehtz9v. Accessed November 30, 2017
  17. 17.↵
    2. Ceccarelli M,
    3. Barthel FP,
    4. Malta TM, et al
    . TCGA Research Network. Molecular profiling reveals biologically discrete subsets and pathways of progression in diffuse glioma. Cell 2016;164:550–63 doi:10.1016/j.cell.2015.12.028 pmid:26824661
    CrossRefPubMed
  18. 18.↵
    2. Menze BH,
    3. Jakab A,
    4. Bauer S, et al
    . The multimodal Brain Tumor Image Segmentation Benchmark (BRATS). IEEE Trans Med Imaging 2015;34:1993–2024 doi:10.1109/TMI.2014.2377694 pmid:25494501
    CrossRefPubMed
  19. 19.↵
    2. Bakas S,
    3. Akbari H,
    4. Sotiras A, et al
    . Advancing The Cancer Genome Atlas glioma MRI collections with expert segmentation labels and radiomic features. Sci Data 2017;4:170117 doi:10.1038/sdata.2017.117 pmid:28872634
    CrossRefPubMed
  20. 20.↵
    2. Yogananda CG,
    3. Shah BR,
    4. Vejdani-Jahromi M, et al
    . A novel fully automated MRI-based deep learning method for classification of IDH mutation status in brain gliomas. Neuro Oncol 2020;22:402–11 doi:10.1093/neuonc/noz199 pmid:31637430
    CrossRefPubMed
  21. 21.↵
    2. Avants BB,
    3. Tustison NJ,
    4. Song G, et al
    . A reproducible evaluation of ANTs similarity metric performance in brain image registration. Neuroimage 2011;54:2033–44 doi:10.1016/j.neuroimage.2010.09.025 pmid:20851191
    CrossRefPubMed
  22. 22.↵
    2. Rohlfing T,
    3. Zahr NM,
    4. Sullivan EV, et al
    . The SRI24 multichannel atlas of normal adult human brain structure. Hum Brain Mapp 2010;31:798–819 doi:10.1002/hbm.20906 pmid:20017133
    CrossRefPubMed
  23. 23.↵
    2. Smith SM
    . Fast robust automated brain extraction. Hum Brain Mapp 2002;17:143–55 doi:10.1002/hbm.10062 pmid:12391568
    CrossRefPubMed
  24. 24.↵
    2. Smith SM,
    3. Jenkinson M,
    4. Woolrich MW, et al
    . Advances in functional and structural MR image analysis and implementation as FSL. Neuroimage 2004;23(Suppl 1):S208–19 doi:10.1016/j.neuroimage.2004.07.05]1 pmid:15501092
    CrossRefPubMed
  25. 25.↵
    2. Woolrich MW,
    3. Jbabdi S,
    4. Patenaude B, et al
    . Bayesian analysis of neuroimaging data in FSL. Neuroimage 2009;45:S173–86 doi:10.1016/j.neuroimage.2008.10.055 pmid:19059349
    CrossRefPubMed
  26. 26.↵
    2. Jenkinson M,
    3. Beckmann CF,
    4. Behrens TE, et al
    . FSL. Neuroimage 2012;62:782–90 doi:10.1016/j.neuroimage.2011.09.015 pmid:21979382
    CrossRefPubMed
  27. 27.↵
    2. Tustison NJ,
    3. Cook PA,
    4. Klein A, et al
    . Large-scale evaluation of ANTs and FreeSurfer cortical thickness measurements. Neuroimage 2014;99:166–79 doi:10.1016/j.neuroimage.2014.05.044 pmid:24879923
    CrossRefPubMed
  28. 28.↵
    2. Petrick N,
    3. Mori K
    2. Wegmayr VA,
    3. Aitharajus S,
    4. Buhmann J, et al
    . Classification of brain MRI with big data and deep 3D convolutional neural networks. In: Petrick N, Mori K, eds. Medical Imaging 2018: Computer-Aided Diagnosis. SPIE Proceedings 2018;1057501 doi:10.1117/12.2293719
    CrossRef
  29. 29.↵
    2. Feng X,
    3. Yang J,
    4. Lipton ZC, et al
    . Deep larning on MRI affirms the prominence of the hippocampal formation in Alzheimer's disease classification. bioRxiv 2018 https://www.biorxiv.org/content/10.1101/456277v1.full.pdf. Accessed April 8, 2020
  30. 30.↵
    2. Abadi M,
    3. Barham P,
    4. Chen J, et al
    . Tensorflow: a system for large-scale machine learning. 12th (USENIX) symposium on operating systems design and implementation (OSDI 16); 2016:265–283. May 27, 2016. https://arxiv.org/a bs/1605.08695. Accessed April 24, 2018
  31. 31.↵
    2. Chollet F. Keras
    : The python deep learning library. Astrophysics Source Code Library 2018. https://keras.io/. Accessed February 10, 2018
  32. 32.↵
    2. Kingma DP,
    3. Ba JL
    . Adam: a method for stochastic optimization. In: Proceedings of the International Conference on Learning Representations, San Diego, California. May 7–9, 2015
  33. 33.↵
    2. Korfiatis P,
    3. Kline TL,
    4. Lachance DH, et al
    . Residual deep convolutional neural network predicts MGMT methylation status. J Digit Imaging 2017;30:622–28 doi:10.1007/s10278-017-0009-z pmid:28785873
    CrossRefPubMed
  34. 34.↵
    2. Han L,
    3. Kamdar MR
    . MRI to MGMT: predicting methylation status in glioblastoma patients using convolutional recurrent neural networks. Pac Symp Biocomput 2018;23:331–42 pmid:29218894
    PubMed
  35. 35.↵
    2. Chang P,
    3. Grinband J,
    4. Weinberg BD, et al
    . Deep-learning convolutional neural networks accurately classify genetic mutations in gliomas. AJNR Am J Neuroradiol 2018;39:1201–07 doi:10.3174/ajnr.A5667 pmid:29748206
    Abstract/FREE Full Text
  36. 36.↵
    2. Jégou S,
    3. Drozdzal M,
    4. Vazquez D, et al
    . The one hundred layers tiramisu: fully convolutional densenets for semantic segmentation. In: Proceedings of the 2017 IEEE Conference on Computer Vision and Pattern Recognition Workshops. 2017:11–19. https://ieeexplore.ieee.org/xpl/conhome/8014302/proceeding. Accessed March 28, 2018
  37. 37.↵
    2. Wang G,
    3. Li W,
    4. Ourselin S, et al
    . Automatic brain tumor segmentation based on cascaded convolutional neural networks with uncertainty estimation. Front Comput Neurosci 2019;13:56 doi:10.3389/fncom.2019.00056 pmid:31456678
    CrossRefPubMed
  38. 38.↵
    2. Everhard S,
    3. Tost J,
    4. El Abdalaoui H, et al
    . Identification of regions correlating MGMT promoter methylation and gene expression in glioblastomas. Neuro Oncol 2009;11:348–56 doi:10.1215/15228517-2009-001 pmid:19224763
    CrossRefPubMed
  39. 39.↵
    2. Cankovic M,
    3. Nikiforova MN,
    4. Snuderl M, et al
    . The role of MGMT testing in clinical practice: a report of the association for molecular pathology. J Mol Diagn 2013;15:539–55 doi:10.1016/j.jmoldx.2013.05.011 pmid:23871769
    CrossRefPubMed
  40. 40.↵
    2. Estival A,
    3. Sanz C,
    4. Ramirez JL, et al
    . Pyrosequencing versus methylation-specific PCR for assessment of MGMT methylation in tumor and blood samples of glioblastoma patients. Sci Rep 2019;9:11125 doi:10.1038/s41598-019-47642-2
    CrossRef
  41. 41.↵
    2. Poulin M,
    3. Zhou JY,
    4. Yan L, et al
    . Pyrosequencing methylation analysis. Methods Mol Biol 2018;1856:283–96 doi:10.1007/978-1-4939-8751-1_17 pmid:30178259
    CrossRefPubMed
  42. 42.↵
    2. Bady P,
    3. Sciuscio D,
    4. Diserens AC, et al
    . MGMT methylation analysis of glioblastoma on the Infinium methylation BeadChip identifies two distinct CpG regions associated with gene silencing and outcome, yielding a prediction model for comparisons across datasets, tumor grades, and CIMP-status. Acta Neuropathol 2012;124:547–60 doi:10.1007/s00401-012-1016-2 pmid:22810491
    CrossRefPubMed
  43. 43.↵
    Cancer Genome Atlas Research Network. Comprehensive genomic characterization defines human glioblastoma genes and core pathways. Nature 2008;455:1061–68 doi:10.1038/nature07385 pmid:18772890
    CrossRefPubMed
  • Received June 24, 2020.
  • Accepted after revision November 21, 2020.
Back to top

In this issue

Advertisement
Print
Download PDF
Email Article
Citation Tools
MRI-Based Deep-Learning Method for Determining Glioma MGMT Promoter Methylation Status
C.G.B. Yogananda, B.R. Shah, S.S. Nalawade, G.K. Murugesan, F.F. Yu, M.C. Pinho, B.C. Wagner, B. Mickey, T.R. Patel, B. Fei, A.J. Madhuranthakam, J.A. Maldjian
American Journal of Neuroradiology May 2021, 42 (5) 845-852; DOI: 10.3174/ajnr.A7029
del.icio.us logo Twitter logo Facebook logo Mendeley logo
Purchase

Jump to section

Related Articles

Cited By...

More in this TOC Section

Similar Articles

Advertisement