Radiomics models based on multisequence MRI for predicting PD-1/PD-L1 expression in hepatocellular carcinoma

Abstract The purpose of this study was to explore the effectiveness of radiomics based on multisequence MRI in predicting the expression of PD-1/PD-L1 in hepatocellular carcinoma (HCC). One hundred and eight patients with HCC who underwent contrast-enhanced MRI 2 weeks before surgical resection were enrolled in this retrospective study. Corresponding paraffin sections were collected for immunohistochemistry to detect the expression of PD-1 and PD-L1. All patients were randomly divided into a training cohort and a validation cohort at a ratio of 7:3. Univariate and multivariate analyses were used to select potential clinical characteristics related to PD-1 and PD-L1 expression. Radiomics features were extracted from the axial fat-suppression T2-weighted imaging (FS-T2WI) images and the arterial phase and portal venous phase images from the axial dynamic contrast-enhanced MRI, and the corresponding feature sets were generated. The least absolute shrinkage and selection operator (LASSO) was used to select the optimal radiomics features for analysis. Logistic regression analysis was performed to construct single-sequence and multisequence radiomics and radiomic-clinical models. The predictive performance was judged by the area under the receiver operating characteristic curve (AUC) in the training and validation cohorts. In the whole cohort, PD-1 expression was positive in 43 patients, and PD-L1 expression was positive in 34 patients. The presence of satellite nodules served as an independent predictor of PD-L1 expression. The AUC values of the FS-T2WI, arterial phase, portal venous phase and multisequence models in predicting the expression of PD-1 were 0.696, 0.843, 0.863, and 0.946 in the training group and 0.669, 0.792, 0.800 and 0.815 in the validation group, respectively. The AUC values of the FS-T2WI, arterial phase, portal venous phase, multisequence and radiomic-clinical models in predicting PD-L1 expression were 0.731, 0.800, 0.800, 0.831 and 0.898 in the training group and 0.621, 0.743, 0.771, 0.810 and 0.779 in the validation group, respectively. The combined models showed better predictive performance. The results of this study suggest that a radiomics model based on multisequence MRI has the potential to predict the preoperative expression of PD-1 and PD-L1 in HCC, which could become an imaging biomarker for immune checkpoint inhibitor (ICI)-based treatment.