Skip to main content

Abstract

Aim: To investigate the factors associated with treatment selection and overall survival for first-line EGFR-tyrosine kinase inhibitors (EGFR-TKIs) therapy among patients with non-small-cell lung cancer. Materials & methods: We conducted a retrospective cohort study of linked administrative health databases in Ontario, Canada. Results: A total of 1011 patients received an EGFR-TKI as first-line therapy. Treatment selection and overall survival associated with these treatments were affected by age, sex, geographical residency, comorbidities and different sites of metastasis. Conclusion: Though recent approval of osimertinib offers a potential new standard of care in the first-line setting, earlier generation TKIs remain pillars in treatment of non-small-cell lung cancer therapeutic armamentarium. Our findings may contribute to optimizing treatment sequencing of EGFR-TKIs to maximize clinical benefits.
Lung cancer is a leading cause of cancer-related death worldwide [1]. In 2018, estimates suggested a total of 1.76 million persons died of lung cancer, accounting for 18.4% of all cancer-related deaths [2]. Approximately 85% of all lung cancer cases are non-small-cell lung cancer (NSCLC) [3]. The 5-year survival rate for lung cancer ranges between 13 and 20% [4,5]. Most patients are diagnosed in advanced stages of NSCLC, which results in poor prognosis [6].
Advancements in our understanding of cancer biology have allowed us to tailor treatment approaches based on patients’ genetic profiles. A half of NSCLC cases are associated with known mutations, and several actionable gene alterations have been identified for targeted treatment [7–9]. EGFR has been one of the most prevailing targets for devising specific treatment algorithms for patients with NSCLC. Approximately 15% of NSCLC cases have an activating mutation in the EGFR genes in exon 18–21 [10]. Those harboring EGFR mutations are eligible to receive EGFR-tyrosine kinase inhibitors (EGFR-TKIs), which have been demonstrated to improve objective response rates and progression-free survival compared with conventional chemotherapy in first-line settings [11–23].
There are currently five EGFR-TKIs approved as first-line treatment for NSCLC in Canada, including the first-generation TKIs erlotinib (Tarceva; Hoffman-La Roche, Basel, Switzerland) and gefitinib (Iressa; AstraZeneca, London, UK), second-generation TKIs afatinib (Giotrif; Boehringer Ingelheim, Ingelheim, Germany) and dacomitinib (Vizimpro; Pfizer, NY, USA), and third-generation TKI osimertinib (Tagrisso; AstraZeneca, London, UK).
Currently, a dearth of clinical evidence exists to suggest whether one EGFR-TKI should be chosen over another in a typical first-line setting [24,25]. This suggests multiple factors could affect treatment selection, though a delineation of these factors in the case of EGFR-TKIs has never been undertaken. Previous research has identified general factors related to prescribing decisions, for example, treatment sequencing, evidence from clinical trials, safety/toxicity profiles associated with each agent, growing familiarity with new agents among practitioners, regional/institutional preference, reimbursement and influence of pharmaceutical companies [26,27], though the applicability of these factors to EGFR-TKIs is unknown.
While several randomized controlled trials and observational studies have estimated the efficacy/effectiveness of prognostic factor-guided EGFR-TKIs in advanced/metastatic EGFR mutation-positive NSCLC population [11–23,28–35], limited information is available on the longitudinal effects of EGFR-TKIs at the population level. Furthermore, evidence regarding the comparative effectiveness of these EGFR-TKIs is inconsistent across the literature.
To our knowledge, no studies thus far have investigated how prescribing decisions for EGFR-TKIs are made and the factors that may affect these prescribing decisions. Using population-based administrative health datasets, we sought to determine which factors influence the receipt of certain EGFR-TKIs in first-line settings and investigate how these are associated with overall survival (OS). Due to the recency of the approval dates and concomitant lack of data for dacomitinib and osimertinib, the present study only pertained to afatinib, erlotinib and gefitinib.

Materials & methods

Study design

This was a retrospective, population-based cohort study of linked health administrative data in Ontario, Canada. The datasets are housed at the Institute for Clinical Evaluative Sciences (ICES), a prescribed entity under Ontario’s Personal Health Information Protection Act. The Act authorizes ICES to draw individual patient-level data from multiple health administrative datasets for researchers to use in secondary analyses. The research was cleared for ethics by the Office of Research Ethics at the University of Waterloo (ORE # 41067).

Study population

The study included all patients diagnosed with NSCLC in Ontario between 01 January 2010 and 31 August 2019. NSCLC cases were identified through the Ontario Cancer Registry (OCR). The International Classification of Diseases for Oncology, Third Edition (ICD-O-3) site codes 34.0–34.9, in combination with relevant histology codes for nonsquamous, squamous and not otherwise specified were used to identify cases of primary lung cancer from the OCR. Inclusion criteria were age ≥18 years at diagnosis, locally advanced or metastatic NSCLC, and receipt of afatinib, erlotinib or gefitinib as first-line treatment. We excluded persons with death dates on or before the date of NSCLC diagnosis, and individuals who received more than one EGFR-TKI as first-line treatment. The dataset did not contain any information on biomarker status; thus, we assumed patients with records of EGFR-TKI prescription in the first-line setting had positive EGFR mutation status.

Data sources

Multiple health administrative datasets were linked using encrypted unique identifiers. The OCR contains information on incident cancer cases and patients who have died of cancer in Ontario since 1964 [36,37]. The OCR includes data on date and stage of NSCLC diagnosis, age, sex, geographical location, rural versus urban residence and date of death. The Registered Persons Database contains demographic information and vital statistics on all residents of Ontario who are eligible for universal healthcare coverage in the province. The Canadian Institute for Health Information – Discharge Abstract Database contains data on diagnoses and procedures for all hospitalizations in Ontario. The Ontario Drug Benefits (ODB) database contains data on all prescription medications dispensed to those eligible for publicly funded drug coverage. These include all persons aged ≥65 years, persons living in homes for special care and long-term care homes, persons receiving professional services through the home and community care service programs, persons receiving social assistance, and persons receiving benefits through Trillium Drug Program, a scheme which helps people with high prescription drug costs relative to their net household income. The ODB does not capture information covered by private insurance and compassionate supplies from manufacturers. The Activity Level Reporting system contains information on all systemic and radiation therapy services and outpatient oncology clinic visits.

Covariates

We searched the literature and consulted expert opinion to identify several sociodemographic and clinical factors that may influence treatment selection and OS [24,38–44]. These factors included: year of diagnosis, age, sex, rurality, neighborhood income quintile, local health integration network (LHIN), clinical stage, histology, Charlson comorbidity index and sites of metastasis (bone, brain, liver and lung).
Neighborhood household income was determined through linkage of postal codes to Canadian census data and stratified into three tertiles, with the first and last tertiles representing neighborhoods with the lowest and highest income status, respectively. Charlson comorbidity index was determined from hospitalization data utilizing a 2-year ‘look-back’ window, with the score from the most recently available hospitalization record applied to each participant [40,41]. We followed Stavem et al.’s approach and considered missing comorbidities to be absent [45]. At the time of data collection, publicly funded healthcare services in Ontario were administered on a regional basis by 14 LHINs, each with its own distinct geographical territory. Recently, Ontario integrated these LHINs into five regions consisting of North, West, Toronto Central, East and South regions. The analyses were conducted reflective of these changes.

Statistical analysis

All variables were categorical and described using frequencies and percentages. To explore the factors associated with treatment selection, we conducted two separate logistic regression analyses comparing afatinib to gefitinib and erlotinib to gefitinib. We used gefitinib as the reference category because it was the most established treatment group among the three EGFR-TKIs; erlotinib and gefitinib have been in use since 2010, while afatinib was approved for use and publicly funded in 2014. Furthermore, while erlotinib is only publicly funded for second- and third-line settings in Ontario, it is indicated for first-line setting as well. The aim of our models was not to predict treatment selection, but to identify which variables may be of importance to clinicians prescribing EGFR-TKIs for first-line treatment. Therefore, our focus was not to identify the most parsimonious model, but rather to build an explanatory model and examine the effects of all relevant covariates on treatment selection.
A priori, we defined sociodemographic and clinical factors that may be important in clinical decision making for treatment selection and its associated outcomes (see ‘Covariates’ above) [24,38–44]. A series of chi-square tests were conducted to test the associations between independent variables and the outcome variable. In addition, all explanatory variables were checked for multicollinearity using variance inflation factors. Discrimination of the models was assessed with the area under the receiver operating characteristic curve (AUC). Calibration of the models was evaluated using the Hosmer–Lemeshow goodness-of-fit tests.
OS was assessed using the Kaplan–Meier method on the overall population and various patient subgroups. The OS was calculated from the date of diagnosis of NSCLC to death (for any reason) or the last day of patient follow-up (censored). Comparisons between groups were performed using the log-rank test. Multivariable Cox proportional hazards models were used to determine adjusted hazard ratios (aHRs) and to evaluate the predictive factors for survival. Statistical significance was set at α = 0.05. All analyses were performed using SAS version 9.4 (SAS Institute, Inc., NC, USA).

Results

A total of 88,179 patients were identified as having a primary diagnosis of lung cancer in the OCR between 2010 and 2019. Of these, 68,334 were NSCLC cases and 33,321 had stage IIIB/IV NSCLC. In total, 1011 patients met the eligibility criteria and were enrolled in the study (Table 1); 67 patients were excluded as they had records of receiving more than one EGFR-TKI in the first-line setting (Figure 1). A total of 110 (10.9%) patients received afatinib, while 482 (47.7%) and 419 (41.4%) received erlotinib and gefitinib, respectively. Male patients constituted 41.8% of the study population. Almost all patients had nonsquamous histology (98.7%) and no patients had squamous cell carcinoma. The majority of patients were at stage IV NSCLC (89.8%) at the time of diagnosis.
Table 1. Baseline patient sociodemographic and clinical characteristics.
VariablesAfatinibErlotinibGefitinibTotalp-value
 n = 110 (%)n = 482 (%)n = 419 (%)n = 1011 (%) 
Year of diagnosis<0.0001
– 2010–20146 (5.4%)425 (88.2%)178 (42.5%)609 (60.2%) 
– 2015–2019104 (94.6%)57 (11.8%)241 (57.5%)402 (39.8%) 
Age, years<0.0001
– 18–5936 (32.7%)126 (26.1%)84 (20.1%)246 (24.3%) 
– 60–6915 (13.6%)85 (17.6%)53 (12.7%)153 (15.1%) 
– 70–7955 (50.0%)247 (51.3%)217 (51.8%)519 (51.3%) 
– 80+4 (3.6%)24 (5.0%)65 (15.5%)93 (9.2%) 
Sex<0.0001
– Males44 (40.0%)245 (50.8%)134 (32.0%)423 (41.8%) 
– Females66 (60.0%)237 (49.2%)285 (68.0%)588 (58.2%) 
Rurality0.20
– Rural15 (13.6%)65 (13.5%)41 (9.8%)121 (12.0%) 
– Urban95 (86.4%)417 (86.5%)378 (90.2%)890 (88.0%) 
Neighborhood income quintile0.27
– 1 (poorest)57 (51.8%)208 (43.2%)167 (39.9%)432 (42.7%) 
– 216 (14.6%)86 (17.8%)78 (18.6%)180 (17.8%) 
– 3 (wealthiest)37 (33.6%)188 (39.0%)174 (41.5%)399 (39.5%) 
LHIN<0.0001
– North11 (10.3%)30 (6.3%)17 (4.1%)58 (5.8%) 
– West27 (25.2%)150 (31.4%)84 (20.4%)261 (26.2%) 
– Toronto9 (8.4%)33 (6.9%)55 (13.3%)97 (9.7%) 
– Central34 (31.8%)89 (18.6%)168 (40.8%)291 (29.2%) 
– East26 (24.3%)176 (36.8%)88 (21.4%)290 (29.1%) 
Clinical stage0.06
– IIIB11 (10.0%)60 (12.4%)32 (7.6%)103 (10.2%) 
– IV99 (90.0%)422 (87.6%)387 (92.4%)908 (89.8%) 
Histology0.33
– Nonsquamous109 (99.1%)478 (99.2%)411 (98.1%)998 (98.7%) 
– Squamous cell0 (0%)0 (0%)0 (0%)0 (0%) 
– NOS1 (0.9%)4 (0.8%)8 (1.9%)13 (1.3%) 
CCI0.25
– No107 (97.3%)452 (93.8%)400 (95.5%)959 (94.9%) 
– Yes3 (2.7%)30 (6.2%)19 (4.5%)52 (5.1%) 
Sites of metastasis 
– Liver16 (14.6%)66 (13.7%)44 (10.5%)126 (12.5%)0.27
– Bone40 (36.4%)150 (31.1%)167 (39.9%)357 (35.3%)0.02
– Brain26 (23.6%)66 (13.7%)98 (23.4%)190 (18.8%)0.0004
– Lung19 (17.3%)121 (25.1%)101 (24.1%)241 (23.8%)0.22
CCI: Charlson comorbidity index; LHIN: Local health integration network; NOS: Not otherwise specified.
Figure 1. Flow diagram of patient selection.
NSCLC: Non-small-cell lung cancer; TKI: Tyrosine kinase inhibitors.

Treatment selection

The results of the logistic regression analyses are presented in Table 2. We found no evidence of multicollinearity between explanatory variables. The goodness-of-fit of the models were confirmed (Hosmer–Lemeshow: chi-square = 7.29 and p = 0.50 for gefitinib vs erlotinib; chi-square = 8.33 and p = 0.40 for gefitinib vs afatinib) and the models exhibited moderate discriminatory capacity (AUC = 0.75 for gefitinib vs erlotinib; AUC = 0.69 for gefitinib vs afatinib).
Table 2. Adjusted odds ratios for prescription of afatinib and erlotinib compared with gefitinib.
VariableAfatinib, OR (95% CI)
(n = 110)
Erlotinib, OR (95% CI)
(n = 482)
Gefitinib (reference)
(n = 419)
Age, years
– 18–59111
– 60–690.83 (0.40–1.73)1.12 (0.69–1.83)1
– 70–790.66 (0.39–1.12)0.73 (0.50–1.06)1
– 80+0.14 (0.04–0.42)0.19 (0.10–0.34)1
Sex
– Males1.43 (0.89–2.27)2.59 (1.90–3.52)1
– Females111
Rurality
– Rural1.10 (0.52–2.31)0.98 (0.60–1.59)1
– Urban111
Neighborhood income quintile
– 1 (poorest)1.53 (0.78–2.97)1.07 (0.71–1.63)1
– 2111
– 3 (wealthiest)0.94 (0.47–1.87)0.98 (0.64–1.48)1
LHIN
– North3.30 (1.02–10.64)2.57 (1.15–5.75)1
– West2.03 (0.86–4.82)2.94 (1.68–5.14)1
– Toronto111
– Central1.03 (0.45–2.34)0.73 (0.42–1.26)1
– East1.66 (0.70–3.95)3.51 (2.01–6.11)1
Clinical stage
– IIIB1.70 (0.75–3.85)1.43 (0.84–2.44)1
– IV111
CCI
– No1.58 (0.44–5.68)0.81 (0.42–1.59)1
– Yes111
Bone metastasis
– Yes0.78 (0.47–1.29)0.58 (0.42–0.80)1
– No/unknown111
Liver metastasis
– Yes1.58 (0.78–3.20)1.74 (1.10–2.78)1
– No/unknown111
Brain metastasis
– Yes1.14 (0.67–1.97)0.53 (0.36–0.78)1
– No/unknown111
Lung metastasis
– Yes0.64 (0.36–1.16)1.37 (0.96–1.95)1
– No/unknown111
CCI: Charlson comorbidity index; LHIN: Local health integration network; OR: Odds ratio.
Age was associated with prescribing choice of EGFR-TKIs, with older patients more likely to be prescribed gefitinib over afatinib and erlotinib. Compared with patients aged 18–59 years, the adjusted odds ratio (aOR) for prescribing afatinib in lieu of gefitinib in patients aged ≥80 years was 0.14 (95% CI: 0.04–0.42), and 0.19 (95% CI: 0.10–0.34) for erlotinib versus gefitinib. A larger proportion of patients aged ≥70 years received gefitinib over the other two drugs. Erlotinib was more commonly prescribed for male patients than gefitinib (aOR: 2.59; 95% CI: 1.90–3.52).
Regional prescribing preferences were evident. Patients residing in LHIN – North region, compared with Toronto Central, were less likely to be prescribed gefitinib than afatinib and erlotinib. The adjusted odds of receiving afatinib was 3.30-times greater than receiving gefitinib (95% CI: 1.02–10.64), while the adjusted odds was 2.57-times greater for erlotinib versus gefitinib (95% CI: 1.15–5.75). In the West and East regions, erlotinib was more commonly prescribed over afatinib and gefitinib. The adjusted odds of being prescribed erlotinib in the West region was 2.94 (95% CI: 1.68–5.14) and 3.51 (95% CI: 2.01–6.11) in the East region.
We found associations between sites of metastasis and prescribing decisions. Patients with metastasis to bone and brain were less likely to be prescribed erlotinib than gefitinib. Among patients with bone metastasis, the aOR for erlotinib prescription was 0.58 (95% CI: 0.42–0.80) and 0.53 (95% CI: 0.36–0.78) for patients with brain metastasis. Erlotinib was more commonly prescribed for patients with liver metastasis than gefitinib (aOR: 1.74; 95% CI: 1.10–2.78).

Survival analysis

The median OS of the overall cohort was 19.53 months (95% CI: 18.38–20.75; Figure 2). Statistically significant differences in OS were observed across the EGFR-TKIs; the median OS were 31.04 (95% CI: 23.41–42.05), 17.36 (95% CI: 16.04–18.48) and 21.63 months (95% CI: 19.27–23.18) for afatinib, erlotinib and gefitinib, respectively (Figure 3A). Significantly shorter OS was observed in patients who were male (median: 17.33 months; 95% CI: 16.01–18.97; p for log-rank test = 0.0017; Figure 3C), had presence of comorbidities (15.81 months; 95% CI: 13.32–20.52; p for log-rank test = 0.026; Figure 3H), and had metastasis to liver (16.27 months; 95% CI: 14.14–18.44; p for log-rank test = 0.0001; Figure 3I) and bone (17.98 months; 95% CI: 16.21–20.38; p for log-rank test = 0.0094; Figure 3J). Furthermore, shorter median OS was observed in older patients, but the difference was not statistically significant.
Figure 2. Overall survival for the entire study population.
Figure 3. Overall survival by patient factors.
Overall survival according to treatment (afatinib, erlotinib or gefitinib) (A); age groups (B); sex (C); income (D); geographical residency (E); rural versus urban (F); clinical stage (G); presence of comorbidity (H); liver metastasis (I); bone metastasis (J); brain metastasis (K); and lung metastasis (L).
A multivariable Cox regression model showed that prescription of erlotinib (aHR: 1.58; 95% CI: 1.33–1.86; p < .0001), aged ≥80 years (aHR: 1.42; 95% CI: 1.07–1.88; p = 0.02), presence of comorbidities (aHR: 1.37; 95% CI: 1.01–1.86; p = 0.04), liver metastasis (aHR: 1.49; 95% CI: 1.22–1.83; p = 0.0001), bone metastasis (aHR: 1.31; 95% CI: 1.13–1.53; p = 0.0004) and brain metastasis (aHR: 1.30; 95% CI: 1.07–1.57; p = 0.007) were inversely associated with OS (Table 3). Prescription of afatinib (aHR: 0.61; 95% CI: 0.45–0.82; p = 0.0013) was positively associated with OS (Table 3).
Table 3. Multivariable Cox regression of overall survival.
VariableHR (95% CI)p-value
EGFR-TKI
– Afatinib0.61 (0.45–0.82)0.0013
– Erlotinib1.58 (1.33–1.86)<0.0001
– GefitinibRef 
Age
– 18–59Ref 
– 60–690.86 (0.68–1.08)0.18
– 70–791.15 (0.96–1.37)0.13
– 80+1.42 (1.07–1.88)0.02
Sex
– Males1.14 (0.99–1.32)0.08
– FemalesRef 
Rurality
– UrbanRef 
– Rural1.05 (0.84–1.30)0.67
Clinical stage
– IIIB1.05 (0.83–1.34)0.68
– IVRef 
CCI
– NoRef 
– Yes1.37 (1.01–1.86)0.04
Income
– PoorestRef 
– Middle1.09 (0.90–1.32)0.40
– Wealthiest1.04 (0.89–1.21)0.67
LHIN
– North1.35 (0.92–1.97)0.12
– West1.13 (0.86–1.48)0.40
– TorontoRef 
– Central1.05 (0.81–1.37)0.72
– East1.15 (0.88–1.50)0.30
Liver metastasis
– NoRef 
– Yes1.49 (1.22–1.83)0.0001
Bone metastasis
– NoRef 
– Yes1.31 (1.13–1.53)0.0004
Brain metastasis
– NoRef 
– Yes1.30 (1.07–1.57)0.007
Lung metastasis
– NoRef 
– Yes0.94 (0.80–1.11)0.48
CCI: Charlson comorbidity index; HR: Hazards ratio; LHIN: Local health integration network; Ref: Reference; TKI: Tyrosine kinase inhibitor.

Discussion

We identified sociodemographic and clinical factors influencing treatment selection and OS among patients who received EGFR-TKIs between 2010 and 2019. Compared with gefitinib, erlotinib was prescribed more frequently for those who were males, residing in certain geographical locations (North, West and East regions), and had liver metastasis. The results for afatinib were similar; while a higher prescription of afatinib was noted among patients with no comorbidities and had brain metastasis, the results were not statistically significant. Compared with afatinib and erlotinib, gefitinib was more commonly prescribed for older patients and those with bone metastasis. We found type of EGFR-TKI therapy, age, presence of comorbidities, and metastasis to liver, bone and brain as potential independent prognostic factors for OS.
We expected to find some of the associations reported in Tables 2 and 3. Previous studies showed that brain metastasis was an independent prognostic factor for progression-free survival and OS [46,47]. Among patients with EGFR mutation, the incidence of brain metastasis can be as large as one in two persons [48]. The use of afatinib is widely accepted for patients with brain metastasis due to improved capacity for blood–brain barrier penetration into the intracranial lesions to initiate antitumor activities compared with reversible EGFR-TKIs [49,50].
In terms of the use of gefitinib for patients with bone metastasis, previous studies suggested that gefitinib may reduce bone metastasis growth through inhibition of EGF signaling pathways in bone stromal cells and improve pathologic fractures [51,52]. Furthermore, the high usage of gefitinib in older patients may partly be explained by the safety/toxicity profile of afatinib over erlotinib and gefitinib. While results from the noninterventional RealGiDo study indicated that adverse events (AEs) with afatinib can be managed with dose adjustments and care measures [53], older patients may not be able to handle the intensity of AEs associated with afatinib. For patients with terminal NSCLC, the intent of treatment would likely focus on health-related quality of life, which involves minimizing treatment-related AEs and managing symptoms.
Another factor to consider is acquired resistance, the most common being the development of a T790M mutation, which occurs in 50–70% of cases [54–56]. Studies have shown that afatinib is able to overcome acquired first-generation EGFR-TKI resistance [57,58]; therefore, afatinib may have been reserved for subsequent use in later lines for patients who fail previous EGFR-TKI therapy. Although afatinib is not publicly funded for second-line settings in Ontario, funding of afatinib for patients who have initiated another EGFR-TKI therapy in the first-line setting, and who have not had disease progression, are considered on a case-by-case basis.
Our median OS estimate stratified by EGFR-TKIs are inconsistent with what was reported in several trials [11–23]. The median OS observed in our study for erlotinib and gefitinib were shorter than what was observed in most published trials. Exception was the results reported in the IPASS trial, where median OS was 18.8 months for gefitinib [11]. However, the median OS of afatinib observed in our study was longer than what was reported in Phase IIb Lux-Lung 7 trial (27.9 months), Phase III Lux-Lung 3 (28.2 months) and Lux-Lung 6 (23.1 months) trials [20,59]. The differences in median OS between EGFR-TKIs were consistent with other observational studies that reported effectiveness of EGFR-TKIs in clinical practice [28–35]. Most studies reported higher median OS for afatinib over erlotinib and gefitinib [29,31,33]. However, Li and colleagues reported higher median OS for patients who received erlotinib (23.2 months) than afatinib (20.7 months) [32]. In addition, Chao and colleagues reported substantially higher median OS for patients who received erlotinib (34.6 months) than gefitinib (19.2 months) [28].
Given the fact that 15% of patients with nonsquamous histology harbor EGFR mutation, we had expected a larger sample for our study. However, the relatively small sample size could be attributed to the initial challenges in implementation of biomarker testing in Ontario in the early 2010s, along with logistical difficulties, for example, delayed turnaround times, which led chemotherapy to be used as the first-line treatment to avoid clinical deterioration [60,61]. A previous study has suggested approximately one in four patients do not undergo biomarker testing [62]. Furthermore, the ODB database captured only information related to publicly-funded medications; therefore, prescription medications covered by private insurance and compassionate supplies from manufacturers could not be considered.
An interesting finding of this study was that only 38.9% of our study cohort initiated a second-line therapy, which questions the notion of reserving therapies for subsequent use (e.g., development of acquired T790M mutation resistance) to maximize the duration of chemotherapy-free treatment. The results from the Phase III FLAURA trial demonstrated superior efficacy and safety profiles associated with osimertinib compared with first-generation TKIs in the first-line setting, regardless of T790M mutational status [63]. In terms of treatment sequencing, patients receiving osimertinib as first-line treatment would not receive any subsequent EGFR-TKIs upon progression and would likely involve treatments with platinum doublets or immune checkpoint inhibitors. Therefore, clinical challenges remain in deciding whether the most effective therapy should be used as first-line treatment or be reserved for later lines to expand treatment options.
The strength of our study was the linkage and use of population-based administrative datasets, which captured all relevant data and complete follow-up for all patients. To our knowledge, this was the first study to systematically assess the factors influencing prescribing decisions associated with EGFR-TKIs using administrative datasets. Nonetheless, our study has limitations. First, the number of patients who received afatinib was relatively small compared with patients who received erlotinib and gefitinib, due to later approval of afatinib. The difference in sample size resulting from the late licensing date may have contributed to selection and length-time bias in our study, which may have caused overestimation of odds ratios and survival estimates in our analyses. The OS data for persons who received erlotinib and gefitinib were more mature compared with persons who received afatinib. In addition, previous studies have demonstrated the importance of Eastern Cooperative Oncology Group performance status and the type of EGFR mutation status, for example, exon 19 deletion (Exon19DEL) and the exon 21 codon 858 point mutation (L858R), as important factors to consider in treatment selection and survival [20,64–70]. However, a lack of data on Eastern Cooperative Oncology Group performance status and EGFR mutation status prohibited us from carrying out any analyses involving these factors. We assumed that patients who received any of the three EGFR-TKIs had positive EGFR mutation, regardless of the type. Although highly unlikely, we cannot rule out the possibility that TKIs may have been prescribed to EGFR wild-type patients. Lastly, since the ODB database did not capture information on private insurance claims and compassionate supplies, we were not able to assess prescribing differences based on types of insurance.

Conclusion

The results of the present study demonstrated that factors including age, sex, geographical residency, and metastasis to bone, liver and brain were independent factors influencing treatment selection of EGFR-TKIs. Presences of comorbidities, in addition to the aforementioned factors, were independent prognostic factors for OS. In clinical practice, there were significant differences in OS between the three EGFR-TKIs. Additional population-based studies are required to compare the clinical effectiveness of EGFR-TKIs stratified by EGFR-mutation status.
Summary Points
Factors affecting treatment selection and overall survival (OS) of first-line EGFR-tyrosine kinase inhibitors remain unclear.
We conducted a retrospective, population-based cohort study linking non-small-cell lung cancer cases identified in the Ontario Cancer Registry between 2010 and 2019 to other administrative health databases.
The outcomes of interest for this study were prescription of one of three EGFR-tyrosine kinase inhibitors (afatinib, erlotinib and gefitinib) and OS.
Logistic regression models were conducted to evaluate the associations between treatment selection and patient characteristics.
The Kaplan–Meier method was used to estimate OS; multivariable Cox regression analyses were performed to test the associations between OS and patient characteristics.
A total of 1011 patients met the inclusion/exclusion criteria.
The median OS were 31.04 (95% CI: 23.41–42.05), 17.36 (95% CI: 16.04–18.48), and 21.63 months (95% CI: 19.27–23.18) for afatinib, erlotinib and gefitinib, respectively.
Age, sex, geographical residency and sites of metastasis affected treatment selection.
Prescription of erlotinib or gefitinib, older age, presence of comorbidities, and metastasis to bone, liver and brain were poor prognostic factors associated with OS.

Author contributions

Y-J Kim, M Oremus, HH Chen and S Horton were responsible for study conception and design. Y-J Kim, M Oremus and D Fearon were responsible for data analysis. Y-J Kim, M Oremus, HH Chen, T McFarlane, D Fearon and S Horton revised and approved final manuscript for publication.

Acknowledgments

The opinions, results and conclusions reported in this paper are those of the authors and are independent from the funding sources. No endorsement by ICES or the Ontario MOHLTC is intended or should be inferred. Parts of this material are based on data and/or information compiled and provided by the Canadian Institute for Health Information (CIHI). However, the analyses, conclusions, opinions and statements expressed herein are those of the author(s), and not necessarily those of CIHI. We thank B Zhang at ICES for her assistance with data linkage.

Financial & competing interests disclosure

This study was supported by ICES, which is funded by an annual grant from the Ontario Ministry of Health and Long-Term Care (MOHLTC). The authors have no other relevant affiliations or financial involvement with any organization or entity with a financial interest in or financial conflict with the subject matter or materials discussed in the manuscript apart from those disclosed.
No writing assistance was utilized in the production of this manuscript.

Data sharing statement

The dataset from this study is held securely in coded form at ICES. While data sharing agreements prohibit ICES from making the data set publicly available, access may be granted to those who meet prespecified criteria for confidential access, available at www.ices.on.ca/DAS. The full dataset creation plan and underlying analytic code are available from the authors upon request, understanding that the computer programs may rely upon coding templates or macros that are unique to ICES and are therefore either inaccessible or may require modification.

References

Papers of special note have been highlighted as: • of interest; •• of considerable interest
1.
Siegel RL, Miller KD, Jemal A. Cancer statistics, 2019. CA Cancer J. Clin. 69(1), 7–34 (2019).
2.
Bray F, Ferlay J, Soerjomataram I et al. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J. Clin. 68(6), 394–424 (2018).
3.
Howlader N, Noone AM, Krapcho M et al. SEER Cancer Statistics Review, 1975–2016. National Cancer Institute, MD, USA (2019). https://seer.cancer.gov/csr/1975_2016/
4.
Allemani C, Matsuda T, Di Carlo V et al. Global surveillance of trends in cancer survival 2000–14 (CONCORD-3): analysis of individual records for 37 513 025 patients diagnosed with one of 18 cancers from 322 population-based registries in 71 countries. Lancet 391(10125), 1023–1075 (2018).
5.
Barta JA, Powell CA, Wisnivesky JP. Global epidemiology of lung cancer. Ann. Glob. Health 85(1), 8 (2019).
6.
Walters S, Maringe C, Coleman MP et al. Lung cancer survival and stage at diagnosis in Australia, Canada, Denmark, Norway, Sweden and the UK: a population-based study, 2004–2007. Thorax 68(6), 551–564 (2013).
7.
Dearden S, Stevens J, Wu Y et al. Mutation incidence and coincidence in non-small-cell lung cancer: meta-analyses by ethnicity and histology (mutMAP). Ann. Oncol. 24(9), 2371–2376 (2013).
8.
Greulich H. The genomics of lung adenocarcinoma. Genes Cancer 1(12), 1200–1210 (2010).
9.
Korpanty G, Graham D, Vincent M et al. Biomarkers that currently affect clinical practice in lung cancer: EGFR, ALK, MET, ROS-1, and KRAS. Front. Oncol. 4, 204 (2014).
• Reviews treatment options specified by biomarkers relevant to non-small-cell lung cancer – EGFR, ALK, Met, ROS-1 and KRAS.
10.
Graham RP, Treece AL, Lindeman NI et al. Worldwide frequency of commonly detected EGFR mutations. Arch. Pathol. Lab. Med. 142(2), 163–167 (2018).
11.
Fukuoka M, Wu YL, Thongprasert S et al. Biomarker analyses and final overall survival results from a Phase III, randomized, open-label, first-line study of gefitinib versus carboplatin/paclitaxel in clinically selected patients with advanced non-small-cell lung cancer in Asia (IPASS). J. Clin. Oncol. 29(21), 2866–2874 (2011).
12.
Han JY, Park K, Kim SW et al. First-SIGNAL: first-line single-agent iressa versus gemcitabine and cisplatin trial in never-smokers with adenocarcinoma of the lung. J. Clin. Oncol. 30(10), 1122–1128 (2012).
13.
Inoue A, Kobayashi K, Maemondo M et al. Updated overall survival results from a randomized Phase III trial comparing gefitinib with carboplatin-paclitaxel for chemo-naïve non-small cell lung cancer with sensitive EGFR gene mutations (NEJ002). Ann. Oncol. 24(1), 54–59 (2013).
14.
Maemondo M, Inoue A, Kobayashi K et al. Gefitinib or chemotherapy for non-small-cell lung cancer with mutated EGFR. N. Engl. J. Med. 362(25), 2380–2388 (2010).
15.
Mitsudomi T, Morita S, Yatabe Y et al. Gefitinib versus cisplatin plus docetaxel in patients with non-small-cell lung cancer harbouring mutations of the epidermal growth factor receptor (WJTOG3405): an open label, randomised Phase III trial. Lancet Oncol. 11(2), 121–128 (2010).
16.
Rosell R, Carcereny E, Gervais R et al. Erlotinib versus standard chemotherapy as first-line treatment for European patients with advanced EGFR mutation-positive non-small-cell lung cancer (EURTAC): a multicentre, open-label, randomised Phase III trial. Lancet Oncol. 13(3), 239–246 (2012).
17.
Sequist LV, Yang JCH, Yamamoto N et al. Phase III study of afatinib or cisplatin plus pemetrexed in patients with metastatic lung adenocarcinoma with EGFR mutations. J. Clin. Oncol. 31(27), 3327–3334 (2013).
18.
Wu YL, Zhou C, Hu CP et al. Afatinib versus cisplatin plus gemcitabine for first-line treatment of Asian patients with advanced non-small-cell lung cancer harbouring EGFR mutations (LUX-Lung 6): an open-label, randomised Phase III trial. Lancet Oncol. 15(2), 213–222 (2014).
19.
Wu YL, Zhou C, Liam CK et al. First-line erlotinib versus gemcitabine/cisplatin in patients with advanced EGFR mutation-positive non-small-cell lung cancer: analyses from the Phase III, randomized, open-label, ENSURE study. Ann. Oncol. 26(9), 1883–1889 (2015).
20.
Yang JCH, Wu YL, Schuler M et al. Afatinib versus cisplatin-based chemotherapy for EGFR mutation-positive lung adenocarcinoma (LUX-Lung 3 and LUX-Lung 6): analysis of overall survival data from two randomised, Phase III trials. Lancet Oncol. 16(2), 141–151 (2015).
21.
Yoshioka H, Shimokawa M, Seto T et al. Final overall survival results of WJTOG3405, a randomized Phase III trial comparing gefitinib versus cisplatin with docetaxel as the first-line treatment for patients with stage IIIB/IV or postoperative recurrent EGFR mutation-positive non-small-cell lung cancer. Ann. Oncol. 30(12), 1978–1984 (2019).
22.
Zhou C, Wu YL, Chen G et al. Erlotinib versus chemotherapy as first-line treatment for patients with advanced EGFR mutation-positive non-small-cell lung cancer (OPTIMAL, CTONG-0802): a multicentre, open-label, randomised, Phase III study. Lancet Oncol. 12(8), 735–742 (2011).
23.
Zhou C, Wu YL, Chen G et al. Final overall survival results from a randomised, Phase III study of erlotinib versus chemotherapy as first-line treatment of EGFR mutation-positive advanced non-small-cell lung cancer (OPTIMAL, CTONG-0802). Ann. Oncol. 26(9), 1877–1883 (2015).
24.
Girard N. Optimizing outcomes and treatment sequences in EGFR mutation-positive non-small-cell lung cancer: recent updates. Future Oncol. 15(25), 2983–2997 (2019).
•• Reviews important clinical factors associated with optimization of treatment sequencing of EGFR-tyrosine kinase inhibitors.
25.
Nan X, Xie C, Yu X et al. EGFR TKI as first-line treatment for patients with advanced EGFR mutation-positive non-small-cell lung cancer. Oncotarget 8(43), 75712–75726 (2017).
26.
Fleischman W, Agrawal S, King M et al. Association between payments from manufacturers of pharmaceuticals to physicians and regional prescribing: cross sectional ecological study. BMJ 354, i4189 (2016).
27.
Schumock GT, Walton SM, Park HY et al. Factors that influence prescribing decisions. Ann. Pharmacother. 38(4), 557–562 (2004).
28.
Chao PI, Cheng G, Zhang L et al. Retrospective analysis of gefitinib and erlotinib in EGFR-mutated non-small-cell lung cancer patients. J. Lung Health Dis. 1(1), 16–24 (2017).
29.
Fujiwara A, Yoshida M, Fujimoto H et al. A retrospective comparison of the clinical efficacy of gefitinib, erlotinib, and afatinib in Japanese patients with non-small-cell lung cancer. Oncol. Res. 26(7), 1031–1036 (2018).
30.
Ho GF, Chai CS, Alip A et al. Real-world experience of first-line afatinib in patients with EGFR-mutant advanced NSCLC: a multicenter observational study. BMC Cancer 19(1), 896 (2019).
31.
Ito K, Murotani K, Kubo A et al. Comparative analysis of overall survival using propensity score between first- and second-generation EGFR-TKI: real world data of 1354 patients with EGFR mutant NSCLC. Ann. Oncol. 29(Suppl. 8), viii493–viii547 (2018).
32.
Li Y, Appius A, Pattipaka T et al. Real-world management of patients with epidermal growth factor receptor (EGFR) mutation-positive non-small-cell lung cancer in the USA. PLoS ONE 14(1), e0209709 (2019).
33.
Lin YT, Chen JS, Liao WY et al. Clinical outcomes and secondary epidermal growth factor receptor (EGFR) T790M mutation among first-line gefitinib, erlotinib, and afatinib-treated non-small-cell lung cancer patients with activating EGFR mutations. Int. J. Cancer 144(11), 2887–2896 (2019).
•• Investigated independent factors for secondary T790M mutation.
34.
Perol M, Morere JF, Fabre E et al. French real-life efficacy of 1st line gefitinib in EGFR mutation-positive NSCLC in the prospective EPIDAURE study: results by EGFR exon 19 Del and L85R mutation subtypes. Ann. Oncol. 27(Suppl. 6), vi429 (2016).
35.
Tokaca N, O'Brien M, Bhosle J et al. Real-world outcomes with first-line afatinib in EGFR mutant NSCLC adenocarcinoma: a single centre experience exploring effects of dose-reduction. J. Thorac. Oncol. 13(Suppl. 4), S1–S139 (2018).
36.
Clarke EA, Marrett LD, Kreiger N. Cancer registration in Ontario: a computer approach. IARC Sci. Publ. 95, 246–257 (1991).
37.
Robles SC, Marrett LD, Clarke EA et al. An application of capture-recapture methods to the estimation of completeness of cancer registration. J. Clin. Epidemiol. 41(5), 495–501 (1988).
38.
Bergqvist M, Christensen HN, Wiklund F et al. Real world utilization of EGFR TKIs and prognostic factors for survival in NSCLC during 2010–2016 in Sweden: a nationwide observational study. Int. J. Cancer 146(9), 2510–2517 (2020).
• This population-based study investigated treatment pattern and survival using Swedish nationwide cancer registry.
39.
Booth CM, Li G, Zhang-Salomons J et al. The impact of socioeconomic status on stage of cancer at diagnosis and survival: a population-based study in Ontario, Canada. Cancer 116(17), 4160–4167 (2010).
• Investigated the impact of patient socioeconomic status on cancer survival.
40.
Charlson ME, Pompei P, Ales KL et al. A new method of classifying prognostic comorbidity in longitudinal studies: development and validation. J. Chronic Dis. 40(5), 373–383 (1987).
41.
Deyo RA, Cherkin DC, Ciol MA. Adapting a clinical comorbidity index for use with ICD-9-CM administrative databases. J. Clin. Epidemiol. 45(6), 613–619 (1992).
42.
Kim SY, Myung JK, Kim HR et al. Factors that predict clinical benefit of EGFR TKI therapy in patients with EGFR wild-type lung adenocarcinoma. Tuberc. Respir. Dis. (Seoul) 82(1), 62–70 (2019).
43.
Lin JH, Lin D, Xu L et al. The association between clinical prognostic factors and epidermal growth factor receptor-tyrosine kinase inhibitor (EGFR-TKI) efficacy in advanced non-small-cell lung cancer patients: a retrospective assessment of 94 cases with EGFR mutations. Oncotarget 8(2), 3412–3421 (2017).
44.
Schrijvers CT, Coebergh JW, Mackenbach JP. Socioeconomic status and comorbidity among newly diagnosed cancer patients. Cancer 80(8), 1482–1488 (1997).
45.
Stavem K, Hoel H, Skjaker SA et al. Charlson comorbidity index derived from chart review or administrative data: agreement and prediction of mortality in intensive care patients. Clin. Epidemiol. 9, 311–320 (2017).
46.
Hendriks LEL, Smit EF, Vosse BAH et al. EGFR mutated non-small-cell lung cancer patients: more prone to development of bone and brain metastases? Lung Cancer 84(1), 86–91 (2014).
47.
Sperduto PW, Yang TJ, Beal K et al. Estimating survival in patients with lung cancer and brain metastases: an update of the Graded Prognostic Assessment for lung cancer using molecular markers (Lung-molGPA). JAMA Oncol. 3(6), 827–831 (2017).
48.
Sorensen JB, Hansen HH, Hansen M et al. Brain metastases in adenocarcinoma of the lung: frequency, risk groups, and prognosis. J. Clin. Oncol. 6(9), 1474–1480 (1988).
49.
Hochmair M. Medical treatment options for patients with epidermal growth factor receptor mutation-positive non-small cell lung cancer suffering from brain metastases and/or leptomeningeal disease. Target. Oncol. 13(3), 269–285 (2018).
•• Discusses EGFR-tyrosine kinase inhibitor treatment options for patients with brain metastases.
50.
Su PL, Wu YL, Chang WY et al. Preventing and treating brain metastases with three first-line EGFR-tyrosine kinase inhibitors in patients with EGFR mutation-positive advanced non-small-cell lung cancer. Ther. Adv. Med. Oncol. 10, 1758835918797589 (2018).
51.
Lu X, Wang Q, Hu G et al. ADAMTS1 and MMP1 proteolytically engage EGF-like ligands in an osteolytic signaling cascade for bone metastasis. Genes Dev. 23(16), 1882–1894 (2009).
52.
Okano Y, Nishio M. Efficacy of gefitinib in treatment of lung cancer patients with bone metastasis. Clin. Calcium 18(4), 527–533 (2008).
53.
Halmos B, Tan EH, Soo RA et al. Impact of afatinib dose modification on safety and effectiveness in patients with EGFR mutation-positive advanced NSCLC: results from a global real-world study (RealGiDo). Lung Cancer 127, 103–111 (2019).
54.
Arcila ME, Oxnard GR, Nafa K et al. Rebiopsy of lung cancer patients with acquired resistance to EGFR inhibitors and enhanced detection of the T790M mutation using a locked nucleic acid-based assay. Clin. Cancer Res. 17(5), 1169–1180 (2011).
55.
Sequist LV, Waltman BA, Dias-Santagata D et al. Genotypic and histological evolution of lung cancers acquiring resistance to EGFR inhibitors. Sci. Transl. Med. 3(75), 75ra26 (2011).
56.
Yang JCH, Ahn MJ, Kim DW et al. Osimertinib in pretreated T790M-positive advanced non-small-cell lung cancer: AURA study Phase II extension component. J. Clin. Oncol. 35(12), 1288–1296 (2017).
57.
Heigener DF, Reck M. Mutations in the epidermal growth factor receptor gene in non-small cell lung cancer: impact on treatment beyond gefitinib and erlotinib. Adv. Ther. 28(2), 126–133 (2011).
58.
Heigener DF, Schumann C, Sebastian M et al. Afatinib in non-small-cell lung cancer harboring uncommon EGFR mutations pretreated with reversible EGFR inhibitors. Oncologist 20(10), 1167–1174 (2015).
59.
Park K, Tan EH, O'Byrne K et al. Afatinib versus gefitinib as first-line treatment of patients with EGFR mutation-positive non-small-cell lung cancer (LUX-Lung 7): a Phase IIB, open-label, randomised controlled trial. Lancet Oncol. 17(5), 577–589 (2016).
60.
Cheema PK, Raphael S, El-Maraghi R et al. Rate of EGFR mutation testing for patients with nonsquamous non-small-cell lung cancer with implementation of reflex testing by pathologists. Curr. Oncol. 24(1), 16–22 (2017).
61.
Ellis PM, Verma S, Sehdev S et al. Challenges to implementation of an epidermal growth factor receptor testing strategy for non-small-cell lung cancer in a publicly funded health care system. J. Thorac. Oncol. 8(9), 1136–1141 (2013).
62.
Spicer J, Tischer B, Peters M. EGFR mutation testing and oncologist treatment choice in advanced NSCLC: global trends and differences. Ann. Oncol. 26(Suppl. 1), i57 (2015).
63.
Soria JC, Ohe Y, Vansteenkiste J et al. Osimertinib in untreated EGFR-mutated advanced non-small-cell lung cancer. N. Engl. J. Med. 378(2), 113–125 (2018).
64.
Cha YK, Lee HY, Ahn MJ et al. Survival outcome assessed according to tumor burden and progression patterns in patients with epidermal growth factor receptor mutant lung adenocarcinoma undergoing epidermal growth factor receptor tyrosine kinase inhibitor therapy. Clin. Lung Cancer 16(3), 228–236 (2015).
65.
Hung MS, Fang YH, Lin YC et al. Survival-associated factors of first-line EGFR-tyrosine kinase inhibitor responders and non-responders in lung adenocarcinoma patients with common EGFR mutations. Mol. Clin. Oncol. 8(3), 421–428 (2018).
66.
Jackman DM, Yeap BY, Sequist LV et al. Exon 19 deletion mutations of epidermal growth factor receptor are associated with prolonged survival in non-small cell lung cancer patients treated with gefitinib or erlotinib. Clin. Cancer Res. 12(13), 3908–3914 (2006).
67.
Pirker R, Pereira JR, Szczesna A et al. Prognostic factors in patients with advanced non-small cell lung cancer: data from the Phase III FLEX study. Lung Cancer 77(2), 376–382 (2012).
68.
Ren S, Su C, Wang Z et al. Epithelial phenotype as a predictive marker for response to EGFR-TKIs in non-small cell lung cancer patients with wild-type EGFR. Int. J. Cancer 135(12), 2962–2971 (2014).
69.
Riely GJ, Pao W, Pham D et al. Clinical course of patients with non-small-cell lung cancer and epidermal growth factor receptor exon 19 and exon 21 mutations treated with gefitinib or erlotinib. Clin. Cancer Res. 12(3 Pt 1), 839–844 (2006).
70.
Tu HY, Ke EE, Yang JJ et al. A comprehensive review of uncommon EGFR mutations in patients with non-small-cell lung cancer. Lung Cancer 114, 96–102 (2017).