Dataset Information

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

ABSTRACT:

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.

SUBMITTER: Liu P

PROVIDER: S-EPMC9259879 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:BackgroundHepatocellular carcinoma (HCC) was the sixth common malignancies characteristic with highly aggressive in the world. It was well established that tumor mutation burden (TMB) act as indicator of immunotherapeutic responsiveness in various tumors. However, the role of TMB in tumor immune microenvironment (TIME) is still obscure.MethodThe mutation data was analyzed by employing "maftools" package. Weighted gene co-expression network analysis (WGCNA) was implemented to determine candidate module and significant genes correlated with TMB value. Differential analysis was performed between different level of TMB subgroups employing R package "limma". Gene ontology (GO) enrichment analysis was implemented with "clusterProfiler", "enrichplot" and "ggplot2" packages. Then risk score signature was developed by systematical bioinformatics analyses. K-M survival curves and receiver operating characteristic (ROC) plot were further analyzed for prognostic validity. To depict comprehensive context of TIME, XCELL, TIMER, QUANTISEQ, MCPcounter, EPIC, CIBERSORT, and CIBERSORT-ABS algorithm were employed. Additionally, the potential role of risk score on immune checkpoint blockade (ICB) immunotherapy was further explored. The quantitative real-time polymerase chain reaction was performed to detect expression of HTRA3.ResultsTMB value was positively correlated with older age, male gender and early T status. A total of 75 intersection genes between TMB-related genes and differentially expressed genes (DEGs) were screened and enriched in extracellular matrix-relevant pathways. Risk score based on three hub genes significantly affected overall survival (OS) time, infiltration of immune cells, and ICB-related hub targets. The prognostic performance of risks score was validated in the external testing group. Risk-clinical nomogram was constructed for clinical application. HTRA3 was demonstrated to be a prognostic factor in HCC in further exploration. Finally, mutation of TP53 was correlated with risk score and do not interfere with risk score-based prognostic prediction.ConclusionCollectively, a comprehensive analysis of TMB might provide novel insights into mutation-driven mechanism of tumorigenesis further contribute to tailored immunotherapy and prognosis prediction of HCC.