Abstract
Background
Studies have demonstrated that propionate metabolism-related genes (PMRGs) are associated with cancer progression. PMRGs are not known to be involved in Hepatocellular carcinoma (HCC).
Methods
In this study, The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases were accessed for HCC-related transcriptome data and clinical information. First, DE-PMRGs were derived by intersecting PMRGs and DEGs between HCC tissues and normal controls. The clusterProfiler R package was then used to enrich DE-PMRGs. In addition, biomarkers of HCC were identified, and a prognostic model was developed. Using functional analysis and tumor microenvironment analysis, new insights were obtained into HCC. The expression of biomarkers was validated using quantitative real-time polymerase chain reaction (qRT-PCR).
Results
132 DE-PMRGs were obtained by intersecting 3690 DEGs and 291 PMRGs. Steroid and organic acid metabolism were associated with these genes. For the construction of the risk model for HCC samples, five biomarkers were identified, including Acyl-CoA dehydrogenase short chain (ACADS), CYP19A1, formiminotransferase cyclodeaminase (FTCD), glucose-6-phosphate dehydrogenase (G6PD), and glutamic-oxaloacetic transaminase (GOT2). ACADS, FTCD, and GOT2 were positive factors, whereas CYP19A1 and G6PD were negative. HCC patients with AUC greater than 0.6 were predicted to survive 1/2/3/4/5 years, indicating decent efficiency of the model. The probability of 1/3/5-survival for HCC was also predicted by the nomogram using the risk score, pathologic T stage, and cancer status. Moreover, functional enrichment analysis revealed the high-risk genes were associated with invasion and epithelial-mesenchymal transition. Significantly, immune cell infiltration and immune checkpoint expression were linked to HCC development.
Conclusion
This study identified five biomarkers of propionate metabolism that can predict HCC prognosis. This finding may provide a deeper understanding of PMRG function in HCC.
Abbreviations
HCC, Hepatocellular carcinoma; BCLC, Barcelona Clinic Liver Cancer; PMRG, Propionate Metabolism Related Genes; DE-PMRGs, differentially expressed PMRGs; TMS, Tumor Microenvironment; TCGA, The Cancer Genome Atlas; GEO, Gene Expression Omnibus; LASSO, Least Absolute Shrinkage and Selection Operator; ROC, Receiver Operating Characteristic; DCA, Decision Curve Analysis; ROC, Receiver Operating Characteristic; GSEA, Gene Set Enrichment Analysis; EMT, Epithelial–Mesenchymal Transition; IPS, Immunophenoscore; TIDS, Tumor Immune Dysfunction and Exclusion; ESTIMATE, Estimation of STromal and Immune cells in MAlignant Tumor tissues using Expression data; qRT-PCR, quantitative Real-Time Polymerase Chain Reaction; ACADS, Acyl-CoA Dehydrogenase Short Chain; FTCD, Formimidoyltransferase Cyclodeaminase; GOT2, glutamic-oxaloacetic transaminase 2; CYP19A1, cytochrome P450 family 19 subfamily A member 1; G6PD, glucose-6-phosphate dehydrogenase.
Patient Consent Statement
We adhere to the 1964 《 Declaration of Helsinki 》 and subsequent amendments or similar ethical standards. And we confirm that all clinical data contained in this study have been anonymous or confidential.
Acknowledgments
We thank Home for Researchers editorial team (www.home-for-researchers.com) for language editing service. We want to thank the whole study team at 3 Hospitals (Huaihe Hospital of Henan University, Henan Cancer Hospital and General Hospital of CPLA) for their continuous support.
Disclosure
The authors of this manuscript declare no relationships with any companies whose products or services may be related to the article’s subject matter.