Open Access Highly Accessed Original research

Bone Scan Index: a prognostic imaging biomarker for high-risk prostate cancer patients receiving primary hormonal therapy

Reza Kaboteh1*, Jan-Erik Damber2, Peter Gjertsson1, Peter Höglund3, Milan Lomsky1, Mattias Ohlsson4 and Lars Edenbrandt15

Author Affiliations

1 Department of Molecular and Clinical Medicine, Clinical Physiology, Sahlgrenska University Hospital, Sahlgrenska Academy at the University of Gothenburg, Gothenburg, 413 45, Sweden

2 Department of Urology, Institute of Clinical Sciences, Gothenburg University, Gothenburg, 413 45, Sweden

3 Competence Centre for Clinical Research, Lund University Hospital, Lund, 221 00, Sweden

4 Department of Theoretical Physics, Lund University, Lund, 221 00, Sweden

5 Department of Clinical Sciences, Malmö, Lund University, Lund, 221 00, Sweden

For all author emails, please log on.

EJNMMI Research 2013, 3:9  doi:10.1186/2191-219X-3-9


The electronic version of this article is the complete one and can be found online at: http://www.ejnmmires.com/content/3/1/9


Received:16 November 2012
Accepted:24 January 2013
Published:6 February 2013

© 2013 Kaboteh et al.; licensee Springer.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

The objective of this study was to explore the prognostic value of the Bone Scan Index (BSI) obtained at the time of diagnosis in a group of high-risk prostate cancer patients receiving primary hormonal therapy.

Methods

This was a retrospective study based on 130 consecutive prostate cancer patients at high risk, based on clinical stage (T2c/T3/T4), Gleason score (8 to 10) and prostate-specific antigen (PSA) (> 20 ng/mL), who had undergone whole-body bone scans < 3 months after diagnosis and who received primary hormonal therapy. BSI was calculated using an automated method. Cox proportional-hazards regression models were used to investigate the association between clinical stage, Gleason score, PSA, BSI and survival. Discrimination between prognostic models was assessed using the concordance index (C-index).

Results

In a multivariate analysis, Gleason score (p = 0.01) and BSI (p < 0.001) were associated with survival, but clinical stage (p = 0.29) and PSA (p = 0.57) were not prognostic. The C-index increased from 0.66 to 0.71 when adding BSI to a model including clinical stage, Gleason score and PSA. The 5-year probability of survival was 55% for patients without metastases, 42% for patients with BSI < 1, 31% for patients with BSI = 1 to 5, and 0% for patients with BSI > 5.

Conclusions

BSI can be used as a complement to PSA to risk-stratify high-risk prostate cancer patients at the time of diagnosis. This imaging biomarker, reflecting the extent of metastatic disease, can be of value both in clinical trials and in patient management when deciding on treatment.

Keywords:
Image analysis; Radionuclide imaging; Bone metastases; Prostate cancer

Background

Life expectancy is a major factor to be considered in the management of prostate cancer patients. Risk stratification schemes based on clinical T stage, Gleason score and prostate-specific antigen (PSA) are widely used to estimate risk in individual patients. The extent of bone metastases is also associated with survival [1,2], but there has not been any clinically useful technique of quantifying the skeletal tumour burden and including this information in the risk assessment. Bone scintigraphy, however, is commonly used to assess skeletal tumour burden in prostate cancer patients, both in clinical routine and in nearly every clinical trial. In order to extract as much clinical information as possible from the bone scans, the Bone Scan Index (BSI) was developed as a quantitative tool to improve the interpretability and clinical relevance of the bone scan [3]. BSI is a method of expressing the tumour burden in the bone as a percentage of the total skeletal mass.

BSI has recently been presented as a response indicator in patients with castration-resistant metastatic prostate cancer who have received chemotherapy [4]. The patients had bone scans at baseline and at 3-month and 6-month follow-up, and BSI changes post-treatment were a significant prognostic factor for survival. A doubling of BSI post-treatment resulted in a 1.9-fold increase in the risk of death. These results showed the feasibility of capturing bone scintigraphy data as a single quantitative measure and thereby allowing bone scans to be explored as imaging biomarkers. Furthermore, changes in BSI post-treatment were significantly associated with survival, but post-treatment changes in PSA were not, while adjusting for changes in BSI.

In order to improve the usefulness of BSI, we recently presented a fully automated method of quantifying BSI [5]. The automated method was evaluated in a group of newly diagnosed prostate cancer patients, and BSI increased the predictive accuracy for death as a result of prostate cancer when added to a model containing pre-treatment clinical stage, Gleason score and PSA. The 5-year probability of dying as a result of prostate cancer was under 6% for patients with BSI < 0.1 and over 78% for patients with BSI > 1.0 at the time of diagnosis. Initial treatment may influence outcomes, and that was not adjusted for in our previous study. Furthermore, low- and intermediate-risk patients (median PSA = 16.8 ng/mL) were included, and many of these patients probably had a specific indication for receiving a bone scan, for example, bone pain. This may have resulted in a selection bias. High-risk patients at the time of diagnosis have been referred for bone scans before treatment decisions as part of the clinical routine. The purpose of this study was to explore the prognostic value of BSI obtained at the time of diagnosis in a group of high-risk prostate cancer patients receiving primary hormonal therapy.

Methods

Patients

All patients with the diagnosis of prostate cancer who, during the period 2002 to 2008, had undergone a whole-body bone scan < 3 months from the time of diagnosis at the Sahlgrenska University Hospital, Gothenburg, Sweden, were retrospectively considered for inclusion in the study. Only patients at high risk who received primary hormonal therapy were included (chemical or surgical castration (n = 97), anti-androgen deprivation as monotherapy (n = 13), total androgen blockade (n = 11), neoadjuvant treatment before curative treatment (n = 6), a combination of external irradiation with anti-androgen deprivation (n = 2), radical prostatectomy (n = 1)). None of the patients had received any type of PCa treatment prior to the bone scan. We defined patients as high-risk if at least one of the following criteria was met [6]:

• T2c/T3/T4

• Gleason score 8 to 10

• PSA level > 20 ng/mL

In total, 130 patients were included in our analysis. Patient characteristics are presented in Table 1. During the study period, patients at high risk at the time of diagnosis were referred for bone scan examination before treatment decisions as part of the clinical routine. The study was approved by the Research Ethics Committee at Gothenburg University.

Table 1. Patient characteristics (n = 130)

Data collection

Overall survival was defined as the time from bone scan examination to death from any cause. Survival data and PSA values collected from the computerised medical records were updated until 24 September 2010. A total of 76 were dead at follow-up with a median survival time of 2.7 years (IQR 1.6 to 4.2) and 54 were alive with a median follow-up time of 5.1 years (IQR 3.9 to 8.5). Data on clinical T stage and Gleason score at the time of diagnosis were collected from the Swedish National Cancer Registry [7].

The percentage of the skeleton affected by tumour mass in a bone scan was measured by calculating the BSI. We have recently presented an automated method [5] based on the clinically validated methodology for manually computing BSI as presented by a group at the Memorial Sloan-Kettering Cancer Center in New York [3,8,9]. The automated method is trained to mimic an expert reader in distinguishing hotspots due to metastases from those caused by factors such as degenerative disease or fractures. Manual correction was required in less than 5% of the patients to exclude hotspots clearly misclassified and representing features such as a very large urinary bladder, a urinary catheter attached to a drainage bag or urine contamination. No other manual intervention was applied. The method is implemented in the commercially available software package EXINI bone™ (EXINI Diagnostics AB, Lund, Sweden).

Statistical methods

Cox proportional-hazards regression models were used to investigate the association between clinical stage, Gleason score, PSA, BSI and survival, both in univariate and multivariable modelling. Hazard ratios were estimated from the Cox models together with 95% confidence intervals. The concordance index (C-index) was used to discriminate between the different survival models [10]. The significance of a difference in C-index between different models was calculated using the method described by Haibe-Kains et al. [11].

Kaplan-Meier estimates of the survival function were used together with the log-rank test to indicate a significant difference between four groups of patients: without metastases (M0), with BSI < 1, BSI = 1 to 5 and BSI > 5. The BSI levels 1 and 5 were used to categorize the patients since that resulted in approximately the same number of patients in each of the metastatic groups.

The survival rate expected in a group of men in the general population with the same age as the M0 group was calculated. This control survival curve of the general Swedish population was generated using the Hakulinen method [12], using data collected from the administrative agency Statistics Sweden. This group is similar to the M0 group with respect to other possible factors affecting survival, except for the patients of the M0 group having prostate cancer.

In the survival analysis, data were censored at a follow-up time of 5 years. All analyses were carried out using the R statistical computing environment.

Results

The 5-year probability of survival was 41% in the total group of high-risk patients. Metastatic disease was present in 58/130 (45%) of the patients. The 5-year survival probabilities were 24% and 55% in patients with (M1) and without (M0) metastases, respectively. Considering the tumour burden, measured using the continuous variable BSI, and not only the presence of metastases (M0 or M1), it was possible to divide the metastatic group. We divided this group into three subgroups in order to have approximately 20 patients in each subgroup: BSI < 1 (n = 22), BSI = 1 to 5 (n = 16) and BSI > 5 (n = 20). Figure 1 shows one patient from each of these subgroups illustrating the extent of disease. The Kaplan-Meier curves for patients M0, BSI < 1, BSI = 1 to 5 and BSI > 5 were significantly different (p < 0.001), and the 5-year survival probabilities were 55%, 42%, 31% and 0%, respectively. The survival curve for the M0 patients is above the corresponding survival curve for the general male population of the same age for the first 4 years (Figure 2). The corresponding value for an age-matched control cohort was 65%.

thumbnailFigure 1. Three patients, one each from the groups BSI < 1, BSI = 1 to 5 and BSI > 5. The actual BSI values were 0.7, 2.2 and 8.1, respectively.

thumbnailFigure 2. Kaplan-Meier curves showing patient-survival probability stratified by BSI categories. The difference between the four groups of patients without metastases (M0) and with metastases BSI < 1, BSI = 1 to 5 and BSI > 5 was statistically significant (p < 0.001). The broken line shows an age-matched control survival curve for the M0 group.

In the univariate analysis clinical stage, Gleason score, PSA and BSI were all associated with survival (Table 2). In a multivariate analysis, Gleason score (p = 0.01) and BSI (p < 0.001) were associated with survival, but clinical stage (p = 0.29) and PSA (p = 0.57) were not prognostic. The C-index increased from 0.66 to 0.71 (p = 0.006) when adding BSI to a model including clinical stage, Gleason score and PSA.

Table 2. Survival analysis demonstrating association between clinical stage, Gleason score, PSA, BSI and survival

Discussion

BSI includes prognostic information in addition to that obtained from clinical T stage, Gleason score and PSA. Patients with metastatic disease and BSI < 1 showed a 5-year probability of survival of 42% compared to 31% for those with BSI = 1 to 5 and 0% for those with BSI > 5. These findings show that a quantitative measurement of tumour burden in the skeleton can be used to risk-stratify patients more efficiently than only using M-staging based on evidence of the presence or absence of metastatic spread.

The added value of BSI was shown as an increase in C-index from 0.66 to 0.71 when adding BSI to a model including clinical T stage, Gleason score and PSA. This is in agreement with the results presented by Ulmert et al. [5], where the corresponding increase when adding BSI was of the same magnitude, but from a higher level (0.77 to 0.82). The higher values in the Ulmert study are likely to be due to the fact that death as a result of prostate cancer was used in that study, while all-cause mortality or overall survival was used in this one. The cause of death for the patients was not available in this retrospective study. Furthermore, the common primary end point in large clinical trials exploring the effects of drugs for prostate cancer patients is ‘overall survival’. This study showed the value of BSI in a homogenous group of prostate cancer patients, with those at high risk receiving primary hormonal therapy. In the Ulmert study, low-risk patients were also included and initial treatment was not adjusted for.

The sensitivity for metastatic disease of bone scans is regarded to be high, that is, the number of false-negative cases is low, but some of our M0 patients are probably falsely classified as M0. These patients might at least partly contribute to the higher mortality in the M0 population after 4 years in comparison with the age-matched control group. A follow-up bone scan to validate the BSI from the time of diagnosis could therefore be of value. Change in BSI over time can also contain important information in the M1 group. Dennis et al. recently showed that on-treatment change in BSI was associated with overall survival for patients receiving chemotherapy [4]. Their study also showed that changes in PSA were not associated with survival while adjusting for changes in BSI.

Recently, another study demonstrated the value of BSI in patients with castration-resistant prostate cancer [13]. In a group of 42 patients who underwent taxane-based chemotherapy, reduction of BSI after treatment was associated with longer overall survival. The findings of recent BSI studies show the value of this imaging biomarker in both at the time of diagnosis and at a later stage of the disease.

Based on the results of the present study, BSI can be considered as an informative predictor of patient survival in men with high-risk disease treated with hormonal therapy. BSI can also be of value in clinical trials as an inclusion criterion so that baseline prognosis in treatment and placebo groups can be matched.

The purpose of this study was to generate evidence regarding how BSI can be used to risk-stratify prostate cancer patients by applying an automated software to bone scans - the established method of evaluating skeletal metastases. Bone scans are widely used in both clinical routine and clinical trials, but the methodology has its shortcomings. The radiotracer locates not tumour cells but regions of the skeleton that are actively undergoing tissue repair - a well-known sign of tumour involvement. Future developments in diagnostic imaging may provide improved and more specific methods of analysing skeletal lesions, including three-dimensional methods such as SPECT/CT and PET/CT. At present, bone scan is the dominant clinical method for evaluation of skeletal metastases. Furthermore, prognostic studies based on large patient groups and long follow-up time are not currently available for the new modalities. This type of bone scan study, showing that imaging can be used as a biomarker, will be important in terms of gaining experience, which can also be of value when new imaging modalities become widely used.

This study was based on a fully automated software for calculation of tumour burden in patients with metastatic disease. With this technique, widespread use in clinical routine and clinical trials is possible. Operator-dependent subjectivity in interpretation of bone scans, which has proved to be substantial [14], is eliminated, and a reproducible analysis can be obtained at different hospitals. The software has proved to be capable of differentiating between malignant lesions and hotspots due to degenerative disease, e.g. by considering whether uptake in the shoulder regions is symmetrical. In this study, only obvious false lesions were manually corrected, e.g. a very large urinary bladder, a urinary catheter attached to a drainage bag or urine contamination. The quality of the software was assessed in a study by Sadik et al. [15]. This study showed that 35 physicians improved their sensitivity to metastatic disease from 78% to 88% using the software for automated analysis of bone scans.

Conclusions

In conclusion, this study showed that BSI is strongly associated with overall survival in patients with high-risk prostate cancer receiving primary hormonal therapy and that BSI includes prognostic information in addition to clinical T stage, Gleason score and PSA. This study builds on previous studies and generates evidence that contributes to qualification of BSI as an imaging biomarker in prostate cancer patients. This biomarker has recently become automated and thus highly reproducible, eliminating operator-dependent subjectivity and providing prognostic information with a processing time < 10 s.

Competing interests

Lars Edenbrandt and Mattias Ohlsson are shareholders in EXINI Diagnostics AB (Lund, Sweden), which provides software for Bone Scan Index calculations. Reza Kaboteh, Jan-Erik Damber, Peter Gjertsson, Peter Höglund and Milan Lomsky indicated no potential conflicts of interest.

Authors’ contributions

RK and LE participated in the design of the study and in the analysis and interpretation of data, and drafted the manuscript. JED, PG and ML participated in the design of the study and in the analysis and interpretation of data. PH and MO carried out the analysis of data. All authors read and approved the final manuscript.

Acknowledgements

This study was supported by grants from the Swedish Research Council (Medicine), and Sahlgrenska Academy ALF grants.

References

  1. Soloway MS, Hardeman SW, Hickey D, Todd B, Soloway S, Raymond J, Moinuddin M: Stratification of patients with metastatic prostate cancer based on extent of disease on initial bone scan.

    Cancer 1988, 61:195-202. PubMed Abstract | Publisher Full Text OpenURL

  2. Noguchi M, Kikuchi H, Ishibashi M, Noda S: Percentage of the positive area of bone metastasis is an independent predictor of disease death in advanced prostate cancer.

    Brit J Cancer 2003, 88:195-201. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  3. Erdi YE, Humm JL, Imbriaco M, Yeung H, Larson SM: Quantitative bone metastases analysis based on image segmentation.

    J Nucl Med 1997, 38:1401-1406. PubMed Abstract | Publisher Full Text OpenURL

  4. Dennis ER, Jia X, Mezheritskiy IS, Stephenson RD, Schoder H, Fox JJ, Heller G, Scher HI, Larson AM, Morris MJ: Bone scan index: a quantitative treatment response biomarker for castration-resistant metastatic prostate cancer.

    JCO 2012, 30:519-524. Publisher Full Text OpenURL

  5. Ulmert D, Kaboteh R, Fox JJ, Savage C, Evans MJ, Lilja H, Abrahamsson PA, Björk T, Gerdtsson A, Bjartell A, Gjertsson P, Höglund P, Lomsky M, Ohlsson M, Richter J, Sadik M, Morris MJ, Scher HI, Sjöstrand K, Yu A, Suurküla M, Edenbrandt L, Larson SM: A novel automated platform for quantifying the extent of skeletal tumor involvement in prostate cancer patients using the Bone Scan Index.

    Eur Urol 2012, 62:78-84. PubMed Abstract | Publisher Full Text OpenURL

  6. D'Amico AV, Desjardin A, Chung A, Chen MH, Schultz D, Whittington R, Malkowicz SB, Wein A, Tomaszewski JE, Renshaw AA, Loughlin K, Richie JP: Assessment of outcome prediction models for patients with localized prostate carcinoma managed with radical prostatectomy or external beam radiation therapy.

    Cancer 1998, 82:1887-1896. PubMed Abstract | Publisher Full Text OpenURL

  7. Stattin P, Holmberg E, Johansson JE, Holmberg L, Adolfsson J, Hugosson J, National Prostate Cancer Register (NPCR) of Sweden: Outcomes in localized prostate cancer: National Prostate Cancer Register of Sweden follow-up study.

    J Natl Cancer Inst 2010, 102:950-958. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  8. Imbriaco M, Larson SM, Yeung HW, Mawlawi OR, Erdi Y, Venkatraman ES, Scher HI: A new parameter for measuring metastatic bone involvement by prostate cancer: the Bone Scan Index.

    Clin Cancer Res 1998, 4:1765-1772. PubMed Abstract | Publisher Full Text OpenURL

  9. Sabbatini P, Larson SM, Kremer A, Zhang ZF, Sun M, Yeung H, Imbriaco M, Horak I, Conolly M, Ding C, Ouyang P, Kelly WK, Scher HI: Prognostic significance of extent of disease in bone in patients with androgen-independent prostate cancer.

    J Clin Oncol 1999, 17:948-957. PubMed Abstract | Publisher Full Text OpenURL

  10. Harrell FE, Lee KL, Mark DB: Multivariable prognostic models: issues in developing models, evaluating assumptions and adequacy, and measuring and reducing errors.

    Stat Med 1996, 15:361-387. PubMed Abstract | Publisher Full Text OpenURL

  11. Haibe-Kains B, Desmedt C, Sotiriou C, Bontempi G: A comparative study of survival models for breast cancer prognostication on microarray data: does a single gene beat them all?

    Bioinformatics 2008, 24:2200-2208. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  12. Hakulinen T: Cancer survival corrected for heterogeneity in patient withdrawal.

    Biometrics 1982, 38:933-942. PubMed Abstract | Publisher Full Text OpenURL

  13. Mitsui Y, Shiina H, Yamamoto Y, Haramoto M, Arichi N, Yasumoto H, Kitagaki H, Igawa M: Prediction of survival benefit using an automated bone scan index in patients with castration-resistant prostate cancer.

    BJU Int 2012. Publisher Full Text OpenURL

  14. Sadik M, Suurkula M, Hoglund P, Jarund A, Edenbrandt L: Quality of planar whole-body bone scan interpretations — a nationwide survey.

    Eur J Nucl Med Mol Imaging 2008, 35:1464-1472. PubMed Abstract | Publisher Full Text OpenURL

  15. Sadik M, Suurkula M, Hoglund P, Jarund A, Edenbrandt L: Improved classification of planar whole-body scan using a computer-assisted diagnosis system: a multicenter, multiple-reader, multiple-case study.

    J Nucl Med 2009, 50:368-375. PubMed Abstract | Publisher Full Text OpenURL