Abstract
BACKGROUND: The tumor microenvironment (TME) has recently been proven to play a crucial role in the development and prognosis of tumors. However, the current knowledge on the potential of the TME in prostate cancer (PCa) remains scarce. PURPOSE: This study aims to elucidate the value of TME-related genes for PCa prognosis by integrative bioinformatics analysis. MATERIALS AND METHODS: We downloaded the immune and stromal scores of PCa samples via the ESTIMATE and correlated these scores to clinicopathological characteristics and recurrence-free survival (RFS) of patients. Based on these scores, the TME-related differentially expressed genes were identified for functional enrichment analysis. Cox regression analyses were performed to identify prognostic genes and establish a predictive risk model. Moreover, gene set enrichment analysis (GSEA) was performed to evaluate the relationship between risk score and immune pathway. RESULTS: The stromal and immune scores were associated with clinicopathological characteristics and RFS in PCa patients. In total, 238 intersecting differentially expressed genes were identified. Functional enrichment analysis further revealed that these genes dramatically participated in the immune-related pathways. The immune-related risk model was built with C-type lectin domain containing 7A (CLEC7A) and collagen type XI alpha 1 chain (COL11A1) using Cox regression analyses. Kaplan-Meier survival analysis showed that the expression levels of CLEC7A and COL11A1 were significantly associated with the RFS. Further, the RFS time in high-risk group was significantly shorter than that in low-risk group. The areas under the curve for the risk model in predicting 3- and 5-year RFS rates were 0.694 and 0.731, respectively. GSEA suggested that immunosuppression existed in high-risk PCa patients. CONCLUSION: CLEC7A and COL11A1 were selected to build a predictive risk model, which may help clinicians to assess the prognosis of PCa patients and select appropriate targets for immunotherapy.