The management of nonmetastatic prostate cancer is complex because of the interplay of multiple considerations, including risk stratification, relative treatment efficacy/toxicity, competing risk of death from cancer versus other causes, and patient preferences.1 Published randomized controlled trials support the use of radiotherapy in the primary2 and postoperative3 settings, with other randomized controlled trials assessing appropriate dosing of external-beam radiotherapy4 and the concomitant use of hormonal therapy.5 Alternatives to dose-escalated external-beam radiotherapy include low-dose-rate (LDR)6 and high-dose-rate (HDR)7 brachytherapy (with or without integrated external-beam radiation therapy), hypofractionated radiation therapy,8 stereotactic body radiation therapy,9 and particle therapy.10 Various clinical trials have been initiated to further define best practices for the use of radiotherapy either alone or in combination with other treatments.
Multiple pretreatment risk stratification systems have been published to support decision-making, interphysician communication, clinical trial stratification/design, and outcome reporting in nonmetastatic prostate cancer. Five unique commonly used systems (not exhaustive) have been described in the literature, including
Harvard11/American Urological Association12/European Association of Urology13
Genitourinary Radiation Oncologists of Cancer (GUROC)14/National Institute for Health and Clinical Excellence15
Cancer of the Prostate Strategic Research Endeavour (CaPSURE)16
European Society of Medical Oncology17
Despite the proliferation of publications on prostate cancer prognostic factors and risk stratification systems, the definition of a universally accepted classification system with superior statistical and clinical decision-making properties is still required.20 In a recent review of pretreatment prostate cancer risk stratification, 4 specific improvements to risk stratification systems were identified, including 1) the definition of a very-low-risk category that would be appropriate for surveillance strategies, 2) reassessment of the interface between intermediate risk and high risk, 3) the definition of a new extreme-risk category, and 4) division of the intermediate-risk category into 2 subgroups.21 The latest NCCN risk stratification system has addressed 3 or 4 of these potential improvements (ie, very-low-risk category, change of the intermediate- and high-risk interface, and the definition of a very-high-risk category) with the creation of a 5-group risk stratification system.18,19 The ability of this new system to categorize patients into distinct clinical outcomes group has not been previously assessed; however, the previous 3-group NCCN system compared favorably with the American Joint Committee on Cancer Staging system in a recent publication.22
This investigation reports on the creation and clinical outcomes related to a pan-Canadian Prostate Cancer Risk Stratification (ProCaRS) radiotherapy database, which was constructed to explore various risk stratification issues in prostate cancer. This database will also explore the statistical properties of the 5-group NCCN risk stratification system for the prediction of clinical prostate cancer outcomes.
Materials and Methods
The GUROC ProCaRS Database
A GUROC ProCaRS database template based on the British Columbia Cancer Agency (BCCA) cancer database23 was created and approved for dissemination to all Canadian investigators with potential databases for inclusion in this national database project. Primary ethics approval was obtained from the Western University Research Ethics Board (London, ON, Canada), with secondary ethics approvals obtained at all participating sites before data exchange. ProCaRS template-compliant data on 7974 patients who underwent radiotherapy were received from 4 institutions (BCCA; n=3771 [47%]), Princess Margaret Hospital (PMH; n=1752 [22%]), McGill (n=194 [3%]), and Laval (n=2257 [28%]), with a total of 7 unique databases, and locked for analysis on October 15, 2011. These 7 databases included a provincial database (BCCA, n=3771 [47%]), 2 LDR brachytherapy databases (PMH, n=1163 [15%] and Laval, n=1630 [20%]), 1 combined HDR plus external-beam radiotherapy database (Laval, n=627 [8%]), and 3 external-beam databases (PMH dose-escalation, n=401 [5%]; PMH neoadjuvant Casodex study #9907, n=188 [2.5%]; and McGill radiotherapy, n=194 [2.5%]).
Database Quality Assurance and Clinical Outcomes
An extensive quality assurance program was initiated to confirm database quality, including various error checking, missing data minimization, and outcome confirmation procedures. Risk stratification categories were calculated from the available information in the ProCaRS database, including GUROC 3-group14 and NCCN 5-group18,19 risk stratification. In terms of GUROC risk stratification, the following risk groups have been defined:
Low-risk: 1997 AJCC T1-T2a, prostate-specific antigen (PSA) ≤10 ng/mL, and Gleason score ≤6
Intermediate-risk: 1997 AJCC T1-T2, PSA ≤20 ng/mL, and Gleason score ≤7 not otherwise low-risk
High-risk: 1997 AJCC T3-T4 or PSA >20 ng/mL, or Gleason score 8-10
The current NCCN risk categories18,19 are defined as:
Very-low-risk: T1c, Gleason score ≤6, PSA <10, <3 biopsy cores positive, ≤50% cancer in each core, and PSA density <0.15 ng/mL/g
Low-risk: T1-T2a, Gleason score 2-6, and PSA <10 not very-low-risk
Intermediate-risk: T2b or T2c, and/or Gleason score 7, and/or PSA 10-20 not low-risk
High-risk: T3a, or PSA >20, or Gleason score 8-10 not very high risk
Very-high-risk: T3b-4
Various clinical outcomes were computed from the ProCaRS database using the date of radiotherapy initiation as the starting date for all outcome calculations. The primary outcome for the analysis was defined as the biochemical failure-free survival (BFFS; biochemical failure or initiation of salvage therapy before biochemical failure) as currently defined by the American Society of Radiation Oncology (ASTRO).24 The secondary end point for this analysis was overall survival. Other end points that have been computed in the ProCaRS database include prostate cancer-specific survival and ASTRO biochemical failure (ie, biochemical failure alone not including initiation of salvage therapy). Additionally, biochemical failure and failure-free survival using the original ASTRO definition (ie, 3 successive increases in PSA, with failure backdated to halfway between 0th and 1st rising PSA) are also available in the database.25 Technical biochemical failures from brachytherapy bounces were adjusted using a quality assurance procedure whereby patients with PSA levels that decreased to 0.5 ng/mL or less without intervention were considered to not have had a biochemical failure.26
Statistical Analyses
Descriptive analysis of all available patient, tumor, and treatment variables was performed. Additionally, a descriptive analysis of various risk stratification categories and available outcomes was also performed. Both univariable and multivariable Cox proportional hazards modeling were performed on the primary (ASTRO BFFS) and secondary (prostate cancer-specific survival, overall survival) outcomes of interest. Kaplan-Meier curves were generated for both the primary and secondary outcomes stratified by GUROC/NCCN risk stratification for all patients, all external-beam-alone patients, all external-beam-alone patients receiving greater than the equivalent of 70 Gy in 35 fractions of radiotherapy in 1.8 to 2.0 Gy/d (70 Gy/35#: 185/904 [20.5%]; 72 Gy/36#: 15/904 [1.7%]; 74 Gy/37#: 211/904 [23.3%]; 75.6 Gy/42#: 82/904 [9.1%]; 76 Gy/38#: 3/904 [0.3%]; 78 Gy/39#: 3/904 [0.3%]; and 79.8 Gy/42#: 405/904 [44.8%]), and all patients treated with brachytherapy (LDR, HDR brachytherapy with or without external-beam). Statistical significance was determined using the log-rank statistic. All statistical analyses were performed using either SAS 9.2 software (SAS institute, Cary NC) or the R software platform (www.r-project.org).
Results
Descriptive Analysis
Mean age was 66.5 years (SD, 7.4). Mean PSA, PSA velocity, and PSA doubling time were 9.19 ng/mL (SD, 11.0), 0.27 ng/mL/y, and 1.17 years, respectively. T-stage distribution was 45.2% for T1, 46.0% for T2, 8.2% for T3, and 0.6% for T4. Gleason grade distribution was 65% for 2 to 6, 30% for 7, and 5% for 8 to 10. Mean percentage core involvement was 43%. GUROC risk stratification was low in 3928 (49.5%), intermediate in 2888 (36.5%), and high in 1097 (14%) patients. NCCN risk stratification was very low in 877 (11%), low in 2776 (35%), intermediate in 3100 (39.5%), high in 800 (10%), and very high in 293 (3.5%) patients. Descriptives broken down according to participating institutional databases are depicted in Table 1.
External-beam radiation therapy was used as monotherapy in 2793 of patients with a median dose fractionation of 68 Gy in 34 fractions (range, 50.0 Gy in 20 fractions to 79.8 Gy in 42 fractions). LDR brachytherapy was used in 4560 (57%) patients with a mean reported dose of 153.9 Gy (SD, 13.2). HDR brachytherapy (in combination with 40-45 Gy in 1.8-2.0 Gy/d external-beam radiation therapy) was used in 737 (9%) patients with a mean reported dose of 17 Gy (SD, 4.0). Collectively, 5297 (66%) patients received some form of brachytherapy as primary treatment of their prostate cancer. Hormonal therapy was used in 2999 patients (38%), with a mean use of 10.4 months. Table 2 depicts hormonal therapy use and duration versus GUROC/NCCN risk stratification grouping.
Table 3 summarizes clinical outcomes for the entire ProCaRS database and its constituent databases. Median follow-up for the entire ProCaRS database was 78.9 months (range, 0-190.7). Median follow-up for the 7 included databases ranged from 28.1 to 94.1 months. A total of 1442 of 7974 patients (19.1%) had ASTRO BFFS, with 1392 (18.5%) having a biochemical failure before any initiation of salvage therapy. A total of 1230 (15.4%) patients died, with 273 of these mortality events (22.0%) confirmed as cancer-related deaths. Additionally, 821 of 1230 patients (67.0%) died of other causes, with another 136 (11.0%) dying without any mortality attribution.
Cox Regression Analyses
Multivariable modeling for both the primary (ASTRO BFFS) and secondary (overall survival) end points showed that age, PSA, T stage, Gleason score grouping, hormone therapy duration, and type of radiation treatment were predictive of both of these clinical outcomes (P<.05; Table 4). Presence or absence of hormonal therapy was found to be significant only on univariable analysis (and not multivariable) for both clinical outcomes.
Baseline Tumor, Patient, and Treatment Characteristics for All Patients and Stratified by Center/Cohort
Kaplan-Meier Analyses
Kaplan-Meier analysis of the existing GUROC and NCCN classification system shows good separation of ASTRO BFFS and overall survival curves for all ProCaRS patients (log-rank P<.0001; Figure 1), all external-beam patients (log-rank P<.0001, curves not shown), all external-beam patients at 70 Gy in 35 fractions equivalent or greater (log-rank P<.0001, curves not shown). GUROC and NCCN classification systems separated groups into distinct ASTRO BFFS groups in the brachytherapy population (log-rank P<.005, curves not shown). Additionally, the NCCN classification was able to categorize brachytherapy patients into distinct overall survival groups (log rank P=.001), whereas the GUROC classification was not able to resolve these groups to the same extent (log rank P=.264).
Discussion
Risk stratification in the management of prostate cancer allows for prediction of the risk associated with various positive or negative clinical outcomes, the directing of appropriate therapy, clinical trial stratification and inclusion/exclusion criteria definition, and a common nomenclature for institutional outcome reporting. No consensus exists for the ideal risk stratification scheme, as evidenced by the multiple different (but related) systems20 that have been adopted by various organizations interested in the management of prostate cancer.21 Other risk stratification systems have been proposed in the literature, but these systems have not gained widespread acceptance.27,28
Summary of Hormone Therapy Use and Hormone Duration Stratified by GUROC and NCCN Risk Classification Systems (N=7974)
In addition to the reporting of important clinical outcomes related to radiotherapy treatment, this newly formed ProCaRS database will allow for the systematic investigation of various aspects of nonmetastatic prostate cancer pretreatment risk stratification/categorization. Other large prostate cancer databases exist in the literature, including CaPSURE,29 SEER and SEER-Medicare,30 and the Memorial Sloan-Kettering Cancer Center,31 which have investigated various research questions, including comparative effectiveness analyses, risk modeling/nomograms, population-based analyses, and treatment use.
This ProCaRS analysis has confirmed previously investigated prognostic factors related to ASTRO BFFS and overall survival, including various patient (age), tumor (PSA, Gleason score, positive core percentage), and treatment (hormonal therapy duration, form of radiation therapy) factors.20 Missing prediagnostic PSA data within the ProCaRS dataset prevented robust investigation into novel prognostic factors, such as PSA velocity and doubling time, in terms of the univariable and multivariable analyses. This analysis strongly confirms the utility of existing risk stratification schema to categorize different radiotherapy patients (external-beam radiation therapy vs brachytherapy) into different prognostic groups. Furthermore, analysis of the ProCaRS dataset suggested that the implementation of more than 3 risk categories may have important prognostic benefits for overall survival. In particular, inspection of clinical outcomes related to the “very-low-risk” and “low-risk” NCCN categories shows a difference between groups in terms of the BFFS end point but not the prostate cancer-specific survival or overall survival end points (Figure 1). This Kaplan-Meier information confirms the high curability of all forms of low-risk prostate cancer, likely because of the underlying indolent nature.
Risk Stratification and Clinical Outcomes for All Patients and Stratified by Center/Cohort
Multivariable Cox Regression Models of ASTRO BFFS and Overall Survival for the ProCaRS Database (N=7974)
This investigation has several limitations. The ProCaRS database is a retrospective entity based on data that were made available by various Canadian clinical investigators. The patients contained within the database may not be completely representative of a great patient population. Additionally, increasing levels of hormonal therapy were used with increasing risk stratification, which may impact the interpretation of the Kaplan-Meier curves. Evidence of this phenomenon exists with the observation of the large low-risk and LDR brachytherapy patient populations contained within ProCaRS. Despite the extensive quality assurance procedures used in the creation and curation of this database, it is important to acknowledge that different investigators at different institutions probably collected the constituent datasets during different periods. Despite these potential issues, the sample size of this dataset allows for high levels of statistical power to test new and novel hypotheses related to risk stratification and/or outcome analyses. Another limitation of this analysis was that it was restricted to common risk stratification systems used in clinical practice; alternative published multivariable systems/models were not considered in this initial investigation.
Subsequent to the establishment of this large prostate cancer radiotherapy database, the GUROC research group initiated several investigations related to the prediction of clinical outcomes. Specifically, they are using recursive partitioning analysis, clinical nomogram creation, and artificial neural network techniques to create de novo risk classification systems for the prediction of ASTRO BFFS and overall survival. Head-to-head comparisons of various alternative multivariable models/systems can also be conducted using the ProCaRS dataset. Additionally, they are investigating the utility of adapting the existing 3-group classification system (ie, the GUROC system) to allow for the splitting of low-, intermediate-, and high-risk categories into clinically relevant subcategories. In this project, up to 6 categories (very-low-risk, low-risk, low-intermediate-risk, high-intermediate-risk, high-risk, and extreme-risk) will be defined and characterized in terms of ASTRO BFFS and overall survival. In addition, the ProCaRS database will be used to perform direct propensity score matched pair analyses of various interventions (eg, brachytherapy vs external-beam radiation therapy).
Conclusions
This investigation has shown that cross-institutional collaboration to obtain robust clinical data to assess prostate cancer radiotherapy outcomes and to explore risk stratification is feasible. This investigation has demonstrated the ability of 2 risk stratification schemes to categorize patients who underwent external-beam radiation therapy or brachytherapy into various risk groups for 2 important clinical outcomes (ASTRO BFFS and overall survival). This work suggests that the definition of additional risk categories (ie, very-low-risk and very-high-risk consistent with the NCCN approach) may further improve patient risk categorization, particularly in brachytherapy patient populations. Further work will assess the importance of novel risk stratification systems, clinical nomograms, artificial neural networks, and direct treatment comparisons.
The authors would like to acknowledge the financial support of the Canadian Association of Radiation Oncology ACURA research award and the Ontario Institute of Cancer Research - High Impact Trials Program (Oncology Research and Methods Training Project).
References
- 1.↑
Kollmeier MA, Zelefsky MJ. How to select the optimal therapy for early-stage prostate cancer. Crit Rev Oncol Hematol 2012;83:225–234.
- 2.↑
Warde P, Mason M, Ding K et al.. Combined androgen deprivation therapy and radiation therapy for locally advanced prostate cancer: a randomized phase 3 trial. Lancet 2011;378:2104–2111.
- 3.↑
Pasquier D, Ballereau C. Adjuvant and salvage radiotherapy after prostatectomy for prostate cancer: a literature review. Int J Radiat Oncol Biol Phys 2008;72:972–979.
- 4.↑
Viani GA, da Silva LG, Stefano EJ. High-dose conformal radiotherapy reduces prostate cancer-specific mortality: results of a meta-analysis. Int J Radiat Oncol Biol Phys 2012;83:e619–625.
- 5.↑
Cuppone F, Bria E, Giannarelli D et al.. Impact of hormonal treatment duration in combination with radiotherapy for locally advanced prostate cancer: meta-analysis of randomized trials. BMC Cancer 2010;10:675.
- 6.↑
Peinemann F, Grouven U, Bartel C et al.. Permanent interstitial low-dose-rate brachytherapy for patients with localised prostate cancer: a systematic review of randomised and nonrandomised controlled clinical trials. Eur Urol 2011;60:881–893.
- 7.↑
Ghilezan M. Role of high dose rate brachytherapy in the treatment of prostate cancer. Cancer Radiother 2012;16:418–422.
- 8.↑
Vavassis P, Nguyen DH, Bahary JP et al.. Hypofractionated radiotherapy in prostate cancer. Expert Rev Anticancer Ther 2012;12:965–972.
- 9.↑
Arcangeli S, Scorsetti M, Alongi F. Will SBRT replace conventional radiotherapy in patients with low-intermediate risk prostate cancer? A review. Crit Rev Oncol Hematol 2012;84:101–108.
- 10.↑
Allen AM, Pawlicki T, Dong L et al.. An evidence based review of proton beam therapy: the report of ASTRO’s emerging technology committee. Radiother Oncol 2012;103:8–11.
- 11.↑
D’Amico AV, Whittington R, Malkowicz SB 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.
- 12.↑
Thompson I, Thrasher JB, Aus G et al.. Guideline for the management of clinically localized prostate cancer: 2007 update. J Urol 2007;177:2106–2131.
- 13.↑
Heidenreich A, Aus G, Bolla M et al.. EAU guidelines on prostate cancer. Eur Urol 2008;53:68–80.
- 14.↑
Lukka H, Warde P, Pickles T et al.. Controversies in prostate cancer radiotherapy: consensus development. Can J Urol 2001;8:1314–1322.
- 15.↑
Graham J, Baker M, Macbeth F et al.. Diagnosis and treatment of prostate cancer: summary of NICE guidance. BMJ 2008;336:610–612.
- 16.↑
Cooperberg MR, Lubeck DP, Mehta SS et al.. Time trends in clinical risk stratification for prostate cancer: implications for outcomes (data from CaPSURE). J Urol 2003;170:S21–25.
- 17.↑
Horwich A, Parker C, Bangma C et al.. Prostate cancer: ESMO Clinical Practice Guidelines for diagnosis, treatment and follow-up. Ann Oncol 2010;21:v129–33.
- 18.↑
Mohler JL, Armstrong AJ, Bahnson RR et al.. Prostate Cancer, Version 3.2012. Featured updates to the NCCN Guidelines. J Natl Compr Canc Netw 2012;10:1081–1087.
- 20.↑
Sutcliffe P, Hummel S, Simpson E et al.. Use of classical and novel biomarkers as prognostic risk factors for localised prostate cancer: a systematic review. Health Technol Assess 2009;13:iii, xi-xiii 1-219.
- 21.↑
Rodrigues G, Warde P, Pickles T et al.. Pre-treatment risk stratification of prostate cancer patients: A critical review. Can Urol Assoc J 2012;6:121–127.
- 22.↑
Zaorsky NG, Li T, Devarajan K, Horwitz EM et al.. Assessment of the American Joint Committee on Cancer staging (sixth and seventh editions) for clinically localized prostate cancer treated with external beam radiotherapy and comparison with the National Comprehensive Cancer Network risk-stratification method. Cancer 2012;118:5535–5543.
- 23.↑
Pickles T, Morris WJ, Kattan MW et al.. Comparative 5-year outcomes of brachytherapy and surgery for prostate cancer. Brachytherapy 2011;10:9–14.
- 24.↑
Roach M 3rd, Hanks G, Thames H Jr et al.. Defining biochemical failure following radiotherapy with or without hormonal therapy in men with clinically localized prostate cancer: recommendations of the RTOG-ASTRO Phoenix Consensus Conference. Int J Radiat Oncol Biol Phys 2006;65:965–974.
- 25.↑
Horwitz EM, Vicini FA, Ziaja EL et al.. The correlation between the ASTRO Consensus Panel definition of biochemical failure and clinical outcome for patients with prostate cancer treated with external beam irradiation. American Society of Therapeutic Radiology and Oncology. Int J Radiat Oncol Biol Phys 1998;41:267–272.
- 26.↑
Thompson A, Keyes M, Pickles T et al.. Evaluating the Phoenix definition of biochemical failure after (125)I prostate brachytherapy: can PSA kinetics distinguish PSA failures from PSA bounces? Int J Radiat Oncol Biol Phys 2010;78:415–421.
- 27.↑
Beasley M, Williams SG, Pickles T et al.. Expanded risk groups help determine which prostate radiotherapy sub-group may benefit from adjuvant androgen deprivation therapy. Radiat Oncol 2008;3:8.
- 28.↑
Roach M 3rd, Weinberg V, Sandler H et al.. Staging for prostate cancer: time to incorporate pretreatment prostate-specific antigen and Gleason score? Cancer 2007;109:213–220.
- 29.↑
Lubeck DP, Litwin MS, Henning JM et al.. The CaPSURE database: a methodology for clinical practice and research in prostate cancer. CaPSURE Research Panel. Cancer of the Prostate Strategic Urologic Research Endeavor. Urology 1996;48:773–777.
- 30.↑
Siegel R, DeSantis C, Virgo K et al.. Cancer treatment and survivorship statistics, 2012. CA Cancer J Clin 2012;62:220–241.
- 31.↑
Zelefsky MJ, Kattan MW, Fearn P et al.. Pretreatment nomogram predicting ten-year biochemical outcome of three-dimensional conformal radiotherapy and intensity-modulated radiotherapy for prostate cancer. Urology 2007;70:283–287.