Machine learning-based radiomics analysis in predicting RAS mutational status using magnetic resonance imaging

Purpose

To assess the efficacy of radiomics features, obtained by magnetic resonance imaging (MRI) with hepatospecific contrast agent, in pre-surgical setting, to predict RAS mutational status in liver metastases.

Methods

Patients with MRI in pre-surgical setting were enrolled in a retrospective study. Manual segmentation was made by means 3D Slicer image computing, and 851 radiomics features were extracted as median values using the PyRadiomics Python package. The features were extracted considering the agreement with the Imaging Biomarker Standardization Initiative (IBSI). Balancing was performed through synthesis of samples for the underrepresented classes using the self-adaptive synthetic oversampling (SASYNO) approach. Inter- and intraclass correlation coefficients (ICC) were calculated to assess the between-observer and within-observer reproducibility of all radiomics characteristics. For continuous variables, nonparametric Wilcoxon–Mann–Whitney test was utilized. Benjamini and Hochberg's false discovery rate (FDR) adjustment for multiple testing was used. Receiver operating characteristics (ROC) analysis with the calculation of area under the ROC curve (AUC), sensitivity (SENS), specificity (SPEC), positive predictive value (PPV), negative predictive value (NPV) and accuracy (ACC) were assessed for each parameter. Linear and non-logistic regression model (LRM and NLRM) and different machine learning-based classifiers including decision tree (DT), k-nearest neighbor (KNN) and support vector machine (SVM) were considered. Moreover, features selection were performed before and after a normalized procedure using two different methods (3-sigma and z-score). McNemar test was used to assess differences statistically significant between dichotomic tables. All statistical procedures were done using MATLAB R2021b Statistics and Machine Toolbox (MathWorks, Natick, MA, USA).

Results

Seven normalized radiomics features, extracted from arterial phase, 11 normalized radiomics features, from portal phase, 12 normalized radiomics features from hepatobiliary phase and 12 normalized features from T2-W SPACE sequence were robust predictors of RAS mutational status. The multivariate analysis increased significantly the accuracy in RAS prediction when a LRM was used, combining 12 robust normalized features extracted by VIBE hepatobiliary phase reaching an accuracy of 99%, a sensitivity 97%, a specificity of 100%, a PPV of 100% and a NPV of 98%. No statistically significant increase was obtained, considering the tested classifiers DT, KNN and SVM, both without normalization and with normalization methods.

Conclusions

Normalized approach in MRI radiomics analysis allows to predict RAS mutational status.