Dataset Information

Multi-omics analysis reveals prognostic value of tumor mutation burden in hepatocellular carcinoma.

ABSTRACT:

Background

Hepatocellular 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.

Method

The 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.

Results

TMB 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.

Conclusion

Collectively, 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.

SUBMITTER: Xu Q

PROVIDER: S-EPMC8254981 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:IntroductionTumor mutation burden (TMB) is an emerging biomarker for immunotherapy of hepatocellular carcinoma (HCC), but its value for clinical application has not been fully revealed.Materials and methodsWe used the Wilcox test to identify the differentially expressed immune-related genes (DEIRGs) in groups with high and low TMB as well as screened the immune gene pairs related to prognosis using univariate Cox regression analysis. A LASSO Cox regression prognostic model was developed by combining The Cancer Genome Atlas Liver Hepatocellular Carcinoma (TCGA-LIHC) with the International Cancer Genome Consortium (ICGC) LIRI-JP cohort, and internal (TCGA, ICGC) and external (GSE14520) validation analyses were conducted on the predictive value of the model. We also explored the relationship between the prognostic model and tumor microenvironment via the ESTIMATE algorithm and performed clinical correlation analysis by the chi-square test, revealing its underlying molecular mechanism with the help of Gene Set Enrichment Analysis (GSEA).ResultsThe prognostic model consisting of 15 immune gene pairs showed high predictive value for short- and long-term survival of HCC in three independent cohorts. Based on univariate multivariate Cox regression analysis, the prognostic model could be used to independently predict the prognosis in each independent cohort. The immune score, stromal score, and estimated score values were lower in the high-risk group than in the low-risk group. As shown by the chi-square test, the prognostic model exhibited an obvious correlation with the tumor stage [American Joint Committee on Cancer tumor-node-metastasis (AJCC-TNM) (p < 0.001), Barcelona Clinic Liver Cancer (BCLC) (p = 0.003)], histopathological grade (p = 0.033), vascular invasion (p = 0.009), maximum tumor diameter (p = 0.013), and background of liver cirrhosis (p < 0.001). GSEA revealed that the high-risk group had an enrichment of many oncology features, including the cell cycle, mismatch repair, DNA replication, RNA degradation, etc.ConclusionOur research developed and validated a reliable prognostic model associated with TMB for HCC, which may help to further enrich the therapeutic targets of HCC.

Project description:BackgroundGlycosylation 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.MethodsGene 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.ResultsIn 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.ConclusionsGlycosylation-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.