Integrated Multi-Omics Data Analysis Reveals Associations Between Glycosylation and Stemness in Hepatocellular Carcinoma

Background: Glycosylation plays an essential role in driving the progression and treatment resistance of hepatocellular carcinoma (HCC). However, its function in regulating the acquisition and maintenance of the cancer stemness-like phenotype in HCC remains largely unknown. There is also very little known about how CAD and other potential glycosylation regulators may influence stemness. This study explores the relationship between glycosylation and stemness in HCC.

Methods: Gene set variance analysis (GSVA) was used to assess the TCGA pan-cancer enrichment in glycosylation-related pathways. Univariate, LASSO, and multivariate COX regression were then used to identify prognostic genes in the TCGA-LIHC and construct a prognostic signature. HCC patients were classified into high- and low-risk subgroups based on the signature. The relationship between gene expression profiles and stemness was confirmed using bulk and single-cell RNA-sequencing data. The role of CAD and other genes in regulating the stemness of HCC was also validated by RT-qPCR, CCK-8, and colony formation assay. Copy number variation (CNV), immune infiltration, and clinical features were further analyzed in different subgroups and subsequent gene expression profiles. Sensitive drugs were also screened.

Results: In the pan-cancer analysis, HCC was shown to have specific glycosylation alterations. Five genes, CAD, SLC51B, LGALS3, B3GAT3, and MT3, identified from 572 glycosylation-related genes, were used to construct a gene signature and predict HCC patient survival in the TCGA cohort. The results demonstrated a significant positive correlation between patients in the high-risk group and both elevated gene expression and HCC dedifferentiation status. A significant reduction in the stemness-related markers, CD24, CD44, CD20, FOXM1, and EpCAM, was found after the knockdown of CAD and other genes in HepG2 and Huh7 cells. Frequent mutations increased CNVs, immune-suppressive responses, and poor prognosis were also associated with the high-risk profile. The ICGC-LIRI-JP cohort confirmed a similar relationship between glycosylation-related subtypes and stemness. Finally, 84 sensitive drugs were screened for abnormal glycosylation of HCC, and carfilzomib was most highly correlated with CAD.

Conclusions: Glycosylation-related molecular subtypes are associated with HCC stemness and disease prognosis. These results provide new directions for further research on the relationship between glycosylation and stemness phenotypes.