- Open Access
The impact of increasing dose on overall survival in prostate cancer
Radiation Oncology volume 10, Article number: 115 (2015)
To assess the impact of increasing dose on overall survival (OS) for prostate cancer patients.
Treatment data were obtained on more than 20,000 patients in the National Oncology Data Alliance®, a proprietary database of merged tumor registries, who were treated for prostate cancer with definitive radiotherapy between 1995 and 2006. Eligible patients had complete data on total dose, T stage, use and timing of androgen deprivation therapy (ADT), and treatment start date (n = 20,028). Patients with prior malignancies were excluded.
On multivariate analysis, dose, T stage, grade, marital status, age, and neoadjuvant ADT were significant predictors of OS. Hazard ratios for OS declined monotonically with increasing dose, reaching 0.63 (95 % Confidence Interval 0.53–0.76) at ≥80 Gy. On subset analysis, neoadjuvant ADT significantly improved OS in high risk patients but was not significant in lower risk patients. The dose response was maintained across all risk groups. Medical comorbidities were balanced across all dose strata and sensitivity analysis demonstrated that other prognostic factors were unlikely to explain the observed dose response.
This study suggests that increasing dose significantly improves OS in prostate cancer patients treated with radiotherapy.
Dose-escalated radiotherapy provides superior freedom from clinical and biochemical failure in men treated for prostate cancer (PCa) [1–5] and is widely used in clinical practice . Combined androgen deprivation therapy (ADT) and radiotherapy have demonstrated an OS advantage in intermediate and high risk patients in multiple randomized clinical trials (RCTs), albeit at doses ≤70 Gy [7–10]. Retrospective series have suggested that the benefit of dose escalation may be greater than that of ADT with conventional doses . With a median follow-up of 7 years, Radiation Therapy Oncology Group (RTOG) 0126 showed significant improvements in the rates of biochemical failure and distant metastases with dose-escalated radiotherapy, but no OS benefit has been observed with higher doses to date .
The optimal radiotherapy dose must balance the risks and benefits of improved local control relative to treatment morbidity, but a treatment benefit in terms of improved OS is generally assumed [13, 14]. Given the potential toxicities of therapy and growing evidence that the treatment benefit in men with indolent PCa is limited [15, 16], it is important to analyze the outcomes of dose-escalated therapy in a population-based cohort to characterize this benefit where randomized data with adequate statistical power have been lacking. Therefore, the purpose of this study is to assess the impact on OS of increasing dose in PCa patients treated with definitive radiotherapy using a large database and to quantify the effects of prognostic factors.
Data on all men diagnosed with PCa from 1995 through 2006 were extracted from the National Oncology Data Alliance® (Elekta/IMPAC Medical Systems, Inc., Sunnyvale, CA and Stockholm, Sweden), a proprietary database of merged tumor registries. The registry is fully compliant with American College of Surgeons (ACOS) regulatory requirements. It captures newly diagnosed cancer cases at more than 150 hospitals in the United States. Data in the NODA registries are exactly the data sent to state tumor registries and to Surveillance, Epidemiology, and End Results (SEER), in regions that participate in SEER .
The study population included all patients treated with definitive radiotherapy with complete data on total dose, T stage, use and timing of ADT, and treatment start dates. Patients who underwent prostatectomy or brachytherapy, had distant disease, or had history of prior malignancy were excluded, leaving 20,028 evaluable patients. Patient characteristics included age, race, marital status, T stage, grade, total dose, ADT use, year of diagnosis, year of treatment, and medical comorbidities. In ACOS registries, Gleason scores are binned into four histologic grades . PSA values were not available. This study was conducted with the approval of the City of Hope Institutional Review Board.
The primary endpoint was OS, measured from the start of treatment to the date of death from any cause. Participating institutions regularly update patient files with vital status. For patients reported alive at last follow-up, the Social Security death index (SSDI) was searched to extend follow-up duration. OS outcomes were updated through November 1, 2011, the last date that state death records were included in the Death Master File by the Social Security Administration as per Section 205r of the Social Security Act. Patients alive as of this date were censored.
Statistical analysis was performed using SPSS v. 18.0 (SPSS Inc., Chicago) after patient-identifiers were removed from the data set. Multivariate proportional hazards analysis (MVA) was used to identify factors associated with OS . The proportional hazards assumption was tested for a non-zero slope in the generalized linear model using scaled Schoenfeld residuals, and no substantial deviations were detected. Survival curves were plotted by Kaplan-Meier and compared using the log-rank test. For survival analysis, dose was considered as a categorical variable with three dose bins (68–71.99 Gy, 72–75.99 Gy, and ≥76 Gy). For MVA, the effect of dose was also analyzed as a categorical variable and was divided into seven dose bins with 68–69.99 Gy as the referent category. Correlation between significant prognostic variables and dose was examined by least squares regression and by plotting studentized residuals. MVA was performed with variates common to both the NODA and the Surveillance, Epidemiology, and End Results (SEER) database to confirm the validity of the data and the findings . MVA yielded similar hazard ratios (HRs) and demonstrated congruence between the two databases (Additional file 1: Appendix 1). To address the potential differences between this retrospective study and RCTs, 5-year OS was calculated for our data and RCTs using high risk patients in both cases.
The observed differences in OS between groups receiving different radiotherapy doses may reflect the effects of an unmeasured variable, like PSA, leading to a false association between dose and OS. For example, assigning patients to lower doses because they have higher PSAs and reduced life expectancy could mimic a dose response. Sensitivity analysis was performed to measure the potential influence of an unmeasured confounder on the hazard ratio (HR) estimate for OS. The association between dose and OS may be affected by the prevalences of the unmeasured confounder in the dose groups and its HR (HR confounder ), which was assumed to be independent of dose. Only sampling bias for an unmeasured confounder with HR > 1 and higher prevalence in the lower dose group can mimic the dose response [20, 21]. Sensitivity analysis estimated the necessary prevalence to mimic the perceived dose response based on the presence of the confounder. Williams et al. showed that patients had significantly worse OS with pretreatment 40 ≥ PSA > 20 (HR = 1.32) and PSA > 40 (HR = 1.91), compared to patients with PSA < 10. . Therefore, we modeled the effect of the putative confounder as a PSA imbalance using these published values.
Following the work of Lin et al., the degree of imbalance needed between dose groups to mimic the observed dose response is given by,
where P low and P high are the prevalences of the confounder in the low and high dose groups, respectively, and HR dose is the observed HR for dose from the MVA (See Additional file 1: Appendix 2).
Patient demographics and tumor characteristics are presented in Table 1. Median age was 71 years. ADT was administered in 40 % of patients as neoadjuvant or concurrent therapy with radiotherapy, defined as beginning within 6 months of the radiotherapy start date. The first notation for the use of intensity-modulated radiotherapy (IMRT) occurred in 2002. The prevalence of significant comorbidities was similar across all dose groups.
Median follow-up was 8.6 years for surviving patients. At the close-out date, 40 % of patients had died. The Kaplan-Meier curves showed a clear and statistically significant dose response relationship for OS (p < 0.0001) (Fig. 1). Median OS was 11.4 years (95 % Confidence Interval [CI] 11.1–11.6) for 68–71.99 Gy, 12.0 years (95 % CI 11.6–12.3) for 72–75.99 Gy, and 12.8 years (95 % CI 10.9–14.7) for ≥76 Gy.
Dose, T stage, grade, age, marital status, and neoadjuvant ADT were significant on MVA (Table 2). Race and treatment year were not significant. When dose was defined as a categorical variable with seven dose bins using 68–69.99 Gy as the referent category, the same variables were significant and no additional factors were significant. Adjusting for the effect of other prognostic factors, HRs for OS declined monotonically with increasing dose, reaching 0.63 (95 % CI 0.53–0.76) at ≥80 Gy, without evidence of a plateau (Fig. 2 ). Even using the 74–76 Gy group (HR = 0.83, 95 % CI 0.77–0.90) as the referent category, treatment with ≥80 Gy maintained its significant OS advantage.
On MVA, patients treated with higher doses also experienced a significant improvement in cause-specific survival (CSS). Using 68–71.99 Gy as the referent category, the HRs for CSS decreased sequentially to 0.62 (95 % CI 0.49–0.78) for 72–75.99 Gy and 0.26 (95 % CI 0.12–0.55) for ≥76 Gy.
As an exploratory analysis, patients were organized into high and low risk groups, approximating the National Comprehensive Cancer Network risk categories. Patients with either T3-T4 stage or Grade 3–4 histology were classified as high risk. In our analysis, T1 and T2, Grade 1 and T1, Grade 2 patients had comparable HRs for OS and were classified as low risk and T2 Grade 2 patients were categorized as intermediate. Dose retained significance on MVA for low, intermediate, and high risk patients, with comparable HRs for dose observed in all three subgroups. On subset analysis, ADT provided a significant OS benefit in high risk patients (HR 0.87, 95 % CI 0.81–0.94), but was not significant in lower risk patients on MVA (p > 0.5).
Dose was not correlated with age, race, marital status, T stage, grade, and neoadjuvant ADT. IMRT was first noted in patient records in 2002. When all patients treated in 2002 or later were excluded, dose remained statistically significant on MVA with small changes in the observed HRs. In the multivariate setting, there was no suggestion or trend of improved OS in later years (p > 0.7). Thus, the dose response does not appear to be an effect of treatment technique.
Treatment outcomes were examined with patients divided into three age subgroups (< 65, 65–74, and ≥75 years) to determine if the benefit associated with high dose treatment was greater in younger patients. HRs for OS were comparable in patients of all ages, adjusting for the effects of increasing age on mortality.
The association between dose and improved OS was evident in all subgroups examined, including stratification by ADT use. The observed OS advantage remained significant when patients who died within 6 months, 1, 2, 3, and 5 years of treatment were excluded, indicating that the benefit of high dose treatment was not related to an imbalance in the number of patients with abbreviated life expectancy consigned to receive lower doses. Reported medical comorbidities were balanced across all dose strata.
Five-year OS estimates for our study and for RCTs enrolling high risk patients treated with external beam radiotherapy and ADT are shown in Fig. 3 [10, 23–26]. All patients in our study had a minimum 5-year follow-up, making 5-year OS a fair metric for comparison. Five-year OS was similar for RCTs and our cohort, suggesting not only that a very similar dose response was operative in both this study and the RCTs, but that the demographics and survival from non-PCa-related deaths were similar for all of these patient populations.
Sensitivity analysis was performed to assess the potential effects of an unmeasured confounder on the estimated HR for dose. The observed dose response could be simulated mathematically if patients in the high dose groups had lower pretreatment PSAs than patients in the lower dose groups. To determine the degree of PSA disparity needed to achieve the observed HR for OS, we assumed that all patients in the ≥80 Gy group had PSA < 10. We then calculated the prevalence of patients with elevated pretreatment PSA in the 68–70 Gy group that must be present to match the observed HR dose = 0.63 (Table 2). To replicate the dose response, the 68–70 Gy group must have either 65 % of patients with PSA > 40 (HR = 1.91) and the remainder with PSA < 10 or 45 % with PSA > 40 and the remainder with 40 ≥ PSA > 20 (HR = 1.32) . Thus, the median PSA must be approximately 40 in the 68–70 Gy group to account for the observed OS difference.
A disparity in the prevalence of a medical comorbidity, such as cardiovascular disease, could also have this effect, but the recorded prevalence of significant comorbidities was similar across all dose groups.
This observational study is the largest on PCa treated with known radiotherapy doses. On MVA, dose was a clear and statistically significant prognostic factor for OS without evidence of plateau up to 82 Gy. On subset analysis, ADT was protective in high risk patients but was not significant in low risk patients, whereas dose was significant in both low and high risk patients. The association between dose and improved OS was evident in all subgroups examined, including stratification by age and treatment with and without ADT. Although a retrospective analysis cannot account for all patient factors, the dose response was durable across a number of clinical scenarios and this analysis accounted for potential disparities in medical comorbidities. Sensitivity analysis also suggested that the presence of other uncoded variables, such as PSA, was unlikely to explain the observed dose response.
Several limitations in our study should be acknowledged. Although this study is the largest to address the impact of increasing dose on OS in patients with localized PCa treated with definitive radiotherapy, the data do not enable the opportunity to examine several important factors, including duration of ADT, use of pelvic nodal radiotherapy, and complications. Second, outcome measures were limited to OS and CSS, and data regarding local and distant failure cannot be extracted. Third, patient selection can influence the results in population-based observational series and must be considered when applying our results to clinical practice.
The completeness of the available data, however, helps to counter the possibility that the observed benefit of high dose radiotherapy is due solely to selection bias. First, the advantage of high dose remained when patients with shortened life expectancy (dying within 6 months to 5 years after treatment) were excluded. Second, the HRs for dose were consistent across all age groups, indicating that the benefit of high dose therapy was not biased by age discrepancies. Third, sensitivity analysis demonstrated that the observed dose response could be mimicked only by an unrealistic disparity in the prevalence of an unmeasured variable such as PSA that severely compromised survival. The agreement of this study with RCTs (Fig. 3 ) strongly imply that patients in this retrospective study were no more likely, at any dose, to have life-threatening comorbidities than patients in RCTs. Consequently, the accuracy of HR estimates may be affected by unmeasured variables, but reason dictates that these findings represent a true biological dose response.
RCTs provide precise measures and comparisons of treatment efficacy, have excellent internal validity, and are less prone to bias. Population-based observational studies cannot replace randomized data as the standard for outcomes research; however, it can provide evidence of the effectiveness of therapies in the general population and it can address questions that have not been adequately evaluated in RCTs . The impact of bias can be minimized in observational studies by accounting for important covariates, including comorbidity, that are controlled in RCTs. In addition, the impact of unmeasured variables, such as PSA, can also be estimated using sensitivity analysis (as we did here). With careful planning, population-based observational studies can help to answer important clinical questions and fill the knowledge gaps left by existing randomized data.
Of note, PSA cannot be reliably extracted from US population-based databases before 2004, and changes in Gleason score measurements over time could result in upgrading of tumors diagnosed in later years .
Finally, increasing dose conferred a clear and statistically significant OS benefit, but the magnitude of benefit was modest. In this study, the median survival for patients treated with ≥76 Gy was 12.8 years compared to 11.4 years for 68–72 Gy. Men diagnosed with PCa have competing causes of mortality, including cardiovascular disease and second cancers [29, 30]. Appropriate patient selection is clearly an important consideration in the decision to treat that is beyond the scope of this study. Once this decision is made, however, our study suggests that treatment with higher dose and ADT in high risk patients confers a significant OS benefit. Only a prospective RCT can provide absolute measures of this advantage.
Five randomized dose escalation trials have shown superior freedom from failure at higher doses, but were underpowered to detect an OS difference [1–5]. In 393 patients receiving radiotherapy followed by proton boost to 70.2 vs. 79.2 GyE, Zietman et al. reported 5-year OS rates of 96 % and 97 %, respectively. Given the 96 % OS for the low dose group, no survival advantage could have been demonstrated in this study . Kuban et al. reported no difference in 8-year OS in 303 patients treated with either 70 or 78 Gy. However, 73 % of patients were censored, salvage ADT was used in an unknown number of cases, and only 10 patients died of PCa . In the largest dose escalation trial, Dearnaley et al. reported 10-year OS of 71 % in 843 patients treated with either 64 or 74 Gy. All patients received neoadjuvant and concurrent ADT, regardless of risk status .
RTOG 0126 demonstrated a significant reduction in biochemical failure and distant metastases with higher doses, but not an improvement in OS after a median follow-up of 7 years. Given the small number of CSS events to date, however, this finding is not unexpected and longer follow-up is required to credibly address this question, even in an adequately powered clinical trial .
In our study, increasing dose was associated with a significant OS benefit. The small incremental benefit of 1.4 years median survival, however, required thousands of patients to detect this difference at a level of statistical significance. A dose response can only be conclusively demonstrated in a clinical trial. To answer this question, two RCTs in the United States that were powered to detect an OS difference are maturing and ongoing. Our findings, however, suggest that although a dose response for OS exists, the magnitude of this advantage may compromise their ability to verify this outcome.
In a large population-based cohort, our study suggests that OS is significantly improved in PCa patients in all risk groups treated with higher doses. The analysis also indicates that medical comorbidities, PSA, and other prognostic factors cannot credibly account for the observed dose response. We conclude that treatment with higher dose should be considered in all PCa patients receiving definitive radiotherapy provided that acceptable toxicity constraints can be achieved.
Zietman AL, DeSilvio ML, Slater JD, Rossi Jr CJ, Miller DW, Adams JA, et al. Comparison of conventional-dose vs. high-dose conformal radiation therapy in clinically localized adenocarcinoma of the prostate. JAMA. 2005;294:1233–39.
Dearnaley DP, Sydes MR, Graham JD, Aird EG, Bottomley D, Cowan RA, et al. Escalated-dose versus standard-dose conformal radiotherapy in prostate cancer: First results from the MRC RT01 randomised controlled trial. Lancet Oncol. 2007;8:475–87.
Al-Mamgani A, van Putten WL, Heemsbergen WD, van Leenders GJ, Slot A, Dielwart MF, et al. Update of Dutch multicenter dose-escalation trial of radiotherapy for localized prostate cancer. Int J Radiat Oncol Biol Phys. 2008;72(4):980–8.
Kuban DA, Tucker SL, Dong L, Starkschall G, Huang EH, Cheung MR, et al. Long-term results of the M. D. Anderson randomized dose-escalation trial for prostate cancer. Int J Radiat Oncol Biol Phys. 2008;70(1):67–74.
Beckendorf V, Guerif S, Le Prisé E, Cosset JM, Bougnoux A, Chauvet B, et al. 70 Gy vs. 80 Gy in localized prostate cancer: 5-year results of GETUG 06 randomized trial. Int J Radiat Oncol Biol Phys. 2011;80(4):1056–63.
Swisher-McClure S, Mitra N, Woo K, Smaldone M, Uzzo R, Bekelman JE. Increasing use of dose-escalated external beam radiation therapy for men with nonmetastatic prostate cancer. Int J Radiat Oncol Biol Phys. 2014;89(1):103–12.
Bolla M, Collette L, Blank L, Warde P, Dubois JB, Mirimanoff RO, et al. Long-term results with immediate androgen suppression and external irradiation in patients with locally advanced prostate cancer (an EORTC study): a phase III randomised trial. Lancet. 2002;360(9327):103–6.
Jones CU, Hunt D, McGowan DG, Amin MB, Chetner MP, Bruner DW, et al. Radiotherapy and short-term androgen deprivation for localized prostate cancer. N Engl J Med. 2011;365(2):107–18.
Widmark A, Klepp O, Solberg A, Damber JE, Angelsen A, Fransson P, et al. Endocrine treatment, with or without radiotherapy, in locally advanced prostate cancer (SPCG-7/SFUO-3): an open randomised phase III trial. Lancet. 2009;373(9660):301–8.
Warde P, Mason M, Ding K, Kirkbride P, Brundage M, Cowan R, et al. Combined androgen deprivation therapy and radiation therapy for locally advanced prostate cancer: a randomised, phase 3 trial. Lancet. 2014;378(9809):2104–11.
Nguyen KH, Horwitz EM, Hanlon AL, Uzzo RG, Pollack A. Does short-term androgen deprivation substitute for radiation dose in the treatment of high-risk prostate cancer? Int J Radiat Oncol Biol Phys. 2003;57:377–83.
Michalski JM, Moughan J, Purdy J, Bosch W, Bahary J-P, Lau HY, et al. A randomized trial of 79.2 Gy versus 70.2 Gy radiation therapy (RT) for localized prostate cancer. J Clin Oncol. 2015;33(suppl 7S; abstr 4):4.
Zelefsky MJ, Fuks Z, Hunt M, Lee HJ, Lombardi D, Ling CC, et al. High dose radiation delivered by intensity modulated conformal radiotherapy improves the outcome of localized prostate cancer. J Urol. 2001;166(3):876–81.
Chauvet B, Alfonsi M, de Rauglaudre G, Reboul F. Prostate cancer: has local radiation treatment had an impact on survival? Cancer Radiother. 2002;6(3):141–6.
Bill-Axelson A, Holmberg L, Ruutu M, Garmo H, Stark JR, Busch C, et al. Radical prostatectomy versus watchful waiting in early prostate cancer. N Engl J Med. 2011;364:1708–17.
Wilt TJ, Brawer MK, Jones KM, Barry MJ, Aronson WJ, Fox S, et al. Radical prostatectomy versus observation for localized prostate cancer. N Engl J Med. 2012;367:203–13.
Surveillance, Epidemiology, and End Results (SEER) Program (www.seer.cancer.gov) SEER*Stat Database: Incidence - SEER 18 Regs Research Data, Nov 2012 Sub (1973–2010) < Katrina/Rita Population Adjustment > − Linked To County Attributes - Total U.S., 1973–2012 Counties, National Cancer Institute, DCCPS, Surveillance Research Program, Surveillance Systems Branch, released September 2013
Bilimoria KY, Stewart AK, Winchester DP, Ko CY. The national cancer data base: a powerful initiative to improve cancer care in the United States. Ann Surg Oncol. 2007;15(3):683–90.
Harrell Jr FE. Regression modeling strategies with applications to linear models, logistic regression, and survival analysis. New York: Springer. 2014;13–15:509–22.
Lin DY, Psaty BM, Kronmal RA. Assessing the sensitivity of regression results to unmeasured confounders in observational studies. Biometrics. 1998;54(3):948–63.
Mitra N, Heitjan DF. Sensitivity of the hazard ratio to nonignorable treatment assignment in an observational study. Statist Med. 2007;26:1398–414.
Williams SG, Duchesne GM, Millar JL, Pratt GR. Both pretreatment prostate-specific antigen level and posttreatment biochemical failure are independent predictors of overall survival after radiotherapy for prostate cancer. Int J Radiat Oncol Biol Phys. 2004;60(4):1082–7.
Armstrong JG, Gillham CM, Dunne MT, Fitzpatrick DA, Finn MA, Cannon ME, et al. A randomized trial (Irish clinical oncology research group 97–01) comparing short versus protracted neoadjuvant hormonal therapy before radiotherapy for localized prostate cancer. Int J Radiat Oncol Biol Phys. 2011;81(1):35–45.
Bolla M, de Reijke TM, Van Tienhoven G, Van den Bergh AC, Oddens J, Poortmans PM, et al. Duration of androgen suppression in the treatment of prostate cancer. N Engl J Med. 2009;360(24):2516–27.
Arcangeli G, Fowler J, Gomellini S, Arcangeli S, Saracino B, Petrongari MG, et al. Acute and late toxicity in a randomized trial of conventional versus hypofractionated three-dimensional conformal radiotherapy for prostate cancer. Int J Radiat Oncol Biol Phys. 2011;79(4):1013–21.
Mendenhall NP, Hoppe BS, Nichols RC, Mendenhall WM, Morris CG, Li Z. Five-year outcomes from 3 prospective trials of image-guided proton therapy for prostate cancer. Int J Radiat Oncol Biol Phys. 2014;88(3):596–602.
Booth CM, Tannock IF. Randomised controlled trials and population-based observational research: partners in the evolution of medical evidence. Br J Cancer. 2014;110(3):551–5.
Jani AB, Master VA, Rossi PJ, Liauw SL, Johnstone PA. Grade migration in prostate cancer: an analysis using the Surveillance, Epidemiology, and End Results registry. Prostate Cancer Prostatic Dis. 2007;10(4):347–51.
Tendulkar RD, Hunter GK, Reddy CA, Stephans KL, Ciezki JP, Abdel-Wahab M, et al. Causes of mortality after dose-escalated radiation therapy and androgen deprivation for high-risk prostate cancer. Int J Radiat Oncol Biol Phys. 2013;87(1):94–9.
Harlan LC, Potosky A, Gilliland FD, Hoffman R, Albertsen PC, Hamilton AS, et al. Factors associated with initial therapy for clinically localized prostate cancer: prostate cancer outcomes study. J Natl Cancer Inst. 2001;93(24):1864–71.
The authors declare that they have no competing interests.
MH participated in the design of the study, performed the statistical analysis, and drafted the manuscript. TS conceived of the study, participated in its design and coordination, performed the statistical analysis, and helped to draft the manuscript. DS participated in the design of the study and performed the statistical analysis. BT participated in the design of the study. JW participated in the design and coordination of the study. All authors read and approved the final manuscript.
About this article
Cite this article
Hall, M.D., Schultheiss, T.E., Smith, D.D. et al. The impact of increasing dose on overall survival in prostate cancer. Radiat Oncol 10, 115 (2015). https://doi.org/10.1186/s13014-015-0419-3
- Radiation therapy
- Prostate cancer
- Dose response
- Overall survival