- Research article
- Open Access
Risk factors, prognostic factors, and nomograms for bone metastasis in patients with newly diagnosed infiltrating duct carcinoma of the breast: a population-based study
BMC Cancer volume 20, Article number: 1145 (2020)
Breast cancer is the most common malignancy in women, and it is also the leading cause of death in female patients; the most common pathological type of BC is infiltrating duct carcinoma (IDC). Some nomograms have been developed to predict bone metastasis (BM) in patients with breast cancer. However, there are no studies on diagnostic and prognostic nomograms for BM in newly diagnosed IDC patients.
IDC patients with newly diagnosed BM from 2010 to 2016 in the Surveillance, Epidemiology and End Results (SEER) database were reviewed. Multivariate logistic regression analysis was used to identify risk factors for BM in patients with IDC. Univariate and multivariate Cox proportional hazards regression analysis were used to explore the prognostic factors of BM in patients with IDC. We then constructed nomograms to predict the risk and prognosis of BM for patients with IDC. The results were validated using bootstrap resampling and retrospective research on 113 IDC patients with BM from 2015 to 2018 at the Affiliated Hospital of Chengde Medical University.
This study included 141,959 patients diagnosed with IDC in the SEER database, of whom 2383 cases were IDC patients with BM. The risk factors for BM in patients with IDC included sex, primary site, grade, T stage, N stage, liver metastasis, race, brain metastasis, breast cancer subtype, lung metastasis, insurance status, and marital status. The independent prognostic factors were brain metastases, race, grade, surgery, chemotherapy, age, liver metastases, breast cancer subtype, insurance status, and marital status. Through calibration, receiver operating characteristic curve and decision curve analyses, we found that the nomogram for predicting the prognosis of IDC patients with BM displayed great performance both internally and externally.
These nomograms are expected to be a precise and personalized tool for predicting the risk and prognosis for BM in patients with IDC. This will help clinicians develop more rational and effective treatment strategies.
Breast cancer (BC) is the most common malignancy and the leading cause of death among all female cancer patients [1, 2]. Globally, there were approximately 2.1 million newly diagnosed female BC cases in 2018 . Recently, with the advancement of early diagnosis and comprehensive treatment, the mortality rate of BC has gradually decreased, and distant metastasis has become the main cause of death for these patients [4, 5]. It has been reported that the incidence of metastases in BC patients ranges from 20 to 30% . More importantly, bone metastasis (BM) accounts for 50% of all distant metastases in these patient . At present, most BC patients with BM receive palliative treatment . Although some patients choose surgery, it is not suitable for patients with multiple metastases or a poor overall health . Some studies have shown that the median survival for patients with breast cancer and BM is only 24–36 months .
The TNM staging system is the most common tool used to predict the prognosis of cancer patients by assessing tumor size and location (T), distant metastasis (M), and regional lymph node metastasis (N) . However, the TNM staging system does not sufficiently cover cancer biology or predict the outcome for all subtypes of BC . In particular, the TNM staging system fails to quantify the risk for patients with distant metastatic malignancies. Therefore, an increasing number of cancer-related nomograms (statistical tools to estimate the probability of survival or a specific result through a simple graphical representation) have been developed for predicting the prognosis of cancer patients . Nomograms have a number of advantages in predicting the prognosis of some malignant tumors compared to the traditional American Joint Committee for Cancer (AJCC) TNM staging system, making them a good alternative.
It is well established that histological subtypes of breast cancer affect prognosis, and the most common pathological type of BC is infiltrating duct carcinoma (IDC) . At present, there are no studies that have focused on diagnostic and prognostic nomograms for BM in newly diagnosed IDC patients. Therefore, it is necessary to fully understand the epidemiological characteristics of IDC patients with BM to identify the risk and prognostic factors for BM. Well-developed clinical nomograms can be used to predict individual outcomes, which is beneficial to both patients and clinicians .
Thus, the aim of this study was to develop a predictive model by analyzing the data of the Surveillance, Epidemiology and End Results (SEER) database to determine the risk and prognosis for BM in patients with IDC.
We included patients with newly diagnosed IDC in the SEER database from 2010 to 2016 in our study. Exclusion criteria were as follows: (1) patients with two or more primary malignancies; (2) patients whose pathological type was not IDC; (3) patients missing important clinical pathological information, including laterality, primary tumor site, grade, TNM stage, or estrogen receptor (ER) or progesterone receptor (PR) status, or HER2 status. Finally, 141,959 patients diagnosed with IDC were included in the present study, of whom 2383 patients (1.68%) had BM, while 139,576 patients (98.32%) did not. In addition, we retrospectively collected data for IDC patients with BM from the Affiliated Hospital of Chengde Medical University (AHOCMU) between 2015 and 2018 as an external validation cohort for our research.
The variables were selected to identify the risk factors of BM in IDC patients are as follows: age at diagnosis, sex, race, tumor site, laterality, grade, T stage, N stage, liver metastasis, brain metastasis, lung metastasis, breast cancer subtype, ER status, PR status, HER2 status, insurance, and marital status. In our research, we also performed the survival analyses to study the prognostic factors of IDC patients with BM. In addition to the above variables, the treatment information, including surgery, radiotherapy, and chemotherapy, were also included to study the prognostic factors. Moreover, patients with overall survival (OS) less than 1 month were also excluded from the survival analyses. In the survival analysis, the main endpoint of our study was OS, which was defined as the date from diagnosis to death (due to any cause) or the date of the last follow-up. Risk of developing metastasis was defined as the risk of bone metastasis when the patient was first diagnosed with IDC of the breast. Survival prognosis was defined as the OS of the patient who was first diagnosed with IDC of the breast. Our study was approved by the Institutional Research Committee from AHOCMU.
Development of a diagnostic nomogram
All statistical analyses in our research were performed in R software (version 3.6.1). To identify the risk factors of BM in IDC patients, univariate analysis was performed. Comparisons of continuous data were performed by independent t-tests, while the chi-square test or the Fisher exact probability method were used for categorical data. Variables with a P value < 0.05 in the univariate analysis were included in the multivariate logistic analysis to identify the risk factors for BM in IDC patients. Based on independent risk factors, the rms package was used to build a nomogram and calculate the individual risk score. Meanwhile, the receiver operating characteristic (ROC) curve was plotted, and the area under the curve (AUC) was used to show the discrimination of the nomogram. Moreover, a calibration curve and decision curve analyses (DCA) were performed to evaluate the nomogram .
Development of a prognostic nomogram
To identify the prognostic factors of IDC patients with BM, 2383 patients were included to perform survival analyses. All BM patients were randomly divided into training (n = 1671) and validation (n = 712) cohorts with a ratio of 7:3. The classification process was completely randomized and it was performed in R software. The best age cutoff values for OS were determined by X-tile software; patients were divided into high, middle, and low groups. We then performed univariate Cox proportional hazards regression analysis to determine the OS-related variables. Afterward, significant variables in the univariate Cox proportional hazards regression analyses were incorporated into the multivariate Cox proportional hazards regression analyses to determine the independent prognostic factors for IDC patients with BM. Then, a nomogram based on the independent prognostic factors was established to predict the OS for IDC patients with BM. Additionally, time-dependent ROC curves of 1, 3, and 5 years were generated, and the corresponding time-dependent AUCs were used to show the discrimination of the nomogram. Calibration curves and DCA of 1, 3, and 5 years were established. To further validate that the nomogram could perform well in an independent cohort, we validated the nomogram with data from the SEER validation cohort and the AHOCMU cohort. Time-dependent ROC curve, calibration curve, and DCA were also performed in the validation cohort. In the present study, a P value< 0.05(two side) was identified as statistical significance.
Baseline characteristics of the study population
Based on our criteria, a total of 141,959 IDC patients from the SEER database were included, and an additional 113 IDC patients with BM were identified from the AHOCMU for this study. Additionally, 1671 patients were included in the training cohort and 712 patients were included in the validation cohort. As shown in Table 1, 99.24% of the patients were female and 80.23% were white. The most common tumor grade of differentiation was grade II (42.01%). The most common primary site location was in the upper-outer quadrant of the breast (39.04%). There was minimal laterality, with left primary site origins accounting for 50.54% of the study group and right primary site origins accounting for 49.46%. The most common T and N stages were T1 (63.60%) and N0 (69.79%). Regarding the classifications of breast cancer subtypes, luminal A (HR+/HER2-) accounted for 71.09%. A total of 1133 (0.80%) patients had lung metastases, 202 (0.14%) patients had brain metastasis and 979 (0.69%) patients had liver metastases. Most patients were insured (98.10%) and married (62.57%). In our study, most patients were positive for PR (71.76%) and ER (81.51%). Regarding therapy, 136,494 (96.15%) of the patients underwent surgery, 61,831 (43.56%) underwent chemotherapy, and 80,424 (56.65%) underwent radiotherapy. Table 2 displays information on the clinical and pathological features for the IDC patients with BM.
Risk factors for IDC patients with BM
As shown in Table 3, variables with a P value < 0.05 in the univariate analysis were included in the multivariate logistic regression analysis to determine the risk factors for BM in IDC patients. The results revealed that sex, primary site, grade, T stage, N stage, brain metastasis, lung metastasis, liver metastasis, breast cancer subtype, race, insurance status, and marital status were independent predictors for BM in IDC patients (Table 4).
Diagnostic nomogram development and validation
A nomogram for predicting the risk of BM in IDC patients was established based on the independent predictors (Fig. 1). ROC analysis showed that the AUCs of the nomogram reached 0.907, demonstrating a better discriminative ability (Fig. 2a). The calibration curve showed high consistency between the observed and predicted results (Fig. 2b). In addition, the DCA indicated that the nomogram had good performance in clinical practice (Fig. 2c).
Prognostic factors for IDC patients with BM
In the training cohort, the univariate Cox proportional hazards regression analysis showed that age, race, primary site, grade, radiotherapy, surgery, chemotherapy, liver metastasis, lung metastasis, brain metastasis, breast cancer subtype, HER2 status, insurance status, and marital status were prognostic factors (all P < 0.05) (Table 5). Then, the multivariate Cox proportional hazards regression analysis was performed. Finally, ten factors, including age, race, grade, surgery, chemotherapy, brain metastases, liver metastases, breast cancer subtypes, insurance status, and marital status, were identified as independent prognostic factors for OS (Table 5).
Prognostic nomogram development and validation
Based on the prognostic factors selected in the training cohort, a nomogram was established to predict the OS for IDC patients with BM (Fig. 3). ROC analysis showed that the AUCs of these nomograms for the 1-, 3-, and 5-year OS reached 0.775, 0.758, and 0.731 in the training cohort; 0.770, 0.773, and 0.753 in the internal validation cohort; and 0.756, 0.764, and 0.767 in the external validation cohort, respectively (Fig. 4a, b, c). The calibration curves of the nomograms showed a strong agreement between actual observations and predictions (Fig. 5). Due to data reasons, the 5-year OS calibration curve for the AHOCMU cohort could not be generated. The clinical application value of the nomogram was evaluated by DCA. As shown in Fig. 6, this nomogram shows a notable positive net benefit over a wide range of death risks, indicating that it has a good clinical utility in predicting the OS for IDC patients with BM. The external validation using the established nomogram in the AHOCMU cohort also demonstrated the high accuracy of the prediction model. Kaplan–Meier survival analysis was performed on the training cohort, internal validation cohort, and external validation cohort, and the results showed that there was an obvious difference in survival rates between the three cohorts (Fig. 7).
Almost all deaths in patients with BC are caused by metastatic disease [4, 5]. Common metastatic sites include bone, lung, liver, and brain, of which bone is the most common [17, 18]. However, unlike the metastases to the lung, liver and brain, BM is generally considered to be less fatal . Once BC patients are diagnosed with BM, the OS decreases dramatically and the median life expectancy decreases to 2–3 years [20, 21]. IDC is the most common pathological type of BC; therefore, it is necessary to identify the risk and prognostic factors for IDC patients with BM to facilitate the early prevention and detection of BM and improve the prognosis for IDC patients with BM.
Currently, there are many studies focused on BM in patients with BC, but there are few studies on IDC patients with BM. Chen et al. reported that in axillary lymph node metastasis, CA125, CA153, ALP, and hemoglobin concentration were independent risk factors for BM in BC patients . Yue Gong et al. determined that age, ethnicity, histology, grade, tumor subtype, extra bone metastasis site, and education level were predictors of BM in BC patients . Other studies have also reported that involvement of more than four axillary lymph nodes at initial diagnosis, primary tumor size greater than 2 cm, estrogen receptor positive and progesterone receptor negative tumors and younger age are risk factors for BM in BC patients [24, 25]. This is similar to the results of our study. In our study, sex, primary site, grade, T stage, N stage, brain metastasis, lung metastasis, liver metastasis, breast cancer subtype, race, insurance status and marital status were significant predictors for BM in IDC patients. Although Zhao et al. established a nomogram model based on gene expression to predict the risk of BM in BC patients, it is not suitable for a wide range of clinical applications and includes all types of BC, which is not conducive to individualized and accurate predictions . To date, no realistic model has been established to predict the risk and prognosis of BM in ICD patients. To address this problem we extracted, screened, and organized specific and relevant prognostic and risk factors of IDC patients with BM and established an intuitive and practical prediction model. This model is beneficial to both the clinician and the individual patient.
It is generally believed that IDC with only metastases to the bone has a better OS prognosis than IDC with bone and visceral metastasis . Previous studies have also found that patients with BM alone had a median survival of approximately two to three times that of patients with additional visceral metastases [28,29,30]. Lobbezoo DJ et al. compared the results of 815 patients with primary or recurrent metastatic BC and found that patients with visceral metastases and patients with multiple metastatic sites had a worse prognosis . Interestingly, our results showed that the presence of brain metastasis and liver metastasis had a significant negative impact on the OS, which is consistent with the above results. In addition, we found that the number of metastatic organ sites also had a significant effect on survival. Previous studies have shown that patients with four metastatic sites are 2.2 times more likely to die than patients with only one metastatic site . We speculate that patients with only BM develop vital organ dysfunction later, so these patients have a higher survival rate than those with both bone and extraosseous metastases. According to previous research, the breast cancer subtype is an independent risk factor for the occurrence of metastasis, and the incidence of BM is highest in BC patients that are HR+/HER2− or HR+/HER2+ [23, 32]. Our results show that patients with HR+/HER2- BC have a higher risk of BM, and patients with Grade 2 BC are more likely to have BM compared to patients with Grades 3 and 4 BC, which is controversial. At present, most people think that once a tumor has distant organ metastasis, it may accelerate the metastasis to other organs, which is consistent with our results . According to our results, chemotherapy had a positive effect on prognosis. Contrary to what we expected, radiotherapy was not a relevant factor for prognosis. Unfortunately, we were unable to compare the effects of different chemotherapy regimens on survival rates because there was no detailed information on chemotherapy strategies in our data.
To facilitate clinical work, we established two nomograms to predict the risk and prognosis for BM in IDC patients. Through calibration, ROC curve and DCA, the nomogram shows great performance, both internally and externally, for predicting the prognosis of IDC patients with BM. These models have better prediction capabilities and higher credibility and can provide references for patient consultations, risk assessment and clinical decision-making. To our knowledge, this is the first population-based model to predict the risk and prognosis of newly diagnosed BM in IDC patients. However, we should acknowledge that this study has some limitations. First, it is a retrospective study and only patients with complete information were included. Therefore, selection bias is likely to exist. Second, some patients with BM have no symptoms, causing the number of newly diagnosed patients with BM to be lower than the actual number. Third, we did not have specific information about systemic treatments, such as endocrine therapy or HER2 targeted therapy. Fourth, since the data in this study were from the SEER database, the nomogram we constructed may not be applicable to IDC patients worldwide.
These nomograms could be used as a supportive graphic tool in IDC to help clinicians distinguish, assess and evaluate the risk and prognosis of IDC with BM. Internal and external validation and application in an independent population demonstrated the satisfactory performance and clinical utility of this predictive model.
Availability of data and materials
The dataset from SEER database generated and/or analyzed during the current study are available in the SEER dataset repository (https://seer.cancer.gov/).
American Joint Committee for Cancer
Infiltrating duct carcinoma
Surveillance, epidemiology and end results
Affiliated Hospital of Chengde Medical University
Receiver operating characteristic
Area under the curve
Decision curve analyses
Siegel RL, Miller KD, Jemal A. Cancer statistics, 2018. CA Cancer J Clin. 2018;68(1):7–30.
Fitzmaurice C, Akinyemiju TF, Al Lami FH, Alam T, Alizadeh-Navaei R, Allen C, Alsharif A-GN, Amini E, Anderson BO, et al. Global, regional, and National Cancer Incidence, mortality, years of life lost, years lived with disability, and disability-adjusted life-years for 29 Cancer groups, 1990 to 2016: a systematic analysis for the global burden of disease study. JAMA Oncol. 2018;4(11):1553–68.
Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68(6):394–424.
DeSantis CE, Fedewa SA, Goding Sauer A, Kramer JL, Smith RA, Jemal A. Breast cancer statistics, 2015: convergence of incidence rates between black and white women. CA Cancer J Clin. 2016;66(1):31–42.
Li Z, Kang Y. Emerging therapeutic targets in metastatic progression: a focus on breast cancer. Pharmacol Ther. 2016;161:79–96.
Yerushalmi R, Tyldesley S, Kennecke H, Speers C, Woods R, Knight B, Gelmon KA. Tumor markers in metastatic breast cancer subtypes: frequency of elevation and correlation with outcome. Ann Oncol. 2012;23(2):338–45.
Kennecke H, Yerushalmi R, Woods R, Cheang MC, Voduc D, Speers CH, Nielsen TO, Gelmon K. Metastatic behavior of breast cancer subtypes. J Clin Oncol Off J Am Soc Clin Oncol. 2010;28(20):3271–7.
Suva LJ, Washam C, Nicholas RW, Griffin RJ. Bone metastasis: mechanisms and therapeutic opportunities. Nat Rev Endocrinol. 2011;7(4):208–18.
McDonald ES, Clark AS, Tchou J, Zhang P, Freedman GM. Clinical diagnosis and Management of Breast Cancer. J Nuclear Med. 2016;57:9S–16S.
Lipton A. Should bisphosphonates be utilized in the adjuvant setting for breast cancer? Breast Cancer Res Treat. 2010;122(3):627–36.
Burke HB. Outcome prediction and the future of the TNM staging system. J Natl Cancer Inst. 2004;96(19):1408–9.
Park YH, Lee SJ, Cho EY, La Choi Y, Lee JE, Nam SJ, Yang JH, Shin JH, Ko EY, Han BK, et al. Clinical relevance of TNM staging system according to breast cancer subtypes. Ann Oncol. 2019;30(12):2011.
Balachandran VP, Gonen M, Smith JJ, DeMatteo RP. Nomograms in oncology: more than meets the eye. Lancet Oncol. 2015;16(4):e173–80.
Molland JG, Donnellan M, Janu NC, Carmalt HL, Kennedy CW, Gillett DJ. Infiltrating lobular carcinoma--a comparison of diagnosis, management and outcome with infiltrating duct carcinoma. Breast. 2004;13(5):389–96.
Liu RZ, Zhao ZR, Ng CS. Statistical modelling for thoracic surgery using a nomogram based on logistic regression. J Thorac Dis. 2016;8(8):E731–6.
Vickers AJ, Elkin EB. Decision curve analysis: a novel method for evaluating prediction models. Med Decis Mak. 2006;26(6):565–74.
Parkin DM, Pisani P, Ferlay J. Estimates of the worldwide incidence of 25 major cancers in 1990. Int J Cancer. 1999;80(6):827–41.
Wu Q, Li J, Zhu S, Wu J, Chen C, Liu Q, Wei W, Zhang Y, Sun S. Breast cancer subtypes predict the preferential site of distant metastases: a SEER based study. Oncotarget. 2017;8(17):27990–6.
Abdel-Rahman O. Population-based validation of the National Cancer Comprehensive Network recommendations for breast cancer staging. Breast Cancer Res Treat. 2018;172(1):231–8.
Li S, Peng Y, Weinhandl ED, Blaes AH, Cetin K, Chia VM, Stryker S, Pinzone JJ, Acquavella JF, Arneson TJ. Estimated number of prevalent cases of metastatic bone disease in the US adult population. Clin Epidemiol. 2012;4:87–93.
Zhang H, Zhu W, Biskup E, Yang W, Yang Z, Wang H, Qiu X, Zhang C, Hu G, Hu G. Incidence, risk factors and prognostic characteristics of bone metastases and skeletal-related events (SREs) in breast cancer patients: a systematic review of the real world data. J Bone Oncol. 2018;11:38–50.
Chen WZ, Shen JF, Zhou Y, Chen XY, Liu JM, Liu ZL. Clinical characteristics and risk factors for developing bone metastases in patients with breast cancer. Sci Rep. 2017;7(1):11325.
Gong Y, Zhang J, Ji P, Ling H, Hu X, Shao ZM. Incidence proportions and prognosis of breast cancer patients with bone metastases at initial diagnosis. Cancer Med. 2018;7(8):4156–69.
Colleoni M, O'Neill A, Goldhirsch A, Gelber RD, Bonetti M, Thürlimann B, Price KN, Castiglione-Gertsch M, Coates AS. Lindtner jet al: identifying breast cancer patients at high risk for bone metastases. J Clin Oncol. 2000;18(23):3925–35.
Wei B, Wang J, Bourne P, Yang Q, Hicks D, Bu H, Tang P. Bone metastasis is strongly associated with estrogen receptor-positive/progesterone receptor-negative breast carcinomas. Hum Pathol. 2008;39(12):1809–15.
Zhao C, Lou Y, Wang Y, Wang D, Tang L, Gao X, Zhang K, Xu W, Liu T, Xiao J. A gene expression signature-based nomogram model in prediction of breast cancer bone metastases. Cancer Med. 2019;8(1):200–8.
Leone BA, Vallejo CT, Romero AO, Machiavelli MR, Pérez JE, Leone J, Leone JP. Prognostic impact of metastatic pattern in stage IV breast cancer at initial diagnosis. Breast Cancer Res Treat. 2017;161(3):537–48.
Harries M, Taylor A, Holmberg L, Agbaje O, Garmo H, Kabilan S, Purushotham A. Incidence of bone metastases and survival after a diagnosis of bone metastases in breast cancer patients. Cancer Epidemiol. 2014;38(4):427–34.
Jacobson AF, Shapiro CL, Van den Abbeele AD, Kaplan WD. Prognostic significance of the number of bone scan abnormalities at the time of initial bone metastatic recurrence in breast carcinoma. Cancer. 2001;91(1):17–24.
Plunkett TA, Smith P, Rubens RD. Risk of complications from bone metastases in breast cancer. implications for management. Eur J Cancer. 2000;36(4):476–82.
Lobbezoo DJ, van Kampen RJ, Voogd AC, Dercksen MW, van den Berkmortel F, Smilde TJ, van de Wouw AJ, Peters FP, van Riel JM, Peters NA, et al. Prognosis of metastatic breast cancer: are there differences between patients with de novo and recurrent metastatic breast cancer? Br J Cancer. 2015;112(9):1445–51.
Metzger-Filho O, Sun Z, Viale G, Price KN, Crivellari D, Snyder RD, Gelber RD, Castiglione-Gertsch M, Coates AS, Goldhirsch A, et al. Patterns of recurrence and outcome according to breast cancer subtypes in lymph node-negative disease: results from international breast cancer study group trials VIII and IX. J Clin Oncol. 2013;31(25):3083–90.
Xie J, Ying YY, Xu B, Li Y, Zhang X, Li C. Metastasis pattern and prognosis of male breast cancer patients in US: a population-based study from SEER database. Ther Adv Med Oncol. 2019;11:1758835919889003.
We would like to thank all the staff in Department of Spine Surgery, Affiliated Hospital of Chengde Medical University for their contribution on our research.
We received no external funding for this study.
Ethics approval and consent to participate
We received permission to access the research data file in the SEER program from the National Cancer Institute, US (reference number 15260-Nov2018). Approval was waived by the local ethics committee, as SEER data is publicly available and de-identified.
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
Huang, Z., Hu, C., Liu, K. et al. Risk factors, prognostic factors, and nomograms for bone metastasis in patients with newly diagnosed infiltrating duct carcinoma of the breast: a population-based study. BMC Cancer 20, 1145 (2020). https://doi.org/10.1186/s12885-020-07635-1
- Breast cancer
- Infiltrating duct carcinoma
- Bone metastasis