The mutual interplay between epigenetic modifications and metabolic rewiring contributes to malignant features of prostate adenocarcinoma (PRAD). This study aimed to uncover the biological roles of deubiquitylase USP35 in PRAD and find effective epigenetic or metabolic targets. Bioinformatic tools or methods revealed that USP35 is upregulated in PRAD samples and correlates with inferior prognosis. The in vitro and in vivo assays suggested that USP35 could enhance malignant features of PRAD cells. Mechanistically, we found that USP35 could directly deubiquitinate and stabilize BRPF1 proteins. USP35 depends on accumulated BRPF1 proteins to accelerate cell growth, stem-like properties, and migration in vitro and in vivo. Interestingly, high BRPF1 could bind to promoter of SREBP2 and activate the SREBP2 transcriptional capacity. Therefore, USP35/BRPF1 aixs could promote expressions of mevalonate (MVA) metabolism signature in a SREBP2-dependent manner. USP35 depends on BRPF1 to maintain the activity of mevalonate metabolism in PRAD cells. Last of all, we observed that targeting BRPF1 or using MVA inhibitor (atorvastatin) are effective to suppress USP35high PRAD in vivo tumor growth. USP35 is an indicator of MVA metabolic signature in PRAD. Collectively, our study highlighted the USP35/BRPF1/SREBP2 axis in modulating MVA metabolism in PRAD, suggesting the significance of BRPF1 or MVA as the potential therapeutic targets for PRAD treatment.

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.