Multiparametric MR radiomics in brain glioma: models comparation to predict biomarker status

Abstract

Abstract Background Genotype status of glioma have important significance to clinical treatment and prognosis. At present, there are few studies on the prediction of multiple genotype status in glioma by method of multi-sequence radiomics. The purpose of the study is to compare the performance of clinical features (age, sex, WHO grade, MRI morphological features etc.), radiomics features from multi MR sequence (T2WI, T1WI, DWI, ADC, CE-MRI (contrast enhancement)), and a combined multiple features model in predicting biomarker status (IDH, MGMT, TERT, 1p/19q of glioma. Methods In this retrospective analysis, 81 glioma patients confirmed by histology were enrolled in this study. Five MRI sequences were used for radiomic feature extraction. Finally, 107 features were extracted from each sequence on Pyradiomics software, separately. These included 18 first-order metrics, such as the mean, standard deviation, skewness, and kurtosis etc., 14 shape features and second-order metrics including 24 grey level run length matrix (GLCM), 16 grey level run length matrix (GLRLM), 16 grey level size zone matrix (GLSZM), 5 neighboring gray tone difference matrix (NGTDM), and 14 grey level dependence matrix (GLDM). Then, Univariate analysis and LASSO (Least absolute shrinkage and selection operator regression model were used to data dimension reduction, feature selection, and radiomics signature building. Significant features (p < 0.05 by multivariate logistic regression were retained to establish clinical model, T1WI model, T2WI model, T1 + C (T1WI contrast enhancement model, DWI model and ADC model, multi sequence model. Clinical features were combined with multi sequence model to establish a combined model. The predictive performance was validated by receiver operating characteristic curve (ROC analysis and decision curve analysis (DCA). Results The combined model showed the better performance in some groups of genotype status among some models (IDH AUC = 0.93, MGMT AUC = 0.88, TERT AUC = 0.76). Multi sequence model performed better than single sequence model in IDH, MGMT, TERT. There was no significant difference among the models in predicting 1p/19q status. Decision curve analysis showed combined model has higher clinical benefit than multi sequence model. Conclusion Multi sequence model is an effective method to identify the genotype status of cerebral glioma. Combined with clinical models can better distinguish genotype status of glioma. Key Points The combined model showed the higher performance compare with other models in predicting genotype status of IDH, MGMT, TERT. Multi sequence model showed a better predictive model than that of a single sequence model. Compared with other models, the combined model and multi sequence model show no advantage in prediction of 1p/19q status.