Identification and validation of an individualized metabolic prognostic signature for predicting the biochemical recurrence of prostate cancer based on the immune microenvironment

Abstract

Background Prostate cancer (PCa) is the most prevalent genitourinary malignancy in men, with a significant proportion of patients developing biochemical recurrence (BCR) after treatment. The immune microenvironment and metabolic alterations have crucial implications for the tumorigenesis and progression of PCa. Therefore, identifying metabolic genes associated with the immune microenvironment holds promise for predicting BCR and improving PCa prognosis. Methods In this study, ssGSEA and hierarchical clustering analysis were first conducted to evaluate and group PCa samples, followed by the use of the ESTIMATE and CIBERSORT algorithms to characterize the immunophenotypes and tumor microenvironment. The differential metabolic genes (MTGs) between groups were utilized to develop a prognostic-related signature. The predictive performance of the signature was assessed by principal component analysis (PCA), receiver operating characteristic (ROC) curve analysis, survival analysis, and the TIDE algorithm. A miRNA-MTGs regulatory network and predictive nomogram were constructed. Moreover, the expression of prognostic MTGs in PCa was detected by RT‒qPCR. Results PCa samples from the TCGA cohort were separated into two groups: the immune-low group and immune-high group. Forty-eight differentially expressed MTGs between the groups were identified, including 37 up-regulated and 11 down-regulated MTGs. Subsequently, CEL, CYP3A4, and PDE6G were identified as the genes most strongly associated with the BCR of PCa patients and these genes were utilized to establish the MTGs-based prognostic signatures. PCA, ROC curves analysis, Kaplan–Meier survival analysis, and the nomogram all showed the good predictive ability of the signature regardless of clinical variables. Furthermore, the MTGs-based signature was indicated as a potential predictive biomarker for immunotherapy response. Nine miRNAs involved in the regulation of prognostic MTGs were determined. In addition to the CEL gene, the PDE6G and CYP3A4 genes were expressed at higher levels in PCa samples. Conclusions The MTGs-based signature represents a novel approach with promising potential for predicting BCR in PCa patients.