Alterations in Genes of the EGFR Signaling Pathway and Their Relationship to EGFR Tyrosine Kinase Inhibitor Sensitivity in Lung Cancer Cell Lines

doi:10.1371/journal.pone.0004576

PLoS ONE. 2009; 4(2): e4576.

Published online 2009 February 24. doi: 10.1371/journal.pone.0004576

PMCID: PMC2642732

Copyright Gandhi et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Alterations in Genes of the EGFR Signaling Pathway and Their Relationship to EGFR Tyrosine Kinase Inhibitor Sensitivity in Lung Cancer Cell Lines

Jeet Gandhi,^#¹ Jianling Zhang,^#¹ Yang Xie,⁵ Junichi Soh,¹ Hisayuki Shigematsu,¹ Wei Zhang,¹ Hiromasa Yamamoto,¹ Michael Peyton,¹ Luc Girard,¹ William W. Lockwood,⁶ Wan L. Lam,⁶ Marileila Varella-Garcia,⁷ John D. Minna,^1,^3,⁴ and Adi F. Gazdar^1,²^*

¹Hamon Center for Therapeutic Oncology Research, University of Texas Southwestern Medical Center at Dallas, Dallas, Texas, United States of America

²Department of Pathology, University of Texas Southwestern Medical Center at Dallas, Dallas, Texas, United States of America

³Department of Internal Medicine, University of Texas Southwestern Medical Center at Dallas, Dallas, Texas, United States of America

⁴Department of Pharmacology, University of Texas Southwestern Medical Center at Dallas, Dallas, Texas, United States of America

⁵Department of Clinical Sciences, University of Texas Southwestern Medical Center at Dallas, Dallas, Texas, United States of America

⁶Department of Cancer Genetics and Developmental Biology, British Columbia Cancer Research Centre, Vancouver, British Columbia, Canada

⁷Department of Internal Medicine, University of Colorado Cancer Center, Aurora, Colorado, United States of America

Alfred Lewin, Editor

University of Florida, United States of America

^#Contributed equally.

* E-mail: adi.gazdar/at/utsouthwestern.edu

Conceived and designed the experiments: JKG JZ JS HS HY MP LG WWL WLL MG JDM AFG. Performed the experiments: JKG JZ JS HS WZ HY MP WWL. Analyzed the data: JKG JZ YX JS LG WWL MG AFG. Contributed reagents/materials/analysis tools: WLL MG JDM. Wrote the paper: JKG JZ YX JS WZ HY WWL WLL MG JDM AFG.

Received September 16, 2008; Accepted December 18, 2008.

This article has been cited by other articles in PMC.

Abstract

Background

Deregulation of EGFR signaling is common in non-small cell lung cancers (NSCLC) and this finding led to the development of tyrosine kinase inhibitors (TKIs) that are highly effective in a subset of NSCLC. Mutations of EGFR (mEGFR) and copy number gains (CNGs) of EGFR (gEGFR) and HER2 (gHER2) have been reported to predict for TKI response. Mutations in KRAS (mKRAS) are associated with primary resistance to TKIs.

Methodology/Principal Findings

We investigated the relationship between mutations, CNGs and response to TKIs in a large panel of NSCLC cell lines. Genes studied were EGFR, HER2, HER3 HER4, KRAS, BRAF and PIK3CA. Mutations were detected by sequencing, while CNGs were determined by quantitative PCR (qPCR), fluorescence in situ hybridization (FISH) and array comparative genomic hybridization (aCGH). IC50 values for the TKIs gefitinib (Iressa) and erlotinib (Tarceva) were determined by MTS assay. For any of the seven genes tested, mutations (39/77, 50.6%), copy number gains (50/77, 64.9%) or either (65/77, 84.4%) were frequent in NSCLC lines. Mutations of EGFR (13%) and KRAS (24.7%) were frequent, while they were less frequent for the other genes. The three techniques for determining CNG were well correlated, and qPCR data were used for further analyses. CNGs were relatively frequent for EGFR and KRAS in adenocarcinomas. While mutations were largely mutually exclusive, CNGs were not. EGFR and KRAS mutant lines frequently demonstrated mutant allele specific imbalance i.e. the mutant form was usually in great excess compared to the wild type form. On a molar basis, sensitivity to gefitinib and erlotinib were highly correlated. Multivariate analyses led to the following results:

1. mEGFR and gEGFR and gHER2 were independent factors related to gefitinib sensitivity, in descending order of importance.

2. mKRAS was associated with increased in vitro resistance to gefitinib.

Conclusions/Significance

Our in vitro studies confirm and extend clinical observations and demonstrate the relative importance of both EGFR mutations and CNGs and HER2 CNGs in the sensitivity to TKIs.

Introduction

Lung cancer is the leading cause of all cancer deaths worldwide [1]. Despite the recent advances in diagnosis and multimodality therapies for lung cancers, the prognosis still remains poor with 5-year survival rates of only 16% for all stages [2].

Lung cancer is characterized by the accumulation of multiple genetic and epigenetic alterations including somatic mutations and gene copy number gains or both which results in the activation of oncogenes or inactivation of tumor suppressor genes [3]. Epidermal growth factor receptor (EGFR) deregulation has been observed in multiple tumor types including non-small cell lung cancers (NSCLCs) [4]. Hirsch et. al. identified frequent EGFR protein over expression (62%) in NSCLCs of squamous cell and adenocarcinoma subtypes [5]. EGFR over expression is often associated with adverse prognosis [6]. The receptor tyrosine kinase (RTK) super-family of cell surface receptors serves as mediators of cell signaling by extra-cellular growth factors. Members of the ERBB family of RTKs including EGFR(HER1/ERBB1), HER2 (ERBB2/EGFR2), HER3 (ERBB3/EGFR3) and HER4 (ERBB4/EGFR4) have received much attention given their strong association with malignant proliferation.

The RAS/MAPK and PI3K/AKT pathways are major signaling networks linking EGFR activation to cell proliferation and survival [7]. As discussed below, EGFR signaling pathway genes have been reported to be mutated in NSCLC. Depending on the geographical location, EGFR and KRAS mutations have been identified in ~10% −30% of NSCLCs [3], [8]. EGFR mutations are independently associated with adenocarcinoma histology, East Asian ethnicity, never smoking status and female gender. Mutations of KRAS also target adenocarcinoma histology, but otherwise differ from EGFR mutations because they are relatively rare in East Asians and occur more frequently in males and smokers [9]. Less commonly, somatic mutations have also been found in other EGFR pathway genes including HER2 (~2%) [10], HER4 (~2%) [11], BRAF (~2%) [12], and PIK3CA (~4%) [13], [14].

Gene copy number gains (CNGs) due to focal amplification or chromosomal polysomy, is one of the other major mechanisms of oncogene activation [15]. Lockwood et al. identified multiple components of the EGFR pathway signaling were frequently amplified and over-expressed in NSCLC. Interestingly, they also found that EGFR pathway gene amplification was more frequent in the adenocarcinoma subtype of NSCLC.

Because of the frequent deregulation of EGFR pathway genes in NSCLC, EGFR became one of the first rationally selected molecules for targeted therapy. While initial targeted approaches utilized monoclonal antibodies which block the ligand-receptor interaction, newer approaches have utilized small molecule reversible tyrosine kinase inhibitors (TKIs). The tyrosine kinase activity of EGFR is required for the biochemical responses induced by this receptor [7], [16], [17]. Two TKIs, gefitinib (Iressa, AstraZeneca) and erlotinib (Tarceva, Genentech) have been widely used in the treatment of advanced or recurrent NSCLC. Responses were noted in subsets, notably East Asian ethnicity, female gender, never smoking status and adenocarcinoma histology [18], [19], [20]. Subsequently, EGFR mutations were identified in the tyrosine kinase domain, and predicted for response to TKI in the same subset of patients [21], [22], [23]. According to a meta-analysis of 1170 patients, more than 70% of NSCLCs with EGFR mutations responded to TKIs, whereas 10% of tumors without EGFR mutation responded [24]. However, further studies indicated that factors other than EGFR mutations may play a role in determining response and survival after TKI therapy. EGFR gene copy number gain was associated with significantly improved TKI sensitivity and survival in a large unselected study with appropriate controls [25]. In addition, other members of the EGFR family i.e. HER2 [26] and EGFR3 [27] may be important factors involved in TKI sensitivity. A further complexity is the clinical observation that somatic mutations of KRAS confer intrinsic resistance to TKIs [28].

To further understand the relationship between TKI sensitivity and deregulation of EGFR pathway genes, we examined the mutation and copy number status of seven of these genes (EGFR, HER2, HER3, HER4, KRAS, BRAF, and PIK3CA) in a large panel of lung cancer cell lines and correlated the data with in vitro sensitivity to TKIs.

Materials and Methods

Tumor Cell Lines

We studied a total of 112 cell lines consisting of 77 NSCLC and 32 small cell lung cancer (SCLC) and three lines from small cell cancers at extrapulmonary sites (extrapulmonary small cell cancers, ExPuSC) [29]. All except three of these cell lines were established by the authors [30] at one of two locations. Cell lines initiated at the NCI have the prefix NCI-H while lines established at UT Southwestern have the prefix HCC. NCI-H3255 was obtained from Dr. Bruce Johnson (Lowe Center for Thoracic Oncology, Dana-Farber Cancer Institute, Boston, MA) [22]. PC-9 (originally from the Tokyo Medical University, Tokyo, Japan) was obtained from Dr. Bert Vogelstein (Johns Hopkins University School of Medicine, Baltimore, MD). Calu-3 was purchased from American Type Culture Center (Manassas, VA). We also investigated eight immortalized human bronchial epithelial cell lines (HBECs, HBEC1KT, HBEC3KT, HBEC4KT, HBEC5KT, HBEC17KT, HBEC30KT, HBEC31KT and HBEC34KT), which were initiated by us [31].

Most of the tumor cell lines were maintained in RPMI1640 supplemented with 5–10% fetal bovine serum (FBS). A few cell lines were grown in ACL4 (for NSCLC lines) and HITES (for SCLC lines) supplemented with 5% FBS. All HBEC cell lines were maintained in Keratinocyte-SFM (Invitrogen, Carlsbad, CA) with bovine pituitary extract (BPE) and recombinant epidermal growth factor (EGF) [31]. All cell lines were incubated at 37°C in a humidified atmosphere with 5% CO₂.

The genetic fingerprint of each cell line was obtained (Powerplex 1.2 system, Promega, Madison, WI) and each cell line had a unique genetic profile which was identical to the profiles obtained from the American Type Culture Collection.

DNA and RNA extraction

Genomic DNA was obtained from cell line pellets by standard phenol-chloroform (1 [ratio]

1) extraction followed by ethanol precipitation [32] or by using DNeasy Tissue Kit (QIAGEN, Valencia, CA). Total RNA was obtained from cell lines using RNeasy Plus Mini Kit (QIAGEN). cDNA was prepared as described previously [14].

Gene Sequencing

DNA sequencing and mutational analysis for EGFR (exons 18–21), HER2 (exons 19 and 20), KRAS (codons 12, 13, and 61), BRAF (exons 11 and 15) and PIK3CA (exons 9 and 20) were done as reported by us previously [10], [14], [32]. HER3 (exons 18–21) and HER4 (exons18–23) mutation status were analyzed by PCR amplification of genomic DNA or cDNA and direct sequencing of the PCR products [10]. The mutations in these genes were determined using the PCR primers as listed (Table S6). All PCRs were performed in 25 µL volumes containing 100 ng of DNA using HotStarTaq DNA polymerase (QIAGEN Inc., Valencia, CA). DNA was amplified for 32 to 34 cycles at 95°C for 30 seconds, 62°C to 68°C for 30 seconds, and 72°C for 30 seconds followed by 7 minutes extension at 72°C. All PCR products were incubated using exonuclease I and shrimp alkaline phosphatase (Amersham Biosciences Co., Piscataway, NJ) and sequenced directly using Applied Biosystems PRISM dye terminator cycle sequencing method (Perkin-Elmer Co., Foster City, CA). All mutations were confirmed by independent sequencing in both directions.

Copy number evaluation

Copy number gains can be evaluated by a number of methodologies. For clinical samples, fluorescence in situ hybridization (FISH) is frequently used as it can discriminate between tumor and non-malignant cells. Laboratory studies of tumor cell lines (which are free of non-malignant cells) often utilize quantitative PCR (qPCR). An alternative method is array comparative genomic hybridization (aCGH). Because direct comparisons of these methods have seldom been reported [33], we utilized all three methods.

Real-time qPCR

Gene copy numbers were determined using real-time quantitative PCR employing the Chromo4 PCR System (Bio-Rad Laboratories, Hercules, CA). To determine the copy number of a target gene, we used control genes located on the same chromosome as the target gene (Table S6). The resultant ratios were compared to similar ratios from diploid cells (the mean values of the eight HBEC cell lines. TaqMan methodology was used except for PIK3CA [14]. Primers and probes were chosen by TaqMan Primer Express™ 1.5 (Applied Biosystem, Foster City, CA). Primers were purchased from Invitrogen and probes from Biosearch Technologies (Novato, CA). The sequences of primers and probes are shown in Table S6. Standard curves were constructed with serial dilutions of specific PCR products. Amplification mixes (25 µl) contained the sample DNA (20 ng), 10× TaqMan buffer (2.5 µl), 200 µM dNTP, 1.25 U Hotstar Taq™ DNA polymerase, 200 nM each primer and 100 nM probe. The thermal cycling conditions comprised 5 min at 95°C and 40 cycles at 95°C for 15 s and 60°C for 30 s. All the samples were analyzed in triplicate. The parameter C_t is defined as the fractional cycle number at which the fluorescence generated by cleavage of the probe passes a fixed threshold above baseline. The target gene copy number in unknown samples is quantified by measuring C_t values and by using a standard curve to determine the copy number [34]. Gene copy number was calculated using the following equation: g

(S_T/S_C)/(N_T/N_C). PIK3CA copy number was assessed as described by us previously [14].

Tiling path aCGH

aCGH was performed as previously described [14], [15]. Genomic imbalances were identified using aCGH-Smooth [35] as previously described [36].

FISH assays

Dual-color FISH assays were performed according to a standard protocol [37], [38]. The probe sets, EGFR/CEP7 and PathVysion and the controls CEP7 were obtained from Abbott Molecular (Des Plaines, Illinois), HER3/CEP12 was obtained from QBiogene; BRAF probe was generated from the bacterial artificial clone (BAC) clone used for aCGH. The copy numbers of the target genes (labeled in red fluorophores) and the CEP probes (labeled in Spectrum Green) were verified in at least 100 interphase cells and 20 metaphase spreads. Images were captured and merged using the CytoVision workstation (Applied Imaging, San Jose, CA).

Sub-cloning

Genomic DNA was isolated from mutant EGFR or KRAS NSCLC cell lines and used as a template for PCR amplification of EGFR or KRAS. The primers and the conditions of PCR reactions were as described previously. PCR products were cloned into pCR2.1-TOPO vector using TOPO TA cloning kit (Invitrogen). About 20 clones (range15–25) were selected for sequencing using either M13 forward primer or corresponding EGFR or KRAS primers. The results are expressed as the percentages of mutant alleles present in the total number cloned.

TKI sensitivity

The MTS colorimetric assay (Promega) was performed as per the manufacturer's instructions. This assay is based on the conversion of MTS into soluble formazan by endogenous dehydrogenase enzymes found in metabolically active cells. Cells were plated at 1×10³–4×10³ cells/well in tissue culture treated 96-well plates. The following day, TKI (gefitinib or erlotinib) was added to each plate in a dilution series across the plate such that eight different concentrations of the drug were tested. On day 5, 20 µl of MTS was added, followed by a 1 hour incubation at 37°C and then the absorbance was read at 490 nm on a plate reader.

96-well plate data were imported into an in-house Database of In VItro drug Sensitivity Assays (DIVISA by Luc Girard, manuscript in preparation) where IC50s are calculated as well as various statistical analyses. To calculate IC50 values, the data were background-subtracted (columns 1 and 12 typically contained media with no cells or drugs and served as background values), and fitted to the DRC model (R package ‘drc’ by Christian Ritz and Jens Streibig, http://www.bioassay.dk) to generate a sigmoidal curve from which the concentration that inhibits 50% of the cells (IC50) was determined.

Statistical Analyses

Fisher's two-tailed exact tests were determined using the Prism 4 software (Graph Pad, San Diego, CA). P values<0.05 were considered significant. Other statistical analyses are discussed under appropriate categories in the Results section.

Results

Mutations (m) and CNGs (g) in lung cancer cell lines

We examined 32 SCLC and three lines from small cell cancers from extrapulmonary sites (extrapulmonary small cell cancers, ExPuSC), for somatic mutations and CNGs of EGFR pathway genes. We found a total of only three mutations in 35 cell lines and all three were PIK3CA mutations (Table S1). CNGs for EGFR pathway genes were infrequently encountered in SCLC (Table S2). Since somatic mutations and CNGs for EGFR pathway genes were rare in SCLC, we limited our further studies to NSCLC.

For any of the seven genes tested, mutations (39/77, 50.6%), copy number gains (50/77, 64.9%) or either (65/77, 84.4%) were frequent in NSCLC lines. These findings are discussed in greater detail below.

Mutations (m) of EGFR pathway genes in NSCLC

We examined the NSCLC lines, for somatic mutations of EGFR pathway genes as listed in the Methods section. We found a total of 40 mutations in 39 (50.6%) cell lines (Fig 1a, Table S3). These mutations consisted of 19 KRAS (24.7%), 10 EGFR (13%), five BRAF (6.5%), four PIK3CA (5.2%), one HER2 (1.3%), one HER4 (2.2%) and none for HER3. mEGFR, mBRAF and mHER2 were present exclusively in adenocarcinomas while mKRAS were more frequent in adenocarcinoma and large cell histologies. mPIK3CA were the only mutations that did not target adenocarcinoma histology (Fig. 1c).

Figure 1

Mutations (m) and Copy number gains (g) of EGFR pathway genes in NSCLC.

Mutations are mutually exclusive

We examined for any association between somatic mutations in individual cell lines. In order to test the hypothesis that EGFR pathway gene mutations are mutually exclusive, we used Monte Carlo simulations to calculate the empirical p-value. The null hypothesis is that mutations will occur independently. We simulated the joint distribution of mutation events among these seven genes using the observed marginal mutation rates and the independent assumption. In each simulation, we noted the number of cell lines with multiple mutations. We repeated the simulations 10,000 times and obtained the empirical distribution of the number of cell lines with multiple mutations; then we compared the observed number of cell lines with multiple mutations with this empirical distribution to calculate the p-value.

For this study, the one-sided p-value is 0.0014; therefore, we rejected the null hypothesis and concluded that gene mutations are mutually exclusive events for these genes in NSCLC cell lines, with one exception (Table 1). Large cell carcinoma line NCI-H460 had both KRAS and PIK3CA mutations.

Table 1

Are mutations of EGFR pathway genes mutually exclusive?

Comparison of methods for determining copy number gains

aCGH, FISH and qPCR were used to test the gene copy numbers for NSCLC cell lines. There is no clear biological threshold value for defining the abnormal copy numbers for NSCLC human cell lines; we selected the threshold values which could achieve the best positive or negative category agreement among the three tests. Specifically, we computed the Kappa statistics [39] between aCGH and qPCR tests over two dimensional cut-off grids as 2, 2.5, 3, 3.5, 4, 4.5 and 5; and then did the same for FISH and qPCR. The cut-off values that achieved the best agreement among three tests were as follows: 4 for qPCR, 3 for CGH and 4.5 for FISH; and we used these values to define the presence of copy number gains in NSCLC cell lines. Using these cut-off values and all available data, there was highly significant concordance (p<0.001) between the three methods as shown in Table 2.

Table 2

Comparison of the three methods used for studying gene copy number.

EGFR pathway gene CNGs

As the qPCR data were complete for all lines (while subsets were tested by the other two methods which are more laborious and expensive), we used qPCR data for all further analyses. Copy number gains were frequent (>10%) for EGFR, HER2, HER3 and KRAS (Fig 1a, Table S4). In particular gains for EGFR were more than twice as frequent (40%) than for any other gene. In general, with the exceptions of gBRAF (limited to adenocarcinoma histology), and gPIK3CA (largely limited to squamous cell histology), CNGs did not show any apparent histology subtype bias (Fig. 1c).

Relationship between Mutations and Copy number gains

Unlike mutations, copy number gains were not mutually exclusive either with other CNGs or with mutations (Fig. 2). We found CNGs for two or more genes in 32.5% (25/77) of cell lines (Table S4). However, we noted a highly significant association between mutations of EGFR or KRAS and CNGs of their respective genes (Fig 3a, 3b).

Figure 2

Copy number gains are not mutually exclusive with either other copy number gains or with mutations.

Figure 3

Mutant allele specific imbalance (MASI) for EGFR and KRAS genes.

Mutant allele specific imbalance (MASI) for EGFR and KRAS genes

As noted above CNGs of EGFR and KRAS were more frequent in cell lines harboring these mutations. Employing the sub-cloning method previously described, we investigated whether the EGFR and KRAS genes preferentially demonstrated mutant allele specific CNGs in cell lines harboring the respective mutations (Fig. 3). The data, for testing whether mutant alleles are amplified in gene mutant cell lines are clustered data with binary outcome. Each cell line is a cluster and the null hypothesis is that the mutation rate is 0.5. Because the number of subclones in each cell line varied, we used the analysis approach for the clustered binary outcome with various cluster sizes [40]. In mutant lines, the mutant allele was almost always in excess compared to the wild type allele, a phenomenon we have termed MASI. The p-value for EGFR mutant cell lines was 0.019 and for KRAS mutant cell line was 0.0003 indicating that the mutant alleles were preferentially dominant in both cases. Most MASI cases were due to increased copy number of the mutant allele. However, in a minority of cases MASI was also noted in cells having diploid or near diploid copy numbers (acquired uniparental disomy). Thus we used the term MASI as opposed to mutant allele specific gains.

Characteristics of EGFR Mutant Cell Lines

The clinico-pathologic and molecular data for the 10 NSCLC cell lines which harbor EGFR mutations are summarized in Table 3. All contained one of the two major mutations in the kinase domain, either deletions of exon 19 (n

7) or L858R mutation in exon 21 (n

3). Mutations were exclusively seen in adenocarcinoma and never smokers or patients with low tobacco exposure (≤10 pack years). One EGFR mutant cell line was developed in Japan and the remaining nine were developed in North America. Of the nine developed in North America, only one was from an individual of Asian ethnicity. In addition, we identified that EGFR mutations were more common in comparatively younger age group (<55 years).

Table 3

Characteristics of EGFR Mutant Cell Lines.

Seven of these cell lines were sensitive to gefitinib. The three resistant lines had a second mutation: either the T790M resistance associated mutation in exon 20 (n

2) or a homozygous deletion of the PTEN gene and absence of its protein [41] (authors' unpublished observations).

Effects of mutations and copy number gains on sensitivity to TKIs

To evaluate the effect of mutations and copy number gains on sensitivity to TKIs, we analyzed the IC50 values of 45 NSCLC cell lines. Because TKI clinical sensitivity preferentially targets adenocarcinoma histology, we included 31 adenocarcinomas. The entire subset included all of the EGFR, HER2 and HER4 mutant cell lines and 12 KRAS and 3 BRAF cell lines. Because PIK3CA mutations favored non-adenocarcinoma histology only one PIK3CA mutant cell line was included. Of the entire subset 17 lines were wild type for all genes tested (Table S5).

We found an excellent concordance between the IC50 values between gefitinib and erlotinib (p<0.0001) (Fig 4, Table S7). Because our data set for gefitinib sensitivity was more extensive we elected to utilize the gefitinib data for further analyses.

Figure 4

Concordance between IC50 values for gefitinib vs erlotinib.

Fig. 5 illustrates the sensitivity patterns of the tested cell lines. The in vitro concentration of 1 µM used in tissue culture correlates to the plasma concentration of gefitinib in patients treated with the standard dose of gefitinib i.e. 250 mg a day. This threshold has been used by researchers in the past to distinguish sensitive from insensitive and/or resistant cell lines [7]. On the basis of the above-mentioned threshold and the shape of the curve, we divided our cell lines into 3 categories: sensitive (≤1 µM), intermediate (>1 µM but ≤10 µM) and resistant (>10 µM) as seen in figure Fig. 5. The IC50 values follow a line that demonstrated a slope change point at approximately 10 µM, and thus demonstrated an apparently bimodal distribution (Fig. 5). As cell lines with values below 1 µM were regarded as sensitive, we arbitrarily categorized values between 1 and 10 µM as being intermediate in sensitivity.

Figure 5

Rank Order of NSCLC Cell Line depending on the Iressa IC50.

There were nine cell lines in the sensitive group and they included seven of the ten mEGFR lines, one gEGFR cell line and one gHER2 cell line. The intermediate category consisted of only two cell lines; a mEGFR cell line having a secondary T790M mutation and a gHER2 cell line. The resistant cell line category was the largest and included two mEGFR cell lines, one having the secondary T790M mutation and the other having homozygous deletion of PTEN gene. The remaining resistant lines included all of the wild type lines and all lines having KRAS, BRAF, HER2, HER4 and PIK3CA mutations.

Using univariate tests, we analyzed the association between gefitinib sensitivity and mutational status or CNGs of the various EGFR pathway genes and found a significant correlation between gefitinib sensitivity with EGFR mutation (p

0.002) and EGFR copy number gains (p

0.001). We tested the hypothesis that KRAS mutations confer intrinsic resistance to TKIs. Because mEGFR was associated with sensitivity to TKIs we excluded these lines from our analyses. We compared the rank order of IC50s of mKRAS cell lines with a) all other cell lines and b) cell lines wild type for all other tested genes. Our univariate analyses demonstrated that indeed KRAS mutations conferred in vitro resistance to gefitinib (Table 5). None of the other gene mutations or copy number gains showed any significant correlation with gefitinib sensitivity.

Table 5

Do KRAS mutations confer resistance to TKIs?

We then used a multivariate linear regression model [42] to explore the association between gefitinib sensitivity and EGFR mutation, EGFR copy number and HER2 copy number (Table 4). The response variable was log IC50 values of gefitinib (as a continuous variable), the predictors are EGFR mutation status (mutation or wild type EGFR), EGFR copy number (measured by qPCR) and HER2 copy number (measured by qPCR) as continuous variables. After adjusting for the effects of EGFR and HER2 copy numbers, the EGFR mutation status was significantly associated with the gefitinib sensitivity (p<0.001). Similarly after adjusting for the effects of the other two variables, the EGFR copy number (p value

0.002) and HER2 copy number (p

0.021) were independently and significantly associated with the gefitinib sensitivity. To summarize the multivariate analysis, all 3 parameters, i.e. mEGFR mutations, gEGFR and gHER2, showed a correlation with TKI sensitivity in decreasing order of importance. We also repeated the analysis using the copy number of EGFR and HER2 as binary variables. In these analyses, EGFR copy number remained significantly associated with gefitinib sensitivity (p value<0.001) while HER2 copy number was no longer significantly associated (p value>0.05; Table 4).

Table 4

Results of Multivariate Regression Analyses.

Discussion

In order to clarify some of the conflicting data from clinical studies, we performed an in vitro study of deregulation of the EGFR signaling pathway in a large panel of lung cancer cell lines and correlated the results with TKI sensitivity. In contrast to SCLC, mutations or CNGs of the studied genes were frequent in NSCLC – 84.4% of NSCLC cell lines had a mutation, one or more CNGs or both. Thus the EGFR pathway is deregulated at high frequency in NSCLC.

Somatic mutations or CNGs of the studied genes were rare in SCLC and further studies were limited to NSCLC. About half of the NSCLC lines had a mutation in one of the genes of the EGFR pathway. The relative distribution of the mutations of EGFR pathway genes in established NSCLC cell lines were very similar to those reported from large clinical studies [3], [8], [10], [12], [14], [43]. EGFR (13%) and KRAS (24.7%) mutations were more frequently observed than other gene mutations. We found a single mutation of HER2 and HER4 and no mutations of EGFR3. Mutations of BRAF and PIK3CA were intermediate in frequency. In general mutations were more common in adenocarcinoma subtype and relatively rare in squamous cell carcinomas except for PIK3CA which was more common in squamous cell carcinomas, as also occurs in tumors [14]. Mutations were mutually exclusive except for a single cell line that had both KRAS and PIK3CA mutations. We have previously reported that mutations of EGFR, KRAS, BRAF and HER2 were mutually exclusive in 691 resected NSCLC tumors indicating that a single activating mutation in the EGFR-RAS-RAF signaling pathway may be sufficient for the pathogenesis of many lung cancers [10], [14]. However, the mutational status of PIK3CA was not mutually exclusive, which is similar to another recent report [14]. We have postulated that the PIK3CA and EGFR signaling pathways closely interact but that PI3K signaling represents a partially independent pathway [14].

Copy number gains are believed to be an important mechanism for activation of the EGFR pathway genes and the downstream signaling network. CNGs have been difficult to interpret due to the differences stemming from the use of different methods. Our comparison of the three methods yielded excellent concordance for binary values after we used kappa statistics to determine the optimal cut off values for increased copy numbers. However, the three methods we compared have their individual advantages and disadvantages. CNGs result from gene over-representation due to chromosomal polysomy or focal amplification. FISH, a widely used technique to assess gene copy number takes can identify gains due to both mechanisms. On the other hand, not only is it expensive and time consuming but it analyzes only a small subset of cells (~100 cells) to assess the copy number, and tandem segmental duplications may not always be detected. Thus intratumoral heterogeneity is difficult to assess. Both aCGH and qPCR analyze DNA from many thousands of cells, but because they compare values for the test gene against an internal standard (a reference locus), they can detect focal amplifications more readily than polysomy. In the case of qPCR, if the amplicon encompasses the reference locus, then amplification may be underestimated. In the case of aCGH, the copy number deduced from the signal ratio can appear lower than those obtained by the other two techniques, for aCGH measures relative and not absolute copy number. It does not take into account changes in ploidy and therefore will dampen ratio shifts for copy number gains in samples with a high DNA index [33]. Thus, the absolute copy number for a given chromosome region is best determined by using a combination of complementary methods [44]. Despite the shortcomings of each method, we found a high concordance among the three techniques when using binary values (normal/increased), similar to findings in a recent comparison [33]. We had most extensive data for gene copy number generated via qPCR and used it for all our further analysis.

Sixty seven percent of the NSCLC lines had CNGs for one or more EGFR pathway genes. CNGs for EGFR (40.3%), HER2 (18.2%) and KRAS (14.3%) were more frequent than the other genes. BRAF CNGs (13%) were more common in adenocarcinoma subtype and PIK3CA CNGs (36.4%) were more frequent in squamous cell carcinoma as compared to the other subtypes. Other genes did not show subtype bias. Unlike mutations, CNGs were not mutually exclusive either with other CNGs or with mutations.

However, CNGs were not completely random. We identified that EGFR and KRAS CNGs were significantly more frequent in EGFR and KRAS mutant cell lines respectively. Using subcloning, we determined that in almost all mutant lines with CNGs, the mutant allele was in great excess in comparison to the wild type allele. In some lines with diploid amounts of the mutant gene, the wild type allele was absent or in minute amounts. This finding for oncogenes (as compared to tumor suppressor genes) has been described, particularly in hematologic malignancies [45], [46], and represents a form of acquired uniparental disomy. Taken together, we term this phenomenon as mutant allele specific imbalance (MASI). Thus MASI may result from specific amplification of the mutant allele, or by loss of the wild type gene, or by a combination of these events. Of interest, wild type KRAS may function as an inhibitor of tumorigenesis and thus represent a form of tumor suppressor gene [47]. These results suggest that the combination of two methods of activating the oncogenes EGFR and KRAS (i.e. mutation and MASI) may confer a greater growth or survival advantage to the malignant cell than a single method. These findings also indicate that mutations must precede CNGs in cells harboring both changes. Other evidence for the combination of mutations and CNGs in tumor cells exists [48], [49], [50], [51], and for the concept that mutations are the initial event [48], [49], [50], [52], [53].

EGFR mutations were exclusively found in adenocarcinoma histology and in never smokers or individuals with low smoke exposure (<15 PYR). In addition EGFR mutations were more common in a comparatively younger age group (<55 years). EGFR mutations did not demonstrate gender bias. Most mutant lines were derived from Caucasians, which is not surprising considering that 90% were established in the USA. In rank order of sensitivity to gefitinib, nine of the EGFR mutant lines were represented in the 12 most sensitive lines, in keeping with the previously discussed importance of mutations in sensitivity to TKIs.

We found an excellent concordance between the IC50 values between gefitinib and erlotinib. Depending on the previously determined in vitro threshold for gefitinib and the bimodal curve we classified our NSCLC cell lines into 3 categories: sensitive, intermediate and resistant. We had nine sensitive cell lines which included seven of the 10 EGFR mutant cell lines, one cell line with EGFR CNG and one with HER2 CNG. Of the remaining 3 EGFR mutant cell lines, two had the resistance associated secondary T790M mutation. The third resistant EGFR mutant line had a deletion of the PTEN gene, a finding associated with TKI resistance in other systems [54].

Pao et al demonstrated that mutations in KRAS are associated with a lack of sensitivity to TKIs [28]. We tested the hypothesis that KRAS mutations confer intrinsic resistance to TKIs and our univariate analyses demonstrated that indeed KRAS mutations conferred in vitro resistance to gefitinib. Thus KRAS mutations are associated with both clinical and in vitro resistance to gefitinib.

None of the other gene mutations or copy number gains showed any significant correlation with gefitinib sensitivity. In the clinical setting, approximately 10–20% of NSCLC patients who do not harbor identifiable EGFR mutations respond partially to gefitinib. Evidently EGFR mutations are one of the most important but not the sole determinant of TKI response.

We evaluated the effect of EGFR mutation, EGFR CNGs and HER2 CNGs on TKI sensitivity using a multivariate regression analysis. To summarize the multivariate analysis, all 3 parameters, i.e. EGFR mutations, EGFR CNG and HER2 CNG, showed a correlation with the TKI sensitivity (used as a continuous variable) in decreasing order of importance. Thus three previously identified factors related to patient response to TKIs are major factors in the in vitro response. However, the most important independent factor is the presence of activating EGFR mutations. Our findings may assist in the prediction of response and the selection of patients for targeted therapies. Another important finding was the frequent presence of selective imbalance of the mutant form (MASI) of oncogenes in tumor cells that harbor mutations of EGFR and KRAS. While MASI may confer a selective advantage to the tumor cell, its effect on clinical course or response to therapy remains to be determined.

Supporting Information

Table S1

(0.01 MB PDF)

Click here for additional data file.^{(10K, pdf)}

Table S2

(0.01 MB PDF)

Click here for additional data file.^{(12K, pdf)}

Table S3

(0.01 MB PDF)

Click here for additional data file.^{(12K, pdf)}

Table S4

(0.02 MB PDF)

Click here for additional data file.^{(15K, pdf)}

Table S5

(0.02 MB PDF)

Click here for additional data file.^{(16K, pdf)}

Table S6

(0.01 MB PDF)

Click here for additional data file.^{(14K, pdf)}

Table S7

(0.02 MB PDF)

Click here for additional data file.^{(16K, pdf)}

Footnotes

Competing Interests: Dr. Gazdar is a paid consultant/lecturer for AstraZeneca PLC. Dr. Garcia recieves Research Funding >10,000 from AstraZeneca, Genentech and OSI; Honorarium <10,000 from Roche. Dr. Minna receives research support from AstraZeneca PLC.

Funding: Funding was recieved from the Specialized Program of Research Excellence in Lung Cancer (P50CA70907) and Early Detection Research Network, National Cancer Institute, Bethesda, Maryland (U01CA084971). Funds from both grants were utilized for salary support and for performance of assays. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

References

1. Sun S, Schiller JH, Gazdar AF. Lung cancer in never smokers–a different disease. Nat Rev Cancer. 2007;7:778–790. [PubMed]

2. Jemal A, Siegel R, Ward E, Hao Y, Xu J, et al. Cancer statistics, 2008. CA Cancer J Clin. 2008;58:71–96. [PubMed]

3. Shigematsu H, Gazdar AF. Somatic mutations of epidermal growth factor receptor signaling pathway in lung cancers. Int J Cancer. 2006;118:257–262. [PubMed]

4. Rowinsky EK. The erbB family: targets for therapeutic development against cancer and therapeutic strategies using monoclonal antibodies and tyrosine kinase inhibitors. Annu Rev Med. 2004;55:433–457. [PubMed]

5. Hirsch FR, Varella-Garcia M, Bunn PA, Jr, Di Maria MV, Veve R, et al. Epidermal growth factor receptor in non-small-cell lung carcinomas: correlation between gene copy number and protein expression and impact on prognosis. J Clin Oncol. 2003;21:3798–3807. [PubMed]

6. Nicholson RI, Gee JM, Harper ME. EGFR and cancer prognosis. Eur J Cancer. 2001;37(Suppl 4):S9–15. [PubMed]

7. Sharma SV, Bell DW, Settleman J, Haber DA. Epidermal growth factor receptor mutations in lung cancer. Nat Rev Cancer. 2007;7:169–181. [PubMed]

8. Eberhard DA, Johnson BE, Amler LC, Goddard AD, Heldens SL, et al. Mutations in the epidermal growth factor receptor and in KRAS are predictive and prognostic indicators in patients with non-small-cell lung cancer treated with chemotherapy alone and in combination with erlotinib. J Clin Oncol. 2005;23:5900–5909. [PubMed]

9. Gazdar AF, Shigematsu H, Herz J, Minna JD. Mutations and addiction to EGFR: the Achilles ‘heal’ of lung cancers? Trends Mol Med. 2004;10:481–486. [PubMed]

10. Shigematsu H, Takahashi T, Nomura M, Majmudar K, Suzuki M, et al. Somatic mutations of the HER2 kinase domain in lung adenocarcinomas. Cancer Res. 2005;65:1642–1646. [PubMed]

11. Soung YH, Lee JW, Kim SY, Wang YP, Jo KH, et al. Somatic mutations of the ERBB4 kinase domain in human cancers. Int J Cancer. 2006;118:1426–1429. [PubMed]

12. Naoki K, Chen TH, Richards WG, Sugarbaker DJ, Meyerson M. Missense mutations of the BRAF gene in human lung adenocarcinoma. Cancer Res. 2002;62:7001–7003. [PubMed]

13. Samuels Y, Wang Z, Bardelli A, Silliman N, Ptak J, et al. High frequency of mutations of the PIK3CA gene in human cancers. Science. 2004;304:554. [PubMed]

14. Yamamoto H, Shigematsu HNM, Lockwood WW, Sato M, et al. PIK3CA mutations and gene copy number in lung cancers. Cancer Res. 2008;68

15. Lockwood WW, Coe BP, Williams AC, MacAulay C, Lam WL. Whole genome tiling path array CGH analysis of segmental copy number alterations in cervical cancer cell lines. Int J Cancer. 2007;120:436–443. [PubMed]

16. Pawson T. Regulation of the Ras signalling pathway by protein-tyrosine kinases. Biochem Soc Trans. 1994;22:455–460. [PubMed]

17. Chen WS, Lazar CS, Poenie M, Tsien RY, Gill GN, et al. Requirement for intrinsic protein tyrosine kinase in the immediate and late actions of the EGF receptor. Nature. 1987;328:820–823. [PubMed]

18. Fukuoka M, Yano S, Giaccone G, Tamura T, Nakagawa K, et al. Multi-institutional randomized phase II trial of gefitinib for previously treated patients with advanced non-small-cell lung cancer (The IDEAL 1 Trial) [corrected]. J Clin Oncol. 2003;21:2237–2246. [PubMed]

19. Kris MG, Natale RB, Herbst RS, Lynch TJ, Jr, Prager D, et al. Efficacy of gefitinib, an inhibitor of the epidermal growth factor receptor tyrosine kinase, in symptomatic patients with non-small cell lung cancer: a randomized trial. JAMA. 2003;290:2149–2158. [PubMed]

20. Shepherd FA, Rodrigues Pereira J, Ciuleanu T, Tan EH, Hirsh V, et al. Erlotinib in previously treated non-small-cell lung cancer. N Engl J Med. 2005;353:123–132. [PubMed]

21. Lynch TJ, Bell DW, Sordella R, Gurubhagavatula S, Okimoto RA, et al. Activating mutations in the epidermal growth factor receptor underlying responsiveness of non-small-cell lung cancer to gefitinib. N Engl J Med. 2004;350:2129–2139. [PubMed]

22. Paez JG, Janne PA, Lee JC, Tracy S, Greulich H, et al. EGFR mutations in lung cancer: correlation with clinical response to gefitinib therapy. Science. 2004;304:1497–1500. [PubMed]

23. Pao W, Miller V, Zakowski M, Doherty J, Politi K, et al. EGF receptor gene mutations are common in lung cancers from “never smokers” and are associated with sensitivity of tumors to gefitinib and erlotinib. Proc Natl Acad Sci U S A. 2004;101:13306–13311. [PubMed]

24. Uramoto H, Mitsudomi T. Which biomarker predicts benefit from EGFR-TKI treatment for patients with lung cancer? Br J Cancer. 2007;96:857–863. [PMC free article] [PubMed]

25. Shepherd FA, Rosell R. Weighing tumor biology in treatment decisions for patients with non-small cell lung cancer. J Thorac Oncol. 2007;2(Suppl 2):S68–76. [PubMed]

26. Cappuzzo F, Varella-Garcia M, Shigematsu H, Domenichini I, Bartolini S, et al. Increased HER2 gene copy number is associated with response to gefitinib therapy in epidermal growth factor receptor-positive non-small-cell lung cancer patients. J Clin Oncol. 2005;23:5007–5018. [PubMed]

27. Cappuzzo F, Toschi L, Domenichini I, Bartolini S, Ceresoli GL, et al. HER3 genomic gain and sensitivity to gefitinib in advanced non-small-cell lung cancer patients. Br J Cancer. 2005;93:1334–1340. [PMC free article] [PubMed]

28. Pao W, Wang TY, Riely GJ, Miller VA, Pan Q, et al. KRAS mutations and primary resistance of lung adenocarcinomas to gefitinib or erlotinib. PLoS Med. 2005;2:e17. [PMC free article] [PubMed]

29. Galanis E, Frytak S, Lloyd RV. Extrapulmonary small cell carcinoma. Cancer. 1997;79:1729–1736. [PubMed]

30. Phelps RM, Johnson BE, Ihde DC, Gazdar AF, Carbone DP, et al. NCI-Navy Medical Oncology Branch cell line data base. J Cell Biochem. 1996;(Suppl 24):32–91. [PubMed]

31. Ramirez RD, Sheridan S, Girard L, Sato M, Kim Y, et al. Immortalization of human bronchial epithelial cells in the absence of viral oncoproteins. Cancer Res. 2004;64:9027–9034. [PubMed]

32. Shigematsu H, Lin L, Takahashi T, Nomura M, Suzuki M, et al. Clinical and biological features associated with epidermal growth factor receptor gene mutations in lung cancers. J Natl Cancer Inst. 2005;97:339–346. [PubMed]

33. Lyng H, Lando M, Brovig RS, Svendsrud DH, Johansen M, et al. GeneCount: genome-wide calculation of absolute tumor DNA copy numbers from array comparative genomic hybridization data. Genome Biol. 2008;9:R86. [PMC free article] [PubMed]

34. Bieche I, Onody P, Laurendeau I, Olivi M, Vidaud D, et al. Real-time reverse transcription-PCR assay for future management of ERBB2-based clinical applications. Clin Chem. 1999;45:1148–1156. [PubMed]

35. Jong K, Marchiori E, Meijer G, Vaart AV, Ylstra B. Breakpoint identification and smoothing of array comparative genomic hybridization data. Bioinformatics. 2004;20:3636–3637. [PubMed]

36. Coe BP, Lockwood WW, Girard L, Chari R, Macaulay C, et al. Differential disruption of cell cycle pathways in small cell and non-small cell lung cancer. Br J Cancer. 2006;94:1927–1935. [PMC free article] [PubMed]

37. Varella-Garcia M. Molecular cytogenetics in solid tumors: laboratorial tool for diagnosis, prognosis, and therapy. Oncologist. 2003;8:45–58. [PubMed]

38. Helfrich BA, Raben D, Varella-Garcia M, Gustafson D, Chan DC, et al. Antitumor activity of the epidermal growth factor receptor (EGFR) tyrosine kinase inhibitor gefitinib (ZD1839, Iressa) in non-small cell lung cancer cell lines correlates with gene copy number and EGFR mutations but not EGFR protein levels. Clin Cancer Res. 2006;12:7117–7125. [PubMed]

39. Everit B. Statistical methods for medical investigations. London, U.K: Edward Arnold; 1994.

40. Jung SH, Kang SH, Ahn C. Sample size calculations for clustered binary data. Stat Med. 2001;20:1971–1982. [PubMed]

41. Guo A, Villen J, Kornhauser J, Lee KA, Stokes MP, et al. Signaling networks assembled by oncogenic EGFR and c-Met. Proc Natl Acad Sci U S A. 2008;105:692–697. [PubMed]

42. Netter J, Kutner M, Nachtsheim C, Wasserman W. Applied linear statistical model. McGraw-Hill/Irwin; 1996.

43. Marks JL, McLellan MD, Zakowski MF, Lash AE, Kasai Y, et al. Mutational analysis of EGFR and related signaling pathway genes in lung Adenocarcinomas identifies a novel somatic kinase domain mutation in FGFR4. PLoS ONE. 2007;2:e426. [PMC free article] [PubMed]

44. Rosenberg C, Schut TB, Mostert MC, Tanke HJ, Raap AK, et al. Comparative genomic hybridization in hypotriploid/hyperdiploid tumors. Cytometry. 1997;29:113–121. [PubMed]

45. Gondek LP, Dunbar AJ, Szpurka H, McDevitt MA, Maciejewski JP. SNP array karyotyping allows for the detection of uniparental disomy and cryptic chromosomal abnormalities in MDS/MPD-U and MPD. PLoS ONE. 2007;2:e1225. [PMC free article] [PubMed]

46. Tyybakinoja A, Elonen E, Vauhkonen H, Saarela J, Knuutila S. Single nucleotide polymorphism microarray analysis of karyotypically normal acute myeloid leukemia reveals frequent copy number neutral loss of heterozygosity. Haematologica. 2008;93:631–632. [PubMed]

47. Zhang Z, Wang Y, Vikis HG, Johnson L, Liu G, et al. Wildtype Kras2 can inhibit lung carcinogenesis in mice. Nat Genet. 2001;29:25–33. [PubMed]

48. Gazdar AF, Minna JD. Deregulated EGFR Signaling during Lung cancer Progression: Mutations, Amplicons and Autocrine loops. Cancer Prevention Research. 2008;1:6. [PMC free article] [PubMed]

49. Nomura M, Shigematsu H, Li L, Suzuki M, Takahashi T, et al. Polymorphisms, Mutations, and Amplification of the EGFR Gene in Non-Small Cell Lung Cancers. PLoS Med. 2007;4:e125. [PMC free article] [PubMed]

50. Soh J, Toyooka S, Ichihara S, Asano H, Kobayashi N, et al. Sequential molecular changes during multistage pathogenesis of small peripheral adenocarcinomas of the lung. J Thorac Oncol. 2008;3:340–347. [PMC free article] [PubMed]

51. Yatabe Y, Takahashi T, Mitsudomi T. Epidermal growth factor receptor gene amplification is acquired in association with tumor progression of EGFR-mutated lung cancer. Cancer Res. 2008;68:2106–2111. [PubMed]

52. Zhang J, Iwanaga K, Choi KC, Wislez M, Raso MG, et al. Intratumoral epiregulin is a marker of advanced disease in non-small cell lung cancer patients and confers invasive properties on EGFR-mutant cells. Cancer Prevention Research. 2008;1:201–207. [PMC free article] [PubMed]

53. Tang X, Varella-Garcia M, Xavier AC, Massarelli E, Ozburn NC, et al. Epidermal growth factor receptor abnormalities in the pathogenesis and progression of lung adenocarcinomas. Cancer Prevention Research. 2008;1:9. [PMC free article] [PubMed]

54. Mellinghoff IK, Cloughesy TF, Mischel PS. PTEN-mediated resistance to epidermal growth factor receptor kinase inhibitors. Clin Cancer Res. 2007;13:378–381. [PubMed]

Articles from PLoS ONE are provided here courtesy of
Public Library of Science