- 1Department of Radiology, Third Affiliated Hospital of Southern Medical University (Academy of Orthopedics Guangdong), Guangzhou, China
- 2Imaging Department of Guangdong 999 Brain Hospital, Guangzhou, China
- 3GE Healthcare, Guangzhou, China
- 4Department of Neuro-Oncology of Guangdong 999 Brain Hospital, Guangzhou, China
Objectives: To explore the feasibility of predicting overall survival (OS) of patients with midline glioma using multi-parameter magnetic resonance imaging (MRI) features.
Methods: Data of 84 patients with midline gliomas were retrospectively collected, including 40 patients with OS > 12 months (28 cases were adults, 14 cases were H3 K27M-mutation) and 44 patients with OS < 12 months (29 cases were adults, 31 cases were H3 K27M-mutation). Features were extracted from the largest slice of tumors, which were manually segmented on T2-weighted (T2w), T2 fluid-attenuated inversion recovery (T2 FLAIR), and contrast-enhanced T1-weighted (T1c) images. Data were randomly divided into training (70%) and test cohorts (30%) and normalized and standardized using Z-scores. Feature dimensionality reduction was performed using the variance method and maximum relevance and minimum redundancy (mRMR) algorithm. We used the logistic regression algorithm to construct three models for T2w, T2 FLAIR, and T1c images as well as one combined model. The test cohort was used to evaluate the models, and receiver operating characteristic (ROC) curves, areas under the curve (AUCs), sensitivity, specificity, and accuracy were calculated. The nomogram of the combined model was built and evaluated using a calibration curve. Decision curve analysis (DCA) was used to evaluate the clinical application value of the four models.
Results: A total of 1,316 features were extracted from T2w, T2 FLAIR, and T1c images, respectively. And then the best non-redundant features were selected from the extracted features using the variance method and mRMR. Finally, five features were extracted each from T2w, T2 FLAIR, and T1c images, and 12 features were extracted for the combined model. Four models were established using the optimal features. In the test cohort, the combined model performed the best out of all models. The AUCs of the T2w, T2 FLAIR, T1c, and combined models were 0.73, 0.78, 0.74, and 0.87, respectively, and accuracies were 0.72, 0.76, 0.72, and 0.84, respectively. The ROC curves and DCA showed that the combined model had the highest efficiency and most favorable clinical benefits.
Conclusion: The combined radiomics model based on multi-parameter MRI features provided a reliable non-invasive method for the prognostic prediction of midline gliomas.
Introduction
Diffuse midline gliomas (DMGs) with an H3 K27M mutation are a group of tumors newly defined in 2016 and amended to DMGs with an H3 K27-altered in 2021 (1, 2). They occur commonly in the brainstem of children, although they are not uncommon in adults. The median survival time of patients with a DMG with H3 K27M mutation is shorter than that of patients without H3 K27M mutation (3). Furthermore, H3 K27M mutation not only have prognostic implications but may also represent a potential new immunotherapeutic target (4). However, because of the high risk of operating on the midline area, most patients only undergo radiotherapy or chemotherapy (5, 6). Without information regarding H3 K27-altered, predicting the survival of patients with midline gliomas is challenging. Therefore, it is of significant clinical value to develop a non-invasive method to evaluate the survival of patients with midline gliomas.
In clinical practice, magnetic resonance imaging (MRI) is one of the most valuable methods for evaluating the survival of patients with midline gliomas before and during treatment. MRI findings of tumors often help clinicians adjust treatment strategies. Furthermore, the World Health Organization's classification of tumors of the central nervous system began integrating molecular and genetic profiling in 2016, with gradual improvements being made continuously. However, the relationship between macroscopic MRI findings and molecular subtyping of tumors remains uncertain. H3 K27M mutant patients has shorter survival time than those of wild type, and several studies have shown no significant correlation between H3 K27M mutant- and wild-type DMGs in terms of tumor necrosis, patterns of enhancement, edema, infiltrative features, or diffusion characteristics (3, 7, 8). This suggests that macroscopic MRI findings, such as tumor location, size, signal intensity, contrast enhancement, and advanced MRI techniques, such as diffusion-weighted imaging, do not currently meet the clinical needs of evaluation of DMG survival. On the other hand, according to the 2021 WHO classification of tumors of the central nervous system, DMGs, H3 K27-altered include H3 wild type with enhancer of zeste homolog inhibitory protein (EZHIP) overexpression and epidermal growth factor receptor (EGFR) mutation in addition to H3 K27M mutation (2). Therefore, the prediction of H3 K27M mutation alone in the past cannot accurately evaluate the prognosis of patients with DMG. Thus, there is an urgent need for better and more advanced tools to explore biological information hidden in traditional MR images, which would have substantial clinical implications, such as avoiding inappropriate treatments.
One such tool is radiomics, which involves the non-invasive extraction of a large amount of quantitative information from medical images that cannot be perceived by human vision. To date, radiomics has been demonstrated to be valuable in predicting the molecular genetic characteristics, tumor grading, differential diagnosis, and prognostic evaluation of gliomas (6, 9–14). Numerous studies have accurately predicted the survival of gliomas using radiomics. Prasanna et al. (15) found that radiomic features from the peritumoral brain parenchyma on routine pre-operative MRI can predict long- vs. short-term survival in glioblastomas. Senders et al. (16) have trained fifteen statistical and machine learning algorithms based on 13 demographic, socioeconomic, clinical, and radiographic features to predict overall survival, 1-year survival status of glioblastomas, and found the accelerated failure time model demonstrated superior performance compared to Cox proportional hazards regression and other machine learning algorithms. Therefore, we investigated the feasibility of predicting the overall survival (OS) of patients with midline gliomas using machine learning to offer an alternative method to non-invasively predict prognoses of midline glioma patients.
Materials and Methods
Patients
This study was approved by the hospital ethics committee, and the requirement for obtaining informed consent was waived. We retrospectively collected data of patients with a midline glioma treated at Guangdong 999 Brain Hospital between January 2017 and October 2020. The inclusion criteria were: (1) glioma located in the midline area (e.g., vermis, brainstem, thalamus, corpus callosum, sellar region, and basal ganglia); (2) patients had undergone surgical resection or biopsy, and the pathological diagnosis was certain; and (3) patients underwent T2-weighted (T2w), T2 fluid-attenuated inversion recovery (T2 FLAIR), and contrast-enhanced T1-weighted (T1c) imaging at our hospital before the operation. Exclusion criteria were: (1) missing T2w, T2 FLAIR, or T1c images; (2) uncertain pathological diagnosis; and (3) poor image quality. Poor image quality is defined as obvious motion artifacts; the vascular pulsation artifact in the posterior fossa was obvious after enhancement; and (4) patients who had received related preoperative treatments, such as radiotherapy or chemotherapy; (5) patients who had been followed up for < 12 months.
Initially, 161 cases were included, of whom 77 were excluded; 20 had received treatment before surgery, 11 had uncertain pathology, 12 had incomplete imaging data, 6 had poor image quality, and 28 were followed up for < 12 months. Finally, a total of 84 cases were included. The OS of 40 cases was >12 months (long-term group) and the OS of 44 cases was <12 months (short-term group) (Figure 1).
T1w, T2w, T2 FLAIR, and T1c Images Acquisition
All patients were scanned using a GE HDxt 3.0T MR (General Electric, Milwaukee, WI, USA) or Philips Intera 1.5T MR (Royal Philips Electronics, Amsterdam, Netherlands) scanner. Imaging sequences included T1w, T2w, T2 FLAIR, and T1c sequences. The 3.0T MRI scan parameters were: repetition time (TR) 1,900 ms, echo time (TE) 24 ms, and inversion time (TI) 780 ms for T1w imaging; TR 4,480 ms and TE 120 ms for T2w imaging; TR 9,480 ms, TE 120 ms, and TI 2,300 ms for T2 FLAIR imaging. The T1w, T2w, and T2 FLAIR imaging sequences on the 3.0T scanner used the following parameters: field of view (FOV): 240 × 240 mm, matrix: 256 × 256, slice thickness: 5.5 mm, slice gap: 1.0 mm, number of excitations (NEX): 1. The 1.5T MRI scan parameters were: TR 488 ms, TE 15 ms, and matrix: 152 × 121 for T1w imaging; TR 3,980 ms, TE 110 ms, and matrix: 230 × 130 for T2w imaging; TR 6,000 ms, TE 120 ms, TI 2,000 ms, and matrix: 192 × 115 for T2 FLAIR imaging. The T1w, T2w, and T2 FLAIR imaging sequences on the 1.5T scanner used the following parameters: FOV: 230 × 182 mm, slice thickness: 5.5 mm, slice gap: 1.0 mm, NEX: 1. T1c imaging on the 1.5T and 3.0T MRI scanners was performed after the administration of Gadolinium-Diethylenetriaminepentaacetic Acid (Gd-DTPA) (Kangchen Company, Guangzhou, China) or Gd diamine (GE Pharmaceuticals, USA) at a dose of 0.1 mmol/kg body weight.
Tumor Segmentation
All images were saved as Digital Imaging and Communications in Medicine (DICOM) and imported into the ITK-SNAP software (version 3.8.0, http://www.itksnap.org) for segmentation. Some studies have shown that one slice (2D) with the largest cross-section of the tumor and the entire tumor volume (3D) for segmentation has comparable diagnostic performance, and there are also some studies that use one slice for tumor segmentation (17–20). So, in this study, the regions of interest (ROIs) were manually segmented (Figure 2) on the axial slice with the largest cross-section of the tumors on T2w, T2 FLAIR, and T1c images by two authors (S Wu and P He, with 4 and 6 years of diagnostic experience in neuroradiology, respectively) who were blinded to histological results.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/866274/fneur-13-866274-HTML/image_m/fneur-13-866274-g002.jpg)
Figure 2. MRI feature extraction and machine learning flow chart. By manually segmenting the tumors, we extracted 1,316 imaging features from masked preoperative T2w, T2 FLAIR, and T1c images. The variance method and maximum relevance and minimum redundancy were used to select the best features, and the four models were established. The test cohort was then used to verify the model.
The final ROIs were calculated using the overlapping segmentations generated by the two authors. If the overlap rate was <90%, the ROI was defined by L Cheng, who had >10 years of diagnostic experience in neuroradiology (16, 21). The ROIs included cystic and necrotic components which may reflect tumor heterogeneity.
Imaging Feature Extraction
A total of 1,316 texture parameters were extracted from T2w, T2 FLAIR, and T1c images using the open-source package Pyradiomics (version 3.0.1, https://github.com/Radiomics/pyradiomics) on Python (version 3.5.6, https://www.python.org) and included first-order statistical features, shape features, gray level co-occurrence matrix, gray level dependence matrix, gray level run length matrix, gray level size zone matrix, and neighboring gray tone dependence matrix (22). The median and block methods were used to replace missing values and outliers.
Data Cohort Division and Preprocessing
Data were randomly divided into training (70%; 59 patients including the long-term group 31 patients and the long-term group 28 patients) and test (30%; 25 patients including the long-term group13 patients and the long-term group 12 patients) cohorts. All features extracted from the ROIs were normalized using Z-scores.
Feature Dimensionality Reduction and Selection of Optimal Features
The variance method was used to calculate the variance of features extracted from imaging. When the variance of the feature was 0, the feature was removed. The maximum relevance and minimum redundancy (mRMR) method was used to select the best non-redundant features from the extracted features.
Radiomics Feature Model and Evaluation
Logistic regression was used to establish the model. According to the best features of the T2w, T2 FLAIR, and T1c images, three models for T2w, T2 FLAIR, and T1c images and one combined model (T2w + T2 FLAIR + T1c) were constructed. Receiver operating characteristic (ROC) curve analysis and decision curve analysis (DCA) were used to assess the performance of the models (Figure 2).
Statistical Analysis
All statistical analyses were performed using R (version 3.5.1, http://www.R-project.org) and Python (version 3.5.6, https://www.python.org/). Chi-squared or Fisher's exact tests were used for nominal variables. Kruskal–Wallis H tests were used for ordinal variables. Student's t-tests were used for continuous variables. The variance method and mRMR were used for feature dimensionality reduction. Logistic regression was used to establish the model and construct a nomogram. ROC curve analysis was used to assess the established models according to accuracy (AC), area under the curve (AUC), sensitivity, and specificity values. DCA was used to evaluate the application value of the models. A p < 0.05 was considered statistically significant.
Results
Patient Characteristics
A total of 84 patients with a midline glioma were enrolled in the study (62 men and 22 women; mean age 27.73 ± 16.36 years; range 1–63 years). The clinical characteristics of the long-term and short-term groups are shown in Table 1, Figure 3. There was a significant difference in radiotherapy between the long-term group and the short-term group (p = 0.002). The patients with H3 K27M mutation in short-term group outnumbered those in long-term group significantly (p = 0.001). There was also a significant difference in OS between the short-term and long-term groups (5.16 ± 5.38 months vs. 25.34 ± 9.46 months, p < 0.001). The clinical characteristics of the training and test cohorts are shown in Table 2. The patients who underwent radiotherapy in the training outnumbered those in the test cohort significantly (p = 0.0028).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/866274/fneur-13-866274-HTML/image_m/fneur-13-866274-g003.jpg)
Figure 3. Images of patients with midline gliomas. T2w (A), T2 FLAIR (B), and T1c (C) images of a 34-year-old man with a pontine diffuse midline glioma with a H3 K27M mutation (WHO 4), and histopathology of diffuse astrocytoma. He underwent a biopsy and was treated with radiotherapy and chemotherapy. The overall survival of the patient was 4.4 months. T2w (D), T2 FLAIR (E), and T1c (F) images of a 33-year-old man with corpus callosum diffuse astrocytoma (WHO 2). He underwent a biopsy and was treated with radiotherapy and chemotherapy. The overall survival of the patient was 39 months (up to October 2020).
Feature Selection and Model Comparison
The extracted texture parameters needed to be normalized and standardized to reduce bias due to the missing values and outliers. The best non-redundant features were then selected from the extracted features using the variance method and mRMR. Finally, five features were extracted each from T2w, T2 FLAIR, and T1c images, and 12 features were extracted for the combined model (Table 3).
The three separate models and the combined model performed well on both the training and test cohorts, shown by the high AC. The performance of each model on the training and test cohorts, respectively, were: T2w: AUC = 0.86 vs. 0.734, AC = 0.84 vs. 0.72; T2 FLAIR: AUC = 0.89 vs. 0.78, AC = 0.90 vs. 0.76; T1c: AUC = 0.84 vs. 0.74, AC = 0.76 vs. 0.72; combined model: AUC = 0.96 vs. 0.87, AC = 0.93 vs. 0.84 (Table 4). The combined model had the best performance. The ROC curves of the T2w, T2 FLAIR, T1c, and combined models are shown in Figure 4. DCA of the four models showed that the combined model offered the largest net benefit (Figure 5). The radiomics nomogram of the combined model applied to the training cohort showed that it directly predicted the probability of OS > 12 months in patients with midline gliomas (Figure 6).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/866274/fneur-13-866274-HTML/image_m/fneur-13-866274-g004.jpg)
Figure 4. The receiver operating characteristic curves of the four models for the training (A) and test (B) cohorts. The combined model (red line) applied to the training and test cohorts had the best performance.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/866274/fneur-13-866274-HTML/image_m/fneur-13-866274-g005.jpg)
Figure 5. Decision-curve analysis of the four models. The x-axis represents the threshold probability, and the y-axis represents the net benefit. The combined model (red line) had a higher net benefit and better application value for predicting the survival time of patients with midline gliomas.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/866274/fneur-13-866274-HTML/image_m/fneur-13-866274-g006.jpg)
Figure 6. Nomogram of the combined model. A vertical line was drawn from each predictor to the point to obtain the score of the predictor. The scores of each predictor are then added, and a vertical line is drawn from the total score to the survival probability axis. The total score corresponds to the probability of >12 months survival.
Discussion
In this study, for the reproducibility and replicability of feature extraction, we used an open-source Pyradiomics 3.0.1 package of Python to calculate the features. To minimize the bias stemming from the use of data acquired at a single center and parameters of a single machine, we collected MRI data from different suppliers and different imaging protocols to recognize the diversity of data. Finally, we obtained optimal radiomics characteristics of T2w, T2 FLAIR, and T1c images to construct four models to non-invasively predict the prognosis of patients with midline gliomas. We compared the four models in terms of ROC curve, AUC, AC, sensitivity, and specificity and found that all models performed well on both the training and test cohorts, showing high AC. Among these, the combined model had the best performance. Moreover, DCA showed that the combined model offered the highest net benefit and application value. Our results demonstrate the feasibility of evaluating prognoses of patients with midline gliomas using traditional MRI radiomics. Furthermore, the fusion of radiomics characteristics obtained from different imaging sequences may improve the performance of the predictive model.
Numerous studies have shown that a single radiomics predictive model performs worse than a fusion radiomics model. Fusion radiomics models perform better and include more useful features, demonstrating the multimodal images play an important role in the prognostic evaluation of gliomas (23, 24). In this study, we constructed a machine learning model based on multi-parameter MRI data to predict the OS of patients with midline gliomas. Our results also showed that the fusion model performed better than the single radiomics model. In our training cohort, the combined model included 12 radiomics features (three features from T1c data, four features from T2 FLAIR data, and five features from T2w data), which were mainly derived from second- and higher-order texture features. The higher-order texture features can capture deeper imaging heterogeneity and provide information on tissue microstructure and local tumor microenvironment, which can help predict patient survival time (25, 26). Higher-order radiomics characteristics have also been shown to be closely related to genetic changes in tumors (6, 27–31). Five features in our combined model come from T2w, indicating that T2w images may provide more valuable texture features than T2 FLAIR and T1c images for evaluating survival time, which may explain why most studies prefer T2w sequences as the first choice for separate or joint analyses (26). However, T2 FLAIR and T1c images can offer additional details to the model and are thus indispensable (6, 27, 28). Therefore, we chose these three sequences. In addition, we applied the 12 features of the combined model to the nomogram, which enabled a more intuitive prediction of survival time.
To the best of our knowledge, the H3 K27M mutation, an independent prognostic factor, is common in midline gliomas. Several studies on the prediction of H3 K27M mutation status using machine learning based on radiomics features extracted from MRI data have been conducted to date. Su et al. (32) retrospectively studied 100 patients with midline gliomas and built 10 models using the Tree-based Pipeline Optimization Tool-based automated machine learning algorithm with radiomics features extracted from T2 FLAIR images. The AC of the best pipeline ranged from 0.788 to 0.867 in the training cohort and from 0.60 to 0.84 in the testing cohort. Pan et al. (14) retrospectively evaluated 151 brainstem gliomas and built a prediction model that incorporated 36 MRI features and three clinical features using a random forest algorithm. The model achieved an AC of 84.44% in the test cohort. Moreover, they constructed a simplified model that achieved an AC of 75.55% in the test cohort. Kandemirli et al. (33) developed an extreme gradient boosting algorithm classifier based on machine learning to predict the H3 K27M mutation using conventional MRI sequences. The model had an AUC of 0.791 in the training set and 0.737 in the test set, and the AC of the model was 73% in the test set. Furthermore, Li et al. (34) explored visually accessible MRI features of 30 patients with DMGs with and without the H3 K27M mutation using principal component analysis based on radiomics and found that T2w sequences may be more valuable and that cystic formation may be a biomarker for diagnosing DMGs with an H3 K27M mutation.
In addition to predicting the mutation status of H3 K27M, several researchers have evaluated the OS of patients with midline glioma using texture analysis. Szychot et al. (35) retrospectively analyzed 32 children with DMGs using a T2w imaging and apparent diffusion coefficient map-based texture analysis to predict outcomes of patients with DMG and found that the best predictor was mean of positive pixels, which could divide patients into poor and good prognosis groups according to a median survival time of 7.5 and 17.5 months, respectively. In this article, we included a larger number of cases to predict the OS of midline glioma patients using machine learning based on multiparametric MRI radiomics features, and our ACs were similar to those of the abovementioned machine learning studies. Thus, our model may offer an alternative method to predict the OS of patients with midline gliomas.
The current study has several limitations. First, a relatively small sample size and lack of an independent external test dataset in this study, which we will take as the key work in the next step. Second, the use of a single center may have introduced some bias. The next step is to conduct a multicenter study to improve the performance of the predictive model using more data. Third, only one method was used to build the model. If multiple methods were used to build and verify the models, a better model could be developed with improved AC. Forth, the patients who underwent radiotherapy in the long-term group are more than those in the short-term group, and some RT patients were treated with chemotherapy and/or targeted therapy at the same time, so it cannot be simply considered that RT can better predict the survival. In the next step, more cases need to be collected for subgroup analysis to detect the impact of RT on predicting the survival of DMG. Fifth, At present, there are many methods for extracting image features, such as local binary pattern (LBP), local ternary pattern (LTP), histogram of oriented gradients (HOG), etc. We used only one conventional method for feature extraction in this study due to the issues with overfitting and related techniques, In the future, we will try to use more methods to extract image features for further research.
Our radiomics model based on traditional MRI sequence data, especially the fusion radiomics model, may offer a reliable and non-invasive method to predict the survival time of patients with midline gliomas. Although there are several limitations, our findings provide a reference for the preoperative prediction and individualized treatment of gliomas as well as for further radiomics research on midline gliomas.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.
Ethics Statement
The studies involving human participants were reviewed and approved by IRB of Guangdong Sanjiu Brain Hospital. Written informed consent from the participants' legal guardian/next of kin was not required to participate in this study in accordance with the national legislation and the institutional requirements. Written informed consent was obtained from the individual(s), and minor(s)' legal guardian/next of kin, for the publication of any potentially identifiable images or data included in this article.
Author Contributions
D-BD and QZ: conceptualization. D-BD: methodology, formal analysis, and writing—original draft preparation. Y-TL: software. W-SW and QZ: validation. J-FZ: investigation. PH: resources. S-NW, PH, and L-NC: data curation. W-SW: visualization. QZ: supervision and project administration. All authors: writing—review and editing, read, and agreed to the submission of the manuscript.
Conflict of Interest
Y-TL was employed by GE Healthcare.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fneur.2022.866274/full#supplementary-material
References
1. Louis DN, Perry A, Reifenberger G, von Deimling A, Figarella-Branger D, Cavenee WK, et al. The 2016 World Health Organization Classification of tumors of the central nervous system: a summary. Acta Neuropathol. (2016) 131:803–20. doi: 10.1007/s00401-016-1545-1
2. Louis DN, Perry A, Wesseling P, Brat DJ, Cree IA, Figarella-Branger D, et al. The 2021 WHO classification of tumors of the central nervous system: a summary. Neuro Oncol. (2021) 23:1231-51. doi: 10.1093/neuonc/noab106
3. Qiu T, Chanchotisatien A, Qin Z, Wu J, Du Z, Zhang X, et al. Imaging Characteristics of Adult H3 K27M-mutant Gliomas. J Neurosurg. (2019) 15:1–9. doi: 10.3171/2020.2.JNS191920a
4. Ochs K, Ott M, Bunse T, Sahm F, Bunse L, Deumelandt K, et al. K27M-mutant Histone-3 as a Novel Target for Glioma Immunotherapy. Oncoimmunology. (2017) 6:e1328340. doi: 10.1080/2162402X.2017.1328340
5. Liu Z, Feng S, Li J, Cao H, Huang J, Fan F, et al. The Epidemiological Characteristics and Prognostic Factors of Low-Grade Brainstem Glioma: a Real-World Study of Pediatric and Adult Patients. Front Oncol. (2020) 10:391. doi: 10.3389/fonc.2020.00391
6. Huang RY, Guenette JP. Non-invasive Diagnosis of H3 K27M Mutant Midline Glioma. Neuro Oncol. (2020) 22:309–10. doi: 10.1093/neuonc/noz240
7. Aboian MS, Solomon DA, Felton E, Mabray MC, Villanueva-Meyer JE, Mueller S, et al. Imaging Characteristics of Pediatric Diffuse Midline Gliomas with Histone H3 K27M Mutation. AJNR Am J Neuroradiol. (2017) 38:795–800. doi: 10.3174/ajnr.A5076
8. Aboian MS, Tong E, Solomon DA, Kline C, Gautam A, Vardapetyan A, et al. Diffusion Characteristics of Pediatric Diffuse Midline Gliomas with Histone H3-K27M Mutation Using Apparent Diffusion Coefficient Histogram Analysis. Am J Neuroradiol. (2019) 40:1804–10. doi: 10.3174/ajnr.A6302
9. Zhang X, Tian Q, Wang L, Liu Y, Li B, Liang Z, et al. Radiomics Strategy for Molecular Subtype Stratification of Lower-Grade Glioma: Detecting IDH and TP53 Mutations Based on Multimodal MRI. J Magn Reson Imaging. (2018) 48:916–26. doi: 10.1002/jmri.25960
10. Peng H, Huo J, Li B, Cui Y, Zhang H, Zhang L, et al. Predicting Isocitrate Dehydrogenase (IDH) Mutation Status in Gliomas Using Multiparameter MRI Radiomics Features. J Magn Reson Imaging. (2021) 53:1399–407. doi: 10.1002/jmri.27434
11. Zhang S, Sun H, Su X, Yang X, Wang W, Wan X, et al. Automated machine learning to predict the co-occurrence of isocitrate dehydrogenase mutations and O6-methylguanine-DNA methyltransferase promoter methylation in patients with gliomas. J Magn Reson Imaging. (2021) 54:197–205. doi: 10.1002/jmri.27498
12. Zhou H, Vallières M, Bai HX, Su C, Tang H, Oldridge D, et al. MRI Features Predict Survival and Molecular Markers in Diffuse Lower-grade Gliomas. Neuro Oncol. (2017) 19:862–70. doi: 10.1093/neuonc/now256
13. Kocak B, Durmaz ES, Ates E, Sel I, Turgut Gunes S, Kaya OK, et al. Radiogenomics of Lower-grade Gliomas: Machine Learning-based MRI Texture Analysis for Predicting 1p/19q Codeletion Status. Eur Radiol. (2020) 30:877–86. doi: 10.1007/s00330-019-06492-2
14. Pan CC, Liu J, Tang J, Chen X, Chen F, Wu YL, et al. A Machine Learning-based Prediction Model of H3K27M Mutations in Brainstem Gliomas Using Conventional MRI and Clinical Features. Radiother Oncol. (2019) 130:172–79. doi: 10.1016/j.radonc.2018.07.011
15. Prasanna P, Patel J, Partovi S, Madabhushi A, Tiwari P. Radiomic features from the peritumoral brain parenchyma on treatment-naïve multi-parametric MR imaging predict long versus short-term survival in glioblastoma multiforme: Preliminary findings. Eur Radiol. (2017) 27:4188–97. doi: 10.1007/s00330-016-4637-3
16. Senders JT, Staples P, Mehrtash A, Cote DJ, Taphoorn MJB, Reardon DA, et al. An Online Calculator for the Prediction of Survival in Glioblastoma Patients Using Classical Statistics and Machine Learning. Neurosurgery. (2020) 86:2:E184–E192. doi: 10.1093/neuros/nyz403
17. Li S, Liu J, Xiong Y, Han Y, Pang P, Luo P, et al. Application Values of 2D and 3D Radiomics Models Based on CT Plain Scan in Differentiating Benign from Malignant Ovarian Tumors. Biomed Res Int. (2022) 2022:5952296. doi: 10.1155/2022/5952296
18. Lubner MG, Stabo N, Lubner SJ, del Rio AM, Song C, Halberg RB, et al. CT textural analysis of hepatic metastatic colorectal cancer: pre-treatment tumor heterogeneity correlates with pathology and clinical outcomes. Abdom Imaging. (2015) 40:2331–7. doi: 10.1007/s00261-015-0438-4
19. Akkus Z, Ali I, Sedlár J, Agrawal JP, Parney IF, Giannini C, et al. Predicting Deletion of Chromosomal Arms 1p/19q in Low-Grade Gliomas from MR Images Using Machine Intelligence. J Digit Imaging. (2017) 30:469-476. doi: 10.1007/s10278-017-9984-3
20. Han Y, Xie Z, Zang Y, Zhang S, Gu D, Zhou M, et al. Non-invasive genotype prediction of chromosome 1p/19q co-deletion by development and validation of an MRI-based radiomics signature in lower-grade gliomas. J Neurooncol. (2018) 140:297–306. doi: 10.1007/s11060-018-2953-y
21. Kniep HC, Madesta F, Schneider T, Hanning U, Schönfeld MH, Schön G, et al. Radiomics of Brain MRI: Utility in Prediction of Metastatic Tumor Type. Radiology (2019) 290:2:479-87. doi: 10.1148/radiol.2018180946
22. van Griethuysen JJM, Fedorov A, Parmar C, Hosny A, Aucoin N, Narayan V, et al. Computational Radiomics System to Decode the Radiographic Phenotype. Cancer Res. (2017) 77:e104–e107. doi: 10.1158/0008-5472.CAN-17-0339
23. Jiang C, Kong Z, Liu S, Feng S, Zhang Y, Zhu R, et al. Fusion Radiomics Features from Conventional MRI Predict MGMT Promoter Methylation Status in Lower Grade Gliomas. Eur J Radiol. (2019) 121:108714. doi: 10.1016/j.ejrad.2019.108714
24. Li ZC, Bai H, Sun Q, Li Q, Liu L, Zou Y, et al. Multiregional Radiomics Features from Multiparametric MRI for Prediction of MGMT Methylation Status in Glioblastoma Multiforme: a Multicentre Study. Eur Radiol. (2018) 28:3640–50. doi: 10.1007/s00330-017-5302-1
25. Aerts HJ, Velazquez ER, Leijenaar RT, Parmar C, Grossmann P, Carvalho S, et al. Decoding Tumour phenotype by Noninvasive Imaging Using a Quantitative Radiomics approach. Nat Commun. (2014) 5:4006. doi: 10.1038/ncomms5006
26. Qian Z, Li Y, Sun Z, Fan X, Xu K, Wang K, et al. Radiogenomics of Lower-grade Gliomas: a Radiomic Signature as a Biological Surrogate for Survival Prediction. Aging. (2018) 10:2884–99. doi: 10.18632/aging.101594
27. Seow P, Wong JHD, Ahmad-Annuar A, Mahajan A, Abdullah NA, Ramli N. Quantitative Magnetic Resonance Imaging and Radiogenomic Biomarkers for Glioma Characterisation: a Systematic Review. Br J Radiol. (2018) 91:20170930. doi: 10.1259/bjr.20170930
28. Sasaki T, Kinoshita M, Fujita K, Fukai J, Hayashi N, Uematsu Y, 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
29. Tian H, Wu H, Wu G, Xu G. Noninvasive Prediction of TERT Promoter Mutations in High-Grade Glioma by Radiomics Analysis Based on Multiparameter MRI. Biomed Res Int. (2020) 2020:3872314. doi: 10.1155/2020/3872314
30. Patel SH, Poisson LM, Brat DJ, Zhou Y, Cooper L, Snuderl M, et al. T2-FLAIR Mismatch, an Imaging Biomarker for IDH and 1p/19q Status in Lower-grade Gliomas: A TCGA/TCIA Project. Clin Cancer Res. (2017) 23:6078-85. doi: 10.1158/1078-0432.CCR-17-0560
31. Li Y, Qian Z, Xu K, Wang K, Fan X, Li S, et al. MRI Features Predict p53 Status in Lower-grade Gliomas via a Machine-learning Approach. Neuroimage Clin. (2017) 17:306–11. doi: 10.1016/j.nicl.2017.10.030
32. Su X, Chen N, Sun H, Liu Y, Yang X, Wang W, et al. Automated Machine Learning Based on Radiomics Features Predicts H3 K27M Mutation in Midline Gliomas of the Brain. Neuro Oncol. (2020) 22:393-401. doi: 10.1093/neuonc/noz184
33. Kandemirli SG, Kocak B, Naganawa S, Ozturk K, Yip SSF, Chopra S, et al. Machine learning-based multiparametric magnetic resonance imaging radiomics for prediction of H3K27M mutation in midline gliomas. World Neurosurg. (2021) 151:e78–85. doi: 10.1016/j.wneu.2021.03.135
34. Li Q, Dong F, Jiang B, Zhang M. Exploring MRI Characteristics of Brain Diffuse Midline Gliomas With the H3 K27M Mutation Using Radiomics. Front Oncol. (2021) 11:646267. doi: 10.3389/fonc.2021.646267
Keywords: machine learning, radiomics, MRI, midline glioma, overall survival
Citation: Deng D-B, Liao Y-T, Zhou J-F, Cheng L-N, He P, Wu S-N, Wang W-S and Zhou Q (2022) Non-Invasive Prediction of Survival Time of Midline Glioma Patients Using Machine Learning on Multiparametric MRI Radiomics Features. Front. Neurol. 13:866274. doi: 10.3389/fneur.2022.866274
Received: 31 January 2022; Accepted: 31 March 2022;
Published: 02 May 2022.
Edited by:
Nico Sollmann, University of California, San Francisco, United StatesReviewed by:
Jun Shen, Sun Yat-sen University, ChinaTing Song, Third Affiliated Hospital of Guangzhou Medical University, China
Jayant Jagtap, Symbiosis Institute of Technology, India
Vera Catharina Keil, VU Medical Center, Netherlands
Xiangrong Yu, Zhuhai People's Hospital, China
Copyright © 2022 Deng, Liao, Zhou, Cheng, He, Wu, Wang and Zhou. 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: Wen-Sheng Wang, wws161616@sina.com; Quan Zhou, zhouquan3777@smu.edu.cn