Abstract
The impact of treatment delays on prostate cancer–specific outcomes remains ill-defined. This study investigates the effect of time to treatment on biochemical disease control after prostatectomy.
This retrospective study includes 1,807 patients who received a prostatectomy as a primary treatment at two large tertiary referral centers from 1987 to 2015. Multivariate cox model with restricted cubic spline was used to identify optimal time to receive treatment and estimate the risk of biochemical recurrence.
Median follow-up time of the study was 46 (interquartile range, 18–86) months. Time to treatment was subcategorized based on multivariate cubic spline cox model. In multivariate spline model, adjusted for all the pertinent pretreatment variables, inflection point in the risk of biochemical recurrence was observed around 3 months, which further increased after 6 months. Based on spline model, time to treatment was then divided into 0 to 3 months (61.5%), >3 to 6 months (31.1%), and 6 months (7.4%). In the adjusted cox model, initial delays up to 6 months did not adversely affect the outcome; however, time to treatment >6 months had significantly higher risk of biochemical recurrence (HR, 1.84; 95% confidence interval, 1.30–2.60; P < 0.01).
The initial delays up to 6 months in prostate cancer primary treatment may be sustainable without adversely affecting the outcome. However, significant delays beyond 6 months can unfavorably affect biochemical disease control.
Time to treatment can aid clinicians in the decision-making of prostate cancer treatment recommendation and educate patients against unintentional treatment delays.
Introduction
Time to receive treatment or treatment delays is an important quality metric in patient-centered care that has been shown to affect outcome among various cancers (1), including breast, head and neck, colorectal, and melanoma (2–6), with patients receiving treatment sooner after initial diagnosis being more likely to experience favorable treatment outcomes. However, the implications of treatment delays on prostate cancer prognosis remain unclear in contemporary literature. Unlike other cancers, prostate cancer usually has slow progression and does not present immediate risk to patient survival (7), rendering both clinicians and patients ample time to make decisions on the appropriate treatment choice. In addition, concerns related with overtreatment of prostate cancer have prompted clinicians to adapt more pragmatic approaches toward treatment delivery which may increase the response time to treatment initiation, leading to increased time to treatment (8).
Few studies have demonstrated the possibility of adverse outcomes related with increasing time to treatment in prostate cancer (9–11); however, the length of delays prior to the initiation of treatment which is sustainable to achieve durable prostate cancer outcomes, without losing window of cure, remains debatable (12, 13). In addition, studies have argued against the use of early intervention in prostate cancer (14) and have shown that delays up to several years do not adversely affect prostate cancer outcomes (15–17). Consequently, a lack of consensus among previously published studies has resulted in significant ambiguity in the association between time to treatment and prostate cancer outcomes (14, 15, 18–21), which warrants further investigation. This study leverages the availability of 25 years’ worth of data from two large tertiary referral centers, with detailed clinicopathologic information, on prostate cancer patients who received prostatectomy as part of their primary treatment. This study seeks to determine the optimal time to start treatment that maximizes favorable prostate cancer outcomes and the clinical implications of treatment delays on biochemical disease control after prostatectomy.
Materials and Methods
Patient population
This retrospective analysis includes a total of 1,891 patients with histologically confirmed localized nonactive surveillance prostate cancer who received prostatectomy as a primary treatment at two large tertiary referral centers—H. Lee Moffitt Cancer Center and Research Institute (MCC) and the University of Pennsylvania Health System (UPHS, Philadelphia, PA)—from 1987 to 2015 (Table 1). Patients who received prostatectomy at MCC were identified from the Health Research and Informatics platform, an enterprise-wide data warehouse. Patients from UPHS were those recruited to the Study of Clinical Outcomes, Risk and Ethnicity (SCORE) between 1990 and 2012. In addition, patients who received neoadjuvant treatment prior to their surgery were excluded from the analysis (n = 84). Therefore, a total 1,807 patients were used in the final analysis. Institutional Review Board approval was obtained before the commencement of the study.
Clinicopathologic/demographic characteristics . | MCC (n = 648) . | UPHS (n = 1,159) . |
---|---|---|
Time to treatment (months) | ||
Median (IQR) | 3 (2–4) | 3 (2–4) |
Time to treatment (months) | ||
0–3 | 435 (67.1) | 676 (58.3) |
>3–6 | 169 (26.1) | 393 (33.9) |
>6 | 44 (6.8) | 90 (7.8) |
Age category | ||
≤50 | 91 (14.0) | 105 (9.1) |
>50–60 | 282 (43.5) | 519 (44.8) |
>60–70 | 252 (38.9) | 507 (43.7) |
≥70 | 23 (3.5) | 28 (2.4) |
Age | ||
Median (IQR) | 59 (53–64) | 60 (55–64) |
Race | ||
AAM | 147 (22.7) | 147 (12.7) |
EAM | 501 (77.3) | 1,012 (87.3) |
NCCN risk group | ||
Low risk | 277 (42.7) | 637 (55.0) |
Intermediate risk | 285 (44.0) | 315 (27.2) |
High risk | 69 (10.6) | 84 (7.2) |
Unknown | 17 (2.6) | 123 (10.6) |
iPSA category | ||
0–6 | 374 (57.7) | 650 (56.1) |
>6–10 | 177 (27.3) | 317 (27.3) |
>10–20 | 77 (11.9) | 151 (13.0) |
>20 | 20 (3.1) | 41 (3.5) |
Clinical GS | ||
Group 1 | 335 (51.7) | 780 (67.3) |
Group 2/3 | 235 (36.3) | 183 (15.8) |
Group 4 | 41 (6.3) | 43 (3.7) |
Group 5 | 5 (0.8) | 12 (1.0) |
Unknown | 32 (4.9) | 141 (12.2) |
Clinical T stage | ||
T1 | 493 (76.1) | 649 (56.0) |
T2 | 127 (19.6) | 265 (22.9) |
T3 | 11 (1.7) | 9 (0.8) |
Unknown | 17 (2.6) | 236 (20.4) |
Pathology Gleason | ||
Group 1 | 240 (37.0) | 598 (51.6) |
Group 2/3 | 370 (57.1) | 479 (41.3) |
Group 4 | 21 (3.2) | 51 (4.4) |
Group 5 | 17 (2.6) | 31 (2.7) |
AJCC pathologic stage | ||
pT2A–pT2B | 517 (79.8) | 814 (70.2) |
pT3A | 104 (16.0) | 281 (24.2) |
pT3B–pT4 | 6 (0.9) | 63 (5.4) |
PX | 21 (3.2) | 1 (0.1) |
Marital status | ||
Married | 569 (85.0) | — |
Single | 100 (15.0) | — |
Insurance status | ||
Private insurance | 423 (63.2) | — |
Public insurance | 182 (27.2) | — |
Uninsured | 64 (9.6) | — |
History tobacco status | ||
Present | 310 (46.3) | — |
Absent | 325 (48.6) | — |
Unknown | 34 (5.1) | — |
Clinicopathologic/demographic characteristics . | MCC (n = 648) . | UPHS (n = 1,159) . |
---|---|---|
Time to treatment (months) | ||
Median (IQR) | 3 (2–4) | 3 (2–4) |
Time to treatment (months) | ||
0–3 | 435 (67.1) | 676 (58.3) |
>3–6 | 169 (26.1) | 393 (33.9) |
>6 | 44 (6.8) | 90 (7.8) |
Age category | ||
≤50 | 91 (14.0) | 105 (9.1) |
>50–60 | 282 (43.5) | 519 (44.8) |
>60–70 | 252 (38.9) | 507 (43.7) |
≥70 | 23 (3.5) | 28 (2.4) |
Age | ||
Median (IQR) | 59 (53–64) | 60 (55–64) |
Race | ||
AAM | 147 (22.7) | 147 (12.7) |
EAM | 501 (77.3) | 1,012 (87.3) |
NCCN risk group | ||
Low risk | 277 (42.7) | 637 (55.0) |
Intermediate risk | 285 (44.0) | 315 (27.2) |
High risk | 69 (10.6) | 84 (7.2) |
Unknown | 17 (2.6) | 123 (10.6) |
iPSA category | ||
0–6 | 374 (57.7) | 650 (56.1) |
>6–10 | 177 (27.3) | 317 (27.3) |
>10–20 | 77 (11.9) | 151 (13.0) |
>20 | 20 (3.1) | 41 (3.5) |
Clinical GS | ||
Group 1 | 335 (51.7) | 780 (67.3) |
Group 2/3 | 235 (36.3) | 183 (15.8) |
Group 4 | 41 (6.3) | 43 (3.7) |
Group 5 | 5 (0.8) | 12 (1.0) |
Unknown | 32 (4.9) | 141 (12.2) |
Clinical T stage | ||
T1 | 493 (76.1) | 649 (56.0) |
T2 | 127 (19.6) | 265 (22.9) |
T3 | 11 (1.7) | 9 (0.8) |
Unknown | 17 (2.6) | 236 (20.4) |
Pathology Gleason | ||
Group 1 | 240 (37.0) | 598 (51.6) |
Group 2/3 | 370 (57.1) | 479 (41.3) |
Group 4 | 21 (3.2) | 51 (4.4) |
Group 5 | 17 (2.6) | 31 (2.7) |
AJCC pathologic stage | ||
pT2A–pT2B | 517 (79.8) | 814 (70.2) |
pT3A | 104 (16.0) | 281 (24.2) |
pT3B–pT4 | 6 (0.9) | 63 (5.4) |
PX | 21 (3.2) | 1 (0.1) |
Marital status | ||
Married | 569 (85.0) | — |
Single | 100 (15.0) | — |
Insurance status | ||
Private insurance | 423 (63.2) | — |
Public insurance | 182 (27.2) | — |
Uninsured | 64 (9.6) | — |
History tobacco status | ||
Present | 310 (46.3) | — |
Absent | 325 (48.6) | — |
Unknown | 34 (5.1) | — |
Abbreviations: AJCC, American Joint Committee on Cancer; iPSA, preoperative PSA.
Baseline covariates
All the patients included in the study received prostatectomy as definite primary treatment for their localized prostate cancer. Detailed clinicopathologic and demographic information including pretreatment variables such as clinical Gleason score (GS), PSA, and National Comprehensive Cancer Network (NCCN) risk grouping on patients from both institutes was recorded (Tables 1 and 2). GS for the cohort was categorized using International Society of Urological Pathology and was divided in group 1 (GS 3+3), Group 2 and 3 (GS 7), group 4 (GS 8), and group 5 (GS 9 and 10; ref. 22). Patients’ pretreatment risk status was determined by NCCN and was categorized as low (PSA ≤ 10 ng/mL and GS 3+3 and clinical T stage ≤ 2A), intermediate (PSA >10–20 ng/mL or GS 7, or clinical T stage 2b–2c), and high risk (PSA > 20 ng/mL or GS ≥ 8 or clinical T stage > 3b; ref. 23). In addition, MCC cohort had additional information on patients’ marital, insurance, and history tobacco exposure status. Date of diagnosis was used as an index date to determine time to treatment (delays between diagnosis and prostatectomy) and was calculated as months elapsed since the date of diagnosis to the date of prostatectomy.
. | Time to treatment in months . | |||
---|---|---|---|---|
Variables . | 0–3 (N = 1,111) . | >3–6 (N = 562) . | >6 (N = 134) . | P value . |
Age category | ||||
≤50 | 126 (11.3) | 56 (10.0) | 14 (10.4) | 0.6 |
>50–60 | 505 (45.4) | 243 (43.2) | 53 (39.5) | |
>60–70 | 450 (40.5) | 245 (43.6) | 64 (47.8) | |
≥70 | 30 (2.7) | 18 (3.2) | 3 (2.2) | |
Age | ||||
Median | 59 | 60 | 60 | 0.09 |
IQR | 54–64 | 55–65 | 55–65 | |
Race | ||||
AAM | 149 (13.4) | 111 (19.7) | 34 (25.4) | <0.001 |
White | 962 (86.6) | 451 (80.2) | 100 (74.6) | |
NCCN risk group | ||||
Low risk | 529 (47.6) | 312 (55.5) | 73 (54.5) | 0.03 |
Intermediate risk | 385 (34.6) | 173 (30.8) | 42 (31.3) | |
High risk | 108 (9.7) | 34 (6.0) | 11 (8.2) | |
Unknown | 89 (8.0) | 43 (7.6) | 8 (6.0) | |
iPSA category | ||||
0–6 | 639 (57.5) | 315 (56.0) | 70 (52.2) | 0.4 |
>6–10 | 288 (25.9) | 166 (29.5) | 40 (29.8) | |
>10–20 | 145 (13.0) | 66 (11.7) | 17 (12.7) | |
>20 | 39 (3.5) | 15 (2.7) | 7 (5.2) | |
Clinical GS | ||||
Group 1 | 661 (59.5) | 371 (66.0) | 83 (61.9) | 0.2 |
Group 2/3 | 269 (24.2) | 119 (21.2) | 30 (22.4) | |
Group 4 | 59 (5.3) | 17 (3.0) | 8 (6.0) | |
Group 5 | 13 (1.2) | 3 (0.5) | 1 (0.7) | |
Unknown | 109 (9.8) | 52 (9.2) | 12 (8.9) | |
Clinical T stage | ||||
T1 | 699 (62.9) | 362 (64.4) | 81 (60.4) | 0.5 |
T2 | 238 (21.4) | 120 (21.3) | 34 (25.4) | |
T3 | 16 (1.4) | 4 (0.7) | 0 (0) | |
Unknown | 158 (14.2) | 76 (13.5) | 19 (14.2) | |
Pathology Gleason | ||||
Group 1 | 498 (44.8) | 278 (49.5) | 62 (46.3) | 0.4 |
Group 2/3 | 535 (48.1) | 253 (45.0) | 61 (45.5) | |
Group 4 | 46 (4.1) | 18 (3.2) | 8 (6.0) | |
Group 5 | 32 (2.9) | 13 (2.3) | 13 (2.2) | |
Extracapsular extension | ||||
Present | 290 (26.1) | 134 (23.8) | 35 (26.1) | 0.5 |
Absent | 821 (73.9) | 428 (76.2) | 99 (73.9) | |
Surgical margins | ||||
Present | 257 (23.1) | 115 (20.5) | 34 (25.4) | 0.3 |
Absent | 854 (76.9) | 447 (79.5) | 100 (74.6) | |
Seminal vesicle invasion | ||||
Yes | 63 (5.7) | 27 (4.8) | 9 (6.7) | 0.6 |
No | 1,048 (94.3) | 535 (95.2) | 125 (93.3) | |
AJCC pathologic stage | ||||
pT2A–pT2B | 810 (72.9) | 424 (75.4) | 97 (72.4) | 0.3 |
pT3A | 250 (22.5) | 108 (19.2) | 27 (20.1) | |
pT3B–pT4 | 39 (3.5) | 24 (4.3) | 6 (4.5) | |
PX | 12 (1.1) | 6 (1.1) | 4 (3.0) | |
Gleason upgrade | ||||
Upgraded | 292 (26.3) | 157 (27.9) | 36 (26.9) | 0.8 |
No upgrade | 711 (64.0) | 357 (63.5) | 88 (65.7) | |
Not determined | 108 (9.7) | 48 (8.5) | 10 (7.5) |
. | Time to treatment in months . | |||
---|---|---|---|---|
Variables . | 0–3 (N = 1,111) . | >3–6 (N = 562) . | >6 (N = 134) . | P value . |
Age category | ||||
≤50 | 126 (11.3) | 56 (10.0) | 14 (10.4) | 0.6 |
>50–60 | 505 (45.4) | 243 (43.2) | 53 (39.5) | |
>60–70 | 450 (40.5) | 245 (43.6) | 64 (47.8) | |
≥70 | 30 (2.7) | 18 (3.2) | 3 (2.2) | |
Age | ||||
Median | 59 | 60 | 60 | 0.09 |
IQR | 54–64 | 55–65 | 55–65 | |
Race | ||||
AAM | 149 (13.4) | 111 (19.7) | 34 (25.4) | <0.001 |
White | 962 (86.6) | 451 (80.2) | 100 (74.6) | |
NCCN risk group | ||||
Low risk | 529 (47.6) | 312 (55.5) | 73 (54.5) | 0.03 |
Intermediate risk | 385 (34.6) | 173 (30.8) | 42 (31.3) | |
High risk | 108 (9.7) | 34 (6.0) | 11 (8.2) | |
Unknown | 89 (8.0) | 43 (7.6) | 8 (6.0) | |
iPSA category | ||||
0–6 | 639 (57.5) | 315 (56.0) | 70 (52.2) | 0.4 |
>6–10 | 288 (25.9) | 166 (29.5) | 40 (29.8) | |
>10–20 | 145 (13.0) | 66 (11.7) | 17 (12.7) | |
>20 | 39 (3.5) | 15 (2.7) | 7 (5.2) | |
Clinical GS | ||||
Group 1 | 661 (59.5) | 371 (66.0) | 83 (61.9) | 0.2 |
Group 2/3 | 269 (24.2) | 119 (21.2) | 30 (22.4) | |
Group 4 | 59 (5.3) | 17 (3.0) | 8 (6.0) | |
Group 5 | 13 (1.2) | 3 (0.5) | 1 (0.7) | |
Unknown | 109 (9.8) | 52 (9.2) | 12 (8.9) | |
Clinical T stage | ||||
T1 | 699 (62.9) | 362 (64.4) | 81 (60.4) | 0.5 |
T2 | 238 (21.4) | 120 (21.3) | 34 (25.4) | |
T3 | 16 (1.4) | 4 (0.7) | 0 (0) | |
Unknown | 158 (14.2) | 76 (13.5) | 19 (14.2) | |
Pathology Gleason | ||||
Group 1 | 498 (44.8) | 278 (49.5) | 62 (46.3) | 0.4 |
Group 2/3 | 535 (48.1) | 253 (45.0) | 61 (45.5) | |
Group 4 | 46 (4.1) | 18 (3.2) | 8 (6.0) | |
Group 5 | 32 (2.9) | 13 (2.3) | 13 (2.2) | |
Extracapsular extension | ||||
Present | 290 (26.1) | 134 (23.8) | 35 (26.1) | 0.5 |
Absent | 821 (73.9) | 428 (76.2) | 99 (73.9) | |
Surgical margins | ||||
Present | 257 (23.1) | 115 (20.5) | 34 (25.4) | 0.3 |
Absent | 854 (76.9) | 447 (79.5) | 100 (74.6) | |
Seminal vesicle invasion | ||||
Yes | 63 (5.7) | 27 (4.8) | 9 (6.7) | 0.6 |
No | 1,048 (94.3) | 535 (95.2) | 125 (93.3) | |
AJCC pathologic stage | ||||
pT2A–pT2B | 810 (72.9) | 424 (75.4) | 97 (72.4) | 0.3 |
pT3A | 250 (22.5) | 108 (19.2) | 27 (20.1) | |
pT3B–pT4 | 39 (3.5) | 24 (4.3) | 6 (4.5) | |
PX | 12 (1.1) | 6 (1.1) | 4 (3.0) | |
Gleason upgrade | ||||
Upgraded | 292 (26.3) | 157 (27.9) | 36 (26.9) | 0.8 |
No upgrade | 711 (64.0) | 357 (63.5) | 88 (65.7) | |
Not determined | 108 (9.7) | 48 (8.5) | 10 (7.5) |
Abbreviations: AJCC, American Joint Committee on Cancer; iPSA, preoperative PSA.
Follow-up and outcome
Biochemical disease control was used as the primary outcome of the study and was defined as biochemical recurrence (BCR) after prostatectomy. To capture BCR, patients were followed after their prostatectomy for postsurgical PSA. BCR was determined by calculating the time between prostatectomy and PSA failure. PSA failure was determined by clinician-documented single PSA ≥ 0.2 ng/mL or two consecutive PSA values of 0.2 ng/mL after prostatectomy (23, 24).
Statistical analysis
In order to identify the optimal cutoff; time to treatment from the date of diagnosis and prostatectomy was introduced as a continuous variable in a multivariate cox regression model. Restricted cubic spline (RCS) function with four knots was used in cox model to allow nonlinear association between HR and time to treatment (2, 25). RCS-based cox model was adjusted using all the pertinent pretreatment clinicodemographic variables including clinical GS, PSA, clinical T stage, race, age, and year of prostatectomy. RCS in cox regression allowed us to identify not only the time point at which the risk of BCR increased (Fig. 1) but also the precise subcategorization of the different time interval for the treatment start time for further comparison. Differences between final subcategories of time to treatment and clinicopathologic/demographic characteristics were analyzed using the χ2 test. Methods of survival analysis including Kaplan–Meier (KM) analysis and cox proportional hazard (CPH) model were used to analyze the association between time to treatment subcategories and biochemical disease control. Both univariate and multivariate CPH models were used to estimate the risk of BCR. In the models, cGS, PSA, time to treatment, and race category were used as categorical variables, whereas age at diagnosis and year of prostatectomy were used as numeric variables. A P value of ≤ 0.05 was considered as statistically significant to make clinically meaningful inference from the analysis. Analysis was completed using SAS 9.4 (26).
Results
Comparison of baseline characteristics from the two institutions is provided in Table 1. Median follow-up time of the study cohort was 46 [interquartile range (IQR), 18–86] months, with an overall median time to treatment of 3 months. In multivariate cox regression model with RCS function, adjusted for all the pretreatment variables (Fig. 1), optimal cutoffs were provided for the time to treatment subcategories. Based on RCS model (Fig. 1), around 3 to 4 months after the date of diagnosis, inflection in the risk of BCR was observed (HR started to increase) which further crossed the reference line (HR, 1) near 6 months. Using this model-based approach, time to treatment was divided as 0 to 3, >3 to 6, and >6 months for further comparisons. A large proportion (61.5 %) of the patients received prostatectomy in 0 to 3 months, followed by 31.1% in >3 to 6 months, and 7.4% after 6 months. Higher proportions of African-American men (AAM) were likely to receive treatment after >6 months (Table 2). Delays beyond 6 months were also more prominent among patients with low-risk NCCN group. There was no apparent association among time to treatment categories and other clinicopathologic characteristics (Table 2).
Overall, 5-year freedom from biochemical failure (FFbF) rates were 78% [95% confidence interval (CI), 75%–81%], 82% (95% CI, 78%–85%), and 69% (95% CI, 59%–77%) for time to treatment of 0 to 3, >3 to 6, and >6 months, respectively (P < 0.001; Fig. 2A). Similarly, the trends in unfavorable BCR rate within the strata of time-to-treatment categories were also preserved when KM curves were stratified by NCCN especially among intermediate risk (Fig. 2B–D). In adjusted cox model, treatment delays after 6 months were significantly associated with an increased risk of BCR compared with those with 0 to 3 months (HR, 1.84; 95% CI, 1.30–2.60; P = < 0.001). Conversely, time to treatment >3 to 6 months did not show a significant effect on BCR. Both unadjusted and adjusted estimates for the risk of BCR are provided in Table 3. To explain the possible heterogeneity within the strata of GS, PSA, race, and NCCN risk groupings, the cross-product interaction between these covariates and time to treatment subcategories was analyzed separately. There was no interaction (likelihood ratio P > 0.05) observed between these covariates and time to treatment (result not presented). We also ran a sensitivity analysis in the MCC cohort for the multivariate cox model (Supplementary Table S1) to test whether inclusion of sociodemographic variables such as marital status, health insurance status, and history of tobacco exposure can affect the association between risk of BCR and time to treatment. Inclusion of these variable to the multivariate model adjusted for pretreatment clinical variables does not affect the hazard estimates as treatment delays beyond 6 month continue to negatively affect BCR (HR, 2.41; 95% CI, 1.39–4.17; P = 0.001).
. | . | . | Univariate model . | . | Multivariate modela . | . |
---|---|---|---|---|---|---|
Variables . | Number of patients . | Number of events . | HR (95% CI) . | P value . | HR (95% CI) . | P value . |
Time to treatment (in months) | ||||||
0–3 | 1,111 | 242 | 1 (Ref.) | 1 (Ref.) | ||
>3–6 | 562 | 99 | 0.83 (0.66–1.05) | 0.1 | 0.90 (0.71–1.14) | 0.4 |
>6 | 134 | 39 | 1.70 (1.21–2.38) | 0.002 | 1.84 (1.30–2.60) | <0.001 |
Race | ||||||
White | 1,513 | 311 | 1 (Ref.) | 1 (Ref.) | ||
AAM | 294 | 69 | 1.10 (0.84–1.43) | 0.4 | 1.04 (0.79–1.38) | 0.7 |
Clinical GS | ||||||
Group 1 | 1,115 | 164 | 1 (Ref.) | 1 (Ref.) | ||
Group 2/3 | 418 | 117 | 2.14 (1.68–2.71) | <0.001 | 2.19 (1.71–2.82) | <0.001 |
Group 4 | 84 | 42 | 4.41 (3.13–6.19) | <0.001 | 3.58 (2.50–5.11) | <0.001 |
Group 5 | 17 | 10 | 4.67 (2.46–8.85) | <0.001 | 3.57 (1.84–6.90) | <0.001 |
iPSA category | ||||||
0–6 | 1,024 | 146 | 1 (Ref.) | 1 (Ref.) | ||
>6–10 | 494 | 121 | 1.76 (1.38–2.24) | <0.001 | 1.61 (1.25–2.06) | <0.001 |
>10–20 | 228 | 82 | 3.10 (2.36–4.07) | <0.001 | 2.91 (2.19–3.84) | <0.001 |
>20 | 61 | 31 | 3.89 (2.64–5.74) | <0.001 | 2.56 (1.67–3.91) | <0.001 |
Clinical T stage | ||||||
T1 | 1,142 | 204 | 1 (Ref.) | 1 (Ref.) | ||
T2 | 392 | 112 | 1.49 (1.18–1.87) | <0.001 | 1.28 (1.0–1.64) | 0.04 |
T3 | 20 | 13 | 5.31 (3.19–8.84) | <0.001 | 3.20 (1.83–5.57) | <0.001 |
Age | 1,807 | 48 | 1.0 (0.99–1.02) | 0.2 | 0.99 (0.97–1.0) | 0.2 |
RP year | 1,807 | 48 | 0.96 (0.94–0.98) | <0.001 | 0.97 (0.95–1.0) | 0.1 |
. | . | . | Univariate model . | . | Multivariate modela . | . |
---|---|---|---|---|---|---|
Variables . | Number of patients . | Number of events . | HR (95% CI) . | P value . | HR (95% CI) . | P value . |
Time to treatment (in months) | ||||||
0–3 | 1,111 | 242 | 1 (Ref.) | 1 (Ref.) | ||
>3–6 | 562 | 99 | 0.83 (0.66–1.05) | 0.1 | 0.90 (0.71–1.14) | 0.4 |
>6 | 134 | 39 | 1.70 (1.21–2.38) | 0.002 | 1.84 (1.30–2.60) | <0.001 |
Race | ||||||
White | 1,513 | 311 | 1 (Ref.) | 1 (Ref.) | ||
AAM | 294 | 69 | 1.10 (0.84–1.43) | 0.4 | 1.04 (0.79–1.38) | 0.7 |
Clinical GS | ||||||
Group 1 | 1,115 | 164 | 1 (Ref.) | 1 (Ref.) | ||
Group 2/3 | 418 | 117 | 2.14 (1.68–2.71) | <0.001 | 2.19 (1.71–2.82) | <0.001 |
Group 4 | 84 | 42 | 4.41 (3.13–6.19) | <0.001 | 3.58 (2.50–5.11) | <0.001 |
Group 5 | 17 | 10 | 4.67 (2.46–8.85) | <0.001 | 3.57 (1.84–6.90) | <0.001 |
iPSA category | ||||||
0–6 | 1,024 | 146 | 1 (Ref.) | 1 (Ref.) | ||
>6–10 | 494 | 121 | 1.76 (1.38–2.24) | <0.001 | 1.61 (1.25–2.06) | <0.001 |
>10–20 | 228 | 82 | 3.10 (2.36–4.07) | <0.001 | 2.91 (2.19–3.84) | <0.001 |
>20 | 61 | 31 | 3.89 (2.64–5.74) | <0.001 | 2.56 (1.67–3.91) | <0.001 |
Clinical T stage | ||||||
T1 | 1,142 | 204 | 1 (Ref.) | 1 (Ref.) | ||
T2 | 392 | 112 | 1.49 (1.18–1.87) | <0.001 | 1.28 (1.0–1.64) | 0.04 |
T3 | 20 | 13 | 5.31 (3.19–8.84) | <0.001 | 3.20 (1.83–5.57) | <0.001 |
Age | 1,807 | 48 | 1.0 (0.99–1.02) | 0.2 | 0.99 (0.97–1.0) | 0.2 |
RP year | 1,807 | 48 | 0.96 (0.94–0.98) | <0.001 | 0.97 (0.95–1.0) | 0.1 |
Abbreviation: iPSA, preoperative PSA.
aMultivariate model was also adjusted for the institutions used in the study (UPHS and MCC) to account for the institution-specific effects.
Discussion
Time to treatment and biochemical disease control
This study optimizes the ideal time to treatment by applying model-based approach and assesses the implications of treatment delays on prostate cancer outcome. In our analysis, we carefully included all the pretreatment or clinical variables in multivariate RCS cox model, that are available to clinicians in order to make treatment decision, to identify the time point beyond which the risk of BCR was higher. Based on our results, delays in treatment after diagnosis for up to 6 months may not adversely affect prostate cancer outcome; however, significant delays beyond 6 months may lead to poor biochemical disease control. These results were consistent with another retrospective study by O’Brien and colleagues who reported a significantly higher risk of BCR among patients with surgical delays beyond 6 months. It is important to note that in our RCS-based model, risk of BCR increased around the same time point as O’Brien and colleagues reported (9). Interestingly, when stratified by NCCN risk groups, rate of BCR was significantly higher among intermediate-risk patients, as there was no statistical difference observed among low-risk patients. Our finding aligns perfectly with a large retrospective study by Abern and colleagues, who reported a high risk of BCR among intermediate-risk patients with delayed prostatectomy, whereas no apparent risk of poor outcomes among low-risk patients (10). Given that higher proportions of low-risk patients were likely to wait for surgery beyond 6 months, observed nonsignificant effect of treatment delays on BCR can be reassuring. However, attempts to perform restaging after 6 months may be necessary to avoid any disease upgrade as delays negatively affect BCR.
Contrary to our findings, there are studies which argue against the association of treatment delays and its implications on prostate cancer outcomes. In a study by Korest and colleagues (14), the authors reported negligible effect of treatment delays on BCR-free survival and prostate cancer pathology. In addition, a review by Van Den Bergh and colleagues reported nonsignificant effect of treatment delays on prostate cancer outcomes and reported that delays up to several years may not inferiorly affect prostate cancer outcome (15). In another study by Morini and colleagues on 908 prostate cancer patients, authors reported nonsignificant association of treatment delays and BCR among low-/intermediate-risk patients (16). However, the study design used in most of these contemporary studies utilized arbitrary or a priori–determined cutoff points to categorize time to treatment and then used these categories to predict prostate cancer outcome (14, 15, 18–21). This approach is susceptible to misclassification bias as it ignores the important factors such as clinical characteristics which have a direct association with treatment planning and may affect time to treatment (2, 10, 27). Furthermore, most studies often fail to acknowledge that prostate cancer progression is not fixed and alterations in the tumor characteristics as a result of biological changes are inevitable due to lack of definite treatment. Vickers and colleagues draw attention to this issue and point to the frequent use of null hypothesis–based approach in reporting that delayed treatment does not lead to any change in prostate cancer pathology (17), which can be misleading because it does not render any clinical utility toward treatment recommendation. While analyzing our data, we addressed this issue differently and applied a model-based approach to optimize time to treatment and provided an optimal time point beyond which risk of BCR starts to increase. In addition, with the inclusion of all the key clinical variables that are used in prostate cancer treatment decision-making, we further increased the clinical relevance of our estimates to better guide prostate cancer treatment planning. In addition, we included pooled data from two large cancer institutes and had a significantly large sample size with higher statistical power to assess the association treatment delays and prostate cancer.
Clinical implications of time to treatment
Our results will aid in the shared decision-making process of prostate cancer treatment for both clinicians and patients. Although the cohort we used in analysis consisted of nonactive surveillance patients, our results may still have implications on the existing practices related to repeat biopsy especially among those who are likely to experience treatment delays. Currently, NCCN guidelines identify active surveillance as viable strategy for very low, low-risk, and a small subset of favorable intermediate-risk patients and recommend repeat biopsy in 12 months (23, 28). However, observed outcome deficits with increasing time to treatment may necessitate the restaging before 12 months especially among intermediate-risk patients. Although the absence of any statistical association of treatment delays and BCR among low-risk patients can be reassuring, attempts to avoid undue delays must be ensured to obviate any possibility of losing window of cure and prevent inferior outcomes. Our results underscore the clinical utility of time to treatment in targeting the subset of patients who are likely to benefit from initiation of treatment soon after diagnosis. This result necessitates the inclusion of time to treatment in the decision-making of prostate cancer treatment recommendation. Although this study does not corroborate whether Gleason progression or pathologic upstaging could have resulted in poor outcome among those with treatment delays beyond 6 months, future studies should focus on understanding the interaction between prostate cancer pathology and treatment delays to explain outcomes’ heterogeneity with increasing time to treatment.
The presence of detailed clinicopathologic information along with longer follow-up from two large tertiary care referral centers adds significant value to the generalizability of our results. In addition, patients across the different surgical era (1987–2015) were included in the analysis, therefore allowing us to capture data spanning over 25 years. Compared with other studies evaluating similar associations (14, 15, 18–21), our study consisted of a higher percentage of a traditionally underrepresented population of AAM (∼16% of the study population), thus making our results more generalizable to AAM who are at risk of experiencing inferior outcomes. In our analysis, application of the model-based approach over arbitrary cut point selection methods, to select optimal time point for time to treatment, reduced the likelihood of selection biases among those with delayed treatment. This study also has few limitations. First, our cohort only consisted of prostatectomy patients, thus lacking the information on those who received or contemplated receiving radiotherapy as a primary treatment. Therefore, our results may lack the generalizability across the treatment spectrum of prostate cancer, and further investigation may be warranted to confirm our finding in patients that received radiotherapy as their primary treatment. Second, in order to significantly increase our sample size and provide generalizability of our findings, the analysis for this study was conducted using combined data from two different tertiary referral centers. Although we controlled for the health centers’ effects in multivariate analysis, unmeasured confounding resulting from the differential treatment practices and diagnostic ascertainment (lack of standardization in pathologist review across the patient population) remains a potential limitation. It is also important to note that UPHS was relatively older institution and contributed a large proportion of older cases in our cohort (diagnosed before 2000). Because PSA screening also began around the same time, treatment delays in the older cases were more prominent among patients from UPHS than MCC. This variation can explain the higher proportion of pT3b-pT4 disease (advance stage) in the UPHS cohort (Table 1). Furthermore, loss to follow-up was another limitation of this study. Likelihood that many patients after receiving their primary treatment at one of the two tertiary centers (MCC and UPHS) would have moved on to seek care in their local facility remains higher and can result in inadequate follow-up. Compared with other categories, there were fewer patients with treatment delays beyond 6 months. Smaller numbers in this category may reduce the precision in the observed risk estimate and their corresponding confidence interval. Lastly, BCR was the only surrogate used in the analysis compared with metastasis and prostate cancer–specific mortality. Use of BCR as a surrogate may potentially limit our ability to associate our finding with long-term survival endpoints. As highlighted by McLaughlin and colleagues (3), ethical issues remain one of the major limitations in understanding the optimal time to initiate treatment, resulting in a lack of interventional prospective studies to determine optimal time to treatment. Continuing research will be required to advance existing guidelines on the delivery of prostate cancer treatment.
Conclusions
Time to treatment >6 months has detrimental effects on biochemical disease control after surgery. In order to achieve sustainable outcomes, definitive treatment for prostate cancer must be delivered within 6 months. This study identifies time to treatment as an independent predictor of prostate cancer outcome, and therefore, careful planning may be warranted while considering treatment delays.
Disclosure of Potential Conflicts of Interest
D.I. Lee received honoraria from the speakers’ bureau of Intuitive Surgical. No potential conflicts of interest were disclosed by the other authors.
Authors' Contributions
Conception and design: S. Awasthi, R. Balkrishnan, P. Lal, T.R. Rebbeck, K. Yamoah
Development of methodology: S. Awasthi, F.A. Asamoah, T.R. Rebbeck, K. Yamoah
Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): S. Awasthi, D.I. Lee, S.B. Malkowicz, P. Lal, T.R. Rebbeck, K. Yamoah
Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): S. Awasthi, T. Gerke, V.L. Williams, S.B. Malkowicz, K. Yamoah
Writing, review, and/or revision of the manuscript: S. Awasthi, T. Gerke, J.Y. Park, F.A. Asamoah, V.L. Williams, A.K. Fink, R. Balkrishnan, D.I. Lee, J. Dhillon, J.M. Pow-Sang, T.R. Rebbeck, K. Yamoah
Administrative, technical, or material support (i.e., reporting or organizing data, constructing databases): S. Awasthi, T. Gerke, K. Yamoah
Study supervision: R. Balkrishnan, K. Yamoah
Others (pathology review and database management): P. Lal
Acknowledgments
This work was supported by the Prostate Cancer Foundation Young Investigator Award (to K. Yamoah), American Cancer Society MRSG-17-108-01-TBG (to K. Yamoah), and V Foundation (to K. Yamoah and J.Y. Park).
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.