Abstract
Renal cell carcinoma (RCC) is the most common type of renal cancer, characterized by the dysregulation of metabolic pathways. RCC is the second highest cause of death among patients with urologic cancers and those with cancer cell metastases have a 5-years survival rate of only 10-15%. Thus, reliable prognostic biomarkers are essential tools to predict RCC patient outcomes. This study identified differentially expressed genes (DEGs) in the gene expression omnibus (GEO) database that are associated with pre-and post-metastases in clear cell renal cell carcinoma (ccRCC) patients and intersected these with metabolism-related genes in the Kyoto encyclopedia of genes and genomes (KEGG) database to identify metabolism-related DEGs (DEMGs). GOplot and ggplot packages for gene ontology (GO) and KEGG pathway enrichment analysis of DEMGs with log (foldchange) (logFC) were used to identify metabolic pathways associated with DEMG. Upregulated risk genes and downregulated protective genes among the DEMGs and seven independent metabolic genes, RRM2, MTHFD2, AGXT2, ALDH6A1, GLDC, HOGA1, and ETNK2, were found using univariate and multivariate Cox regression analysis, intersection, and Lasso-Cox regression analysis to establish a metabolic risk score signature (MRSS). Kaplan-Meier survival curve of Overall Survival (OS) showed that the low-risk group had a significantly better prognosis than the high-risk group in both the training cohort (p < 0.001; HR = 2.73, 95% CI = 1.97-3.79) and the validation cohort (p = 0.001; HR = 2.84, 95% CI = 1.50-5.38). The nomogram combined with multiple clinical information and MRSS was more effective at predicting patient outcomes than a single independent prognostic factor. The impact of metabolism on ccRCC was also assessed, and seven metabolism-related genes were established and validated as biomarkers to predict patient outcomes effectively.