Deep dissection of stemness-related hierarchies in hepatocellular carcinoma

Abstract

Background

Increasing evidence suggests that hepatocellular carcinoma (HCC) stem cells (LCSCs) play an essential part in HCC recurrence, metastasis, and chemotherapy and radiotherapy resistance. Multiple studies have demonstrated that stemness-related genes facilitate the progression of tumors. However, the mechanism by which stemness-related genes contribute to HCC is not well understood. Here, we aim to construct a stemness-related score (SRscores) model for deeper analysis of stemness-related genes, assisting with the prognosis and individualized treatment of HCC patients.Further, we found that the gene LPCAT1 was highly expressed in tumor tissues by immunohistochemistry, and sphere-forming assay revealed that knockdown of LPCAT1 inhibited the sphere-forming ability of hepatocellular carcinoma cells.

Methods

We used the TCGA-LIHC dataset to screen stemness-related genes of HCC from the MSigDB database. Prognosis, tumor microenvironment, immunological checkpoints, tumor immune dysfunction, rejection, treatment sensitivity, and putative biological pathways were examined. Random forest created the SRscores model. The anti-PD-1/anti-CTLA4 immunotherapy, tumor mutational burden, medication sensitivity, and cancer stem cell index were compared between the high- and low-risk score groups. We also examined risk scores for different cell types using single-cell RNA sequencing data and correlated transcription factor activity in cancer stem cells with SRscores genes. Finally, we tested core marker expression and biological functions.

Results

Patients can be divided into two subtypes (Cluster1 and Cluster2) based on the TCGA-LIHC dataset's identification of 11 stemness-related genes. Additionally, a SRscores was developed based on subtypes. Cluster2 and the group with the lowest SRscores had superior survival and immunotherapy response than Cluster1 and the group with the highest SRscores. The group with a high SRscores was significantly more enriched in classical tumor pathways than the group with a low SRscores. Multiple transcription factors and SRscores genes are correlated. The core gene LPCAT1 is highly expressed in rat liver cancer tissues and promotes tumor cell sphere formation.

Conclusion

A SRscores model can be utilized to predict the prognosis of HCC patients as well as their response to immunotherapy.