Background: The existing metabolic gene signatures for predicting breast cancer outcomes only focus on gene expression data without considering clinical characteristics. Therefore, this study aimed to establish a predictive risk model combining metabolic enzyme genes and clinicopathological characteristics to predict the overall survival in patients with breast cancer. Methods: Transcriptomics and corresponding clinical data for patients with breast cancer were downloaded from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases. Differentially expressed metabolic genes between tumors and normal tissues were identified in the TCGA dataset (training dataset). A prognostic model was then built using univariate and multifactorial Cox proportional hazards regression analyses in the training dataset. The capability of the predictive model was then assessed using the receiver operating characteristic in both datasets. Pathway enrichment analysis and immune cell infiltration were performed using Kyoto Encyclopedia of Genes and Genomes (KEGG)/Gene Ontology (GO) enrichment and CIBERSORT algorithm, respectively. Results: In breast cancer and normal tissues, 212 metabolic enzyme genes were differentially expressed. The predictive model included four factors: age, stage, and expression of SLC35A2 and PLA2G10. Patients with breast cancer were classified into high- and low-risk groups based on the model; the high-risk group had a significantly poorer overall survival rate than the low-risk group. Furthermore, the two risk groups showed different activation of pathways and alterations in the properties of tumor microenvironment-infiltrating immune cells. Conclusion: We developed a powerful model to predict prognosis in patients with breast cancer by combining the gene expression of metabolic enzymes with clinicopathological characteristics.
Keywords: TCGA; breast cancer; metabolism-related gene; prognosis; survival analysis.
Copyright © 2022 Lu, Liu and Zhang.