Overall survival according to immunotherapy and radiation treatment for metastatic non-small-cell lung cancer: a National Cancer Database analysis

Background Preclinical studies suggest enhanced anti-tumor activity with combined radioimmunotherapy. We hypothesized that radiation (RT) + immunotherapy would associate with improved overall survival (OS) compared to immunotherapy or chemotherapy alone for patients with newly diagnosed metastatic non-small-cell lung cancer (NSCLC). Methods The National Cancer Database was queried for patients with stage IV NSCLC receiving chemotherapy or immunotherapy from 2013 to 2014. RT modality was classified as stereotactic radiotherapy (SRT) to intra- and/or extracranial sites or non-SRT external beam RT (EBRT). OS was analyzed using the Kaplan-Meier method and Cox proportional hazards models. Results In total, 44,498 patients were included (13% immunotherapy, 46.8% EBRT, and 4.7% SRT). On multivariate analysis, immunotherapy (hazard ratio [HR]:0.81, 95% confidence interval [CI]:0.78–0.83) and SRT (HR:0.78, 95%CI:0.70–0.78) independently associated with improved OS; however, the interaction term for SRT + immunotherapy was insignificant (p = 0.89). For immunotherapy patients, the median OS for no RT, EBRT, and SRT was 14.5, 10.9, and 18.2 months, respectively (p < 0.0001), and EBRT (HR:1.37, 95%CI:1.29–1.46) and SRT (HR:0.78, 95%CI:0.66–0.93) associated with OS on multivariate analysis. In the SRT subset, median OS for immunotherapy and chemotherapy was 18.2 and 14.3 months, respectively (p = 0.004), with immunotherapy (HR:0.82, 95%CI:0.69–0.98) associating with OS on multivariate analysis. Furthermore, for patients receiving SRT, biologically effective dose (BED) > 60 Gy was independently associated with improved OS (HR:0.79, 95%CI:0.70–0.90, p < 0.0001) on multivariate analysis with a significant interaction between BED and systemic treatment (p = 0.008). Conclusions Treatment with SRT associated with improved OS for patients with metastatic NSCLC irrespective of systemic treatment. The high survival for patients receiving SRT + immunotherapy strongly argues for evaluation in randomized trials.


Background
In the United States, lung cancer is the most common cause of cancer-related mortality with 234,000 new cases and 154,000 deaths from this disease expected in 2018 [1]. Of these new cases, approximately 85% will be non-small-cell lung cancer (NSCLC) [2] which is most commonly diagnosed after distant metastasis has already occurred [3]. Unfortunately, the 5-year overall survival (OS) for patients with stage IV NSCLC from 2007 to 2013 was just 5.2% [3]. Therefore, it is not surprising that there has been significant interest in the application of new treatment modalities and therapeutic combinations to improve prognosis for this group of patients, including the use of immunotherapy [2] and consolidative radiotherapy (RT) for select oligometastatic cases [4].
Immune checkpoint blockade targeting the PD-1/ PD-L1 axis was approved for metastatic NSCLC based on the groundbreaking CheckMate 017 and CheckMate 057 trials [5,6]. Since that time, the indications for immune checkpoint blockade in NSCLC have expanded to include first-line therapy for patients with metastatic disease and tumor PD-L1 expression ≥50% [7]. Additionally, an OS benefit for all PD-L1 categories receiving pembrolizumab plus cytotoxic chemotherapy compared to cytotoxic chemotherapy alone has been observed in the setting of metastatic NSCLC without EGFR or ALK mutations [8,9]. Moreover, patients with any history of predominantly intracranial (61%) RT receiving pembrolizumab for advanced NSCLC had improved OS in a secondary analysis of the phase I KEYNOTE-001 trial [10]. These data combined with promising preclinical observations [11,12] and the theoretical advantages of combined radioimmunotherapy [13] have promoted enthusiasm for concurrent or sequential RT plus immunotherapy in advanced NSCLC. Therefore, we hypothesized that RT plus immunotherapy would associate with improved OS for patients with stage IV NSCLC in the National Cancer Database (NCDB) and further examined the influence of RT technique and dose on OS.

Data source
Data were procured from Health Insurance Portability and Accountability Act-compliant, de-identified participant user files extracted from the NCDB in an institutional-review board exempt study. The NCDB is a joint project of the Commission on Cancer (CoC) of the American College of Surgeons and the American Cancer Society. It was established in 1989 and includes demographic and oncologic outcomes data from over 1500 CoC-accredited cancer programs. Notably, the NCDB captures over 70% of incident cancer diagnoses nationally with more than 34 million historical records to date [14,15]. While the data used in this study are derived from a de-identified NCDB file, the American College of Surgeons and the CoC have not verified and are not responsible for the analytic or statistical methodology employed, or the conclusions drawn from these data by the investigators.

Study population
The NCDB was queried for patients with histologically proven stage IV NSCLC (American Joint Committee on Cancer staging manual, 7th Edition) who had received chemotherapy or immunotherapy as a first course of treatment at a CoC-accredited center from 2013 to 2014. This timeframe was chosen to ensure consistency in immunotherapy coding as some agents previously coded as chemotherapy were reclassified as immunotherapy in the NCDB effective January 1, 2013 [16]. Additionally, restricting the cohort to those receiving systemic therapy prior to the end of 2014 allowed for sufficient follow-up for correlation with OS.
In total, 44,498 patients met these inclusion criteria with n = 38,691 (86.9%) receiving chemotherapy and n = 5807 (13.1%) receiving immunotherapy. Patients were then further classified as receiving no RT, stereotactic radiotherapy (SRT), or non-SRT external beam radiotherapy (EBRT). RT modality was defined using NCDB codes whereby SRT included any patient specifically coded as receiving "stereotactic radiosurgery, NOS," "linac radiosurgery," or "Gamma Knife." Thus, patients defined as receiving SRT necessarily included individuals receiving RT to intraand/or extracranial sites. All other patients receiving RT and not labeled with at least 1 of these SRT-defining codes were considered to have received EBRT regardless of prescription dose or RT fractionation.

Statistical considerations
Statistical analyses were performed with SPSS version 24 (IBM Inc., Armonk, NY) and JMP statistical software version 13.0 (SAS Institute, Cary, NC). The chi-square test was used for comparisons of patient-, tumor-, and treatment-related factors that were grouped categorically for analytic purposes. Variables that were specifically compared between chemotherapy and immunotherapy patient subsets included sex, age, modified Charlson-Deyo comorbidity index [17,18], tumor histology, race (white vs. non-white), insurance status, treating facility type, and RT modality. Cumulative RT dose in Gray (Gy) and the number of RT fractions were used to calculate biologically effective dose (BED) using the formula BED = n · d · (1 + d/ 10) where n = number of fractions, d = dose per fraction, and 10 is the assumed alpha/beta ratio for NSCLC tumors. BED was compared between groups using univariate logistic regression.
Kaplan-Meier methods were used for univariate survival analyses while the log-rank test was used to make OS comparisons between and/or among groups based on treatment modality. Multivariate survival analyses with resultant hazard ratios (HRs) and 95% confidence intervals (CIs) were performed using Cox proportional hazards regression. Covariates for multivariate analyses were identified a priori and included factors known to influence OS for metastatic NSCLC as well as factors significantly differing between chemotherapy and immunotherapy groups. A test for interaction between SRT and immunotherapy was performed by comparing HRs for the SRT and no RT groups as a function of systemic therapy using a multiplicative interaction model. A similar test for interaction between BED and systemic therapy was performed for patients receiving SRT with available dose/fractionation data. All statistical tests used two-sided hypothesis testing with a type I error < 0.05 defining statistical significance.

Population characteristics
Patient-, tumor-, and treatment-related characteristics for the entire cohort and for the immunotherapy vs. chemotherapy subpopulations are displayed in Table 1. A total of 44,498 patients with stage IV NSCLC met inclusion criteria with n = 5807 (13.1%) receiving immunotherapy and n = 38,691 (86.9%) receiving chemotherapy. Compared to the chemotherapy group, the immunotherapy group contained a larger proportion of female patients, younger patients, and patients with adenocarcinoma histology. There was no significant difference in the proportion of patients receiving treatment at academic or non-academic centers as a function of systemic therapy (p = 0.48).
Regarding RT, n = 21,593 (48.5%) did not receive RT, n = 20,821 (46.8%) received EBRT (34.6% intracranial), and n = 2048 (4.7%) received SRT (81.3% intracranial). Of these, n = 20,031 (87.6%) had evaluable RT dose/fractionation information for BED calculations. The median dose for patients receiving EBRT was 30 Gy in a median of 10 fractions (interquartile range [IQR]: 10-15) while the median dose for patients receiving SRT was 22 Gy in a median of a single fraction (IQR: 1-4). The median BED for the SRT group was 60 Gy (IQR: 43.2-70.4) and was significantly higher than the median BED of 39 Gy (IQR: 39.0-50.7) for the EBRT group (p < 0.0001). Of patients in the SRT group with known dose/fractionation information, 34.5% received a BED > 60 Gy with this subset prescribed a median RT dose of 24 Gy in a median of a single fraction. The most common RT dose/ fractionation prescriptions for patients receiving SRT to a BED ≤60 Gy vs. > 60 Gy are displayed in Table 2.

Survival
The median survival of patients receiving immunotherapy was 17.3 months vs. 14.4 months in the chemotherapy group (p = 0.03). Kaplan-Meier plots investigating the relationship between combined systemic and RT modalities for the immunotherapy (n = 5807) and SRT (n = 2084) subpopulations are displayed in panels A and B of Fig. 1 Fig. 3 and also found a significant association with OS (p < 0.0001) with patients receiving extracranial SRT and immunotherapy having the best prognosis. The results of a multivariate analysis for OS for the entire cohort are displayed in Table 3. The independent predictors with the largest impact on OS included Charlson-Deyo comorbidity score (HR:1.23, 95% CI:1.19-1.28 for scores ≥2 vs. 0, p < 0.0001), systemic therapy (HR: 0.81, 95% CI: 0.78-0.84 for immunotherapy vs. chemotherapy, p < 0.0001), and RT modality (HR: 0.78, 95% CI: 0.69-0.78 for SRT vs. no RT; HR: 1.16, 95% CI: 1.13-1.18 for EBRT vs. no RT, p < 0.0001). This same multivariate analysis performed while adjusting for intra-vs. extracranial RT continued to identify systemic therapy (HR: 0.87, 95% CI: 0.83-0.91 for immunotherapy vs. chemotherapy, p < 0.0001) and RT modality (HR: 0.60, 95% CI: 0.57-0.63 for SRT vs. EBRT, p < 0.0001) as independently associated with OS. Since RT site was also significantly associated with OS (HR: 1.12, 95% CI: 1.08-1.15 for intracranial vs. extracranial RT, p < 0.0001), results of the multivariate analysis were separately reported for the intra-and extracranial RT subgroups as displayed in Table 4. RT modality as well as systemic therapy continued to be independently associated with OS for both the intracranial and extracranial subgroups. Moreover, the association of combined immunotherapy and SRT with improved OS was persistent on multivariate analyses for the immunotherapy and SRT subgroups as respectively displayed in Tables 5 and 6. Specifically, for patients receiving immunotherapy, EBRT (HR: 1.37, 95% CI: 1.29-1.46) was associated with significantly reduced OS while SRT (HR: 0.78, 95% CI: 0.66-0.93) was associated with significantly improved OS (p < 0.0001 comparing no RT, EBRT, and SRT). Furthermore, for patients undergoing SRT, immunotherapy (HR: 0.82, 95% CI: 0.69-0.98, p = 0.031) associated with higher OS compared to chemotherapy. Despite the prolonged OS for patients receiving combined immunotherapy and SRT compared to patients receiving alternative combinations of systemic and RT modalities, the interaction term for SRT and immunotherapy did not reach significance (p = 0.89).
The influence of BED ("low BED" ≤60 Gy vs. "high BED" > 60 Gy) [19] and systemic therapy on OS for patients receiving SRT with available dose/fractionation information (n = 1739) is displayed in panel B of Fig. 2. Notably, patients receiving high-BED SRT with immunotherapy had a median OS of 32.1 months compared to 16.6 months for high-BED SRT with chemotherapy, 15.3 months for low-BED SRT with immunotherapy, and 13.7 months for low-BED SRT with chemotherapy (p < 0.0001 when comparing all groups). A similar Kaplan-Meier plot showing OS using the same stratification for patients receiving specifically extracranial SRT is displayed in panel B of Fig. 3 and once again demonstrated the highest OS for patients receiving high-BED SRT and immunotherapy (p < 0.0001). On multivariate analysis including all SRT patients with available dose/fractionation data as displayed in Table 7, BED had the largest impact on OS of all covariates (HR:0.79, 95% CI:0.70-0.90, p < 0.0001), and systemic therapy was no longer significantly associated with OS after including BED in the  presented phase III KEYNOTE-042 trial [20]. This similarity in survival outcomes adds validity to these results and demonstrates the efficacy of immunotherapy in patients treated outside of a randomized phase III study.
On subset analysis, patients undergoing predominantly intracranial SRT had improved OS compared to those who did not regardless of systemic therapy choice. Interestingly, among patients receiving immunotherapy, SRT associated with improved OS while non-SRT EBRT associated with These data add to a growing body of literature attempting to clarify the potential role of SRT in the setting of metastatic NSCLC [4,21,22]. Specifically, our finding that patients receiving SRT experienced improved OS should be considered alongside the previously reported high progression-free survival (PFS) of 14.7 months and median OS of 20.4 months in a single-arm, phase II protocol of SRT with concurrent erlotinib for patients with stage IV NSCLC involving ≤6 sites after progression through first-or subsequent-line chemotherapy [22]. Moreover, compared to an expected median OS of just 11 months for stage IV NSCLC treated with a first-line platinum doublet [23], the median OS of 14.3 months for SRT + chemotherapy and 18.2 months for SRT + immunotherapy as observed in our report are quite promising. Although analyses investigating differences in OS based on the existence of oligometastatic versus polymetastatic disease are not possible using the NCDB, Gomez et al. have previously reported a high median PFS of 11.9 months after local consolidative therapy including SRT for patients with oligometastatic NSCLC involving ≤3 sites that had not progressed after first-line chemotherapy [4]. This again suggests that select patients with metastatic NSCLC may be appropriate candidates for an aggressive treatment paradigm incorporating local ablative therapy.
Preclinical studies have suggested enhanced anti-tumor activity for combined radioimmunotherapy including intracranial RT in mouse glioblastoma models [11,12]. This observation in the laboratory was bolstered by a secondary analysis of the KEYNOTE-001 phase I trial wherein patients with advanced NSCLC receiving any prior RT (61% intracranial) followed by pembrolizumab had a significantly higher median OS of 10.7 months vs. 5.3 months for those receiving pembrolizumab without previous RT [10]. These data are encouraging when viewed as a foundation to support the use of RT as a strategy to induce response to immune checkpoint blockade for patients with metastatic NSCLC as the unselected overall response rate with single agent PD-1/PD-L1 blockade is only 17-19% [24,25]. This possibility is further suggested by the results of a randomized phase II study of SRT and sequential pembrolizumab vs. pembrolizumab alone for patients with advanced NSCLC with the experimental arm having an overall response rate of 41% and a median PFS of 6.4 months compared to 19% and 1.8 months, respectively, for the control arm [26]. Our observed significantly prolonged median OS of 18.2 months for combined SRT and immunotherapy in the NCDB is       consistent with these previously published data; however, evidence of a synergistic interaction between SRT and immunotherapy did not reach significance. RT dose and fractionation may influence the complex relationship between RT-induced immunomodulation and anti-tumor response [13,19]. Standard RT fractionation ranges from 1.5-2.2 Gy/fraction; however, SRT delivering substantially higher doses in a small (e.g., ≤5) number of fractions and "hypofractionated" RT typically delivering 40-50 Gy over 5-10 fractions have been increasingly used for aggressive local therapy for select cases of metastatic NSCLC [4,22,26]. This shift towards higher RT doses per fraction for select patients with stage IV NSCLC has been supported by mouse models wherein higher RT doses in fewer fractions were associated with increased tumor CD8+ T-cell infiltration [27,28]. Although these preclinical data are somewhat inconsistent with the results of Shaverdian et al.'s secondary analysis of KEYNOTE-001 showing improved oncologic outcomes for predominantly non-SRT radiation (n = 49/70, 70%) followed by pembrolizumab, they offer a potential mechanistic explanation for several observations in our data including (1) reduced OS for EBRT +immunotherapy, (2) an impressive median OS of 32.1 months for high-BED SRT + immunotherapy, and (3) a significant interaction for BED and systemic therapy among patients receiving SRT. Of note, it is possible that the relative influence of dose and fractionation on the biological interactions between RT and immunotherapy may differ as a function of treatment setting since conventionally fractionated EBRT followed by durvalumab has been associated with improved OS for patients with locally advanced rather than metastatic NSCLC [29].
Areas for future investigation include elucidation of RT's immunomodulatory effects on the tumor microenvironment, determination of clinically optimal radioimmunotherapy combinations/sequencing for metastatic NSCLC, and investigation of the toxicity of combined radioimmunotherapy. To this end, inclusion of molecular correlative studies on future prospective trials of combined radioimmunotherapy for metastatic NSCLC is imperative. A growing number of early-phase protocols [30][31][32][33][34][35][36] including diverse radioimmunotherapy strategies will inevitably work towards this goal and lay the foundations for eventual standard-of-care defining phase III trials. Although reports have suggested that combined thoracic radioimmunotherapy is generally safe [37], continued prospective evaluation is recommended to determine the influence of RT dose/sequencing and patient-related factors on the incidence of grade 3+ toxicity with this treatment approach.
Limitations to this study include the lack of information on biomarkers including PD-L1 status, possibility for coding errors in NCDB data, the inability to examine alternative endpoints of interest including toxicity or PFS, heterogeneity of agents coded as immunotherapy in the NCDB, and selection bias regarding treatment received. It is likely that patient selection factors influenced some associations between aggressive local therapy and improved OS as patients with a lower disease burden including oligometastases may have preferentially received SRT and/or higher-BED SRT regimens. Furthermore, the inclusion of patients receiving heterogeneous RT doses within the EBRT subset may limit generalizations for this group. Despite these caveats, our data demonstrate that SRT associates with improved OS for patients with stage IV NSCLC in the NCDB regardless of systemic treatment.

Conclusions
Overall, the hypothesis-generating observation of superior OS for patients receiving high-BED SRT and immunotherapy, the reduction in OS for patients receiving EBRT and immunotherapy, and the significant interaction between BED and systemic therapy suggest that RT technique may influence the efficacy of immunotherapy in this setting. These findings strongly support efforts to evaluate optimal multimodality radiation and immunotherapy strategies in prospective randomized trials.

Funding
This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.

Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Authors' contributions CCF and MK were responsible for formal data analysis and interpretation. MK was responsible for conceptualization, data curation, and study supervision. CCF was responsible for writing the original draft of the manuscript. All authors were responsible for review and editing of the manuscript and approval of the final version for submission.

Ethics approval and consent to participate
This study was institutional-review board exempt and performed with de-identified participant user files from the National Cancer Database.

Consent for publication
Individual data is not contained in this manuscript.
Competing interests RRW provides consulting and/or advising services for Aettis, Inc., AstraZeneca, Genus, ImmunoVir LLC, Merck Serono S.A., Nano Proteagen, Reflexion Pharmaceuticals, RiMO, and Shuttle Pharmaceuticals. RRW participated in guest speaking with sponsored industry travel for Boehringer Ingelheim Ltd. RRW also has equity and/or intellectual property rights with Boost