Abstract
Objective
Patients with bone metastasis (BM) of small cell lung cancer (SCLC) have a poor prognosis. We aimed to identify predictors and prognostic factors in patients with BM of SCLC and construct nomograms to predict BM.
Methods
We retrospectively analyzed 18,187 cases from the Surveillance, Epidemiology, and End Results database reported between 2010 and 2016. Differences in overall survival (OS) and cancer-specific survival (CSS) were evaluated after propensity score matching. Independent predictors for BM and prognostic factors for patients with BM of SCLC were determined using univariate and multivariate regression analyses. Two nomograms were constructed and evaluated using C-statistics.
Results
BM was observed in 4014 (22.07%) patients. Kaplan–Meier survival analysis revealed significant differences between BM and non-BM groups. The median OS for patients with and without BM was 6 and 7 months, respectively. The median CSS for patients with and without BM was 9 and 13 months, respectively. Age, sex, tumor size, N stage, chemotherapy, surgery, radiotherapy, and liver/brain/lung metastases were related to BM and independent prognostic factors for OS and CSS. Diagnostic and prognostic nomograms were generated.
Conclusion
Our nomograms predicted the incidence of BM and the 5-month survival rate of patients with SCLC and BM.
Keywords
Introduction
Cancer of the lung and bronchus ranks second among all malignancies. Although the mortality associated with this cancer has declined in recent years, it remains the leading cause of death among all cancer types. According to the 2020 Cancer Statistics, the American Cancer Society estimated that the number of lung and bronchus cancer cases in the United States in 2020 increased by 116,300 in men and 112,520 in women. 1 Previous reports indicated that lung cancer is the most commonly diagnosed malignancy in both sexes (combined), accounting for 11.6% of all cancer cases. It is also the leading cause of cancer death (18.4% of the total cancer deaths) worldwide. 2 Among the various subtypes of lung cancer, small cell lung cancer (SCLC) is an aggressive pathological type. Although it is highly sensitive to initial chemotherapy and radiotherapy treatment, most patients have poor outcomes due to recurrence and disease progression. 3
Bone metastasis (BM) is a common disease progression outcome and occurs in a variety of malignant tumors, including lung cancer, prostate cancer, breast cancer, and others. It has been previously reported that the deterioration of BM results in hypercalcemia, severe bone pain, pathological fracture, spinal cord compression, and other skeletal-related events (SREs). 4 Severe SREs further lead to loss of bone function, thereby reducing the quality of life of patients. Additionally, BM significantly reduces the survival time of patients and is almost always incurable. 5
Early diagnosis and treatment of BM can significantly reduce the occurrence and development of SREs, subsequently improving the survival rate of patients.6,7 Approximately 40% of patients with lung cancer develop BM at an advanced stage of the disease. 8 However, in accordance with the screening guidelines of the National Comprehensive Cancer Network, invasive examination or skeletal imaging is not routinely recommended for asymptomatic patients. Therefore, it is particularly important to identify predictors of BM for lung cancer and construct an early screening model. Thus, the present study analyzed epidemiological data for BM in patients with SCLC from the Surveillance, Epidemiology, and End Results (SEER) database, the largest publicly available cancer database, to identify factors affecting the occurrence of BM and prognosis of patients with BM of SCLC. An epidemiological prediction score model and survival model of BM were generated to predict the incidence of BM in patients with SCLC and the 5-month survival rate of patients with BM of SCLC.
Materials and methods
Patient selection
The data from 1973 to 2016 included in the present study were downloaded from SEER*Stat software (Version 8.3.6.1 https: //seer.cancer.gov/data/access.html). Because this was a retrospective analysis of data from the SEER database, medical ethics review and patient consent were not applicable. All procedures performed in the studies were in line with the Helsinki Declaration and its later amendments and other comparable ethical standards. The reporting of the present study conforms to the TRIPOD guidelines. 9
We defined the inclusion criteria of patients with lung cancer as follows: (1) patients who were histologically diagnosed with SCLC (The International Classification of Diseases for Oncology third edition (ICD-O-3) was used to identify SCLC by site codes [8002, 8041, 8043, 8144, 8145]) from 2010 to 2016; (2) the staging of lymph nodes followed the 7th edition of the American Joint Committee on Cancer; (3) all relevant data were available, such as age at the initial diagnosis, race, sex, N-stage, primary site, histological type (ICD-O-3), bone/liver/lung metastasis, chemotherapy recode, radiation recode, operation information, survival duration (months), cause-specific death classification, and survival status. The following conditions were used as the exclusion criteria: (1) diagnosed with autopsies or death certificates and (2) patients with incomplete information. Finally, only eligible patients were included in this cohort study. Patients with SCLC and BM were defined as the BM group, and those without BM were defined as the non-BM group.
Statistical analyses and nomogram generation
All statistical analyses were conducted with IBM SPSS Statistics for Windows, Version 26.0 (IBM Corp., Armonk, NY, USA) and R (www.r-project.org). Categorical variables in the BM and non-BM groups were compared by the chi-square test. Univariate logistic analysis was applied to determine the factors related to BM. Significant variables were included in the multivariate logistic regression analysis to determine the independent risk factors of BM in patients with SCLC, and the results of multivariate logistic regression were presented in forest plots using R software. Regarding prognostic factors, the significant variables in the univariate Cox proportion hazard regression analysis were applied to construct a multivariate Cox proportional hazards model to determine the independent prognostic factors of patients with BM of SCLC. In this study, a two-sided p < 0.05 was considered to indicate statistical significance.
Propensity score matching (PSM) analyses were used to balance differences in the clinical baseline characteristics between the BM and non-BM groups. PSM was performed using IBM SPSS Statistics for Windows, Version 26.0. The two groups were matched based on a ratio of 1:1. The overall survival (OS) and cancer-specific survival (CSS) were analyzed using Kaplan–Meier curves, and the differences were compared with the log-rank test.
A nomogram was used to predict the incidence of BM and the OS of patients with BM of SCLC. Meaningful and commonly used indicators in multivariate analyses were incorporated into the construction of nomograms using R. The bootstrap method was used for internal validation, and the c-index for each nomogram was calculated.
Results
Characteristics of patients with SCLC
A total of 268,251 patients were newly diagnosed with SCLC between 2010 and 2016. Based on the selection criteria, a total of 18,187 patients were included in this study. The characteristics of the patients before PSM are presented in Table 1. Regarding treatment methods, only 4.27% of all patients underwent surgery, whereas 72.58% and 49.82% received chemotherapy and radiotherapy, respectively. In terms of metastatic sites, BM was reported in 22.07% of patients, whereas 15.70% exhibited brain metastasis.
Characteristics of patients with SCLC before PSM.
BM, bone metastasis; SCLC, small cell lung cancer; PSM, propensity score matching.
PSM analysis
To avoid any bias from other factors, such as age and sex, PSM was used to balance 12 pairs of clinical and histopathological characteristics during the evaluation of prognosis. After 1:1 matching, a total of 7920 cases were matched in the BM and non-BM groups. As shown in Table 2, the differences in the clinical and histopathological characteristics between the two groups after PSM were balanced (p > 0.05). OS and CSS were calculated using the data collected before and after PSM. Before PSM, the median OS for patients with and without BM was 6 months [95% confidence interval (CI): 5.7–6.3] and 9 months (95% CI: 8.8–9.2), respectively. The median CSS for patients with and without BM was 6 months (95% CI: 5.7–6.3) and 9 months (95% CI: 8.7–9.3), respectively. Importantly, the OS and CSS of patients with BM were significantly different from patients without BM after PSM (p < 0.001). The median OS for patients with BM and without BM was 6 (95% CI: 5.7–6.3) and 7 (95% CI: 6.7–7.3) months, respectively. After PSM, the median CSS for patients with BM and without BM was 9 (95% CI: 8.6–9.4) and 13 (95% CI: 12.3–13.6) months, respectively. No, statistically significant differences were observed in OS between the patients with BM before and after PSM. The Kaplan–Meier curves for OS and CSS after PSM are shown in Figure 1.
Characteristics of patients with SCLC after PSM.
BM, bone metastasis; SCLC, small cell lung cancer; PSM, propensity score matching.

Kaplan–Meier curves of overall survival (a) and cancer-specific survival (b) for patients stratified by the presence (N = 4014) or absence (N = 14,173) of BM.
Predictors of BM in patients with SCLC
For univariate logistic regression analysis, age, race, sex, tumor size, N stage, surgery, chemotherapy, radiation, brain metastasis, liver metastasis, and lung metastasis were identified to be closely related to the risk of BM. The related factors were further used for multivariate analysis, which indicated no differences between the two groups in patients 55 to 74 years old (Table 3). Importantly, age ≥85 years (p < 0.001), Black race (p = 0.018), other races (p = 0.041), men (p < 0.001), larger tumor size (p < 0.001), higher N stage (p < 0.001), surgery (p < 0.001), chemotherapy (p < 0.001), radiation (p < 0.001), brain metastasis (p < 0.001), liver metastasis (p < 0.001), and lung metastasis (p < 0.001) were identified as independent predictors of BM in newly diagnosed patients with SCLC. For better presentation of the results, forest plots were used (Figure 2).
Univariate logistic regression analyses of factors associated with BM.
*Other: American Indian/Alaska Native, Asian/Pacific Islander.
BM, bone metastasis; HR, hazard ratio; CI, confidence interval.

Multivariable logistic regression analyses of factors associated with BM.
Prognosis of patients with BM of SCLC
The prognostic factors of OS and CSS were predicted using univariate and multivariate Cox proportional hazard regression analyses (Table 4). Importantly, sex, tumor size, N stage, surgery performed, chemotherapy, radiation therapy, and bone/brain/liver/lung metastases were found to be statistically significant in the univariate analysis for both OS and CSS (p < 0.05). As expected, age >65 years was statistically significant in the univariate analysis of OS. In the case of CSS, patients aged 65 to 74 years and those older than 85 years exhibited a statistically significant difference (p < 0.05). In contrast, race was identified as a risk factor for CSS but not OS.
Cox-regression of univariate and multivariate analyses associated with the OS and CSS of patients with BM of SCLC.
*Other: American Indian/Alaska Native, Asian/Pacific Islander.
BM, bone metastasis; SCLC, small cell lung cancer; HR, hazard ratio; CI, confidence interval; CSS, cancer-specific survival; OS, overall survival.
The multivariate Cox proportional hazard regression analysis revealed that factors, including men, Black, other races, N2 stage, N3 stage, N4 stage, surgery, chemotherapy, radiation, brain metastasis, liver metastasis, and lung metastasis were independent prognostic factors for both OS and CSS. However, age >65 years was identified as an independent prognostic factor in the OS group, whereas Black and other races were independent prognostic factors in the CSS group.
Diagnostic and prognostic nomograms for patients with BM of SCLC
The nomograms for diagnosis and prognosis were established based on the results of the multivariate logistic regression analysis and Cox proportional hazards regression model (Figure 3 and Figure 4). These included multiple factors, such as age, race, sex, tumor size, N stage, surgery, chemotherapy, radiation, brain metastasis, liver metastasis, and lung metastasis. In the internal validation set, the c-index of the diagnosis nomogram was 0.745, and the c-index of the prognosis nomogram was 0.565.

A nomogram for predicting the risk of BM in patients with SCLC. Other includes American Indians/Alaska Natives and Asian/Pacific Islanders. The size of the cyan square and grey area under the curve on the total points axis represents the sample size. The population distribution of patients with SCLC are shown.

A nomogram for predicting the prognosis of patients with BM of SCLC. Other includes American Indians/Alaska Natives and Asian/Pacific Islanders. The size of the cyan square and grey area under the curve on the total points axis represents the sample size. The population distribution of patients with SCLC are shown.
Using the diagnosis nomogram, we included an example of an 86-year-old white woman who presented with a 4-cm tumor and a positive supraclavicular lymph node. The patient had not received radiotherapy, chemotherapy, or surgical treatment but exhibited liver metastasis and bilateral lung lesions. The total points were 372, and the probability of BM was approximately 48.55%. Similarly, using the prognosis nomogram, we included an example of a 76-year-old Black man with BM of SCLC who exhibited a 6-cm tumor and stage N2. The patient presented with no previous history of surgical treatment and radiotherapy. However, chemotherapy was previously administered. He showed no signs of liver, brain, or lung metastases. Here, the total score was ∼859, and the probability that the patient would survive for >5 months (median OS) was calculated to be 66.4%.
Discussion
The present study focused on the malignant progression and prognosis of lung cancer. Although several studies have been conducted on non-small cell lung cancer (NSCLC) in recent years, and the treatment efficacy of patients with NSCLC has improved significantly, the prognosis of patients with SCLC remains poor, mainly due to limited treatment options. SCLC accounts for ∼15% of total lung cancer cases 10 and is a particularly aggressive cancer type with a high propensity to metastasize. Generally, patients are diagnosed at an advanced stage, which is a negative independent prognostic factor for patients with SCLC. 11
Bone is the third most common metastatic site for a variety of solid tumors. 4 Bone metastases are further categorized as simultaneous bone metastases (SBM) and metachronous bone metastases (MBM). In lung cancer, SBM and MBM might exhibit different clinicopathological features, treatment sensitivities, and prognoses. 12 Because the criteria for synchronous and heterochronic bone metastases of cancer remain unclear, the survival differences in patients with SCLC have not been clearly elucidated, but patients with SBM usually exhibit a significant tumor burden and widespread organ destruction. Additionally, these patients often suffer from emotional and financial stress and poorer prognoses.13,14 However, there are limited studies on the predictive factors and prognosis of BM in patients with SCLC. Therefore, the identification of accurate prognostic factors and development of a convenient and practical prediction model is important to facilitate personalized treatment.
The incidence of BM from SCLC and its associated mortality rate are high. According to a previous report, BM was observed in 27% to 41% of patients with SCLC. 12 In the present study, 22.7% of patients with SCLC exhibited BM, which was similar to previous studies. The mortality rates associated with BM were 90.49%, 90.3%, and 91.04% in 4, 6, and 12 months, respectively, in patients with SCLC. The selection of related risk factors was based on the study conducted by Zhang et al., in which the relationship between different pathological types of lung cancer and BM was explored. In particular, the study included several factors, including population characteristics (sex and race), tumor characteristics (stage and metastatic site), and information regarding treatment methods. 15 Similar to the results reported by Li et al., 16 older age, sex, race, and N stage were found to be related to the occurrence of BM. Compared with patients with a tumor size ≤3 cm, patients with tumors 3 to 7 cm in size exhibited a higher risk of BM. Interestingly, both studies from our laboratory reported no correlation between a larger tumor (tumor size ≥7 cm) and a higher risk of metastasis. Previous studies showed that BM occurs in three stages. Specifically, metastases escape from the primary tumor, enter the circulation, and finally colonize the bone. 17 Li et al. hypothesized that under the regulation of cell adhesion factors, tumors ≥ 7 cm in size would not easily detach from the primary tumor to begin subsequent metastasis. 16 However, based on the discussion of the present study, we predict that the adhesion between tumors might have an impact. Nevertheless, other factors might also affect this particular result. Wang et al. demonstrated that a larger tumor size and multiple lymph node metastases were associated with the development of SCLC. 18 Thus, during the analysis of the effect of tumor size on the occurrence of BM, the status of lymph nodes must be considered. In particular, more precise analyses of several factors might better explain this result. Thus, clinicians should pay more attention to larger tumors and re-examine emission computed tomography or positron emission tomography-computed tomography scans in a timely manner. As shown in the forest plot (Figure 3), the hazard ratio of liver metastasis reached 4.025. This indicated that the presence of liver metastasis would substantially increase the risk of BM. The liver was previously reported to be the most prevalent site of metastasis (61.9%), and liver metastasis was the most common type of hematogenous metastasis in extensive-stage SCLC. 19 Therefore, in patients with liver metastases, closer attention should be paid to the condition of their bones.
Because diverse factors affect the survival time of patients, PSM is used to balance factors that might affect survival analysis. PSM is a commonly used method to eliminate the effect of bias during the statistical analysis of observational data. 20 To the best of our knowledge, the present study is the first to analyze the CSS and OS of patients with BM of SCLC using PSM. This further highlights the importance of determining the influence of BM on the prognosis of SCLC. After PSM, a significant difference was observed in OS and CSS between the BM group and non-BM group. Similarly, the median survival time was different (OS: 6 versus 7 months; CSS: 9 versus 13 months). The results of OS were consistent with a previous study conducted in 2016. 21 This is important as the data included in the present study were collected up to 2016. Furthermore, the survival times of patients with BM before and after PSM were compared. Importantly, no significant difference was detected in the survival time before and after PSM, which further demonstrated that this method did not affect the raw data. For prognostic factors, the results of multivariate Cox proportional hazard regression analysis identified men, larger tumor size, higher N stage, and brain/liver/lung metastases as risk factors for both OS and CSS in patients with BM. In addition, age >65 years was identified as a risk factor for OS but not CSS. One possible explanation is an age-related increase in the probability of death from other causes, whereas age might not be as closely related to CSS. Importantly, Black or other races were identified as risk factors for CSS but not OS. This result is not surprising because African Americans in the United States bear a disproportionate share of the cancer burden, with the highest death rate and shortest survival compared with any other racial or ethnic group for most cancers. 22 Therefore, it is necessary to further the progress aimed at the elimination of these racial disparities. Finally, two novel nomograms were generated for the diagnosis and prognosis of SCLC. Compared with previous studies with a broad scope, the nomogram developed in the present study was found to be more accurate, which might increase the convenience and provide important clues for clinical diagnosis and treatment.
The present study also had certain limitations. First, the study did not include an independent external cohort to validate the model, which is an important focus of a future study. This study presented a retrospective analysis that might have led to bias. Thus, prospective research is required to verify the conclusions of this study. Additionally, some indicators were not provided by the SEER database, such as driver genes, programmed death-ligand 1 expression levels, Eastern Cooperative Oncology Group scores, type of operation, radiotherapy dose, and chemotherapy regimen. Including this information may increase the accuracy of the diagnosis and prognosis model.
For the diagnosis and treatment of patients with cancer, a multidisciplinary team approach is necessary. In fact, oncology, orthopedics, pathology, imaging, and nuclear medicine departments are indispensable. The multidisciplinary diagnosis and treatment of lung cancer have been reported to shorten the time interval between diagnosis and oncology assessment and treatment, further ensuring that patients receive individualized and effective treatment at an earlier stage. This may help further improve the prognosis of patients to a certain extent. 23
Conclusions
The present study showed that multiple factors, such as sex, tumor size, N-stage, surgery, radiotherapy, and liver/brain/lung metastases, might be related to the occurrence of BM in patients with SCLC. Additionally, these factors were identified as independent prognostic factors in patients with BM of SCLC. The OS and CSS of these patients were determined to be poor. The two generated nomograms may be conveniently applied in clinical work to predict the incidence and survival rate for BM from SCLC. In conclusion, the BM of SCLC remains a significant challenge, and additional studies are required to explore and develop novel treatment methods to improve the survival for these patients.
