ORIGINAL RESEARCH article
Sec. Cancer Imaging and Image-directed Interventions
Volume 9 - 2019 | https://doi.org/10.3389/fonc.2019.01164
Radiomics-Based Machine Learning Technology Enables Better Differentiation Between Glioblastoma and Anaplastic Oligodendroglioma
- 1Department of Ophthalmology, West China Hospital, Sichuan University, Chengdu, China
- 2State Key Laboratory of Biotherapy and Cancer Center, West China Hospital, Sichuan University, Chengdu, China
- 3Department of Neurosurgery, West China Hospital, Sichuan University, Chengdu, China
- 4Department of Radiology, West China Second University Hospital, Sichuan University, Chengdu, China
- 5School of Computer Science, Nanjing University of Science and Technology, Nanjing, China
- 6Department of Biotherapy, Cancer Center, West China Hospital, Sichuan University, Chengdu, China
- 7State Key Laboratory of Biotherapy and Cancer Center, West China Hospital, Sichuan University, and Collaborative Innovation Center for Biotherapy, Chengdu, China
Purpose: The aim of this study was to test whether radiomics-based machine learning can enable the better differentiation between glioblastoma (GBM) and anaplastic oligodendroglioma (AO).
Methods: This retrospective study involved 126 patients histologically diagnosed as GBM (n = 76) or AO (n = 50) in our institution from January 2015 to December 2018. A total number of 40 three-dimensional texture features were extracted from contrast-enhanced T1-weighted images using LIFEx package. Six diagnostic models were established with selection methods and classifiers. The optimal radiomics features were separately selected into three datasets with three feature selection methods [distance correlation, least absolute shrinkage and selection operator (LASSO), and gradient boosting decision tree (GBDT)]. Then datasets were separately adopted into linear discriminant analysis (LDA) and support vector machine (SVM) classifiers. Specificity, sensitivity, accuracy, and area under curve (AUC) of each model were calculated to evaluate their diagnostic performances.
Results: The diagnostic performance of machine learning models was superior to human readers. Both classifiers showed promising ability in discrimination with AUC more than 0.900 when combined with suitable feature selection method. For LDA-based models, the AUC of models were 0.986, 0.994, and 0.970 in the testing group, respectively. For the SVM-based models, the AUC of models were 0.923, 0.817, and 0.500 in the testing group, respectively. The over-fitting model was GBDT + SVM, suggesting that this model was too volatile that unsuitable for classification.
Conclusion: This study indicates radiomics-based machine learning has the potential to be utilized in clinically discriminating GBM from AO.
High-grade gliomas, the most common malignant solidary brain tumors in adults, are traditionally classified into anaplastic oligodendroglioma (AO), anaplastic astrocytoma (AA), and Glioblastoma (GBM) (1, 2). According to the WHO classification, AO was ascribed in Grade III while GBM was ascribed in Grade IV based on their histology characteristics. The early diagnosis of GBM and AO is clinically challenging but necessary due to their different treatment choice as well as the therapeutic responsiveness and patient survival (3). As for GBM extended resection is recommended to increase patient survival, whereas for AO this strategy lacks solid evidence (4–6). The treatment after surgery is also different as well. For GBM, standardized therapy after surgery recommended by NCCN guidelines is standard brain radiation therapy (RT) + concurrent temozolomide (TMZ) followed by adjuvant chemotherapy (7). While for AO, it is recommended to use fractionated external beam RT together with neoadjuvant or adjuvant PCV (procarbazine, lomustine, and vincristine) regarding the specific condition of each patient (7, 8).
A glioma-specific blood biomarker for glioma has not been identified yet. Therefore, the radiology examination is critical for tumor detection and lesion localization. Brain magnetic resonance imaging (MRI) plays a key role in the preoperative diagnostic of gliomas with high image resolution on tumor tissue. However, in some cases, MRI may be unable to provide enough information for differentiation between GBM and AO. The MRI characteristics of two tumors are pretty similar when GBM is characterized by perilesional vasogenic edema and ring-like enhancement (9–11); while AO also shows peritumoral edema and heterogeneous enhancement (11, 12). In this regard, the urgency of new radiological method has been highlighted.
Given that texture analysis on images provides a more objective information beyond naked eye assessment, quantitative descriptions of tumor characteristics could be an option for clinical diagnosis (13–16). Moreover, with digital parameters, new technology, such as machine learning, can be introduced for further statistical analysis. Machine learning, a hotspot in the field of artificial intelligence, enables the extraction of meaningful patterns from massive datasets and thereby achieving precise predictions with the model built (17). Machine learning has demonstrated outstanding performance in previous research including segmentation of the tumor, classification of certain types of tumor, and prediction of survival or genotype (18–23). Although the differentiation between GBM and AO is of high clinical relevance, the machine learning approach has never been explored yet. In this study, we investigated the feasibility of radiomics-based machine learning to differentiate GBM and AO.
Materials and Methods
In this retrospective single-center research, we viewed medical records in neurosurgery department to initially search for patients histologically diagnosed with GBM or AO from January 2015 to December 2018. The medical records were reviewed by two researchers to enroll the potentially qualified patients and to collect relevant clinical information for our research. The inclusion criteria for patients were: (1) with pathological diagnosis of GBM or AO in intraoperative freezing biopsy, and (2) with available high-quality pre-treatment MR scan performed at our institution before surgical resection. Then the pre-surgical MRI images of patients were exported from radiological department though Picture Archiving and Communication Systems (PACS) with uniform standard.
For patients before 2016, we made correction on their pathological diagnoses based on the new World Health Organization 2016 classification of gliomas by a senior neuropathologist with working experience of 10 years (24). The new standards required the presence of both IDH-mt and 1p19q co-deletion for the diagnosis of AO, otherwise it could only be regarded as NOS (Not Otherwise Specified) (24). Therefore, we excluded patients based on new classification who were with incomplete gene reports or with absent presence of both gene expression.
Seventy-nine consecutive patients with GBM and 56 consecutive patients with AO fulfilled inclusion criteria in the initial selection. Three patients with GBM and six patients with AO were excluded in the following evaluation according to the exclusion criteria, which were: (1) presence of motion artifacts on MRI, (2) previous history of brain surgery or biopsy, (3) previous history of intracranial diseases, such as subarachnoid hemorrhage, cerebral infarction, etc. Based on this strategy, a study cohort was built consisting of 76 GBM patients (mean age: 46.5 years) and 50 AO patients (mean age: 47.1 years).
All procedures performed in studies involving human participants were in accordance with the ethical standards of the institutional and/or national research committee and with the 1964 Helsinki declaration and its later amendments or comparable ethical standards. The institutional review board approved this retrospective study. The written informed consent was obtained from participants enrolled in this study. The written informed consent was necessary before radiological examination (written informed consent for patients <16 years old was signed by parents or guardians) for each patient. The patients agreed to undertake examination when needed and were informed that the statistics (including MR image), which could be used for academic purpose in the future, would be stored in our institutional database. The Ethics Committee of Sichuan University and radiology department of our institution have approved for statistics export and utilization for this study.
MR Image Acquisition
The current study focused on the conventional MR sequences. The suitable sequence should be chosen first for two reasons, that the descriptions on features boundary were vague in some sequences, and that features from all sequences would bring too much burden on classifiers. After initial evaluation on images and consultation with senior radiologists, the contrast-enhanced MRI sequence was the only one used to perform texture analysis in this study.
The MR scans were performed in the radiology department of institution. The contrast-enhanced MRI sequences were obtained with a 3.0T Siemens Trio Scanner using a MPRAGE sequence with the following imaging parameters: TR/TE/TI = 1900/2.26/900 ms, Flip angle = 9°, slice thickness = 1 mm, axial FOV = 25.6 × 25.6 cm2 and data matrix = 256 × 256. Intravenous injection of gadopentetate dimeglumine (0.1 mmol/Kg) was taken as contrast agent for patients. Multi-directional data for contrast-enhanced MRI were collected during the interval time of 90-250s. Figure 1 shows two examples of contrast-enhanced MRI images.
Figure 1. Two examples of contrast-enhanced MRI images. (A–C) Patient with GBM in (A) parasagittal, (B) axial, and (C) coronal view. (D–F) Patient with AO in (D) parasagittal, (E) axial, and (F) coronal view. GBM, glioblastoma; AO, anaplastic oligodendroglioma.
Human Readers Assessment
To test whether machine learning could outperform human readers, the diagnostic performance of them was compared. A senior neurosurgeon and a senior radiologist independently made diagnosis based on contrast-enhanced T1-weighted images, which were presented randomly, regarding classification as GBM or AO. Both readers were blinded to patient information and pathology reports. Then, the accuracy, sensitivity, and specificity were calculated for further analysis.
Texture Feature Extraction
The texture features of tumor tissue were extracted by two researchers using the LIFEx package under the supervision of a senior radiologist (25). Disagreements between researchers were recorded and adjudicated by consulting senior radiologists and neurosurgeons. The volume of interest (VOI) was drawn on T1C images by contouring the outer margin of tumor tissues slice by slice. The peritumoral edema band and adjacent structure invasion were separated from the primary tumor with the difference in contrast enhancement. For the lesions with multiple (>2) enhancement foci, ROI was only performed on the biggest one for those with clear boundary, and on tumor-confirmed area for those with vague boundary. After the ROI delineation, texture features were calculated automatically with default setting.
A total of 40 three-dimensional (3D) texture features were calculated from two orders. In the first order, texture features were calculated from shape histogram-based matrix and histogram-based matrix. In the second or higher order, features were calculated from gray-level co-occurrence matrix (GLCM), gray-level zone length matrix (GLZLM), neighborhood gray-level dependence matrix (NGLDM), and gray-level run length matrix (GLRLM). To avoid the interference of the lower image matrix resolution, texture analysis performed only on the VOIs with more than 64 voxels by default setting. All original data about extracted features were shown in Supplementary Material 1.
To ensure the validity and reproducibility of the extraction, the procedure was performed twice, and the difference between two sets was examined with Manny-Whitney U-test. We adjusted the q < 0.01 as significant (before was p < 0.05) to avoid the interference of false-positive errors rising from a large number of texture features. The results suggested that none of the features were significantly different, implying that the results could be considered reliable and reproducible (shown in Supplementary Material 2).
Classification Model Establishment
The purpose of machine learning was to train the models to predict whether each tumor was a GBM or AO with radiomics parameters extracted from the tumor tissue image. However, feature selection was necessary to eliminate statistically insignificant features and to avoid overfitting, which contributes to decreased running time and increased accuracy of the resulting models (26–28). In this study, we employed three selection methods with different selection mechanisms: distance correlation as representative of filter models, least absolute shrinkage and selection operator (LASSO) and gradient boosting decision tree (GBDT) as representatives of embedded models. Then, three datasets were generated with three different selection methods, which were each classified separately. The list of features selected with three different methods are shown in Supplementary Material 3, and the explanation of the features are summarized in Supplementary Material 4.
The next step was to choose suitable classifiers. Since linear classifier and non-linear classifier represent the state-of-the-art in pattern recognition, we adopted linear discriminant analysis (LDA) and support vector machine (SVM) classification algorithms in the current study as representatives of two classifier types (29). This way, overall six diagnostic models were established based on three selection methods and two classifiers.
As for the algorithm deployment, the study cohort was randomly divided into two subsets as training group and validation group on a proportion of 4:1. When the training on classifiers finished, the validation group was fed to evaluate the diagnostic performance of the models. Sensitivity, specificity, accuracy, and area under receiver operating characteristic curve (AUC) were calculated for both the training and validation group. To appraise the robustness of the methods, the procedure was repeated for 100 cycles with different and independent case assignments. The schematic workflow from image processing to machine learning is shown in Figure 2.
Among 126 patients were enrolled in the current study, 76 patients were diagnosed with GBM, and 50 patients with AO. The sex ratio, mean age, and time between MR scan and pathological diagnosis were summarized in Table 1. As for the human reader assessment, the accuracy for the neurosurgeon was 63.49%, and for the radiologist was 66.77%. Based on the results, a strong tendency on misdiagnosis of AO could be observed.
Diagnostic Performance of Models
The classification models exhibited promising discriminative ability when combined with suitable selection methods. For LDA-based models, all three models presented feasible performance with the AUC in the validation groups of 0.986, 0.994, and 0.970, respectively. For the SVM-based models, the models showed feasible performance with the AUC in the training groups of 0.923, 0.817, and 0.500. Overfitting was observed in one SVM-based model (SVM + GBDT), suggesting this model was volatile in application. The value of average sensitivity, specificity, accuracy, and AUC of training group and testing group are summarized in Table 2.
Table 2. Results of the discriminative model in distinguishing GBM from AO in the training and validation group.
Figure 3 represents the two-dimensional projection of the LDA-based models, illustrating that the GBM and AO formed distinctive clusters in the space defined by discriminant functions 1 and 2 generated by LASSO + LDA. Figure 4 shows the examples of performance of LDA-based models in terms of the distribution of the canonical functions in the 100 independent training cycles in the MRI analysis. A clear negative-values shift of the LDA function can be observed for AO, and all positive-values shift for GBM. ROCs of all models are shown in Supplementary Material 5.
Figure 3. Relationships between the canonical discriminant functions for GBM, AO and the group centroids. The two distinctive clusters formed by GBM and AO suggest three LDA-based models have excellent discriminant ability for GBM and AO. GBM, glioblastoma; AO, anaplastic oligodendroglioma; LDA, linear discriminant analysis; LASSO, least absolute shrinkage and selection operator; GBDT, gradient boosting decision tree; SVM, support vector machine. (A) Canonical discriminant functions for Distance Correlation + LDA. (B) Canonical discriminant functions for LASSO + LDA. (C) Canonical discriminant functions for GBDT + LDA.
Figure 4. Examples of performance of the LDA-based models in terms of the distribution of the canonical functions determined for the GBM and AO for two of the 100 independent training cycles in the MRI analysis. Minimal overlap is observed, suggesting high differential ability of the models. (A) Distance Correlation + LDA, (B) LASSO + LDA, (C) GBDT + LDA. GBM, glioblastoma; AO, anaplastic oligodendroglioma; AUC, area under curve; LDA, linear discriminant analysis; LASSO, least absolute shrinkage and selection operator; GBDT, gradient boosting decision tree.
For patients with high-grade gliomas, accurate tumor classification is clinically important because of its close relation with treatment strategy as well as therapeutic responsiveness and prognosis (3). In this study, we applied radiomics-based machine learning to pre-surgically differentiate between GBM and AO. Six models based on three selection methods (distance correlation, LASSO, and GBDT) and two classifiers (LDA and SVM) were built and evaluated. Our results demonstrated that machine learning approaches can be utilized and are clearly superior to human reader diagnosis.
Previous studies have explored the possibility of using machine learning for classification of brain tumor types (18, 21, 30). In the setting of gliomas, several studies have proved the utility of machine learning to differentiate between high-grade and low-grade glioma with high accuracy (0.80 and 0.945) (31, 32). In the setting of differentiation among specific histological subtypes of gliomas, a computer-aided diagnosis system was proposed and evaluated in a previous study to distinguish GBM from lower-grade gliomas, with positive results (33). A multicenter investigation also confirmed the feasibility of using 3D texture analysis for pediatric glioma subtype classification (medulloblastoma, pilocytic astrocytoma, and ependymoma) with an overall accuracy of 0.87 (34). The current study investigated a subject that has never been explored before, that the feasibility of radiomics-based machine learning in discriminating GBM from AO. Diagnostic performance of six models was assessed in the current study built on three selection methods (distance correlation, LASSO, and GBDT) coupled with two classifiers (LDA and SVM). In general, both classifiers showed high diagnostic performance with AUC more than 0.900 when combined with a suitable selection method. Nevertheless, when comparing between two classifiers, LDA-based models had slightly better diagnostic performance than that of SVM-based models.
The diagnostical models were established based on two types of classifies which differ in computing mechanism considering the performance of a certain classifier may be various in the settings of different tumors. LDA is a representative of the linear classifier which uses a straight line (a vector) to separates two classes (GBM and AO in this case), while SVM, a representative of the non-linear classifier, uses so-called support vectors to define a polynomic hyperplane to separate classes (35). In the settings of differentiating GBM and AO, our results showed LDA-based models had slightly better diagnostic performance than that of SVM-based models. However, the difference between the models was too slight to select the superior one, specifically given that all models investigated seemed to perform quite comparably and variance in AUC might be partially attributed due to the small statistical group. Therefore, limited by the small study cohort and relatively complicated methods, our results could only be regarded as hypothesis generation for future larger studies.
The results also implied that feature selection methods have impacts on diagnostic performance, especially for SVM-based models. Current feature selection methods can be categorized into three types depending on their selection mechanism: (1) Filter models select features by ranking them based on certain general characteristics such as correlation to remove irrelevant features without using any machine-learning algorithms. (2) Wrapper Models utilize a specific classifier to evaluate the quality of selected features, and offer a simple and powerful way to address the problem of feature selection, regardless of the chosen learning machine. (3) Embedded models are similar to wrapper models but embeds feature selection with classifier construction. Such models have the advantages of wrapper models-they include the interaction with the classification model, while embedded models are far less computationally intensive than wrapper models (28). In this study, we employed three selection methods as representatives of different selection mechanisms: distance correlation as representative of filter models, LASSO, and GBDT as representatives of embedded models. There was a common set of features selected by all three selection methods or two of the methods, which suggested these features might be important for the classification. For other features, it is hard to tell what extent they influenced the algorithms, since the AUCs showed minimal difference. However, even with feature selection, overfitting was still observed in one model (GBDT + SVM). We are unable to provide the exact reason but hypothesis that this model might be overly complex to be used as a discriminative tool to differentiate between GBM and AO.
Besides the comparison between machine learning models, we also performed comparison between machine and human readers in this study. Two readers unaware of the information on the exact number of patients were asked to make diagnosis on GBM or AO based on MR images. The readers were chosen from neurosurgery department or radiology department to ensure the convincing and reliable conclusions. The results were unexpected, considering there were only two options, that the diagnostical accuracy on AO was even lower than 0.500. Specifically, AO was easily misdiagnosed as GBM in human readers' radiological assessments. The explanations from the readers were the same that GBM and AO usually represented similar patterns on MR images, and they prefer to choose GBM rather than AO in these cases due to the epidemiological reason that the incidence of GBM is much higher than AO. As we mentioned before, the accurate pre-surgical diagnosis for two types of tumor is clinical important given the differences in surgical strategy. Therefore, it is reasonable to draw the conclusion that the patients will benefit from better treatment with machine learning clinical assistances. Machine intelligence will urge the radiological practice to change dramatically. However, we should also realize that the current machine technology is far from replacing human readers, and a combination of radiologist and machine might be the best choice for the foreseeable future. Radiologists still lead the central role in diagnosis while machine only act as assistance. This combination virtually eliminates simple blunders, increases play level, and provides better insight into the decision process (36).
Our study has several limitations. Firstly, it was a retrospective single-center investigation, which may lead to a patient selection bias and limited sample size. However, at present stage, it is still unknowable how much data is required to establish a predictive model, which may be answered through empirical investigation. The number of patients enrolled in previous studies focusing on similar topic ranged from 25 to 534 (31–34, 37–39). Secondly, we did not perform subgroup analysis regarding the IDH mutation status of GBM patients. Recent studies reported machine-learning based MRI texture analysis could be used as a new method for prediction of IDH mutational status, which suggested that IDH mutational status might have bearing on texture features (37, 39, 40). Thirdly, we used conventional contrast-enhanced MRI images only and did not use other sequences or advanced imaging tools such as magnetic resonance spectroscopy (MRS). Contrast-enhanced MRI sequence was chosen in this study for the clear delineation of tumor boundaries. The combined use of other sequences or imaging tools may enable better diagnostic ability. Fourthly, models built in current study were not externally validated. Since medical centers use different MRI scanners, imaging parameters and contrast, radiomic features may change accordingly. Therefore, the efficacy of machine learning-based models in this study cannot be guaranteed for external datasets. Nevertheless, we used the open-source package to perform the image processing and texture analysis, which allows others to reproduce the texture analysis with other datasets.
In conclusion, radiomics-based machine learning enables differentiation between glioblastoma and anaplastic oligodendroglioma. Our data indicate that the performance of this approach is superior to a human reader. This method may be a valuable addition to routine clinical practice to improve GBM and AO differentiation. However, multicenter investigations including larger patient cohorts and analysis combined with other MRI sequences or imaging techniques are warranted so that this non-invasive approach can be introduced into routine clinical practice.
Data Availability Statement
We are pleased to share our data to any qualified researchers without undue reservation. Please contact corresponding author if there is anything they need.
The studies involving human participants were reviewed and approved by Ethics Committee of Sichuan University. Written informed consent to participate in this study was provided by the participants' legal guardian/next of kin. 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.
YF, CC, JX, and XM contributed conception and design of the study. CC and ZT enrolled eligible patients, obtained medical records, and MRI images of each patient. FZ and CC independently made diagnosis based on contrast-enhanced T1-weighted images. YF and ZT did texture analysis. JW established the models and performed other statistical analysis. YF wrote the first draft of the manuscript. CC wrote sections of the manuscript. All authors contributed to manuscript revision, read, and approved the submitted version.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2019.01164/full#supplementary-material
Supplementary Material 1. Original data of extracted features. LASSO, least absolute shrinkage and selection operator; GBDT, gradient boosting decision tree.
Supplementary Material 2. Results of Manny-Whitney U-test to examine the difference between two sets of data extraction.
Supplementary Material 3. Texture features selected with three different methods. Features selected by all three methods are labeled red. Features selected by two methods are labeled in yellow. LASSO, least absolute shrinkage and selection operator; GBDT, gradient boosting decision tree.
Supplementary Material 4. Explanation of each texture features.
Supplementary Material 5. ROCs of six models to differentiate between GBM and AO. GBM, glioblastoma; AO, anaplastic oligodendroglioma; LASSO, least absolute shrinkage and selection operator; GBDT, gradient boosting decision tree; LDA, linear discriminant analysis; SVM, support vector machine.
AO, anaplastic oligodendroglioma; AUC, area under curve; GBDT, gradient boosting decision tree; GBM, glioblastoma; LASSO, least absolute shrinkage and selection operator; LDA, linear discriminant analysis; MRI, magnetic resonance imaging; SVM, support vector machine; T1C image, contrast-enhanced T1-weighted image; VOI, volume of interest.
2. Ostrom QT, Gittleman H, Truitt G, Boscia A, Kruchko C, Barnholtz-Sloan JS. CBTRUS statistical report: primary brain and other central nervous system tumors diagnosed in the United States in 2011-2015. Neuro Oncol. (2018) 20 (Suppl. 4):iv1–86. doi: 10.1093/neuonc/noy131
4. Lacroix M, Abi-Said D, Fourney DR, Gokaslan ZL, Shi W, DeMonte F, et al. A multivariate analysis of 416 patients with glioblastoma multiforme: prognosis, extent of resection, and survival. J Neurosurg. (2001) 95:190–8. doi: 10.3171/jns.2001.95.2.0190
6. Simpson JR, Horton J, Scott C, Curran WJ, Rubin P, Fischbach J, et al. Influence of location and extent of surgical resection on survival of patients with glioblastoma multiforme: results of three consecutive Radiation Therapy Oncology Group (RTOG) clinical trials. Int J Radiat Oncol Biol Phys. (1993) 26:239–44.
7. Nabors LB, Portnow J, Ammirati M, Baehring J, Brem H, Butowski N, et al. NCCN guidelines insights: central nervous system cancers, version 1.2019. J Natl Compr Canc Netw. (2019) 15:1331–45. doi: 10.6004/jnccn.2017.0166
8. Thomas AA, Abrey LE, Terziev R, Raizer J, Martinez NL, Forsyth P, et al. Multicenter phase II study of temozolomide and myeloablative chemotherapy with autologous stem cell transplant for newly diagnosed anaplastic oligodendroglioma. Neuro Oncol. (2017) 19:1380–90. doi: 10.1093/neuonc/nox086
13. Hatt M, Tixier F, Pierce L, Kinahan PE, Le Rest CC, Visvikis D. Characterization of PET/CT images using texture analysis: the past, the present…any future? Eur J Nucl Med Mol Imaging. (2017) 44:151–65. doi: 10.1007/s00259-016-3427-0
15. Ahn SY, Park CM, Park SJ, Kim HJ, Song C, Lee SM, et al. Prognostic value of computed tomography texture features in non-small cell lung cancers treated with definitive concomitant chemoradiotherapy. Invest Radiol. (2015) 50:719–25. doi: 10.1097/RLI.0000000000000174
16. Pickles MD, Lowry M, Gibbs P. Pretreatment prognostic value of dynamic contrast-enhanced magnetic resonance imaging vascular, texture, shape, and size parameters compared with traditional survival indicators obtained from locally advanced breast cancer patients. Invest Radiol. (2016) 51:177–85. doi: 10.1097/RLI.0000000000000222
18. Blumenthal DT, Artzi M, Liberman G, Bokstein F, Aizenstein O, Ben Bashat D. Classification of high-grade glioma into tumor and nontumor components using support vector machine. AJNR Am J Neuroradiol. (2017) 38:908–14. doi: 10.3174/ajnr.A5127
19. Citak-Er F, Firat Z, Kovanlikaya I, Ture U, Ozturk-Isik E. Machine-learning in grading of gliomas based on multi-parametric magnetic resonance imaging at 3T. Comput Biol Med. (2018) 99:154–60. doi: 10.1016/j.compbiomed.2018.06.009
21. Kunimatsu A, Kunimatsu N, Yasaka K, Akai H, Kamiya K, Watadani T, et al. Machine learning-based texture analysis of contrast-enhanced MR imaging to differentiate between glioblastoma and primary central nervous system lymphoma. Magn Reson Med Sci. (2019) 18:44–52. doi: 10.2463/mrms.mp.2017-0178
22. Papp L, Pötsch N, Grahovac M, Schmidbauer V, Woehrer A, Preusser M, et al. Glioma survival prediction with combined analysis of in vivo (11)C-MET PET features, ex vivo features, and patient features by supervised machine learning. J Nucl Med. (2018) 59:892–9. doi: 10.2967/jnumed.117.202267
23. Ranjith G, Parvathy R, Vikas V, Chandrasekharan K, Nair S. Machine learning methods for the classification of gliomas: initial results using features extracted from MR spectroscopy. Neuroradiol J. (2015) 28:106–11. doi: 10.1177/1971400915576637
25. Nioche C, Orlhac F, Boughdad S, Reuzé S, Goya-Outi J, Robert C, et al. LIFEx: a freeware for radiomic feature calculation in multimodality imaging to accelerate advances in the characterization of tumor heterogeneity. Cancer Res. (2018) 78:4786–9. doi: 10.1158/0008-5472.CAN-18-0125
29. Dellacasa Bellingegni A, Gruppioni E, Colazzo G, Davalli A, Sacchetti R, Guglielmelli E, et al. NLR, MLP, SVM, and LDA: a comparative analysis on EMG data from people with trans-radial amputation. J Neuroeng Rehabil. (2017) 14:82. doi: 10.1186/s12984-017-0290-6
31. Inano R, Oishi N, Kunieda T, Arakawa Y, Yamao Y, Shibata S, et al. Voxel-based clustered imaging by multiparameter diffusion tensor images for glioma grading. Neuroimage Clin. (2014) 5:396–407. doi: 10.1016/j.nicl.2014.08.001
32. Zhang X, Yan LF, Hu YC, Li G, Yang Y, Han Y, et al. Optimizing a machine learning based glioma grading system using multi-parametric MRI histogram and texture features. Oncotarget. (2017) 8:47816–30. doi: 10.18632/oncotarget.18001
34. Fetit AE, Novak J, Rodriguez D, Auer DP, Clark CA, Grundy RG, et al. Radiomics in paediatric neuro-oncology: a multicentre study on MRI texture analysis. NMR Biomed. (2018) 31:e3781. doi: 10.1002/nbm.3781
35. Gerhardt N, Schwolow S, Rohn S, Pérez-Cacho PR, Galán-Soldevilla H, Arce L, et al. Quality assessment of olive oils based on temperature-ramped HS-GC-IMS and sensory evaluation: comparison of different processing approaches by LDA, kNN, and SVM. Food Chem. (2019) 278:720–8. doi: 10.1016/j.foodchem.2018.11.095
37. Bisdas S, Shen H, Thust S, Katsaros V, Stranjalis G, Boskos C, et al. Texture analysis- and support vector machine-assisted diffusional kurtosis imaging may allow in vivo gliomas grading and IDH-mutation status prediction: a preliminary study. Sci Rep. (2018) 8:6108. doi: 10.1038/s41598-018-24438-4
38. Hu LS, Ning S, Eschbacher JM, Gaw N, Dueck AC, Smith KA, et al. Multi-parametric MRI and texture analysis to visualize spatial histologic heterogeneity and tumor extent in glioblastoma. PLoS ONE. (2015) 10:e0141506. doi: 10.1371/journal.pone.0141506
39. Jakola AS, Zhang YH, Skjulsvik AJ, Solheim O, Bø HK, Berntsen EM, et al. Quantitative texture analysis in the prediction of IDH status in low-grade gliomas. Clin Neurol Neurosurg. (2018) 164:114–20. doi: 10.1016/j.clineuro.2017.12.007
Keywords: machine learning, magnetic resonance imaging, glioblastoma, anaplastic oligodendroglioma, texture analysis
Citation: Fan Y, Chen C, Zhao F, Tian Z, Wang J, Ma X and Xu J (2019) Radiomics-Based Machine Learning Technology Enables Better Differentiation Between Glioblastoma and Anaplastic Oligodendroglioma. Front. Oncol. 9:1164. doi: 10.3389/fonc.2019.01164
Received: 14 June 2019; Accepted: 17 October 2019;
Published: 05 November 2019.
Edited by:Yin Li, University of Wisconsin-Madison, United States
Reviewed by:Andre Bongers, University of New South Wales, Australia
Yanwei Miao, Dalian Medical University, China
Copyright © 2019 Fan, Chen, Zhao, Tian, Wang, Ma and Xu. 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.
†These authors have contributed equally to this work