Abstract
PURPOSE: The response to transarterial chemoembolization (TACE) varies among individuals with hepatocellular carcinoma (HCC). This study aimed to identify a biomarker for predicting TACE response in HCC patients and to investigate its correlations with the tumor microenvironment and pre-TACE radiomics features. PATIENTS AND METHODS: GSE104580 data were obtained from the Gene Expression Omnibus (GEO) database. Differentially expressed gene analysis and machine learning algorithms were used to identify genes for constructing the TACE failure signature (TFS). TFS scores were then calculated for HCC patients in The Cancer Genome Atlas (TCGA) cohort. After obtaining images from The Cancer Imaging Archive (TCIA), tumor labeling and radiomics feature extraction, the Rad-score model was generated. Correlation analysis was performed between the TFS score and the Rad-score. CIBERSORT, ssGSEA and TME analysis were performed to explore differences in the immune landscape among distinct risk groups. The immunotherapy response was compared between different groups. RESULTS: ADH1C, CXCL11, EMCN, SPARCL1 and LIN28B were selected and incorporated into the TFS, which demonstrated satisfactory performance in predicting TACE response. Patients in the high TFS score group had poorer overall survival (OS) than those in the low TFS score group. The Rad-score model was constructed using six radiomics features, and the Rad-score was significantly correlated with hub gene expression and the TFS score. The high-TFS group was also characterized by an immunosuppressive tumor microenvironment and exhibited unfavorable responses to immunotherapy with PD-1 and CTLA-4 checkpoint inhibitors. CONCLUSION: This study established a transcriptomic biomarker for predicting the efficacy of TACE that correlates with radiomics features on pretreatment imaging, tumor immune microenvironment characteristics, and the efficacy of immunotherapy and targeted therapy in HCC patients.