Non-Small-Cell Lung Cancer Patients with a High Predicted Risk of Irradical Resection: Can Chemoradiotherapy Offer Similar Survival?

Purpose Irradical resection of non-small-cell lung cancer (NSCLC) is a detrimental prognostic factor. Recently, Rasing et al. presented an internationally validated risk score for pre-treatment prediction of irradical resection. We hypothesized that chemoradiation therapy (CRT) could serve as an alternative approach in patients with a high risk score and compared overall survival (OS) outcomes between surgery and CRT. Methods Patients from a population-based cohort with stage IIB–III NSCLC between 2015 and 2018 in The Netherlands were selected. Patients with a ‘Rasing score’ > 4 who underwent surgery were matched with patients who underwent CRT using 1:1 nearest-neighbor propensity score matching. The primary endpoint of OS was compared using a Kaplan–Meier analysis. Results In total, 2582 CRT and 638 surgery patients were eligible. After matching, 523 well-balanced pairs remained. Median OS in the CRT group was 27.5 months, compared with 45.6 months in the surgery group (HR 1.44, 95% CI 1.23–1.70, p < 0.001). The 114 surgical patients who underwent an R1–2 resection (21.8%) had a worse median OS than the CRT group (20.2 versus 27.5 months, HR 0.77, 95% CI 0.61–0.99, p = 0.039). Conclusion In NSCLC patients at high predicted risk of irradical resection, CRT appears to yield inferior survival compared with surgery. Therefore, choosing CRT instead of surgery cannot solely be based on the Rasing score. Since patients receiving an R1–2 resection do have detrimental outcomes compared with primary CRT, the treatment decision should be based on additional information, such as imaging features, comorbidities, patient preference, and the surgeon’s confidence in achieving an R0 resection. Supplementary Information The online version contains supplementary material available at 10.1245/s10434-021-10982-3.

Recently, a prediction score for irradical resection has been developed in the Netherlands and validated in the United States of America, by Rasing et al. 8 This model was developed using data from the Netherlands Cancer Registry, which resulted in a multivariable logistic regression model for the prediction of irradical resection with a good discriminative performance (external c-statistic 0.71). Predictive parameters included histology, clinical T-stage, clinical N-stage, planned extent of resection (e.g., lobectomy, pneumonectomy), and surgical approach (thoracoscopic, open). Patients with a Rasing score [ 4 were deemed at high risk, as their individual predicted and observed probability of an R1-2 resection was[13% in all these patients, and 19% on average. 8 We hypothesized that CRT could serve as an alternative approach in patients with a high risk of an irradical resection as predicted by the Rasing score. 8 Therefore, the primary aim of this study was to compare overall survival (OS) between surgery and CRT in those patients. A secondary aim was to identify potentially varying OS outcomes for CRT versus surgery across different patient subgroups.

MATERIALS AND METHODS
This study was a population-based retrospective cohort study. Institutional review board approval was obtained and the need for written informed consent was waived, as all patient information in the database used was anonymized.

Study Population
Surgery Group The surgery group used in this study was formed using a group of patients described in Rasing et al. 8 These patients were identified using the national database of the Netherlands Cancer Registry (managed by the Netherlands Comprehensive Cancer Organization). All newly diagnosed cancer cases in the Netherlands are registered in this database by independent data managers. The database contains in-depth patient, tumor, and treatment characteristics, which are gathered based on information from electronic medical records, the Pathological Anatomical National Automated Archive and the national registry of hospital discharge diagnoses and diagnosis-treatment combinations (DBC). All patients diagnosed with NSCLC between January 1, 2015, and December 31, 2018 and treated surgically were included for analysis.
The predicted probability (p) of an R1-2 resection was calculated using the following formula: 8 Chemoradiation Therapy Group The same data set and time period were used to identify patients who underwent CRT. Patients were considered to have received CRT if they received both chemo-and radiotherapy, with chemotherapy being administered first or within 10 days of the first radiotherapy fraction, and with \ 120 days between the first dosage of chemotherapy and radiotherapy. Exclusion criteria were age \ 18 years, cM1 stage, cN3 stage, cTis, cT0, overall clinical stage I or IIA, carcinoid histology, a time interval from diagnosis to the beginning of treatment of 0 days or [ 90 days, and a duration of radiotherapy of \ 26 or [ 50 days.

Study Parameters
Variables analyzed for both groups include age, sex, WHO performance status, previous history of malignancy, localization and lateralization of the tumor, tumor histology, cT and cN stages, and overall staging. For the surgery group, completeness of resection (R0 or R1-2), extent of surgery, surgical approach and whether adjuvant chemotherapy was administered were additional variables. For surgical approach, the intended approach was recorded (with conversion to open surgery counting as scopic approach). Type of chemoradiation therapy (sequential or concurrent) was an additional variable in the chemoradiation group. Outcome variables were vital status at the end of follow-up and time duration in days of follow-up from diagnosis. As the database used does not register the ethnicity or race of patients, no data regarding the distribution of ethnicity or race of patients was available for analysis.

Statistical Analysis
Continuous variables are presented with mean ± standard deviation (SD). Differences in baseline between the two groups in continuous variables were assessed with an unpaired T-test or Mann-Whitney U test, depending on the normality of the data distribution. Differences in nominal and ordinal categorical variables were assessed with a v 2 test and the Mann-Whitney U test, respectively.
Missing data in the chemoradiation group were assumed to be missing at random, and handled using multiple imputation by chained equations, creating 20 new data sets. The most representative imputation set was used, determined via survival modeling of all imputation sets separately and comparing them with the pooled OS model for all 20 imputed data sets.
The imputed data set of CRT patients was matched to the group of high-risk patients who had undergone surgery identified in the Rasing study using propensity score matching. A logistic regression model was performed to determine a propensity score for each patient in which the variables age, sex, WHO performance score, history of malignancy, lateralization and localization of the tumor, histology, and cT-and cN-staging were accounted for. Patients from both groups were matched (1:1) according to nearest-neighbor matching without replacement. Withinpair difference was minimized by setting a caliper of 0.1 of the standard deviation of the logit of the propensity score. Standardized mean differences (SMD) were calculated to check the balance of the match, with a threshold of 0.1 considered well balanced. Additionally, the groups were compared on all variables accounted for in the matching process, similarly to the baseline comparison, to further check the balance of the match. Imputation of missing data and propensity score matching were performed with ''mice'' and ''MatchIt'' packages in R 4.0.3 software (The R Foundation for Statistical Computing, Vienna, Austria).
To compare the matched groups, OS analysis was performed with Kaplan-Meier survival analysis. In addition to the main survival analysis, subgroup analyses were performed using Cox regression models to identify potential differences in hazard ratios across subgroups. To compare OS for CRT with those receiving irradical resection, Kaplan-Meier survival analysis was performed. Finally, to assess whether results would be different when a different Rasing score threshold for high risk was chosen, a sensitivity analysis was performed, using a risk score of [ 5 as cut-off for patient selection. Analysis was performed using SPSS version 25.0 (IBM Corp, IBM SPSS Statistics for Windows, Armonk, NY). A p value \ 0.05 was considered statistically significant.

RESULTS
Among 7156 surgically treated patients for stage I-III NSCLC between 2015 and 2018, 928 (13.0%) had a high predicted risk (Rasing score [ 4). Of these, 182 (19.6%) underwent an irradical resection. For this study, 290 patients were excluded based on the exclusion criteria. Additionally, a total of 4876 patients treated with CRT for stage I-III NSCLC between 2015 and 2018 were identified. Of these, 2294 were excluded based on the exclusion criteria (Fig. 1).
Before propensity score matching, major statistical differences across most baseline variables were observed between both groups (Table 1). After propensity score matching, 523 pairs of patients remained, and were included for further analysis. SMDs were below the 0.1 threshold, indicating a well-balanced match. On average, patients were 67.5 years old, the majority were male and most had a WHO performance status of 0-1. The majority (60.4%) of patients had a cT4 tumor, and 30.4% and 17.1% were staged cN1 and cN2, respectively. Two thirds of patients had an overall clinical stage IIIA, with 21.4% and 12.5% having stage IIB and IIIB disease, respectively. Of the surgical patients, 114 (21.8%) underwent an R1-2 resection, whereas an R0 resection was achieved in 409 (78.2%). Most surgical patients (49.3%) underwent a lobectomy or bilobectomy, whereas pneumonectomy, sleeve lobectomy, and sublobar resections were performed in 37.3%, 10.7%, and 2.7%, respectively. The majority (71.5%) of surgical patients underwent an open procedure. Adjuvant chemotherapy was administered in 206 (39.4%) of surgical patients. Of the CRT patients, 144 (27.5%) and 379 (72.5%) received sequential and concurrent CRT, respectively.
After a median follow-up of 28.2 months, the median OS among patients treated with surgery was 45.6 months, compared with 27.5 months in the CRT group (log-rank p \ 0.001; Fig. 2). Cox regression analysis revealed an HR of 1.44 (95% CI 1.23-1.70). The CRT group had better OS than those in the surgery group who ended up receiving an R1 or R2 resection (n = 114), with an estimated median OS of the latter group of 20.2 months (p = 0.039; Fig. 3). The HR between the CRT group and the patients who received an R1 or R2 resection was 0.77 (95% CI 0.61-0.99).
The HRs for OS across subgroups are presented in Fig. 4. A subgroup of patients with a histology other than adenocarcinoma or squamous cell carcinoma had a significantly different (interaction p = 0.003) HR for CRT versus surgery, with no significant superiority of surgery nor CRT in that subgroup. Other subgroups, including those based on age, WHO performance status, and clinical staging, did not modify the overall observed effect of CRT versus surgery on OS. Sensitivity analysis using a Rasing score cut-off for high risk of [ 5 (instead of [ 4) yielded similar outcomes regarding the HR of CRT versus surgery on OS (HR 1.36, 95% CI 1.09-1.71; Supplemental Table 1 and Supplemental Fig. 1).

DISCUSSION
To improve patient selection for NSCLC surgery, the Rasing score was developed in 7156 patients and validated in 82,235 patients for individualized pre-treatment prediction of an irradical (R1-2) resection. 8 In the patient group with a high predicted risk of such irradical resection (score [ 4), the current study demonstrates that CRT as an alternative treatment approach resulted in diminished OS compared with surgery. Specifically, the median survival for surgical treatment was found to be 45.6 months, compared with 27.5 months for CRT. However, the subgroup who indeed underwent R1 or R2 resection (21.8%) were found to have a significantly worse survival compared with CRT. These findings support the current standard of surgery as the preferred treatment, even in high-risk cases, provided that an R0 resection is achieved. CRT remains a good alternative for patients in whom an R1 or R2 resection is expected. The Rasing score 8 may aid in estimating this risk, but additional clinical information and judgement (based on parameters not measured in the Rasing score), such as imaging features, comorbidities, assessment by the surgeon, and patient preference remain important for final treatment decision-making.
The current treatment for NSCLC mainly relies on surgery, especially for stage I and II tumors. 9,10 Five-year OS rates for surgery in early stages of the disease have been reported to be 73-78%. 11 For stage III NSCLC a median survival of about 2 years has been achieved. 12,13 The surgical group in the current study had a favorable median OS of 45.6 months (3.8 years). A likely explanation for this superior survival outcome is stricter selection of patients compared with other studies; for example, patients with earlier stage disease, better overall patient condition, and more patients who received adjuvant therapy. Unfortunately, irradical resection remains a considerable impediment, resulting in significantly worse survival, as found in many reports, including the current study. 4,5,[14][15][16][17][18][19] When tumors are deemed irresectable due to local advancement or a patient has contra-indications for surgical treatment, radiation-based therapy can be used as an alternative. 9,10,20 Over the years, this treatment approach has improved from simple thoracic radiotherapy to advanced image-guided radiotherapy and platinum-based chemotherapy to improve survival. Radiation-based therapy when used as an alternative for surgery in clinical stage I patients can achieve a local tumor-control rate of more Preoperatively estimating the chance of an irradical resection may be challenging. In an aim to aid in the clinical decision-making, the Rasing score was developed and validated using large population-based databases. 8 Predictive parameters for an R1-2 resection included histology, clinical T-stage, clinical N-stage, planned extent of resection (e.g., lobectomy, pneumonectomy), and surgical approach (i.e., thoracoscopic, open). Patients with a Rasing score[4 were deemed at high risk and were the population of interest in the current study. When using the Rasing score, it is of crucial importance to realize that no causal   relationship can be inferred from the prediction model (or any other prediction model). For this reason, the risk score should not be inappropriately used to plan the type of surgery. Rather, the model can be appropriately used after staging has been performed and a surgical plan has been made, to then predict the individual risk of an R1-2 resection. The team of physicians can then use the individually calculated risk of an R1-2 resection to optimize or reconsider the proposed treatment approach. Comparison of surgery versus CRT was not previously performed for the specific group with a high predicted risk of irradical resection. This study has compared both treatment methods via propensity score matching of patients identified as high-risk patients in the study by Rasing et al. to patients treated with chemoradiation therapy. The overall median survival found in both groups supports the use of surgery in most cases, with a hazard ratio of 1.44 (95% CI 1.23-1.70). The subgroup analysis did not find any subgroups with a clear benefit of CRT. On the other hand, the overall superiority of surgery was not evident in some subgroups of patients where the confidence intervals around the subgroup HRs contained 1.0. This included those with WHO performance scores 2-3 and patients with cN2 or stage IIIB disease.
The results of this study should be interpreted with consideration of their limitations. First, due to the limited data available, this study could only assess overall survival, and not ascertain other outcomes, such as progression-free survival. Other outcomes, such as burden of treatment or quality of life, were not reported either, and therefore information on the comparison of treatments on those grounds is still lacking. Secondly, as this was not a randomized controlled trial, a risk of selection bias remains. The known confounders have been corrected using propensity score matching, but residual unknown confounders could remain. Potential confounders could include tumor characteristics on imaging, surgeon and/or team experience level, hospital volume, and patient comorbidities. The current study also has considerable strengths, including the relatively large sample size, the population-based design resulting in higher external validity, and the possible generalization of the results to the international context due to the international validation of the Rasing score.
In conclusion, CRT for selected patients with stage IIB-III NSCLC at a high predicted risk for irradical resection using the internationally validated Rasing risk score 8 resulted in worse OS compared with surgery. Consequently, the treatment decision of surgery versus CRT cannot solely be based on that risk score. Since patients who end up receiving an R1 or R2 resection do have inferior outcomes compared with primary CRT, the treatment decision should be based on additional information not covered by the Rasing score, such as imaging features, comorbidities, patient preference, and surgeon's confidence in achieving an R0 resection. ACKNOWLEDGMENT None.
FUNDING No external funding source was involved in this investigation.

FIG. 4
Forest plot comparing hazard ratios of survival after chemoradiation therapy versus surgery in NSCLC patients with a high predicted risk of an irradical (R1-2) resection, organized per subgroup and compared using Cox regression indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons. org/licenses/by/4.0/.