Identification of a Novel Prognostic Signature Based on N-Linked Glycosylation and Its Correlation with Immunotherapy Response in Hepatocellular Carcinoma

BACKGROUND: The complex tumor microenvironment of hepatocellular carcinoma (HCC) has led to a low response to immune checkpoints inhibitors (ICIs) and a poor prognosis. PD-L1, as one of the indications for ICIs, is rich in glycosylation modifications, which result in untimely ICIs. Our study constructed a prognostic model based on N-linked glycosylation related genes for predicting the prognosis and the response to ICIs. METHODS: The list of N-linked glycosylation related genes is from the AmiGO2 database. The patients in The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) cohorts were enrolled. The Cox regression was performed to develop a prognostic model and patients were divided into a low- and high-risk subgroups. The role of signature in HCC was well investigated by prognostic analysis, gene set enrichment analysis, and immune infiltration analysis. 21 recurrent HCC patients who received postoperative adjuvant ICIs were recruited to evaluate the relationship between immunotherapy response and the signature. In vitro studies were conducted to investigate the oncogenic effects of DDOST, STT3A and TMEM165 in HCC. RESULTS: 59 N-linked glycosylation related differentially expressed genes were screened from HCC and normal tissues in the TCGA cohort. The prognostic model was developed with DDOST, STT3A and TMEM165. The risk score could be an independent prognostic factor. Patients in the high-risk subgroup showed a worse prognosis than patients in the low-risk one. ssGSEA showed that patients in the low-risk subgroup tended to be in the immune-activated state, with higher levels of B cell and macrophage cell infiltrations and lower levels of regulatory T cell (Treg) infiltrations in both TCGC and GEO cohorts. Immunohistochemistry studies showed that DDOST, STT3A and TMEM165 are highly expressed in tumor tissues and patients with a high-risk score correlated with poor progression free survival and worse immunotherapeutic response. Furthermore, the proliferation of HCC cells was reduced after the knockdown of DDOST, as well as upon the knockdown of STT3A and TMEM165. CONCLUSION: In this study, we establish that the risk model based on N-linked glycosylation related genes could efficiently predict the prognosis and tumor microenvironment immune state of HCC patients, and the risk score could serve as a novel indicator of immunotherapy.