The novel transcriptomic signature of angiogenesis predicts clinical outcome, tumor microenvironment and treatment response for prostate adenocarcinoma

Abstract

Angiogenesis plays the critical roles in promoting tumor progression, aggressiveness, and metastasis. Although few studies have revealed some angiogenesis-related genes (ARGs) could serve as prognosis-related biomarkers for the prostate cancer (PCa), the integrated role of ARGs has not been systematically studied. The RNA-sequencing data and clinical information of prostate adenocarcinoma (PRAD) were downloaded from The Cancer Genome Atlas (TCGA) as discovery dataset. Twenty-three ARGs in total were identified to be correlated with prognosis of PRAD by the univariate Cox regression analysis, and a 19-ARG signature was further developed with significant correlation with the disease-free survival (DFS) of PRAD by the least absolute shrinkage and selection operator (LASSO) Cox regression with tenfold cross-validation. The signature stratified PRAD patients into high- and low-ARGs signature score groups, and those with high ARGs signature score were associated with significantly poorer outcomes (median DFS: 62.71 months vs unreached, p < 0.0001). The predicting ability of ARGs signature was subsequently validated in two independent cohorts of GSE40272 & PRAD_MSKCC. Notably, the 19-ARG signature outperformed the typical clinical features or each involved ARG in predicting the DFS of PRAD. Furthermore, a prognostic nomogram was constructed with three independent prognostic factors, including the ARGs signature, T stage and Gleason score. The predicted results from the nomogram (C-index = 0.799, 95%CI = 0.744–0.854) matched well with the observed outcomes, which was verified by the calibration curves. The values of area under receiver operating characteristic curve (AUC) for DFS at 1-, 3-, 5-year for the nomogram were 0.82, 0.83, and 0.83, respectively, indicating the performance of nomogram model is of reasonably high accuracy and robustness. Moreover, functional enrichment analysis demonstrated the potential targets of E2F targets, G2M checkpoint pathways, and cell cycle pathways to suppress the PRAD progression. Of note, the high-risk PRAD patients were more sensitive to immune therapies, but Treg might hinder benefits from immunotherapies. Additionally, this established tool also could predict response to neoadjuvant androgen deprivation therapy (ADT) and some chemotherapy drugs, such as cisplatin, paclitaxel, and docetaxel, etc. The novel ARGs signature, with prognostic significance, can further promote the application of targeted therapies in different stratifications of PCa patients.