Decoding Tumor–Immune Interactions in Hepatocellular Carcinoma Through Network-Centered Identification of CXCR2

Abstract

Hepatocellular carcinoma (HCC) is one of the most prevalent cancers worldwide and exhibits considerable biological heterogeneity in both molecular and clinical characteristics. The diverse molecular alterations and clinical manifestations of HCC indicate substantial heterogeneity across patient subgroups. This study aimed to identify novel therapeutic targets and predictive biomarkers associated with HCC using an integrative bioinformatics approach. High-throughput genomic datasets were obtained from the UCSC Xena browser to retrieve mRNA HTSeq-count data from the TCGA-HCC cohort. Gene co-expression network (GCN), protein–protein interaction network (PPIN), and enrichment analyses were performed to identify key dysregulated genes and their biological significance. Integrated network analyses identified three dysregulated hub genes, namely CXCR2, TLR2, and TLR4. Genomic alterations in these genes were further evaluated across tumor samples in the TCGA-HCC cohort. Kaplan–Meier (KM) survival analysis demonstrated that lower CXCR2 mRNA expression was significantly associated with poorer overall survival (OS) and recurrence-free survival (RFS). Furthermore, TIMER and UALCAN analyses revealed significant associations between CXCR2 expression and tumor purity, as well as immune cell infiltration levels, including T cells, macrophages, dendritic cells (DCs), and neutrophils. These findings suggest that CXCR2 is significantly associated with the immune microenvironment of HCC and represents a potential prognostic biomarker whose biological role warrants further mechanistic investigation.