Utilization of Biopsy-Based Genomic Classifier to Predict Distant Metastasis after Definitive Radiation and Short-Course ADT - Full Text Article

BACKGROUND: We examined the ability of a biopsy-based 22-marker genomic classifier (GC) to predict for distant metastases after radiation and a median of 6 months of androgen deprivation therapy (ADT). 

METHODS: We studied 100 patients with intermediate-risk (55%) and high-risk (45%) prostate cancer who received definitive radiation plus a median of 6 months of ADT (range 3–39 months) from 2001–2013 at a single center and had available biopsy tissue. Six to ten 4 micron sections of the needle biopsy core with the highest Gleason score and percentage of tumor involvement were macro dissected for RNA extraction. GC scores (range, 0.04–0.92) were determined. The primary end point of the study was time to distant metastasis. Median follow-up was 5.1 years. There were 18 metastases during the study period. 

RESULTS: On univariable analysis (UVA), each 0.1 unit increase in GC score was significantly associated with time to distant metastasis (hazard ratio: 1.40 (1.10–1.84), P = 0.006) and remained significant after adjusting for clinical variables on multivariable analysis (MVA) (adjusted hazard ratio: 1.36 (1.04–1.83), P = 0.024). The c-index for 5-year distant metastasis was 0.45 (95% confidence interval: 0.27–0.64) for Cancer of the Prostate Risk Assessment score, 0.63 (0.40–0.78) for National Comprehensive Cancer Network (NCCN) risk groups, and 0.76 (0.57–0.89) for the GC score. Using pre-specified GC risk categories, the cumulative incidence of metastasis for GC40.6 reached 20% at 5 years after radiation (P = 0.02). 

CONCLUSIONS: We believe this is the first demonstration of the ability of the biopsy-based GC score to predict for distant metastases after definitive radiation and ADT for intermediate- and high-risk prostate cancer. Patients with the highest GC risk (GC40.6) had high rates of metastasis despite multi-modal therapy suggesting that they could potentially be candidates for treatment intensification and/or enrollment in clinical trials of novel therapy.

Prostate Cancer and Prostatic Diseases (2017) 20, 186–192; doi:10.1038/pcan.2016.58; published online 24 January 2017

INTRODUCTION: Radiation and androgen deprivation therapy (ADT) is a standard therapy for contemporary patients with intermediate and high-risk prostate cancer.1 While many men will be cured with this treatment, there remains a proportion of men who will progress after therapy and develop metastatic disease. For these men, intensification of therapy beyond standard duration ADT and radiation may be required to further reduce the risk of metastasis.

Because intensification of therapy using longer durations of ADT, second-generation anti-androgens, chemotherapy or novel agents carries risk of additional toxicity, it is critically important to correctly identify the subgroup of patients who may be in need of such intensification. While clinical factors such as PSA, T-category and Gleason score have traditionally been used to risk-stratify patients, their accuracy has been significantly enhanced in the last few years by genomic-based tests. However, most genomic tests for prostate cancer have been discovered or validated in either the post-prostatectomy setting or more recently the active surveillance setting, and it remains unknown whether they can predict for a clinically meaningful end point such as distant metastasis after definitive radiation and ADT.

For example, the Oncotype Dx test from Genomic Health is used on biopsy tissue from National Comprehensive Cancer Network (NCCN) very low through intermediate risk prostate cancer and gives the risk of harboring adverse pathology (pT3 or Gleason ⩾ 4 +3) at prostatectomy.2 The Prolaris test from Myriad Genetics can predict 10-year prostate cancer specific mortality with conservative management from biopsy tissue or biochemical recurrence from radical prostatectomy tissue.3 Finally, the Decipher test from GenomeDx Biosciences can be used on both biopsy and prostatectomy specimens to predict the risk of distant metastasis after prostatectomy.4–6

In this study, we examined whether the biopsy-based 22-gene Decipher genomic classifier (GC) can be used to accurately predict for the risk of distant metastasis in a contemporary cohort with both intermediate and high-risk prostate cancer who all received radiation and ADT.

MATERIALS AND METHODS - Patient selection and treatment: Our initial cohort included 153 patients who received radiation and ADT at the Dana-Farber/Brigham and Women’s Cancer Center from 2001–2013 for NCCN intermediate or high-risk prostate cancer and had archived formalin-fixed paraffin-embedded tissue available. Twenty-six patients did not have adequate tumor tissue for RNA extraction while 27 patients did not pass prespecified microarray quality control thresholds and were thus excluded from analysis. GC scores were calculated for the remaining 100 patients, consisting of 55 intermediate risk and 45 high-risk patients. Of note, the clinical characteristics between those with adequate vs inadequate tissue for analysis was similar except for a higher proportion of positive biopsy cores in those with adequate tissue (median 50% vs 33%, Po0.001; Supplementary Table 1). The duration of ADT was at the discretion of the treating physician, and most (68%) received 6 months of ADT with a range of 3–39 months. The ADT was typically started 2 months before the initiation of radiation. Radiation was external beam alone for 97% of the patients to a median dose of 72 Gy (range: 68.4–81.0 Gy). Included patients had all consented to a prospective tissue collection protocol, and this study was approved by the Dana-Farber/Harvard Cancer Center Institutional Review Board.

Specimen collection and processing Six to ten 4 micron sections of the needle biopsy core with the highest Gleason score and percentage of tumor involvement were macrodissected for RNA extraction. GC scores were determined from the Decipher prostate cancer classifier assay (GenomeDx Biosciences Laboratory, San Diego, CA, USA) as previously described.4,5 Cancer of the Prostate Risk Assessment (CAPRA) was calculated as previously described5 while substituting pre-operative PSA with pre-radiation PSA. Risk categorization of GC and CAPRA were based on prior publications.5,7,8 The primary end point of the study was distant metastasis following radiation therapy. Secondary end points evaluated were biochemical failure (nadir plus 2 definition) and castrate resistance (defined as any rise in PSA despite being on salvage ADT for biochemical failure).

Statistical analysis Duration of ADT (⩽6 months vs 46 months), biopsy Gleason (⩽3+4 vs ⩾4+3) and clinical stage (oT2b vs ⩾T2b) were treated as categorical variables while age at radiation therapy (RT), pretreatment PSA (log 2 transformed), and percent positive cores were modeled as continuous covariates. In time to event analyses, event times were defined as the time from initiation of RT to metastases or date of last follow-up.

The performance of GC, CAPRA and NCCN risk categories were compared and contrasted by measuring their ability to (1) independently predict metastases using univariable (UVA) and multivariable (MVA) penalized Cox regression;9 (2) discriminate risk among patients using survival receiver operating characteristic curves and their respective c-indices;10 (3) exhibit clinical usefulness by plotting out decision curves adapted to survival data;11 and (4) stratify metastatic risk using cumulative incidence curves12 based on their previously published cutpoints. Least absolute shrinkage and selection operator (LASSO) regression was a secondary penalized method used to assess the relative importance of the variables of interest in predicting metastasis. Furthermore, it has been demonstrated that LASSO regression without penalization on the exposure variable produces an estimate of the coefficient of the exposure variable of minimal bias with as few as three events per variable.13 This method was also applied to the data. Confidence intervals for survival c-indices were computed via the bootstrap. The c-indices were considered statistically significant if the lower bound of the 95% confidence interval exceeded 0.50. The significance level was set at 0.05 for all tests while analyses were performed in R v3.0. With an anticipated event rate of 0.13, this study had 80% power to detect a hazard ratio of 1.54 per 0.1 unit increase in continuous genomic biomarker score.

RESULTS - Patient Characteristics
Patient characteristics are provided in Table 1. The median age of patients in the study was 67 years (range, 45–87) and 16% of patients were African–American. Thirty percent of patients had a biopsy Gleason ⩽3+4, 36% had Gleason 4+3, and 34% of patients had Gleason 8 or 9 disease. The median pretreatment PSA was 7.3 ng ml − 1 (interquartile range: 4.7–-14.9 ng ml − 1 ) while the median time between biopsy and radiation therapy was 4.5 months. Sixty-eight percent of patients received 6 months of ADT, 1% received 3 months of ADT, and the remaining 31% received 12–39 months of ADT. Eighty-seven percent received both a gonadotropin-releasing hormone agonist and an anti-androgen. The median follow-up on censored patients was 5.1 years (interquartile range: 3.4-6.3) and 18 patients developed metastasis during study follow-up. During this same period, 28 patients had biochemical failure while 12 developed castrate resistant disease.

Demographic and clinical characteristics of eligible patients

Calculation of GC and CAPRA scores The distributions of GC and CAPRA scores are presented in Figure 1. The median GC score for this biopsy-based cohort was 0.39 (interquartile range: 0.22–0.61). The median CAPRA score for these patients was 5 (interquartile range: 4–6). We observed a significant positive correlation between GC and biopsy Gleason score (r = 0.38, Po0.001; Supplementary Figure 1) as well as CAPRA scores (r = 0.37, Po0.001).

CAPRA b genomic classifier risk scores
Figure 1. Distributions of the study cohort by (a) Cancer of the Prostate Risk Assessment (CAPRA), (b) genomic classifier risk scores.

Results of cox proportional hazards


Ability of GC to predict for metastases and secondary end points On UVA of baseline clinical and genomic risk factors we found only GC was a significant predictor of metastasis, associated with a 40% increase in the risk of metastasis per 0.1 unit increase in score (HR: 1.40; 95% confidence interval (CI): 1.10–1.84; P = 0.006; Table 2). The CAPRA score had a hazard ratio of 1.15 per unit increase but was not statistically significant (P = 0.271). Likewise, NCCN high-risk (ref: intermediate risk) was not significantly predictive of metastasis (HR: 2.00; 95% CI: 0.78-5.35; P = 0.147).

1 Validating biopsy decipher post radiation and ADT
Figure 2. (a) Least absolute shrinkage and selection operator (LASSO) coefficient path demonstrating the order of importance of genomic classifier (GC) and clinical variables in predicting metastasis. Moving from right to left, the order of nonzero hazards coefficients represents the order of variable importance. (b) LASSO coefficient path without penalization on GC using a cross-validated penalty parameter of 0.049, represented by a vertical dashed line. Only GC, percent of positive cores and clinical stage have nonzero hazards coefficients at this level of penalization. This model estimates a less biased hazard ratio for GC when the number of events per variable is low. (c) Survival c-indices at 5 years following radiation therapy (RT) for GC, Cancer of the Prostate Risk Assessment (CAPRA) and National Comprehensive Cancer Network (NCCN) risk. (d) Decision curve analysis comparing net benefit at 5 years post-RT of GC and CAPRA across various threshold probabilities. Compared with ‘treat none’ and ‘treat all’ scenarios (in which no risk prediction model is employed) to make treatment decisions, across a range of threshold probabilities GC had the highest net benefit compared with the clinical-only CAPRA risk model. The net benefit is defined as a measure of the relative value of benefits from identifying higher risk men that should for example, receive more intensive therapy (for example, longer duration hormonal suppression) and harms (for example, morbidity of long-term androgen deprivation therapy (ADT)) associated with the GC and CAPRA risk models.

2 Validating biopsy decipher post radiation and ADT
Figure 3. (a) Cumulative incidence curves in which patients are stratified by National Comprehensive Cancer Network (NCCN) risk categories. (b) Cumulative incidence curves in which patients are stratified by Cancer of the Prostate Risk Assessment (CAPRA) risk categories. (c) Cumulative incidence curves in which patients are stratified by genomic classifier (GC) risk categories. (d) Cumulative incidence curves in which patients are stratified by GC risk using an exploratory cutoff of 0.2. RT, radiation therapy.

When adjusting for relevant clinical variables in MVA including stage, Gleason, PSA, percent of positive cores, duration of ADT and year of treatment, we observed a small reduction in the GC hazard ratio (HR: 1.36; 95% CI: 1.04–1.83) but it remained a significant predictor of metastasis (P = 0.024). LASSO regression established

that GC was the most important variable in predicting metastasis as it was the first variable to enter the model followed by percent of positive biopsy cores with age at RT being the last and, therefore, least important variable (Table 3; Figure 2a). LASSO regression without penalizing GC, but with a cross-validated penalty parameter of 0.049 on all other variables, produced a hazard ratio for GC of 1.44 (Table 3; Figure 2b), closely corroborating the results from MVA.

In a separate MVA that included both GC and CAPRA, as well as one that included GC and NCCN risk, GC’s independent predictive capability was consistent across models with HR’s of 1.44 (p = 0.012) and 1.37 (p = 0.014), respectively. Neither CAPRA nor NCCN were significant predictors in these models and each saw a large variance in their hazard ratios when compared with their UVA results. Similar results were observed when GC was modeled with biopsy Gleason score (Table 2).

Cox regression analysis on secondary end points demonstrated that only CAPRA was a significant predictor of biochemical failure on UVA (HR: 1.23; 95% CI: 1.01-1.49; P = 0.042; Supplementary Table 2) but neither of the models were significant predictors on MVA. With respect to the development of castrate resistant disease, GC was a significant predictor on UVA as well as after adjusting for NCCN risk (HR for GC: 1.43; 95% CI: 1.01-2.09; P = 0.044) and for CAPRA (HR for GC: 1.48; 95% CI: 1.00–2.45; P = 0.049; Supplementary Table 3).

We next evaluated the discrimination performance of GC, which had a survival c-index at 5 years following RT for predicting metastasis of 0.76 (95% CI: 0.57–0.89; Figure 2c). This compared favorably to a c-index of 0.45 for CAPRA (95% CI: 0.27–0.64) and 0.63 for NCCN (95% CI: 0.40–0.78), as only the GC’s c-index confidence interval did not contain 0.5. The c-index showed consistent discrimination by GC over time, as the c-index of GC for predicting metastasis at 10 years following RT was 0.78 (95% CI: 0.60–0.87; Supplemental Figure 2). Decision curve analysis also demonstrated that the net benefit of using GC for treatment decisions was generally higher than basing clinical decisions on either the CAPRA risk model or naively choosing to either treat all patients or to treat none (Figure 2d). Therefore, irrespective of the threshold probability for making a treatment decision based on the GC risk model, decision curve analysis shows that using GC will improve decision-making.

Metastatic risk stratification When studying risk stratification, there was not a significant difference in the cumulative incidence of metastasis between NCCN high and intermediate patients (P = 0.238; Figure 3). CAPRA appeared to stratify risk better but the differences were not statistically significant (P = 0.215). However, pre-determined GC risk groups showed a significant ability to stratify patient risk of metastasis, with the 26% of patients in the high genomic risk group (GC40.6) experiencing a 5-year cumulative risk of metastasis of 20% after radiation and ADT (P = 0.02). The low and intermediate genomic risk group curves overlapped beyond 5 years as the number of patients at risk within each group dropped substantially thereafter. In exploratory analysis, we found that patients with a GC score ⩽0.2 had a 0% cumulative incidence of metastasis throughout the period of follow-up (P = 0.07).

DISCUSSION: To the best of our knowledge, this study of contemporary intermediate and high-risk patients treated with radiation and a median of 6 months of ADT is the first to demonstrate that a biopsy-based 22-gene GC can predict for distant metastases after definitive radiation and ADT. We found that the GC outperformed clinical variables and accurately predicted the 5-year risk of distant metastases with a c-Index of 0.76. In addition, while this study contained mostly men with intermediate risk disease, the 26% with a high GC score (40.6) had a 5-year risk of distant metastasis of 20% despite treatment with both radiation and ADT.

The clinical implication of this study is that the GC, which was developed on men treated with prostatectomy, may be used to accurately risk-stratify men for metastatic failure after radiation and ADT. For the 26 percent of men with a high GC score, there was a high-risk of metastases at 5 years despite radiation and ADT, and these men may need to consider treatment intensification, such as with longer duration ADT which has been shown to improve survival for locally advanced disease,14,15 or with docetaxel which has shown early promise in the initial report of RTOG 05-21 (ref. 16), or with entry into a trial of novel therapies for aggressive prostate cancer. Dose was not predictive of metastasis in this cohort, and whether further radiation dose intensification could benefit men with a high GC score remains unknown. Conversely, for the 28% of men who had a GC score ⩽ 0.2, the risk of metastases appeared to be very low. While this cutpoint was exploratory and will need validation before being considered for clinical use, it raises the possibility that sometime in the future, given the lack of development of metastases in these men, and the increasing recognition of the harms of ADT,17,18 men with a GC score ⩽0.2 may be reasonable candidates for forgoing ADT and choosing dose-escalated radiation alone, which is currently an NCCN-endorsed option for some men with intermediate and even select men with high-risk disease.

Although further studies in larger data sets will be needed to validate the risk-adapted treatment paradigm described above, this study highlights the significant value that GCs may have for patients who choose radiation as primary management of their prostate cancer. The GC outperformed clinical variables in predicting for distant metastasis, and this may reflect the fact that the classifier is able to measure multiple pathways at once, including information about androgen signaling, cell-cycle progression and chromosome structure maintenance.19 Another advantage of the GCs is that they are reproducible lab-based assays that do not depend on institutional expertize, as opposed to Gleason grading, which can be highly variable depending on the expertize of the reading pathologist.20,21 Finally, recent data suggests that GCs may also have the advantage of being less sensitive to sampling error.2

Work by Freedland et al.22 has also highlighted the potential value of genomics in radiation-managed men. In their study of 147 radiation-managed men of whom about half received ADT and half did not, a 33 gene classifier based on cell-cycle progression (Myriad) was associated with the risk of biochemical failure after radiation. It was also associated with the risk of prostate-cancer mortality, although the number of events was small at 6. Our study, which uses a different GC, adds to this literature by focusing exclusively on men treated with both radiation and ADT, and emphasizing the highly clinically significant end point of distant metastasis. This has allowed us to identify a subgroup representing 26% of the men who remain at high-risk of distant metastasis despite treatment both with radiation and ADT and, therefore, need consideration of even further treatment intensification.

A limitation of the current study is the size of the data set and length of follow-up, and so further studies will be needed to validate these findings in larger data sets with longer follow-up. In addition, the hypothesis that patients with very low GC score ⩽ 0.2 may be able to omit ADT requires further testing in previously treated cohorts as well as in prospective studies, which are currently being planned. Finally, while the GC score was prognostic for distant metastasis, it did not have a significant association with biochemical recurrence. This may reflect the fact that the test was originally developed to specifically predict for distant metastasis and generally only a minority of biochemical recurrences will lead to distant metastasis. Another consideration is that none of the patients in this study received a multiparametric magnetic resonance imaging, and it has previously been shown that magnetic resonance imaging can add some prognostic information to clinical variables through upstaging or through detecting potentially higher grade lesions.23 Future studies should evaluate how GC adds to prognostic value when multiparametric magnetic resonance imaging has also been performed,24 although studies with the cell-cycle progression score suggest that the magnetic resonance imaging and the genomic information are capturing different types of prognostic information.25

CONCLUSION: In summary, this is the first report of the ability of a 22-gene biopsy-based GC to accurately predict distant metastasis in men with intermediate or high-risk prostate cancer treated with radiation and a median of 6 months of ADT. Patients with the highest GC risk (GC40.6) had high rates of metastasis despite multi-modal therapy and could potentially be candidates for intensification of therapy and/or entry into clinical trials of novel therapy.

CONFLICT OF INTEREST: PLN has consulted for Medivation, GenomeDx and Ferring. FYF has consulted for Medivation and GenomeDx. NEM has received personal fees from Via Oncology. VC, BP-A, TK, HS and ED are all employees of GenomeDx Biosciences. The remaining authors declare no conflicts of interest.

ACKNOWLEDGEMENTS: This work was funded by a grant from an Anonymous Family Foundation, the Prostate Cancer Foundation, Fitz’s Cancer Warriors, Cynthia and David Chapin, Hugh Simons in Honor of Frank and Anne Simons, The Gina Ventre and Scotty Forbes Fund, The Campbell Family in Honor of Joan Campbell and GenomeDx Biosciences.

Read More: A Commentary from the Associate Editor of PCAN - Andrew Armstrong, MD, ScM, FACP

Written By: 1, NE Martin1, V Choeurng2, B Palmer-Aronsten2, T Kolisnik2, CJ Beard1, PF Orio1, MD Nezolosky1, Y-W Chen1, H Shin2, E Davicioni2 and FY Feng3

1. Department of Radiation Oncology, Dana-Farber/Brigham and Women’s Cancer Center and Harvard Medical School, Boston, MA, USA;
2. GenomeDx Biosciences, Vancouver, BC, Canada and 3 Department of Radiation Oncology, University of California at San Francisco, San Francisco, CA, USA.
Correspondence: Dr PL Nguyen, Department of Radiation Oncology, Dana-Farber/Brigham and Women’s Cancer Center and Harvard Medical School, Boston, MA, USA


1. Mohler JL, Armstrong AJ, Bahnson RR, D’Amico AV, Davis BJ, Eastham JA et al. Prostate cancer, version 1. 2016. J Natl Compr Canc Netw 2016; 14: 19–30. 

2. Klein EA, Cooperberg MR, Magi-Galluzzi C, Simko JP, Falzarano SM, Maddala T 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. 

3. Cuzick J, Swanson GP, Fisher G, Brothman AR, Berney DM, Reid JE et al. Prognostic value of an RNA expression signature derived from cell cycle proliferation genes in patients with prostate cancer: a retrospective study. Lancet Oncol 2011; 12: 245–255. 

4. Karnes RJ, Bergstralh EJ, Davicioni E, Ghadessi M, Buerki C, Mitra AP et al. Validation of a genomic classifier that predicts metastasis following radical prostatectomy in an at risk patient population. J Urol 2013; 190: 2047–2053. 

5. Klein EA, Haddad Z, Yousefi K, Lam LLC, Wang Q, Choeurng V et al. Decipher genomic classifier measured on prostate biopsy predicts metastasis risk. Urology 2016; 90: 148–152. 

6. Knudsen BS, Kim HL, Erho N, Shin H, Alshalalfa M, Lam LLC et al. Application of a clinical whole-transcriptome assay for staging and prognosis of prostate cancer diagnosed in needle core biopsy specimens. J Mol Diagn 2016; 18: 395–406. 

7. Cooperberg MR, Broering JM, Carroll PR. Risk assessment for prostate cancer metastasis and mortality at the time of diagnosis. J Natl Cancer Inst 2009; 101: 878–887. 

8. Ross AE, Yousefi K, Davicioni E, Ghadessi M, Johnson MH, Sundi D et al. Utility of risk models in decision making after radical prostatectomy: lessons from a natural history cohort of intermediate- and high-risk men. Eur Urol 2015; 69: 496–504. 

9. Firth D. Bias reduction of maximum likelihood estimates. Biometrika 1993; 80: 27. 

10 Heagerty PJ, Lumley T, Pepe MS. Time-dependent ROC curves for censored survival data and a diagnostic marker. Biometrics 2000; 56: 337–344. 

11. Vickers AJ, Cronin AM, Elkin EB, Gonen M. Extensions to decision curve analysis, a novel method for evaluating diagnostic tests, prediction models and molecular markers. BMC Med Inform Decis Mak 2008; 8: 53. 

12. Fine JP, Gray RJ. A proportional hazards model for the subdistribution of a competing risk. J Am Stat Assoc 1999; 94: 496. 

13. Chen Q, Nian H, Zhu Y, Talbot HK, Griffin MR, Harrell FE. Too many covariates and too few cases? - a comparative study. Stat Med 2016; 35: 4546–4558. 

14. Horwitz EM, Bae K, Hanks GE, Porter A, Grignon DJ, Brereton HD et al. Ten-year follow-up of radiation therapy oncology group protocol 92-02: a phase III trial of the duration of elective androgen deprivation in locally advanced prostate cancer. J Clin Oncol 2008; 26: 2497–2504. 

15. Bolla M, de Reijke TM, Van Tienhoven G, Van den Bergh ACM, Oddens J, Poortmans PMP et al. Duration of androgen suppression in the treatment of prostate cancer. N Engl J Med 2009; 360: 2516–2527.

16. Sandler H, Hu C, Rosenthal S, Sartor O, Gomella LG, Amin MB et al. A phase III protocol of androgen suppression (AS) and 3DCRT/IMRT versus AS and 3DCRT/ IMRT followed by chemotherapy (CT) with docetaxel and prednisone for localized, high-risk prostate cancer (RTOG 0521). In: ASCO Genitourinary Cancer Symposium. Chicago, IL, 2015. http://abstracts.asco.org/156/AbstView_ 156_144388.html. 

17. Nguyen PL, Alibhai SMH, Basaria S, D’Amico AV, Kantoff PW, Keating NL et al. Adverse effects of androgen deprivation therapy and strategies to mitigate them. Eur Urol 2014; 67: 1–12. 

18. Dinh KT, Reznor G, Muralidhar V, Mahal BA, Nezolosky MD, Choueiri TK et al. Association of androgen deprivation therapy with depression in localized prostate cancer. J Clin Oncol 2016; 67: 825–836. 

19. Alshalalfa M, Schliekelman M, Shin H, Erho N, Davicioni E. Evolving transcriptomic fingerprint based on genome-wide data as prognostic tools in prostate cancer. Biol Cell 2015; 107: 232–244. 

20. Burchardt M, Engers R, Müller M, Burchardt T, Willers R, Epstein JI et al. Interobserver reproducibility of Gleason grading: evaluation using prostate cancer tissue microarrays. J Cancer Res Clin Oncol 2008; 134: 1071–1078. 

21. Nguyen PL, Schultz D, Renshaw AA, Vollmer RT, Welch WR, Cote K et al. The impact of pathology review on treatment recommendations for patients with adenocarcinoma of the prostate. Urol Oncol 22: 295–299. 

22. Freedland SJ, Gerber L, Reid J, Welbourn W, Tikishvili E, Park J et al. Prognostic utility of cell cycle progression score in men with prostate cancer after primary external beam radiation therapy. Int J Radiat Oncol Biol Phys 2013; 86: 848–853. 

23. Ho R, Siddiqui MM, George AK, Frye T, Kilchevsky A, Fascelli M et al. Preoperative multiparametric magnetic resonance imaging predicts biochemical recurrence in prostate cancer after radical prostatectomy. PLoS One 2016; 11: e0157313. 

24. Porpiglia F, Cantiello F, De Luca S, Manfredi M, Veltri A, Russo F et al. In-parallel comparative evaluation between multiparametric magnetic resonance imaging, prostate cancer antigen 3 and the prostate health index in predicting pathologically confirmed significant prostate cancer in men eligible for active surveillance. BJU Int 2015; 118: 527–534. 

25. Renard-Penna R, Cancel-Tassin G, Comperat E, Varinot J, Léon P, Roupret M et al. Multiparametric magnetic resonance imaging predicts postoperative pathology but misses aggressive prostate cancers as assessed by cell cycle progression score. J Urol 2015; 194: 1617–1623.