Skip to main content

Transcriptional changes associated with resistance to inhibitors of epidermal growth factor receptor revealed using metaanalysis

Abstract

Background

EGFR is important in maintaining metabolic homeostasis in healthy cells, but in tumors it activates downstream signaling pathways, causing proliferation, angiogenesis, invasion and metastasis. Consequently, EGFR is targeted in cancers using reversible, irreversible or antibody inhibitors. Unfortunately, tumors develop inhibitor resistance by mutations or overexpressing EGFR, or its ligand, or activating secondary, EGFR-independent pathways.

Methods

Here we present a global metaanalysis comparing transcriptional profiles from matched pairs of EGFR inhibitor-sensitive vs. -resistant cell lines, using 15 datasets comprising 274 microarrays. We also analyzed separately pairs of cell lines derived using reversible, irreversible or antibody inhibitors.

Results

The metaanalysis identifies commonalities in cell lines resistant to EGFR inhibitors: in sensitive cell lines, the ontological categories involving the ErbB receptors pathways, cell adhesion and lipid metabolism are overexpressed; however, resistance to EGFR inhibitors is associated with overexpression of genes for ErbB receptors-independent oncogenic pathways, regulation of cell motility, energy metabolism, immunity especially inflammatory cytokines biosynthesis, cell cycle and responses to exogenous and endogenous stimuli. Specifically in Gefitinib-resistant cell lines, the immunity-associated genes are overexpressed, whereas in Erlotinib-resistant ones so are the mitochondrial genes and processes. Unexpectedly, lines selected using EGFR-targeting antibodies overexpress different gene ontologies from ones selected using kinase inhibitors. Specifically, they have reduced expression of genes for proliferation, chemotaxis, immunity and angiogenesis.

Conclusions

This metaanalysis suggests that ‘combination therapies’ can improve cancer treatment outcomes. Potentially, use of mitochondrial blockers with Erlotinib, immunity blockers with Gefitinib, tyrosine kinase inhibitors with antibody inhibitors, may have better chance of avoiding development of resistance.

Peer Review reports

Background

Cancer is principally caused by changes in three types of genes i.e. oncogenes, tumor suppressor genes and DNA stability genes; common environmental factors contributing to these changes could be smoking, oncogenic viruses, occupational and environmental carcinogens and predisposing genetic polymorphisms [1,2]. Massive research efforts are ongoing to find treatments for cancer, still an unresolved problem and of course a heavy burden on health care. Targeting specific pathways and modulating the immune system are key strategies to control cancer progression and increase effectiveness of treatment [3].

EGFR, important growth factor receptor implicated in many cancers, is one of the targets for chemotherapy. EGFR belongs to ErbB family of tyrosine kinase receptors [4]. In tumors these receptors are activated by increased expression or structural rearrangement of receptor gene, mutation in ligand binding or tyrosine kinase domain or by the production of autocrine and paracrine ligands. On activation, EGFR dimerizes and triggers Ras-RAF-MEK-ERK-MAPK, JAK-STAT and other signaling cascades [4,5]. These pathways activate transcription factors, ultimately resulting in activation of cellular processes including proliferation, carbohydrate utilization, protein synthesis, angiogenesis, cell growth and cell survival [4,6]. In cancer cells, EGFR activation is important in maintaining the metabolic homeostasis and stimulates proliferation, invasion, angiogenesis, survival, decreased apoptosis, migration, differentiation and adhesion. Because of its central signaling position, EGFR is targeted in number of malignancies e.g. lung, colorectal, pancreatic, head and neck cancers, glioblastomas etc. [7].

To treat malignancies, EGFR activity is targeted with reversible, irreversible or antibody inhibitors and their combinations. Reversible inhibitors, e.g., Gefitinib and Erlotinib, compete for the intracellular ATP binding site of the kinase; the irreversible inhibitors, e.g., PF299804 and WZ4002, block the ATP binding site by covalent interaction with Cys773 in EGFR [8]. Antibody inhibitors block the extracellular ligand binding domain of EGFR thereby preventing ligand binding and receptor dimerization. Different types of inhibitors generate different transcriptional responses in EGFR-targeted cells [9].

Unfortunately, after a period, tumors inexorably develop resistance to inhibitors by e.g., overexpression of EGFR ligand, activating mutations in the tyrosine kinase or the ligand-binding domain, or by mutations in downstream or parallel signaling pathways, e.g., Axl or IGFIR [10-14]. Multiple studies focused on defining the secondary mutations that cause resistance to EGFR inhibitors [13-15]. However, much less attention has been paid to the transcriptional and metabolic changes that distinguish the resistant cells from the sensitive ones [12]. Therefore, we decided to explore the fundamental functional changes that distinguish the resistant cells from the sensitive ones, using transcriptional profiling. To cast a wide net, we used metaanalysis approach to find the differential gene expression between matched pairs of cell lines sensitive and resistant to EGFR inhibitors. Here we included eight studies with 15 distinct datasets directly comparing the transcriptional profiles in EGFR inhibitor-sensitive vs. resistant cell lines. The cell lines included non-small cell lung cancer, head and neck cancer, and epidermoid carcinoma cell lines. The inhibitors included both reversible and irreversible kinase inhibitors, as well as antibodies.

We found that in EGFR inhibitor-sensitive cell lines characteristically overexpressed gene ontologies are adhesion, negative regulation of cell proliferation, lipid metabolism and oncogenic processes involving ErbB receptors. But when cells become resistant, ontological categories associated with energy metabolism, immunity involving overexpressing inflammatory cytokines, responses to external and internal stimuli, proliferation and ErbB-independent oncogenic pathways are overexpressed. The specific resistance to Gefitinib apparently develops by overexpressing immunomodulatory genes; resistance to Erlotinib by energy producing mitochondrial pathways; resistance to irreversible inhibitors by overexpressing EGFR ligands, whereas resistance to antibody inhibitors develops differently from the resistance to tyrosine kinase inhibitors.

Methods

Downloading the data files

The overall flowchart of our methodology is graphically represented in Additional file 1: Figure S1.

Different microarray platforms used for transcriptional profiling produced different, characteristic data files, which were worked up separately and then synchronized. The CEL or TXT files deposited in these studies were first downloaded and unzipped. For each study, data obtained from sensitive and resistant cell lines were saved in different columns of excel spread sheets. Datasets obtained from Affymetrix studies were combined and analyzed using RMAExpress for quality control [16,17]. For non-Affymetrix studies, where we could not run RMAExpress quality control, we downloaded already normalized, _RAW.tar files and used these without further modifications, as submitted by the original authors.

Grouping studies for analysis using RankProd software

RankProd package analyses gene expression microarray data specifically to identify differentially expressed genes. RankProd uses non-parametric rank product method to detect genes that are consistently found among the most strongly upregulated ones and the most strongly downregulated ones in a number of replicate experiments, comparing two different condition [18]. We have combined into a single spreadsheet microarray data for sensitive and resistant cell lines with 20552 common genes in all datasets using data-loader [17]. Five datasets comprising 214 microarrays and 28235 genes for Gefitinib-sensitive and resistant cell lines were combined into a single excel spreadsheet and analyzed using RankProd. Differentially expressed genes in each of the class were recorded. Microarray data for the seven datasets comprising forty Erlotinib-sensitive and resistant microarrays, having 32062 common genes were combined for analysis using RankProd software [17].

We have pooled and compared the microarray data for EGFR irreversible inhibitors from two datasets, fourteen microarrays and 21631 common genes. For studying EGFR antibody inhibitors responses we found a single study with 3 microarrays from Cetuximab-sensitive and 3 from resistant cell lines, with 48607 genes.

We used the RankProd Software to find out the genes differentially expressed in EGFR inhibitor-sensitive and resistant cell lines with p-values better than 10−4. For each analysis we derived two tables, one representing the ontological categories over expressed in sensitive cell lines and second table with gene ontologies overregulated in resistant cell lines [18]. The results of the RankProd analysis are presented in Additional file 2: Figure S2.

DAVID analysis

We used online Database for Annotation, Visualization and Integrated Discovery (DAVID) software as described before [17,19]. We also generated clusters, which reduced overlaps and redundancies in regulated ontological categories, for Erlotinib-, irreversible inhibitor- and Cetuximab-sensitive and resistant cell lines. These are provided in Additional files 3,4,5,6,7, and 8.

Results

Datasets characterization

Inhibitors have been studied principally to find out their mechanism of action in different cancer types. We here aim to study the differential gene expression in EGFR inhibitor-sensitive vs. resistant cancer cell lines. We searched GEO Datasets using key term “EGFR Resistance” and, limiting our choices to human cell lines, selected studies that directly compare transcriptional profiles of matched EGFR inhibitor-sensitive vs. resistant datasets (Table 1). We found 8 appropriate studies comprising 15 datasets and 274 microarrays. In seven datasets EGFR inhibitor-sensitive vs. resistant cell lines were compared without any other treatment. In one study comprising four datasets (GSE38310), two different types of Erlotinib-resistant cell lines were originated from Erlotinib-sensitive cell lines and then treated either with DMSO or Erlotinib. We compared DMSO-treated sensitive cell lines with DMSO-treated resistant cell lines and Erlotinib-treated sensitive cell lines with Erlotinib-treated resistant cell lines. In a Gefitinib study comprising three datasets [10], sensitive and resistant cell lines were treated with Gefitinib and EGF separately or in combination (GSE34228). In one study Erlotinib-sensitive and resistant cell lines were treated with miR-7 (GSE40130). In another study, GSE38404, resistant cell lines were produced by exposing the sensitive cell lines to inhibitor for different time periods. We have combined data from all the resistant cell lines for metaanalysis.

Table 1 EGFR inhibitors-sensitive versus resistant cell lines datasets used in metaanalysis of microarrays

Illumina microarrays were used in majority of datasets, comprising 3 studies and 7 datasets. In four studies, comparing four datasets Affymetrix microarrays were used. Agilent microarrays were used in one study comparing four datasets. Different studies used different types of cell lines for example Non-Small Cell Lung Cancer, Head and Neck Cancer, A431 and Skin cancer cell lines; therefore, we have also individually compared datasets from individual studies (data not shown).

Information about the type of sensitive and resistant cells, pretreatment and type of EGFR inhibitor, labeling accession number, platform, number of datasets, number of sensitive and resistant microarray chips in each dataset is summarized in Table 1.

We have compared these datasets using different approaches; 1) datasets comparing all EGFR inhibitor-sensitive vs. resistant cell lines; 2) Gefitinib- and 3) Erlotinib-sensitive vs. resistant cell lines; 4) irreversible inhibitor-sensitive vs. resistant cell lines; 5) Cetuximab-sensitive vs. resistant cell lines [11]. We used RankProd software to select differentially expressed genes in these individual groups i.e. sensitive vs. resistant, with results graphically presented in (Additional file 2: Figure S2).

Ontological categories overrepresented in lists of differentially expressed genes were identified using DAVID [19]; complete lists are given in Additional files 3, 4, 5, 6 and 7. To condense redundancies, we have also identified clusters for differentially expressed ontological categories in the sensitive vs. Erlotinib-, Irreversible- and Cetuximab-resistant groups (see below).

Global comparison of all EGFR inhibitors-sensitive vs. resistant cell lines

We have compared overexpressed ontological categories in all sensitive vs. resistant cell lines. Table 2a contains the 10 categories with the best p-values; in Table 2b, we selected the characteristically different categories between sensitive and resistant cells with p-values better than 10−4.

Table 2 Global comparison of the ontological categories differentially expressed in EGFR inhibitors-sensitive vs. resistant cell lines

In Table 2a, we find that in the sensitive cell lines the genes related to systems involved in transport of materials across the membranes are overexpressed, whereas in the resistant cell lines the genes related to the metabolism and macromolecule biosynthesis, especially of proteins, are overexpressed. Although we have found genes associated with regulation of the protein translation in both groups, this category is more significant in the resistant cell lines. Both EGFR inhibitor-sensitive and resistant cell lines overexpress cell proliferation genes, however the sensitive cell lines tend to overexpress apoptosis and cancer-related genes as well. Interestingly, the resistant cell lines seem significantly more responsive to various endogenous and exogenous stimuli than the sensitive ones.

Oncogenesis-related genes and pathways are up-regulated in both sensitive and resistant cell lines. As reported in previous studies, in sensitive cell lines tumor growth seems dependent on the oncogene activation through ErbB receptor kinases [15]. In contrast, in resistant cell lines often the Ras pathway is activated independently of receptors (Table 2b). We note that while Ras is considered a downstream target of EGFR signaling in noncancerous cells [20], the EGFR- and Ras-associated genes comprise widely different groups (Table 2c).

Interestingly, we have observed that in the sensitive cell lines lipids are preferentially metabolized as the source of energy. But as cells become resistant, both carbohydrates and lipids are metabolized to provide energy. In addition, resistant cell lines are inducing the expression of energy generating genes, including the oxidoreductases. This observation suggests that energy production is an important matter for development of resistant cell lines.

Genes involved in responses to various extracellular and intracellular stimuli, for example steroid hormone, hydrogen peroxide and stress, are over expressed significantly in the resistant cell lines compared to sensitive ones. These responses may be related either to cells’ protection from toxins or overall cell survival, as stress and hydrogen peroxide responses prepare the cell for a toxic environment and steroids stimulate proliferation in certain cancers [21,22]. Those cancer cells that increase expression of genes responsive to these stimuli tend to adapt and survive [23].

We have also observed that certain processes related to immunity are overexpressed in the sensitive as well as in resistant cell lines, but we found that ontological categories for these processes have significantly better p-values in the resistant ones. Further, detailed study of the genes involved in immunity revealed that the innate immunity involving the complement component system is active in sensitive cells. In contrast, the resistant cells utilize both innate and adaptive immune systems, especially involving the cytokines and chemokines (Table 2d).

In EGFR inhibitor-sensitive cell lines, cell death and proliferation are relatively suppressed. But in resistant cell lines proliferation is relatively increased, while differentiation is decreased, thus favoring cancer cells growth and persistence. Another major distinction is that sensitive cell lines strongly express adhesion-related genes, whereas in resistant cells the genes related to cell movement are overexpressed. This suggests that the resistant cells have increased tendency to metastasize, perhaps via EMT [13,24].

Comparison of Erlotinib-sensitive vs. resistant cell lines

Surprisingly, protein biosynthesis was the only ontological category in the top ten differentially expressed ones in Erlotinib-sensitive cell lines (Table 3a). In the resistant cells, we have also found mitochondria, immunity and cytoskeleton genes in the top ten ontological categories. In Erlotinib-sensitive cell lines, carbohydrate and protein metabolism genes are overexpressed with very good p-values, 10−21 (Table 3b). These lines seem to utilize the maximum of their energy reservoirs from glucose and protein molecules. Conversely, in the Erlotinib-resistant cell lines, glycolysis and gluconeogenesis are suppressed and, importantly, genes related to mitochondria and mitochondrial processes are remarkably boosted. These results suggest that increased production of energy to support cellular metabolic processes and survival are very important in the development of Erlotinib resistance. From these observations we suggest that energy level could be the limiting factor for the tumor cells survival in the presence of EGFR inhibitors.

Table 3 Comparison of ontological categories differentially expressed in Erlotinib-sensitive vs. resistant cell lines, using DAVID program

Comparison of Gefitinib-sensitive vs. resistant cell lines

In the Table 4a, we show that the cancer-related genes and pathways are preferentially activated in the sensitive cell lines. But in the cells resistant to Gefitinib, the ontological categories related to macromolecule biosynthesis, specifically of proteins, are most significantly overexpressed. Apparently, the sensitive cell lines have higher tendency for attachment, as the adhesion related genes are upregulated in these, whereas the resistant cell lines are significantly more responsive to endogenous and exogenous stimuli. In the sensitive cell lines ontological categories related to cell cycle, metabolism of proteins and signaling are also up regulated, but the p-values are not as significant as in the resistant cell lines. Tumor related pathways and cell cycle-related genes are also relatively upregulated in resistant cell lines.

Table 4 Comparison of ontological categories differentially expressed in Gefitinib-sensitive vs . resistant cell lines

Interestingly, we have not found any ontological category related to immune system development overexpressed in the sensitive cell lines. In contrast, in the resistant cell lines immunity systems, specifically biosynthesis of cytokines genes, are boosted (Table 4b). This suggests that Gefitinib-resistant cell lines use cytokines of the immune systems to develop resistance against Gefitinib and thus maintain tumor cells growth and progression [25].

In resistant cells the ontological categories related to responses to various stimuli and angiogenesis were robustly and consistently overexpressed. We have also found increased expression of genes responsive to reactive oxygen species, which ultimately tends to activate apoptosis and immunity pathways in these cells. Expression of genes related to hypoxia has also been seen in the resistant cell lines, which may increase their drug resistance and ultimately their growth rate [26].

Both in Gefitinib-sensitive and in resistant cell lines proliferation and apoptosis are regulated. However, in the resistant cell lines, the cell death processes are negatively regulated and the cell division is enhanced, suggesting that resistant cell lines have higher propensity to survive. Meanwhile, responses to various endogenous and exogenous stimuli are enhanced, maybe to cope with stresses from environment. In addition, immunity genes and cytokine biosynthesis are upregulated to help cancer cells to survive and grow.

Interestingly, we have observed that lipid metabolism genes are prominent in the sensitive cell lines, whereas oxidoreductases, as well as both glucose and lipid metabolism, i.e. genes related to production of energy, are more prominent in the resistant cell lines. Energy requirements for the tumor cells are met preferentially by formation of new blood vessels but also by glucose and lipid metabolism.

Comparison of sensitive vs. resistant cell lines obtained using irreversible EGFR inhibitors

We have combined the microarray data from the cell lines selected as resistant to irreversible EGFR inhibitors WZ4002 and PF299804. Unexpectedly, in the top ten categories from the sensitive cell lines, genes related to membrane systems are found (Table 5a). In contrast, the extracellular region-related ontological categories are upregulated in the resistant cell lines. Cell cycle and lipid metabolism categories are relatively upregulated in the sensitive cell lines, similar to what was observed in reversible EGFR inhibitors-sensitive cell lines (Table 4b). Ontological categories related to cell division were found in both groups, however, in the resistant cell lines, cell division is positively regulated with better p-value than in sensitive cell lines. Importantly, the immunity genes are relatively overexpressed in resistant cell lines, reinforcing our observation that cytokine overexpression is one of the survival strategies gained by resistant cell lines (Table 5b).

Table 5 Comparison of ontological categories differentially expressed in irreversible inhibitor-sensitive vs . resistant cell lines

Surprisingly, we found increased expression of EGFR ligands in cells selected using irreversible inhibitors (Table 5b). Perhaps these cells are still dependent on EGFR-mediated signaling cascades. This differs from observations in cell lines selected for resistance to reversible inhibitors, which use alternative pathways for oncogenes activation (Tables 2b and 3b).

Comparison of Cetuximab-sensitive vs. resistant cell lines

For the analysis of cells selected as resistant to antibody inhibitor Cetuximab, we could find only a single study including six microarrays [27], which means that the statistical robustness of the results is reduced. Interestingly, we found that the differential gene expression in Cetuximab-sensitive and resistant cell lines is quite unlike the differential gene expression seen in the cell lines obtained using tyrosine kinase inhibitors. This agrees well with our previous study documenting that the transcriptional responses to EGFR antibody inhibitors are different from those to kinase inhibitors [9]. In sensitive cell lines, ontological categories related to immunity, cell proliferation, cell migration and response to external stimulus are overexpressed (Table 6a). In contrast, in the resistant cell lines epithelium development- and differentiation-associated gene ontologies are frequently and significantly overexpressed. Intriguing, the ontological category for wound response was overexpressed in resistant cells as well, although it had relatively less significant p-value.

Table 6 Comparison of ontological categories differentially expressed in Cetuximab-sensitive vs . resistant cell lines

In Table 6b, the overregulated processes in Cetuximab-sensitive cell line comprise immunity responses, signal transduction involving JAK-STAT pathway and angiogenesis. In contrast, in the resistant cell line genes for cell cycle, cell structure and reduction of transcription are increased.

Clustering of ontological categories in sensitive vs. resistant cell lines

The charts comparing overexpressed ontological categories, although very informative, contain many redundant and overlapping ontological categories (Tables 2,3,4,5 and 6). To get around this problem we have clustered the ontological categories using DAVID software. Cluster outputs for Erlotinib-, irreversible inhibitor and Cetuximab-sensitive versus resistant cell lines are presented in Additional file 8.

In Erlotinib-sensitive cell lines, transcription, translation and protein transport processes, marked as ‘T’ in Additional file 8a, are statistically more significant. In contrast, we found mitochondria-related genes, designated as ‘E’ for ‘energy’, the most frequent ontological category in resistant cell lines. Ontologies related to apoptosis, ‘A’, are the second most upregulated category in sensitive cells. In resistant cell lines genes related to cell cycle, ‘CC’, and cytoskeleton, ‘CY’, are the most expressed categories after energy production. In sensitive cell lines apparently energy is spent for synthesizing proteins and regulating cell death processes, whereas in resistant cells energy is mainly produced and utilized in cell growth maintenance. Thus, if this hypothesis is confirmed by the laboratory approaches, energy starvation strategies could be used to treat cancer cells.

Unexpectedly, categories ‘M’, describing motion of cells, and ‘V’, denoting the vascularization and angiogenesis, were found in sensitive cell lines, whereas adhesion-related genes, ‘AD’, were found relatively overexpressed in resistant cell lines. We have also found that certain categories related to energy metabolism; marked as ‘M-c’ for metabolism of carbohydrates and ‘E’ for energy-related processes are also up regulated in sensitive cells (Additional file 8a).

Although extracellular region ‘EC’ is the top regulated category in irreversible inhibitor-resistant cell lines, important EGFR-activated cellular processes, including immunity, ‘IM’, angiogenesis, ‘V’, and cytoskeleton, ‘CY’, were upregulated (Additional file 8b). In sensitive cell lines membrane systems, ‘MM’, response to hormones, ‘R’, and cell cycle, ‘CC’, were upregulated.

In Cetuximab antibody-resistant cell lines, many ontological categories are suppressed. For example, ontological category ‘IM’, immunity, is statistically the most prominent, biosynthesis of inflammatory cytokine and responses to these cytokines are specifically downregulated. We have observed that gene ontologies related to migration of the cells, marked ‘M’, are prominent in Additional file 8c. Cell adhesion categories, designated as ‘AD’, are also present there but the enrichment score and p-values are comparatively less. Blood vessels development processes, ‘V’, are overexpressed in Cetuximab-sensitive cell lines.

We have also found increased expression of signaling, ‘S’, apoptosis inhibition, ‘A’, and homeostasis-related genes, ‘H’, in sensitive cell lines, but not in resistant cell lines. Some of the overregulated clusters were also related to cell cycle, ‘CC’, transcription, ‘T’, and responses to external stimuli, ‘R’. Ectoderm development category was overexpressed with higher p-value in antibody-resistant cell lines, reflecting epidermal cell origin. Other ontological categories were extracellular matrix ‘EC’ and adhesion (Additional file 8c).

As observed above, (Tables 3,4,5 and 6), we have seen that responses of antibody-resistant cell line are reverse of those observed in tyrosine kinase inhibitors-resistant cell lines (comparing Additional files 8a and b with Table 6c). All ontological categories observed in Cetuximab-sensitive cell lines seem responsible for tumor cell survival and progression and were previously seen overexpressed in tyrosine kinase inhibitor-resistant cell lines.

Discussion

Global metaanalysis of EGFR inhibitor-resistant vs. sensitive cell lines presented here identifies the transcriptional and metabolic differences that allow tumor cells to evade EGFR-targeted therapies. Specifically, we find that the most acute problem created by EGFR inhibition is to provide a sufficient energy supply to the proliferating and metabolizing transformed cells. This metabolic predicament for tumors was one of the first ones identified, already by Warburg in 1927 [28]. The energy deficit caused by inhibition of EGFR is overcome by several mechanisms, e.g., by major enhancement of the mitochondrial oxidative phosphorylation in Erlotinib-inhibited cells, or by promoting vascularization and angiogenesis in Gefitinib-inhibited tumors.

In aerobic glycolysis, tumor cells produce lactate by glycolysis in the presence of oxygen and decrease the oxidative phosphorylation due to presence of an isoform of pyruvate kinase. If pyruvate kinase is in its native form, then increased oxygen uptake and oxidative phosphorylation occur and lactate formation is decreased [29]. In metaanalysis of Erlotinib-resistant cells we found that mitochondrial processes involving oxidative phosphorylation were upregulated. Moreover, the resistant cells recruit supplementary energy sources, inducing enzymes that catabolize carbohydrates in addition to those catabolizing lipids.

Additionally, our study proposes that EGFR tyrosine kinase inhibitors-resistant cancer cell lines express specific inflammatory cytokines and angiogenesis signals to promote vascularization and perhaps autocrine immunity stimulation as important strategy to combat cancer treatment [13,30]. In-depth perusal of the processes enhanced in the resistant cell lines revealed involvement of interleukins and chemokines with C-X-C motif. These observations are in line with previously known facts that cancer cells use immune cells and cytokines in various ways to maintain their own proliferation [25]. This mode of overcoming EGFR inhibition is particularly prominent in cells selected for Gefitinib resistance.

Perhaps expectedly, while the sensitive cell lines depend on the ErbB-dependent signaling, resistant ones circumvent this pathway and, at least in some cases, rely on the Ras pathway. Parenthetically, the mutations that confer EGFR inhibitor resistance generally occur in EGFR, seldom in Ras. In targeted therapy, as is case of EGFR inhibition, small molecules targeting tyrosine kinase domain of EGFR or monoclonal antibodies specifically blocking the ligand binding site are used. In some tumors, ErbB-independent signaling is sustained by activation of downstream molecules, for example Ras, PKI3CA and BRAF [15,21,31-34]. Correspondingly, we find that in EGFR inhibitor-resistant cell lines, Ras-dependent oncogenic pathways are preferentially expressed as a strategy to overcome the EGFR inhibition (Table 2b). Detailed analysis of the ErbB and Ras pathway genes enhanced in the sensitive and resistant cell lines respectively, revealed that genes involved in these ontological categories are quite dissimilar (Table 2d).

Interestingly, we found that irreversible inhibitors of the EGFR signaling pathway increased in resistant cell lines the expression of multiple EGFR ligands (Table 5b,c). This seems to be an important alternative to other modes of overcoming EGFR inhibition [11,35].

Hypoxia stimulates various processes in cancer cells including vascularization, growth factor signaling, genetic instability, cell survival, metastasis, cell death and antitumor drug resistance, in general favoring tumor survival and propagation [36]. The cells in surrounding areas of low oxygen are considered resistant to antitumor drugs mainly due to non-availability of drug, change in responses to drug, induction of drug-resistant genes, increase in mutations, angiogenesis and metastasis genes enhancement [26]. Metaanalysis comparing data from Gefitinib-sensitive versus resistant cell lines revealed that gene ontologies responsive to lower oxygen levels are specifically upregulated in resistant cell lines thus favoring survival. This finding further supports our observation that use of different inhibitors results in different, specific strategies to resist the antitumor targeted therapy.

An important caveat in this analysis stems from the fact that few of the studies, namely GSE34228 and GSE38310 for Gefitinib and Erlotinib, respectively, dominate the metaanalysis by overwhelming numbers of microarrays. Therefore, we cannot claim at this point that it is the resistance to Gefitinib that specifically produces the described changes in all cell lines selected for resistance to Gefitinib, and mutatis mutandis, the same for Erlotinib. However, we note that both GSE34228 and GSE38310 compared Non-Small Cell Lung Cancer lines, so the characteristic differences between Gefitinib- vs. Erlotinib-selected resistance cannot be ascribed to different cell types. Additional data are needed to corroborate or refute inhibitor-specific characteristics of resistant cells. That said, our previous metaanalysis demonstrated that different inhibitors, although working by ostensibly same mechanisms, cause identifiably different transcriptional responses, which allows for inhibitor-specific developments of resistance as well.

Conclusions

This metaanalysis of transcriptional and metabolic differences between EGFR inhibitor-resistant vs. sensitive cell lines identifies changes that allow tumor cells to evade EGFR-inhibition. The use of different inhibitors results in different, specific strategies to resist the antitumor therapy. Common pathways are upregulated in cell lines resistant to inhibitors targeting the kinase domain of EGFR, however, there are certain processes uniquely expressed against some of the inhibitors but, apparently, not others. The development of resistance to antibody inhibitors can vary significantly.

We found that the most acute problem created by EGFR inhibition is to provide a sufficient energy supply to the cells. The energy deficit can be overcome by several routes, e.g., by boosting mitochondrial oxidative phosphorylation, or by promoting vascularization and angiogenesis.

Further study of sensitive and resistant cancer cell lines responses to additional EGFR inhibitors will improve our understanding of drug resistance development and thus lead to improved anticancer treatment strategies. For example, use of mitochondrial blockers with Erlotinib, immunity blockers with Gefitinib, or combining tyrosine kinase inhibitors with antibody inhibitors, may avoid development of resistance to EGFR inhibitors.

References

  1. Vogelstein B, Kinzler KW. Cancer genes and the pathways they control. Nat Med. 2004;10(8):789–99.

    Article  CAS  PubMed  Google Scholar 

  2. Peto J, Houlston RS. Genetics and the common cancers. Eur J Cancer. 2001;37 Suppl 8:S88–96.

    Article  CAS  PubMed  Google Scholar 

  3. Vanneman M, Dranoff G. Combining immunotherapy and targeted therapies in cancer treatment. Nat Rev Cancer. 2012;12(4):237–51. doi:210.1038/nrc3237.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Hynes NE, Lane HA. ERBB receptors and cancer: the complexity of targeted inhibitors. Nat Rev Cancer. 2005;5(5):341–54.

    Article  CAS  PubMed  Google Scholar 

  5. Schlessinger J. Common and distinct elements in cellular signaling via EGF and FGF receptors. Science. 2004;306(5701):1506–7.

    Article  CAS  PubMed  Google Scholar 

  6. Blumenberg M. Profiling and metaanalysis of epidermal keratinocytes responses to epidermal growth factor. BMC Genomics. 2013;14:85. (doi):10.1186/1471-2164-1114-1185.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Wheeler DL, Dunn EF, Harari PM. Understanding resistance to EGFR inhibitors-impact on future treatment strategies. Nat Rev Clin Oncol. 2010;7(9):493–507. doi: 410.1038/nrclinonc.2010.1097. Epub 2010 Jun 1015.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Li SN, Li HQ. Epidermal growth factor receptor inhibitors: a patent review (2010 - present). Expert Opin Ther Pat. 2014;24(3):309–21. doi:310.1517/13543776.13542014.13871527. Epub 13542014 Jan 13543773.

    Article  CAS  PubMed  Google Scholar 

  9. Blumenberg M. Differential Transcriptional Effects of EGFR Inhibitors. PLoS One. 2014;9(9):e102466. LID - 102410.101371/journal.pone.0102466 [doi].

    Article  PubMed  PubMed Central  Google Scholar 

  10. Nakata A, Gotoh N. Recent understanding of the molecular mechanisms for the efficacy and resistance of EGF receptor-specific tyrosine kinase inhibitors in non-small cell lung cancer. Expert Opin Ther Targets. 2012;16(8):771–81. doi:710.1517/14728222.14722012.14697155. Epub 14722012 Jul 14728225.

    Article  CAS  PubMed  Google Scholar 

  11. Hatakeyama H, Cheng H, Wirth P, Counsell A, Marcrom SR, Wood CB, et al. Regulation of heparin-binding EGF-like growth factor by miR-212 and acquired cetuximab-resistance in head and neck squamous cell carcinoma. PLoS One. 2010;5(9):e12702. doi:12710.11371/journal.pone.0012702.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Cortot AB, Repellin CE, Shimamura T, Capelletti M, Zejnullahu K, Ercan D, et al. Resistance to irreversible EGF receptor tyrosine kinase inhibitors through a multistep mechanism involving the IGF1R pathway. Cancer Res. 2013;73(2):834–43. doi:810.1158/0008-5472.CAN-1112-2066. Epub 2012 Nov 1119.

    Article  CAS  PubMed  Google Scholar 

  13. Zhang Z, Lee JC, Lin L, Olivas V, Au V, LaFramboise T, et al. Activation of the AXL kinase causes resistance to EGFR-targeted therapy in lung cancer. Nat Genet. 2012;44(8):852–60. doi:810.1038/ng.2330.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Giles KM, Kalinowski FC, Candy PA, Epis MR, Zhang PM, Redfern AD, et al. Axl mediates acquired resistance of head and neck cancer cells to the epidermal growth factor receptor inhibitor erlotinib. Mol Cancer Ther. 2013;12(11):2541–58. doi:2510.1158/1535-7163.MCT-2513-0170. Epub 2013 Sep 2511.

    Article  CAS  PubMed  Google Scholar 

  15. Chong CR, Janne PA. The quest to overcome resistance to EGFR-targeted therapies in cancer. Nat Med. 2013;19(11):1389–400. doi:1310.1038/nm.3388. Epub 2013 Nov 1387.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Gautier L, Cope L, Bolstad BM. Irizarry RA: affy–analysis of Affymetrix GeneChip data at the probe level. Bioinformatics. 2004;20(3):307–15.

    Article  CAS  PubMed  Google Scholar 

  17. Mimoso C, Lee DD, Zavadil J, Tomic-Canic M, Blumenberg M. Analysis and meta-analysis of transcriptional profiling in human epidermis. Methods Mol Biol. 2014;1195:61–97. (doi):10.1007/7651_2013_1060.

    Article  CAS  PubMed  Google Scholar 

  18. Hong F, Breitling R, McEntee CW, Wittner BS, Nemhauser JL, Chory J, et al. RankProd: a bioconductor package for detecting differentially expressed genes in meta-analysis. Bioinformatics. 2006;22(22):2825–7. Epub 2006 Sep 2818.

    Article  CAS  PubMed  Google Scholar 

  19. Dennis Jr G, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, et al. DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biol. 2003;4(5):P3. Epub 2003 Apr 2003.

    Article  PubMed  Google Scholar 

  20. Ohashi K, Sequist LV, Arcila ME, Moran T, Chmielecki J, Lin YL, et al. Lung cancers with acquired resistance to EGFR inhibitors occasionally harbor BRAF gene mutations but lack mutations in KRAS, NRAS, or MEK1. Proc Natl Acad Sci U S A. 2012;109(31):E2127–33.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Fulda S, Gorman AM, Hori O, Samali A. Cellular stress responses: cell survival and cell death. Int J Cell Biol. 2010;2010:214074.

    PubMed  PubMed Central  Google Scholar 

  22. Ulrich CM, Bigler J, Potter JD. Non-steroidal anti-inflammatory drugs for cancer prevention: promise, perils and pharmacogenetics. Nat Rev Cancer. 2006;6(2):130–40.

    Article  CAS  PubMed  Google Scholar 

  23. Yap TA, Sandhu SK, Carden CP, de Bono JS. Poly (ADP-ribose) polymerase (PARP) inhibitors: Exploiting a synthetic lethal strategy in the clinic. CA Cancer J Clin. 2011;61(1):31–49.

    Article  PubMed  Google Scholar 

  24. Mimeault M, Batra SK. Functions of tumorigenic and migrating cancer progenitor cells in cancer progression and metastasis and their therapeutic implications. Cancer Metastasis Rev. 2007;26(1):203–14.

    Article  CAS  PubMed  Google Scholar 

  25. Dranoff G. Cytokines in cancer pathogenesis and cancer therapy. Nat Rev Cancer. 2004;4(1):11–22.

    Article  CAS  PubMed  Google Scholar 

  26. Brown JM, Wilson WR. Exploiting tumour hypoxia in cancer treatment. Nat Rev Cancer. 2004;4(6):437–47.

    Article  CAS  PubMed  Google Scholar 

  27. Boeckx C, Op de Beeck K, Wouters A, Deschoolmeester V, Limame R, Zwaenepoel K, et al. Overcoming cetuximab resistance in HNSCC: The role of AURKB and DUSP proteins. Cancer Lett. 2014;354(2):365–77.

    Article  CAS  PubMed  Google Scholar 

  28. Warburg O, Wind F, Negelein E. The metabolism of tumors in the body. J Gen Physiol. 1927;8(6):519–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Christofk HR, Vander Heiden MG, Harris MH, Ramanathan A, Gerszten RE, Wei R, et al. The M2 splice isoform of pyruvate kinase is important for cancer metabolism and tumour growth. Nature. 2008;452(7184):230–3.

    Article  CAS  PubMed  Google Scholar 

  30. Luo J, Solimini NL, Elledge SJ. Principles of cancer therapy: oncogene and non-oncogene addiction. Cell. 2009;136(5):823–37.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Diaz Jr LA, Williams RT, Wu J, Kinde I, Hecht JR, Berlin J, et al. The molecular evolution of acquired resistance to targeted EGFR blockade in colorectal cancers. Nature. 2012;486(7404):537–40.

    CAS  PubMed  PubMed Central  Google Scholar 

  32. Misale S, Yaeger R, Hobor S, Scala E, Janakiraman M, Liska D, et al. Emergence of KRAS mutations and acquired resistance to anti-EGFR therapy in colorectal cancer. Nature. 2012;486(7404):532–6.

    CAS  PubMed  PubMed Central  Google Scholar 

  33. Sartore-Bianchi A, Martini M, Molinari F, Veronese S, Nichelatti M, Artale S, et al. PIK3CA mutations in colorectal cancer are associated with clinical resistance to EGFR-targeted monoclonal antibodies. Cancer Res. 2009;69(5):1851–7.

    Article  CAS  PubMed  Google Scholar 

  34. Lievre A, Bachet JB, Le Corre D, Boige V, Landi B, Emile JF, et al. KRAS mutation status is predictive of response to cetuximab therapy in colorectal cancer. Cancer Res. 2006;66(8):3992–5.

    Article  CAS  PubMed  Google Scholar 

  35. Rajput A, Koterba AP, Kreisberg JI, Foster JM, Willson JK, Brattain MG. A novel mechanism of resistance to epidermal growth factor receptor antagonism in vivo. Cancer Res. 2007;67(2):665–73.

    Article  CAS  PubMed  Google Scholar 

  36. Harris AL. Hypoxia–a key regulatory factor in tumour growth. Nat Rev Cancer. 2002;2(1):38–47.

    Article  CAS  PubMed  Google Scholar 

Download references

Financial support

IRSIP, Higher Education Commission, Islamabad, Pakistan, and the R.O.Perelman Department of Dermatology, NYU Medical School, New York, USA.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Miroslav Blumenberg.

Additional information

Competing interests

The authors declared that they have no competing interests.

Authors’ contributions

SY carried out the data acquisition and analysis and drafted the manuscript. QJ participated in the study design and helped to draft the manuscript. MB conceived of the project, designed the study and wrote the manuscript. All authors read and approved the final manuscript.

Additional files

Additional file 1: Figure S1.

Flow chart of metaanalysis procedure.

Additional file 2: Figure S2.

Differentially expressed genes in sensitive cell lines vs. resistant cell lines. Top panels genes overexpressed in resistant cell lines, bottom panels genes overexpressed in sensitive cell lines. a) All studies. b) Gefitinib. c) Erlotinib. d) Irreversible inhibitors. e) Cetuximab.

Additional file 3:

Ontological categories characteristically expressed in EGFR inhibitors-sensitive vs. resistant cell lines.

Additional file 4:

Ontological categories characteristically expressed in Erlotinib-sensitive vs. resistant cell lines.

Additional file 5:

Ontological categories characteristically expressed in Gefitinib-sensitive vs. resistant cell lines.

Additional file 6:

Ontological categories characteristically expressed in irreversible inhibitors-sensitive vs. resistant cell lines.

Additional file 7:

Ontological categories characteristically expressed in Cetuximab-sensitive vs. resistant cell lines.

Additional file 8:

a. Cluster analysis for Erlotinib Sensitivevs.regulated analysis using DAVID software. For complete lists see Additional file 4. Additional files 8b and 8c Clusters of differentially expressed gene ontologies in irreversible inhibitor- and Cetuximab-sensitive vs. resistant cell lines, respectively. For complete lists see Additional files 6 and 7.

Rights and permissions

Open Access  This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.

The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.

To view a copy of this licence, visit https://creativecommons.org/licenses/by/4.0/.

The Creative Commons Public Domain Dedication waiver (https://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Younis, S., Javed, Q. & Blumenberg, M. Transcriptional changes associated with resistance to inhibitors of epidermal growth factor receptor revealed using metaanalysis. BMC Cancer 15, 369 (2015). https://doi.org/10.1186/s12885-015-1337-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12885-015-1337-3

Keywords