Abstract
Purpose
We aimed to establish nomograms to predict the survival in patients aged ≥45 years with lung squamous cell carcinoma and brain metastasis.
Methods
We collected patients diagnosed as lung squamous cell carcinoma with brain metastasis aged ≥45 years between 2010 and 2019 from the Surveillance, Epidemiology, and End Results database. Prognostic factors were determined by the univariate and multivariate Cox regression analysis, and then the nomogram was constructed to predict cancer-specific survival and overall survival. Nomograms were evaluated by decision curve analysis, the area under the receiver operating characteristic curve, calibration plot, concordance index, and risk group stratification.
Results
In total, 2437 patients were included, with 1706 and 731 in the cohorts of training and validation, respectively. The age, N stage, T stage, liver metastasis, chemotherapy, bone metastasis, along with radiotherapy were significant in predicting the survival, and adopted for the establishment of nomograms. In the training and validation sets, the concordance index were .713(95%CI:0.699–.728) & .700(95%CI:0.677–.722) in predicting cancer-specific survival and .715(95%CI:0.701–.729) & .712(95%CI:0.690–.735) in predicting overall survival, respectively. Besides, the area under the receiver operating characteristic curve for predicting cancer-specific survival and overall survival in the training set were all >.7 at 1-, 2-, and 3- years. Calibration plots proved the survival predicted by nomograms were consistent with the actual values. decision curve analysis revealed better clinical validity of the nomogram in predicting cancer-specific survival and overall survival at 1-year than TNM staging. Patients were stratified into the high-/low-risk groups according to the optimal cutoff value of 100.21 for cancer-specific survival and 91.98 for overall survival. A web-based probability calculator was constructed finally.
Conclusion
Two nomograms were developed for the prognostic prediction of lung squamous cell carcinoma patients with brain metastasis aged ≥45 years, providing guidance for decision-making in clinical practice.
Keywords
Introduction
As one of the most common malignancy, the non-small cell lung cancer (NSCLC) accounts for nearly 85% of patients with lung cancer (LC), which is the leading cause of deaths in the world. 1 Brain metastasis (BM) occur in about 10%–25% of NSCLC patients at the first visit and about 50% of patients during treatment. 2 The proportion of lung squamous cell carcinoma (LUSC) in all the NSCLCs is about 25%–30%, 3 especially more likely to occur in middle-aged and older patients. In contrast to non-squamous cell carcinoma, incidence of BM has been found to be much lower in LUSC. 4 Nevertheless, there have been few studies on BM of LUSC.
LUSC patients with BM usually present with a poor prognosis regarding of its insidious onset as well as the lack of targeted therapy. Approximately 70% of patients are already at the intermediate/advanced stage at diagnosis, losing opportunities for complete surgical resection. 5 Targetable driver mutations, including epidermal growth factor receptor (EGFR), anaplastic lymphoma kinase (ALK), or ROS proto-oncogene 1 (ROS1) alterations, are exclusively found in the histological samples of adenocarcinoma but not in LUSC. 6 Considering the limitations of treatment methods, the survival of LUSC patients with BM is as short as 3-6 months in median, 7 bringing a heavy burden on their families and the whole society.
One study has reported on the prognostic influences of metastatic NSCLC in the elderly. 8 There is also a study reporting on prognostic influences in LC younger than 45 years of age. 9 In the study, LC patients younger than 45 years of age were defined as having early onset LC. These patients accounted for approximately 5% of all LC patients. 10 Hereditary cancer factors are considered to be the main cause of early onset LC. 11 The World Health Organization (WHO) classifies age as follows: younger people under the age of 44, middle-aged people between the ages of 45 and 59, and older people over the age of 60. LUSC is most often seen in middle-aged and older patients. To date, no studies have reported on the prognostic factors influencing the prognosis of middle-aged and elderly patients with LUSC with BM.
Early diagnosis together with timely therapy is potential for the improvement of the survival of LUSC patients with BM. Since the WHO classification of thoracic tumors in 2021 didn’t incorporate histopathologic criteria for the prognosis of LUSC, 12 TNM staging has become the most frequently used tool for clinical staging and prognostic assessment of LUSC. Unfortunately, the prognostic factors such as the sex, age, histologic grade, or treatments in these patients were neglected in TNM staging. 13 Thus, efficient tools are urgently needed to provide a reference for the diagnosis as well as prognostic assessment of BM in LUSC cases.
Currently, nomograms have been recognized as a unique and reliable method for the prognostic prediction of various tumor in the past decade, including gastric cancer, breast cancer, and testicular cancer,14-16 but few in LUSC patients combined with BM. NSCLC can be classified into LUSC, lung adenocarcinoma (LUAD), adenosquamous carcinoma (ADSQC), and large-cell lung carcinoma (LCLC) according to histopathological classification, and the prognosis of patients with different pathological types varies significantly. Patients with LC with BM have a poor prognosis, and untreated patients have a short survival. Therefore, more attention needs to be paid to this problem and the factors affecting prognosis need to be studied. Prognostic factors influencing LC with BM have been reported in studies. 17 Other scholars have reported the prognostic factors of NSCLC with BM and small cell lung cancer (SCLC) with BM.18,19 In addition, other scholars have reported the prognostic influencing factors of metastatic LUSC and metastatic LUAD, respectively.20,21 To date, no prognostic model for LUSC with BM has been reported in the literature. Thus, this objective of this study was to develop appropriate nomograms to evaluate the prognosis of LUSC patients combined with BM using the Surveillance, Epidemiology and End Results (SEER) database, facilitating the decision-making in clinical practice. We present the following article in accordance with the RECORD reporting checklist.
Patients and Methods
Data Source
Data used in this study were obtained from the Surveillance, Epidemiology, and End Results (SEER) database, covering approximately 30% of the whole U.S. population. Besides, SEER*Stat (version 8.4.0.1) was used (http://seer.cancer.gov/SEERSTAT/) to access the database. The study was carried out complying with the Declaration of Helsinki (as revised in 2013). Local ethics approval or statements were not required because the clinical data used in this study were obtained from the public-access SEER database and thus, the requirement for informed consent was waived.
Patient Selection
Patients with a pathological diagnosis of LUSC (histologic type code: 8070-8074, 8083) between 2010 and 2019 from the SEER database were included. The exclusion criteria included (1) patients with non-brain metastasis (n = 94 561); (2) patients less than 45 years old (n = 56); (3) patients with non-pathological diagnosis (n = 623); (4) patients with non-first primary cancer (n = 834); (5) patients with unknown Race (n = 8); (6) patients with unknown Marital (n = 147); (7) patients without left or right laterality (n = 129); (8) patients with T0\Tx cancer staging (n = 298); (9) patients with Nx cancer staging (n = 90); (10) patients with unknown bone metastasis (n = 36); (11) patients with unknown liver metastasis (n = 27); (12) patients with unknown lung metastasis (n = 23); (13) patients with unknown surgery (n = 2); (14) patients with the survival <1 month (n = 311). The flowchart for inclusion and exclusion is presented in Figure 1. Flowchart for inclusion and exclusion of LUSC patients with BM aged ≥45 years.
Data Collection
The demographic information of the patients (such as the age), the clinicopathological information of the tumor (such as the primary location) as well as results of follow-up (such as the survival months) were collected from the SEER database.
The included patients were divided into 3 groups (45–60 years old, 61–70 years old, and ≥71 years old) using the optimal cutoff value of the age calculated by X-tile (version 3.6.1, Yale University School of Medicine, US). Figure S1 shows estimation of the cutoff value for the age determined by X-tile software. Besides, to avoid the error caused by incongruous stage system, we completed the transformation of TNM staging and evaluated the patient data according to the 7th edition of TNM.
Endpoints included cancer-specific survival (CSS) and overall survival (OS). CSS referred to the duration between the diagnosis and a specific cancer‐related death; while OS referred to the duration between the diagnosis and the all-cause death/the final follow‐up.
Statistical Analysis
All statistical analyses were carried out via SPSS (26.0, IBM Corporation, NY, USA) together with R (4.1.0, R Foundation for statistical computing, Vienna, Austria). Patients were allocated into either the training (70%) or validation (30%) cohort. Categorical data were compared between cohorts using the chi-square test. Besides, the prognostic factors for the training set were determined through univariate & multivariate Cox proportional regression, and then a nomogram was constructed for prognostic prediction. The concordance index (C-index) together with the area under the receiver operating characteristic curve (AUROC) was adopted to quantify the discriminative power of nomograms. In addition, calibration plots (1000 bootstrap replicates) were used to calibrate the nomograms. Decision curve analysis (DCA) was conducted for the evaluation of their clinical validity. Total points of individuals were then computed using the nomograms. Finally, all patients were stratified into groups with high/low risks using the optimal cutoff value of receiver operating characteristic (ROC) curves, and their survival were compared via Log-rank test as well as Kaplan–Meier (K-M) curves. A bilateral P value less than .05 was indicative of statistical significance.
Results
Analysis of Demographic as Well as Clinicopathological Characteristics
In total, there were 2437 patients with LUSC and BM between 2010 and 2019 included, and then randomly allocated to the training cohort (N = 1706) along with the validation cohort (N = 731). There were 1562 (64.10%) males and 875 (35.90%) females. Besides, their age ranged within 45–96 years, including 602 (24.70%), 973 (39.93%), and 862 (35.37%) in the groups of 45-60 years, 61-70 years, and ≥71 years, respectively, suggesting the predominance of LUSC with BM in elderly patients. The majority of patients (80.10%) were white. More than half of patients (52.03%) were married, the 27.99% were separated, divorced or widowed (SDW), and 19.98% were unmarried.
There were 57.90% of patients with LUSC in the right lung and 42.10% in the left lung. The lesions occurred mainly in the upper lobe, accounting for 54.37% of all. Patients at Grade III/IV, Grade I/II, and unknown Grade accounted for 13.46%, 29.63%, and 56.91%, respectively. According to the 7th edition of TNM staging of American Joint Committee on Cancer (AJCC), the primary tumors were mainly in T4, T3, T2, and T1 stages in 865 (35.49%), 628 (25.77%), 760 (31.19%), and 184 (7.55%) patients, respectively. In addition, the patients tended to have later T stages (25.77% in T3 and 35.49% in T4) and later N stages (50.47% in N2 and 17.32% in N3). Bone metastasis accounted for 27.37%, lung metastasis accounted for 21.99% and liver metastasis accounted for 17.85%. In terms of the treatments, radiotherapy, and chemotherapy were the predominant ones, accounting for 77.10% and 52.56%, respectively, only 4.51% of patients underwent surgery.
Clinicopathological Characteristics of Patients Aged 45 Years or Older With LUSC With BM.
SDW: separated, divorced, or widowed.
Univariate as well as Multivariate Cox Regression for Predictive factors of Prognosis in the Training Cohort
Univariate and Multivariate Analyses of Risk Factors for CSS in Training Cohort.
REF: reference; HR, hazard ratio; SDW: separated, divorced, or widowed.
Univariate and Multivariate Analyses of Risk Factors for OS in Training Cohort.
REF: reference; HR, hazard ratio; SDW: separated, divorced or widowed.
Construction and Calibration of Nomograms
We constructed nomograms for the prediction of the CSS (Figure 2(A)) as well as OS (Figure 2(B)) at 1-/2-/3-year among LUSC patients with BM using the 7 significant prognostic factors found in multivariate Cox regression models. The nomograms for predicting 1-, 2-, 3-year CSS and OS in the training set of LUSC patients with BM aged ≥45 years. (A) The nomogram in predicting CSS. (B) The nomogram in predicting OS.
We further calibrated the nomograms using calibration plots, and the predicted probabilities for CSS (Figure 3(A) and (B)) and OS (Figure 3(C) and (D)) by nomograms were found to be consistent with the actual observed values in the training and validation sets, suggesting a good accuracy of the nomograms. Calibration plots of the nomograms in predicting 1-, 2-, 3-year CSS and OS in LUSC patients with BM aged ≥45 years. (A) Calibration plots of the nomograms in predicting 1-, 2-, 3-year CSS in the training set. (B) Calibration plots of the nomograms in predicting 1-, 2-, and 3-year CSS in the validation set. (C) Calibration plots of the nomograms in predicting 1-, 2-, 3-year OS in the training set. (D) Calibration plots of the nomograms in predicting 1-, 2-, 3-year OS in the validation set. The horizontal axis is the predicted value in the nomogram, and the vertical axis is the observed value.
Validation of Nomograms Using C-Index, Area Under The Receiver Operating Characteristic Curve and Decision Curve Analysis
Internal validation was used to test discriminative power of nomograms. For the training and validation sets, the C-index was .713(.699–.728) and .700(.677–.722) in predicting CSS, and .715(.701–.729) and .712(.690–.735) in predicting OS, respectively.
Besides, the AUROC at one-, two-, and three-years in predicting CSS was found to be .759, .717, and .737 for the training set as well as .725, .651, and .647 for the validation set, respectively (Figure 4(A) and (B)); whereas the AUROC for the prediction of the OS at one-, two-, and three-years was .761, .722, and .741 for training set as well as .737, .669, and .687 for the validation set, respectively (Figure 4(C) and (D)). The AUROC for predicting 1-, 2-, and 3-year CSS and OS in LUSC patients with BM aged ≥45 years. (A) The AUROC at 1-, 2-, and 3-year for the prediction of CSS in the training set was .759, .717, and .737, respectively. (B) The AUROC at 1-, 2-, and 3-year for the prediction of CSS in the validation set was .725, .651, and .647, respectively. (C) The AUROC at 1-, 2-, and 3-year for the prediction of OS in the training set was .761, .722, and .741, respectively. (D) The AUROC at 1-, 2-, and 3-year for the prediction of OS in the validation set was .737, .669, and .687, respectively.
DCAs for the nomograms as well as TNM staging systems are shown in Figure 5. Results suggested that the nomograms for the prediction of CSS was more beneficial at one-year in either the training or validation set in contrast to TNM staging; In comparison, the nomogram for the prediction of CSS at 2-, 3-year in either the training or validation set showed no apparent advantages over TNM staging (Figure 5(A) and (B)). Meanwhile, the nomogram for the prediction of OS was more beneficial at one-year in either the training or validation set in contrast to TNM staging; in contrast, the nomogram for the prediction of OS at two- and three-year in either the training or validation set showed no apparent advantages over TNM staging (Figure 5(C) and (D)). DCA of the nomograms for predicting CSS and OS in LUSC patients with BM aged ≥45 years. (A) In the training set, the nomogram in predicting CSS at 1 year was more beneficial compared to TNM staging, while the nomogram for CSS at 2, 3-year had no apparent advantages over TNM staging. (B) In the validation set, the nomogram in predicting CSS at the 1-year was more beneficial compared to TNM staging, while the nomogram for CSS at 2, 3-year had no apparent advantages over TNM staging system. (C) In the training set, the nomogram in predicting OS at 1-year was more beneficial compared to TNM staging, while the nomogram in predicting OS at 2, 3-year had no apparent advantages over TNM staging. (D) In the validation set, the nomogram in predicting OS at 1-year was more beneficial compared to TNM staging, while the nomogram for OS at 2, 3-year had no apparent advantages over TNM staging system.
Based on the above findings, it could be inferred that the nomograms had a favorable discriminative capability.
Risk Stratification of Patients Based on Scores of Nomograms
Total score of each patient was computed using the nomograms. Besides, the optimal cutoff value was determined by the receiver operating characteristic (ROC) curve. Based on the optimal cutoff values (100.21 for the prediction of CSS; 91.98 for the prediction of OS), patients were then further stratified into the groups with high or low risks. In addition, K-M curve revealed that compared to the group with low risks, the rate of CSS or OS were lower in the group with high risks in both cohorts (Figure 6) with statistical significance. The Kaplan-Meier curves of LUSC patients with BM aged ≥45 years in the low-risk and high-risk groups. The K-M curve showed that the CSS rate of the patients in the high-risk group was significantly lower than that in the low-risk group both in the training set (A) and validation set (B). The K-M curve showed that the OS rate of the patients in the high-risk group was significantly lower than that in the low-risk group both in the training set (C) and validation set (D).
Construction of Web-Based Probability Calculator
To make our nomograms more convenient for the application in clinical practice, we constructed a web-based calculator for the probability calculation of OS and CSS (https://fengyang.shinyapps.io/DynNomapp_OS/ and https://fengyang.shinyapps.io/DynNomapp_CSS/) for LUSC patients with BM aged ≥45 years. The survival probability of patients could be directly displayed once the variables of individual patients were inputed into the calculator.
Discussion
There have been several clinical models established for LUSC patients.20,22 In the present study, two nomograms were developed and specifically predicted the OS and CSS in LUSC patients with BM aged ≥45 years based on SEER database for the initial time. The age, N stage, T stage, liver metastasis, chemotherapy, bone metastasis, as well as radiotherapy were identified as predictive factors for the prognosis and included in the nomograms. Subsequently, patients were further stratified into groups with high-/low-risks according to the scores of nomograms. The two nomograms may have a potential for the convenient prediction of OS & CSS of individual patients, providing significant implications for decision-making in clinical practice.
The role of age in the prognosis of LUSC patients has been controversial. Hu et al 23 found the age not to be a predictor of the survival among NSCLC patients combined with pericardial effusion using the SEER database. However, Gu et al has demonstrated that increasing age was correlated with a poor lung cancer-specific survival (LCSS) for LUSC patients, 24 consistent with the findings in our study. Therefore, more studies would be needed to further evaluate the impact of age on the survival of LUSC patients. Previous literature has reported that marital status has been shown to significantly affect the prognosis of patients with LUAD, 25 as divorced patients are less likely to enjoy a relatively harmonious family environment. In our study, marital status was included as one of the basic demographic information. Our study found through statistical analysis that marriage was not an independent factor affecting prognosis, and therefore, marriage was not included in the final prediction model. Further future studies are needed to assess the role of marital status on prognosis.
TNM staging is an internationally recognized prognostic factor for LC. The 8th edition of AJCC staging system emphasized the importance of primary tumor size in the treatment and prognosis of LC. 26 The diameter of the primary tumor has been proved to be related with the number of clonal cells, which is indicative of differentiation of cells, the malignancy of tumors, as well as the possibility of distant metastases. 27 Several studies have revealed a significant association of tumor size with the prognosis of NSCLC patients after surgery, radiotherapy, and chemotherapy,28-30 which was similar to the findings in our study.
Researches have proved bone/liver metastases to be associated with the prognosis of LC patients, 31 which is consistent with the results of our study. Bone metastasis has been found in 30%–40% of NSCLC patients, among which approximately 6% developed within one year after the diagnosis of LC, characterized by the pain and pathological fractures; Patients with bone metastases showed 4.4 months shorter OS compared to those without bone metastases. 32 The 1-year survival rate in patients without bone metastasis was about 37%, which was more than 3 times of that in those with bone metastasis. 33 Compared to brain and bone metastases, liver metastases presented the worst prognosis among single-organ metastases among the patients with LC, with the OS of only 4 months in median. 31 An analysis based on the US oncology database between 2010 and 2015 found that the 1-year survival rate for all patients with liver metastases was 15.1%, which was obviously lower than 24.0% for patients without liver metastases. 34
With the advent of genomics, targeted therapies have brought significant benefits for LC patients. However, chemotherapy has remained the preferred treatment for patients with advanced LUSC in clinical practice due to the low rate of genetic mutations in them. Some researchers found that chemotherapy significantly prolonged their survival, including elderly patients. 22 Radiotherapy is an essential treatment for patients with stage IV LUSC, mainly including local and systemic radiotherapy. An analysis based on SEER database between 2010 and 2015 found that radiotherapy was significantly associated with OS rates in LUSC. 35 Consistent with the above findings, our study also demonstrated the significant roles of chemotherapy and radiotherapy in LUSC patients.
In this study, there were still several limitations. First, selection bias would be inevitable since it was retrospective study and composed of only patients having complete data. Second, SEER database included no information of some specific prognostic factors, such as smoking history, family cancer history, genetic mutations, as well as the assessment of physical status. Third, the information on targeted therapies and immunotherapy were unavailable in the SEER database. Due to the lack of relevant information on tyrosine kinase inhibitors (TKIs) and immune checkpoint inhibitors treatment in the SEER database, the two factors were not included in the study, which is one of the limitations of our study. In the future, it is necessary to include TKIs and immune checkpoint inhibitors in clinical studies. Besides, there was no information on systemic therapies, especially the type of operations, dose for radiation, together with the choice of chemotherapeutic agents. Fourth, the SEER database only contained the data of the single-organ metastases at diagnosis but not during follow-up. Finally, only internal validation was performed in this study. More studies would be necessary for external validation in different populations to evaluate the predictive value of our nomograms in the future.
Conclusion
Two nomograms were successfully established for the prediction of the CSS/OS among LUSC patients combined with BM aged ≥45 years based on the significant prognostic factors, including the age, N stage, T stage, liver metastasis, chemotherapy, bone metastasis, as well as radiotherapy. The two nomograms presented with a good accuracy and a discriminative power, potentially making the clinical decision-making more convenient in the real world.
Supplemental Material
Supplemental Material - Development and Validation of Prognostic Nomograms for Lung Squamous Cell Carcinoma With Brain Metastasis in Patients Aged 45 Years or Older: A Population-Based Study
Supplemental Material for Development and Validation of Prognostic Nomograms for Lung Squamous Cell Carcinoma With Brain Metastasis in Patients Aged 45 Years or Older: A Population-Based Study by Feng Yang, Lianjun Gao, Qimin Wang, and Wei Gao in Cancer Control
Footnotes
Acknowledgments
We thank the SEER database for providing clinical information, and we also thank CureEdit for English language editing.
Declaration of Conflicting Interests
All authors have completed the ICMJE uniform disclosure form. The authors have no conflicts of interest to declare.
Funding
The author(s) received no financial support for the research, authorship, and/or publication of this article.
Ethical Approval
The authors are accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved. The SEER database is an open-access cancer database that only contains de-identified patient data. Therefore, this study was exempted from the approval by the Institutional Review Board of China Rehabilitation Research Center.
Supplemental Material
Supplemental material for this article is available online.
Appendix
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
