Next Article in Journal
A Toolbox for Functional Analysis and the Systematic Identification of Diagnostic and Prognostic Gene Expression Signatures Combining Meta-Analysis and Machine Learning
Next Article in Special Issue
Characterization of HMGB1/2 Interactome in Prostate Cancer by Yeast Two Hybrid Approach: Potential Pathobiological Implications
Previous Article in Journal
The Microenvironment of Pituitary Tumors—Biological and Therapeutic Implications
Previous Article in Special Issue
Neuroendocrine Differentiation of Prostate Cancer—An Intriguing Example of Tumor Evolution at Play
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Cancers
Volume 11
Issue 10
10.3390/cancers11101603
cancers-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Novel Predictor Tool of Biochemical Recurrence after Radical Prostatectomy Based on a Five-MicroRNA Tissue Signature

by
Zhongwei Zhao
1,2,
Sabine Weickmann
1,
Monika Jung
1,
Michael Lein
3,4,
Ergin Kilic
5,
Carsten Stephan
1,4,
Andreas Erbersdobler
6,
Annika Fendler
7,† and
Klaus Jung
1,4,*,†
1
Department of Urology, Charité-Universitätsmedizin Berlin, 10117 Berlin, Germany
2
Department of Urology, Qilu Hospital of Shandong University, Jinan 250012, Shandong, China
3
Department of Urology, Sana Medical Center Offenbach, 63069 Offenbach am Main, Germany
4
Berlin Institute for Urologic Research, 10115 Berlin, Germany
5
Institute of Pathology, Hospital Leverkusen, 51375 Leverkusen, Germany
6
Institute of Pathology, University of Rostock, 18055 Rostock, Germany
7
Max Delbrueck Center for Molecular Medicine in the Helmholtz Association, Cancer Research Program, 13125 Berlin, Germany
*
Author to whom correspondence should be addressed.
†
These authors share senior authorship.
Cancers 2019, 11(10), 1603; https://0-doi-org.brum.beds.ac.uk/10.3390/cancers11101603
Submission received: 11 October 2019 / Accepted: 17 October 2019 / Published: 21 October 2019
(This article belongs to the Special Issue Prostate Cancer: Past, Present, and Future)
Browse Figures
Review Reports Versions Notes

Abstract

:
Within five to ten years after radical prostatectomy (RP), approximately 15–34% of prostate cancer (PCa) patients experience biochemical recurrence (BCR), which is defined as recurrence of serum levels of prostate-specific antigen >0.2 µg/L, indicating probable cancer recurrence. Models using clinicopathological variables for predicting this risk for patients lack accuracy. There is hope that new molecular biomarkers, like microRNAs (miRNAs), could be potential candidates to improve risk prediction. Therefore, we evaluated the BCR prognostic capability of 20 miRNAs, which were selected by a systematic literature review. MiRNA expressions were measured in formalin-fixed, paraffin-embedded (FFPE) tissue RP samples of 206 PCa patients by RT-qPCR. Univariate and multivariate Cox regression analyses were performed, to assess the independent prognostic potential of miRNAs. Internal validation was performed, using bootstrapping and the split-sample method. Five miRNAs (miR-30c-5p/31-5p/141-3p/148a-3p/miR-221-3p) were finally validated as independent prognostic biomarkers. Their prognostic ability and accuracy were evaluated using C-statistics of the obtained prognostic indices in the Cox regression, time-dependent receiver-operating characteristics, and decision curve analyses. Models of miRNAs, combined with relevant clinicopathological factors, were built. The five-miRNA-panel outperformed clinically established BCR scoring systems, while their combination significantly improved predictive power, based on clinicopathological factors alone. We conclude that this miRNA-based-predictor panel will be worth to be including in future studies.

1. Introduction

Prostate cancer (PCa) is the second most common cancer among men worldwide [1]. Six recent studies, including more than 1000 patients each, showed that approximately 15–34% of patients suffer from a biochemical recurrence (BCR), 5 to 10 years after radical prostatectomy (RP) [2,3,4,5,6,7]. The guidelines of the European Association of Urology (EAU) and the American Urological Association define BCR as a re-increase in prostate specific antigen (PSA) above 0.2 µg/L, which is confirmed by two consecutive elevated values [8,9]. This definition of BCR is associated with 90% probability of subsequently increasing PSA values [10]. At the moment, this PSA threshold is the most frequently used BCR definition in studies, and generally accepted as standard BCR value [11]. A current systematic review, based on fourteen studies which compared patients with and without BCR, confirmed BCR as decisive risk factor for distant metastasis, as well as PCa-specific and overall mortality [12]. The BCR cannot be equated with clinical recurrence after RP [9,12,13], but the clinician considers an increasing PSA level as the first sign of local recurrence or metastasis [13]. Without a secondary therapy following BCR, approximately 30% of patients experienced clinically manifested distant metastasis, and 19–27% of patients who experience BCR may suffer PCa-specific mortality within 10 years [10,14]. Early and more reliable prediction of PCa patients with high risk of BCR would aid the decision-making for adjuvant therapy, or more frequent monitoring during follow-up.
Clinicopathological variables, such as pathological tumor stage (pT stage), Gleason score and PSA, are important predictors, known to correlate with the survival of PCa patients. Currently, predictive scoring systems including these variables, are widely used for recurrence risk prediction. In 1998, D’Amico et al. [15] suggested a new risk classification method for BCR prediction in PCa patients. Subsequently, the Cancer of the Prostate Risk Assessment Postsurgical Score (CAPRAS) was used to predict BCR, based on pathological information from RP [16]. In 2005, Stephenson et al. proposed a postoperative nomogram including a number of clinicopathological variables to predict the 10-year probability of BCR after RP [17]. In addition, the National Comprehensive Cancer Network (NCCN) guidelines also suggested a novel risk classification system based on clinicopathological variables [18]. Even though clinicopathological factors certainly have a prognostic value, predictions based on them have a limited accuracy [19,20,21]. Despite the high value of the Gleason score for prognosis, there is considerable evidence that patients with the same Gleason score may experience vastly different clinical outcomes [22,23,24]. As one of the most powerful predictive scoring systems, CAPRAS was found to lack sufficient prediction accuracy in two meta-analyses [25,26]. Therefore, it is hoped that new tissue-based molecular biomarkers could be potential candidates to improve the disease recurrence prediction [27,28,29,30].
Efforts have been made to find genetic variables such as microRNAs (miRNA, miR) that may improve overall prediction accuracy [31,32,33]. We have recently reviewed all the related literature on miRNA-based BCR markers [34]. We found that more than 80% of previous studies focused only on single miRNA or multiple-miRNA-based models, while classical clinicopathological variables with known predictive value have seldom been taken into consideration during model building. However, miRNA-based predictors do not always outperform clinicopathological predictors [35]. A clinico-genomic model, combining well-defined prognostic models based on clinicopathological variables and miRNAs, should be a possibility of further improving accuracy [36,37]. Based on previous studies, it seems that models with multiple miRNAs (>2 miRNAs) perform better than a single miRNA predictor in PCa prognosis [33,38].
In this study, we aimed to (a) evaluate and compare miRNAs in prostatectomy specimens with the predictive potential of BCR reported in previous studies and filtered by a systematic literature search, (b) validate the predictive capability of these selected miRNAs by reverse-transcription quantitative real-time polymerase chain reaction (RT-qPCR), and using different evaluation approaches such as Cox regression analysis with internal validation, time-dependent receiver-operating characteristic curve analysis (ROC), and decision curve analysis, and (c) try to build a model, considering relevant clinicopathological variables combined with a multiple miRNA panel, based on the mentioned validation analyses.

2. Results

2.1. Patient Characteristics and Flow Chart of the Study

This retrospective study was performed on formalin-fixed, paraffin-embedded (FFPE) tissue specimens collected from 206 PCa patients, who underwent radical prostatectomy between 2001 and 2007 at the Charité–University Hospital. Patients were followed-up until October 2018. An alternate selection of available samples, based on patients with BCR, was conducted, so that the cohort comprised 98 (48%) patients with BCR and 108 (52%) patients without BCR. The flow diagram (Figure 1) outlines the steps performed in this study. The study was performed in accordance with the “Reporting Recommendations for Tumor Marker Prognostic Studies (REMARK)” [39].
Clinicopathological characteristics of the patient cohort, and the follow-up and recurrence free survival times, are summarized in Table 1. Statistically significant differences in most variables (age, PSA, pT status, ISUP Grade groups, and surgical margin; details in Table 1) were observed between recurrence-related PCa subgroups.

2.2. Selection of BCR-Related miRNAs and Their Differential Expression between Tumor and Normal Tissue and between the Samples from Patients with and without BCR

We decided to include the 20 miRNAs listed in Table 2 in our study. This selection was based on a comprehensive literature search, that we summarized in a systematic review on potential miRNA candidates of BCR prediction [34]. That review was based on 53 studies published between 2008 and 2017, and resulted in 58 distinct miRNAs as potential BCR markers. Sixteen out of these 58 miRNAs were recommended at least in two studies, or in a study with external validation, as well as partly supported by other study data when we carried out our analyses. The selection of four miRNAs (miR-31-5p, miR-34a-5p, miR-204-5p, miR-494-3p, and miR-939-5p), out of the 20 miRNAs listed in Table 2, was based on the subsequent literature data after the publication of the review [34], own data of previous microarray analyses, and/or their functional evidence regarding prostate cancer progression.
Expressions of the selected 20 miRNAs were then analyzed by RT-qPCR in 206 PCa specimens and 29 adjacent normal tissue samples. We identified 15 miRNAs, out of the 20 examined miRNAs, with differential expression in PCa patients, compared to adjacent normal controls (Figure S1; significances additionally indicated in Table 2). Of those, the expression of 10 miRNAs were found to be increased in PCa, whereas five miRNAs showed reduced levels. The expression data between the PCa specimen and adjacent normal tissue samples served only for controls and will not be further discussed. The comparison of miRNA expression in tumor tissue samples from patients with and without BCR is shown in Figure 2, and additionally included in Table 2. Significantly differential expression (p < 0.05) was observed for eight miRNAs, with levels always being downregulated in tissue samples from BCR patients, compared with expression levels in tissue samples from patients without BCR.

2.3. Correlations of Clinicopathological Variables with miRNAs and Correlations between miRNAs

Correlations between the various clinicopathological variables and the miRNAs were calculated by Spearman rank correlation analysis (Table 3).
PSA, surgical margin status, pN status, pT stage, and ISUP grade were significantly correlated with five, six, six, two, and nine different miRNAs, out of 20 examined miRNAs, respectively, while no correlations were found between age and DRE status (Table 3). For the dichotomous variables DRE, surgical margin, and pN status, and the ordinal scale variables pT stage and ISUP grade, additional calculations using the Mann–Whitney test and the Kruskal–Wallis test were compiled (Tables S2–S5). In addition, Figures S2 and S3 illustrate the behavior of the different miRNAs, depending on the pT and ISUP scale. However, all statistically significant correlations between clinicopathological variables and miRNAs are small and did not exceed the low correlation coefficient level of 0.30, with one exception (miR-1-3p and ISUP; rs = −0.351). In addition, significant associations were not observed between all clinicopathological factors and differentially expressed miRNAs with regard to BCR status (Fisher’s exact test, p-values between 0.161 and 1.000).
Correlations between the 20 analyzed miRNAs are summarized in Table S6. Approximately 25% of all Spearman rank correlation coefficients were higher than 0.50. However, all analyzed miRNAs were not clustered to other miRNAs forming the analyzed miRNA panel in this study, and were not related to other miRNAs of the miRNA gene family (Table S1).

2.4. Prognostic Potential of miRNAs Predicting BCR

Given the mostly weak or absent correlation of the differentially expressed miRNAs with clinicopathological variables, we hypothesized that our panel of miRNA might function as useful additional markers to predict BCR. According to particular recommendation of the REMARK guidelines [39], the continuous miRNA expression data and uncategorized data were used for all further miRNA analyses to avoid reduced power in detecting associations between miRNAs and BCR prediction. Univariate Cox regression analysis was performed to evaluate the predictive capability of all 20 examined miRNAs in tumor samples (Table 4). Those 16 miRNAs with p-values < 0.2 were selected for multivariate Cox regression analysis to avoid type II errors (Table 4). MiRNAs that were verified as independent factors (p < 0.05) in subsequent multivariate Cox regression analysis in a full model, and after stepwise backward elimination, are shown in Table 4. Six miRNAs (miR-30c-5p, miR-30d-5p, miR-31-5p, miR-141-3p, miR-148a-3p, and miR-221-3p) with p < 0.05, except for miR-30d-5p with p = 0.075, remained in the final model after backward elimination. In the following, this miRNA signature was named six-miR-panel.
Furthermore, we used the split-sample method as additional internal validation method splitting the patient cohort into a training (n = 140) and validation (n = 66) set. Multivariate Cox regression analysis, with stepwise backward elimination in the training set, resulted in a model with the same five miRNAs mentioned above, without miR-30d-5p, and was named five-miR-panel (Table S7). Using the Cox regression, we tested the predictive capacity of this five-miR-panel and the above-mentioned six-miR-panel in the total cohort, as well in the training and the test set (Table 5). C-statistics of the prognostic indices from both panels did not differ significantly, confirming that miR-30d-5p does not add a significant prediction advantage (Table 5). This was also confirmed in the decision curve analysis, showing similar curves for the five- and six-miR-panel (Figure S4). For practical reasons, we considered the five-miR-panel in further analyses.
Only two miRNAs (miR-31-5p and miR-221-3p) of this final five-miR-panel overlap with the eight miRNAs that showed significantly differential expression according to BCR status (p < 0.05, Table 2). Taking into account a p-value threshold of <0.100 with regard to the differential expression according to BCR status, four miRNAs (miR-31-5p, miR-221-3p, miR-30c-5p, and miR-141-3p) of this panel coincide with the respective 10 miRNAs (Table 2).

2.5. 5-miR-Panel Outperformed Models Based only on Clinicopathological Variables in Predicting BCR

Models based on clinicopathological variables are the conventional tools that are used for BCR in clinical practice to date. We therefore built models containing both clinicopathological information and miRNAs, to evaluate whether miRNAs might add a significant benefit to the prediction models based on clinicopathological factors.
In a first step, univariate Cox regression analysis of all clinicopathological variables was performed, resulting in six variables with p < 0.2, which were selected for further multivariate analyses (Table 6). The “full model” with all six clinicopathological factors and a “reduced model”, after stepwise backward elimination, finally including only pT stage and ISUP group grade as remaining independent predictors in the model, were established. C-statistics between the reduced and full model did not differ significantly (0.712 vs. 0.723, p = 0.230; Table 7).
C-statistics of the obtained prognostic indices were performed, to compare their discriminative abilities with those of other established predictive reference models mentioned in the Introduction, like CAPRAS [16], NCCN [18], according to D’Amico et al. [15] and Stephenson et al. [17] (Table 7).
As shown in Table 7, the CAPRAS model demonstrated the best AUC value among the four reference models, but its AUC value did not significantly differ from the values of our full and reduced models, (p = 0.187 and 0.618, respectively). We then added the five-miR-panel to all six clinicopathological models, to test whether this miRNA panel was able to improve BCR prediction significantly (Table 7). The discriminative ability of all six models based on clinicopathological variables could be significantly improved by the addition of the five-miR-panel (p-values of <0.0001 to 0.017).
To confirm the predictive benefit of all models, time-dependent ROC curve analyses were performed based on different postoperative time points (Figure 3). Depending on the median recurrence-free survival time of the total cohort after RP, the postoperative time points of 1, 2, 3, 4, and 5 years were applied. As shown in Figure 3a, the AUCs of the five-miR-panel and six-miR-panel were consistently higher than those of all six clinicopathological variable-based prediction models. The AUC lines of our full and reduced model, as well as the CAPRAS model, were quite similar, and consistently higher than those of the three remaining reference models at all time points. However, as shown in Figure 3b, inclusion of the five-miR-panel into the clinicopathological variable-based models improved all AUCs, which is consistent with the results in Table 7.
We also applied decision curve analysis as recommended most informative metrics [78] to demonstrate the incremental prognostic value of the five-miR-panel. It is exemplarily shown in the curves of Figure 4 for the reference models according to D’Amico et al. [15] and CAPRAS [16], as well as for our “full model”, in comparison to those of the five-miR-panel and the combined models. The curves indicate the benefit of the combined models in comparison to both single models, whereas C-statistics is not able to detect differences between the results of the combined models and the five-miR-panel (Table S8).

2.6. Functional Links Between miRNAs of the 5-miR-Panel and Prostate Cancer

Validated functional links between the miRNAs of the five-miR-panel and PCa are given in Table 8. Both tumor suppressive and oncogenic effects were observed for miR-141-3p, miR-148a-3p, and miR-221-3p. The differences of the five miRNAs, with regard to their expression levels (Table 2), hazard ratios (Table 4), and tumor suppressive/oncogenic functions (Table 8), show that associations between the molecular mechanisms and the outcome results are not evident. This is not surprising, because every miRNA can affect numerous target genes in different ways, whereas outcome results, here as predicted BCRs, are related to the integration of several data.

3. Discussion

In this retrospective study, we evaluated the potential of 20 individual miRNAs alone, or in combinations with clinicopathological variables, to predict BCR after RP. As shown in the present study, five miRNAs (miR-30c-5p, miR-31-5p, miR-141-3p, miR-148a-3p, and miR-221-3p) were finally validated as independent prognostic panel of BCR. Thus, the five-miR panel, together with standard clinicopathological factors, is a novel promising predictor of BCR, and will aid in decision-making and treatment management of PCa patients at risk of recurrence after RP.
As briefly outlined in the Introduction, the present study was designed to avoid the numerous shortcomings of previous studies, that we characterized in a systematic review of studies that used microRNAs as BCR predictors after RP [34]. For example, 44 out of the evaluated 53 studies (83%) measured expression levels of only one or two miRNAs; 18 studies (34%) calculated the predictive ability of miRNAs only in univariate analyses, without considering clinicopathological variables; and only nine studies (17%) showed an improved BCR prediction ability of miRNAs, in comparison to the prediction ability with clinicopathological variables; only 10 studies (19%) included a BCR group with more than 30 patients; and internal/external validation of the predictive models using miRNAs was reported in only eight studies (15%) [34]. Therefore, the present study was prepared with a total of 206 patients, including 98 patients suffering from BCR, and determination of 20 miRNAs. As these miRNAs were already proven in other studies as potential BCR predictors (Table 2), this approach provided a good opportunity to overcome the deficiencies of the previously reviewed miRNA-oriented BCR prediction studies. In this respect, two objectives were outlined in the present study, namely to (i) generate a robust multiple miRNA-based classifier and (ii) combine this classifier with relevant standard clinicopathological variables to a reliable model [36,91]. Nam et al. [38] described a panel of five miRNAs as a novel and promising predictor, after adjustment and combination with known clinicopathological variables. In addition, Kristensen et al. [33] showed that the prognostic benefit of a model with clinicopathological variables could be improved by the inclusion of three miRNAs. Both studies provide good examples for the establishment of clinico-genomic prediction models, based on miRNA classifier with multiple miRNAs.
By multivariate Cox regression analyses, we developed a six-miR-panel and five-miR-panel as effective BCR predictor classifiers. Two independent internal validation approaches, bootstrapping and the split-sample method, confirmed the robustness of these models, in line with the meaning of a “fit-for purpose” method [92]. C-statistics and decision curve analysis of both panels indicated similar performance for both panels (Table 5, Figure S4). We therefore decided to use the five-miR-panel as the final classifier for further analyses, in accordance with the principle “do not utilize more laboratory tests than absolutely necessary” [93].
The predictive potential of the individual miRNAs included in our final five-miR-panel was reported in previous studies (Table 2). For instance, downregulated expression of miR-30c-5p was an independent BCR marker in multivariate regression analysis including clinicopathological variables, [50] and reduced levels of miR-30c-5p were also found to be significantly associated with Gleason score and pathological stage [50], an association that was not confirmed in our study. In addition, downregulated miR-30c-5p expression, combined with expression of its target gene BCL9, was considered an independent predictor in multivariate regression analysis in another study [51]. Five previous studies focused on the role of miR-221-3p in predicting BCR [33,45,46,59,67]. Three of these studies confirmed reduced miR-221-3p levels to be a promising predictive biomarker, using multivariate Cox regression analysis, with adjustment of clinicopathological variables [33,45,67]. Furthermore, a strong correlation was found in these three studies between miR-221-3p expression, pathological stage, Gleason score, and pre-operative PSA levels, which was also not confirmed in our study. Two studies failed to detect the predictive potential of miR-221-3p in PCa patients in multivariate analyses, which may have resulted from the short follow-up time (median <2 years) as well as analytical reasons [46,59]. In two other studies [45,47], downregulated miR-141-3p was found to be associated with BCR, but validation as an independent predictor failed in multivariate regression analysis. Reduced levels of miR-148a-3p were observed to correlate with BCR in the study of Lichner et al. [47]. MiR-31-5p was already described as a deregulated miRNA in PCa [54,94] and reported in relationship to BCR [55], but was confirmed as a predictive biomarker of BCR in our multivariate prediction model for the first time.
We found important particularities between our study and the studies by Nam et al. [38] and Kristensen et al. [33]. Nam et al. [38] used five miRNAs in their signature, including miR-301a, which was most strongly associated with PCa recurrence. This miRNA was also included in our initial 20 miR measurement panel, but was not statistically significant when associated with BCR in univariate Cox regression analysis (p = 0.246, Table 4). The multiple miRNA signature by Kristensen et al. [33] consists of three miRNAs: miR-185-5p, miR-221-3p, and miR-326. All three miRNAs were also used in our initial panel and were significantly associated with BCR (Table 4), but only miR-221-3p remained as an independent marker in the final panel of multivariate Cox regression analysis after stepwise backward elimination (Table 4). Thus, it is worth noting that significant miRNAs in the other signatures failed in their inclusion in our final panel. After we had finished the experimental work in our study, three other panels for prediction of PCa outcome after RP based on multiple miRNA have been published. The working group of Sorensen et al. [95,96], with the above-mentioned panel of three miRNAs, introduced two novel panels. These are a four-miRNA ratio (miR-23a-3p x miR-10b-5p/miR-133a x miR-374b-5p), named MICaP [95], and a panel based on nine miRNAs, among them miR-185-5p and miR-221-3p, combined with three methylation markers, named miMe [96]. Nam et al. [97] introduced a new five-miRNA panel associated with metastasis after RP, while no miRNA was identical with those miRNAs recommended in their preceding report [38]. All these new panels showed promising predictive capacities, but the panels differed and were not compared within the same working group. Consistently good prediction accuracy data under conditions of differing miRNAs in the panels may be the result of different starting conditions, in particular the initial measurements of miRNAs, and the multivariate assessment of data. Ultimately, it is scarcely surprising that our five-miR-panel contains only two miRNAs with significantly different expression levels (p < 0.05) between the two BCR groups and three other miRNAs with p-values > 0.05. This panel is the result of multivariate data assessment that considers possible interactions between the various factors. On the other hand, it additionally confirms our selection approach of including all miRNAs in multivariate Cox regression analysis that had p-values < 0.200 in univariate Cox regressions. This procedure enables us to avoid type II errors, because all five miRNAs show characteristics of independent factors (Table 4).
The missing or weak correlations/associations of individual miRNAs of the final five-miR-panel, with conventional clinicopathological variables and among each other, characterize their potential as orthogonal predictors (Table 3, Table S5) [98]. Orthogonal biomarkers are characterized by the fact that they can improve the prediction accuracy of BCR, due to their additional information, independently from other variables. Associations of miRNAs with clinicopathological variables, for example, ISUP grading with miR-31-5p, miR-141-3p, and miR-148a-3p, did not exclude their potential as common independent predictive BCR biomarkers (Table 3, Table S5) [47].
Until now, the BCR prediction and other outcomes in clinical practice were generally based on prediction tools using standard clinicopathological variables (Table 7). In addition, a current systematic review of the prognostic impact of BCR on oncological outcome considered only clinical factors and tumor characteristics [12]. The limited prognostic accuracy of different approaches has been critically discussed [19,20,21,26]. Recently, several tissue-based tests known as Decipher, OncoType DX Prostate, Prolaris, and ProMark, using multigene or multiprotein expression data, have been introduced into clinical practice for PCa risk stratification and prognosis (reviewed in [99]). The NCCN Prostate Cancer Guidelines Panel stated in their 2019 guidelines that these tests might be considered for an initial risk stratification, based on biopsy samples and RP specimens [99]. Other current reports support this view [27,28,29,30,100,101,102]. Yet, the particularly different prognostic outcome results obtained by a head-to-head comparison of three tests (namely Decipher, Prolaris, and Oncotype), obviously complicate the routine implementation of these tests [103]. MiRNA-based panels, like the here presented five-miR-panel, have a great advantage in comparison to the mentioned multigene tests, because miRNAs are more robust molecular analytes than mRNAs [104]. This is particularly important under clinical routine conditions to procure stable tissue samples.
We would also classify the analysis of molecular biomarkers (specifically miRNA expression) in tissue specimens, either FFPE or fresh-frozen samples, from RP, as a promising tool to improve the prognostic accuracy of clinicopathological factors only. C-statistics and decision curve analysis clearly indicates that the five-miR-panel significantly improved the predictive accuracy of all clinicopathological models. Even the CAPRAS tool, which outperformed all other clinical models, was significantly improved. Our results indicate that there is considerable potential for improvement of the predictive ability and accuracy of the currently existing prognostic models (Table 7, Figure 3 and Figure 4). The mentioned studies by Nam et al. [38] and Kristensen et al. [33] support this view. We want to add that the sole focus on molecular markers to improve the prognostic accuracy contradicts the necessary global view in translational medicine to consider all aspects of disease associated features appropriately [105].
Besides identifying novel biomarkers, the identification of dysregulated miRNAs can also help improve understanding of the mechanisms of prostate cancer recurrence, which may provide a basis for novel therapeutic strategies or a better understanding of current treatment response. MiRNAs can function as both oncogenes and tumor suppressors [106]. The complexity of interpreting the functional role of the miRNAs that are relevant markers for BCR prediction can be illustrated by miR-141-3p and 148-3p, both included in the five-miR-panel. MiR-141-3p was characterized as a tumor suppressor in PCa, and inhibits PCa cell proliferation and migration by repressing genes involved in extracellular matrix-mediated pathways [47]. On the other hand, miR-141-3p was also described to act as an oncogene in untreated PCa and castration-resistant PCa, by enhancing PCa cell proliferation, even though no target genes have been described [107]. MiR-148-3p was shown to act as an oncogene by promoting prostate cancer cell growth due to repression of its target CAND1 [87], but it was also described as a tumor suppressor, by inhibiting the growth of androgen-refractory PCa cells through repression of mitogen- and stress-activated protein kinase 1 [88]. Regarding the different types of prostate cancer cells used in these studies, we assume that the expression of both miR-141 and miR-148a might be continuously altered during PCa progression. Both examples highlight that it might not be meaningful to characterize specific effects of individual miRNAs without considering possible interactions between the different molecular components. MiRNAs are part of a large network of competing endogenous RNAs, which further consists of mRNAs and long non-coding RNAs, that are regulated via miRNA responsive elements [108]. Changes in this network could explain the different functional outcome of an expressed miRNA in a specific setting, like in forthcoming BCRs. Moreover, the functional impact of a biomarker might explain the biological rationale of its effect, but does not reflect its clinical applicability and validity [109]. This view is supported by the fact that a concordance between the expression levels of the BCR and non-BCR samples, the hazard ratios, and the functional mechanisms of the miRNAs of the five-miR-panel, is not evident (Table 2, Table 4 and Table 8).
Despite our effort to avoid deficiencies observed in previous studies, as discussed above, some limitations of our study should be noted. These are the retrospective nature of the study, the sample collection from only a single center, and the lack of external validation. Another limitation is the focus on BCR as an endpoint, without considering clinical endpoints like metastasis-free survival or cancer-specific-free survival. Because of partially unclear or unavailable data, corresponding analyses were not possible, but were also not primarily intended. Future studies should focus on validating the predictive power of the five-miRNA panel in a prospective and multi-institutional setting. It would be worthwhile to include a head-to-head comparison with other multiple-miRNA-based tools, using a currently published novel EAU BCR risk stratification system as outcome endpoint [110].

4. Materials and Methods

4.1. Patient Selection and Data Collection

This retrospective study was approved by the local Ethics Committee of the Charité–University Hospital (EK-CCM-2004-09-14, approval date: September 20, 2004; EA1/153/07, approval date: October 22, 2007; EA1/134/12, approval date: June 22, 2012) and informed patient consent was obtained. The study was carried out in accordance with the Declaration of Helsinki. FFPE tissue specimens were collected from 206 PCa patients who underwent RP between 2001 and 2007 at Charité–University Hospital. Patients were treatment-naive before RP and were followed-up until October 2018. Follow-up data were based on medical records, telephone contacts with the patients’ urologists and patients or family members. Sample size was determined by a power-adapted calculation (α = 5%, power = 80%; Supplementary Information S1). The cohort comprised 98 (48%) patients with BCR, and 108 (52%) patients without BCR, in the follow-up after RP (Table 1). According to the guidelines of the European Association of Urology (EAU) and the American Urological Association (AUA), BCR was defined as a re-increase of PSA above 0.2 µg/L, which is confirmed by consecutive elevated values [8,9]. Twenty-nine samples from adjacent normal tissues were used as controls.
Histopathological evaluation, grading, staging, and determination of margin status are critical parameters for prognostic considerations in PCa. Pathological specimens were reviewed by two independent pathologists (E.K., A.E), in order to provide an unbiased and unanimous pathological diagnosis for all cases. Criteria of the International Union Against Cancer and the World Health Organization/International Society of Urological Pathology (ISUP) were applied.

4.2. RNA Extraction and Reverse Transcription Quantitative Real-Time Polymerase Chain Reaction (RT-qPCR)

All procedures were carried out as documented in detail in our previous publications [32,54,104,111,112], taking into account the items of “Minimum Information for Publication of Quantitative Real-time PCR Experiments” (MIQE guidelines) [113]. Briefly, total RNA was extracted from dissected FFPE tissue using the miRNeasy FFPE kit (Qiagen, Hilden, Germany), according to the manufacturer’s protocol. The RNA concentration and purity were measured on a NanoDrop ND-1000 spectrophotometer (NanoDrop Technologies, Wilmington, DE, USA). RNA extracts showed absorbance ratios of 260 nm to 280 nm and of 260 nm to 230 nm between 1.84 to 2.04 and 1.71 to 1.92, respectively, confirming their appropriate purity for further measurements. The selection of the 20 potential miRNA candidates for BCR prediction in this study has been explained in the above-mentioned Section 2.2. RT-qPCR analyses were carried out on the Light Cycler 480 Instrument (Roche Diagnostics, Mannheim, Germany) using TaqMan miRNA assays (Applied Biosystems, Foster City, CA, USA; Table S1). Expression values of all examined miRNAs were normalized to the geometric mean of let-7g-5p and miR-103a-3p, using the software qBasePLUS, v.3.0 software (Biogazelle, Zwijnaarde, Belgium). The suitability of both miRNAs as stable reference miRNAs in the normalization process is shown in Figure S5.

4.3. Statistical Analysis

IBM SPSS Statistics for Windows, version 25 (IBM Corp., Armonk, NY, USA) with bootstrap module, GraphPad Prism 8.20 (GraphPad Software, La Jolla, CA, USA), and MedCalc 19.06 (MedCalc Software, Ostend, Belgium) were used for statistical analysis. Nonparametric statistical tests used in our study included the Mann–Whitney U test, Chi-square or Fisher’s exact test, and Spearman rank correlation test. Univariate and multivariate Cox regression analyses were carried out for the prediction of BCR. For the internal validation of the models, bias-corrected and accelerated bootstrap calculations with the total cohort and the split-sample method (total cohort was randomly divided into a training (n = 140) and test (n = 66) set) were performed. Output data of Cox regression analyses defined as prognostic indices were utilized for C-statistics and decision curve analyses, as previously described [114,115]. AUCs under the time-dependent ROC curves of all models were additionally computed based on a cumulative case/dynamic control approach [77] at postoperative 1, 2, 3, 4, and 5 years, respectively. A p-value of < 0.05 (two-sided) was regarded as statistically significant.

5. Conclusions

In conclusion, we have identified miRNAs that are significantly associated with BCR after RP in prostate cancer, namely miR-30c-5p, miR-31-5p, miR-141-3p, miR-148a-3p and miR-221-3p. The inclusion of this five-miRNA-panel significantly improved the predictive accuracy of known clinicopathological models. Our study provides a good example of the benefit of incorporating molecular biomarkers, like miRNAs, into previously established risk prediction models based only on clinicopathological factors. Such an approach would be useful to identify patients at risk, set up personalized surveillance protocols, counsel patients, and select patients for adjuvant treatment trials. However, future prospective studies are needed, and comparison with other genomic classifiers, to implement a robust and reliable tool into the decision-making process of clinical practice.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2072-6694/11/10/1603/s1, Supplementary Information S1: Sample size and power calculations, Table S1: TaqMan MicroRNA assays used for RT-qPCR, Table S2: MicroRNAs in tumor samples in relationship to DRE, 3: MicroRNAs in tumor samples in relationship to the surgical margin status, Table S4: MicroRNAs in tumor samples in relationship to the lymph node status, Table S5: Associations between clinicopathological variables and microRNAs, Table S6: Spearman rank correlation coefficients among the microRNAs, Table S7: miRNA-based predictive BCR classifier using sample splitting approach, Table S8: C-statistics of three predictive BCR models, Figure S1: Expression of the 20 microRNAs, Figure S2: Expression of the 20 microRNAs analyzed in PCa tissue samples categorized according to the pT stages, Figure S3: Expression of the 20 microRNAs analyzed in PCa tissue samples categorized according to the ISUP grades, Figure S4: ROC curve and decision curve analysis of the 5-miR-panel and the 6-miR-panel, Figure S5: Stability of let-7g-5p and miR-103a-3p as suitability criterion for their use as normalizers.

Author Contributions

Conceptualization, C.S., A.E., A.F. and K.J.; Data curation, Z.Z., S.W. and K.J.; Formal analysis, Z.Z., C.S., A.F. and K.J.; Funding acquisition, M.L., A.E. and A.F.; Investigation, Z.Z., S.W., M.J., E.K. and A.E.; Methodology, Z.Z., S.W., M.J. and E.K.; Project administration, M.L., C.S., A.E. and K.J.; Resources, M.L., C.S., A.E. and A.F.; Supervision, C.S., A.E., A.F. and K.J.; Validation, Z.Z., S.W., M.J., E.K. and A.E.; Visualization, Z.Z., S.W., A.F. and K.J.; Writing—original draft, Z.Z., M.J. and K.J. Writing—review & editing, M.L., C.S., A.E., A.F. and K.J. All authors approved the final manuscript.

Funding

This research was partly funded by Foundation Urological Research (Grant-No. 17/02) and Wilhelm Sander-Stiftung (Grant No. 2010.111.1).

Acknowledgments

The authors thank Siegrun Blauhut for her valuable technical assistance.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Bray, F.; Ferlay, J.; Soerjomataram, I.; Siegel, R.L.; Torre, L.A.; Jemal, A. Global cancer statistics 2018: Globocan estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J. Clin. 2018, 68, 394–424. [Google Scholar] [CrossRef] [PubMed]
  2. Hu, X.; Cammann, H.; Meyer, H.-A.; Jung, K.; Lu, H.; Leva, N.; Magheli, A.; Stephan, C.; Busch, J. Risk prediction models for biochemical recurrence after radical prostatectomy using prostate-specific anitigen and Gleason score. Asian J. Androl. 2014, 16, 897–901. [Google Scholar] [PubMed]
  3. Tilki, D.; Mandel, P.; Schlomm, T.; Chun, F.K.; Tennstedt, P.; Pehrke, D.; Haese, A.; Huland, H.; Graefen, M.; Salomon, G. External validation of the CAPRA-S score to predict biochemical recurrence, metastasis and mortality after radical prostatectomy in a European cohort. J. Urol. 2015, 193, 1970–1975. [Google Scholar] [CrossRef] [PubMed]
  4. Liesenfeld, L.; Kron, M.; Gschwend, J.E.; Herkommer, K. Prognostic Factors for Biochemical Recurrence More than 10 Years after Radical Prostatectomy. J. Urol. 2017, 197, 143–148. [Google Scholar] [CrossRef]
  5. Ludwig, W.W.; Feng, Z.; Trock, B.J.; Humphreys, E.; Walsh, P.C. Prostate specific antigen testing after radical prostatectomy—Can we stop at 20 years? J. Urol. 2018, 199, 114–119. [Google Scholar] [CrossRef]
  6. Pak, S.; You, D.; Jeong, I.G.; Kim, Y.S.; Hong, J.H.; Kim, C.S.; Ahn, H. Time to biochemical relapse after radical prostatectomy and efficacy of salvage radiotherapy in patients with prostate cancer. Int. J. Clin. Oncol. 2019, 24, 1238–1246. [Google Scholar] [CrossRef]
  7. Loeb, S.; Zhu, X.; Schröder, F.H.; Roobol, M.J. Long-term radical prostatectomy outcomes among participants from the European Randomized Study of Screening for Prostate Cancer (ERSPC) Rotterdam. BJU Int. 2012, 110, 1678–1683. [Google Scholar] [CrossRef]
  8. Cornford, P.; Bellmunt, J.; Bolla, M.; Briers, E.; De Santis, M.; Gross, T.; Henry, A.M.; Joniau, S.; Lam, T.B.; Mason, M.D.; et al. EAU-ESTRO-SIOG Guidelines on Prostate Cancer. Part II: Treatment of relapsing, metastatic, and castration-resistant prostate cancer. Eur. Urol. 2017, 71, 630–642. [Google Scholar] [CrossRef]
  9. Cookson, M.S.; Aus, G.; Burnett, A.L.; Canby-Hagino, E.D.; D’Amico, A.V.; Dmochowski, R.R.; Eton, D.T.; Forman, J.D.; Goldenberg, S.L.; Hernandez, J.; et al. Variation in the definition of biochemical recurrence in patients treated for localized prostate cancer: The American Urological Association Prostate Guidelines for Localized Prostate Cancer Update Panel report and recommendations for a standard in the reporting of surgical outcomes. J. Urol. 2007, 177, 540–545. [Google Scholar]
  10. Stephenson, A.J.; Kattan, M.W.; Eastham, J.A.; Dotan, Z.A.; Bianco, F.J.; Lilja, H.; Scardino, P.T. Defining biochemical recurrence of prostate cancer after radical prostatectomy: A proposal for a standardized definition. J. Clin. Oncol. 2006, 24, 3973–3978. [Google Scholar] [CrossRef]
  11. Tourinho-Barbosa, R.; Srougi, V.; Nunes-Silva, I.; Baghdadi, M.; Rembeyo, G.; Eiffel, S.S.; Barret, E.; Rozet, F.; Galiano, M.; Cathelineau, X.; et al. Biochemical recurrence after radical prostatectomy: What does it mean? Int. Braz. J. Urol 2018, 44, 14–21. [Google Scholar] [CrossRef] [PubMed]
  12. Van den Broeck, T.; van den Bergh, R.C.N.; Arfi, N.; Gross, T.; Moris, L.; Briers, E.; Cumberbatch, M.; De, S.M.; Tilki, D.; Fanti, S.; et al. Prognostic value of biochemical recurrence following treatment with curative intent for prostate cancer: A systematic review. Eur. Urol. 2019, 75, 967–987. [Google Scholar] [CrossRef] [PubMed]
  13. Pound, C.R.; Partin, A.W.; Eisenberger, M.A.; Chan, D.W.; Pearson, J.D.; Walsh, P.C. Natural history of progression after PSA elevation following radical prostatectomy. JAMA 1999, 281, 1591–1597. [Google Scholar] [CrossRef] [PubMed]
  14. Brockman, J.A.; Alanee, S.; Vickers, A.J.; Scardino, P.T.; Wood, D.P.; Kibel, A.S.; Lin, D.W.; Bianco, F.J.; Rabah, D.M.; Klein, E.A.; et al. Nomogram predicting prostate cancer-specific mortality for men with biochemical recurrence after radical prostatectomy. Eur. Urol. 2015, 67, 1160–1167. [Google Scholar] [CrossRef] [PubMed]
  15. D’Amico, A.V.; Whittington, R.; Malkowicz, S.B.; Schultz, D.; Blank, K.; Broderick, G.A.; Tomaszewski, J.E.; Renshaw, A.A.; Kaplan, I.; Beard, C.J.; et al. Biochemical outcome after radical prostatectomy, external beam radiation therapy, or interstitial radiation therapy for clinically localized prostate cancer. JAMA 1998, 280, 969–974. [Google Scholar] [CrossRef]
  16. Cooperberg, M.R.; Pasta, D.J.; Elkin, E.P.; Litwin, M.S.; Latini, D.M.; Du Chane, J.; Carroll, P.R. The University of California, San Francisco Cancer of the Prostate Risk Assessment score: A straightforward and reliable preoperative predictor of disease recurrence after radical prostatectomy. J. Urol. 2005, 173, 1938–1942. [Google Scholar] [CrossRef]
  17. Stephenson, A.J.; Scardino, P.T.; Eastham, J.A.; Bianco, F.J.; Dotan, Z.A.; DiBlasio, C.J.; Reuther, A.; Klein, E.A.; Kattan, M.W. Postoperative nomogram predicting the 10-year probability of prostate cancer recurrence after radical prostatectomy. J. Clin. Oncol. 2005, 23, 7005–7012. [Google Scholar] [CrossRef]
  18. Mohler, J.L. The 2010 NCCN clinical practice guidelines in oncology on prostate cancer. J. Natl. Compr. Cancer Netw. 2010, 8, 145. [Google Scholar] [CrossRef]
  19. Remmers, S.; Verbeek, J.F.M.; Noeboer, D.; van der Kwast, T.; Roobol, M.J. Predicting biochemical recurrence and prostate cancer specific mortality after radical prostatectomy: Comparison of six prediction models in a cohort of screen and clinically detected prostate cancer patients. BJU Int. 2019, 124, 635–642. [Google Scholar] [CrossRef]
  20. Shariat, S.F.; Karakiewicz, P.I.; Suardi, N.; Kattan, M.W. Comparison of nomograms with other methods for predicting outcomes in prostate cancer: A critical analysis of the literature. Clin. Cancer Res. 2008, 14, 4400–4407. [Google Scholar] [CrossRef]
  21. Kang, M.; Jeong, C.W.; Choi, W.S.; Park, Y.H.; Cho, S.Y.; Lee, S.; Lee, S.B.; Ku, J.H.; Hong, S.K.; Byun, S.S.; et al. Pre- and post-operative nomograms to predict recurrence-free probability in korean men with clinically localized prostate cancer. PLoS ONE 2014, 9, e100053. [Google Scholar] [CrossRef] [PubMed]
  22. Goodman, M.; Ward, K.C.; Osunkoya, A.O.; Datta, M.W.; Luthringer, D.; Young, A.N.; Marks, K.; Cohen, V.; Kennedy, J.C.; Haber, M.J.; et al. Frequency and determinants of disagreement and error in gleason scores: A population-based study of prostate cancer. Prostate 2012, 72, 1389–1398. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Jung, J.W.; Lee, J.K.; Hong, S.K.; Byun, S.S.; Lee, S.E. Stratification of patients with intermediate-risk prostate cancer. BJU Int. 2015, 115, 907–912. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Siadat, F.; Sykes, J.; Zlotta, A.R.; Aldaoud, N.; Egawa, S.; Pushkar, D.; Kuk, C.; Bristow, R.G.; Montironi, R.; van der Kwast, T. Not all gleason pattern 4 prostate cancers are created equal: A study of latent prostatic carcinomas in a cystoprostatectomy and autopsy series. Prostate 2015, 75, 1277–1284. [Google Scholar] [CrossRef] [PubMed]
  25. Lorent, M.; Maalmi, H.; Tessier, P.; Supiot, S.; Dantan, E.; Foucher, Y. Meta-analysis of predictive models to assess the clinical validity and utility for patient-centered medical decision making: Application to the CAncer of the Prostate Risk Assessment (CAPRA). BMC Med. Inform. Decis. Mak. 2019, 19, 2. [Google Scholar] [CrossRef] [PubMed]
  26. Meurs, P.; Galvin, R.; Fanning, D.M.; Fahey, T. Prognostic value of the CAPRA clinical prediction rule: A systematic review and meta-analysis. BJU Int. 2013, 111, 427–436. [Google Scholar] [CrossRef] [PubMed]
  27. Hiser, W.M.; Sangiorgio, V.; Bollito, E.; Esnakula, A.; Feely, M.; Falzarano, S.M. Tissue-based multigene expression tests for pretreatment prostate cancer risk assessment: Current status and future perspectives. Future Oncol. 2018, 14, 3073–3083. [Google Scholar] [CrossRef]
  28. Spratt, D.E.; Zhang, J.; Santiago-Jimenez, M.; Dess, R.T.; Davis, J.W.; Den, R.B.; Dicker, A.P.; Kane, C.J.; Pollack, A.; Stoyanova, R.; et al. Development and validation of a novel integrated clinical-genomic risk group classification for localized prostate cancer. J. Clin. Oncol. 2018, 36, 581–590. [Google Scholar] [CrossRef]
  29. Ross, A.E.; Johnson, M.H.; Yousefi, K.; Davicioni, E.; Netto, G.J.; Marchionni, L.; Fedor, H.L.; Glavaris, S.; Choeurng, V.; Buerki, C.; et al. Tissue-based genomics augments post-prostatectomy risk stratification in a natural history cohort of intermediate- and high-risk men. Eur. Urol. 2016, 69, 157–165. [Google Scholar] [CrossRef]
  30. Fine, N.D.; LaPolla, F.; Epstein, M.; Loeb, S.; Dani, H. Genomic classifiers for treatment selection in newly diagnosed prostate cancer. BJU Int. 2019, 124, 574–586. [Google Scholar] [CrossRef]
  31. Bell, E.H.; Kirste, S.; Fleming, J.L.; Stegmaier, P.; Drendel, V.; Mo, X.; Ling, S.; Fabian, D.; Manring, I.; Jilg, C.A.; et al. A novel miRNA-based predictive model for biochemical failure following post-prostatectomy salvage radiation therapy. PLoS ONE 2015, 10, e0118745. [Google Scholar] [CrossRef] [PubMed]
  32. Fendler, A.; Jung, M.; Stephan, C.; Honey, R.J.; Stewart, R.J.; Pace, K.T.; Erbersdobler, A.; Samaan, S.; Jung, K.; Yousef, G.M. miRNAs can predict prostate cancer biochemical relapse and are involved in tumor progression. Int. J. Oncol. 2011, 39, 1183–1192. [Google Scholar] [PubMed]
  33. Kristensen, H.; Thomsen, A.R.; Haldrup, C.; Dyrskjot, L.; Hoyer, S.; Borre, M.; Mouritzen, P.; Orntoft, T.F.; Sorensen, K.D. Novel diagnostic and prognostic classifiers for prostate cancer identified by genome-wide microRNA profiling. Oncotarget 2016, 7, 30760–30771. [Google Scholar] [CrossRef] [PubMed]
  34. Zhao, Z.W.; Stephan, C.; Weickmann, S.; Jung, M.; Kristiansen, G.; Jung, K. Tissue-based microRNAs as predictors of biochemical recurrence after radical prostatectomy: What can we learn from past studies? Int. J. Mol. Sci. 2017, 18, 2023. [Google Scholar] [CrossRef]
  35. Karatas, O.F.; Guzel, E.; Suer, I.; Ekici, I.D.; Caskurlu, T.; Creighton, C.J.; Ittmann, M.; Ozen, M. miR-1 and miR-133b are differentially expressed in patients with recurrent prostate cancer. PLoS ONE 2014, 9, e98675. [Google Scholar] [CrossRef]
  36. Bovelstad, H.M.; Nygard, S.; Borgan, O. Survival prediction from clinico-genomic models—A comparative study. BMC Bioinform. 2009, 10, 413. [Google Scholar] [CrossRef]
  37. Cooperberg, M.R.; Davicioni, E.; Crisan, A.; Jenkins, R.B.; Ghadessi, M.; Karnes, R.J. Combined value of validated clinical and genomic risk stratification tools for predicting prostate cancer mortality in a high-risk prostatectomy cohort. Eur. Urol. 2015, 67, 326–333. [Google Scholar] [CrossRef]
  38. Nam, R.K.; Amemiya, Y.; Benatar, T.; Wallis, C.J.; Stojcic-Bendavid, J.; Bacopulos, S.; Sherman, C.; Sugar, L.; Naeim, M.; Yang, W.; et al. Five microRNA signature predictive of prostate cancer recurrence and metastasis: A cohort study. J. Cancer 2015, 6, 1160–1171. [Google Scholar] [CrossRef]
  39. Sauerbrei, W.; Taube, S.E.; McShane, L.M.; Cavenagh, M.M.; Altman, D.G. Reporting Recommendations for Tumor Marker Prognostic Studies (REMARK): An abridged explanation and elaboration. J. Natl. Cancer Inst. 2018, 110, 803–811. [Google Scholar] [CrossRef]
  40. Wei, W.; Leng, J.; Shao, H.; Wang, W. MiR-1, a potential predictive biomarker for recurrence in prostate cancer after radical prostatectomy. Am. J. Med. Sci. 2017, 353, 315–319. [Google Scholar] [CrossRef]
  41. Hudson, R.S.; Yi, M.; Esposito, D.; Watkins, S.K.; Hurwitz, A.A.; Yfantis, H.G.; Lee, D.H.; Borin, J.F.; Naslund, M.J.; Alexander, R.B.; et al. MicroRNA-1 is a candidate tumor suppressor and prognostic marker in human prostate cancer. Nucleic Acids Res. 2012, 40, 3689–3703. [Google Scholar] [CrossRef] [PubMed]
  42. Leite, K.R.; Reis, S.T.; Viana, N.; Morais, D.R.; Moura, C.M.; Silva, I.A.; Pontes, J.; Katz, B.; Srougi, M. Controlling RECK miR21 promotes tumor cell invasion and is related to biochemical recurrence in prostate cancer. J. Cancer 2015, 6, 292–301. [Google Scholar] [CrossRef] [PubMed]
  43. Melbo-Jorgensen, C.; Ness, N.; Andersen, S.; Valkov, A.; Donnem, T.; Al-Saad, S.; Kiselev, Y.; Berg, T.; Nordby, Y.; Bremnes, R.M.; et al. Stromal expression of MiR-21 predicts biochemical failure in prostate cancer patients with Gleason score 6. PLoS ONE 2014, 9, e113039. [Google Scholar] [CrossRef] [PubMed]
  44. Li, T.; Li, R.S.; Li, Y.H.; Zhong, S.; Chen, Y.Y.; Zhang, C.M.; Hu, M.M.; Shen, Z.J. miR-21 as an independent biochemical recurrence predictor and potential therapeutic target for prostate cancer. J. Urol. 2012, 187, 1466–1472. [Google Scholar] [CrossRef]
  45. Zheng, Q.; Peskoe, S.B.; Ribas, J.; Rafiqi, F.; Kudrolli, T.; Meeker, A.K.; De Marzo, A.M.; Platz, E.A.; Lupold, S.E. Investigation of miR-21, miR-141, and miR-221 expression levels in prostate adenocarcinoma for associated risk of recurrence after radical prostatectomy. Prostate 2014, 74, 1655–1662. [Google Scholar] [CrossRef] [Green Version]
  46. Amankwah, E.K.; Anegbe, E.; Park, H.; Pow-Sang, J.; Hakam, A.; Park, J.Y. miR-21, miR-221 and miR-222 expression and prostate cancer recurrence among obese and non-obese cases. Asian J. Androl. 2013, 15, 226–230. [Google Scholar] [CrossRef] [Green Version]
  47. Lichner, Z.; Ding, Q.; Samaan, S.; Saleh, C.; Nasser, A.; Al-Haddad, S.; Samuel, J.N.; Fleshner, N.E.; Stephan, C.; Jung, K.; et al. miRNAs dysregulated in association with Gleason grade regulate extracellular matrix, cytoskeleton and androgen receptor pathways. J. Pathol. 2015, 237, 226–237. [Google Scholar] [CrossRef]
  48. Li, J.; Fu, F.; Wan, X.; Huang, S.; Wu, D.; Li, Y. Up-regulated miR-29c inhibits cell proliferation and glycolysis by inhibiting SLC2A3 expression in prostate cancer. Gene 2018, 665, 26–34. [Google Scholar] [CrossRef]
  49. Kato, M.; Kurozumi, A.; Goto, Y.; Matsushita, R.; Okato, A.; Nishikawa, R.; Fukumoto, I.; Koshizuka, K.; Ichikawa, T.; Seki, N. Regulation of metastasis-promoting LOXL2 gene expression by antitumor microRNAs in prostate cancer. J. Hum. Genet. 2017, 62, 123–132. [Google Scholar] [CrossRef]
  50. Ling, X.H.; Han, Z.D.; Xia, D.; He, H.C.; Jiang, F.N.; Lin, Z.Y.; Fu, X.; Deng, Y.H.; Dai, Q.S.; Cai, C.; et al. MicroRNA-30c serves as an independent biochemical recurrence predictor and potential tumor suppressor for prostate cancer. Mol. Biol. Rep. 2014, 41, 2779–2788. [Google Scholar] [CrossRef]
  51. Ling, X.H.; Chen, Z.Y.; Luo, H.W.; Liu, Z.Z.; Liang, Y.K.; Chen, G.X.; Jiang, F.N.; Zhong, W.D. BCL9, a coactivator for Wnt/-catenin transcription, is targeted by miR-30c and is associated with prostate cancer progression. Oncol. Lett. 2016, 11, 2001–2008. [Google Scholar] [CrossRef] [PubMed]
  52. Kobayashi, N.; Uemura, H.; Nagahama, K.; Okudela, K.; Furuya, M.; Ino, Y.; Ito, Y.; Hirano, H.; Inayama, Y.; Aoki, I.; et al. Identification of miR-30d as a novel prognostic maker of prostate cancer. Oncotarget 2012, 3, 1455–1471. [Google Scholar] [CrossRef] [PubMed]
  53. Lin, Z.Y.; Chen, G.; Zhang, Y.Q.; He, H.C.; Liang, Y.X.; Ye, J.H.; Liang, Y.K.; Mo, R.J.; Lu, J.M.; Zhuo, Y.J.; et al. MicroRNA-30d promotes angiogenesis and tumor growth via MYPT1/c-JUN/VEGFA pathway and predicts aggressive outcome in prostate cancer. Mol. Cancer 2017, 16, 48. [Google Scholar] [CrossRef] [PubMed]
  54. Schaefer, A.; Jung, M.; Mollenkopf, H.-J.; Wagner, I.; Stephan, C.; Jentzmik, F.; Miller, K.; Lein, M.; Kristiansen, G.; Jung, K. Diagnostic and prognostic implications of microRNA profiling in prostate carcinoma. Int. J. Cancer 2010, 126, 1166–1176. [Google Scholar] [CrossRef]
  55. Daniunaite, K.; Dubikaityte, M.; Gibas, P.; Bakavicius, A.; Rimantas, L.J.; Ulys, A.; Jankevicius, F.; Jarmalaite, S. Clinical significance of miRNA host gene promoter methylation in prostate cancer. Hum. Mol. Genet. 2017, 26, 2451–2461. [Google Scholar] [CrossRef]
  56. Liu, X.; Luo, X.; Wu, Y.; Xia, D.; Chen, W.; Fang, Z.; Deng, J.; Hao, Y.; Yang, X.; Zhang, T.; et al. MicroRNA-34a attenuates paclitaxel resistance in prostate cancer cells via direct suppression of JAG1/Notch1 Axis. Cell Physiol. Biochem. 2018, 50, 261–276. [Google Scholar] [CrossRef]
  57. Song, C.J.; Chen, H.; Chen, L.Z.; Ru, G.M.; Guo, J.J.; Ding, Q.N. The potential of microRNAs as human prostate cancer biomarkers: A meta-analysis of related studies. J. Cell Biochem. 2018, 119, 2763–2786. [Google Scholar] [CrossRef]
  58. Avgeris, M.; Stravodimos, K.; Fragoulis, E.G.; Scorilas, A. The loss of the tumour-suppressor miR-145 results in the shorter disease-free survival of prostate cancer patients. Br. J. Cancer 2013, 108, 2573–2581. [Google Scholar] [CrossRef] [Green Version]
  59. Kang, S.G.; Ha, Y.R.; Kim, S.J.; Kang, S.H.; Park, H.S.; Lee, J.G.; Cheon, J.; Kim, C.H. Do microRNA 96, 145 and 221 expressions really aid in the prognosis of prostate carcinoma? Asian J. Androl. 2012, 14, 752–757. [Google Scholar] [CrossRef] [Green Version]
  60. Larne, O.; Martens-Uzunova, E.; Hagman, Z.; Edsjo, A.; Lippolis, G.; den Berg, M.S.; Bjartell, A.; Jenster, G.; Ceder, Y. miQ—A novel microRNA based diagnostic and prognostic tool for prostate cancer. Int. J. Cancer 2013, 132, 2867–2875. [Google Scholar] [CrossRef]
  61. Leite, K.R.; Tomiyama, A.; Reis, S.T.; Sousa-Canavez, J.M.; Sanudo, A.; Dall’Oglio, M.F.; Camara-Lopes, L.H.; Srougi, M. MicroRNA-100 expression is independently related to biochemical recurrence of prostate cancer. J. Urol. 2011, 185, 1118–1122. [Google Scholar] [CrossRef] [PubMed]
  62. Zhang, S.; Li, Z.; Zhang, L.; Xu, Z. MEF2 activated long noncoding RNA PCGEM1 promotes cell proliferation in hormonerefractory prostate cancer through downregulation of miR148a. Mol. Med. Rep. 2018, 18, 202–208. [Google Scholar] [PubMed]
  63. Guo, J.; Wang, M.; Liu, X. MicroRNA-195 suppresses tumor cell proliferation and metastasis by directly targeting BCOX1 in prostate carcinoma. J. Exp. Clin. Cancer Res. 2015, 34, 91. [Google Scholar] [CrossRef] [PubMed]
  64. Cai, C.; Chen, Q.B.; Han, Z.D.; Zhang, Y.Q.; He, H.C.; Chen, J.H.; Chen, Y.R.; Yang, S.B.; Wu, Y.D.; Zeng, Y.R.; et al. miR-195 inhibits tumor progression by targeting RPS6KB1 in human prostate cancer. Clin. Cancer Res. 2015, 21, 4922–4934. [Google Scholar] [CrossRef] [PubMed]
  65. Todorova, K.; Metodiev, M.V.; Metodieva, G.; Zasheva, D.; Mincheff, M.; Hayrabedyan, S. miR-204 is dysregulated in metastatic prostate cancer in vitro. Mol. Carcinog. 2016, 55, 131–147. [Google Scholar] [CrossRef] [PubMed]
  66. Lin, Y.; Chen, F.; Shen, L.; Tang, X.; Du, C.; Sun, Z.; Ding, H.; Chen, J.; Shen, B. Biomarker microRNAs for prostate cancer metastasis: Screened with a network vulnerability analysis model. J. Transl. Med. 2018, 16, 134. [Google Scholar] [CrossRef]
  67. Spahn, M.; Kneitz, S.; Scholz, C.J.; Stenger, N.; Rudiger, T.; Strobel, P.; Riedmiller, H.; Kneitz, B. Expression of microRNA-221 is progressively reduced in aggressive prostate cancer and metastasis and predicts clinical recurrence. Int. J. Cancer 2010, 127, 394–403. [Google Scholar]
  68. Lin, Z.Y.; Huang, Y.Q.; Zhang, Y.Q.; Han, Z.D.; He, H.C.; Ling, X.H.; Fu, X.; Dai, Q.S.; Cai, C.; Chen, J.H.; et al. MicroRNA-224 inhibits progression of human prostate cancer by downregulating TRIB1. Int. J. Cancer 2014, 135, 541–550. [Google Scholar] [CrossRef]
  69. Wan, Y.; Zeng, Z.C.; Xi, M.; Wan, S.; Hua, W.; Liu, Y.L.; Zhou, Y.L.; Luo, H.W.; Jiang, F.N.; Zhong, W.D. Dysregulated microRNA-224/apelin axis associated with aggressive progression and poor prognosis in patients with prostate cancer. Hum. Pathol. 2015, 46, 295–303. [Google Scholar] [CrossRef]
  70. Nam, R.K.; Benatar, T.; Wallis, C.J.; Amemiya, Y.; Yang, W.; Garbens, A.; Naeim, M.; Sherman, C.; Sugar, L.; Seth, A. MiR-301a regulates E-cadherin expression and is predictive of prostate cancer recurrence. Prostate 2016, 76, 869–884. [Google Scholar] [CrossRef]
  71. He, H.C.; Han, Z.D.; Dai, Q.S.; Ling, X.H.; Fu, X.; Lin, Z.Y.; Deng, Y.H.; Qin, G.Q.; Cai, C.; Chen, J.H.; et al. Global analysis of the differentially expressed miRNAs of prostate cancer in Chinese patients. BMC Genom. 2013, 14, 757. [Google Scholar] [CrossRef] [PubMed]
  72. Bian, H.; Zhou, Y.; Zhou, D.; Zhang, Y.; Shang, D.; Qi, J. The latest progress on miR-374 and its functional implications in physiological and pathological processes. J. Cell. Mol. Med. 2019, 23, 3063–3076. [Google Scholar] [CrossRef] [PubMed]
  73. Zhu, J.; Sun, C.; Wang, L.; Xu, M.; Zang, Y.; Zhou, Y.; Liu, X.; Tao, W.; Xue, B.; Shan, Y.; et al. Targeting survivin using a combination of miR494 and survivin shRNA has synergistic effects on the suppression of prostate cancer growth. Mol. Med. Rep. 2016, 13, 1602–1610. [Google Scholar] [CrossRef] [PubMed]
  74. Shen, P.F.; Chen, X.Q.; Liao, Y.C.; Chen, N.; Zhou, Q.; Wei, Q.; Li, X.; Wang, J.; Zeng, H. MicroRNA-494-3p targets CXCR4 to suppress the proliferation, invasion, and migration of prostate cancer. Prostate 2014, 74, 756–767. [Google Scholar] [CrossRef] [PubMed]
  75. Guo, K.; Liang, Z.; Li, F.; Wang, H. Comparison of miRNA and gene expression profiles between metastatic and primary prostate cancer. Oncol. Lett. 2017, 14, 6085–6090. [Google Scholar] [CrossRef] [PubMed]
  76. Pashaei, E.; Pashaei, E.; Ahmady, M.; Ozen, M.; Aydin, N. Meta-analysis of miRNA expression profiles for prostate cancer recurrence following radical prostatectomy. PLoS ONE 2017, 12, e0179543. [Google Scholar] [CrossRef] [PubMed]
  77. Heagerty, P.J.; Zheng, Y. Survival model predictive accuracy and ROC curves. Biometrics 2005, 61, 92–105. [Google Scholar] [CrossRef]
  78. Steyerberg, E.W.; Pencina, M.J.; Lingsma, H.F.; Kattan, M.W.; Vickers, A.J.; Van, C.B. Assessing the incremental value of diagnostic and prognostic markers: A review and illustration. Eur. J. Clin. Investig. 2012, 42, 216–228. [Google Scholar] [CrossRef]
  79. Huang, Y.Q.; Ling, X.H.; Yuan, R.Q.; Chen, Z.Y.; Yang, S.B.; Huang, H.X.; Zhong, W.D.; Qiu, S.P. miR30c suppresses prostate cancer survival by targeting the ASF/SF2 splicing factor oncoprotein. Mol. Med. Rep. 2017, 16, 2431–2438. [Google Scholar] [CrossRef]
  80. Zhang, J.; Wang, X.; Wang, Y.; Peng, R.; Lin, Z.; Wang, Y.; Hu, B.; Wang, J.; Shi, G. Low expression of microRNA-30c promotes prostate cancer cells invasion involved in downregulation of KRAS protein. Oncol. Lett. 2017, 14, 363–368. [Google Scholar] [CrossRef]
  81. Zhang, Q.; Sun, M.; Zhou, S.; Guo, B. Class I HDAC inhibitor mocetinostat induces apoptosis by activation of miR-31 expression and suppression of E2F6. Cell Death Discov. 2016, 2, 16036. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  82. Coarfa, C.; Fiskus, W.; Eedunuri, V.K.; Rajapakshe, K.; Foley, C.; Chew, S.A.; Shah, S.S.; Geng, C.; Shou, J.; Mohamed, J.S.; et al. Comprehensive proteomic profiling identifies the androgen receptor axis and other signaling pathways as targets of microRNAs suppressed in metastatic prostate cancer. Oncogene 2016, 35, 2345–2356. [Google Scholar] [CrossRef] [PubMed]
  83. Zhang, L.; Qi, M.; Feng, T.; Hu, J.; Wang, L.; Li, X.; Gao, W.; Liu, H.; Jiao, M.; Wu, Z.; et al. IDH1R132H promotes malignant transformation of benign prostatic epithelium by dysregulating microRNAs: Involvement of IGF1R-AKT/STAT3 signaling pathway. Neoplasia 2018, 20, 207–217. [Google Scholar] [CrossRef] [PubMed]
  84. Huang, S.; Wa, Q.; Pan, J.; Peng, X.; Ren, D.; Huang, Y.; Chen, X.; Tang, Y. Downregulation of miR-141-3p promotes bone metastasis via activating NF-kappaB signaling in prostate cancer. J. Exp. Clin. Cancer Res. 2017, 36, 173. [Google Scholar] [CrossRef] [PubMed]
  85. Li, J.Z.; Li, J.; Wang, H.Q.; Li, X.; Wen, B.; Wang, Y.J. MiR-141-3p promotes prostate cancer cell proliferation through inhibiting kruppel-like factor-9 expression. Biochem. Biophys. Res. Commun. 2017, 482, 1381–1386. [Google Scholar] [CrossRef]
  86. Zhao, H.; Su, W.; Zhu, C.; Zeng, T.; Yang, S.; Wu, W.; Wang, D. Cell fate regulation by reticulon-4 in human prostate cancers. J. Cell. Physiol. 2019, 234, 10372–10385. [Google Scholar] [CrossRef]
  87. Murata, T.; Takayama, K.; Katayama, S.; Urano, T.; Horie-Inoue, K.; Ikeda, K.; Takahashi, S.; Kawazu, C.; Hasegawa, A.; Ouchi, Y.; et al. miR-148a is an androgen-responsive microRNA that promotes LNCaP prostate cell growth by repressing its target CAND1 expression. Prostate Cancer Prostatic Dis. 2010, 13, 356–361. [Google Scholar] [CrossRef]
  88. Fujita, Y.; Kojima, K.; Ohhashi, R.; Hamada, N.; Nozawa, Y.; Kitamoto, A.; Sato, A.; Kondo, S.; Kojima, T.; Deguchi, T.; et al. MiR-148a attenuates paclitaxel resistance of hormone-refractory, drug-resistant prostate cancer PC3 cells by regulating MSK1 expression. J. Biol. Chem. 2010, 285, 19076–19084. [Google Scholar] [CrossRef]
  89. Farina, N.H.; Zingiryan, A.; Akech, J.A.; Callahan, C.J.; Lu, H.; Stein, J.L.; Languino, L.R.; Stein, G.S.; Lian, J.B. A microRNA/Runx1/Runx2 network regulates prostate tumor progression from onset to adenocarcinoma in TRAMP mice. Oncotarget 2016, 7, 70462–70474. [Google Scholar] [CrossRef] [Green Version]
  90. Song, C.; Chen, H.; Wang, T.; Zhang, W.; Ru, G.; Lang, J. Expression profile analysis of microRNAs in prostate cancer by next-generation sequencing. Prostate 2015, 75, 500–516. [Google Scholar] [CrossRef]
  91. Klein, E.A.; Cooperberg, M.R.; Magi-Galluzzi, C.; Simko, J.P.; Falzarano, S.M.; Maddala, T.; Chan, J.M.; Li, J.; Cowan, J.E.; Tsiatis, A.C.; et al. A 17-gene assay to predict prostate cancer aggressiveness in the context of Gleason grade heterogeneity, tumor multifocality, and biopsy undersampling. Eur. Urol. 2014, 66, 550–560. [Google Scholar] [CrossRef] [PubMed]
  92. Lee, J.W.; Devanarayan, V.; Barrett, Y.C.; Weiner, R.; Allinson, J.; Fountain, S.; Keller, S.; Weinryb, I.; Green, M.; Duan, L.; et al. Fit-for-purpose method development and validation for successful biomarker measurement. Pharm. Res. 2006, 23, 312–328. [Google Scholar] [CrossRef] [PubMed]
  93. Lu, H.; Busch, J.; Jung, M.; Rabenhorst, S.; Ralla, B.; Kilic, E.; Mergemeier, S.; Budach, N.; Fendler, A.; Jung, K. Diagnostic and prognostic potential of circulating cell-free genomic and mitochondrial DNA fragments in clear cell renal cell carcinoma patients. Clin. Chim. Acta 2016, 452, 109–119. [Google Scholar] [CrossRef]
  94. Bhatnagar, N.; Li, X.; Padi, S.K.; Zhang, Q.; Tang, M.S.; Guo, B. Downregulation of miR-205 and miR-31 confers resistance to chemotherapy-induced apoptosis in prostate cancer cells. Cell Death Dis. 2010, 1, e105. [Google Scholar] [CrossRef] [PubMed]
  95. Schmidt, L.; Fredsoe, J.; Kristensen, H.; Strand, S.H.; Rasmussen, A.; Hoyer, S.; Borre, M.; Mouritzen, P.; Orntoft, T.; Sorensen, K.D. Training and validation of a novel 4-miRNA ratio model (MiCaP) for prediction of postoperative outcome in prostate cancer patients. Ann. Oncol. 2018, 29, 2003–2009. [Google Scholar] [CrossRef] [PubMed]
  96. Strand, S.H.; Bavafaye-Haghighi, E.; Kristensen, H.; Rasmussen, A.K.; Hoyer, S.; Borre, M.; Mouritzen, P.; Besenbacher, S.; Orntoft, T.F.; Sorensen, K.D. A novel combined miRNA and methylation marker panel (miMe) for prediction of prostate cancer outcome after radical prostatectomy. Int. J. Cancer 2019. [Google Scholar] [CrossRef] [PubMed]
  97. Nam, R.K.; Wallis, C.J.D.; Amemiya, Y.; Benatar, T.; Seth, A. Identification of a novel microRNA panel associated with metastasis following radical prostatectomy for prostate cancer. Anticancer Res. 2018, 38, 5027–5034. [Google Scholar] [CrossRef]
  98. Gerszten, R.E.; Wang, T.J. The search for new cardiovascular biomarkers. Nature 2008, 451, 949–952. [Google Scholar] [CrossRef]
  99. Mohler, J.L.; Antonarakis, E.S.; Armstrong, A.J.; D’Amico, A.V.; Davis, B.J.; Dorff, T.; Eastham, J.A.; Enke, C.A.; Farrington, T.A.; Higano, C.S.; et al. Prostate Cancer, Version 2.2019, NCCN Clinical Practice Guidelines in Oncology. J. Natl. Compr. Cancer Netw. 2019, 17, 479–505. [Google Scholar] [CrossRef] [Green Version]
  100. Loeb, S.; Ross, A.E. Genomic testing for localized prostate cancer: Where do we go from here? Curr. Opin. Urol. 2017, 27, 495–499. [Google Scholar] [CrossRef]
  101. Glass, A.G.; Leo, M.C.; Haddad, Z.; Yousefi, K.; du, P.M.; Chen, C.; Choeurng, V.; Abdollah, F.; Robbins, B.; Ra, S.; et al. Validation of a genomic classifier for predicting post-prostatectomy recurrence in a community based health care setting. J. Urol. 2016, 195, 1748–1753. [Google Scholar] [CrossRef] [PubMed]
  102. Tilki, D.; Evans, C.P. The Decipher Genomic Classifier independently improves prognostication for patients after prostatectomy. Eur. Urol. 2018, 73, 176–177. [Google Scholar] [CrossRef] [PubMed]
  103. Alam, S.; Tortora, J.; Staff, I.; McLaughlin, T.; Wagner, J. Prostate cancer genomics: Comparing results from three molecular assays. Can. J. Urol. 2019, 26, 9758–9762. [Google Scholar] [PubMed]
  104. Jung, M.; Schaefer, A.; Steiner, I.; Kempkensteffen, C.; Stephan, C.; Erbersdobler, A.; Jung, K. Robust microRNA stability in degraded RNA preparations from human tissue and cell samples. Clin. Chem. 2010, 56, 998–1006. [Google Scholar] [CrossRef] [PubMed]
  105. Littman, B.H.; Di, M.L.; Plebani, M.; Marincola, F.M. What’s next in translational medicine? Clin. Sci. Lond. 2007, 112, 217–227. [Google Scholar] [CrossRef] [PubMed]
  106. Fendler, A.; Stephan, C.; Yousef, G.M.; Jung, K. MicroRNAs as regulators of signal transduction in urological tumors. Clin. Chem. 2011, 57, 954–968. [Google Scholar] [CrossRef] [PubMed]
  107. Waltering, K.K.; Porkka, K.P.; Jalava, S.E.; Urbanucci, A.; Kohonen, P.J.; Latonen, L.M.; Kallioniemi, O.P.; Jenster, G.; Visakorpi, T. Androgen regulation of micro-RNAs in prostate cancer. Prostate 2011, 71, 604–614. [Google Scholar] [CrossRef]
  108. Salmena, L.; Poliseno, L.; Tay, Y.; Kats, L.; Pandolfi, P.P. A ceRNA hypothesis: The Rosetta Stone of a hidden RNA language? Cell 2011, 146, 353–358. [Google Scholar] [CrossRef]
  109. Burke, H.B. Predicting clinical outcomes using molecular biomarkers. Biomark. Cancer 2016, 8, 89–99. [Google Scholar] [CrossRef]
  110. Van den Broeck, T.; van den Bergh, R.C.N.; Briers, E.; Cornford, P.; Cumberbatch, M.; Tilki, D.; De, S.M.; Fanti, S.; Fossati, N.; Gillessen, S.; et al. Biochemical recurrence in prostate cancer: The European Association of Urology Prostate Cancer Guidelines Panel Recommendations. Eur. Urol. Focus 2019, in press. [Google Scholar] [CrossRef]
  111. Fendler, A.; Jung, M.; Stephan, C.; Erbersdobler, A.; Jung, K.; Yousef, G.M. The antiapoptotic function of miR-96 in prostate cancer by inhibition of FOXO1. PLoS ONE 2013, 8, e80807. [Google Scholar] [CrossRef] [PubMed]
  112. Schaefer, A.; Jung, M.; Stephan, C.; Miller, K.; Lein, M.; Kristiansen, G.; Erbersdobler, A.; Jung, K. Suitable reference genes for relative quantification of miRNA expression in prostate cancer. Exp. Mol. Med. 2010, 42, 749–758. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  113. Bustin, S.A.; Benes, V.; Garson, J.A.; Hellemans, J.; Huggett, J.; Kubista, M.; Mueller, R.; Nolan, T.; Pfaffl, M.W.; Shipley, G.L.; et al. The MIQE guidelines: Minimum information for publication of quantitative real-time PCR experiments. Clin. Chem. 2009, 55, 611–622. [Google Scholar] [CrossRef] [PubMed]
  114. Stephan, C.; Jung, K.; Semjonow, A.; Schulze-Forster, K.; Cammann, H.; Hu, X.; Meyer, H.A.; Bogemann, M.; Miller, K.; Friedersdorff, F. Comparative assessment of urinary prostate cancer antigen 3 and TMPRSS2:ERG gene fusion with the serum [−2] proprostate-specific antigen-based prostate health index for detection of prostate cancer. Clin. Chem. 2013, 59, 280–288. [Google Scholar] [CrossRef]
  115. Vickers, A.J.; Elkin, E.B. Decision curve analysis: A novel method for evaluating prediction models. Med. Decis. Mak. 2006, 26, 565–574. [Google Scholar] [CrossRef]
Cancers 11 01603 g001 550
Figure 1. Flowchart of the study. Fourteen and 15 out of 29 adjacent normal samples were from patients without and with BCR, respectively. Between these two sample groups, expression of all microRNAs of interest in this study shown in Figure 2 did not significantly differ (Mann-Whitney U-test, p-values between 0.198 and 0.983) and the combined 29 samples were therefore used as “Normal adjacent tissue”. Abbreviations: FFPE, formalin-fixed, paraffin-embedded; PCa, prostate carcinoma; BCR, biochemical recurrence; RT-qPCR, reverse transcription real time quantitative polymerase chain reaction; MIQE, Minimum Information for Publication of Quantitative Real-time PCR Experiments.
Figure 1. Flowchart of the study. Fourteen and 15 out of 29 adjacent normal samples were from patients without and with BCR, respectively. Between these two sample groups, expression of all microRNAs of interest in this study shown in Figure 2 did not significantly differ (Mann-Whitney U-test, p-values between 0.198 and 0.983) and the combined 29 samples were therefore used as “Normal adjacent tissue”. Abbreviations: FFPE, formalin-fixed, paraffin-embedded; PCa, prostate carcinoma; BCR, biochemical recurrence; RT-qPCR, reverse transcription real time quantitative polymerase chain reaction; MIQE, Minimum Information for Publication of Quantitative Real-time PCR Experiments.
Cancers 11 01603 g001
Cancers 11 01603 g002a 550Cancers 11 01603 g002b 550
Figure 2. Expression of 20 microRNAs examined in adjacent normal tissue samples and prostate cancer (PCa) tissue samples from patients with and without biochemical recurrence (BCR) are presented in subfigures (a–t). Normalized miRNA expression data are given as box- and whisker plots. Boxes represent the lower and upper quartiles with medians; whiskers illustrate the 10 and 90 percentiles of the cohorts. Kruskal–Wallis test with Dunn’s post-hoc was performed. The p-value of the total Kruskal–Wallis test is indicated in the upper left-hand corner of the figure. Abbreviations: N, adjacent normal tissue (n = 29) as explained in legend of Figure 1; R (−), PCa tissue samples from patients without BCR (n = 108); R (+), PCa tissue samples from patients with BCR (n = 98).
Figure 2. Expression of 20 microRNAs examined in adjacent normal tissue samples and prostate cancer (PCa) tissue samples from patients with and without biochemical recurrence (BCR) are presented in subfigures (a–t). Normalized miRNA expression data are given as box- and whisker plots. Boxes represent the lower and upper quartiles with medians; whiskers illustrate the 10 and 90 percentiles of the cohorts. Kruskal–Wallis test with Dunn’s post-hoc was performed. The p-value of the total Kruskal–Wallis test is indicated in the upper left-hand corner of the figure. Abbreviations: N, adjacent normal tissue (n = 29) as explained in legend of Figure 1; R (−), PCa tissue samples from patients without BCR (n = 108); R (+), PCa tissue samples from patients with BCR (n = 98).
Cancers 11 01603 g002aCancers 11 01603 g002b
Cancers 11 01603 g003 550
Figure 3. Time-dependent receiver-operating characteristics (ROC) curve analyses of (a) separate and (b) combined miRNA panel-based and clinicopathological factors-based models at different postoperative time points. Areas under the time-dependent ROC curve (AUCs) of all models were computed based on a cumulative case/dynamic control approach [77] at postoperative 1, 2, 3, 4, and 5 years, respectively. The models of miRNA-panels are described in Table 4 and Table S7. The four reference clinical models [15,16,17,18] and our full model are explained in Table 6 and Table 7. Abbreviations: CAPRAS, Cancer of the Prostate Risk Assessment Postsurgical Score, calculated according to [16]; STEPH, calculated according to Stephenson et al. [17]; NCCN, National Comprehensive Cancer Network, calculated according to [18]; D’AMICO, calculated according to D’Amico et al. [15]; Full model, according to the Cox regression model described in Table 6 with all clinicopathological factors except for digital rectal examination.
Figure 3. Time-dependent receiver-operating characteristics (ROC) curve analyses of (a) separate and (b) combined miRNA panel-based and clinicopathological factors-based models at different postoperative time points. Areas under the time-dependent ROC curve (AUCs) of all models were computed based on a cumulative case/dynamic control approach [77] at postoperative 1, 2, 3, 4, and 5 years, respectively. The models of miRNA-panels are described in Table 4 and Table S7. The four reference clinical models [15,16,17,18] and our full model are explained in Table 6 and Table 7. Abbreviations: CAPRAS, Cancer of the Prostate Risk Assessment Postsurgical Score, calculated according to [16]; STEPH, calculated according to Stephenson et al. [17]; NCCN, National Comprehensive Cancer Network, calculated according to [18]; D’AMICO, calculated according to D’Amico et al. [15]; Full model, according to the Cox regression model described in Table 6 with all clinicopathological factors except for digital rectal examination.
Cancers 11 01603 g003
Cancers 11 01603 g004 550
Figure 4. Decision curve analysis to demonstrate the benefit of the inclusion of the five-miR-panel into only clinicopathological variables-based biochemical recurrence predictive tools, according to (a) D’Amico et al. [15], (b) Cancer of the Prostate Risk Assessment Postsurgical Score (CAPRAS) [16], and (c) the “full model” used in this study (Table 6).
Figure 4. Decision curve analysis to demonstrate the benefit of the inclusion of the five-miR-panel into only clinicopathological variables-based biochemical recurrence predictive tools, according to (a) D’Amico et al. [15], (b) Cancer of the Prostate Risk Assessment Postsurgical Score (CAPRAS) [16], and (c) the “full model” used in this study (Table 6).
Cancers 11 01603 g004
Table
Table 1. Clinicopathological characteristics of the study group.
Table 1. Clinicopathological characteristics of the study group.
CharacteristicsAll PatientsPatients with Biochemical RecurrencePatients without Biochemical Recurrencep-value a
Patients, no. (%)206 (100)98 (48)108 (52)
Age, median years (range)63 (47–74)64 (51–74)62 (47–73)0.014
PSA, median µg/L (range)8.4 (1.3–50.6)9.3 (1.4–50.6)7.5 (1.3–32.9)0.003
Prostate volume, median cm3 (range)33 (14–130)33 (15–130)32 (14–120)0.465
DRE, no. (%) 0.260
Non-suspicious122 (59)54 (55)68 (63)
Suspicious84 (41)44 (45)40 (37)
pT status, no. (%) <0.0001
pT2a18 (9)5 (5)13 (12)
pT2b28 (14)6 (6)22 (20)
pT2c76 (37)31 (32)45 (42)
pT3a62 (30)40 (41) 22 (20)
pT3b21 (10)15 (15)6 (6)
pT41 (0.5)1 (1)0
ISUP Grade groups, no. (%)
152 (25)14 (14)38 (35)0.0003
268 (33)30 (31)38 (35)
327 (13)19 (19)8 (7)
429 (14)16 (16)13 (12)
530 (15)19 (19)11 (10)
pN status, no. (%) b
pN0113 (55)57 (58)56 (52)0.207
pN110 (5)7 (7)3 (3)
pNx83 (40)34 (35)49 (45)
Surgical margin, no. (%)
Negative133 (65)54 (55)79 (73)0.009
Positive73 (35)44 (45)29 (27)
Follow-up after surgery
Median months (range)108 (17–180)101 (19–160)121 (17–180)<0.0001
Recurrence free survival
Median months (95% CI)52 (46–60)16 (11-21)80 (75–93)<0.0001
Abbreviations: CI, confidence interval; DRE, digital rectal examination; ISUP Grade groups, histopathological grade system based on Gleason score according to the International Society of Urologic Pathology; pN, lymph node status; PSA, prostate specific antigen; pT, pathological tumor classification. a p-values (Mann-Whitney U test; Chi-square Cochran-Armitage test for trend or Fisher’s exact test) indicate the association of the patients with and without biochemical recurrence with the clinicopathological variables. b In the following evaluation of data, with regard to the pN status, patients with pNx (with no lymph node dissection) were considered node negative, and were combined with pN0 patients into one group.
Table
Table 2. MicroRNAs analyzed in this study. The selection of microRNAs was based on a previous systematic review of studies with microRNAs as potential BCR predictors after RP [34] and/or on other study data after the publication of the review.
Table 2. MicroRNAs analyzed in this study. The selection of microRNAs was based on a previous systematic review of studies with microRNAs as potential BCR predictors after RP [34] and/or on other study data after the publication of the review.
miRBase ID
Release 22		miRBase Accession No. aReference to BCR-Related miRNADifferential Expression: Tumor vs. Normal Tissue bDifferential Expression: Recurrence vs. Non-Recurrence c
hsa-miR-1-3pMIMAT0000416[35,40,41]n.s., p = 0.177↓, p = 0.0004
hsa-miR-21-5pMIMAT0000076[42,43,44,45,46]↑, p = 0.006n.s., p = 0.271
hsa-miR-29c-3pMIMAT0000681[47,48,49]n.s., p = 0.272n.s., p = 0.390
hsa-miR-30c-5pMIMAT0000244[50,51]↑, p = 0.011n.s., p = 0.056
hsa-miR-30d-5pMIMAT0000245[52,53]↑, p < 0.0001n.s., p = 0.297
hsa-miR-31-5pMIMAT0000089[54,55]↓, p < 0.0001↓, p = 0.0007
hsa-miR-34a-5pMIMAT0000255[47,56,57]n.s., p = 0.155↓, p = 0.0004
hsa-miR-141-3pMIMAT0000432[45,47]n.s., p = 0.219n.s., p = 0.099
hsa-miR-145-5pMIMAT0000437[58,59,60,61]↓, p = 0.011↓, p = 0.0004
hsa-miR-148a-3pMIMAT0000243[47,62]↑, p = 0.018n.s., p = 0.809
hsa-miR-185-5pMIMAT0000455[33]↑, p = 0.015n.s., p = 0.420
hsa-miR-195-5pMIMAT0000461[63,64]↑, p < 0.0001n.s., p = 0.559
hsa-miR-204-5pMIMAT0000265[65,66]↓, p = 0.0001↓, p = 0.045
hsa-miR-221-3pMIMAT0000278[33,45,59,60,67]↓, p < 0.0001↓, p = 0.0001
hsa-miR-224-5p dMIMAT0000281[68,69]↓, p < 0.0001↓, p = 0.0001
hsa-miR-301a-3pMIMAT0000688[38,70]↑, p = 0.0002n.s., p = 0.674
hsa-miR-326MIMAT0000756[33]↑, p = 0.009↓, p = 0.033
hsa-miR-374b-5pMIMAT0004955[71,72]↑, p < 0.0001n.s., p = 0.141
hsa-miR-494-3pMIMAT0002816[73,74,75]n.s., p = 0.721n.s., p = 0.821
hsa-miR-939-5pMIMAT0004982[76]↑, p = 0.028n.s., p = 0.784
Abbreviations: BCR, biochemical recurrence; RP, radical prostatectomy; HR, hazard ratio; ↑, upregulated; ↓, downregulated; n.s., not significant. a Further details about the characteristics of the analyzed microRNAs (location on chromosomes, miRNA families, clustering with other miRNA, and assays) are summarized in Table S1. b Expression difference (Mann–Whitney test) in all tumor samples (n = 206) vs. controls of adjacent normal tissue samples as controls (n = 29), see also Figure S1. c Expression difference (Kruskal–Wallis test with Dunn’s post-hoc test, see also Figure 2) in tissue samples from BCR patients vs. those from non-BCR patients. d Until miRBase 21; however, this miRNA has now the miRBase ID eca-miR-224 (accession no. MIMAT0013206).
Table
Table 3. Spearman rank correlation coefficients between various clinicopathological factors and the examined microRNAs in all tumor samples (n = 206). The eight significantly differentially expressed miRNAs, according to biochemical recurrence status, and correlations with p < 0.05, are highlighted in gray.
Table 3. Spearman rank correlation coefficients between various clinicopathological factors and the examined microRNAs in all tumor samples (n = 206). The eight significantly differentially expressed miRNAs, according to biochemical recurrence status, and correlations with p < 0.05, are highlighted in gray.
(A) miR-1 to miR-1 to miR148a
miR-1miR-21miR-29cmiR-30cmiR-30dmiR-31miR-34amiR-141miR-145miR-148a
Agers−0.059−0.0870.1040.0090.073−0.0940.0700.095−0.0360.042
p-value0.4030.2140.1390.9010.2940.1780.3180.1740.6040.546
PSArs−0.2180.0070.030−0.1070.019−0.086−0.1190.156−0.2110.068
p-value0.0020.9220.6680.1250.7830.2200.0900.0250.0020.334
DRErs−0.1190.132−0.035−0.1270.021−0.053−0.031−0.022−0.0620.005
p-value0.0870.0590.6220.0680.7650.4460.6570.7570.3760.945
Marginrs−0.146−0.1290.056−0.139−0.064−0.198−0.1140.112−0.1840.184
p-value0.0370.0650.4250.0470.3580.0040.1020.1080.0080.008
pN statusrs−0.1100.022−0.147−0.122−0.091−0.129−0.142−0.046−0.1100.005
p-value0.1160.7530.0350.0800.1940.0640.0410.5080.1140.946
pT stagers−0.2180.1460.029−0.0860.015−0.085−0.0610.153−0.1790.144
p-value0.0020.0360.6770.2220.8280.2250.3850.0280.0100.038
ISUPrs−0.3510.1980.188−0.0800.095−0.232−0.0140.298−0.2650.249
p-value0.00010.0040.0070.2560.1730.0010.8440.00010.00010.0003
(B) miR-185 to miR-939
miR-185miR-195miR-204miR-221miR-224miR-301amiR-326miR-374bmiR-494miR-939
Agers−0.0020.0110.1190.001−0.088−0.0440.045−0.0360.0380.000
p-value0.9760.8700.0900.9830.2080.5280.5170.6110.5840.999
PSArs0.063−0.137−0.220−0.080−0.137−0.0300.048−0.0630.005−0.027
p-value0.3710.0500.0020.2530.0490.6700.4950.3670.9450.702
DRErs0.0120.033−0.062−0.082−0.0380.0010.000−0.0340.071−0.003
p-value0.8600.6340.3730.2440.5880.9900.9940.6320.3110.966
Marginrs−0.096−0.214−0.124−0.032−0.175−0.080−0.092−0.0960.127−0.030
p-value0.1690.0020.0760.6520.0120.2540.1860.1680.0690.672
pN statusrs−0.0200.031−0.110−0.048−0.1060.034−0.080−0.097−0.031−0.130
p-value0.7740.6590.1160.4960.1310.6310.2560.1660.6590.062
pT stagers0.0880.032−0.100−0.083−0.1660.121−0.023−0.064−0.049−0.065
p-value0.2110.6480.1540.2350.0170.0840.7380.3610.4850.352
ISUPrs0.092−0.042−0.256−0.015−0.1460.0510.066−0.0590.056−0.069
p-value0.1860.5510.00020.8250.0360.4670.3460.4020.4230.323
Abbreviations: miR, microRNA in its abbreviated form to facilitate the readability of the table; full annotation in Table 2; rs, Spearman rank correlation coefficient; PSA, prostate specific antigen; DRE, digital rectal examination result; pN status, pathological lymph node status, positive/negative; pT stage, pathological tumor classification, see Table 1; ISUP, histopathological grade system, see Table 1.
Table
Table 4. Construction of a miRNA-based predictive classifier for biochemical recurrence, using a bootstrapping approach of Cox regression analysis with all 206 tumor samples. a
Table 4. Construction of a miRNA-based predictive classifier for biochemical recurrence, using a bootstrapping approach of Cox regression analysis with all 206 tumor samples. a
miRNAUnivariate Cox Regression for All miRsMultivariate Cox Regression with Significant Univariate miRs b
Full ModelBackward Elimination
HR (95% CI)p-valueHR (95% CI)p-valueHR (95% CI)p-value
miR-1-3p0.67 (0.52–0.87)0.0030.82 (0.58–1.16)0.262
miR-21-5p1.12 (1.01–1.24)0.0491.35 (0.86–2.12)0.188
miR-29c-3p1.06 (0.81–1.39)0.660
miR-30c-5p0.56 (0.34–0.92)0.0230.46 (0.23–0.93)0.0310.49 (0.28–0.85)0.011
miR-30d-5p1.22 (1.05–1.42)0.0091.79 (1.12–2.85)0.0151.27 (0.98–1.66)0.075
miR-31-5p0.78 (0.68–0.90)0.0010.83 (0.69–0.98)0.0300.78 (0.67–0.91)0.001
miR-34a-5p1.09 (0.96–1.23)0.1740.61 (0.36–1.03)0.063
miR-141-3p1.25 (1.10–1.41)0.0011.96 (1.25–3.07)0.0031.92 (1.32–2.79)0.001
miR-145-5p0.68 (0.51–0.90)0.0081.17 (0.65–2.11)0.604
miR-148a-3p1.15 (0.99–1.34)0.0640.58 (0.35–0.95)0.0310.60 (0.44–0.81)0.001
miR-185-5p1.06 (1.01–1.12)0.0240.99 (0.68–1.43)0.955
miR-195-5p0.91 (0.68–1.23)0.553
miR-204-5p0.76 (0.59–0.98)0.0331.31 (0.92–1.88)0.138
miR-221-3p0.81 (0.67–0.98)0.0330.67 (0.48-0.95)0.0240.74 (0.61–0.90)0.002
miR-224-5p0.67 (0.50–0.89)0.0060.79 (0.51–1.22)0.280
miR-301a-3p0.83 (0.61–1.14)0.241
miR-3261.03 (1.01–1.05)0.0080.97 (0.80–1.18)0.749
miR-374b-5p0.85 (0.52–1,38)0.509
miR-494-3p1.00 (1.00–1.01)0.0231.01 (0.99–1.04)0.447
miR-939-5p1.02 (1.00–1.05)0.0210.98 (0.80–1.21)0.869
Abbreviations: HR, hazard ratio; CI, confidence interval. a Ninety-nine patients with biochemical recurrence and 108 patients without biochemical recurrence. b The 16 miRNAs with p-values < 0.200 in univariate Cox regression analysis were used in the Full model of multivariate Cox regression analysis, and the six miRNAs with p-values < 0.05 in this Full model were subsequently used in a Backward elimination approach.
Table
Table 5. Comparison of the C-statistics of prognostic indices for BCR prediction using two miRNA panels in Cox regression analyses.
Table 5. Comparison of the C-statistics of prognostic indices for BCR prediction using two miRNA panels in Cox regression analyses.
miRNA
Panel		All Samples
(n = 206)		Training Set
(n = 140)		Test Set
(n = 66)		p-value a
AUC (SE)AUC (SE)AUC (SE)
five-miR-panel 0.774 (0.040)0.735 (0.069)0.625
0.745 (0.034)
six-miR-panel 0.779 (0.040)0.741 (0.058)0.587
0.749 (0.034)
p-value b0.6230.9300.947
Abbreviations: AUC, area under the receiver operating characteristics curve; SE, standard error of the mean. a Values indicate that there are no significant differences obtained both with the five-miR-panel and six-miR-panel, using the training or the test set of samples. b Values indicate that there are no significant differences between the five-miR-panel and six-miR-panel using all samples, the training or test set of samples.
Table
Table 6. Construction of a predictive BCR classifier Cox regression in a bootstrapping approach with clinicopathological variables in the 206 patients.
Table 6. Construction of a predictive BCR classifier Cox regression in a bootstrapping approach with clinicopathological variables in the 206 patients.
Variable aUnivariate Cox RegressionMultivariate Cox Regression with Significant Variables
Full Model bReduced Model after Backward Elimination c
HR (95% CI)p-valueHR (95% CI)p-valueHR (95% CI)p-value
Age1.04 (1.00–1.08)0.0291.02 (0.98–1.07)0.231
PSA1.04 (1.01–1.07)0.0041.02 (0.99–1.05)0.207
DRE1.24 (0.83–1.95)0.286
Margin1.72 (1.16–2.56)0.0081.42 (0.64–3.11)0.548
pN status2.66 (0.88–4.10)0.1031.14 (0.74–1.75)0.396
pT stage1.12 (1.08–1.17)<0.00011.08 (1.03–1.14)0.0021.10 (1.04–1.15)0.001
ISUP Group1.37 (1.19–1.57)<0.00011.19 (1.02–1.40)0.0271.23 (1.06–1.43)0.007
a Abbreviations and stratifications of the variables as indicated in Table 1; HR, hazard ratio. b The full model included all variables of the univariate Cox regression with hazard ratios and p-values < 0.200. c Reduced model after backward elimination with entry p < 0.05 and removal p > 0.100.
Table
Table 7. Improved prediction of biochemical recurrence after radical prostatectomy, using clinicopathological-based tools in combination with the five-miR based panel.
Table 7. Improved prediction of biochemical recurrence after radical prostatectomy, using clinicopathological-based tools in combination with the five-miR based panel.
Prediction ToolClinicopathological-Based ToolClinicopathological-Based Tool Combined with the five-miR-Panelp-value
AUC (95% CI)AUC (95% CI)
Reference models
D’Amico et al. [15] 0.590 (0.519–0.657)0.759 (0.695–0.816)<0.0001
CAPRAS0.692 (0.624–0.754)0.769 (0.706–0.825)0.008
NCCN0.642 (0.572–0.707)0.757 (0.693–0.814)0.0005
Stephenson et al. [17]0.664 (0.595–0.728)0.747 (0.682–0.805)0.017
Present study
Full model0.723 (0.657–0.783)0.793 (0.731–0.8465)0.007
Reduced model0.712 (0.638–0.677)0.783 (0.720–0.837)0.007
Abbreviations: AUC, area under the receiver operating characteristics curve; CI, confidence interval; CAPRAS, Cancer of the Prostate Risk Assessment Postsurgical Score, calculated according to [16]; NCCN, National Comprehensive Cancer Network, calculated according to [18]; Full model, according to the Cox regression model described in Table 6, with all clinicopathological factors except for digital rectal examination; Reduced model, according to the Cox regression model described in Table 6 after backward elimination and finally including only the variables pT stage and ISUP Group grade.
Table
Table 8. Functional links between miRNAs of the five-miR-panel and prostate cancer.
Table 8. Functional links between miRNAs of the five-miR-panel and prostate cancer.
miRNA
[Reference]		Tumor Suppressor
Oncogene		Target GeneMolecular Mechanism
miR-30c-5p [79]Tumor suppressorASF/SF2Inhibition of tumor cell proliferation, promotion of apoptosis through the inhibition of ASF/SF2
[80]Tumor suppressorKRASInhibition of tumor cell proliferation, migration, and invasion
[51]Tumor suppressorBCL9Correlation of disease progression via BCL9/Wingles-type signaling
miR-31-5p [81]Tumor suppressorE2F6Promotion of apoptosis, reduced prostate cancer growth
[82]Tumor suppressorAndrogen receptorInhibition of cell proliferation and cell cycle
miR-141-3p [83]Tumor suppressorIGF1Rinhibit malignant transformation of benign prostate epithelial cell via IGF1R-AKT/STAT3 signaling pathway
[84]Tumor suppressorTRAF5, TRAF6Suppression of invasion and migration of PCa cells via inhibiting activation of NF-kB
[85]OncogeneKLF9Regulation of proliferation, spheroid formation, and expression of the stemness factors OCT-4, SOX9, and CCND1
miR-148a-3p [86]Tumor suppressorRTN4Inhibition of cell proliferation and cell cycle blocked in G2/M phase
[87]OncogeneCAND1Increased cell proliferation
[88]Tumor suppressorMSK1Inhibition of cell growth, migration, and invasion
miR-221-3p [89]Tumor suppressorRunx1, Runx2Inhibition of prostate tumorigenesis
[82]Tumor suppressorAndrogen receptor and its receptor coactivatorsInhibition of cell proliferation and targeted key oncogenic pathways, including cell cycle
[90]OncogenePI3K/AKT, P53Associated with signal transduction and cell communication
Abbreviations: ASF/SF2, serine and arginine rich splicing factor 1; KRAS, KRAS proto-oncogene, GTPase; BCL9, BCL9 transcription coactivator; E2F6, E2F transcription factor 6; IGF1R, insulin like growth factor 1 receptor; AKT/STAT3, serine-threonine kinase/ signal transducer and activator of transcription; TRAF5, TNF receptor associated factor 5; TRAF6, TNF receptor associated factor 6; NF-kB, nuclear factor kappa B subunit 1; KLF9, Kruppel like factor 9; OCT-4, organic cation/carnitine transporter4; SOX9, SRY-box transcription factor 9; CCND1, cyclin D1; RTN4, reticulon 4; CAND1, cullin associated and neddylation dissociated 1; MSK1, mitogen- and stress-activated kinase 1; Runx1, runt related transcription factor 1; Runx2, runt related transcription factor 1; P13K/AKT, phosphatidylinositol 3-kinase/serine-threonine kinase; P53, tumor protein 53.

Share and Cite

MDPI and ACS Style

Zhao, Z.; Weickmann, S.; Jung, M.; Lein, M.; Kilic, E.; Stephan, C.; Erbersdobler, A.; Fendler, A.; Jung, K. A Novel Predictor Tool of Biochemical Recurrence after Radical Prostatectomy Based on a Five-MicroRNA Tissue Signature. Cancers 2019, 11, 1603. https://0-doi-org.brum.beds.ac.uk/10.3390/cancers11101603

AMA Style

Zhao Z, Weickmann S, Jung M, Lein M, Kilic E, Stephan C, Erbersdobler A, Fendler A, Jung K. A Novel Predictor Tool of Biochemical Recurrence after Radical Prostatectomy Based on a Five-MicroRNA Tissue Signature. Cancers. 2019; 11(10):1603. https://0-doi-org.brum.beds.ac.uk/10.3390/cancers11101603

Chicago/Turabian Style

Zhao, Zhongwei, Sabine Weickmann, Monika Jung, Michael Lein, Ergin Kilic, Carsten Stephan, Andreas Erbersdobler, Annika Fendler, and Klaus Jung. 2019. "A Novel Predictor Tool of Biochemical Recurrence after Radical Prostatectomy Based on a Five-MicroRNA Tissue Signature" Cancers 11, no. 10: 1603. https://0-doi-org.brum.beds.ac.uk/10.3390/cancers11101603

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop