A novel stemness-hypoxia-related signature for prognostic stratification and immunotherapy response in hepatocellular carcinoma

Background The specific differentiation potential, unlimited proliferation, and self-renewal capacity of cancer stem cells (CSCs) are closely related to the occurrence, recurrence, and drug resistance of hepatocellular carcinoma (HCC), as well as hypoxia. Therefore, an in-depth analysis of the relationship between HCC stemness, oxygenation status, and the effectiveness of immunotherapy is necessary to improve the poor prognosis of HCC patients. Methods The weighted gene co-expression network analysis (WGCNA) was utilized to find hypoxia-related genes, and the stemness index (mRNAsi) was evaluated using the one-class logistic regression (OCLR) technique. Based on stemness-hypoxia-related genes (SHRGs), population subgroup categorization using NMF cluster analysis was carried out. The relationship between SHRGs and survival outcomes was determined using univariate Cox regression. The LASSO-Cox regression strategy was performed to investigate the quality and establish the classifier associated with prognosis. The main effect of risk scores on the tumor microenvironment (TME) and its response to immune checkpoint drugs was also examined. Finally, qRT-PCR was performed to explore the expression and prognostic value of the signature in clinical samples. Results After identifying tumor stemness- and hypoxia-related genes through a series of bioinformatics analyses, we constructed a prognostic stratification model based on these SHRGs, which can be effectively applied to the prognostic classification of HCC patients and the prediction of immune checkpoint inhibitors (ICIs) efficacy. Independent validation of the model in the ICGC cohort yielded good results. In addition, we also constructed hypoxic cell models in Herp3B and Huh7 cells to verify the expression of genes in the prognostic model and found that C7, CLEC1B, and CXCL6 were not only related to the tumor stemness but also related to hypoxia. Finally, we found that the constructed signature had a good prognostic value in the clinical sample. Conclusions We constructed and validated a stemness-hypoxia-related prognostic signature that can be used to predict the efficacy of ICIs therapy. We also verified that C7, CLEC1B, and CXCL6 are indeed associated with stemness and hypoxia through a hypoxic cell model, which may provide new ideas for individualized immunotherapy.