Abstract
The time and speed of biochemical recurrence (BCR) of prostate cancer (PCa) after radical prostatectomy (RP) is highly variable. Stratification methods based on TNM staging and Gleason score (GS) do not allow the identification of patients at risk of BCR following RP. Therefore, the aim of the present study was to identify molecular signatures that can predict BCR risk effectively and facilitate treatment-related decisions for patients with PCa. RNA sequencing data and corresponding clinical data were downloaded from The Cancer Genome Atlas (TCGA) and Oncomine databases. Bioinformatics analysis was performed to identify differentially expressed genes in patients with GS=6 and GS ≥7. Cox regression models were used to determine the PCa signature (PCasig) and a clinical nomogram for the prediction of BCR. The performance of nomograms was assessed using time-dependent receiver operating characteristic curves and the concordance index (C-index). A PCasig comprising 10 genes, including SEMG2, KCNJ16, TFAP2B, SYCE1, KCNU1, AFP, GUCY1B2, GRIA4, NXPH1 and SOX11, was significantly associated with BCR, which was identified in TCGA cohort [hazard ratio (HR), 5.18; 95% CI, 3.241-8.272; C-index, 0.777] and validated in the Oncomine cohort (HR, 2.78; 95% CI, 1.39-5.54; C-index, 0.66). The expression levels of SEMG2, KCNJ16 and TFAP2B were downregulated in patients with GS ≥7. The expression levels of SYCE1, KCNU1, AFP, GUCY1B2, GRIA4, NXPH1 and SOX11 were upregulated in patients with GS ≥7. The clinical nomogram was constructed based on the GS and pathologic T stage (HR, 4.15; 95% CI, 1.39-5.54; C-index, 0.713). The addition of the PCasig to the clinical nomogram significantly improved prognostic value (HR, 7.25; 95% CI, 4.54-11.56; C-index, 0.782) with an net reclassification improvement of 75.3% (95% CI, 46.8-104.6%). Furthermore, the endogenous expression of each gene in the PCasig was measured in five PCa cell lines and in normal prostate cells, and these genes exhibited different expression levels relative to one another. In conclusion, an PCasig was identified by mining TCGA and successfully validated in an Oncomine cohort. This PCasig was an independent prognostic factor with a greater prognostic value for all patients regardless of GS than traditional clinical variables, which can improve the performance of clinical nomograms in predicting BCR of patients with GS ≥7.