Development and validation of a preoperative nomogram for predicting survival of patients with locally advanced prostate cancer after radical prostatectomy
BMC Cancer volume 20, Article number: 97 (2020)
For selected locally advanced prostate cancer (PCa) patients, radical prostatectomy (RP) is one of the first-line treatments. We aimed to develop a preoperative nomogram to identify what kinds of patients can get the most survival benefits after RP.
We conducted analyses with data from the Surveillance, Epidemiology, and End Results (SEER) database. Covariates used for analyses included age at diagnosis, marital status, race, American Joint Committee on Cancer (AJCC) 7th TNM stage, Prostate specific antigen, Gleason biopsy score (GS), percent of positive cores. We estimated the cumulative incidence function for cause-specific death. The Fine and Gray’s proportional subdistribution hazard approach was used to perform multivariable competing risk analyses and reveal prognostic factors. A nomogram was built by these factors (including GS, percent of positive cores and N stage) and validated by concordance index and calibration curves. Risk stratification was established based on the nomogram.
We studied 14,185 patients. N stage, GS, and percent of positive cores were the independent prognostic factors used to construct the nomogram. For validating, in the training cohort, the C-index was 0.779 (95% CI 0.736–0.822), and in the validation cohort, the C-index was 0.773 (95% CI 0.710–0.836). Calibration curves showed that the predicted survival and actual survival were very close. The nomogram performed better over the AJCC staging system (C-index 0.779 versus 0.764 for training cohort, and 0.773 versus 0.744 for validation cohort). The new stratification of risk groups based on the nomogram also showed better discrimination than the AJCC staging system.
The preoperative nomogram can provide favorable prognosis stratification ability to help clinicians identify patients who are suitable for surgery.
Prostate cancer (PCa) is one of the most common malignant tumors of males around the world, and it is estimated that there will be 174,650 new male cases in the United States in 2019 . The European Association of Urology defined PCa with cT3–4 or N+ as locally advanced PCa . Patients with locally advanced PCa have an increased risk of disease progression and cancer-specific death . Due to many factors like the American national guidelines advising against PSA screening, increasing applications for active surveillance, etc., the proportion of men diagnosed with locally advanced PCa has increased [3, 4]. Nowadays, the standard treatment of locally advanced PCa is still unclear due to the absence of level 1 evidence. For selected locally advanced PCa patients, radical prostatectomy (RP) is one of the first-line treatments recommended by the guidelines [2, 5]. However, it is still controversial as to which types of patients could get the most survival benefits from RP.
Some studies about the risk stratification of patients with high-risk PCa (including locally advanced PCa) after RP had been reported. Previous researchers performed a retrospective analysis of 315 high-risk PCa patients after RP in their hospital. They selected biochemical progression as the primary endpoint, and reported that the risk of biochemical progression of high-risk PCa after RP could be stratified by Gleason Score (GS) at biopsy (≥ 8) and % positive core (≥ 30%) . Another research team studied 813 high-risk patients undergoing RP, and found 3 preoperative criteria including stage cT2 or greater, prostate specific antigen (PSA) > 20 ng/mL, and GS > 8. In their study, the number of meeting preoperative criteria was significantly predictive for recurrence-free survival (RFS), and overall survival (OS) . Although these studies had proposed their own risk stratification models for high-risk PCa patients undergone RP, the patients involved in these studies were not fully compliant with the current more common definition of locally advanced PCa, the impact of prognostic factors on patient survival outcome was not quantitative, the weight between the various prognostic factors was not clear enough and the cohorts were not large enough.
To circumvent these defects, we evaluated the prostate cancer-specific survival (CSS) at a large cohort to assess potential preoperative prognostic factors, and developed a preoperative nomogram. This nomogram could be used to predict the CSS of patients with locally advanced PCa after RP and to identify what kinds of patients can get the most survival benefits from RP.
The Surveillance, Epidemiology, and End Results (SEER) database is a free dataset made up of 18 population-based cancer registries. It releases cancer patients’ general information annually and has almost covered 25% population of the United States .
From the SEER database, patients with a diagnosis of adenocarcinoma of the prostate (International Classification of diseases-O-3 code: C61.9) between 2010 and 2016 were selected. The TNM stages were assessed by the 7th edition of American Joint Committee on Cancer [AJCC] Cancer Staging Manual . Inclusion and exclusion criteria were shown in the flowchart in detail (Fig. 1). All the included patients were randomly divided into the training cohort and validation cohort with a ratio of 7:3.
Variables and endpoint
Data for each patient were extracted from the SEER database including age at diagnosis, marital status, race, AJCC 7th TNM stage, PSA, Gleason biopsy score (GS), percent of positive cores at biopsy (% positive core), and follow-up information. Age was categorized as ≤49 years, 50–59 years, 60–69 years, and ≥ 70 years. PSA was categorized as ≤10 ng/mL, 10–20 ng/mL (not including 10), and > 20 ng/mL. GS was classified as ≤6, 7(3 + 4), 7(4 + 3), 8, and ≥ 9. % positive core was classified as 0–25, 25–50% (not including 25%), 50–75% (not including 50%), and 75–100% (not including 75%). In addition to follow-up information, all clinical factors were obtained before undergoing RP. Due to the limitation of the SEER database, the information about detail procedures used to determine the lymph node positive was not been provided .
CSS was used as the primary endpoint. CSS was measured by all deaths caused by prostate cancer, complications of treatments, or unknown processes in patients with active prostate cancer. Follow-up time was defined as the time between the first treatment and the patient’s death or last follow-up.
A χ 2 test was used to expose the difference between the training cohort and validation cohort in categorical variables, and the results were presented as the frequency with its proportion. Kaplan–Meier method and log-rank test were used to expose each potential prognostic variable. Variables in the univariate analysis with p-value < 0.05 were chosen for multivariate Cox proportional hazard regression to identify the independent prognostic factors. In addition, many patients with locally advanced PCa are the elderly with many comorbidities. Thus, there is a high risk of non-cancer competing risk events leading to death in the long-term survival . In order to make the research more precise, we further implemented competing risk analyses. We estimated the cumulative incidence function (CIF) for cause-specific death, and tested the differences by Gray’s test to expose each prognostic variable. The Fine and Gray’s proportional subdistribution hazard approach was used to conduct competing risk multivariable analyses. Based on clinical significance and p-value < 0.05 in competing risk multivariable analyses, these independent prognostic factors were selected to construct a nomogram to predict 5-year cancer-specific mortality (CSM) probabilities of patients.
The performance of the nomogram was validated by measuring discrimination and calibration in both training set and validation set. The discrimination of the nomogram was assessed by C-index (concordance index). The value of C-index ranged from 0.5 (no discrimination) to 1 (perfect discrimination), and higher C-index showed better discrimination of the prognostic model . By comparing the plot of predicted probabilities from the nomogram with the plot of actual probabilities of CSS, the calibration curve was used to compare the predicted survival using nomogram with the actual observed survival, with 1000 resamples of bootstrapping. In addition, we conducted a risk stratification of the nomogram total risk score. The cut-off values were determined using X-tile software and by the minimal p-value approach . The discrimination abilities between the nomogram and AJCC staging system were compared by C-index and Kaplan–Meier curve.
The statistical software packages R (The R Foundation) and R packages cmprsk, rms, survival, mstate, and regplot were used in the above statistical analyses. A p-value < 0.05 was considered statistically significant.
A total of 14,185 eligible patients were included in the current analysis as the primary cohort. Among all patients, 9932 patients were placed within the training cohort, and 4253 patients were placed within the validation cohort randomly. Detail patients’ characteristics were shown in Table 1. The median follow-up was 43 months (95% CI 42–44 months). There were no statistically significant differences in baseline characteristics between the training cohort and the validation cohort.
Determination of independent preoperative prognostic factors
The Fine and Gray’s proportional subdistribution hazard approach was performed to demonstrate the independent preoperative prognostic factors in the training cohort. Table 2 showed the detail univariate and multivariate analysis results. In the univariate analyses, age, N stage, PSA, GS, % positive core biopsy were variables with significant impact on CSS.
In the Fine and Gray’s proportional subdistribution hazard approach, N stage, GS, % positive core biopsy remained significantly, while age and PSA did not show significant impact on CSS. These significant variables were thought as independent prognostic factors of CSS for locally advanced PCa patients after RP.
Construction and validation of the nomogram
Using the independent prognostic factors including N stage, GS, and %positive core biopsy, a nomogram was constructed to predict each locally advanced PCa patient’s probability of CSM at 5 years after RP (Fig. 2).
The nomograms were validated using C-index and calibration curve. In the training cohort, the C-index was 0.779 (95% CI 0.736–0.822). In the validation cohort, the C-index was 0.773 (95% CI 0.710–0.836). The C-index values indicated good discrimination of the nomogram (> 0.7). The calibration curves showed good agreement between prediction by nomogram and actual 5-year CSS in both training cohort and validation cohort (Fig. 3).
Definition of nomogram risk group stratification
The score corresponding to each nomogram variable was listed in detail in Table 3. By summing the score for each nomogram variable, we got the nomogram total risk score for each patient in both the training cohort and the validation cohort. The patients were divided into three nomogram risk groups for CSS: a low-risk group with 0–77 points, a middle-risk group with 79–108 points, a high-risk group with no less than 112 points.
Comparison of nomogram with AJCC staging system
We compared the nomogram with the 8th edition AJCC staging system. In the training cohort, the C-index of AJCC staging system was 0.764 (95% CI 0.719–0.809), and the C-index of the nomogram was 0.779 (95% CI 0.736–0.822). In the validation cohort, the C-index of AJCC staging system was 0.744 (95% CI 0.707–0.817), and the C-index of the nomogram was 0.773 (95% CI 0.710–0.836). In both cohorts, the nomogram performed better in C-index, indicating that the nomogram had a better discrimination ability than AJCC staging system.
The AJCC staging could divide the whole eligible patients in current study into 3 groups including IIIB, IIIC, and IVA. Relatively, as mentioned above, the risk groups based on the nomogram included low-risk group, middle-risk group, high-risk group. For risk groups based on the nomogram, the degree of separation of CIF curves of CSM between groups was more obvious than AJCC staging system in both training cohort and validation cohort (Fig. 4). The CIF curves still demonstrated the new stratification of risk groups based on the nomogram had better prognostic discrimination than AJCC staging system.
Using a large cohort of 14,185 patients from the SEER database, we successfully used 3 independent preoperative prognostic factors to established a nomogram for predicting 5-year CSM for patients with locally advanced PCa. To our knowledge, this is the first preoperative predictive nomogram and nomogram-based risk group stratification built for patients with locally advanced PCa. At the same time, the patient samples are the latest relatively (2010–2016), which is conducive to reducing the impact of defects in pathology, surgery and laboratory testing techniques on the results. This nomogram and nomogram-based risk group stratification is not inferior or even better than the current AJCC staging system in terms of discrimination, and is more quantitative and intuitive, which is convenient for clinicians to use.
Nomograms use a variety of biological and clinical variables to graphically depict the probability that a clinical event may occur for each individual . One of their unique abilities is to estimate individualized risk based on patient and disease characteristics in a graphical and user-friendly form. Nomograms can include many continuous variables and key factors of disease into the prognosis, and taking weights of each variable into account, to make predictive models more practical [14, 15]. Compared with the current prognostic prediction system like AJCC staging system, the nomogram showed advantages in many studies on different cancers [16, 17]. They can also establish individualized risk stratification and help clinicians identify suitable patients for optimal managements, and good results had been obtained in many studies of different tumor types [18,19,20].
In our study, the statistical analysis highlighted 3 preoperative prognostic factors including clinical N stage, GS, and % positive core. Among them, GS had the greatest contribution to 5-year cancer-specific death (GS shared the most proportion of the nomogram total risk score.). GS is one of the most universal prognostic factors for PCa. International Society of Urological Pathology (ISUP) proposed in 2014, GS can be divided into five groups (2–6, 7(3 + 4), 7(4 + 3), 8, ≥9) according to the differences of prognosis . In our study, the classification of GS also referred to the ISUP grading system, and the risk trend shown by the nomogram was consistent with this grading system. Some previous studies had shown the relationship between the prognosis of PCa and GS [6, 22, 23]. These studies had confirmed the increased possibility of adverse clinical events including death and biochemical recurrence in PCa patients with increasing GS. But these studies did not show what the weight of GS compared to other variables on affecting patients was. % positive core showed the obvious impact on CSS for patients in the current study, only inferior to GS. % positive core refers to the proportion of positive needles in the needle biopsy to the total number of needles. Many studies had revealed the prognostic value of % positive core for PCa [6, 22, 24, 25]. In our nomogram, the risk score of group 75–100% increased significantly from 7 (group 50–75%) to 30, suggesting that 75–100% positive core can more significantly affect the patient’s prognosis. This finding was similar to the results of a study of 195 patients with high-risk PCa who undergoing RP. In that study, 70% positive was shown to have significant prognostic discrimination . From the nomogram, regional lymph node metastasis (N1) was another important independent prognostic factor for CSS of locally advanced PCa patients. Previous studies had demonstrated that lymph node metastasis in most patients with PCa was associated with progressive disease, which was important for the prognosis of patients [26, 27]. A retrospective study of 229 patients indicated that the 5-year disease-free survival of rate of patients with regional lymph node metastasis could be decreased from 85% to approximately 50% .
Relatively, the 8th AJCC staging system is currently a general PCa risk stratification system, not aim at specific disease stages and treatment therapies . We compared our nomogram risk group stratification with AJCC staging system, and the results showed that our stratification had advantages in C-index (0.779 versus 0.763 for training cohort, and 0.773 versus 0.745 for validation cohort). At the same time, the nomogram risk stratification’s CIF curve showed a better separation trend in both the training cohort and the validation cohort. There may be several reasons. First, our nomogram risk stratification considered the effect of % positive core on CSS of locally advanced PCa patients after RP, compared to AJCC staging system. AJCC staging system only took GS and clinical N stage into account, but this study and other studies mentioned above confirmed the importance of % positive core. Second, AJCC staging system was not quantitative, and the difference between the effects of GS and clinical N stage on prognosis was not reflected. In our nomogram risk stratification, we used different risk scores to reflect the weight of different prognostic factors, showing the obvious importance of GS among the three prognostic factors. At the same time, such quantitative and graphical tool is also easy for clinicians to use as manuals. In addition, the patients included in our study were more targeted than those in AJCC staging system after filtered by locally advanced PCa and RP. By classifying the prognosis of each patient after surgery, clinicians could better choose whether to perform surgery for each patient.
The main significance of the risk stratification proposed by us is to reveal the potential value of the three prognostic factors for the selection of locally advanced patients to receive RP and the influence weight of each factor on the prognosis of these patients. According to the latest guidelines, radical treatments including RP and RT are still the first-line treatments for locally advanced PCa. Only those patients who are unwilling or unable to receive any form of local treatment, or whose PSA growth is fast, or whose tumor differentiation is poor, are offered androgen deprivation therapy (ADT) monotherapy. The guidelines suggest that RP should be used in highly selective locally advanced PCa patients . But at present, there is no high-level evidence on the appropriate selection criteria of RP or RT and the survival benefits comparison between RP and RT for locally advanced PCa patients. In our study, the risk stratification proposed by the nomogram provided a reference for the selection criteria of RP. Considering the difference of postoperative CSS, patients in high-risk group may be not suitable for RP, while patients in low-risk group may be the appropriate population for RP. For those who are not suitable for RP, the RT + ADT therapy may bring more survival benefits. In fact, previous studies have shown that RT + ADT therapy can bring significant benefits to locally advanced PCa patients with pN+ [31, 32]. However, they didn’t compare RP and RT directly. More clinical studies are needed to explore the most reasonable criteria for the selection of locally advanced PCa patients receiving RP or RT, especially the double arm randomized controlled studies of RP and RT, such as the SPCG-15 currently in progress .
There are still several limitations to our study. First, our research is constructed by retrospective data. Therefore, there may be some undetected potential bias factors in the study. Second, due to the SEER database’s limitations, we lack functional data and disease progression data for further research. Third, we also lack additional independent external validation sets, and this is our important work goal in the future.
By analyzing a large cohort of 14,185 patients, we successfully revealed 3 independent preoperative prognostic factors for CSS of locally advanced PCa patients after RP. Using these 3 factors, we constructed and validated a preoperative prognostic nomogram that could predict the probability of CSM at 5 years of this type of patients after undergoing RP. Based on the nomogram, a risk stratification system was established to help clinicians better identify patients suitable for surgery.
Availability of data and materials
All analyzed data are publicly available at the SEER website (http://www.seer.cancer.gov), and should be requested under the approval of the SEER Program administration.
Androgen deprivation therapy
American Joint Committee on Cancer
Cumulative incidence function
Prostate specific antigen
The Surveillance, Epidemiology, and End Results
Siegel RL, Miller KD, Jemal A. Cancer statistics, 2019. CA Cancer J Clin. 2019;1:7–34.
Mottet N, Bellmunt J, Bolla M, et al. EAU-ESTRO-SIOG guidelines on prostate cancer. Part 1: screening, diagnosis, and local treatment with curative intent. Eur Urol. 2017;4:618–29.
Huland H, Graefen M. Changing trends in surgical management of prostate cancer: the end of overtreatment? Eur Urol. 2015;2:175–8.
Fletcher SA, von Landenberg N, Cole AP, et al. Contemporary national trends in prostate cancer risk profile at diagnosis. Prostate Cancer Prostatic Dis. 2019. https://doi.org/10.1038/s41391-019-0157-y.
Mohler JL, Antonarakis ES, Armstrong AJ, et al. Prostate cancer, version 2.2019, NCCN clinical practice guidelines in oncology. J Natl Compr Cancer Netw. 2019;5:479–505.
Nagao K, Matsuyama H, Matsumoto H, et al. Identification of curable high-risk prostate cancer using radical prostatectomy alone: who are the good candidates for undergoing radical prostatectomy among patients with high-risk prostate cancer? Int J Clin Oncol. 2018;4:757–64.
Ploussard G, Masson-Lecomte A, Beauval JB, et al. Radical prostatectomy for high-risk prostate cancer defined by preoperative criteria: oncologic follow-up in national multicenter study in 813 patients and assessment of easy-to-use prognostic substratification. Urology. 2011;3:607–13.
Cronin KA, Ries LA, Edwards BK. The surveillance, epidemiology, and end results (SEER) program of the national cancer institute. Cancer. 2014;120(Suppl 23):3755-7.
Cuccurullo V. AJCC cancer staging handbook: from the AJCC cancer staging manual (7Th edition). Eur J Nucl Med Mol Imaging. 2011;2:408.
Tward JD, Kokeny KE, Shrieve DC. Radiation therapy for clinically node-positive prostate adenocarcinoma is correlated with improved overall and prostate cancer-specific survival. Pract Radiat Oncol. 2013;3:234–40.
Vagnoni V, Bianchi L, Borghesi M, et al. Adverse features and competing risk mortality in patients with high-risk prostate cancer. Clin Genitourin Cancer. 2017;2:e239–48.
Weiss A, Chavez-Macgregor M, Lichtensztajn DY, et al. Validation study of the American joint committee on cancer eighth edition prognostic stage compared with the anatomic stage in breast cancer. JAMA Oncol. 2018;4:203-9.
Camp RL, Marisa DF, Rimm DL. X-tile: a new bio-informatics tool for biomarker assessment and outcome-based cut-point optimization. Clin Cancer Res. 2004;21:7252–9.
Balachandran VP, Gonen M, Smith JJ, DeMatteo RP. Nomograms in oncology: more than meets the eye. Lancet Oncol. 2015;4:e173–80.
Shouval R, Labopin M, Gorin NC, et al. Individualized prediction of leukemia-free survival after autologous stem cell transplantation in acute myeloid leukemia. Cancer. 2019;125:3566-73.
Cho CSM, Gonen MP, Shia JM, et al. A novel prognostic nomogram is more accurate than conventional staging Systems for predicting survival after resection of hepatocellular carcinoma. J Am Coll Surgeons. 2008;2:281–91.
Wong SL, Kattan MW, McMasters KM, Coit DG. A Nomogram that predicts the presence of sentinel node metastasis in melanoma with better discrimination than the American joint committee on cancer staging system. Ann Surg Oncol. 2005;4:282–8.
Wu Y, Meyers JP, Shi G, et al. A nomogram for predicting survival and retroperitoneal lymph node dissection treatment in patients with resected testicular germ cell tumors. J Surg Oncol. 2019;3:508–17.
Kim BH, Kim K, Chie EK, et al. Risk stratification and prognostic nomogram for post-recurrence overall survival in patients with recurrent extrahepatic cholangiocarcinoma. Hpb. 2017;5:421–8.
Kim Y, Park HC, Yoon SM, et al. Prognostic group stratification and nomogram for predicting overall survival in patients who received radiotherapy for abdominal lymph node metastasis from hepatocellular carcinoma: a multi-institutional retrospective study (KROG 15-02). Oncotarget. 2017;55:94450–61.
Epstein JI, Egevad L, Amin MB, Delahunt B, Srigley JR, Humphrey PA. The 2014 International Society of Urological Pathology (ISUP) consensus conference on Gleason grading of prostatic carcinoma: definition of grading patterns and proposal for a new grading system. Am J Surg Pathol. 2016;2:244–52.
Sundi D, Wang V, Pierorazio PM, et al. Identification of men with the highest risk of early disease recurrence after radical prostatectomy. Prostate. 2014;6:628–36.
Kim TH, Jeon HG, Jeong BC, et al. Development of a new nomogram to predict insignificant prostate cancer in patients undergoing radical prostatectomy. Scand J Urol. 2017;1:27–32.
Grossklaus DJ, Coffey CS, Shappell SB, Jack GS, Chang SS, Cookson MS. Percent of cancer in the biopsy set predicts pathological findings after prostatectomy. J Urol. 2002;5:2032–6.
Hamada R, Nakashima J, Ohori M, et al. Preoperative predictive factors and further risk stratification of biochemical recurrence in clinically localized high-risk prostate cancer. Int J Clin Oncol. 2016;3:595–600.
Evangelista L, Guttilla A, Zattoni F, Muzzio PC, Zattoni F. Utility of choline positron emission tomography/computed tomography for lymph node involvement identification in intermediate- to high-risk prostate cancer: a systematic literature review and meta-analysis. Eur Urol. 2013;6:1040–8.
Sheng W, Zhang H, Lu Y. Survival outcomes of locally advanced prostate cancer in patients aged < 50 years after local therapy in the contemporary US population. Int Urol Nephrol. 2018;8:1435–44.
Danella JF, DeKernion JB, Smith RB, Steckel J. The contemporary incidence of lymph node metastases in prostate cancer: implications for laparoscopic lymph node dissection. J Urol. 1993;6:1488–91.
Amin MB, Greene FL, Edge SB, et al. The eighth edition AJCC Cancer staging manual: continuing to build a bridge from a population-based to a more “personalized” approach to cancer staging. CA Cancer J Clin. 2017;2:93–9.
Mottet N, van den Bergh RCN, Briers E, et al. EAU Guidelines. https://uroweb.org/guidelines/2019. Accessed 10 Jan 2020.
Pilepich MV, Winter K, Lawton CA, et al. Androgen suppression adjuvant to definitive radiotherapy in prostate carcinoma--long-term results of phase III RTOG 85-31. Int J Radiat Oncol Biol Phys. 2005;5:1285–90.
James ND, Spears MR, Clarke NW, et al. Failure-free survival and radiotherapy in patients with newly diagnosed nonmetastatic prostate cancer: data from patients in the control arm of the STAMPEDE trial. JAMA Oncol. 2016;3:348–57.
Stranne J, Brasso K, Brennhovd B, et al. SPCG-15: a prospective randomized study comparing primary radical prostatectomy and primary radiotherapy plus androgen deprivation therapy for locally advanced prostate cancer. Scand J Urol. 2018;5–6:313–20.
The authors acknowledge the efforts of the Surveillance Research Program in the establishment and maintenance of the SEER database.
This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.
Ethics approval and consent to participate
This study used a general dataset from the Surveillance, Epidemiology, and End Results database built by a public health program, and therefore did not require institutional review board approval and the patients’ data were atomized.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Zhou, X., Ning, Q., Jin, K. et al. Development and validation of a preoperative nomogram for predicting survival of patients with locally advanced prostate cancer after radical prostatectomy. BMC Cancer 20, 97 (2020). https://doi.org/10.1186/s12885-020-6565-5