Next Article in Journal
Protein Quality Control Activation and Microtubule Remodeling in Hypertrophic Cardiomyopathy
Next Article in Special Issue
Application of Highly Immunocompromised Mice for the Establishment of Patient-Derived Xenograft (PDX) Models
Previous Article in Journal
Anti-Neuroinflammatory Effect of Alantolactone through the Suppression of the NF-κB and MAPK Signaling Pathways
Previous Article in Special Issue
Impact of Tumour Hypoxia on Evofosfamide Sensitivity in Head and Neck Squamous Cell Carcinoma Patient-Derived Xenograft Models
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Cells
Volume 8
Issue 7
10.3390/cells8070740
cells-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Induction of Acquired Resistance towards EGFR Inhibitor Gefitinib in a Patient-Derived Xenograft Model of Non-Small Cell Lung Cancer and Subsequent Molecular Characterization

by
Julia Schueler
1,*,
Cordula Tschuch
1,
Kerstin Klingner
1,
Daniel Bug
2,
Anne-Lise Peille
1,
Leanne de Koning
3,
Eva Oswald
1,
Hagen Klett
1 and
Wolfgang Sommergruber
4
1
Charles River Discovery Research Services Germany GmbH, Am Flughafen 12-14, 79108 Freiburg, Germany
2
Lehrstuhl für Bildverarbeitung, RWTH-Aachen University, 52074 Aachen, Germany
3
Institut Curie, PSL Research University, Department of Translational Research, 75005 Paris, France
4
University of Applied Sciences, FH Campus Wien, 1200 Vienna, Austria
*
Author to whom correspondence should be addressed.
Cells 2019, 8(7), 740; https://doi.org/10.3390/cells8070740
Submission received: 29 May 2019 / Revised: 4 July 2019 / Accepted: 17 July 2019 / Published: 18 July 2019
(This article belongs to the Special Issue Patient-Derived Cancer Models from Basic Study to Clinical Application)
Browse Figures
Review Reports Versions Notes

Abstract

:
In up to 30% of non-small cell lung cancer (NSCLC) patients, the oncogenic driver of tumor growth is a constitutively activated epidermal growth factor receptor (EGFR). Although these patients gain great benefit from treatment with EGFR tyrosine kinase inhibitors, the development of resistance is inevitable. To model the emergence of drug resistance, an EGFR-driven, patient-derived xenograft (PDX) NSCLC model was treated continuously with Gefitinib in vivo. Over a period of more than three months, three separate clones developed and were subsequently analyzed: Whole exome sequencing and reverse phase protein arrays (RPPAs) were performed to identify the mechanism of resistance. In total, 13 genes were identified, which were mutated in all three resistant lines. Amongst them the mutations in NOMO2, ARHGEF5 and SMTNL2 were predicted as deleterious. The 53 mutated genes specific for at least two of the resistant lines were mainly involved in cell cycle activities or the Fanconi anemia pathway. On a protein level, total EGFR, total Axl, phospho-NFκB, and phospho-Stat1 were upregulated. Stat1, Stat3, MEK1/2, and NFκB displayed enhanced activation in the resistant clones determined by the phosphorylated vs. total protein ratio. In summary, we developed an NSCLC PDX line modelling possible escape mechanism under EGFR treatment. We identified three genes that have not been described before to be involved in an acquired EGFR resistance. Further functional studies are needed to decipher the underlying pathway regulation.

1. Introduction

Lung cancer is still the leading cause of cancer related deaths worldwide, accounting for 1.76 million deaths in 2018 (http://www.webcitation.org/77V92AxNh). Non-small cell lung cancer (NSCLC) represents 85% of all lung cancer cases. Besides surgery, chemo-, radio-, and in recent years, immunotherapy-targeted approaches are one of many treatment options. By targeting cancer-specific molecules and signaling pathways, the systemic treatment options display markedly less nonspecific toxicities. In a high proportion (~30%) of non-small cell lung cancer (NSCLC) patients, the oncogenic driver of tumor growth is the constitutively activated epidermal growth factor receptor (EGFR, HER-1/ErbB1), a receptor tyrosine kinase of the ErbB family, which has been identified as an anticancer target [1]. Apart from EGFR-targeting monoclonal antibodies (mAbs), small-molecule tyrosine kinase inhibitors (TKIs) are in clinical use and approved for fist line therapy in NSCLC. In general, tyrosine kinases are central targets as they not only play an important role in the modulation of growth factor signaling in the initial phase of tumor establishment but also during other stages of tumorigenesis, such as escape from treatment or epithelial-mesenchymal transition (EMT), including prognostic power [1,2,3,4,5,6]. Despite the overall good response towards EGFR, TKI’s acquired resistance occurs in most of the patients after 9 to 14 months of therapy [7]. Numerous mechanisms of acquired resistance towards EGFR-TKIs are reported, ranging from secondary mutations in the EGFR gene to amplification or overexpression of molecules within the EGFR signaling cascade, like MET or hepatocyte growth factor (HGF) [7,8,9,10]. Beside the well-described T790M mutation, novel mutations in the EGFR gene were determined recently. On top of that, non-EGFR related mutations, for example, TP53, were described to be enriched in patients developing an acquired resistance towards EGFR TKIs [11]. To overcome those acquired resistances, preclinical tools to study their development as well as to test new drugs overcoming those mechanisms are urgently needed. In recent times, there have been multiple efforts in the industry as well as academia to establish large panels of well-characterized patient-derived xenograft (PDX) models covering a wide range of different tumor types. Indeed, these collections are becoming the preferred research tool to optimize the drug development process at multiple steps, in particular for target validation, pharmacology, and translational studies [12,13,14,15]. Currently, these collections represent the complexity of tumor heterogeneity and the molecular diversity of human cancers. In our facility, we established a panel of 85 NSCLC PDX models, representing the molecular landscape of NSCLC. In the present study, we derived three Gefitinib-resistant sublines of an NSCLC PDX model (LXFA 677) that was originally sensitive towards EGFR targeted treatment. The PDX model was derived from a patient who received first line Cisplatin therapy and was EGFR as well as KRAS wt, which classified it for second line EGFR treatment. The sublines were established by constant treatment with Gefitinib over a period of a minimum of three months. We characterized the emerging resistant sublines thoroughly on the genetic as well as protein level to decipher the biological difference among them and in comparison, to their parental line. These data led to a better understanding of the evolution of resistance under EGFR TKI treatment. Furthermore, the sublines will serve as research tools to develop next generation compounds, improving the life expectancy of NSCLC patients with acquired resistance.

2. Materials and Methods

2.1. PDX Establishment

This study was carried out in strict accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals of the Society of Laboratory Animals (GV SOLAS). All animal experiments were approved by the Committee on the Ethics of Animal Experiments of the regional council (Permit Numbers: G-09/58, G-13/13 and G13/43). After written informed consent, tumor tissue from NSCLC patients undergoing surgery was placed in a storage solution and transported within a few hours to Charles River. Incoming material of every donor patient received a chronological unique number, which was subsequently used to identify the corresponding PDX model. To facilitate the overview of the PDX models, each model name starts with a three to four letter code identifying the tumor type. For NSCLC, three different subtypes were defined: LXFA, lung cancer Freiburg adeno carcinoma; LXFE, lung cancer Freiburg epidermoid carcinoma; and LXFL, lung cancer Freiburg large cell carcinoma. Four to six week old female NMRI nu/nu mice (Charles River, Germany) placed under isoflurane anesthesia received tumor implants subcutaneously in both flanks. During the first passages, mice were monitored for tumor growth for up to 12 months. When stable tumor growth could be determined, mice were sacrificed and tumor material was implanted into new recipient mice. In addition, xenograft material was stored in liquid nitrogen for future implantation or fixed in formalin and stored in liquid nitrogen for subsequent analyses. A PDX was defined as established when stable growth over at least three passages and regrowth from cultures stored in liquid nitrogen could be observed. The percentage of tumor implants displaying stable growth (take rate) and passage time were recorded for every model and every individual passage. Tumor growth was determined by a two-dimensional measurement with calipers weekly or biweekly depending on the growth characteristics of the respective PDX model. Tumor volumes were calculated according to the following equation: Tumor Vol (mm3) = a (mm) × b2 (mm2) × 0.5, where “a” is the largest diameter and “b” is the perpendicular diameter of the tumor representing an idealized ellipsoid. Animals had to be sacrificed when tumor volume reached 1800 mm3. The NSCLC PDX model, LXFA 677, was derived from a biopsy from a 62 years old male suffering from non-pre-treated adenocarcinoma of the lung. The metadata of the donor patient are listed in Table 1. The tumor growth behavior of the established PDX model LXFA 677 is depicted in Supplemental Figure S1.

2.2. Treatment Experiments In Vivo

Implantation was performed similar to PDX establishment, except that animals received unilateral tumor implants. Animals and tumor implants were monitored daily until the maximum number of implants showed clear signs of beginning solid tumor growth. At randomization, the volume of growing tumors was initially determined. Animals bearing 50 to 250 mm3 tumors, preferably 80 to 200 mm3, were distributed into experimental groups, with comparable median and mean tumor volumes. The day of randomization was designated as day 0 of an experiment and was also the first day of dosing. Gefitinib (#G-4408-1G, LC Laboratories, USA) was applied once a day via oral gavage in doses ranging from 6.25 to 100 mg/kg/d. Details are depicted in Table 2. The relative volume of an individual tumor on day X (RTVx) was calculated by dividing the absolute volume (mm3) of the respective tumor on day X (Tx) by the absolute volume of the same tumor on the day of randomization, i.e., on day 0 (T0), multiplied by 100, as shown by the following equation:
RTVx [%] = Tx/T10∗100.
Tumor inhibition on a particular day (T/Cx) was calculated from the median RTV of a test group and the median RTV of a control group multiplied by 100, as shown by the following equation:
T/Cx [%] = median RTVx treated group/median RTVx control group∗100.
The minimum T/C (%) value recorded for a particular group during an experiment represented the maximum anti-tumor activity for the respective compound (=optimal T/C). All experiments were performed according to the relevant animal welfare guidelines published by FELASA and GV-SOLAS in an AAALAC accredited animal facility. Tumor volume and body weight were determined twice per week. A flow chart indicates the number of mice and passaging intervals during the generation of the three sublines (Supplemental Figure S2).

2.3. Formalin-Fixed Paraffin-Embedded (FFPE) Samples, Tissue Micro Array (TMA), and Immune-Histochemistry (IHC)

FFPE: Tumors were collected immediately after euthanasia of the donor animal. FFPE fixation was performed in 10% neutral buffered formalin for 24 h followed by routine processing and embedding into paraffin
TMA: Whole tumor sections (4 μm) were cut and stained with Hematoxylin-Eosin (H&E). H&E sections of the xenografts were studied by light microscopy and representative areas marked on the slides. Xenograft biopsies, 1 mm in diameter, were taken from the corresponding area in the paraffin block and arrayed in duplicates into a new recipient block as described.
IHC: After antigen retrieval, 5 μM FFPE tissue sections were incubated with anti-human EGFR Antibody (1:36; Dako Cat# M7239, Lot 20055023) overnight at 4 °C, followed by DAB staining and hematoxylin counterstaining. Image analysis: Digitalized images of the IHC slides were evaluated to determine the percentage of EGFR positive areas using OSANO software. A computerized analysis for digitized whole-slide images of the samples was used to quantify the EGFR expression using color classification and morphological image processing techniques.

2.4. RNA Isolation

Total RNA was isolated from the snap frozen tumor xenograft samples with the “mirVana miRNA Isolation kit” (Ambion, Carlsbad, CA, USA) according to the manufacturer’s instructions. The genomic DNA was digested using “RNase-free DNase Set” (Qiagen, Hilden, Germany). The quality and quantity of the RNA preparations were controlled using the Bioanalyzer (Agilent Technologies, Palo Alto, CA, USA) and the NanoDrop 2000 (Thermo Scientific, Waltham, MA, USA).

2.5. DNA Isolation

DNA was extracted from snap frozen tumor xenograft samples. Tumors were digested with proteinase K at 55 °C overnight and the lysates were digested with DNase-free RNase (Qiagen). DNA was extracted with phenol:chloroform:isoamylalcohol and precipitated with ethanol. DNA pellets were washed and resuspended in TElow (Tris 10 mM pH8, EDTA 0.1 mM pH8). The integrity of each DNA preparation was checked on a 1.3% agarose gel and the purity analyzed using the NanoDrop 2000 (Thermo Scientific, Waltham, MA, USA).

2.6. Reverse Transcription

RNA samples were reverse transcribed in a final volume of 20 µL (200 ng total RNA, 10 mmol/L DTT, 1 µg hexamer primers, 2 U MMLV Reverse Transcriptase (Invitrogen, Carlsbad, CA), 40 U of RNasin (Promega, Madison, WI), 0.5 mmol/L each dNTP (Promega), 1X reaction buffer). After cDNA synthesis, samples were adjusted to a final volume of 200 µL and stored at −20 °C until use. Then, 5 µL cDNA (corresponding to a quantity of 5 ng reverse-transcribed RNA) was used for each PCR.

2.7. qRT-PCR

qRT-PCR assays for human c-MET, ERBB2, ERBB3, HGF, PTEN, and AXL were developed and used on DNA and cDNA to quantify the gene copy number and the mRNA expression of these genes in the resistant LXFA 677 models, respectively. A qRT-PCR GAPDH assay and a qRT-PCR 18S ribosomal RNA assay were developed and used to detect these housekeeping genes on DNA and cDNA, respectively. qRT-PCR primers were designed with the online program Primer3, and the primer sequences are indicated in the table below (Table 3).
qRT-PCR was performed in 40 cycles on a StepOnePlus™ system (Applied Biosystems) using 2x SYBR Green I Master Mix (Kapa Biosystems). The c-MET, ERBB3, PTEN, HGF, and AXL data contained on cDNA was normalized using the 18S ribosomal RNA content of the sample. The c-MET and ERBB2 data obtained on DNA was normalized using the 18S ribosomal DNA content of the sample. The gene copy number and expression results are expressed in arbitrary units (AU) and were calculated as described below:
Gene of interest expression or copy number = 2 (Cthousekeeping gene − Ctgene of interest).
If no amplification was obtained or an unspecific product was amplified, the absence of a Ct value was replaced by 40, the maximum possible Ct value. In addition, the specificity of PCR products was assessed by melting curve analysis

2.8. Mutational Analysis (Whole Exome Sequencing)

After DNA extraction, material was sequenced with 126 bp paired-end reads using the ILLUMINA HiSeq-2500 platform and the Agilent V5 50MB enrichment kit, with a coverage of >160x Next, human reads were isolated using Xenome [16], aligned to the GRCh38 reference genome, and mutations were called using a workflow based on GATK (Genome Analysis Toolkit version 4) best practices. Using the variant effect predictor (VEP) [17], candidate mutations were annotated and filtered considering only variants with moderate or high protein impact and those being rare in healthy populations (<1% in gnomAD) [18]. Furthermore, mutations shared by all three resistant clones were annotated with protein functions using SIFT and Polyphen predictions from SNPnexus [19].

2.9. Reverse Phase Protein Array (RPPA)

In total, 150 xenograft samples (30 NSCLC models in five replicates) were provided to Institut Curie, France to perform a Reverse Phase Protein Array (RRPA). The samples were printed onto nitrocellulose covered microscope slides, using an Aushon 2470 precision printer. The samples were prepared as described in [20] and printed onto nitrocellulose covered slides (Supernova, Grace Biolabs) using a dedicated arrayer (2470 arrayer, Aushon Biosystems) in five serial dilutions (2000 to 125 μg/mL) and two replicates per dilution. Arrays were labeled with 33 specific antibodies (see Table 4). Read-out was done using IRDye 800CW (LiCOR) on an Innoscan 710-AL scanner (Innopsys). For the staining of total protein, arrays were incubated 30 min in Super G blocking buffer (Grace Biolabs), rinsed in water, incubated 5 min in 5.10−6% Fast green FCF (Sigma), and rinsed again in water. Raw data were normalized using Normacurve [21]. Next, each RPPA slide was median centered and scaled (divided by median absolute deviation). We then corrected for remaining sample loadings effects individually for each array by correcting the dependency of the data for individual arrays on the median value of each sample over all the arrays using a linear regression. Median normalized data were used to compare expression levels between groups of samples. The ratio between phosphorylated and total protein was calculated by calculating the difference between the log-transformed phospho-protein expression and the log-transformed total protein expression.

2.10. Protein–Protein Association Network Analysis

The Search Tool for Retrieval of Interacting Genes/Proteins (STRING, https://string-db.org/) was used to assign mutated genes as well as activated proteins to a protein interaction network. The enrichment analysis was based on gene ontology (GO) terms and KEGG pathways [22].

2.11. Statistical Analysis

Moderated t-tests [23] (followed by Benjamini–Hochberg multiple testing correction) were applied for statistical analysis of RPPA data. For the evaluation of the statistical significance in all other cases, the Mann–Whitney test (two-tailed) was used unless otherwise indicated. By convention, p-values ≤ 0.05 indicate significant differences. Statistical calculations were performed using GraphPad Prism bio-analytic software (version 6.02 for Windows, GraphPad Software, San Diego, CA, USA, www.graphpad.com) and R.

3. Results

3.1. Patient-Derived NSCLC Xenograft LXFA 677 Showed EGFR Dependency and Dose Dependent Sensitivity towards Gefitinib

To identify the best suited PDX model for this study, we screened our NSCLC PDX collection for EGFR protein expression by IHC. We therefore created a TMA comprising all 85 models in duplicates (Figure 1a and Supplemental Figure S3). To verify these results and to retrieve data for p-EGFR expression, we performed an RPPA with a subset of these models (n = 27) in triplicates (Figure 1b). The lung adenocarcinoma model LXFA 677 displayed high expression levels of EGFR and medium levels of p-EGFR as indicated by the median corrected normalized values (Figure 1a,b, black arrow). The EGFR and p-EGFR protein expression of this models was consistently high in both assays. To ensure that not only is the target expressed but is also druggable in vivo, we tested 47 NSCLC PDX towards their sensitivity against nine different EGFR targeting agents. Each model was tested against one to five compounds in total. LXFA 677 displayed a mean optimal T/C value of 23% (Supplemental Figure S4). In detail, it was sensitive towards treatment with four different EGFR inhibitors, TKIs, and mAb. The multikinase inhibitor, Sorafenib, targeting not exclusively EGFR but also other kinases, such as VEGFR and PDGFR, was only moderately active (Figure 1c). In a more detailed experimental setting, the model showed a dose-dependent sensitivity towards Gefitinib (Figure 1d). Taken together, LXFA 677 was sensitive towards EGFR inhibition, which underlined the essential role of the EGFR pathway, qualifying it as an ideal candidate for the experimental set-up (Figure 1d).

3.2. Gefitinib-Resistant Tumors Emerged under Continuous Treatment

To induce acquired resistance, we started treatment with a clinically relevant dose of 50 mg/kg/day and a suboptimal dose of 40 mg/kg/day per oral gavage, respectively (Figure 2a). Once the tumors went into remission, the dosing was reduced to 80% and 50% of the initial concentration. Marked reduction of the tumor volume occurred 7 to 14 days after start of treatment depending on the treatment scheme. When tumors started regrowth, the Gefitinib dosing was adapted accordingly with increasing amounts of compound up to 50 mg/kg/day. The emerging tumors were designated as individual sublines, LXFA 677res 1–3. Although the number of mice per treatment arm was small (two for the high dose and four for the low dose regimen), we found that starting with the higher dose level of 50 mg/kg/day induced resistance more consistently than the second approach, where 80% of the clinically relevant dose was applied to the treatment-naïve tumors. In the latter, only one out of four tumors showed regrowth under Gefitinib treatments whereas two out of two tumors emerged as EGFR inhibition in the first experimental set-up.
When tumors reached a volume of ≥1000 mm3, the tissue was harvested and re-implanted in four recipient mice per donor tumor. Treatment started when subcutaneously implanted tumors showed signs of successful engraftment (tumor volume ≥50 mm3). The four mice bearing one subline were divided into one treatment arm receiving different doses of Gefitinib and one control arm receiving the control vehicle. All mice allocated to the Gefitinib groups received 60% of the clinically relevant dose (=30 mg/kg/day). Once a tumor volume of >200 mm3 was reached, the daily Gefitinib dose was enhanced to 40 mg/kg/day. This treatment regimen was maintained until individual tumors reached at least 1500 mm3. Subsequently, tumor tissue was harvested and viably frozen down for further use (Figure 2b). LXFA 677res1 showed similar or even faster growth kinetics in treated tumors compared to untreated tumors. LXFA 677res2 depicted the same tumor growth rate in all mice irrespective of their medication. Both sublines were developed with the high dose induction protocol. In contrast, LXFA 677res3 derived from the low dose induction protocol showed inconsistent tumor growth in the different settings. One tumor from the control arm did not grow at all, whereas the other one depicted a growth kinetic like the two tumors under continuous Gefitinib treatment. The doubling time of the different sublines are shown in Table 5. LXFA 677res1 grows faster than the original tumor models whereas the other lines indicate a similar growth behavior. Due to the small number of animals per group, statistical analysis was not performed. To examine if the acquired resistance was reversible, we transplanted tumor material from one animal of the control arms each into four recipient mice. These mice were then stratified into the control and treatment arm and the tumor volume was monitored over time. Despite the fact that the donor tumor material was without Gefitinib treatment for one passage, the re-implanted tumors showed resistance towards Gefitinib to a similar extent as the original resistant sublines, suggesting that the acquired resistance was not reversible (Supplemental Figure S5). In a first attempt to decipher the mechanism behind the resistance, the copy number variation and mRNA expression for a subset of genes was determined (Figure 2c). The copy number of cMET was enhanced in all three sublines, whereas the copy number of HER2 was increased in sublines LXFA 677res1 and -res3 and reduced in -res2. Consistently, the mRNA expression levels of all three sublines were enhanced compared to the original model (relative expression to LXFA 677 = 1). The same observation was made for AXL, but to a much lower extent. HER3, HGF, and pTEN were only slightly modulated in the different sublines with a trend towards downregulation.

3.3. Whole Exome Sequencing Analysis Revealed High Similarities between the Resistant and the Treatment Naïve Clones of NSCLC PDX LXFA 677

To reveal the possible mechanism(s) of resistance, whole exome sequencing analysis was performed with the three resistant sub-clones of the NSCLC PDX and compared to the respective data set of the original tumor model. In all four lines, no mutations, in particular the T790M mutations, were detected in the EGFR gene. Besides that, they shared 85% of their mutated genes (=504) in total (Figure 3a). The total number of mutated genes was also very similar among all lines, displaying only minor differences (Supplemental Table S1). A smaller subset of mutations was shared by at least two of the lines or was even unique for one specific subline (= not common; 26 to 42 mutated genes shared by at least two lines; 26 to 53 unique mutated genes). Hierarchical clustering of genes with differing but not unique mutated genes (at least2; n = 53) showed that the original treatment naïve clone clustered separately from the three other lines (Figure 3b). Within the resistant clones, LXFA 677res2 and LXFA 677res3 clustered together. The mutations shared by all three resistant clones (14 mutations in 13 genes) were annotated with protein functions using SIFT and Polyphen predictions from SNPnexus. Functional annotation was available for missense variants (n = 10), whereas it was missing for insertions and a missense variant on the X chromosome (n = 4). Of note, seven of these mutations were also identified in an NSCLC PDX from a patient with an acquired resistance towards Erlotinib (LXFE 2478). The missense mutation found in NOMO2, ARHGEF5, and SMTNL2 were predicted as deleterious based on the SIFT prediction. Moreover, NOMO2 was predicted as being possibly damaging by the Polyphen algorithm. Using the STRING platform, we investigated the mutated genes that got acquired by at least two resistant models in comparison to those in the naïve treatment (Supplemental Figure S6a). Tolerating a maximum of five interactors and a medium confidence level (=0.400), the Fanconi anemia pathway (KEGG), cell cycle checkpoints (reactome), and DNA topoisomerase type I activity (GO) showed up as the most important interactors with the lowest false discovery rate and the highest observed gene count (Supplemental Table S2).

3.4. The NSCLC PDX Lines Emerging under Gefitinib Treatment Displayed Specific Proteomic Profiles

To elucidate how the tumor cells developed a resistance mechanism against Gefitinib, a Reverse Phase Protein Array (RPPA) was performed on the three resistant sublines and the treatment naïve parental clone. Hereby, it was possible to study protein expression and activation (=phosphorylation) in a large panel of samples and proteins in a quantitative manner. The cluster analysis exhibited three distinct groups of proteins: EGFR and AXL showed high expression levels in all four lines. Medium expression levels were determined for phospho−p38 MAPK, phospho−p44/42 MAPK, phospho−PTEN, STAT3, Akt, STAT1, pTEN, NFkB, IGF−I receptor B, MEK1/2, mTOR, RelB, MET, HER3, and phospho−eIF4B. Low expression levels were detected for phospho−mTOR, phospho−Akt, phospho−c−Jun, phospho−STAT1, p44/42 MAPK, p38 MAPK, phospho−Her3, c−Jun, phospho−MEK1/2, eIF4B, IkB alpha, phospho−NFkB, HER2, phospho−EGFR, phospho−STAT3, and phospho−MET (Figure 4a). As the three resistant sublines displayed high similarities, we determined the logarithmic fold change (log FC) by comparing the values of all three resistant lines vs. the original treatment naïve line. Four proteins were upregulated in the treatment resistant lines with a logFC > 0.5. As the sample number was relatively low, n = 5 and n = 15, respectively, only phospho-STAT1 revealed a statistically significant difference (moderated t-test; multiple testing adjusted p-value = 0.03). A functional interaction analysis using the STRING platform was performed investigating the interaction of the upregulated proteins (Supplemental Figure S6B). The most prominently involved pathways for the upregulated proteins were the viral process and ERBB2 signaling (GO; Supplemental Table S3) and EGFR tyrosine kinase inhibitor resistance as well as prostate cancer (KEGG; Supplemental Table S4). The protein data corroborated the transcriptomic results, with cMET being the only exception. While cMET showed enhanced levels of copy numbers as well as an increase in mRNA levels, it did not translate to a higher protein level compared to the original PDX model.
A better understanding of the actual activity of the expressed proteins was achieved by comparing the ratio between phosphorylated and total protein (n = 14). The comparison of the ratio between phosphorylated vs. total protein displayed a moderate, albeit not significant, upregulation in the resistant clones for four proteins: MEK1/2, STAT3, STAT1, and NFkB (log FC > 0.5). In contrast, Akt was down-regulated in its activation status (logFC < −0.5, Figure 4b). The STRING platform was used again and potential interaction partners of the up-regulated activated proteins were investigated (Supplemental Figure S6c). The enrichment analysis revealed a highly significant positive correlation with the regulation of signal transduction (GO-term) and EGFR TKI resistance, different solid cancer types, namely pancreatic, prostate, and NSCLC as well as ErbB and HIF-1 signaling (KEGG pathway, Supplemental Tables S6 and S7).

4. Discussion

Acquired resistance prevents patients from long term benefits of targeted therapies in many different cancer types. This specifically holds true for anti-EGFR treatment in NSCLC patients, where multiple generations of TKIs all started displaying a loss of activity after a certain period of time [24]. It is therefore even more important to create preclinical models mimicking this phenomenon as close as possible. Studies using tissue culture induced resistance do not take into account that it is not only the tumor cell involved in the development of the resistance mechanism. The complex cross-talk between tumor cells and their micro-environment does not only influence the phenotype of the tumor but its sensitivity towards anticancer agents as well [25]. These aspects can, at least for now, only be mimicked in an animal model. The present study used in vivo PDX of NSCLC to model the onset of resistance against EGFR TKI treatment. A PDX model does recapitulate some aspects of the patient tumor more closely than conventional cell lines, which makes this model an ideal candidate for such a complex task [13].
The PDX model selected in this study represents a typical late stage NSCLC. The fact that the original tumor was EGFR and KRAS wt classified it for second line EGFR treatment [26]. Indeed, it mirrored precisely the clinical situation, with an emerging resistance under constant treatment with Gefitinib.
A preliminary screening for five genes based on the copy number and qPCR led to first conclusions but did not reveal the complete mechanism behind the acquired resistance. Despite the upregulation of cMET on the RNA level, all other gene expression modulations were reflected as well on a protein level. A possible explanation could be the subtle down-regulation of activated eIF4B (logFC < −0.3) being involved in the translation of MET mRNA [27]. Based on the performed whole exome sequence analysis, we identified a set of 13 genes with shared mutations in all three Gefitinib-resistant LXFA 677 subclones (res1 to res3). Interestingly, none of these have yet been directly linked with resistance, except for NOMO2 (NODAL Modulator 2), the mutation which was identified in the current study as prognostically unfavorable by two algorithms (SIFT and Polyphen). NOMO2 together with three other genes has recently linked to radio-resistance and suggested as diagnostic biomarkers in radio-resistant human H460 lung cancer stem-like cells [28]. Functionally, NOMO2 is reported to block Nodal signaling [29]. Furthermore, high levels of Nodal are bad prognostic markers in melanoma, breast, and pancreatic cancer [30]. Thus, it can be hypothesized that the knock-down of the blocking signal from NOMO2 induces an increase of Nodal in the resistant clones as a possible escape mechanism. Of course, this hypothesis maintains to be verified. The two other mutations predicted as deleterious, ARHGEF5 and SMTNL2, are linked to a poor prognosis in lung cancer patients, epithelial-mesenchymal transition (EMT), and Rho-GTPase signaling [31,32]. The biological functions of some identified genes are still not elucidated, which holds true for CCDC74A (coiled-coil domain containing 74A) and OR2T35 (olfactory receptor family 2 subfamily T member 35), for which no publication is available as of Apri 2019. Two other genes, ZNF417 and KDM6B, have been associated as part of a gene signature with cancer aggressiveness [32,33,34,35,36,37,38,39,40]. In a gene enrichment analysis, the Fanconi anemia (FA) pathway was indicated in addition to the above mentioned biological processes. Interestingly, this pathway was reported in connection with EGFR resistance by other groups [41,42,43,44]. Pfaffle et al. explicitly described an FA phenotype in EGFR resistant cells and a possible susceptibility of these cells towards PARP inhibition [43]. This assumption was based inter alia on the work of Garcia-Higuera et al., who described the co-localization of proteins of the FA pathway with BRCA1 [45]. A future study investigating the sensitivity towards PARP inhibitors in the original as well as the resistant clones will help to further characterize this possible link and treatment option.
Although the above discussed genes were found to be mutated in all three Gefitinib resistant clones, it remains to be determined if and how these mutations may contribute to resistance. For most of those genes, only a very vague knowledge about their biologic function in general exists. However, the comparison with the mutational landscape of another NSCLC PDX of our collection, which was derived from a patient with an acquired resistance to Erlotinib, revealed and overlap of 7 out of 14 mutations (in 13 genes) with the current PDX models, where the resistance was acquired in the tumor bearing mice.
The possible trend towards EMT seen in the whole exome sequencing analyses was substantiated by the fact that NFκB activation was enhanced in the resistant clones as determined by RPPA [46]. A biological link between EGFR expression and NFκB activation was described for different solid cancer types, like NSCLC, but also for renal cancer [47,48]. As our PDX model is carrying a p53 mutation, it thereby qualifies as a possible pre-clinical model for the development of NFκB inhibitors. The fact that IκBα was slightly enhanced (logFC > 0.1) in the resistant clones, indicates that the classical NFκB activation pathways seemed to be involved [49]. The activation of EGFR via different upstream proteins, like RET [50] or HER2 [51], was described before. In the present case, the exact activating mechanism of EGFR remains to be elucidated. In any case, when EGFR is overexpressed, it becomes phosphorylated on tyrosine residues located in the C-terminal part and constitutively activated, forming a complex with TBK1 and IRF3 [52]. Additionally, it has been shown that STAT1 gene expression is upregulated by nuclear EGFR and HER2 via cooperation with STAT3 [53]. This observation is supported by our data. The upregulation of AXL determined in all three resistant lines of LXFA 677 was observed by other groups in vitro as well as in patient samples of NSCLC and other solid cancer types [54,55,56,57]. Tian et al. described the possibility to restore sensitivity towards Gefitinib in NSCLC cell lines by silencing AXL [58]. The observed upregulation of activated MEK1/2 is in line with the comments of other groups that MEK1/2 inhibition can prevent or circumvent EGFR TKI resistance in monotherapy or in combination with PI3K/mTOR inhibitors [59]. The MEK1/2 inhibition prevented the EMT emerging in the EGFR resistant NSCLC cells [60]. As a marked proportion of our data is suggesting EMT as one possible escape mechanism in the current model, an upregulation of MEK1/2 is in line with the published data. In contrast, the activation of STAT1 as well as of STAT3 was up to now, not linked to EGFR resistance. Interestingly, the activation of STAT3 is linked to EMT whereas the role of STAT1 is discussed controversially [61,62,63]. Apart from this, specifically, the upregulation of pSTAT1 under treatment with EGFRi correlated with worse survival in patients undergoing combination therapy with anti-EGFR and Cisplatin [64]. This observation was in line with the trend towards a poor prognosis and aggressiveness of tumor growth elucidated in the whole exome sequencing analysis. Cheng et al. described a correlation between a reduction of PDL-1 under EGFRi treatment via STAT1. The upregulation of activated STAT1 in the acquired resistance model can be interpreted as a possible immune escape mechanism of the tumor cells to upregulate PDL-1 and evade the CD8+ tumor infiltrating lymphocytes [65]. The downregulation of activated Akt has not been described yet for acquired resistance in EGFR TKI treatment. Nevertheless, Akt downregulation via inactivation of Ets1- function was described by two independent groups as an innate drug resistance mechanism in EGFR mutated NSCLC cells [66,67]. The possibility that a similar escape mechanism can evolve under treatment with EGFR TKIs in EGFR wt NSCLC lines needs to be investigated. In order to obtain a more robust understanding about the complex escape mechanisms, it will be of great importance to generate and analyze additional resistant clones with regard to their mutation status, their transcriptional profile, and their phospho-proteomics status.
In summary, the generation of sister tumor cell lines derived from a pre- and post-acquired resistance status makes this PDX panel a valuable tool for the exploration of the biological mechanism behind this phenomenon. The further analysis and utilization of this PDX model may facilitate research and development regarding the prevention and treatment of drug resistance and will certainly contribute towards continued improved treatment and long-term survival in NSCLC patients.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4409/8/7/740/s1.

Author Contributions

Conceptualization, J.S., C.T., H.K., and W.S.; methodology, C.T., K.K., L.K., A.-L.P., and E.O.; software, D.B. and H.K.; validation, J.S., E.O., L.K., H.K., and W.S.; formal analysis, J.S., L.K., and H.K.; investigation, C.T., K.K., E.O., and A.-L.P; data curation, C.T., K.K., A-L.P., and L.K.; writing—original draft preparation, J.S.; writing—review and editing, H.K., L.K., and W.S.; visualization, J.S., W.S., and H.K.; supervision, W.S.; project administration, J.S.; funding acquisition, J.S and W.S.

Funding

This work has been carried out as a part of the IMI-funded EU-consortium PREDECT (www.predect.eu).

Acknowledgments

We thank Anke Behnke, Stefanie Klingel, Volker Knauff, Dorothée Lenhardt, Bérengère Ouine, Yana Raeva, and Sabine Rajkumar for excellent technical support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wang, Z. ErbB Receptors and Cancer. Methods Mol. Biol. 2017, 1652, 3–35. [Google Scholar] [PubMed]
  2. Antony, J.; Thiery, J.P.; Huang, R.Y. Epithelial-to-mesenchymal transition: Lessons from development, insights into cancer and the potential of EMT-subtype based therapeutic intervention. Phys. Biol. 2019. [Google Scholar] [CrossRef] [PubMed]
  3. Siu, M.K.; Abou-Kheir, W.; Yin, J.J.; Chang, Y.S.; Barrett, B.; Suau, F.; Casey, O.; Chen, W.Y.; Fang, L.; Hynes, P.; et al. Loss of EGFR signaling regulated miR-203 promotes prostate cancer bone metastasis and tyrosine kinase inhibitors resistance. Oncotarget 2014, 5, 3770–3784. [Google Scholar] [CrossRef] [PubMed]
  4. Tulchinsky, E.; Demidov, O.; Kriajevska, M.; Barlev, N.A.; Imyanitov, E. EMT: A mechanism for escape from EGFR-targeted therapy in lung cancer. Biochim. Et Biophys. Acta 2019, 1871, 29–39. [Google Scholar] [CrossRef] [PubMed]
  5. Xu, J.; Wang, J.; Zhang, S. Mechanisms of resistance to irreversible epidermal growth factor receptor tyrosine kinase inhibitors and therapeutic strategies in non-small cell lung cancer. Oncotarget 2017, 8, 90557–90578. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Zhu, X.; Bao, Y.; Guo, Y.; Yang, W. Proline-Rich Protein Tyrosine Kinase 2 in Inflammation and Cancer. Cancers 2018, 10, 139. [Google Scholar] [CrossRef]
  7. Takeda, M.; Nakagawa, K. First- and Second-Generation EGFR-TKIs Are All Replaced to Osimertinib in Chemo-Naive EGFR Mutation-Positive Non-Small Cell Lung Cancer? Int. J. Mol. Sci. 2019, 20, 146. [Google Scholar] [CrossRef]
  8. Kobayashi, S.; Boggon, T.J.; Dayaram, T.; Jänne, P.A.; Kocher, O.; Meyerson, M.; Johnson, B.E.; Eck, M.J.; Tenen, D.G.; Halmos, B. EGFR Mutation and Resistance of Non–Small-Cell Lung Cancer to Gefitinib. N. Engl. J. Med. 2005, 352, 786–792. [Google Scholar] [CrossRef]
  9. Ohashi, K.; Maruvka, Y.E.; Michor, F.; Pao, W. Epidermal Growth Factor Receptor Tyrosine Kinase Inhibitor–Resistant Disease. J. Clin. Oncol. 2013, 31, 1070–1080. [Google Scholar] [CrossRef]
  10. Sequist, L.V.; Waltman, B.A.; Dias-Santagata, D.; Digumarthy, S.; Turke, A.B.; Fidias, P.; Bergethon, K.; Shaw, A.T.; Gettinger, S.; Cosper, A.K.; et al. Genotypic and Histological Evolution of Lung Cancers Acquiring Resistance to EGFR Inhibitors. Sci. Transl. Med. 2011, 3, ra26–ra75. [Google Scholar] [CrossRef]
  11. Lee, C.K.; Kim, S.; Lee, J.S.; Lee, J.E.; Kim, S.M.; Yang, I.S.; Kim, H.R.; Lee, J.H.; Kim, S.; Cho, B.C. Next-generation sequencing reveals novel resistance mechanisms and molecular heterogeneity in EGFR-mutant non-small cell lung cancer with acquired resistance to EGFR-TKIs. Lung Cancer 2017, 113, 106–114. [Google Scholar] [CrossRef] [PubMed]
  12. Herter-Sprie, G.S.; Kung, A.L.; Wong, K.K. New cast for a new era: Preclinical cancer drug development revisited. J. Clin. Investig. 2013, 123, 3639–3645. [Google Scholar] [CrossRef] [PubMed]
  13. Hidalgo, M.; Amant, F.; Biankin, A.V.; Budinska, E.; Byrne, A.T.; Caldas, C.; Clarke, R.B.; de Jong, S.; Jonkers, J.; Maelandsmo, G.M.; et al. Patient-derived xenograft models: An emerging platform for translational cancer research. Cancer Discov. 2014, 4, 998–1013. [Google Scholar] [CrossRef] [PubMed]
  14. Teicher, B.A. Tumor models for efficacy determination. Mol. Cancer Ther. 2006, 5, 2435–2443. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Tentler, J.J.; Tan, A.C.; Weekes, C.D.; Jimeno, A.; Leong, S.; Pitts, T.M.; Arcaroli, J.J.; Messersmith, W.A.; Eckhardt, S.G. Patient-derived tumour xenografts as models for oncology drug development. Nat. Reviews. Clin. Oncol. 2012, 9, 338–350. [Google Scholar] [CrossRef] [Green Version]
  16. Conway, T.; Wazny, J.; Bromage, A.; Tymms, M.; Sooraj, D.; Williams, E.D.; Beresford-Smith, B. Xenome—a tool for classifying reads from xenograft samples. Bioinformatics 2012, 28, i172–i178. [Google Scholar] [CrossRef] [PubMed]
  17. McLaren, W.; Gil, L.; Hunt, S.E.; Riat, H.S.; Ritchie, G.R.S.; Thormann, A.; Flicek, P.; Cunningham, F. The Ensembl Variant Effect Predictor. Genome Biol. 2016, 17, 122. [Google Scholar] [CrossRef]
  18. Sherry, S.T.; Ward, M.H.; Kholodov, M.; Baker, J.; Phan, L.; Smigielski, E.M.; Sirotkin, K. dbSNP: The NCBI database of genetic variation. Nucleic Acids Res. 2001, 29, 308–311. [Google Scholar] [CrossRef]
  19. Dayem Ullah, A.Z.; Oscanoa, J.; Wang, J.; Nagano, A.; Lemoine, N.R.; Chelala, C. SNPnexus: Assessing the functional relevance of genetic variation to facilitate the promise of precision medicine. Nucleic Acids Res. 2018, 46, W109–W113. [Google Scholar] [CrossRef]
  20. Meseure, D.; Vacher, S.; Lallemand, F.; Alsibai, K.D.; Hatem, R.; Chemlali, W.; Nicolas, A.; De Koning, L.; Pasmant, E.; Callens, C.; et al. Prognostic value of a newly identified MALAT1 alternatively spliced transcript in breast cancer. Br. J. Cancer 2016, 114, 1395–1404. [Google Scholar] [CrossRef] [Green Version]
  21. Troncale, S.; Barbet, A.; Coulibaly, L.; Henry, E.; He, B.; Barillot, E.; Dubois, T.; Hupe, P.; de Koning, L. NormaCurve: A SuperCurve-based method that simultaneously quantifies and normalizes reverse phase protein array data. PLoS ONE 2012, 7, e38686. [Google Scholar] [CrossRef]
  22. Szklarczyk, D.; Gable, A.L.; Lyon, D.; Junge, A.; Wyder, S.; Huerta-Cepas, J.; Simonovic, M.; Doncheva, N.T.; Morris, J.H.; Bork, P.; et al. STRING v11: Protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019, 47, D607–D613. [Google Scholar] [CrossRef]
  23. Ritchie, M.E.; Phipson, B.; Wu, D.; Hu, Y.; Law, C.W.; Shi, W.; Smyth, G.K. Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015, 43, e47. [Google Scholar] [CrossRef]
  24. Morgillo, F.; Della Corte, C.M.; Fasano, M.; Ciardiello, F. Mechanisms of resistance to EGFR-targeted drugs: Lung cancer. ESMO Open 2016, 1, e000060. [Google Scholar] [CrossRef]
  25. Sebens, S.; Schafer, H. The tumor stroma as mediator of drug resistance—A potential target to improve cancer therapy? Curr. Pharm. Biotechnol. 2012, 13, 2259–2272. [Google Scholar] [CrossRef]
  26. Campos-Parra, A.D.; Zuloaga, C.; Manriquez, M.E.; Aviles, A.; Borbolla-Escoboza, J.; Cardona, A.; Meneses, A.; Arrieta, O. KRAS mutation as the biomarker of response to chemotherapy and EGFR-TKIs in patients with advanced non-small cell lung cancer: Clues for its potential use in second-line therapy decision making. Am. J. Clin. Oncol. 2015, 38, 33–40. [Google Scholar] [CrossRef]
  27. Zhang, J.; Babic, A. Regulation of the MET oncogene: Molecular mechanisms. Carcinogenesis 2016, 37, 345–355. [Google Scholar] [CrossRef]
  28. Yun, H.S.; Baek, J.H.; Yim, J.H.; Um, H.D.; Park, J.K.; Song, J.Y.; Park, I.C.; Kim, J.S.; Lee, S.J.; Lee, C.W.; et al. Radiotherapy diagnostic biomarkers in radioresistant human H460 lung cancer stem-like cells. Cancer Biol. Ther. 2016, 17, 208–218. [Google Scholar] [CrossRef] [Green Version]
  29. Haffner, C.; Frauli, M.; Topp, S.; Irmler, M.; Hofmann, K.; Regula, J.T.; Bally-Cuif, L.; Haass, C. Nicalin and its binding partner Nomo are novel Nodal signaling antagonists. EMBO J. 2004, 23, 3041–3050. [Google Scholar] [CrossRef] [Green Version]
  30. Kalyan, A.; Carneiro, B.A.; Chandra, S.; Kaplan, J.; Chae, Y.K.; Matsangou, M.; Hendrix, M.J.C.; Giles, F. Nodal Signaling as a Developmental Therapeutics Target in Oncology. Mol. Cancer Ther. 2017, 16, 787–792. [Google Scholar] [CrossRef] [Green Version]
  31. He, P.; Wu, W.; Wang, H.; Liao, K.; Zhang, W.; Xiong, G.; Wu, F.; Meng, G.; Yang, K. Co-expression of Rho guanine nucleotide exchange factor 5 and Src associates with poor prognosis of patients with resected non-small cell lung cancer. Oncol. Rep. 2013, 30, 2864–2870. [Google Scholar] [CrossRef]
  32. Komiya, Y.; Onodera, Y.; Kuroiwa, M.; Nomimura, S.; Kubo, Y.; Nam, J.M.; Kajiwara, K.; Nada, S.; Oneyama, C.; Sabe, H.; et al. The Rho guanine nucleotide exchange factor ARHGEF5 promotes tumor malignancy via epithelial-mesenchymal transition. Oncogenesis 2016, 5, e258. [Google Scholar] [CrossRef]
  33. DiBardino, D.M.; Rawson, D.W.; Saqi, A.; Heymann, J.J.; Pagan, C.A.; Bulman, W.A. Next-generation sequencing of non-small cell lung cancer using a customized, targeted sequencing panel: Emphasis on small biopsy and cytology. CytoJournal 2017, 14, 7. [Google Scholar] [CrossRef]
  34. Gadea, G.; Blangy, A. Dock-family exchange factors in cell migration and disease. Eur. J. Cell Biol. 2014, 93, 466–477. [Google Scholar] [CrossRef]
  35. Gordon, E.A.; Whisenant, T.C.; Zeller, M.; Kaake, R.M.; Gordon, W.M.; Krotee, P.; Patel, V.; Huang, L.; Baldi, P.; Bardwell, L. Combining docking site and phosphosite predictions to find new substrates: Identification of smoothelin-like-2 (SMTNL2) as a c-Jun N-terminal kinase (JNK) substrate. Cell. Signal. 2013, 25, 2518–2529. [Google Scholar] [CrossRef] [Green Version]
  36. Kim, B.H.; Shenoy, A.R.; Kumar, P.; Das, R.; Tiwari, S.; MacMicking, J.D. A family of IFN-gamma-inducible 65-kD GTPases protects against bacterial infection. Science 2011, 332, 717–721. [Google Scholar] [CrossRef]
  37. Miura, K. ERK2-binding domain is required for phosphorylation of EBITEIN1, a potential downstream interactor of ERK2. Biochem. Biophys. Res. Commun. 2008, 375, 367–371. [Google Scholar] [CrossRef]
  38. Silva, M.P.; Barros-Silva, J.D.; Vieira, J.; Lisboa, S.; Torres, L.; Correia, C.; Vieira-Coimbra, M.; Martins, A.T.; Jeronimo, C.; Henrique, R.; et al. NCOA2 is a candidate target gene of 8q gain associated with clinically aggressive prostate cancer. Geneschromosomes Cancer 2016, 55, 365–374. [Google Scholar] [CrossRef]
  39. Xiao, Z.G.; Shen, J.; Zhang, L.; Li, L.F.; Li, M.X.; Hu, W.; Li, Z.J.; Cho, C.H. The Roles of Histone Demethylase UTX and JMJD3 (KDM6B) in Cancers: Current Progress and Future Perspectives. Curr. Med. Chem. 2016, 23, 3687–3696. [Google Scholar] [CrossRef]
  40. Xu, G.; Zhang, M.; Zhu, H.; Xu, J. A 15-gene signature for prediction of colon cancer recurrence and prognosis based on SVM. Gene 2017, 604, 33–40. [Google Scholar] [CrossRef]
  41. Ali, G.; Bruno, R.; Poma, A.M.; Affinito, O.; Monticelli, A.; Piaggi, P.; Ricciardi, S.; Lucchi, M.; Melfi, F.; Chella, A.; et al. Whole transcriptome targeted gene quantification provides new insights on pulmonary sarcomatoid carcinomas. Sci. Rep. 2019, 9, 3536. [Google Scholar] [CrossRef] [Green Version]
  42. Cathcart-Rake, E.; Lopez-Chavez, A. Young male with fanconi anemia and EGFR-mutant non-small-cell lung cancer. J. Thorac. Oncol. 2014, 9, e83–e85. [Google Scholar] [CrossRef]
  43. Pfaffle, H.N.; Wang, M.; Gheorghiu, L.; Ferraiolo, N.; Greninger, P.; Borgmann, K.; Settleman, J.; Benes, C.H.; Sequist, L.V.; Zou, L.; et al. EGFR-activating mutations correlate with a Fanconi anemia-like cellular phenotype that includes PARP inhibitor sensitivity. Cancer Res. 2013, 73, 6254–6263. [Google Scholar] [CrossRef]
  44. Ye, H.; Zhang, X.; Chen, Y.; Liu, Q.; Wei, J. Ranking novel cancer driving synthetic lethal gene pairs using TCGA data. Oncotarget 2016, 7, 55352–55367. [Google Scholar] [CrossRef] [Green Version]
  45. Garcia-Higuera, I.; Taniguchi, T.; Ganesan, S.; Meyn, M.S.; Timmers, C.; Hejna, J.; Grompe, M.; D’Andrea, A.D. Interaction of the Fanconi anemia proteins and BRCA1 in a common pathway. Mol. Cell 2001, 7, 249–262. [Google Scholar] [CrossRef]
  46. Li, L.; Gu, X.; Yue, J.; Zhao, Q.; Lv, D.; Chen, H.; Xu, L. Acquisition of EGFR TKI resistance and EMT phenotype is linked with activation of IGF1R/NF-kappaB pathway in EGFR-mutant NSCLC. Oncotarget 2017, 8, 92240–92253. [Google Scholar]
  47. Matusan-Ilijas, K.; Damante, G.; Fabbro, D.; Dordevic, G.; Hadzisejdic, I.; Grahovac, M.; Avirovic, M.; Grahovac, B.; Jonjic, N.; Lucin, K. EGFR expression is linked to osteopontin and Nf-kappaB signaling in clear cell renal cell carcinoma. Clin. Transl. Oncol. 2013, 15, 65–71. [Google Scholar] [CrossRef]
  48. Meylan, E.; Dooley, A.L.; Feldser, D.M.; Shen, L.; Turk, E.; Ouyang, C.; Jacks, T. Requirement for NF-kappaB signalling in a mouse model of lung adenocarcinoma. Nature 2009, 462, 104–107. [Google Scholar] [CrossRef]
  49. Baud, V.; Collares, D. Post-Translational Modifications of RelB NF-kappaB Subunit and Associated Functions. Cells 2016, 5, 22. [Google Scholar]
  50. Lin, C.; Lu, W.; Ren, Z.; Tang, Y.; Zhang, C.; Yang, R.; Chen, Y.; Cao, W.; Wang, L.; Wang, X.; et al. Elevated RET expression enhances EGFR activation and mediates EGFR inhibitor resistance in head and neck squamous cell carcinoma. Cancer Lett. 2016, 377, 1–10. [Google Scholar] [CrossRef]
  51. Wang, S.E.; Narasanna, A.; Perez-Torres, M.; Xiang, B.; Wu, F.Y.; Yang, S.; Carpenter, G.; Gazdar, A.F.; Muthuswamy, S.K.; Arteaga, C.L. HER2 kinase domain mutation results in constitutive phosphorylation and activation of HER2 and EGFR and resistance to EGFR tyrosine kinase inhibitors. Cancer Cell 2006, 10, 25–38. [Google Scholar] [CrossRef] [Green Version]
  52. Guo, G.; Gong, K.; Wohlfeld, B.; Hatanpaa, K.J.; Zhao, D.; Habib, A.A. Ligand-Independent EGFR Signaling. Cancer Res. 2015, 75, 3436–3441. [Google Scholar] [CrossRef] [Green Version]
  53. Han, W.; Carpenter, R.L.; Cao, X.; Lo, H.W. STAT1 gene expression is enhanced by nuclear EGFR and HER2 via cooperation with STAT3. Mol. Carcinog. 2013, 52, 959–969. [Google Scholar] [CrossRef]
  54. Byers, L.A.; Diao, L.; Wang, J.; Saintigny, P.; Girard, L.; Peyton, M.; Shen, L.; Fan, Y.; Giri, U.; Tumula, P.K.; et al. An epithelial-mesenchymal transition gene signature predicts resistance to EGFR and PI3K inhibitors and identifies Axl as a therapeutic target for overcoming EGFR inhibitor resistance. Clin. Cancer Res. 2013, 19, 279–290. [Google Scholar] [CrossRef]
  55. Elkabets, M.; Pazarentzos, E.; Juric, D.; Sheng, Q.; Pelossof, R.A.; Brook, S.; Benzaken, A.O.; Rodon, J.; Morse, N.; Yan, J.J.; et al. AXL mediates resistance to PI3Kalpha inhibition by activating the EGFR/PKC/mTOR axis in head and neck and esophageal squamous cell carcinomas. Cancer Cell 2015, 27, 533–546. [Google Scholar] [CrossRef]
  56. Ghiso, E.; Migliore, C.; Ciciriello, V.; Morando, E.; Petrelli, A.; Corso, S.; De Luca, E.; Gatti, G.; Volante, M.; Giordano, S. YAP-Dependent AXL Overexpression Mediates Resistance to EGFR Inhibitors in NSCLC. Neoplasia 2017, 19, 1012–1021. [Google Scholar] [CrossRef]
  57. Guo, G.; Gong, K.; Ali, S.; Ali, N.; Shallwani, S.; Hatanpaa, K.J.; Pan, E.; Mickey, B.; Burma, S.; Wang, D.H.; et al. A TNF-JNK-Axl-ERK signaling axis mediates primary resistance to EGFR inhibition in glioblastoma. Nat. Neurosci. 2017, 20, 1074–1084. [Google Scholar] [CrossRef]
  58. Tian, Y.; Zhang, Z.; Miao, L.; Yang, Z.; Yang, J.; Wang, Y.; Qian, D.; Cai, H.; Wang, Y. Anexelekto (AXL) Increases Resistance to EGFR-TKI and Activation of AKT and ERK1/2 in Non-Small Cell Lung Cancer Cells. Oncol. Res. 2016, 24, 295–303. [Google Scholar] [CrossRef]
  59. Qu, Y.; Wu, X.; Yin, Y.; Yang, Y.; Ma, D.; Li, H. Antitumor activity of selective MEK1/2 inhibitor AZD6244 in combination with PI3K/mTOR inhibitor BEZ235 in gefitinib-resistant NSCLC xenograft models. J. Exp. Clin. Cancer Res. 2014, 33, 52. [Google Scholar] [CrossRef]
  60. Buonato, J.M.; Lazzara, M.J. ERK1/2 blockade prevents epithelial-mesenchymal transition in lung cancer cells and promotes their sensitivity to EGFR inhibition. Cancer Res. 2014, 74, 309–319. [Google Scholar] [CrossRef]
  61. Kachroo, P.; Lee, M.H.; Zhang, L.; Baratelli, F.; Lee, G.; Srivastava, M.K.; Wang, G.; Walser, T.C.; Krysan, K.; Sharma, S.; et al. IL-27 inhibits epithelial-mesenchymal transition and angiogenic factor production in a STAT1-dominant pathway in human non-small cell lung cancer. J. Exp. Clin. Cancer Res. 2013, 32, 97. [Google Scholar] [CrossRef]
  62. Kaowinn, S.; Kaewpiboon, C.; Koh, S.S.; Kramer, O.H.; Chung, Y.H. STAT1HDAC4 signaling induces epithelialmesenchymal transition and sphere formation of cancer cells overexpressing the oncogene, CUG2. Oncol. Rep. 2018, 40, 2619–2627. [Google Scholar]
  63. Li, B.; Huang, C. Regulation of EMT by STAT3 in gastrointestinal cancer (Review). Int. J. Oncol. 2017, 50, 753–767. [Google Scholar] [CrossRef]
  64. Schmitt, N.C.; Trivedi, S.; Ferris, R.L. STAT1 Activation Is Enhanced by Cisplatin and Variably Affected by EGFR Inhibition in HNSCC Cells. Mol. Cancer Ther. 2015, 14, 2103–2111. [Google Scholar] [CrossRef]
  65. Cheng, C.C.; Lin, H.C.; Tsai, K.J.; Chiang, Y.W.; Lim, K.H.; Chen, C.G.; Su, Y.W.; Peng, C.L.; Ho, A.S.; Huang, L.; et al. Epidermal growth factor induces STAT1 expression to exacerbate the IFNr-mediated PD-L1 axis in epidermal growth factor receptor-positive cancers. Mol. Carcinog. 2018, 57, 1588–1598. [Google Scholar] [CrossRef]
  66. Phuchareon, J.; McCormick, F.; Eisele, D.W.; Tetsu, O. EGFR inhibition evokes innate drug resistance in lung cancer cells by preventing Akt activity and thus inactivating Ets-1 function. Proc. Natl. Acad. Sci. USA 2015, 112, E3855–E3863. [Google Scholar] [CrossRef] [Green Version]
  67. Tetsu, O.; Eisele, D.W.; McCormick, F. Resistance to EGFR-targeted therapy by Ets-1 inactivation. Cell Cycle 2015, 14, 3211–3212. [Google Scholar] [CrossRef] [Green Version]
Cells 08 00740 g001 550
Figure 1. Quantification of EGFR and p-EGFR expression in a panel of NSCLC PDX models. (a) EGFR expression was determined by IHC on tissue micro arrays (TMAs) of the NSCLC PDX panel: sections were incubated with anti-human EGFR antibody (1:36) overnight at 4 °C, followed by DAB staining and hematoxylin counterstaining. Digitalized images of the IHC slides were evaluated to determine the percentage of EGFR positive areas using an inhouse software. The computerized analysis was used to quantify the EGFR expression using color classification and morphological image processing techniques. (b) The expression of EGFR and p-EGFR was determined by RPPA in a subset of the NSCLC PDX panel (n = 27). The samples were printed onto nitrocellulose covered microscope slides in five serial solutions and two replicates per dilution. Arrays were labeled with specific antibodies listed in Table 4. Median normalized data were used to compare expression levels between groups of samples. The ratio between phosphorylated and total protein was calculated by calculating the difference between the log-transformed phospho-protein expression and the log-transformed total protein expression. LXFA 677 showed high levels of EGFR and mean levels of p-EGFR. (c) LXFA 677 was implanted subcutaneously in NMRI nude mice and treatment started when a median tumor volume of 250 mm3 was achieved. Animals were assigned to the respective treatment arms. Dosing and schedule of the compounds are shown in Table 2. Treatment duration is indicated as lines on top of the figure. LXFA 677 showed EGFR dependent growth due to its sensitivity towards anti-EGFR monoclonal antibody, Cetuximab, as well as selective tyrosine kinase inhibitor (TKI), Erlotinib. Of note, the multi-TKI Sorafenib exhibited only marginal activity. (d) LXFA 677 was implanted subcutaneously in NMRI nude mice and treatment started when a median tumor volume of 250 mm3 was achieved. Animals were assigned to the respective treatment arms. Dosing and schedule of the compounds are shown in Table 2. Treatment duration is indicated as lines on top of the figure. LXFA 677 displayed a dose-dependent sensitivity towards treatment with Gefitinib.
Figure 1. Quantification of EGFR and p-EGFR expression in a panel of NSCLC PDX models. (a) EGFR expression was determined by IHC on tissue micro arrays (TMAs) of the NSCLC PDX panel: sections were incubated with anti-human EGFR antibody (1:36) overnight at 4 °C, followed by DAB staining and hematoxylin counterstaining. Digitalized images of the IHC slides were evaluated to determine the percentage of EGFR positive areas using an inhouse software. The computerized analysis was used to quantify the EGFR expression using color classification and morphological image processing techniques. (b) The expression of EGFR and p-EGFR was determined by RPPA in a subset of the NSCLC PDX panel (n = 27). The samples were printed onto nitrocellulose covered microscope slides in five serial solutions and two replicates per dilution. Arrays were labeled with specific antibodies listed in Table 4. Median normalized data were used to compare expression levels between groups of samples. The ratio between phosphorylated and total protein was calculated by calculating the difference between the log-transformed phospho-protein expression and the log-transformed total protein expression. LXFA 677 showed high levels of EGFR and mean levels of p-EGFR. (c) LXFA 677 was implanted subcutaneously in NMRI nude mice and treatment started when a median tumor volume of 250 mm3 was achieved. Animals were assigned to the respective treatment arms. Dosing and schedule of the compounds are shown in Table 2. Treatment duration is indicated as lines on top of the figure. LXFA 677 showed EGFR dependent growth due to its sensitivity towards anti-EGFR monoclonal antibody, Cetuximab, as well as selective tyrosine kinase inhibitor (TKI), Erlotinib. Of note, the multi-TKI Sorafenib exhibited only marginal activity. (d) LXFA 677 was implanted subcutaneously in NMRI nude mice and treatment started when a median tumor volume of 250 mm3 was achieved. Animals were assigned to the respective treatment arms. Dosing and schedule of the compounds are shown in Table 2. Treatment duration is indicated as lines on top of the figure. LXFA 677 displayed a dose-dependent sensitivity towards treatment with Gefitinib.
Cells 08 00740 g001
Cells 08 00740 g002 550
Figure 2. Induction of resistance in NSCLC PDX LXFA 677 in vivo. (a): Tumor growth curves for individual tumors under different doses of Gefitinib. The orange line is indicating the duration of the treatment. The respective dose per day is shown above the line. (b): Tumor growth curves for individual tumors after re-implantation of tumors emerging under constant Gefitinib treatment. The orange line is indicating the duration of the treatment. The respective dose per day is shown above the line. (c): Determination of copy number variation (left diagram) and mRNA expression level (right diagram) was determined for three resistant sublines and the treatment naïve NSCLC PDX LXFA 677. The relative expression was determined by using 2.0 (copy number) and 1.0 (mRNA expression level) as default for the treatment naïve line.
Figure 2. Induction of resistance in NSCLC PDX LXFA 677 in vivo. (a): Tumor growth curves for individual tumors under different doses of Gefitinib. The orange line is indicating the duration of the treatment. The respective dose per day is shown above the line. (b): Tumor growth curves for individual tumors after re-implantation of tumors emerging under constant Gefitinib treatment. The orange line is indicating the duration of the treatment. The respective dose per day is shown above the line. (c): Determination of copy number variation (left diagram) and mRNA expression level (right diagram) was determined for three resistant sublines and the treatment naïve NSCLC PDX LXFA 677. The relative expression was determined by using 2.0 (copy number) and 1.0 (mRNA expression level) as default for the treatment naïve line.
Cells 08 00740 g002
Cells 08 00740 g003 550
Figure 3. Results of whole exome sequencing analysis of four different LXFA 677 sublines. (a): Mutational landscape of the four LXFA 677 sublines with more than 85% (n = 504) of genes with mutations being shared across clones. (b): Hierarchical clustering of genes with differing (=at least two; n = 53) mutation status. The whole exome sequencing was performed using the variant effect predictor (VEP). Candidate mutations were annotated and filtered considering only variants with moderate or high protein impact and those being rare in healthy populations (<1% in gnomAD). Mutations shared by all three resistant clones were annotated with protein functions using SIFT and Polyphen predictions from SNPnexus as depicted in detail in Table S2.
Figure 3. Results of whole exome sequencing analysis of four different LXFA 677 sublines. (a): Mutational landscape of the four LXFA 677 sublines with more than 85% (n = 504) of genes with mutations being shared across clones. (b): Hierarchical clustering of genes with differing (=at least two; n = 53) mutation status. The whole exome sequencing was performed using the variant effect predictor (VEP). Candidate mutations were annotated and filtered considering only variants with moderate or high protein impact and those being rare in healthy populations (<1% in gnomAD). Mutations shared by all three resistant clones were annotated with protein functions using SIFT and Polyphen predictions from SNPnexus as depicted in detail in Table S2.
Cells 08 00740 g003
Cells 08 00740 g004 550
Figure 4. The proteomic landscape was determined by RPPA analysis of four different sublines of LXFA 677. The median corrected normalized values were used for all subsequent data analyses. The logarithmic fold change (logFC) was calculated comparing the sum of all resistant lines with the treatment naïve tumor. (a): Hierarchical clustering of average protein expression across five dilutions. (b): Hierarchical clustering of average phosphor- to total protein ratios across five dilutions.
Figure 4. The proteomic landscape was determined by RPPA analysis of four different sublines of LXFA 677. The median corrected normalized values were used for all subsequent data analyses. The logarithmic fold change (logFC) was calculated comparing the sum of all resistant lines with the treatment naïve tumor. (a): Hierarchical clustering of average protein expression across five dilutions. (b): Hierarchical clustering of average phosphor- to total protein ratios across five dilutions.
Cells 08 00740 g004
Table
Table 1. Metadata of the donor patient.
Table 1. Metadata of the donor patient.
DiagnosisBiopsyPre-TreatmentStagePts AgePts GenderEthnicityHLA-Type
undifferentiated adenocarcinoma of the lungprimary tumorCisplatinT1N1Mx62 yearsmaleCaucasianA*31:01
B*56:01
DQB1*05:01
DRB1*10:01
DQB1*05:01
DRB1*01:01
A*11:01
B*37:01
C*06:02 no
C*01:02 no
Table
Table 2. Doses and schedules of compounds applied to tumor-bearing animals.
Table 2. Doses and schedules of compounds applied to tumor-bearing animals.
CompoundSupplierDaily Dose [mg/kg/d]RouteSchedule
Sorafenib# S8502, LC Laboratories, USA200per os0–11
Lapatinib#L-4899, LC Laboratories, USA150per os0–21
Osimertinib#O-7200, LC Laboratories, USA5per os0–21
Afatinib#A-8644, LC Laboratories, USA12.5per os0–21
Necitumumab#LY3012211, local pharmacy20intraperitonealtwice weekly for 3 weeks
Erlotinib# E-4007-25G, LC Laboratories, USA100per os0–21
CetuximabPZN 00493528, Merck Serono, USA20intravenous0, 7, 14
Sunitinib#S-8803, LC Laboratories, USA30per os0–21
Gefitinib#G-4408-1G, LC Laboratories, USAup to 50per os0—end of experiment
Table
Table 3. Primers designed for gene copy number and mRNA expression of selected genes.
Table 3. Primers designed for gene copy number and mRNA expression of selected genes.
TargetForward PrimerReverse PrimerUsed for
GAPDHCAA′ATT′CCA′TGG′CAC′CGT′CGCC′ACA′CCA′TCC′TAG′TTG′CGene copy number (housekeeping gene)
ERBB2CTGAACTGGTGTATGCAGATTGCTTCCGAGCGGCCAAGTCGene copy number
cMET GACATTTCCAGTCCTGCAGTCACTCCGATCGCACACATTTGTGene copy number/mRNA expression
ERBB3AATAAAAGGGCTATGAGGCGATACTAGCTTCCTTAGCTCTGTCTCTTTGAmRNA expression
PTENCAGTAAGCGTTTTTTTTCTTTGAAGATGTGTAAGGTCAAAAGGGTGGAAmRNA expression
HGFGAATACTGCAGACCAATGTGCTTTGCAAGTGAATGGAAGTCCTmRNA expression
18S CTACCACATCCAAGGAAGGCATTTTTCGTCACTACCTCCCCGmRNA expression (housekeeping gene)
Table
Table 4. Antibodies used for RPPA assay.
Table 4. Antibodies used for RPPA assay.
IndexDesignationSourceMono/Polyclonal
1Phospho-Akt (Ser473) (193H12)RabbitMonoclonal
2Phospho-PTEN (ser380/Thr382/383)RabbitPolyclonal
3p44/42 MAPKRabbitPolyclonal
4Stat3RabbitPolyclonal
5Phospho-p38 MAPK (Thr180/Tyr182)RabbitMonoclonal
6Phospho-Stat3 (Ser727)RabbitPolyclonal
7Stat1RabbitPolyclonal
8Phospho-MEK1/2 (Ser217/221)RabbitMonoclonal
9Phospho-p44/42 MAPK (Thr202/Tyr204)RabbitMonoclonal
10Phospho-c-jun (Ser73) (D47G9)RabbitMonoclonal
11NF-kB p65RabbitMonoclonal
12Phospho-NF-kB p65 (Ser536)RabbitMonoclonal
13Phospho-Met (Tyr1234/1235)RabbitMonoclonal
14p38 MAPKRabbitMonoclonal
15phospho-mTOR (Ser2448)RabbitMonoclonal
16Phospho-Her3/Erbb3 (tyr1289)RabbitMonoclonal
17HER3/ErbB3 (c-17)RabbitPolyclonal
18mTORRabbitPolyclonal
19c-Jun (N-term)RabbitMonoclonal
20Phospho-eIF4B (Ser422)RabbitPolyclonal
21eIF4BRabbitMonoclonal
22MEK1/2RabbitPolyclonal
23Akt (pan) (C67E7)RabbitMonoclonal
24PTEN (D4.3) XPRabbitMonoclonal
25Phospho-Stat1 (Y701)RabbitMonoclonal
26EGFR (D38B1)RabbitMonoclonal
27RelBRabbitMonoclonal
28Met (D1C2) XPRabbitMonoclonal
29Axl (C2B12)RabbitMonoclonal
30IkB alpha (44D4)RabbitMonoclonal
31HER2/ErbB2MouseMonoclonal
32IGF-I receptor BRabbitMonoclonal
33Phospho-EGFR (Tyr1068) (D7A5)RabbitMonoclonal
Table
Table 5. Doubling times of different sublines of NSCLC PDX LXFA 677.
Table 5. Doubling times of different sublines of NSCLC PDX LXFA 677.
LXFA 677 UntreatedLXFA 677 TreatedLXFA 677res1 UntreatedLXFA 677res1 TreatedLXFA 677res2 UntreatedLXFA 677res2 TreatedLXFA 677res3 UntreatedLXFA 677res3 Treated
Mean [d]11.8522.987.545.3715.0314.4313.6814.38
Std. Deviation [d]0.077.620.860.232.363.780.002.31
Nr. of values [n]22222212

Share and Cite

MDPI and ACS Style

Schueler, J.; Tschuch, C.; Klingner, K.; Bug, D.; Peille, A.-L.; de Koning, L.; Oswald, E.; Klett, H.; Sommergruber, W. Induction of Acquired Resistance towards EGFR Inhibitor Gefitinib in a Patient-Derived Xenograft Model of Non-Small Cell Lung Cancer and Subsequent Molecular Characterization. Cells 2019, 8, 740. https://doi.org/10.3390/cells8070740

AMA Style

Schueler J, Tschuch C, Klingner K, Bug D, Peille A-L, de Koning L, Oswald E, Klett H, Sommergruber W. Induction of Acquired Resistance towards EGFR Inhibitor Gefitinib in a Patient-Derived Xenograft Model of Non-Small Cell Lung Cancer and Subsequent Molecular Characterization. Cells. 2019; 8(7):740. https://doi.org/10.3390/cells8070740

Chicago/Turabian Style

Schueler, Julia, Cordula Tschuch, Kerstin Klingner, Daniel Bug, Anne-Lise Peille, Leanne de Koning, Eva Oswald, Hagen Klett, and Wolfgang Sommergruber. 2019. "Induction of Acquired Resistance towards EGFR Inhibitor Gefitinib in a Patient-Derived Xenograft Model of Non-Small Cell Lung Cancer and Subsequent Molecular Characterization" Cells 8, no. 7: 740. https://doi.org/10.3390/cells8070740

APA Style

Schueler, J., Tschuch, C., Klingner, K., Bug, D., Peille, A. -L., de Koning, L., Oswald, E., Klett, H., & Sommergruber, W. (2019). Induction of Acquired Resistance towards EGFR Inhibitor Gefitinib in a Patient-Derived Xenograft Model of Non-Small Cell Lung Cancer and Subsequent Molecular Characterization. Cells, 8(7), 740. https://doi.org/10.3390/cells8070740

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop