Prognostic value of the bone scan index using a computer-aided diagnosis system for bone scans in hormone-naive prostate cancer patients with bone metastases
BMC Cancer volume 16, Article number: 128 (2016)
The bone scan index (BSI) using a computer-aided diagnosis system for bone scans is expected to be an objective and quantitative clinical tool for evaluating bone metastatic prostate cancer. This study aimed to evaluate the pretreatment BSI as a prognostic factor in hormone-naive prostate cancer patients with bone metastases.
The study included 60 patients with hormone-naive, bone metastatic prostate cancer that was initially treated with combined androgen blockade therapy. The BONENAVI system was used for calculating the BSI. We evaluated the correlation between overall survival (OS) and pretreatment clinicopathological characteristics, including patients’ age, initial prostate-specific antigen (PSA) value, Gleason scores, clinical TNM stage, and the BSI. Cox proportional hazards regression models were used for statistical analysis.
The median follow-up duration was 21.4 months. Clinical or PSA progression occurred in 37 (61.7 %) patients and 18 (30.0 %) received docetaxel. Death occurred in 16 (26.7 %) patients. Of these deaths, 15 (25.0 %) were due to prostate cancer. The median OS was not reached. In multivariate analysis, age and the BSI were independent prognostic factors for OS. We evaluated the discriminatory ability of our models, including or excluding BSI by quantifying the C-index. The BSI improved the C-index from 0.751 to 0.801 for OS. Median OS was not reached in patients with a BSI ≤1.9 and median OS was 34.8 months in patients with a BSI >1.9 (p = 0.039).
The pretreatment BSI and patients’ age are independent prognostic factors for patients with hormone-naive, bone metastatic prostate cancer.
Prostate cancer is the most common noncutaneous cancer, and the second most frequent cause of death from cancer among men in the USA. In Japan, 11,507 men were estimated to die of prostate cancer in 2014, making this disease the sixth leading cause of death from cancer . The incidence of prostate cancer is lower in Japan than in the USA and other Western countries. However, this incidence has been gradually increasing in Japan in recent years . Huggins and Hodges  reported the efficacy of androgen deprivation therapy for advanced prostate cancer in 1941. Although 80–90 % of prostate cancers with metastasis respond to initial androgen ablation therapy, most patients finally develop castration-resistant prostate cancer (CRPC) [3, 4]. Patients with CRPC show progression of systemic symptoms and local complications. One report showed that the median survival time among patients with advanced prostate cancer was 29 to 34 months from initial treatment , and another study reported a 5-year survival rate of 20–30 % . Because these reports showed a wide range of survival probability, more accurate information on patients’ characteristics related to survival is required.
Several groups have reported prognostic models for survival of patients with progressive disease [7–10]. Most of these reports were of prognostic models for patients with CRPC. There are few reports on prognostic models for patients with metastatic prostate cancer before treatment. A large study on the prognosis of patients with pre-hormonal therapy prostate cancer was reported in Japan and the USA . However, in this previous study, the endpoint was not survival, but recurrence. Previously, we reported a nomogram for overall survival (OS) of patients with bone-metastatic, hormone-naive prostate cancer . This nomogram comprised five pretreatment prognostic factors (patient’s age, clinical T stage, classification of bone metastasis extension [extent of disease on bone scan, EOD scores] , Gleason scores, and PSA) selected by multivariate analysis. Among these prognostic factors, only EOD scores are subjective and semi-quantitative. Therefore, an objective and quantitative scoring system for evaluation of bone metastasis might be ideal and warranted. Recently, a computer-aided diagnosis system (BONENAVI) for bone scans has been developed. The BONENAVI system can calculate the “bone scan index (BSI)”, which provides an objective and quantitative measure of the percentage of the skeleton involved by bone metastases . It is anticipated that a BSI that uses a computer-aided diagnosis system for bone scans will become an objective and quantitative clinical tool for evaluating bone metastatic prostate cancer. The BSI has been reported as being useful as a survival predictor among men with prostate cancer with various conditions such as hormone-naive prostate cancer or CRPC [13, 14].
In this study, we analyzed the relationship between the prognosis of prostate cancer and pretreatment clinical factors, including the BSI as calculated by BONENAVI for OS of patients with bone metastasis. This study aimed to determine whether the BSI is useful as a prognostic marker of hormone-naive, pretreatment prostate cancer with bone metastasis.
Patients and treatments
From 2010–2014, 60 consecutive patients with bone-metastatic prostate cancer were treated at Yokohama City University Hospital and associated hospitals. All of the patients already had metastasis at the time of diagnosis, and none of the patients had been previously treated. This study was a retrospective cohort study. This study was approved by institutional review board of Yokohama City University Medical Center and all subjects signed a written informed consent form.
All of the patients had adenocarcinoma of the prostate, confirmed histologically, with distant metastasis (any T, any N, M1b-c). The 2009 TNM clinical staging system and 2005 International Society of Urologic Pathology Gleason grading system were used. In all of the patients, the clinical stage was evaluated by chest and body computed tomography and bone scans. Based on the number or extent of metastases, the scans were divided into the following five grades according to the EOD on the bone scan : 0, normal or abnormal due to benign bone disease; 1, fewer than six bony metastases, each of which is less than 50 % of the size of a vertebral body (one lesion approximately the size of a vertebral body was counted as two lesions); 2, from 6–20 bone metastases, sized as described above; 3, more than 20 metastases, but fewer than those observed in a “superscan”; and 4, a “superscan” or its equivalent (i.e., more than 75 % of the ribs, vertebrae, and pelvic bones).
Each hospital used the same treatment protocol. All of the patients were initially treated with androgen deprivation therapy (medical or surgical castration with or without anti-androgen). None of the patients received prostatectomies or radiation therapy as initial treatments. After failed initial androgen ablation therapy, almost all of the patients subsequently underwent substitution treatment comprising anti-androgen therapy, anti-androgen withdrawal therapy, and/or oral low-dose steroid therapy. Some patients received bisphosphonate and cytotoxic therapy, such as docetaxel or estramustine, after development of CRPC. None of the patients received abiraterone, enzalutamide, sipuleucel-T, or cabazitaxel because abiraterone, enzalutamide, and cabazitaxel were not approved until 2013, and sipuleucel-T has not been approved yet in Japan. In the terminal state, palliative therapy and pain control with morphine, palliative external beam radiation, and strontium were used as appropriate. Progression was defined according to recommendations of the Cancer Clinical Trials Working Group 2 (PCWG2) .
Bone Scan Index (BSI)
The automated method for analysis of anterior and posterior whole-body bone scan images has been described previously . Each individual hot spot was classified as metastasis or no metastasis by an artificial neural network. The BSI was calculated as the percentage of the sum of all hot spots classified as bone metastases by the artificial neural network values. For calculation of the BSI, we used BONENAVI version 2 (Fujifilm RI Pharma, Co. Ltd., Tokyo, Japan; Exini Bone, Exini Diagnostics, Lund, Sweden) .
Spearman’s rank correlation coefficient was used for analysis of the correlation between EOD scores and the BSI. If there was a strong correlation, these two factors were not analyzed simultaneously because of multicollinearity.
We investigated the usefulness of a BSI as a predictor of OS. Univariate and multivariate analysis of OS was performed using a Cox proportional hazards regression model with stepwise regression analysis. Relative risks and 95 % confidence intervals (CIs) were derived. In univariate or multivariate analysis, all variables (age, PSA, Gleason scores, clinical TNM stage, and the BSI) were analyzed as binary variables. The cutoff point for all variables were determined using the median value of each variable. The C-index was used for discriminatory ability of our models.
The Kaplan–Meier product-limit estimator was used to estimate the distribution of survival. The log-rank test was used for analysis of the survival differences between patients with a BSI over the cut-off value and patients with a BSI under the cut-off value.
All tests were two-sided, and the significance level was fixed at alpha = 0.05. All analyses were conducted with IBM SPSS Statistics for Windows, ver. 22 (IBM Corp., Armonk, NY) and the R stats package (R Foundation for Statistical Computing, Vienna, Austria).
The characteristics of the patients are listed in Table 1. The median age of the patients was 72 years and the median follow-up duration was 21.4 months. The median initial PSA value was 247.0 ng/ml. The median BSI was 1.9 % (range: 0.0–13.2 %). Zoledronic acid was used for 28 (46.7 %) patients. Clinical or PSA progression occurred in 37 (61.7 %) patients and 18 (30.0 %) received docetaxel.
The correlation of EOD scores and the BSI is shown in Fig. 1. The median BSI stratified by EOD 1, 2, 3, and 4 was 0.38, 1.60, 7.38, and 9.29, respectively. There was a significant correlation between EOD and the BSI (rs = 0.769). Because of the strong correlation, we did not analyze the BSI and EOD in multivariate analysis simultaneously owing to multicollinearity.
Death occurred in 16 (26.7 %) patients. Of these patients, 15 (25.0 %) deaths were due to prostate cancer. The median OS was not reached.
In univariate analysis, age (≤72 years old vs >72 years old, hazard ratio (HR): 4.136, 95 % confidential interval (CI): 1.332–12.843, p = 0.014) and BSI (>1.9 % vs ≤1.9 %, HR: 2.921, 95 % CI: 1.006–8.483, p = 0.049) were prognostic factors for OS (Table 2).
In multivariate analysis, age (≤72 years old vs >72 years old, HR: 5.811, 95 % CI: 1.656–20.386, p = 0.006) and BSI (>1.9 % vs ≤1.9 %, HR: 4.676, 95 % CI: 1.238–17.661, p = 0.023) were independent prognostic factors for OS. (Table 2).
We evaluated the discriminatory ability of our models by quantifying the C-index. The C-index of our model including BSI was 0.801 for OS. We also analyzed the discriminatory ability of the model excluding BSI and the model including EOD instead of BSI by quantifying the C-index. The C-index of the models excluding BSI was 0.751 for OS. The C-index of the models including EOD instead of BSI was 0.811.
We compared the survival probability by BSI categories. Median OS was not reached in patients with a BSI ≤1.9 and median OS was 34.8 months in patients with a BSI >1.9 (p = 0.039) (Fig. 2).
Several groups have reported prognostic models for survival of patients with progressive prostate cancer [7–10]. Almost all reports were of prognostic models of patients with CRPC. However, only a few reports have described prognostic models for hormone-naive progressive prostate cancer before treatment [10, 18–20].
Hussain et al.  reported various risk predictors of OS (SWOG9346). They identified baseline variables, such as bone pain, performance status, Gleason sum, weight change, positive lymph node metastasis, pre-study PSA increments, and PSA levels after treatment as strong prognostic factors for OS. Coopeberg et al.  also reported a large study on prostate cancer prognosis in hormone-naive patients in Japan and the USA. They assessed 13,740 US men and 19,265 Japanese men with prostate cancer and developed the Japan Cancer of the Prostate Risk Assessment (J-CAPRA). The CAPRA score, which ranges from 0–12, and is based on the Gleason sum, serum PSA levels at initial treatment, and clinical stage, can predict progression-free survival after primary androgen deprivation therapy. Although the endpoint of the J-CAPRA is progression-free survival, we are interested in developing an OS prognostic model for patients with hormone-naive metastatic prostate cancer. Progression-free survival is predictive of OS in men with CRPC , although the association between progression-free survival and OS is relatively weak. Some reports have indicated improvement in OS without an increase in progression-free survival  or improvement in progression-free survival without an increase in survival . Accurate prediction models for prostate cancer survival probability would be helpful for counseling of patients.
In previous studies, EOD was reported to be a strong predictor of survival [7, 8]. EOD classification is used for evaluation of bone metastases. Because EOD is a subjective and semi-quantitative factor, more accurate and objective parameters for evaluation of spread of bone metastasis are warranted. Recently, the BSI, using a computer-aided diagnosis system for bone scans, has been expected to be an objective and quantitative clinical tool for evaluation of bone metastatic prostate cancer. The BSI is useful for prediction of survival in men with prostate cancer in various situations, such as hormone-naive cancer or CRPC [11, 13, 14, 25–29]. In this study, we measured BSI using BONENAVI version 2, which is a computer-assisted diagnosis system for bone scanning. BONENAVI version 2 software is based on a multi-center database and it showed excellent sensitivity and specificity for diagnosis of skeletal metastasis [17, 30].
In the hormone-naive setting, usefulness of the BSI as prognostic marker was reported by Koboteh and Poulsen [13, 27]. Koboteh and colleagues  reported that the BSI and Gleason scores were independent predictive factors of survival in men with a high risk of prostate cancer who were treated by androgen ablation therapy, although clinical stage and pre-treatment PSA levels were not prognostic factors. Poulsen et al.  also reported that only the BSI was an independent prognostic factor for survival of men with metastatic prostate cancer, although PSA and Gleason scores were not.
In our study, we analyzed the correlations between survival and clinical factors, including age, pre-treatment PSA levels, Gleason scores, clinical TNM classifications, and the BSI using the Cox proportional hazard regression model. EOD scores were significantly associated with the BSI value. EOD scores were not included in multivariate analysis because of multicollinearity. In results, patients’ age and the BSI were independent prognostic factor for OS of men with metastatic, hormone-naive prostate cancer. Because the BSI is continuous variables, the BSI is more ideal and accurate for prediction of individual patient’s prognosis compared with EOD scores. Moreover, we evaluated the discriminatory ability of our models by quantifying the C-index. The C-index of our model including BSI was 0.801 for OS. We also analyzed the discriminatory ability of the model excluding BSI and the model including EOD instead of BSI by quantifying the C-index. The C-index of the model excluding BSI was 0.751 for OS. The C-index of the models including EOD instead of BSI was 0.811.
A limitation of this study is the fact that patients who were enrolled in the study had various health statuses and complications. Our models did not consider health status or patients’ complications that may affect prostate cancer treatment outcomes [21, 31, 32]. Patients with prostate cancer are much older than those with other malignancies. Health status and complications should be classified in the rating score and included as predictive factors in the prognostic model. Our study did not include blood data. In previous studies, blood data, such as serum hemoglobin, alkaline phosphatase, lactate dehydrogenase, and C-reactive protein were reported as survival predictive factors for patients with prostate cancer [19, 20, 33]. Bone pain at diagnosis is also a strong predictor of survival . Unfortunately, data regarding pain at baseline were not available in this study. Lack of data about PSA kinetics was also limitation of our study. Finally, our population was small, observation periods were relatively short, and number of event was also small. Thus, some of the estimates have wide confidential intervals. More large populations and long observations are required for establishing the usefulness of the BSI as a prognostic factor. Although some study limitations exist, the BSI before treatment might be useful for a prognostic biomarker of hormone-naive, bone metastatic prostate cancer.
The pretreatment BSI and patients’ age are independent prognostic factors for patients with hormone-naive, bone metastatic prostate cancer. Compared with EOD, the BSI is a quantitative and objective factor for assessing bone metastases, and a more ideal factor for predicting prognosis in men with metastatic, hormone-naive prostate cancer.
bone scan index
castration-resistant prostate cancer
extent of disease
Japan Cancer of the Prostate Risk Assessment
Matsuda A, Matsuda T, Shibata A, Katanoda K, Sobue T, Nishimoto H, et al. Cancer incidence and incidence rates in Japan in 2008: a study of 25 population-based cancer registries for the Monitoring of Cancer Incidence in Japan (MCIJ) project. Jpn J Clin Oncol. 2014;44(4):388–96.
Huggins C, Hodges CV. Studies on prostatic cancer I: the effect of castration, of estrogen, and of androgen injection on serum phosphatases in metastatic carcinoma of the prostate. Cancer Res. 1941;1:293–7.
Small EJ, Vogelzang NJ. Second-line hormonal therapy for advanced prostate cancer: a shifting paradigm. J Clin Oncol. 1997;15(1):382–8.
Oh WK, Kantoff PW. Management of hormone refractory prostate cancer: current standards and future prospects. J Urol. 1998;160(4):1220–9.
Eisenberger MA, Blumenstein BA, Crawford ED, Miller G, McLeod DG, Loehrer PJ, et al. Bilateral orchiectomy with or without flutamide for metastatic prostate cancer. N Engl J Med. 1998;339(15):1036–42.
Maximum androgen blockade in advanced prostate cancer: an overview of the randomised trials. Prostate Cancer Trialists' Collaborative Group. Lancet 2000, 355(9214):1491–1498.
Miyoshi Y, Noguchi K, Yanagisawa M, Taguri M, Morita S, Ikeda I, et al. Nomogram for overall survival of Japanese patients with bone-metastatic prostate cancer. BMC Cancer. 2015;15(1):338.
Soloway MS, Hardeman SW, Hickey D, Raymond J, Todd B, Soloway S, et al. Stratification of patients with metastatic prostate cancer based on extent of disease on initial bone scan. Cancer. 1988;61(1):195–202.
Cooperberg MR, Hinotsu S, Namiki M, Ito K, Broering J, Carroll PR, et al. Risk assessment among prostate cancer patients receiving primary androgen deprivation therapy. J Clin Oncol. 2009;27(26):4306–13.
Armstrong AJ, Garrett-Mayer ES, Yang YC, de Wit R, Tannock IF, Eisenberger M. A contemporary prognostic nomogram for men with hormone-refractory metastatic prostate cancer: a TAX327 study analysis. Clin Cancer Res. 2007;13(21):6396–403.
Armstrong AJ, Kaboteh R, Carducci MA, Damber JE, Stadler WM, Hansen M, et al. Assessment of the bone scan index in a randomized placebo-controlled trial of tasquinimod in men with metastatic castration-resistant prostate cancer (mCRPC). Urol Oncol. 2014;32(8):1308–16.
Wakabayashi H, Nakajima K, Mizokami A, Namiki M, Inaki A, Taki J, et al. Bone scintigraphy as a new imaging biomarker: the relationship between bone scan index and bone metabolic markers in prostate cancer patients with bone metastases. Ann Nucl Med. 2013;27(9):802–7.
Poulsen MH, Rasmussen J, Edenbrandt L, Hoilund-Carlsen PF, Gerke O, Johansen A, et al. Bone Scan Index predicts outcome in patients with metastatic hormone sensitive prostate cancer. BJU Int. 2015.
Mitsui Y, Shiina H, Yamamoto Y, Haramoto M, Arichi N, Yasumoto H, et al. Prediction of survival benefit using an automated bone scan index in patients with castration-resistant prostate cancer. BJU Int. 2012;110(11 Pt B):E628–634.
Scher HI, Halabi S, Tannock I, Morris M, Sternberg CN, Carducci MA, et al. Design and end points of clinical trials for patients with progressive prostate cancer and castrate levels of testosterone: recommendations of the Prostate Cancer Clinical Trials Working Group. J Clin Oncol. 2008;26(7):1148–59.
Ulmert D, Kaboteh R, Fox JJ, Savage C, Evans MJ, Lilja H, et al. A novel automated platform for quantifying the extent of skeletal tumour involvement in prostate cancer patients using the Bone Scan Index. Eur Urol. 2012;62(1):78–84.
Koizumi M, Wagatsuma K, Miyaji N, Murata T, Miwa K, Takiguchi T, et al. Evaluation of a computer-assisted diagnosis system, BONENAVI version 2, for bone scintigraphy in cancer patients in a routine clinical setting. Ann Nucl Med. 2015;29(2):138–48.
Armstrong AJ, Garrett-Mayer E, de Wit R, Tannock I, Eisenberger M. Prediction of survival following first-line chemotherapy in men with castration-resistant metastatic prostate cancer. Clin Cancer Res. 2010;16(1):203–11.
Smaletz O, Scher HI, Small EJ, Verbel DA, McMillan A, Regan K, et al. Nomogram for overall survival of patients with progressive metastatic prostate cancer after castration. J Clin Oncol. 2002;20(19):3972–82.
Halabi S, Small EJ, Kantoff PW, Kattan MW, Kaplan EB, Dawson NA, et al. Prognostic model for predicting survival in men with hormone-refractory metastatic prostate cancer. J Clin Oncol. 2003;21(7):1232–7.
Iasonos A, Schrag D, Raj GV, Panageas KS. How to build and interpret a nomogram for cancer prognosis. J Clin Oncol. 2008;26(8):1364–70.
Halabi S, Vogelzang NJ, Ou SS, Owzar K, Archer L, Small EJ. Progression-free survival as a predictor of overall survival in men with castrate-resistant prostate cancer. J Clin Oncol. 2009;27(17):2766–71.
Kantoff PW, Higano CS, Shore ND, Berger ER, Small EJ, Penson DF, et al. Sipuleucel-T immunotherapy for castration-resistant prostate cancer. N Engl J Med. 2010;363(5):411–22.
Sternberg CN, Petrylak DP, Sartor O, Witjes JA, Demkow T, Ferrero JM, et al. Multinational, double-blind, phase III study of prednisone and either satraplatin or placebo in patients with castrate-refractory prostate cancer progressing after prior chemotherapy: the SPARC trial. J Clin Oncol. 2009;27(32):5431–8.
Dennis ER, Jia X, Mezheritskiy IS, Stephenson RD, Schoder H, Fox JJ, et al. Bone scan index: a quantitative treatment response biomarker for castration-resistant metastatic prostate cancer. J Clin Oncol. 2012;30(5):519–24.
Horikoshi H, Kikuchi A, Onoguchi M, Sjostrand K, Edenbrandt L. Computer-aided diagnosis system for bone scintigrams from Japanese patients: importance of training database. Ann Nucl Med. 2012;26(8):622–6.
Kaboteh R, Damber JE, Gjertsson P, Hoglund P, Lomsky M, Ohlsson M, et al. Bone Scan Index: a prognostic imaging biomarker for high-risk prostate cancer patients receiving primary hormonal therapy. EJNMMI Res. 2013;3(1):9.
Kaboteh R, Gjertsson P, Leek H, Lomsky M, Ohlsson M, Sjostrand K, et al. Progression of bone metastases in patients with prostate cancer - automated detection of new lesions and calculation of bone scan index. EJNMMI Res. 2013;3(1):64.
Kalderstam J, Sadik M, Edenbrandt L, Ohlsson M. Analysis of regional bone scan index measurements for the survival of patients with prostate cancer. BMC Med Imaging. 2014;14:24.
Nakajima K, Nakajima Y, Horikoshi H, Ueno M, Wakabayashi H, Shiga T, et al. Enhanced diagnostic accuracy for quantitative bone scan using an artificial neural network system: a Japanese multi-center database project. EJNMMI Res. 2013;3(1):83.
Tewari A, Raman JD, Chang P, Rao S, Divine G, Menon M. Long-term survival probability in men with clinically localized prostate cancer treated either conservatively or with definitive treatment (radiotherapy or radical prostatectomy). Urology. 2006;68(6):1268–74.
Droz JP, Balducci L, Bolla M, Emberton M, Fitzpatrick JM, Joniau S, et al. Management of prostate cancer in older men: recommendations of a working group of the International Society of Geriatric Oncology. BJU Int. 2010;106(4):462–9.
Mikkola A, Aro J, Rannikko S, Ruutu M, Finnprostate G. Prognostic grouping of metastatic prostate cancer using conventional pretreatment prognostic factors. Scand J Urol Nephrol. 2009;43(4):265–70.
This work was supported by a MEXT/JSPS KAKENHI grant.
The authors declare that they have no competing interests.
KN and YM designed the study. YM contributed to the statistical analysis. YM and SY wrote the manuscript. TK, YH, JT, KK, MM, ST, YY, MY, and HU supported the statistical analysis and contributed to collecting the clinical data. All authors read and approved the final manuscript.
About this article
Cite this article
Miyoshi, Y., Yoneyama, S., Kawahara, T. et al. Prognostic value of the bone scan index using a computer-aided diagnosis system for bone scans in hormone-naive prostate cancer patients with bone metastases. BMC Cancer 16, 128 (2016). https://doi.org/10.1186/s12885-016-2176-6