MRI radiomics to monitor therapeutic outcome of sorafenib plus IHA transcatheter NK cell combination therapy in hepatocellular carcinoma

Background Hepatocellular carcinoma (HCC) is a common liver malignancy with limited treatment options. Previous studies expressed the potential synergy of sorafenib and NK cell immunotherapy as a promising approach against HCC. MRI is commonly used to assess response of HCC to therapy. However, traditional MRI-based metrics for treatment efficacy are inadequate for capturing complex changes in the tumor microenvironment, especially with immunotherapy. In this study, we investigated potent MRI radiomics analysis to non-invasively assess early responses to combined sorafenib and NK cell therapy in a HCC rat model, aiming to predict multiple treatment outcomes and optimize HCC treatment evaluations. Methods Sprague Dawley (SD) rats underwent tumor implantation with the N1-S1 cell line. Tumor progression and treatment efficacy were assessed using MRI following NK cell immunotherapy and sorafenib administration. Radiomics features were extracted, processed, and selected from both T1w and T2w MRI images. The quantitative models were developed to predict treatment outcomes and their performances were evaluated with area under the receiver operating characteristic (AUROC) curve. Additionally, multivariable linear regression models were constructed to determine the correlation between MRI radiomics and histology, aiming for a noninvasive evaluation of tumor biomarkers. These models were evaluated using root-mean-squared-error (RMSE) and the Spearman correlation coefficient. Results A total of 743 radiomics features were extracted from T1w and T2w MRI data separately. Subsequently, a feature selection process was conducted to identify a subset of five features for modeling. For therapeutic prediction, four classification models were developed. Support vector machine (SVM) model, utilizing combined T1w + T2w MRI data, achieved 96% accuracy and an AUROC of 1.00 in differentiating the control and treatment groups. For multi-class treatment outcome prediction, Linear regression model attained 85% accuracy and an AUC of 0.93. Histological analysis showed that combination therapy of NK cell and sorafenib had the lowest tumor cell viability and the highest NK cell activity. Correlation analyses between MRI features and histological biomarkers indicated robust relationships (r = 0.94). Conclusions Our study underscored the significant potential of texture-based MRI imaging features in the early assessment of multiple HCC treatment outcomes.