Identifying optimal candidates for primary tumor resection among metastatic non-small cell lung cancer patients: a population-based predictive model.

BACKGROUND: A survival benefit was observed in metastatic non-small cell lung cancer (NSCLC) patients that underwent surgical resection of the primary tumor. We developed a model testing the hypothesis that only certain stage IV patients would benefit from surgery and the potential benefit would vary based on primary tumor characteristics.
METHODS: Patients with stage IV NSCLC were identified in the Surveillance, Epidemiology and End Results (SEER) database and then divided into surgery and non-surgery groups. A 1:1 Propensity score matching (PSM) was performed to balance characters. We assumed that patients received primary tumor surgery that lived longer than median cancer specific survival (CSS) time of those who didn't underwent surgery could benefit from the operation. Multivariable Cox model was used to explore the independent factors of CSS in two groups (beneficial and non-beneficial group). Logistic regression was used to build a nomogram based on the significant predictive factors.
RESULTS: A total of 30,342 patients with stage IV NSCLC were identified; 8.03% (2,436) received primary tumor surgery. After PSM, surgical intervention was independently correlated with longer median CSS time (19 vs. 9 months, P<0.001). Among the surgery cohort, 1,374 (56.40%) patients lived longer than 9 months (beneficial group). Differentiated characters (beneficial and non-beneficial group) included age, gender, TNM stage, histologic type, tumor position and differentiation grade, which were integrated as predictors to build a nomogram.
CONCLUSIONS: A practical predictive model was created and might be used to identify the optimal candidates for surgical resection of the primary tumor among stage IV NSCLC patients. 2021 Translational Lung Cancer Research. All rights reserved.

Introduction

Lung cancer is one of the most common cancers worldwide as well as the leading cause of cancer death (1). Non-small cell lung cancer (NSCLC) comprises over 80% of lung cancer cases (2). Approximately half of all NSCLC patients present with metastatic disease at diagnosis, therapeutic scheme are rather limited and life expectancy is greatly compromised. In these patients, chemotherapy, targeted therapy and immunotherapy represent the primary choice of treatment (3).

Local resection to primary tumor site is not usually taken into account as conventional treatment in stage IV NSCLC patients (4). The primary focus and goals in such patients is disease control and palliation but not radical cure (5). Many studies have demonstrated the primary tumor resection was associated with improved survival outcomes. However, surgical intervention whether to primary or metastatic tumor site for advanced-stage NSCLC has been proved to be associated with improvements in survival (6-8) in advanced NSCLC patients (9,10), especially in M1a patients without N3 or pericardial effusion disease (10). From these studies, it can be concluded that not all stage IV patients prolong their lives from primary tumor surgery, the specific group of patients who do benefit from the surgery are still not fully understood.

We aimed to develop a predictive model using a prospective national database to test the hypothesis that only certain patients would benefit from surgery and that the potential benefit would vary depend on primary tumor characteristics. We present the following article in accordance with the Transparent Reporting of a multivariable prediction model for Individual Prognosis Or Diagnosis (TRIPOD) reporting checklist (11) (available at http://dx.doi.org/10.21037/tlcr-20-709).

Methods

Patient selection in SEER database

The Surveillance, Epidemiology and End Results (SEER) database (1973–2014, November 2016 submission) is a national population-based reporting system that collects tumor-related data, including the incidence, treatment, mortality, and other demographics, covering approximately 28% of the US population (12). The SEER data contain no identifiers and are publicly available for studies of cancer-based epidemiology and survival analysis. We had got the permission to access these research data (SEER-Stat username: lianghr). The study was conducted in accordance with the Declaration of Helsinki (as revised in 2013).

Cases of lung cancer (C34.0-34.9) diagnosed from 2004 to 2014 were extracted from the SEER database (SEER-Stat 8.3.5) according to the site code classifications. We selected this range because American Joint Committee on Cancer (AJCC) TMN stage and Collaborative Stage (CS) information was available since 2004 and patients diagnosed after 2014 were excluded to ensure an adequate follow-up time. We reclassified the TNM stage according to AJCC 8th edition. Patients diagnosed as stage IV and histologically confirmed as NSCLC were enrolled. Patients were excluded if the surgery to primary site record was unknown. Other exclusion criteria were as follows: less than 18 years old, unknown TNM stage, unknown survival months, unknown treatment modality, not the first tumor and not only one tumor.

Statistical analysis

The study sample was divided according to treatment, surgery versus non-surgery treatment of the primary tumor. Propensity-score matching (PSM) was calculated to account for confounding variables and facilitate matching patients in the two treatment groups (R software version 2.15.1, https://www.rproject.org/). Variables that could influence the outcomes of treatment were used to generate a propensity score by logistic regression, including age, gender, TNM stage, differentiation grade, primary site, laterality, histological type, radiation therapy, chemotherapy and surgery to metastasis site. Patients in the two groups (surgery and non-surgery) were 1:1 matched using the nearest propensity score on the logit scale. The caliper was set at 0.01. After PSM, differences in continuous (age) and categorical clinical characteristics (others except for age) were tested for significance by student t test and chi-square tests, respectively.

Overall survival (OS), cancer specific survival (CSS) and survival month were extracted from SEER database. OS was the time from diagnosis to death from any cause; living patients were excluded at the time of last recording. CSS was calculated from the date of diagnosis to the date of cancer specific death. Deaths attributed to the lung cancer were treated as events and deaths from other causes were treated as censored observations. OS and CSS were estimated by the Kaplan-Meier (K-M) method and compared with the log-rank test. Multivariate Cox proportional hazard regression was used to determine independent prognostic factors. Hazard ratios (HRs) were calculated with 95% confidence intervals (CIs). Statistical analysis was performed with SPSS 24.0 (IBM Corp., Armonk, NY, USA), all statistical tests were two-sided, and P<0.05 was considered statistically significant.

Construction and validation of nomogram

We assumed that patients received primary tumor resection that lived longer than median CSS time of those who didn’t underwent surgery could benefit from the operation. According to this assumption, eligible patients who had undergone surgery to primary site were randomly divided into training and validation cohorts by 7:3. Both cohorts were then separately divided into two groups, beneficial group and non-beneficial group, according to the median CSS time of the non-surgery group (9 months, result from the cohort after PSM).

Factors that independently affect the CSS in multivariate Cox analysis and can be accessed before surgery were collected into the training set, including age, gender, TNM stage, histologic type, tumor position (over lapping lesion, bronchus tumor, lobe tumor) and differentiation grade. Based on the multivariable logistic analysis, we established the nomogram to predict the stage IV NSCLC candidates who can indeed benefit from primary site surgery in the training cohort. A nomogram was formulated by R software (version 2.15.1) with the rms package.

The areas under the ROC curve (AUC) and a calibration plots were derived based on the validation cohort by the regression analysis to assess the discriminative ability and accuracy of the nomograms, respectively.

According to the result of our predictive model, we divided all stage IV NSCLC patients in SEER database after PSM into three groups, one is beneficial & surgery group whose assuming probability of benefit is over 50%, and another is non-beneficial & surgery group with less than 50% assuming probability of benefit, the other is non-surgery group. Kaplan-Meier analyses and log-rank test were conducted to test whether this model would distinguish patients that could benefit from primary tumor resection.

Results

Patient characteristics before and after PSM

A total of 565,424 lung cancer patients were identified in SEER database from 2004 to 2014 and 30,342 stage IV NSCLC patients met the criteria (). Of these eligible patients, 2,436 (8.03%) underwent surgical treatment to primary tumor site. Obvious differences in age, gender, race, TNM stage, location, lesion laterality, histology, differentiation, radiation and distant surgery were noted between the two treatment groups before PSM (). Specifically, the surgery intervention was associated with small amount of older patient and higher TNM stage patients. This indicated that the baseline characteristics of the two groups (surgery and non-surgery) were not balanced. Variables that could influence the outcomes of treatment were included in 1: 1 PSM, including age, gender, TNM stage, differentiation grade, primary site, laterality, histological type, radiation therapy, chemotherapy and surgery to metastasis site.

Figure 1

Flow chart of predictive model construction. SEER, Surveillance, Epidemiology and End Results; PSM, propensity score matching.

Table 1

Demographic information for patients with stage IV NSCLC before PSM

	All patients (n=30,342)	Surgery to primary site (n=2,436), n (%)	Non-surgery to primary site (n=27,906), n (%)	P
Age				<0.001
Mean	66.67	63.97	66.90
SD	11.48	10.88	11.15
Gender				<0.001
Male	16,976	1,261 (51.77)	15,715 (56.31)
Female	13,366	1,175 (48.23)	12,191 (43.69)
Race				<0.001
White	23,814	1,995 (81.90)	21,819 (78.19)
Black	4,077	268 (11.00)	3,809 (13.65)
Other	2,451	173 (7.10)	2,278 (8.16)
T stage				<0.001
T1	6,309	1,006 (41.30)	5,303 (19.00)
T2	7,908	610 (25.04)	7,298 (26.15)
T3	6,189	380 (15.60)	5,809 (20.82)
T4	9,936	440 (18.06)	9,496 (34.03)
N stage				<0.001
N0	7,622	1,168 (47.95)	6,454 (23.13)
N1	2,788	383 (15.72)	2,405 (8.62)
N2	14,616	765 (31.40)	13,851 (49.63)
N3	5,316	120 (4.93)	5,196 (18.62)
M stage				<0.001
M1a	6,662	774 (31.77)	5,888 (21.10)
M1b	23,680	1,662 (68.23)	22,018 (78.90)
Tumor position				<0.001
Right bronchus	897	44 (1.81)	853 (3.06)
Right upper lobe	4,951	452 (18.56)	4,499 (16.12)
Right middle lobe	10,134	749 (30.75)	9,385 (33.63)
Right lower lobe	1,375	133 (5.46)	1,242 (4.45)
Left bronchus	623	24 (0.99)	599 (2.15)
Left upper lobe	8,113	650 (26.68)	7,463 (26.74)
Left lower lobe	3,831	321 (13.18)	3,510 (12.58)
Over lapping lesion	418	63 (2.59)	355 (1.27)
Histology				<0.001
Large cell	1,333	115 (4.72)	1,218 (4.36)
Adenocarcinoma	18,573	1,649 (67.69)	16,924 (60.65)
Squamous cell	9,130	560 (22.99)	8,570 (30.71)
Adeno-squamous	517	80 (3.28)	437 (1.57)
Neuroendocrine	789	32 (1.31)	757 (2.71)
Grade				<0.001
I	1,789	195 (8.00)	1,594 (5.71)
II	9,184	904 (37.11)	8,280 (29.67)
III	18,164	1,240 (50.90)	16,924 (60.65)
IV	1,205	97 (3.98)	1,108 (3.97)
Radiation				<0.001
Yes	15,450	1,087 (44.62)	14,363 (51.47)
None/unknown	14,892	1,349 (55.38)	13,543 (48.53)
Chemotherapy				0.910
Yes	16,832	1,354 (55.58)	15,478 (55.46)
None/unknown	13,510	1,082 (44.42)	12,428 (44.54)
Surgery to other sites				<0.001
Yes	1,574	374 (15.35)	1,200 (4.30)
None/unknown	28,768	2,062 (84.65)	26,706 (95.70)

NSCLC, non-small cell lung cancer; PSM, propensity score matching; SD, standard deviation.

After the 1:1 PSM, 4,872 stage IV NSCLC patients treated with or without primary site surgery were enrolled in the following survival analysis. Baseline characteristics, including age, gender, race, TNM stage, primary site, laterality, histology, differentiation, chemotherapy, radiation therapy and surgery to distant tumor, were all balanced (P>0.05), the detail is presented in .

Table 2

Demographic information for patients with stage IV NSCLC after PSM

	Surgery to primary site (n=2,436), n (%)	Non-surgery to primary site (n=2,436), n (%)	P
Age			0.877
Mean	63.97	63.56
SD	10.88	11.02
Gender			0.566
Male	1,261 (51.77)	1,247 (51.19)
Female	1,175 (48.23)	1,189 (48.81)
Race			0.994
White	1,995 (81.90)	1,994 (81.80)
Black	268 (11.00)	273 (11.20)
Other	173 (7.10)	169 (7.00)
T stage			0.673
T1	1,006 (41.30)	946 (38.80)
T2	610 (25.04)	624 (25.60)
T3	380 (15.60)	415 (17.00)
T4	440 (18.06)	451 (18.60)
N stage			0.142
N0	1,168 (47.95)	1,158 (47.60)
N1	383 (15.72)	345 (14.20)
N2	765 (31.40)	785 (32.20)
N3	120 (4.93)	148 (6.00)
M stage			0.758
M1a	774 (31.77)	771 (31.6)
M1b	1,662 (68.23)	1,665 (68.4)
Tumor position			0.929
Right bronchus	44 (1.81)	43 (1.77)
Right upper lobe	452 (18.56)	460 (18.88)
Right middle lobe	749 (30.75)	744 (30.54)
Right lower lobe	133 (5.46)	129 (5.30)
Left bronchus	24 (0.99)	26 (1.10)
Left upper lobe	650 (26.68)	642 (26.35)
Left lower lobe	321 (13.18)	330 (13.55)
Over lapping lesion	63 (2.59)	56 (2.20)
Histology			0.878
Large cell	115 (4.72)	105 (4.40)
Adenocarcinoma	1,649 (67.69)	1,663 (68.20)
Squamous cell	560 (22.99)	567 (23.20)
Adeno-squamous	80 (3.28)	66 (2.80)
Neuroendocrine	32 (1.31)	35 (1.40)
Grade			0.492
I	195 (8.00)	173 (7.20)
II	904 (37.11)	885 (36.40)
III	1,240 (50.90)	1,281 (52.60)
IV	97 (3.98)	97 (4.00)
Radiation			0.686
Yes	1,087 (44.62)	1,346 (55.20)
None/unknown	1,349 (55.38)	1,090 (44.80)
Chemotherapy			0.437
Yes	1,354 (55.58)	1,080 (44.40)
None/unknown	1,082 (44.42)	1,356 (56.60)
Surgery to other sites			0.469
Yes	374 (15.35)	355 (14.60)
None/unknown	2,062 (84.65)	2,436 (85.40)

NSCLC, non-small cell lung cancer; PSM, propensity score matching; SD, Standard deviation.

Impact of primary tumor resection on survival outcomes in stage IV NSCLC patients

In the K-M analyses and log-rank test in matched population, as seen in and , patients who underwent primary tumor resection enjoyed longer OS and CSS. The median CSS time was 19 months (95% CI, 17.96–21.04) for patients who underwent primary tumor resection while only 9 months (95% CI, 8.54–10.46) for patients who didn’t receive surgery after PSM (HR =0.47, 95% CI, 0.44–0.51, P<0.001). The 1-year CSS rate was 31.1% and 50.9% for non-surgery and surgery groups, respectively. The 3-year CSS rate was 7.1% and 20.3% for non-surgery and surgery patients after PSM (Table S1).

Table 3

Multivariate Cox analysis for OS and CSS among PSM population

	CSS			OS
	Adjust HR (95% HR)	P value		Adjust HR (95% HR)	P value
Age
≤65	Reference			Reference
65-75	1.21 (1.12–1.30)	<0.001		1.24 (1.15–1.33)	<0.001
≥75	1.42 (1.28–1.57)	<0.001		1.48 (1.35–1.64)	<0.001
Gender
Male	Reference			Reference
Female	0.82 (0.77–0.88)	<0.001		0.81 (0.76–0.87)	<0.001
Race
White	Reference			Reference
Black	1.01 (0.91–1.13)	0.855		1.06 (0.96–1.11)	0.281
Other	0.66 (0.57–0.76)	<0.001		0.68 (0.60–0.78)	<0.001
T stage
T1	Reference			Reference
T2	1.06 (0.97–1.15)	0.216		1.06 (0.98–1.15)	0.156
T3	1.32 (1.20–1.45)	<0.001		1.30 (1.18–1.42)	<0.001
T4	1.55 (1.41–1.71)	<0.001		1.53 (1.40–1.68)	<0.001
N stage
N0	Reference			Reference
N1	1.37 (1.24–1.51)	<0.001		1.34 (1.22–1.48)	<0.001
N2	1.42 (1.32–1.54)	<0.001		1.40 (1.30–1.51)	<0.001
N3	1.56 (1.35–1.81)	<0.001		1.56 (1.35–1.79)	<0.001
M stage
M1a	Reference			Reference
M1b	1.56 (1.44–1.68)	<0.001		1.48 (1.37–1.59)	<0.001
Tumor position
Bronchus	Reference			Reference
Lobe	1.01 (0.83–1.23)	0.945		1.03 (0.85–1.24)	0.763
Over lapping lesion	1.31 (0.99–1.74)	0.061		1.33 (1.02–1.75)	0.038
Differentiation
Well	Reference			Reference
Moderately	1.26 (1.09–1.45)	0.002		1.28 (1.12–1.47)	<0.001
Poorly	1.58 (1.38–1.82)	<0.001		1.60 (1.40–1.83)	<0.001
Undifferentiated	1.96 (1.59–2.40)	<0.001		1.97 (1.61–2.40)	<0.001
Chemotherapy
None/unknown	Reference			Reference
Yes	0.57 (0.53–0.61)	<0.001		0.55 (0.51–0.58)	<0.001
Surgery to other sites
None/unknown	Reference			Reference
Yes	0.77 (0.69–0.84)	<0.001		0.76 (0.70–0.83)	<0.001
Surgery to primary site
None/unknown	Reference			Reference
Yes	0.49 (0.46–0.52)	<0.001		0.51 (0.47–0.54)	<0.001

NSCLC, non-small cell lung cancer; PSM, propensity score matching; SD, standard deviation; OS, overall survival; CSS, cancer specific survival; HR, hazard ratio.

Figure 2

Kaplan-Meier plot of stage IV NSCLC patients according to treatment. NSCLC, non-small cell lung cancer; SEER, Surveillance, Epidemiology and End Results; PSM, propensity score matching; CSS, cancer specific survival; OS, overall survival; HR, hazard ratio.

Primary tumor resection as an independent prognostic factor for survival in stage IV NSCLC patients

In multivariate Cox analysis (), primary tumor resection persisted to be independently associated with better OS (HR =0.51, 95% CI, 0.47–0.55, P<0.001) and CSS (HR =0.49, 95% CI, 0.46–0.52, P<0.001) in matched population. Besides, age, gender, race, TNM stage, differentiation grade, primary site, chemotherapy and surgery to distant site are all independent factor for survival of stage IV NSCLC patients. Radiation and histology had no significant effect on overall survival and cancer specific survival.

Nomogram to identify optimal candidates for primary tumor resection

We assumed that patients received primary tumor surgery that lived longer than median CSS time (9 month) of these who didn’t underwent surgery could benefit from the operation. Among the surgery cohort, 1,374 (56.40%) patients lived longer than 9 months (beneficial group). According to this assumption, surgery receivers who survived longer than 9 months were classified into beneficial group and those who had shorter median CSS time (less than 9 months) were classified into non-beneficial group.

Based on beneficial and non-beneficial group, factors that independently affect the CSS in multivariate Cox analysis and can be accessed before surgery were collected into the training set, including six index: age, gender, TNM stage, histologic type, tumor position and differentiation grade. Based on the multivariate logistic regression, we established a nomogram to predict the stage IV NSCLC candidates who can indeed benefit from primary site surgery in the training cohort. A nomogram was then print (). This model revealed that tumor position, age and differentiation contributed most to prognosis, followed by the T stage, N stage, histology, gender and M stage. By adding the scores of each selected variable, the likelihood of survival of the individual patient can be easily calculated.

Figure 3

A nomogram to predict optimal candidates for primary tumor resection. ASC, adeno-squamous carcinoma; LCC, Large cell carcinoma; SCC, Squamous cell carcinoma; AC, adenocarcinoma; NC, Neuroendocrine.

Validation of the predictive model

Internal validation via the training set indicated that the AUC index for the nomogram to predict CSS were 0.62 (95% CI, 0.603–0.645), which was in agreement with actual CSS. External validation using the validation set, the AUC for the nomogram to predict CSS was 0.60 (95% CI, 0.586–0.621) (Figure S1). The internal and external calibration plots of OS nomogram revealed a good correlation between prediction by the nomogram and actual observation (Figure S2).

We then validated the distinguish ability of the model in validation set. According to Kaplan-Meier analyses and log-rank test (), beneficial & surgery group is significantly live longer than non-beneficial & surgery group (HR =0.57, 95% CI, 0.39–0.79, P<0.001) or non-surgery group (HR =0.54, 95% CI, 0.40–0.74, P<0.001), but there no difference was observed between non-beneficial & surgery and non-surgery group (HR =0.97, 95% CI, 0.87–1.07, P=0.694).

Figure 4

Kaplan-Meier plot to differentiate beneficial group according to our model in SEER database after PSM. SEER, Surveillance, Epidemiology and End Results; PSM, propensity score matching.

Clinical use

To use the nomogram, first draw a vertical line to the top points row to assign points for each variable; then, add the points from each variable together and drop a vertical line from the total points row to obtain the probability of benefit from primary tumor resection (). For example, for a 68-year-old male patient with poorly differentiated bronchus adenocarcinoma which is classified as T3N2M1a stage, the total score is 220 points and the model predicts that the patients would benefit from primary tumor resection with the probability of 30%.

Because the preoperative small specimen could not represent the whole tumor and the neoadjuvant therapy might alter tumor characteristics, not all cases represented accurate differentiation information. As a consequence, we constructed another model without differentiation information, the nomogram is presented in Figure S3.

Discussion

These results confirmed our hypothesis and showed several important findings. Firstly, our study is in line with the previous reports (5,7) that suggested a survival benefit of primary tumor resection for the stage IV in NSCLC patients. Secondly, different from the other study (13) which showed local consolidative therapy for patients with metastases NSCLC improved progression-free survival, we aimed to identify optimal candidates for surgery of the primary tumor and also try to comprehensively describe their clinical characters and power. Further, we created a visualized nomogram based on a logistic regression selective model to classify the optimal patients. Last, both internal and external validation demonstrated our model a useful and stable one.

In our model, lobe tumor position, younger age were the top 2 strongest predictors of patient who would be benefit from primary tumor resection. As is described before, smaller tumor lesion, less node invasion, and younger patients were associated with better outcomes for stage IV NSCLC patients who underwent primary tumor surgery (14), these facts indicated that self-conditions are crucial for surgery in advanced stage patients, which might owing to relative lower difficulty for surgery and thus with a longer life expectancy.

We came up with a new process to establish the decision-making model that pick the optimal patients who could indeed prolong their survival time via primary tumor surgery. Firstly, patients who underwent surgery and not were collected and matched according to baseline character to eliminate the selection bias. Secondly, median CSS time of surgery group and non-surgery group were calculated, we found after PSM, surgery group still own longer CSS time (19 vs. 9 months). Third, we divided the surgery patients into two group by the median CSS time of non-surgery patients (9 months), through which we believe that patients received primary tumor surgery that lived longer than 9 months could really benefit from the operation. Last, we compare the pre-operative baseline characters between surgery performed patients who live longer than 9 month (beneficial group) and who not (non-beneficial group), and then logistic model was used to construct the model. After such screening process, metastatic NSCLC patients who were really suitable for surgery were identified.

In 2017, Elizabeth A. David and colleagues develop a quantitative model predictive of selection for surgical treatment for patients with advanced-stage NSCLC (15). Using detailed clinical variables, readily available at the time of treatment decision making for advanced-stage NSCLC patients, they developed the SSS (a score system) that predicts the use of surgical intervention for these patients and found evidence that a therapeutic effect of resection remains after accounting for this SSS. Overall, AJCC metastatic status, AJCC nodal status, and age were the strongest predictors of patient selection for an operation in the SSS. While in our model, tumor position, age, differentiation and T stage were orderly indicated to be the most factor that could decide whether one could benefit from the primary site resection. The difference here demonstrated that patients who underwent primary resection might not be the optimal candidates, and also indicated the limitation of surgeons’ recommendations.

In 2017, updates to the TNM staging system have reclassified metastatic NSCLC into M1a (separate tumor nodule in a contralateral lobe; tumor with pleural or pericardial nodules; or malignant pleural or pericardial effusion); M1b (single extra-thoracic metastasis in a single organ); or M1c (multiple extra-thoracic metastases in one or several organs) (16). We reclassified stage according to AJCC 8th edition, but materials in SEER database were not enough to differentiate M1c. A study suggested that both M1b and M1c patients would benefit from local surgery (17). For patients with only one metastatic lesion, 5-year survival was approximately 50% after surgery, compared to only 17% for patients undergoing surgery with five metastatic lesions (P<0.001) (4). In our analysis, M1a patients seemed would gain more from primary tumor resection than M1b, thus we have the reason to believe that M1c patients is not as suitable as M1b patients to receive surgery, systemic therapy was considered as priority in M1c patients.

We acknowledge some limitations in our study. First, the SEER database is lack of information on comorbidities that might lead to a selection bias for treatment choice. Second, the database also does not contain the number of metastatic sites, which have been shown to be an important predictors of surgery in IV stage patients (18). Third, data on systemic therapies were unavailable in SEER data base. The lack of information regarding target therapies and immunotherapy represents a major limitation of the current study. We are unable to answer whether before or after the systemic treatment would the primary tumor resection bring more survival profit. However, data from single center are limited and also lack of representative to evaluate such clinical question. SEER is the only comprehensive population-based database and open access to worldwide represents an ideal approach to study the survival of patients in such setting.

Conclusions

We put forward a novel process to screening stage IV NSCLC patients who would really benefit from surgery. Our predictive model demonstrated that primary tumor resection had beneficial impact on metastatic NSCLC patients, and only specific patients would benefit from surgery and that the potential benefit would vary based on primary tumor characteristics. To be specified, in stage IV NSCLC patients, <60 years old, female, adenocarcinoma, well differentiation, tumor in lobe, T1-2, N0 and M1a NSCLC patients would be potentially benefit more from primary tumor resection. Once validated in a prospective cohort, our predictive model can be of great help in selecting the optimal candidates for primary tumor resection among metastatic NSCLC patients.

The article’s supplementary files as