Introduction
Lung cancer is one of the most common cancers in the world, and its mortality of 25% is the highest among all cancers [1]. Non–small cell lung cancer (NSCLC) accounts for approximately 85% of lung cancer cases and surgery is the mainstay treatment for patients with early-stage NSCLC [2]. Surgery is generally performed in patients with stage I to IIIA NSCLC, when the tumor is resectable in patients who are suitable for surgery [3]. The 5-year overall survival (OS) rate for NSCLC is approximately 50% for resected lung cancers [1] and approximately 30%–55% of patients with resected NSCLC eventually show relapse during follow-up [4]. There are several risk factors affecting OS or relapse-free survival (RFS) in patients with resected NSCLC. Tumor-related factors, including tumor-node-metastasis (TNM) staging, cell types, and host-related factors (i.e., sex, age, smoking history) are predictive of OS in patients with NSCLC [5]. In addition, TNM staging, tumor markers, and performance status are reported as significant factors in predicting RFS after surgical resection in patients with NSCLC [5].
In patients who have undergone NSCLC surgery, events such as relapse or metastasis are called intermediate events because they are not fatal, whereas death is called a terminal event. Studies have typically focused on whether or not a patient survives without experiencing relapse or death when analyzing cancer-related survival data. This concept is called disease-free survival (DFS). In DFS studies, the intermediate event and terminal event are considered to have a competing relationship, survival time is defined as the time of occurrence of the first of the two events, and censoring is defined as whether either of the two events has occurred or not. On the other hand, the illness-death model (IDM) [6–8] considers the intermediate event and terminal event to have a semi-competing relationship, so it has the advantage of observing the survival time from the occurrence of the intermediate event to the occurrence of the terminal event, which cannot be observed in DFS.
The purpose of this study was to identify risk factors affecting relapse and death via the IDM among NSCLC patients and evaluate personalized predictive probabilities of disease outcomes to help in establishing follow-up plans or treatment strategies.
Materials and Methods
1. Study population
A total of 917 patients who underwent curative surgery for primary lung cancer from 2010 to 2018 were initially included in this study. Exclusion criteria were as follows: (1) sublobar resection, (2) the presence of synchronous or metachronous lung cancer, (3) an unknown origin of cancer recurrence in patients with double primary cancer, (4) the presence of distant metastases, (5) other than adenocarcinoma or squamous cell carcinoma, and (6) pre-invasive cancers. A total of 612 consecutive patients with resected lung adenocarcinoma or squamous cell carcinoma were ultimately included in the study (Fig. 1).
2. Data processing
Patient demographics, history of malignancy other than lung cancer, family history of lung cancer, smoking history (never smokers vs. previous or current smokers), cancer location (upper lobes, including the middle lobe vs. lower lobes), cancer type (adenocarcinoma vs. squamous cell carcinoma), TNM stage, operation method (lobectomy vs. bi-lobectomy or pneumonectomy), treatment method, and the presence of relapse were recorded from electronic medical records and radiology reports. TNM stage was determined according to the eighth edition of the American Joint Commission on Cancer staging system for lung cancer [9]. The T category was determined based on pathologic reports following lobectomy, and the N category was determined by lymph node dissection or endobronchial ultrasound-guided biopsy.
Time to relapse was defined as the time from the date of operation to the date of the first recorded evidence of intra-thoracic or distant metastasis as confirmed by imaging or histology. Loss to follow-up was defined as a case where imaging follow-up was not performed after December 2018 (1 year before the end of the study). The time of censoring was determined as the date of the last imaging evaluation. Date of death and cause of death data were obtained via linkage to the Korean Statistical Information Service.
3. Statistical analysis
For analyses of the DFS model and the IDM, proportional hazard models were used to examine the associations of covariates (age, sex, smoking history, family history of lung cancer, history of other malignancies, nodule location, operation method, cancer type, treatment method, interaction of cancer types and treatment methods, TNM stage) with relapse or death. For the analysis of DFS, the primary endpoint was defined as the time of relapse or death. The IDM was fitted to estimate the intensities in transitions from surgery to relapse (0→1), from surgery to death (0→2), and from relapse to death (1→2), as illustrated in Fig. 2. These three transitions could occur at any time until the end of follow-up. The most widely used R package for analyzing DFS models is survival, but this package can only handle right-censored data. Therefore, the icensReg [10] R package was used instead of survival because the event of interest, such as relapse, was interval censored in the NSCLC data. In addition, the most widely used R package to analyze multi-state models is msm, but this package can only analyze right-censored data with intermediate events. Therefore, the SmoothHazard [11] R package was used to deal with interval-censored intermediate events.
In both the DFS and ID models, interval censoring was considered as the assessment of relapse was not continuous; the exact date of relapse could have occurred in the interval between two consecutive clinical examinations. As for baseline transition intensities, a Weibull distribution was assumed in both the DFS and ID models. Multivariable analyses in both the DFS and ID models were conducted to select the most parsimonious model that would reduce overfitting. Risk factors included in initial multivariable analyses were selected based on a p-value < 0.2 in TNM-adjusted univariable analyses because TNM stage is considered a primary risk factor. The Akaike information criterion (AIC) was used to decide on final models (i.e., models with the smallest AIC value), thus indicating the best model through a backward elimination strategy.
The IDM can discern four different predictive probabilities for patients who underwent surgery, while the DFS model includes only two predictions. For the former, the probabilities of being alive after relapse, being dead after relapse, being dead without relapse, and being alive without relapse are predictable at each elapsed time since surgery; for the latter, we calculated the probabilities of being relapsed or dead and being alive without experiencing any event. As an illustration, we examined a patient with the following risk factors: age, 65 years; cancer type, adenocarcinoma; sex, male; treatment method, operation only; history of malignancy other than lung cancer, no; smoking history, previous or current smoker; nodule location, upper lobes, and middle lobes; and operation method, lobectomy. Those values correspond to the sample median or mode of each risk factor. We estimated the predictive probabilities of survivorship and relapse by TNM stage (1A, 2A, or 3A) in order to investigate characteristics within each stage and to compare changes in trends according to TNM stage.
Results
1. Patient characteristics
Among the 612 included patients with resected lung adeno-carcinoma or squamous cell carcinoma, 69.3% (424/612) were men (mean age±standard deviation, 66.3±9.6 years) and 30.7% (188/612) were women (mean age±standard deviation, 65.1±8.9 years). Sixteen percent (99/612) had a history of malignancy other than lung cancer, 4% (26/612) had a family history of lung cancer, and 56% (345/612) were previous or current smokers. Forty percent (246/612) of the lung tumors were located in the lower lobes. With respect to TNM staging, 48% (294/612) of cases were 1A, 12.3% (75/612) were 1B, 5.9% (36/612) were 2A, 15.8% (97/612) were 2B, 14.4% (88/612) were 3A, 3.4% (21/612) were 3B, and 0.2% (1/612) were 3C. A combined TNM stage (3B and 3C) was used in the analysis because there was only one patient with stage 3C cancer. Sixty-two percent (382/612) of the patients had adenocarcinoma, and 38% (230/612) had squamous cell carcinoma. Fifty-three percent (323/612) underwent surgery without additional treatment modalities, and 47% (289/612) received chemotherapy or radiation therapy in addition to surgery. With respect to relapse, 51 of the 612 patients (8.3%) were lost to follow-up, 68% (415/612) did not experience relapse, and 24% did experience relapse (146/612). Seventy percent (426/612) were alive, 25% (155/612) of the included patients died from lung cancer, and 5% (31/612) died from causes other than lung cancer during the course of the study. Detailed results are described in Table 1.
2. Comparative risks: DFS and ID models
In univariable analyses of the DFS and ID models adjusted for TNM stage, the p-values of all variables except family history of lung cancer were < 0.20. Those variables were considered for multivariable analyses. S1 Table shows the detailed results of univariable analyses.
The DFS model accounting for age, smoking history, history of other malignancies, operation method, TNM stage, cancer type, treatment method, and the interaction of cancer type and treatment method had the lowest AIC value (AIC, 980.4) (Table 2).
The IDM (AIC, 1,868.3) included age, smoking history, history of other malignancies, nodule location, operation method, TNM stage, cancer type, treatment method, and the interaction of cancer type and treatment method (Table 3). Pneumonectomy against lobectomy had a statistically significantly higher risk within all three transitions (relapse, 0→1: hazard ratio [HR], 1.97; p=0.007; death without relapse, 0→2: HR, 2.58; p=0.009; death following relapse, 1→2: HR, 1.84; p=0.039). Risks for relapse and death without relapse increased substantially with an increase in TNM stage. Smoking history was a statistically significant predictor of death without relapse, with an HR of 4.05 (p < 0.001), though there was no statistically significant association with relapse or death after relapse. History of other malignancies was observed as a statistically significant predictor of relapse, though not of other transitions. There was a statistically significant interaction effect between cancer type and treatment method for relapse (p < 0.001) and death without relapse (p < 0.001), as in the multivariable DFS analysis (p < 0.001). For death following relapse, adenocarcinoma showed substantially less risk compared to squamous cell carcinoma (HR, 0.52; p=0.006). Both the DFS and ID models indicated that family history of lung cancer and sex were not statistically significant predictors of disease outcomes. The IDM revealed a statistically significant relative risk of nodule location on disease outcomes that was not observed in the DFS model.
3. Predictions from the DFS and ID models
Fig. 3 displays predictive probabilities by TNM stage, as estimated from both the ID and DFS models. For the IDM, the height of the sum of the two bottom regions represents the probability of relapse, and the sum of the two middle regions represents the probability of death. Moreover, the height of the sum of the colored regions in the IDM corresponds to the probability of death or relapse in the DFS model.
As expected, a patient diagnosed with a TNM stage of 1A had the lowest estimated probability of being alive after relapse (5-year probability, 5.1%) as well as of death after relapse (2.7%) and death without relapse (5.5%). Probabilities were highest for patients with a TNM stage of 3A (with 5-year probabilities of 14.4%, 9.1%, and 42.8%, respectively). Moreover, the sum of three 5-year probabilities in the IDM gave similar results for the probability of death or relapse as compared to the DFS model (1A, 13.3 vs. 14.2%; 2A, 32.1 vs. 38.5%; 3A, 66.3 vs. 67.8%).
Discussion
In our study, we demonstrated that it is possible to evaluate detailed predictive probabilities of disease outcomes via the IDM in NSCLC patients who received surgical resection. The main reason for this finding is that the IDM evaluates the occurrences of relapse and death simultaneously as well as in conjunction with interactions between these outcomes. As compared to the DFS model, the IDM can identify statistically significant risk factors for disease outcomes according to the patient’s disease course, including relapse, death without relapse, and death after relapse. To our knowledge, this is the first study to apply the IDM for an epidemiological investigation of lung cancer. Our findings facilitate an understanding of multi-state processes in the post-surgery follow-up period for lung cancer patients.
The IDM proposed in our study showed several strengths as compared to the DFS model. First, history of malignancy other than lung cancer (HR, 1.44; p=0.049) and smoking history (HR, 1.68; p=0.001) were statistically significant risk factors related to decreased DFS; the IDM demonstrated that a history of malignancy other than lung cancer was related to relapse (HR, 1.67; p=0.023) and smoking history was related to death without relapse (HR, 4.05; p < 0.001). Second, the location of the lung tumor was not a statistically significant risk factor in the DFS model; however, lower lobe location was positively associated with death without relapse (HR, 1.97; p=0.011) and negatively associated with relapse (HR, 0.63; p=0.009) according to the IDM. This finding suggests that the DFS model offsets these effects. Third, we observed a statistically significant interaction between cancer type and treatment methods in both the DFS and the IDM. However, the IDM could additionally differentiate the magnitude of interaction between cancer type and treatment methods in death without relapse from death after relapse.
Various events can occur in the disease course of a patient who has undergone surgery for NSCLC. Usually, researchers perform separate analyses for each event setting as a primary endpoint, considering the other endpoints as censoring. For example, studies of OS aim to estimate the probability of death by considering relapse as censoring. Conversely, evaluating RFS aims to estimate the probability of relapse by considering death as censoring. However, those separate analyses do not provide completely satisfying results due to a failure to reveal associations between relapse and death [12,13]. In contrast, the IDM simultaneously deals with endpoints focusing on transitioning from one state to another. Moreover, the IDM provides predictions of patient clinical prognoses at specific points in their relapse or death process. Several medical articles have estimated event-related predictive probabilities via the IDM in patients with breast cancer [12,14,15], ovarian cancer [16], or colon cancer [13]. However, to our knowledge, the IDM has not been implemented in patients with lung cancer.
As mentioned above, the IDM allowed us to explore the effects of risk factors on relapse leading to death. Regarding the NSCLC patients who experienced relapse in our study, we found that pneumonectomy and squamous cell carcinoma cancer subtype were unfavorably associated with death after relapse. Sekihara et al. [17] reported that female gender, adenocarcinoma histology, and absence of distant metastasis were favorably associated with post-recurrence survival (PRS), and Shimada et al. [18] revealed that adenocarcinoma showed a favorable PRS; however, an unfavorable PRS was observed among patients with pneumonectomy as well as adjuvant therapy. Our study identified risk factors from the IDM that were similar to previous results obtained with the PRS model, though these two models have different starting points. Specifically, the PRS model examines survival probability only among patients who experienced relapse, and the survival rate is calculated from the time of relapse [19]. In contrast, the IDM evaluates the probability of death after relapse starting from the time of lung cancer surgery. We believe that the IDM shows a more comprehensive view of risk factors for mortality as well as other patient outcomes experienced during lung cancer follow-up.
As shown in Fig. 3, the predictive probability of relapse or death obtained from the DFS model can be decomposed into three parts in the IDM (i.e., the predictive probabilities of being alive after relapse and of having died with or without relapse). This decomposition enables us to investigate the evolving process of events such as relapse and death [20–22]. For example, for a patient diagnosed with lung cancer with a TNM stage of 3A, there is a 67.8% probability of cancer relapse or death within five years of a lung cancer operation based on the results of the DFS model. We report the following separate five-year probabilities for survival after relapse, death after relapse, and death without relapse via the IDM: 14.4%, 9.1%, and 42.8%, respectively. Moreover, the sum of these probabilities corresponds to the predictive probability of relapse (23.5%) and the two mortality probabilities (51.9%) within 5 years following the operation.
Our study has several limitations. First, this was a retrospective single-center study and our results were evaluated within a small number of included patients. Our primary focus was to introduce the IDM and to compare it with the typical DFS model among patients with surgically resected NSCLC. As a further step, implementing the IDM within a larger multi-institutional cohort will help us to better under-stand the association between relapse and death among NSCLC patients. Second, our analysis was limited to patients with resected adenocarcinoma and squamous cell carcinoma, which are two of the most common histologic NSCLC types. Further studies including other cancer subtypes and advanced-stage lung cancers are warranted to gain a comprehensive understanding of this issue. Third, loss to follow-up was defined as a case where imaging follow-up was not performed as of a year before the end of the study, and patients who died of causes other than lung cancer were regarded as being censored at the time of death. Fourth, because the follow-up interval is longer in patients undergoing surgery at an earlier stage, there is a possibility that the date of relapse may not be accurate for early-stage patients, which had a total follow-up period of 8 years. Fifth, although treatment methods change before and after relapse in lung cancer patients, we assumed models in which the effects of treatment methods do not change during follow-up. However, if there are risk factors that may show altered effects before and after relapse, such as treatment methods, a model that considers the effects of these risk factors as time-varying would be more appropriate.
In conclusion, our results demonstrate the IDM can be used as a complementary tool to simultaneously evaluate the predictive probabilities of relapse and death during the post-operation follow-up of NSCLC patients. Furthermore, the IDM may help in establishing follow-up plans or treatment strategies according to individual risk prediction. We believe that the IDM provides a more comprehensive picture of risk factors as compared to the DFS model and may facilitate an improved understanding of the multi-state disease processes occurring among lung cancer patients.
PDF Links
PubReader
ePub Link
XML Download
Full text via DOI
Download Citation 
CrossRef TDM
Supplement
Print 
