Identification and Validation of a Tumor Microenvironment-Related Gene Signature in Hepatocellular Carcinoma Prognosis

Hepatocellular carcinoma (HCC) is a typical inflammatory-related malignant tumor with complex immune tolerance microenvironment and poor prognosis. In this study, we aimed to construct a novel immune-related gene signature for the prognosis of HCC patients, exploring tumor microenvironment (TME) cell infiltration characterization and potential mechanisms.

In this study, we constructed and validated a robust three-gene signature with close immune correlation in HCC, which presented a reliable performance in the prediction of HCC patients' survival.

A total of 364 HCC samples with follow-up information in the TCGA-LIHC dataset were analyzed for the training of the prognostic signature. The Least Absolute Shrinkage and Selector Operation (LASSO) regression based on the IRGs was conducted to identify the prognostic genes and establish an immune risk signature. The immune cell infiltration in TME was estimated via the CIBERSORT method. Gene Set Variation Analysis (GSVA) was conducted to compare the biological pathways involved in the low-risk and high-risk groups. Furthermore, paraffin sections of HCC tissue microarrays containing 77 patients from Fudan University Shanghai Cancer Center were used for IHC staining. The clinical characteristics of the 77 HCC patients were collected and summarized for survival analysis validation via the Kaplan-Meier (KM) method.

Three-gene signature with close immune correlation (Risk score = EPO * 0.02838 + BIRC5 * 0.02477 + SPP1 * 0.0002044) was constructed eventually and proven to be an effective prognostic factor for HCC patients. The patients were divided into a high-risk and a low-risk group according to the optimal cutoff, and the survival analysis revealed that HCC samples with high-risk immuno-score had significantly poorer outcomes than the low-risk group (p < 0.0001). The results of CIBERSORT suggested that the immune cell activation was relatively higher in the low-risk group with better prognosis. Besides, GSVA analysis showed multiple signaling differences between the high- and low-risk group, indicating that the three-gene prognostic model can affect the prognosis of patients by affecting immune-related mechanisms. Tissue microarray (TMA) results further confirmed that the expression of three genes in HCC tissues was closely related to the prognosis of patients, respectively.