- Research article
- Open Access
- Open Peer Review
Nomogram based on homogeneous and heterogeneous associated factors for predicting bone metastases in patients with different histological types of lung cancer
BMC Cancervolume 19, Article number: 238 (2019)
The purpose of the present study was to characterize the prevalence, associated factors, and to construct a nomogram for predicting bone metastasis (BM) with different histological types of lung cancer.
Patients and methods
This study was a descriptive study that basing on the invasive lung cancer patients diagnosed between 2010 and 2014 in Surveillance, Epidemiology, and End Results program. A total of 125,652 adult patients were retrieved. Logistic regression analysis was conducted to investigate homogeneous and heterogeneous factors for BM occurrence. Nomogram was constructed to predict the risk for developing BM and the performance was evaluated by the receiver operating characteristics curve (ROC) and the calibration curve. The overall survival of the patients with BM was analyzed using the Kaplan–Meier method and the survival differences were tested by the log-rank test.
A total of 25,645 (20.9%) were reported to have BM, and the prevalence in adenocarcinoma, squamous cell carcinoma, small cell lung cancer (SCLC), large cell lung cancer (LCLC), and non-small cell lung cancer/not otherwise specified lung cancer (NSCLC/NOS) were 24.4, 12.5, 24.7, 19.5 and 19.4%, respectively, with significant difference (P < 0.001). Male gender, more metastatic sites and lymphatic metastasis were positively associated with BM in all lung cancer subtypes. Larger tumor size was positively associated with BM in all the lung cancer subtypes except for NSCLC/NOS. Poorly differentiated histology was positively associated with adenocarcinoma, squamous cell carcinoma and NSCLC/NOS. The calibration curve and ROC curve exhibited good performance for predicting BM. The median survival of the bone metastatic lung cancer patients was 4.00 (95%CI: 3.89–4.11) months. With the increased number of the other metastatic sites (brain, lung and liver metastasis), the survival significantly decreased (p < 0.001).
Different lung cancer histological subtypes exhibited distinct prevalence and homogeneity and heterogeneity associated factors for BM. The nomogram has good calibration and discrimination for predicting BM of lung cancer.
Lung cancer is the most common carcinoma and the leading cause of cancer death globally . Patients with metastatic disease continue to exhibit a poor prognosis . Bone metastasis (BM) was reported to occur in 15–40% of lung cancer patients. The median survival time of patients with BM was reported to be less than one year [3,4,5].
Early diagnosis and intervention in patients with BM could significantly influence the survival rate of patients [6, 7]. However, the Lung Cancer National Comprehensive Cancer Network (NCCN) screening guidelines do not recommend performing routine assessment or continued reassessment for BM by skeletal imaging in asymptomatic patients . Previous studies reported a series of associated factors for BM, which provided the basis for predicting the BM risk [9, 10]. To promptly perform metastatic screening, a predictive nomogram based on the clinicopathologic features of lung cancer is warranted.
With the development of tumor homogeneity/heterogeneity theory, different histological types of lung cancer are now recognized to exhibit distinct prevalence of BM [4, 11]. We also believe that different histological types of lung cancer are associated with different factors for BM occurrence. Further study identifying the homogeneous and heterogeneous BM associated factors could help physicians to specifically identify the BM risk for different types of lung cancer and tailor targeted preventive treatment strategies.
The purpose of the present study was to characterize the prevalence and associated factors for BM in patients with different histological types of lung cancer using the Surveillance, Epidemiology, and End Results (SEER) database. Meanwhile, the clinical factor based nomogram was built to predict the BM risk and potentially guide the BM screening.
Cancer is a reportable disease in every state of the United States. The data in the SEER database does not require informed patient consent. The present study was complied with the 1964 Helsinki Declaration and its later amendments or comparable ethical standards. This study used previously collected deidentified data, which was deemed exempt from review by the Ethics Board of the Tianjin Medical University Cancer Institute and Hospital.
It was a population based descriptive study and the data were abstracted from the SEER 18 registries research database, comprising approximately 30% of the total US population. As the data for metastatic sites of bone, liver, lung and brain were not collected until 2010, lung cancer patients who were diagnosed between 2010 and 2014 were included in the present study. SEER*Stat Software version 8.3.4 (https://seer.cancer.gov/seerstat/) (Information Management Service, Inc. Calverton, MD, USA) was used to generate the case listing.
The inclusion criteria were as follows: (1) aged 18–79 years; (2) diagnosed as the first and only malignant cancer; (3) only one primary site; (4) American Joint Committee on Cancer (7th edition) stage I-III; and (5) diagnosis not obtained from a death certificate or an autopsy. Patients diagnosed before 2010 were excluded because the SEER did not record BM data until 2010. The lung cancer was histologically classified as adenocarcinoma, squamous cell carcinoma, small cell lung cancer (SCLC), large cell lung cancer (LCLC) and non-small cell lung cancer or not otherwise specified lung cancer (NSCLC/NOS) based on the International Classification of Diseases for Oncology, 3rd Edition codes. The flow-chart for the population selection was shown in Fig. 1.
Quantitative data were described as mean ± standard deviation (SD) and the difference between groups were analysed by student’s t-test. Categorical data was presented as number and the percentage (N, %). Pearson chi-square test was used to evaluate the difference between categorical variables. The univariable and multivariable logistic regression model were conducted to determine the associated factors of BM by different histological types of lung cancer. Factors with a P-value less than 0.05 were incorporated into the multivariable regression model. The final model selection was performed by a backward step-down selection process using the Akaike information criterion.
A nomogram was also formulated based on the results of multivariable logistic analysis using the rms package in R version 3.4.1 (R Foundation for Statistical Computing, Vienna, Austria; www.r-project.org). The performance of the nomogram was evaluated by the receiver operating characteristics (ROC).
To evaluate the calibration of the nomogram, a regression smoothing method was used to produce the calibration plots by bootstrapping with 1000 resamples, where the relationship between the observed and predicted probabilities of BM was described graphically.
Randomly splitting and the temporal splitting method were used to evaluate the stability of the nomogram. An ROC curve was constructed to evaluate the performance of the construction and validation model, and the difference in the aura under the curve (AUC) was tested by DeLong’s test. Statistically significant levels were two-tailed and set at P < .05. Statistical analyses were performed using the statistical Package for the Social Sciences (SPSS) version 23.0 software package for Windows (SPSS, Inc.).
The overall survival of the lung cancer patients with bone metastasis were analyzed using the Kaplan–Meier method and the survival difference between different metastatic sites, treatment regimens of surgery, radiotherapy and chemotherapy were tested by the log-rank test.
The study consisted of 125,652 patients, including 54,957 (43.7%) cases of adenocarcinoma, 27,435 (21.8%) of squamous cell carcinoma, 18,795 (15.0%) of SCLC, 2308 (1.8%) of LCLC and 22,157 (17.6%) of NSCLC/NOS lung cancer patients. Among them, 67,216 (53.5%) were male, and 58,436 (46.5%) were female, mean age was 64.84 ± 9.17 years (Table 1).
Prevalence of BM
After excluding the patients with unknown BM information, 25,645 (20.9%) were reported to exhibit BM. When stratified by histological subtype, the prevalence of BM in adenocarcinoma, squamous cell carcinoma, SCLC, LCLC, and NSCLC/NOS were 24.4, 12.5, 24.7, 19.5 and 19.4%, respectively. The prevalence rates of BM in adenocarcinoma and SCLC were higher than those in the other lung cancer histological types (P < 0.001), and the prevalence of BM in squamous cell lung cancer was lowest among all types (P < 0.001) (Additional file 1: Table S1).
Univariable logistic regression analysis
Advanced age [odds ratio (OR) = 0.81; 95% confidence interval (CI): 0.79–0.83; P < 0.001] and insured status (OR = 0.79; 95% CI: 0.74–0.85; P < 0.001) were negatively associated with BM. However, male gender (OR = 1.27; 95% CI: 1.23–1.31; P < 0.001), race (OR = 1.06; 95% CI: 1.03–1.08; P < 0.001), married status (OR = 0.95; 95% CI: 0.92–0.99; P = 0.005), household income (OR = 1.07; 95% CI: 1.05–1.10; P < 0.001), more metastatic sites (OR = 2.74; 95% CI: 2.69–2.80; P < 0.001), larger tumor size (OR = 1.29; 95% CI: 1.27–1.32; P < 0.001), poor histological differentiation (OR = 1.70; 95% CI: 1.65–1.75; P < 0.001), and more lymphatic metastasis (OR = 1.61; 95% CI: 1.59–1.63; P < 0.001) were positively associated with BM.
Additionally, diverse histological subtypes of the lung cancer exhibit a differential risk for developing BM. Compared with adenocarcinoma, squamous cell carcinoma (OR = 0.44; 95% CI: 0.43–0.46), LCLC (OR = 0.75; 95% CI: 0.68–0.84) and NSCLC/NOS (OR = 0.75; 95% CI: 0.72–0.78) were negatively associated with BM, whereas SCLC exhibited no difference in the risk of BM relative to adenocarcinoma (OR = 1.02; 95% CI: 0.98–1.06). Moreover, subgroup analysis suggested that the associated factors for bone metastases were not consistent across all the histological subtypes of lung cancer (Additional file 2: Table S2).
Multivariable logistic regression analysis
When conducting multivariable logistic regression analysis, the patients with unknown or missing information concerning sex, metastatic site, tumor size, histological type and lymphatic metastasis were excluded. Finally, 68,161 patients with lung cancer were included (Fig. 1). For all the lung cancer patients, multivariable logistic regression suggested male gender, histological differentiation, more lymphatic metastasis, and lung cancer subtype were independent factors associated with BM (Additional file 3: Table S3). When stratified by the histological subtypes of lung cancer, results revealed that male gender, more metastatic sites, larger tumor size, poor histological differentiation and more lymphatic metastasis were all positively associated with BM in adenocarcinoma. The associated factors and OR with 95% CI for squamous cell carcinoma, SCLC, LCLC and NSCLC/NOS lung cancer are presented in Table 2.
Different lung cancer histological subtypes exhibited homogeneity and heterogeneity for the factors associated with bone metastases. Male gender, more metastatic sites and more lymphatic metastasis were positively associated with bone metastases among all lung cancer histological subtypes. However, larger tumor size was not associated with BM in NSCLC/NOS. Poorly differentiated histology was positively associated with adenocarcinoma, squamous cell carcinoma and NSCLC/NOS lung cancer but not with SCLC or LCLC (Fig. 2).
Performance of the nomogram for predicting BM
The prediction nomogram that integrated all significant factors for BM in different lung cancer histologic types is presented in Fig. 3. The calibration curve revealed good agreement between the predicted and observed probabilities for BM in different histological types of lung cancer. Moreover, the ROC curve of the nomogram exhibited good discrimination for predicting BM, and the AUC of the nomogram in adenocarcinoma, squamous cell carcinoma, SCLC, LCLC and NSCLC/NOS lung cancer were 80.3% (95% CI: 79.6–80.9%), 78.1% (95% CI: 76.8–79.4%), 70.8% (95% CI: 69.8–71.8%), 75.1% (95% CI: 72.3–77.8%) and 80.2% (95% CI: 79.2–81.3%), respectively.
Validation of the nomogram
For adenocarcinoma of lung cancer, the random splitting method revealed that the AUC values for the construction and validation model were 80.2% (95% CI: 79.4–81.1%) and 80.2% (95% CI: 79.0–81.4%), respectively, with no significant difference (D = 0.03; P = 0.97) (Fig. 4a). The nomogram for predicting BM was also stable in squamous cell carcinoma (D = 0.67; P = 0.50), SCLC (D = -0.37; P = 0.71), LCLC (D = 1.16; P = 0.25) and NSCLC/NOS (D = 1.14; P = 0.25) (Fig. 4b-e).
Temporal splitting analysis suggested that the AUC in the adenocarcinoma model to predict BM was 79.7% (95% CI: 78.9–80.5%), which was significantly lower than 280.9–83.8%) for patients diagnosed in 2014 (D = -3.21; P = 0.001) (Fig. 4f). However, there was no difference between the constructed and validation model for predicting BM in squamous cell carcinoma (D = 1.23; P = 0.22), SCLC (D = -1.26; P = 0.21), LCLC (D = -0.32; P = 0.75) or NSCLC/NOS (D = -1.77; P = 0.08), suggesting the stability of the predictive nomogram (Fig. 4g-j).
Survival analysis for patients with BM
The median survival of the bone metastatic lung cancer patients was 4.00 (95%CI: 3.89–4.11) months. When stratified by different histological types, the median survival time for adenocarcinoma, squamous cell carcinoma, small cell lung cancer (SCLC), large cell lung cancer (LCLC), and non-small cell lung cancer/not otherwise specified lung cancer (NSCLC/NOS) were 5.00 (95%CI: 4.83–5.17), 3.00 (95%CI: 2.80–3.20), 6.00(95%CI: 5.72–6.29), 3.0 (95%CI: 2.47–3.53) and 3.0(95%CI: 2.80–3.21) months, respectively with significant difference (p < 0.001) (Fig. 5a). Additionally, results also showed with the increased number of the other metastatic sites (brain, lung and liver metastasis), the survival significantly decreased, and the median survival for 0, 1 2 and 3 metastatic sites were 5.00 (95% CI: 4.82–5.18), 4.00(95% CI: 3.84–4.16), 3.00 (95% CI: 2.77–3.23) and 3.00 (95% CI: 2.64–3.37) months, respectively (Fig. 5b).
The present study utilized the SEER database and determined the prevalence and associated factors for BM in patients with different histological types of lung cancer. Meanwhile, the quantitative prediction nomogram for different histological types of lung cancer were built.
BM was less frequent in squamous cell carcinoma. However, adenocarcinoma and SCLC exhibited an association with BM, and adenocarcinoma accounted for more than 50% of all the lung cancer patients with BM, a finding that is partly consistent with a hospital-based study . Based on the large population data set provided by SEER, the present study further concluded that small cell lung cancer and lung adenocarcinoma exhibit higher risks of BM than do the other histological subtypes. Previous studies reported similar results. After evaluating 413 patients who were diagnosed with lung cancer, Oliveira et al. found that adenocarcinoma was associated with a higher risk for developing BM, whereas squamous cell cancer was associated with a lower risk of BM . Another study also reported that the most common histological type accompanied by BM was adenocarcinoma . A high prevalence of BM in lung cancer and different BM prevalence rates across various histological subtypes may partly reflect the homogeneity and heterogeneity of lung cancer. Additionally, a cancer registry-based study conducted in Sweden reported that the prevalence of BM in adenocarcinoma was approximately 39%, which was higher than our results. Moreover, although BM was most prevalent among adenocarcinoma, SCLC exhibited a significantly lower prevalence of BM (25%) than did the other lung cancer subtypes, a finding that was different from ours . Accordingly, the homogeneity and heterogeneity of lung cancer may differ according to ethnicity.
In the present study, we also found that different lung cancer histological subtypes exhibited homogeneity and heterogeneity with respect to the factors associated with BM. Regarding homogeneity, we identified three metastatic associated factors for all types of lung cancer: male gender, extrapulmonary metastatic site and lymphatic metastasis.
Among the present cohort, male gender was an associated factor for BM, a finding that was independent of the clinical features of lung cancer. To our knowledge, this is the first report describing sex as a risk factor for BM in lung cancer patients. Lung cancer, unlike breast cancer, is a non-hormone-dependent tumor. The reason that sex may significantly influence metastasis in non-hormone-dependent tumors remains unknown. Further research to determine potential explanations is needed.
The number of metastatic sites was also reported to be one of the metastatic associated factors for all types of lung cancer. According to our analysis, an increase in the number of metastatic sites was associated with an increased prevalence of BM in lung cancer patients. At the same time, for all histological subtypes, a higher grade of lymphatic involvement was associated with an increasing prevalence of BM in lung cancer patients. The aforementioned homogeneous associated factors may help in the surveillance of BM in lung cancer patients. Thus, physicians should focus on their lung cancer patients with these associated factors.
The advantages of early metastasis detection are as follows: 1) less toxic therapy [14, 15]; 2) SREs can be prevented through timely bone-targeted therapy ; and 3) improved performance status after therapy . Thus, prompt metastatic screening is necessary. To make an early diagnosis and improve the survival rate of cancer patients, efforts have been made to identify the optimal screening time and method. In the latest studies, bone biomarkers and breast osteoblast-like cells were identified as potential strategies for clinical metastatic screening and early diagnosis [18, 19]. However, the aforementioned strategies require extra techniques and equipment support, which may decelerate the process of clinical screening application. Imaging remains a reliable and accepted strategy for metastatic screening and early diagnosis . To prevent unnecessary radiation exposure and cost, in the present study, a histological type-based prediction nomogram was constructed. A quantitative metastatic risk could be generated utilizing the patient’s sex, lymph node metastasis, extrapulmonary metastatic site number, histological differentiation, and tumor size. We concluded that physicians could perform metastatic screening for their lung cancer patients with a high BM risk.
Inevitably, the present study has several limitations. First, in the present study, only the presence/absence of BM based on the initial diagnosis was analyzed. Disease recurrence or subsequent sites of disease were not provided in the SEER database. Thus, the actual rate of BM in patients with lung cancer might be underestimated. Second, the information on performance status, smoking status, blood type, and body mass index were not provided in the SEER database, which may affect the precision of the predictive nomogram. Third, the present predictive nomogram lacked external validation, more studies are needed in future. Fourth, the site of BM was not recorded in the SEER database, thus it cannot be further identified or analyzed in the study.
Despite the limitations, the present study provided insight into the epidemiology of BM in patients with newly diagnosed lung cancer, as recorded by the SEER database. The prevalence of BM in lung cancer was 20.9% and different lung cancer histological types showed homogeneous associated factors (male gender, more metastatic sites and more lymphatic metastasis) and heterogeneous associated factors (tumor size and histology differentiation) for BM. The median survival of the bone metastatic lung cancer patients was relatively low. The histological types significantly affect the median survival time of the lung cancer patients with BM. The nomogram has good performance for predicting the BM development in different histological types of lung cancer and the imaging of the skeletal system should be considered to lung cancer patients with a high BM risk.
aura under the curve
Large cell lung cancer
National Comprehensive Cancer Network
Non-small cell lung cancer/not otherwise specified lung cancer
Receiver operating characteristics curve
Small cell lung cancer
Surveillance, Epidemiology, and End Results
Malvezzi M, Bertuccio P, Rosso T, et al. European cancer mortality predictions for the year 2015: does lung cancer have the highest death rate in EU women? Ann Oncol. 2015;26:779–6.
Australian Institute of Health and Welfare. Cancer in Australia: actual incidence data from 1982 to 2013 and mortality data from 1982 to 2014 with projections to 2017. Asia Pac J Clin Oncol. 2018;14:5–15.
Al Husaini H, Wheatley-Price P, Clemons M, Shepherd FA. Prevention and management of bone metastases in lung cancer: a review. J Thorac Oncol. 2009;4:251–9.
Zhou Y, Yu QF, Peng AF, Tong WL, Liu JM, Liu ZL. The risk factors of bone metastases in patients with lung cancer. Sci Rep. 2017;7:8970.
Kuchuk M, Kuchuk I, Sabri E, Hutton B, Clemons M, Wheatley-Price P. The incidence and clinical impact of bone metastases in non-small cell lung cancer. Lung Cancer. 2015;89:197–202.
Bhatia R, Ravulapati S, Befeler A, Dombrowski J, Gadani S, Poddar N. Hepatocellular carcinoma with bone metastases: incidence, prognostic significance, and management-single-center experience. J Gastrointest Cancer. 2017;48:321–5.
Makino H, Nishio S, Tsubamoto H, et al. Treatment and prognosis of bone metastasis from cervical cancer (KCOG-G1202s). J Obstet Gynaecol Res. 2016;42:701–6.
Network NCC. NCCN Clinical practice guidelines in oncology. Lung Cancer. 2016; https: //www.nccn.org/professionals/physician_gls/pdf /lung.pdf.
Oliveira MB, Mello FC, Paschoal ME. The relationship between lung cancer histology and the clinicopathological characteristics of bone metastases. Lung Cancer. 2016;96:19–24.
Conen K, Hagmann R, Hess V, Zippelius A, Rothschild SI. Incidence and predictors of bone metastases (BM) and skeletal-related events (SREs) in small cell lung Cancer (SCLC): a Swiss patient cohort. J Cancer. 2016;7:2110–6.
Zhang L, Gong Z. Clinical characteristics and prognostic factors in bone metastases from lung Cancer. Med Sci Monit. 2017;23:4087–94.
Sugiura H, Yamada K, Sugiura T, Hida T, Mitsudomi T. Predictors of survival in patients with bone metastasis of lung cancer. Clin Orthop Relat Res. 2008;466:729–36.
Riihimäki M, Hemminki A, Fallah M, et al. Metastatic sites and survival in lung cancer. Lung Cancer. 2014;86:78–84.
Curigliano G, Mayer EL, Burstein HJ, Winer EP, Goldhirsch A. Cardiac toxicity from systemic cancer therapy: a comprehensive review. Prog Cardiovasc Dis. 2010;53:94–104.
Barth RJ Jr, Gibson GR, Carney PA, Mott LA, Becher RD, Poplack SP. Detection of breast cancer on screening mammography allows patients to be treated with less-toxic therapy. AJR Am J Roentgenol. 2005;184:324–9.
LeVasseur N, Clemons M, Hutton B, Shorr R, Jacobs C. Bone-targeted therapy use in patients with bone metastases from lung cancer: a systematic review of randomized controlled trials. Cancer Treat Rev. 2016;50:183–93.
Abrahm JL1, Banffy MB, Harris MB. Spinal cord compression in patients with advanced metastatic cancer: "all I care about is walking and living my life". JAMA. 2008;299:937–946.
Brown J, Rathbone E, Hinsley S, et al. Associations between serum bone biomarkers in early breast Cancer and development of bone metastasis: results from the AZURE (BIG01/04) trial. J Natl Cancer Inst. 2018;110(8):871–9.
Scimeca M, Antonacci C, Toschi N, et al. Breast osteoblast-like cells: a reliable early marker for bone metastases from breast Cancer. Clin Breast Cancer. 2017. https://doi.org/10.1016/j.clbc.2017.11.020 Epub ahead of print.
Guo X, Zhang C, Guo Q, et al. The homogeneous and heterogeneous risk factors for the morbidity and prognosis of bone metastasis in patients with prostate cancer. Cancer Manag Res. 2018;10:1639–46.
The present study was sponsored by the Natural Science Foundation of China (81702161, 81602363), the Russian Foundation of Basic Research (15-29-01338), the Natural Science Foundation of the Tianjin Science and Technology Committee of China (17JCQNJC11000), the Natural Science Foundation of Tianjin Medical University (2016KYZQ10), the Doctor Start-up Grant of Tianjin Medical University Cancer Institute and Hospital (B1711) and Top talent training program of the first affiliated hospital of PLA Army Medical University [SWH2018BJKJ-12].
Availability of data and materials
The data were abstracted from the Surveillance, Epidemiology, and End Results (SEER) database. This is an open database. (https://seer.cancer.gov).
Ethics approval and consent to participate
This study used previously collected deidentified data, which was deemed exempt from review by the Ethics Board of the Tianjin Medical University Cancer Institute and Hospital.
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.
Table S1. Difference in the prevalence of bone metastasis in different histological subtypes of lung cancer. (PDF 252 kb)
Table S2. Univariate logistic regression for the presence of bone metastases at diagnosis of different subtypes of lung cancer. (PDF 298 kb)
Table S3. Multivariable logistic regression for the presence of bone metastases at diagnosis of lung cancer. (PDF 228 kb)