Machine Learning and Radiomics Analysis for Tumor Budding Prediction in Colorectal Liver Metastases Magnetic Resonance Imaging Assessment

Abstract

Purpose: We aimed to assess the efficacy of machine learning and radiomics analysis using magnetic resonance imaging (MRI) with a hepatospecific contrast agent, in a pre-surgical setting, to predict tumor budding in liver metastases. Methods: Patients with MRI in a pre-surgical setting were retrospectively enrolled. Manual segmentation was made by means 3D Slicer image computing, and 851 radiomics features were extracted as median values using the PyRadiomics Python package. Balancing was performed and inter- and intraclass correlation coefficients were calculated to assess the between observer and within observer reproducibility of all radiomics extracted features. A Wilcoxon–Mann–Whitney nonparametric test and receiver operating characteristics (ROC) analysis were carried out. Balancing and feature selection procedures were performed. Linear and non-logistic regression models (LRM and NLRM) and different machine learning-based classifiers including decision tree (DT), k-nearest neighbor (KNN) and support vector machine (SVM) were considered. Results: The internal training set included 49 patients and 119 liver metastases. The validation cohort consisted of a total of 28 single lesion patients. The best single predictor to classify tumor budding was original_glcm_Idn obtained in the T1-W VIBE sequence arterial phase with an accuracy of 84%; wavelet_LLH_firstorder_10Percentile was obtained in the T1-W VIBE sequence portal phase with an accuracy of 92%; wavelet_HHL_glcm_MaximumProbability was obtained in the T1-W VIBE sequence hepatobiliary excretion phase with an accuracy of 88%; and wavelet_LLH_glcm_Imc1 was obtained in T2-W SPACE sequences with an accuracy of 88%. Considering the linear regression analysis, a statistically significant increase in accuracy to 96% was obtained using a linear weighted combination of 13 radiomic features extracted from the T1-W VIBE sequence arterial phase. Moreover, the best classifier was a KNN trained with the 13 radiomic features extracted from the arterial phase of the T1-W VIBE sequence, obtaining an accuracy of 95% and an AUC of 0.96. The validation set reached an accuracy of 94%, a sensitivity of 86% and a specificity of 95%. Conclusions: Machine learning and radiomics analysis are promising tools in predicting tumor budding. Considering the linear regression analysis, there was a statistically significant increase in accuracy to 96% using a weighted linear combination of 13 radiomics features extracted from the arterial phase compared to a single radiomics feature.