Single-Cell RNA Sequencing Reveals the Role of Phosphorylation-Related Genes in Hepatocellular Carcinoma Stem Cells

Abnormal activation of protein kinases and phosphatases is implicated in various tumorigenesis, including hepatocellular carcinoma (HCC). Advanced HCC patients are treated with systemic therapy, including tyrosine kinase inhibitors, which extend overall survival. Investigation of the underlying mechanism of protein kinase signaling will help to improve the efficacy of HCC therapy. Combining single-cell RNA sequencing data and TCGA RNA-seq data, we profiled the protein kinases, phosphatases, and other phosphorylation-related genes (PRGs) of HCC patients in this study. We found nine protein kinases and PRGs with high expression levels that were mainly detected in HCC cancer stem cells, including POLR2G, PPP2R1A, POLR2L, PRC1, ITBG1BP1, MARCKSL1, EZH2, DTYMK, and AURKA. Survival analysis with the TCGA dataset showed that these genes were associated with poor prognosis of HCC patients. Further correlation analysis showed that these genes were involved in cell cycle-related pathways that may contribute to the development of HCC. Among them, AURKA and EZH2 were identified as two hub genes by Ingenuity Pathway Analysis. Treatment with an AURKA inhibitor (alisertib) and an EZH2 inhibitor (gambogenic) inhibited HCC cell proliferation, migration, and invasion. We also found that both AURKA and EZH2 were highly expressed in TP53-mutant HCC samples. Our comprehensive analysis of PRGs contributes to illustrating the mechanisms underlying HCC progression and identifying potential therapeutic targets for future clinical trials.