Bioinformatics Identification of Regulatory Genes and Mechanism Related to Hypoxia-Induced PD-L1 Inhibitor Resistance in Hepatocellular Carcinoma

The combination of a PD-L1 inhibitor and an anti-angiogenic agent has become the new reference standard in the first-line treatment of non-excisable hepatocellular carcinoma (HCC) due to the survival advantage, but its objective response rate remains low at 36%. Evidence shows that PD-L1 inhibitor resistance is attributed to hypoxic tumor microenvironment. In this study, we performed bioinformatics analysis to identify genes and the underlying mechanisms that improve the efficacy of PD-L1 inhibition. Two public datasets of gene expression profiles, (1) HCC tumor versus adjacent normal tissue (N = 214) and (2) normoxia versus anoxia of HepG2 cells (N = 6), were collected from Gene Expression Omnibus (GEO) database. We identified HCC-signature and hypoxia-related genes, using differential expression analysis, and their 52 overlapping genes. Of these 52 genes, 14 PD-L1 regulator genes were further identified through the multiple regression analysis of TCGA-LIHC dataset (N = 371), and 10 hub genes were indicated in the protein-protein interaction (PPI) network. It was found that POLE2, GABARAPL1, PIK3R1, NDC80, and TPX2 play critical roles in the response and overall survival in cancer patients under PD-L1 inhibitor treatment. Our study provides new insights and potential biomarkers to enhance the immunotherapeutic role of PD-L1 inhibitors in HCC, which can help in exploring new therapeutic strategies.