Abstract

With the increase of long-term primary lung cancer survivors, studies focused on metachronous second primary lung cancer (SPLC) have become very urgent. This study aimed to develop a prognostic nomogram and determine therapeutic options of cancer-specific death for patients with metachronous SPLC with and without the competing risk of other-specific death. Study population came from the SEER-18 database between 2006 and 2016. According to the clinical practice guideline of SPLC, the interval time of IPLC and metachronous SPLC was set to 4 years. We constructed nomograms with Lasso + Cox regression model and competing risk model to predict the prognosis and identify therapeutic options of metachronous SPLC patients with the assessment of model performance by the C-index, calibration plot, and decision curve analysis. In addition, two subgroup analyses stratified by histology and tumor size were used to better select therapeutic options for a certain population. 1300 patients with metachronous SPLC were incorporated in this study with 50.1% of the 5-year cumulative incidence in cancer-specific death. Compared with Lasso + Cox regression analysis, competing risk analysis harbored a higher C-index (0.811 vs. 0.76) and better net benefit in predicting cancer-specific death of metachronous SPLC. Two statistical analyses suggested that surgery alone was a preferentially therapeutic option of metachronous SPLC, whereas the effect of surgery + radiation in treating metachronous SPLC was similar to radiation alone. Subgroup analyses indicated that patients with metachronous SPLC were considered receiving different therapeutic options in different histology and tumor size but preferred to receive surgical treatment as the first choice. For primary lung cancer survivors, aggressive surgical treatment was the first-line selection of metachronous SPLC, followed by radiation alone, surgery + radiation, and no surgery + radiation.

1. Introduction

Because of air pollution and tobacco use in the past decades, lung cancer has become one of the most common carcinomas with high morbidity and mortality worldwide [1, 2]. However, the survival time of lung cancer is gradually prolonged due to the advancements in cancer treatment and the CT screen of early lung cancer. Low volume CT screening was associated with significantly lower risk of lung cancer death compared with chest radiography [3] and no screening [4]. Therefore, studies on metachronous second primary lung cancer (SPLC) of long-term lung cancer survivors are urgently needed. Current studies reported approximately 1% to 2% annual incidence of lung cancer among patients with initial primary lung cancer (IPLC), which was fourfold to sixfold higher than those with no history of lung cancer [57]. The cumulative risk of SPLC can reach 20.2% at 10 years postsurgery (95% confidence interval (CI): 15.3–23.2), and 25.2% (95% CI: 19.4–31.3) at 14 years postsurgery [8]. Prognostic risk factors of SPLC need to become the focus of clinicians’ research.

Previous studies reported some independent risk factors of overall survival and cancer-specific death among patients with SPLC, including tumor stage [9], tumor size [10, 11], lymph node metastasis [12], and surgery [1315]. Compared with IPLC, SPLC treatment needs to take into full consideration the influence of surgery on IPLC. Most studies about overall survival and cancer-specific death among patients with SPLC are based on the Kaplan–Meier method and multivariate Cox proportional regression analysis [915]. For old patients with SPLC, the competing risks of cardiovascular and cerebrovascular diseases are inevitable. Competing risk regression analysis is more available to determine the independent risk factors of cancer-specific death among the patients with SPLC in the presence of other-specific death. In addition, a competing risk nomogram can provide a graph of the mathematical model to predict special endpoints on the basis of various biological and clinical variables with relation to cancer-specific death. There is no published study about competing risk analysis and nomogram of cancer-specific death among patients with metachronous SPLC as yet.

Our first study aim is to obtain the 1-year and 5-year cumulative incidence of cancer-specific death and other-specific death among patients with metachronous SPLC by using competing risk regression analysis. Our second study aim is to determine independent prognostic factors of cancer-specific death with and without the competing risk of other-specific death. Subsequently, we construct nomograms with the Cox regression model and competing risk model to predict the prognosis of metachronous SPLC patients and assess model performance by the C-index, calibration plot, and decision curve analysis. In addition, subgroup analysis is performed to better select therapeutic options for a certain population.

2. Methods

2.1. Study Population

Surveillance, Epidemiology, and End Results (SEER) database is a representative population-based database of the National Cancer Institute, which covers 18 states across the United States and approximately 28% of the U.S. population with the strictest data-quality indicators and consistent criteria [16]. In our study, we extracted the data about metachronous SPLC between 2006 and 2016 from the SEER-18 database(SEER ID: huzg) released in April 2019 using the SEER Stat software (Version 8.3.4, https://seer.cancer/gov/seerstat/). The definition of metachronous SPLC was proposed by Martini and Melamed [17]: (1) interval time between IPLC and SPLC ≥2 years; (2) different histology between IPLC and SPLC; (3) When histology between IPLC and SPLC was the same, IPLC and SPLC must origin in situ or locate in different lobes or segments, with no positive intervening lymph nodes and no evidence of metastasis. Because an interval time of 2 to 4 years was considered a gray area [18], the interval time of IPLC and SPLC was set to 4 years in our study. Patients with unknown histological types in IPLC and SPLC were excluded. According to the International Classification of Diseases of Oncology (3rd edition), the histology of IPLC and SPLC was crudely classified into five groups: small cell carcinoma, large cell carcinoma, squamous cell carcinoma, adenocarcinoma, and others. Because partial patients came from 2004 to 2010, American Joint Committee on Cancer sixth Edition staging system was used to perform tumor stage of IPLC and SPLC in our study. No other lung cancer was diagnosed before the IPLC and after SPLC. Detailed surgical options of IPLC and SPLC were not included in the recently released SEER-18 database. Therefore, surgical options of metachronous SPLC were crudely classified into four groups (no vs. surgery alone vs. surgery + radiation vs. radiation alone) according to” Radiation sequence with surgery”, “Radiation recode,” and “Reason no cancer-directed surgery.” Tumor size (≤2 cm) was found to be an independent risk factor of metachronous SPLC [10, 13, 15]; thus we have taken the same classification in terms of tumor size (≤2 cm vs > 2 cm). The survival status of patients with metachronous SPLC was classified into alive, cancer-specific death, and other-specific death with the final follow-up time of Dec 2016. In addition, we collected some patients’ demographics, basic characteristics of IPLC and metachronous SPLC as potential confounding factors, including age (SPLC), sex, race, marital status at diagnosis (SPLC), laterality (IPLC), laterality (SPLC), primary Site (IPLC), primary Site (SPLC), AJCC T 6th ed (SPLC), AJCC N 6th ed (SPLC), AJCC M 6th ed (SPLC), chemotherapy record (IPLC), chemotherapy record (SPLC), SEER Combined Summary Stage (SPLC), and survival time (SPLC).

2.2. Ethics

There were no patients involved in the recruitment and conduct of our study. Therefore, our study was deemed exempt for review by the Institutional Review Board at China, Three Gorges University.

2.3. Statistical Analysis

Our study population was divided into three groups according to survival status. Demographic and clinical characteristics of the three groups were summarized with counts and frequencies, and tested with chi-square tests. The normality of distribution of the data was tested by chi-square goodness-of-fit and Kolmogorov–Smirnov tests. Cumulative incidence function (CIF) was used to estimate the probability of cancer-specific death and other-specific death during follow-up by using Fine and Gray’s competing risk regression analysis. 23 variables were regarded as the potential risk factors of cancer-specific death in our study. For high-dimensional data, the least absolute shrinkage and selection operator (LASSO) method with 10-fold cross-validation is more available to select the useful predictive features from the primary data than conventional univariate regression analysis [19, 20]. Subsequently, multivariate Cox regression analysis with proportional hazard model was used to determine risk factors and develop prognostic nomogram of cancer-specific death. We checked the proportional hazards assumption by using statistical tests and graphical diagnostics on the base of the scaled Schoenfeld residuals. We also compared the difference of cancer-specific death stratified by each potential risk factor through a competing risk analysis with proportional subdistribution hazard model. Before multivariate competing risk analysis, we evaluated the problem of collinearity of all potential risk factors by testing tolerance and variance inflation factor. When the tolerance of each variable was less than 0.1 and variance inflation factor was greater than 5, this variable was considered omitted from this study [21]. Based on multivariate competing risk regression analysis with proportional subdistribution hazard model [22], we examined the independent risk factors of cancer-specific death and constructed a competing risk nomogram to predict 1-year and 5-year probabilities of cancer-specific death. Akaike information criterion and the Bayesian information criterion (BIC) were used for competing risk regression analysis. We estimated the discrimination performance of prognostic nomogram by Concordance index (C-index) values and graphical calibration by plotting the observed rates against the nomogram-predicted probabilities through a bootstrap cross-validation approach with 1000 resamples [22]. The decision curve analysis for survival data with competing risk was used to evaluate the prediction model from the perspective of clinical consequences by calculating the net benefit [23]. When the prediction model produced a larger net benefit, the clinical use of decision curve analysis could effectively discriminate false-positive and false-negative results. In addition, we performed two subgroup analyses of tumor size and histology to better characterize the prognostic differences between chemotherapy and radiation/surgery. The subdistribution hazard ratio (sHR) with 95% CIs was estimated to evaluate the differences, and a two-tailed was considered statistically significant.

All statistical analyses were performed by using R version 3.4.2 software and Stata 14. Cumulative incidence function and the construction of competing risk nomogram were based on the following R packages: “cmprsk”, “rms”, and “mstate”. R package “riskRegression” and “pec” were used to construct calibration plot [22].

3. Results

3.1. Patient Characteristics

A total of 1300 metachronous SPLC patients with IPLC and an interval time of ≥4 years were incorporated in the study according to our selection criteria. The median age at diagnosis of metachronous SPLC was 73 years (interquartile range (IQR): 67–99 years). The majority of the study population was female (58.2%) and White (81.5%). The upper lobe (52.8%) was the most common site of metachronous SPLC, followed by the lower lobe (36.8%) and other areas (10.4%). Similar to IPLC, the main histology types of metachronous SPLC were adenocarcinoma (55.3%) and squamous cell carcinoma (24.5%). Among the total number of patients, 75.2% showed different histologies for IPLC and metachronous SPLC. There was 62.4% of metachronous SPLC with stage I. A total of 1014 patients (78%) had no positive lymph node involvement. The median interval time between IPLC and metachronous SPLC was 69 months (IQR: 57–86 months). The most common treatment of patients with metachronous SPLC was surgery alone (39.6%), followed by radiation alone (30.3%), no surgery/radiation (20.3%), and surgery + radiation (9.8%). Until Dec 2016, 582 patients of metachronous SPLC were alive, with 538 cancer-specific deaths and 180 other-specific deaths. The median follow-up time of individuals with metachronous SPLC was 30 months (IQR: 19–46 months) in the alive group, followed by the other-specific death group (17 months, IQR:7.0–31.2 months), and the cancer-specific death group (12 months, IQR:5.0–24 months). Detailed demographic and clinical characteristics of three survival status groups were shown in Table 1. The overall 1-year, 3-year, and 5-year mortality of metachronous SPLC were 25.2%, 52.7%, and 68.2%, respectively. When considering the competing risk of other-specific death, the 1-year, 3-year, and 5-year cumulative incidence of cancer-specific death were 19.8%, 40.3%, and 50.1%, respectively (see Figure 1).

3.2. The Construction and Assessment of Prognostic Nomogram of Cancer-Specific Death without Regard to Competing Risk

Lasso analysis with Cox regression model observed 16 potential risk factors of cancer-specific death with 0.84 of the area under the receiver operating characteristic curve (see Figure S1). Multivariate Cox regression analysis further found that the independent risk factors of cancer-specific death comprised seven variables, including age, sex, marital status, histology (SPLC), AJCC N (6 eth, SPLC), AJCC T (6 eth, SPLC), and radiation/surgery (SPLC) (see Table S1). Unmarried old male with small cell carcinoma and mediastinal lymph node involvement harbored worse clinical outcomes. Based on these risk factors, we constructed a prognostic nomogram of cancer-specific death with the Cox regression model (see Figure S2). Histology (SPLC) and AJCC N (6 eth, SPLC) were two main risk factors of cancer-specific death. The C-index of this nomogram was 0.736 (95% CI: 0.691–0.781) with an acceptable calibrated plot (see Figure S3). We also performed decision curve analysis to assess the clinical net benefit of the prognostic nomogram, which demonstrated a wide range of thresholds for predicting 1-year probability (9%–76%) and 5-year probability (29%–91%) of cancer-specific death (see Figure S4).

3.3. The Construction and Assessment of Prognostic Nomogram of Cancer-Specific Death in the Present of Competing Risk

Univariate and multivariate competing risk regression analyses (see Table 2) suggested the independent risk factors of cancer-specific death as follow: age, sex, marital status, tumor size, histology (SPLC), AJCC N (6 eth, SPLC), tumor stage (SPLC), and radiation/surgery (SPLC). Surgery alone (no radiation/surgery vs surgery alone, sHR = 0.41, 95%CI: 0.3–0.55, ) and radiation alone (no radiation/surgery vs radiation alone, sHR = 0.60, 95% CI: 0.46–0.78, ) were associated with significantly lower risk of cancer-specific death compared with no radiation/surgery. We integrated eight independent risk factors to develop the prognostic competing risk nomogram of 1-year and 5-year probability of cancer-specific death. The C-index value was 0.811 (95% CI: 0.805–0.817) with an acceptable calibrated plot (see Figure S5), which seemed to be higher than that of prognostic nomogram with Cox regression model (C-index = 0.736, 95% CI: 0.691–0.781). Through this competing risk nomogram, we visually estimated the contribution of each included variable in cancer-specific death (see Figure 2). Of all eight variables, histology of metachronous SPLC harbored the biggest contribution to cancer-specific death, especially small cell carcinoma. In terms of tumor stage, stage IV metachronous SPLC was associated with the highest risk of cancer-specific death, followed by stage III, stage II, and stage I. At the aspect of AJCC N, 6th ed (SPLC), the mediastinal involvement (N2) seemingly had the significantly bigger contribution to cancer-specific death compared with N0. By calculating the sum of points of the eight variables, we predicted the 1-year and 5-year probabilities of cancer-specific death among patients with metachronous SPLC. Moreover, decision curve analysis demonstrated that our competing risk nomogram produced clinical net benefit in all ranges for predicting 1-year probability of cancer-specific death and in a wide range of thresholds (27%–93%) for predicting 5-year probability (see Figure 3).

3.4. Subgroup Analyses of Histology and Tumor Size in Competing Risk Analyses

For nonsmall cell carcinoma, surgery alone, radiation alone, and radiation + surgery had significant therapeutic effects in cancer-specific death compared with no radiation/surgery (see Table 3). However, patients with small cell carcinoma might consider receiving surgery alone (no radiation/surgery vs. surgery alone, sHR = 0.28, 95% CI: 0.11–0.76, ) and radiation alone (no radiation/surgery vs radiation alone, sHR = 0.41, 95% CI: 0.18–0.90, ), but not surgery + radiation (no radiation/surgery vs. surgery + radiation, sHR = 0.85, 95% CI: 0.28–2.55, ).

The authors also stratified the patients by tumor size of metachronous SPLC (≤2 cm vs >2 cm). Chemotherapy had no significant impact on cancer-specific death whether tumor size was less than 2 cm or larger than 2 cm. When tumor size was less than 2 cm, only surgery alone was associated with a lower risk of cancer-specific death than no surgery + radiation (no radiation/surgery vs. surgery alone, sHR = 0.52, 95% CI: 0.31–0.87, ). When tumor size was more than 2 cm, the clinical benefit of cancer-specific death was significant in the surgery alone (no radiation/surgery vs. surgery alone, sHR = 0.33, 95% CI: 0.23–0.48, ) and radiation alone (no radiation/surgery vs. radiation alone, sHR = 0.63, 95% CI: 0.47–0.86, ) groups compared with the no radiation/surgery group.

4. Discussions

This study estimated the cumulative incidence of cancer-specific death in patients with metachronous SPLC, differentiated the influence of other-specific death, and determined potential risk factors of cancer-specific death. The 5-year cumulative incidence of cancer-specific death and other-specific death were 50.1% and 18.1%, respectively. Our study also found that the C-index of competing risk model (0.811, 95% CI:0.805–0.817) was superior to that of Cox regression model (0.736, 95% CI: 0.691–0.781). In addition, our subgroup analyses indicated that the patients of metachronous SPLC with small cell carcinoma and nonsmall cell carcinoma initially preferred to receive surgery alone. When tumor size was less than 2 cm, surgery alone was considered as the main therapeutic option of metachronous SPLC. When tumor size was more than 2 cm, surgery alone and radiation alone were also proposed as therapeutic options for metachronous SPLC along with the preferred method of surgery alone.

Through multivariate competing risk analysis and nomogram points, the authors found some independent risk factors of cancer-specific death among patients with metachronous SPLC. In terms of histology (SPLC), small cell carcinoma seemingly had the biggest contribution to cancer-specific death, followed by large cell carcinoma, squamous cell carcinoma, other, and adenocarcinoma. The vast majority of existing studies focused on nonsmall cell carcinoma of metachronous SPLC [13, 24] and the comparison of operational modes of early SPLC [15, 25, 26]. Our study suggested that surgery alone and radiation alone might be therapeutic options for metachronous SPLC with small cell carcinoma. For metachronous SPLC with nonsmall cell carcinoma, the effect of surgery alone on cancer-specific death seemed to be superior to those of other therapeutic options, especially when tumor size was ≥2 cm. Our study also showed that radiation alone had significant clinical benefit compared with no radiation/surgery (no radiation/surgery vs radiation alone, sHR = 0.63, 95% CI: 0.47–0.86, ) when tumor size was larger than 2 cm. Taioli and the colleagues conducted the first study on the effect of radiation and surgery on SPLC survival [27]. The results obtained by multivariate Cox proportional hazards model and propensity score analysis showed that patients with SPLC who underwent surgery alone demonstrated significantly longer survival than those who underwent radiation alone and no radiation/surgery. Another propensity score analysis from a Chinese population showed that the 5-year overall and progression-free survival of metachronous SPLC with surgical treatment were significantly better than those without surgery [14]. No published study has determined the effect of radiation + surgery on cancer-specific death by using competing risk analysis. Two other important risk factors of CSD, namely, tumor stage and lymph node metastasis, were mentioned in our study (see Figure 2). In multivariate competing risk analysis, stage I was associated with significantly lower cancer-specific mortality than stage III (stage I vs stage III, sHR = 1.81, 95% CI: 1.03–3.19, ) and stage IV (stage I vs stage IV, sHR = 2.17, 95% CI: 1.23–4.18, ). Three previous studies [9, 10, 28] also reported that TNM stage was an independent risk factor of survival in metachronous SPLC, especially stage I. Lee and the colleagues even declared that TNM stage is the most important determinant of survival in metachronous lung cancer [9]. Our study found that mediastinal involvement (N2) of metachronous SPLC was associated with significantly worse survival than N0 (N0 vs. N2, sHR = 1.75, 95% CI: 1.15–2.87, ), which was similar to the study results of Riquet et al. [12] and Yang et al. [26]. By multivariate Cox regression analysis, tumor size (≤2 cm) was regarded as an important risk factor of overall and cancer-specific survival in some studies [10, 15, 24, 26]. In addition, tumor size (≤2 cm) was an important reference for operation modes and can be a basis for the performance of sublobar resection [24, 26]. Our competing risk analysis suggested that a smaller tumor size was associated with a lower risk of cancer-specific death (≤2 cm vs >2 cm, sHR = 1.39, 95% CI: 1.12–1.73, ). In subgroup analysis on the basis of tumor size, only surgery alone harbored a significant effect on cancer-specific mortality when tumor size was less than 2 cm. When tumor size was larger than 2 cm, surgery alone and radiation alone could obtain significant clinical benefits of cancer-specific death compared with no radiation/surgery. There are no published studies about the comparison between the effects of radiation and surgery on cancer-specific death of metachronous SPLC stratified by tumor size. The competing risk nomagram indicated that the other three risk factors (sex, age, marital status) had relatively little contributions to cancer-specific death. Our study also observed that histology between IPLC and metachronous SPLC (different vs same) and surgery/radiation of IPLC (no vs. yes) had no significant effect on cancer-specific death, which were similar to the previous studies [12, 16, 26].

Our study had the following strengths. We included 1300 patients from 18 states across the United States; this was approximately 28% of the U.S. population. Our sample size effectively decreased the selective bias of a small sample size and single-center study. We examined the independent risk factors of cancer-specific death by simultaneously using Cox regression analysis and competing risk regression analysis. We also developed and compared prognostic nomograms to guide clinicians’ risk evaluations for cancer-specific mortality among patients with metachronous SPLC. In addition, our subgroup competing risk analyses provided important references for therapeutic options in cases with different histologies and for selecting between radiation and surgery in metachronous SPLC cases with different tumor sizes.

Certainly, our study also had some limitations. Firstly, some important details, such as smoking status and cardiopulmonary function, were not included in the SEER database, which potentially affected our study results. Patients’ clinical backgrounds might affect the selection of therapeutic options. Surgery selects the fittest patients. Radiation selects intermediary fit. Undocumented patient’s characteristics are likely to result in selection bias. Secondly, the newly released SEER-18 did not reveal surgical treatments. The authors were not able to assess the advantages and disadvantages of surgical treatments in specific populations. Thirdly, the lack of external validation set potentially affected the confidence of our study results. Although competing risk analysis with proportional subdistribution hazards regression can obtain unbiased estimates of the risk of cancer-specific death, potential bias was still inevitable due to the retrospective nature of the study.

5. Conclusion

Our study provided the cumulative incidence of cancer-specific death among metachronous SPLC on the basis of a large-sample study population with considering the competing risk of other-cancer death. A prognostic nomogram with a competing risk model was constructed to predict the 1-year and 5-year probabilities of cancer-specific death with acceptable performance and decision curve analysis. Our subgroup competing risk analyses provided evidence-based references for therapeutic options of metachronous SPLC in specific populations. Future prospective studies and random controlled studies should be conducted to determine the choice of radiation and surgery in certain patients with metachronous SPLC.

Data Availability

The data underlying this study were obtained from open SEER databases. All relevant data are within the paper and its Supplementary Materials.

Ethical Approval

This study was deemed exempt for review by the Institutional Review Board at China, Three Gorges University.

Conflicts of Interest

This was not an industry supported study. The authors declare that there are no conflicts of interest.

Authors’ Contributions

(I) Conception and design: Jiahui Wang and Xinyu Song. (II) Administrative support: Zhen Chen. (III) Provision of study materials or patients: Daokui Xia and Zhigang Hu. (IV) Collection and assembly of data: Jiahui Wang and Zhigang Hu. (V) Data analysis and interpretation: Xinyu Song. (VI) Manuscript writing: all authors. (VII) Final approval of manuscript: all authors.

Acknowledgments

The work was supported by the Natural Science Foundation of Hubei Province CN (No.2018CFB643).

Supplementary Materials

Table S1: risk factors of cancer cause death in prognostic nomogram with Cox regression model among patients with metachronous second primary lung cancer. S1 Figure: the screen of independent risk factors of cancer-specific death using the least absolute shrinkage and selection operator (LASSO) analysis with Cox regression model. (a) Tuning parameter selection in the LASSO model used 10-fold cross-validation via minimum criteria. (b) LASSO coefficient profiles of potential risk factors. (c) The area under the receiver operating characteristic curve (AUC = 0.84). S2 Figure: the prognostic nomogram for predicting the 1-year and 5-year probability of cancer-specific death among metachronous SPLC by using Cox regressiong model. Total point values were independently calculated for each cause of death and then applied to the corresponding probability scale at the bottom of each. S3 Figure: decision curve analysis for cancer-specific death model among metachronous SPLC by using Cox regression model. (a) 1-year probability of cancer-specific death; (b) 5-year probability of cancer-specific death. S4 Figure: calibration plot for cancer-specific death model among metachronous SPLC by using Cox regression model. S5 Figure: calibration plot for cancer-specific death model among metachronous SPLC by using competing risk model. The solid line represents equality between the predicted and observed probabilities. With the dots close to the solid line, the plots reveal excellent agreement between the nomogram-predicted probabilities and actual observations. (Supplementary Materials)