Suppression of RAF/MEK or PI3K synergizes cytotoxicity of receptor tyrosine kinase inhibitors in glioma tumor-initiating cells
- Takashi Shingu1, 2Email author,
- Lindsay Holmes1, 8,
- Verlene Henry1,
- Qianghu Wang3, 4,
- Khatri Latha1,
- Anupama E. Gururaj1,
- Laura A. Gibson1, 2,
- Tiffany Doucette1,
- Frederick F. Lang1,
- Ganesh Rao1,
- Liang Yuan2,
- Erik P. Sulman4,
- Nicholas P. Farrell5,
- Waldemar Priebe6,
- Kenneth R. Hess7,
- Yaoqi A. Wang2,
- Jian Hu2 and
- Oliver Bögler1, 9Email author
© Shingu et al. 2016
Received: 11 October 2015
Accepted: 26 January 2016
Published: 9 February 2016
The majority of glioblastomas have aberrant receptor tyrosine kinase (RTK)/RAS/phosphoinositide 3 kinase (PI3K) signaling pathways and malignant glioma cells are thought to be addicted to these signaling pathways for their survival and proliferation. However, recent studies suggest that monotherapies or inappropriate combination therapies using the molecular targeted drugs have limited efficacy possibly because of tumor heterogeneities, signaling redundancy and crosstalk in intracellular signaling network, indicating necessity of rationale and methods for efficient personalized combination treatments. Here, we evaluated the growth of colonies obtained from glioma tumor-initiating cells (GICs) derived from glioma sphere culture (GSC) in agarose and examined the effects of combination treatments on GICs using targeted drugs that affect the signaling pathways to which most glioma cells are addicted.
Human GICs were cultured in agarose and treated with inhibitors of RTKs, non-receptor kinases or transcription factors. The colony number and volume were analyzed using a colony counter, and Chou-Talalay combination indices were evaluated. Autophagy and apoptosis were also analyzed. Phosphorylation of proteins was evaluated by reverse phase protein array and immunoblotting.
Increases of colony number and volume in agarose correlated with the Gompertz function. GICs showed diverse drug sensitivity, but inhibitions of RTK and RAF/MEK or PI3K by combinations such as EGFR inhibitor and MEK inhibitor, sorafenib and U0126, erlotinib and BKM120, and EGFR inhibitor and sorafenib showed synergy in different subtypes of GICs. Combination of erlotinib and sorafenib, synergistic in GSC11, induced apoptosis and autophagic cell death associated with suppressed Akt and ERK signaling pathways and decreased nuclear PKM2 and β-catenin in vitro, and tended to improve survival of nude mice bearing GSC11 brain tumor. Reverse phase protein array analysis of the synergistic treatment indicated involvement of not only MEK and PI3K signaling pathways but also others associated with glucose metabolism, fatty acid metabolism, gene transcription, histone methylation, iron transport, stress response, cell cycle, and apoptosis.
Inhibiting RTK and RAF/MEK or PI3K could induce synergistic cytotoxicity but personalization is necessary. Examining colonies in agarose initiated by GICs from each patient may be useful for drug sensitivity testing in personalized cancer therapy.
KeywordsGlioma Stem cells 3-D culture Combination therapy Synergistic Autophagy Receptor tyrosine kinase Gompertz
Malignant gliomas, including anaplastic astrocytoma and glioblastoma, are common primary tumors of the central nervous system, and characterized by aggressive cellular proliferation, diffuse infiltration, and resistance to cell death . Diffuse infiltration of malignant glioma cells into adjacent brain parenchyma makes gross total removal by surgery difficult. Glioma cells that are resistant to cell death-inducing stimuli cause recurrence after radiation therapy and chemotherapy. For these reasons, the outcome of patients with malignant gliomas is poor . The tumorigenesis of malignant gliomas has been associated with a number of alterations in related genes causing functional loss of tumor suppressors including PTEN, p53, NF1, RB1 or p16 proteins, and hyperfunction of oncogenic proteins such as EGFR, cyclin D1/3, E3 ubiquitin-protein ligase Mdm2, Mdm4, Met, Bcl-2, cyclin dependent kinase 4/6, PDGFRA and PI3K . Among them, alterations of EGFR including gene amplification and activation mutations are observed in 40–57 % of malignant gliomas [2, 3]. The Cancer Genome Atlas (TCGA) Study also indicates that the vast majority of glioblastomas have aberrant RTK/RAS/PI3K signaling pathways . Therefore, malignant glioma cells are thought to be addictive to these aberrant signaling pathways for their survival and proliferation, and thus targeting these signaling pathways has been thought to be effective on therapy of malignant gliomas . However, a large number of clinical trials in which the targeted therapies were examined for malignant glioma patients indicate that monotherapies using single anticancer agents would have limited efficacy . Recent studies indicate that inter- and intratumoral heterogeneities are major obstacles to the efficient targeted therapies for solid tumors [2–5]. In addition, key signaling nodes to which multiple oncogenic signaling pathways converge might not be inhibited sufficiently by monotherapies or even by combination therapies because of unexpected redundancy or alterations of feedback in the signaling pathway network with complex crosstalk [6–8]. Therefore, it is necessary to find rationale of effective combination treatments and to optimize therapies to individual patients.
Glioma tumor-initiating cells (GICs) are thought to play a critical role in initiation, regrowth and recurrence of the tumor [1, 9]. Therefore, GICs need to be considered a target of glioma therapy, and it might be essential to evaluate effects of therapies on GICs from individual patients in order to predict the clinical relevance of the treatments and optimize therapy to the patients. In vitro screening of anticancer therapy has been done mainly by clonogenic assay because the effect of the therapy on clonogenicity of the tumor cells is thought to be associated with the clinical therapeutic efficacy . However, clonogenic assay using GICs has been a challenge because GICs aggregate in the stem cell culture media, and evaluation of the accurate tumor neurosphere/colony number requires single cell culture system or semi-solid matrix to prevent cell/colony aggregation. Single cell culture systems need large numbers of wells/plates and are not well suited for high-throughput screening of combination therapies . Although colony formation assays of GICs or neural stem cells using gels have been reported, the growth of the colonies initiated by these cells in soft agar has not yet been well characterized [12–15]. In addition, a recent study suggested that proliferating cells with limited self-renewal capacity are more tumorigenic than glioma stem-like cells and thus therapeutic effects on these proliferating cells might be a better predictor for the in vivo efficacy . Therefore, in drug sensitivity testing of gliomas, method by which we can evaluate both clonogenicity of GICs and cell proliferation of GICs and their descendant cells may be useful.
In this study, we cultured GICs in agarose and evaluated the number and volume of the colonies that reflect clonogenicity and cell proliferation, respectively, using a colony counter GelCount. With this method, we examined efficiency of combination treatments using RTK inhibitors, non-receptor kinase inhibitors and transcription factor inhibitors that affect the signaling pathways to which most glioma cells are thought to be addicted.
Antibodies and reagents
Erlotinib, lapatinib and sorafenib were purchased from LC laboratories (Woburn, MA), BKM120 was from Novartis (Basel, Switzerland), PD98059 and PP2 were from Selleck Chemicals (Houston, TX), U0126 and 3-methyladenine (3-MA) were from Sigma-Aldrich (St. Louis, MO), c-Myc inhibitor II was from EMD Millipore Corporation (Billerica, MA). Imatinib mesylate was generously provided from Novartis. A polynuclear platinum BBR3610 was synthesized by Dr. Nicholas P Farrelle (Virginia Commonwealth University) . WP1066, an inhibitor of tyrosine phosphorylated STAT3 and STAT5 was synthesized by Dr. Waldemar Priebe (The University of Texas MD Anderson Cancer Center) . These reagents except for 3-MA, BBR3610 and imatinib were dissolved in DMSO (Sigma-Aldrich). 3-MA was dissolved in culture media, and imatinib and BBR3610 were dissolved in PBS. Antibodies for Akt, AMPK, Atg5, Bad, c-Myc, EGFR, ERK, Met, poly-ADP ribose polymerase (PARP), pyruvate kinase isozyme M2 (PKM2), and ribosomal protein S6, or phosphorylated forms of Akt (Ser473), AMPK (Thr172), Bad (Ser136), EGFR (Tyr1173), ERK (Thr202/Tyr204), Met (Tyr1234/1235), and S6 (Ser235/236) were obtained from Cell Signaling Technology, Inc. (Danvers, MA). Antibodies for Bcl-2, Bcl-XL, β-catenin, Mcl-1, p53, and PTEN were obtained from Santa Cruz Biotechnology, Inc. (Santa Cruz, CA). Anti-LC3B antibody was obtained from Novus Biologicals, Inc. (Littleton, CO). Antibody for CD133 was obtained from Abcam plc (Cambridge, UK). Antibodies for lamin B and nestin were obtained from EMD Millipore. Antibodies for β-actin and vinculin were from Sigma-Aldrich.
Human malignant glioma cell line U87-MG was from American Type Culture Collection (Manassas, VA), and human malignant glioma cell lines LNZ308 and LN428 were from Dr. Nicolas de Tribolet (Lausanne, Switzerland), during 2003–2007. These cells were maintained in DMEM (Life Technologies Corporation, Carlsbad, CA) containing 10 % FBS (Life Technologies), 100 units/ml penicillin G, 100 μg/ml streptomycin, and 50 mg/ml l-glutamine (Life Technologies) in a humidified atmosphere containing 7 % CO2 at 37 °C. Human glioma tumor-initiating cells derived from glioma sphere culture were established by Dr. Frederick F Lang (The University of Texas MD Anderson Cancer Center) as described previously, in which GSC11 and GSC47 cells are non-mesenchymal subtype and GSC2 and GSC20 cells are mesenchymal subtype [19, 20]. These cells are maintained in DMEM/Ham’s F12 medium (F12) (Life Technologies) containing 1× B27 supplement (Life Technologies), 20 ng/ml EGF and 20 ng/ml basic FGF (Life Technologies), 100 units/ml penicillin G, 100 μg/ml streptomycin, and 50 mg/ml l-glutamine in a humidified atmosphere containing 7 % CO2 at 37 °C. Although the authentication of the purchased cell lines was not done by authors, DNA fingerprinting was done for testing cell line contamination using GenomeLab Human STR Primer Set according to the manufacturer’s instruction (Beckman Coulter, Brea, CA). Infection by mycoplasma was examined by Mycoalert (Lonza, Basel, Switzerland) and infected cells were treated with BM cyclin (Roche Diagnostics, Indianapolis, IN).
3D cell culture in agarose gel
For analyses of number and volume (total biomass) of tumor spheres in agarose gel, culture media containing 0.7 % agarose (Thermo Fisher Scientific, Waltham, MA) was used as the bottom layer, and culture media containing 0.4–0.5 % agarose and tumor cells were placed in the middle layer, and then, culture media were added to the well in 96 or 24 well plates (Corning Incorporated, NY) . Cells were seeded at 100, 150 and 200 cells in 96 well plate as triplicate, or 250 and 500 cells in 24 well plate as duplicate. The final volumes of the bottom, middle and top layers for 96 well plate and 24 well plate were 100, 100 and 50, and 500, 500, and 200 μl, respectively. After 4 h of seeding cells, culture media containing indicated reagents were added and cells were further incubated. Plates were scanned at the indicated time points in time course experiments or 12–14 days after addition of drugs in colony formation analysis, and number and total biomass of colonies larger than 40 μm in diameter were evaluated using GelCount colony counter system (Oxford Optronix Inc., UK).
Analysis of tumor growth
Surviving fraction in colony formation analysis and total biomass in volume analysis were calculated by the methods reported by Franken and Kajiwara, respectively [10, 21]. In vivo tumor growth in xenograft models can be described using an equation of Gompertz model as: V = V0 exp[(A/α)(1–exp[–αt])] where V0 = initial volume, A and α = constant, V = volume at any time t . In this study, colony number and total biomass were plotted against time and examined by non-linear regression analysis against Gompertz growth curve where V0 = initial colony number and initial total biomass, respectively. Coefficients of determination (R2) were calculated using GraphPad Prism software (GraphPad Software Inc., San Diego, CA). Combination indices of drugs in colony formation analysis were calculated by modified method of Chou and Talalay , and combination effect was thought to be synergistic when the combination index was less than 0.9. Tumor growth delay was calculated using the formula (Ti−C)/C where Ti and C were time (days) for treated and control tumors to reach double of total biomass at treatment, respectively .
After treatment with reagents, cells were harvested, washed with PBS, lysed in 50 mM Tris–HCl (pH 7.5), 150 mM NaCl, 1 % Igepal-630, 1 mM EDTA, aprotinin, leupeptin, protease inhibitor and phosphatase inhibitor (Roche Diagnostics), and 1 mM PMSF and centrifuged at 10,000×g at 4 °C for 20 min. The supernatant was used as the cell lysate. Extraction of nuclear and cytoplasmic fractions of cells was done using Nuclear Extraction Kit (Sigma-Aldrich) according to the manufacture’s instruction. Ten-fifteen microgram of proteins were subjected to standard western blotting and the blots were quantified with ImageJ software (NIH, Bethesda, MD).
Survival analysis of nude mice with intracranial tumor
The intracranial tumor model was made using male nude mice (Experimental Radiation Oncology at M. D. Anderson Cancer Center) as described previously . Through the screw guide, at a depth of 3.5 mm from the skull, 5 μl aliquots of 5 × 105 GICs in DMEM/F12 were inoculated at a rate of 1.0 μl/min. Mice were randomized into four groups (4–8 mice per each group) and treatment was initiated 1 week after tumor implantation. Erlotinib and sorafenib were dissolved in solvent, OraPlus (Paddock Laboratories Inc., Minneapolis, MN), and administered at doses of 50 mg/kg of erlotinib and 50 mg/kg of sorafenib orally three times a week for a total of twelve treatments. Mice of control group were treated with the solvent. All mouse studies were performed in the veterinary facilities of The University of Texas M. D. Anderson Cancer Center in accordance with IACUC, state, federal, and ethical regulations for experimental animal care.
Cell viability assay
For trypan blue dye exclusion assay, tumor cells were seeded at 2 × 104 cells per well in 24-well flat-bottomed plates and treated with reagents for indicated duration. Cells were dissociated by using accutase (Life Technologies) and the number of viable or dead cells (total >300 cells) was counted in at least three different fields.
Quantification of apoptotic or non-apoptotic cell death with annexin V and propidium iodide
After treatment with indicated reagents, cells were stained with FITC-conjugated annexin V and propidium iodide for 15 min according to the manufacturer’s instruction (BD Biosciences, San Jose, CA), and then analyzed with FACSCalibur (Becton, Dickinson and Company, Franklin Lakes, NJ) using Cell Quest Pro software (Becton, Dickinson and Company).
Analysis of LC3 localization
A GFP-tagged LC3 (GFP-LC3) expression vector was kindly provided by Drs. N Mizushima (The University of Tokyo, Tokyo, Japan) and T Yoshimori (Osaka University, Suita, Japan) . Tumor cells were seeded in 24 well plates and transfected with GFP-LC3 expression vector using FuGENE HD Transfection Reagent (Roche Diagnostics). Cells were incubated overnight and treated with indicated reagents for 8 h. Cells were dissociated by using accutase, fixed with 4 % paraformaldehyde, mounted onto glass slides with VECTASHIELD Mounting Medium (Vector Laboratories Inc., Burlingame, CA), and then examined using fluorescence microscope. The number of cells with or without GFP-LC3 punctate pattern (total >300 cells) was counted in at least three different fields.
siRNA was synthesized by Thermo Fisher Scientific (Dharmacon). The sequence of the siRNA targeting Atg5 was CAACUUGUUUCACGCUAUAdTdT. siCONTROL Non-Targeting siRNA#2 was purchased from Thermo Fisher Scientific and used as control siRNA. Cells were transfected with siRNAs (3 μg) by electroporation using Nucleofector and a Cell Line Nucleofector® Kit T (Amaxa Biosystems, Gaithersburg, MD) according to the manufacturer’s instruction. After incubated overnight, cells were harvested, dissociated by using accutase, and seeded in 24 well plates for cell viability assay or 6 well plates for immunoblotting. After 3–4 h at 37 °C, cells were treated with reagents and subjected to the indicated assays.
Reverse phase protein array
Reverse phase protein array (RPPA) study was performed as described previously . Briefly, serially diluted lysates were arrayed on nitrocellulose-coated slides (Grace Bio-Labs, Inc., Bend, OR) by Aushon 2470 Arrayer (Aushon BioSystems, Billerica, MA). Total 5808 array spots were arranged on each slide including the spots corresponding to positive and negative controls prepared from mixed cell lysates and dilution buffer, respectively. Each slide was probed with a validated primary antibody plus a biotin-conjugated secondary antibody. A panel of 171 antibodies with a Pearson correlation coefficient between RPPA and western blotting of greater than 0.7 were used in RPPA study. The signal obtained was amplified using a Dako Cytomation–catalyzed system (Dako, Glostrup, Denmark) and visualized by DAB colorimetric reaction. The slides were scanned, analyzed, and quantified using a customerized-software Microvigene (VigeneTech Inc., Carlisle, MA) to generate spot intensity. Each dilution curve was fitted with a logistic model (“Supercurve Fitting” developed by the Department of Bioinformatics and Computational Biology in The University of Texas MD Anderson Cancer Center, “http://bioinformatics.mdanderson.org/OOMPA”). Results were normalized, transformed to base 2 logarithms and median centered, and the data from RPPA with three replicate samples were analyzed using two-way analysis of variance with an interaction term to assess whether the combined effect was significantly different from the sum of the individual effects.
Gene expression analysis
Paired-end whole transcriptome sequencing of GIC cultures was performed on the Illumina HiSeq platform after random priming and rRNA reduction. Each sample generated about 50 million paired-ends; each end was 75 bp in size. Short transcript reads were mapped to 20,010 human protein coding genes in Ensembl reference transcriptome (ENSEMBL version 64). Downstream data analyses and RPKM (reads per kilobase per million reads) values were generated using pipeline for RNA sequencing data analysis (PRADA, PMID: 24695405) by involving Burroughs-Wheeler alignment, Samtools, and Genome Analysis Toolkit.
Paired Student’s t test (two-tailed) was used for comparison between two groups, and ANOVA with Tukey’s post hoc test was used for four groups. For the animal experiment, pairs of Kaplan–Meier survival curves were compared by the log-rank Mantel-Cox test using GraphPad Prism software. Differences were considered statistically significant when provided P was less than 0.05.
Glioma cell growth in agarose
Growth suppression of glioma cells by monotherapies and combination therapies
The TCGA study has revealed that 90 % of glioblastomas have aberrant signaling in RTK/RAS/PI3K signaling pathways and thus molecules in these signaling pathways have been thought to be therapeutic targets . However, successful treatments targeting these signaling pathways remain to be clarified [5–8]. Considering multiplicity of cancerous signaling pathways with complex crosstalk in a tumor cell, targeting relatively upstream molecule such as RTK that could affect multiple oncogenic signaling pathways might be reasonable. To overcome signaling redundancy, combined inhibition of multiple molecules that upregulate the same signaling pathway might be necessary. Recent study suggests that inhibition of relatively distal downstream molecule such as mTOR causes activation of the upstream oncogenic molecule such as Akt by decreasing feedback signaling and thus reduces efficacy of erlotinib in the combination therapy . Therefore, we hypothesized that effective combination therapies consist of inhibitors of relatively upstream molecules including cellular membrane receptors and their proximal downstream molecules, and then tested combination effects of RTK inhibitors and RAF/MEK or PI3K inhibitors.
IC50 values of drugs in glioma tumor-initiating cells and U87 cultured in agarose
Receptor tyrosine kinase inhibitors
0.11 ± 0.03
0.49 ± 0.15
0.82 ± 0.34
3.08 ± 0.55
3.72 ± 0.60
0.64 ± 0.22
1.86 ± 0.16
4.39 ± 1.08
3.04 ± 0.99
7.01 ± 2.20
2.59 ± 0.50
1.09 ± 0.44
2.21 ± 0.16
3.51 ± 0.25
5.24 ± 0.01
6.37 ± 0.90
11.0 ± 1.57
6.28 ± 2.24
3.24 ± 2.10
Non-receptor kinase inhibitors
14.0 ± 4.72
18.0 ± 3.25
4.69 ± 1.33
11.0 ± 0.25
11.1 ± 5.35
0.41 ± 0.05
0.46 ± 0.21
0.36 ± 0.05
0.49 ± 0.10
0.59 ± 0.21
1.93 ± 0.50
6.39 ± 2.46
4.33 ± 0.34
NA > 20
Transcription factor inhibitors
2.46 ± 0.62
1.55 ± 0.13
1.96 ± 0.13
1.59 ± 0.40
Myc II (μM)
22.3 ± 3.57
13.3 ± 1.50
19.0 ± 3.47
27.0 ± 5.68
DNA cross-linking agent
46.6 ± 23.0
87.3 ± 33.5
117.3 ± 71
1.46 ± 1.26
Synergistic combinations of targeted drugs in glioma tumor-initiating cells
Erlotinib + Sorafenib
0.85 ± 0.23
1.30 ± 0.16
1.21 ± 0.26
1.04 ± 0.24
Lapatinib + Sorafenib
0.82 ± 0.27
1.17 ± 0.40
1.18 ± 0.05
1.32 ± 0.06
Erlotinib + U0126
0.98 ± 0.24
2.14 ± 0.75
0.73 ± 0.08
0.52 ± 0.31
Erlotinib + PD98059
0.79 ± 0.17
Lapatinib + U0126
0.86 ± 0.12
Erlotinib + BKM120
1.16 ± 0.19
1.04 ± 0.24
0.71 ± 0.05
Sorafenib + U0126
0.67 ± 0.13
0.73 ± 0.09
1.20 ± 0.59
0.56 ± 0.30
Combinations of EGFR inhibitor with sorafenib were synergistic in GSC11 cells, which are of the non-mesenchymal subtype, but not synergistic in GSC47 cells, which are also of the non-mesenchymal subtype (Table 2). In addition, GSC2 and GSC20, which are both of the mesenchymal subtype, showed different sensitivities to combination therapies (Table 2). These results indicate that GICs show diverse sensitivity to therapies and that personalization is necessary for therapy even in the same subtype of glioma.
Next we evaluated the combined effect of erlotinib and sorafenib on the total biomass in time course study. Both tumor growth delay analysis and combination index at 1 week after treatment indicated synergistic effects in GSC11 cells, while tumor growth delay analysis showed synergy but combination index indicated only additivity in GSC20 cells (Fig. 1c, d; Additional file 1: Table S3). Considering these findings in total biomass analysis together with those from the colony formation study (Table 2), we concluded that GSC11 cells are sensitive to combination of erlotinib with sorafenib, and evaluated the effects of the combination therapy on nude mice bearing brain tumor (Fig. 1e). Although analyses of the survival curves and the median survival time did not indicate statistically significant difference, the median survival time of the combination group (83d) was remarkably improved by the combination treatment compared with those of control (74.5d), erlotinib (74d) and sorafenib (70d) groups, even with the doses that did not improve the survival in the monotherapies.
Combination treatment-induced apoptotic and autophagic cell death
Combination treatment-induced alterations in signaling pathways
Among several mathematical representations of tumor growth, the Gompertz model has been used in studies of experimental therapies because it has been shown to correlate not only with in vitro cell number and tumor sphere volume but also with in vivo tumor weight and volume in previous studies using conventional cell lines [22, 27, 33–36]. Our study is, as far as we know, the first showing that both number and total volume of colonies of GICs cultured in agarose correlate with the Gompertz function. Correlation of the Gompertz curve with the in vitro tumor growth curves in our study might indicate that tumor growth pattern in this culture system resembles in vivo tumor growth, though the time scale would be different. One of the advantages of 3D culture with gel is that aggregation of cells or colonies can be prevented and thus we can evaluate both accurate colony number that is thought to reflect clonogenicity of GICs and volume of the colonies formed not by cell aggregation but by cell proliferation. In addition, approximately 20–60 colonies are formed in 1 well of 96-well plate that might enable us to analyze the clonogenicity with lower number of wells and plates but higher statistical power than single cell culture system. We found that combination therapy of erlotinib and sorafenib in GSC11 cells was synergistic in both analyses of colony number and total biomass in vitro using the 3D culture method, and that the in vivo combination therapy showed tendency of improvement of survival by the combination therapy, indicating correlation in the therapeutic efficacy between this in vitro assay system and in vivo experimental model. Although we used minimal doses of erlotinib and sorafenib that did not improve the median survival as monotherapies, the combination treatment tended to improve the survival, suggesting supra-additive effect by the combination therapy in vivo. Since the maximum tolerated doses of erlotinib and sorafenib in combination therapy in mouse are higher, efficiency of the combination therapy might be improved by increasing doses of the drugs .
Based on the TCGA study that showed possible addiction of glioblastoma cells to several signaling pathways, we evaluated the efficiency of single and combination treatments using RTK inhibitors, non-receptor kinase inhibitors, and transcription factor inhibitors . Since clinical trials of targeted therapies for malignant glioma patients have shown that monotherapies and inappropriate combination therapies have limited efficacy, finding rationale of efficient combination treatments is necessary [4, 8]. This is the first report showing synergistic effects of combinations of EGFR inhibitor with sorafenib, sorafenib with U0126, EGFR inhibitor with MEK inhibitor, and erlotinib with BKM120 on GICs beyond their subtype. Among these synergistic treatments, combination of erlotinib with sorafenib was reported to be synergistic in colorectal and lung cancer cells [37, 38]. Although U0126 was reported to be protective against sorafenib-induced cytotoxicity in hepatocellular carcinoma (HCC) cells, combination of sorafenib with another MEK inhibitor CI-1040 was synergistic associated with increased Bim in HCC cells [39, 40]. Although synergistic effect by erlotinib and U0126 has not been reported, a MEK inhibitor was shown to decrease resistance of non-small cell lung cancer (NSCLC) cells to EGFR inhibitors . Similarly, although synergistic effect by erlotinib and BKM120 has not been reported, combinations of erlotinib and PI3K/mTOR inhibitors showed synergistic cytotoxicity in cells of prostate cancer, NSCLC and ovarian cancer in previous studies [42–44]. Findings in these sporadic studies are consistent with those in our systematic evaluation that show significance of combined inhibition of RTK with PI3K or RAF/MEK. Importantly, inhibition of MEK and PI3K did not induce synergistic effect in our study. Instead, RPPA study indicates that synergistic combination therapy could affect multiple signaling pathways including not only RAS/ERK and PI3K/Akt but also other pathways that are critical for tumor cell survival and proliferation. Therefore, the rationale for the effective combination treatments of GICs might be suppression of both an RTK that activates multiple signaling pathways which form complex signaling network and RAF/MEK or PI3K that could play critical roles in the redundancy of the key signaling pathways.
However, we found significant variations in drug sensitivity and intracellular signaling pathways between GICs. Even GSC2 and GSC20 cells, which are classified into the same subtype of gliomas, showed different sensitivity to most of the mono- and combination therapies. MEK/ERK/NFκB signaling induces translocation of PKM2 into nucleus where PKM2 binds to β-catenin that activates transcription of genes related with cell proliferation including c-Myc and CDK1. [45–47]. Nuclear PKM2 also phosphorylates spindle checkpoint protein Bub3 resulting in mitosis progression . Nuclear β-catenin regulates transcription of Bcl-2 and Mcl-1. . Therefore, alterations of these molecules in GSC11 cells induced by combination of erlotinib with sorafenib were consistent with the increased cytotoxicity by the combination therapy. However, GSC20 cells responded to combination of erlotinib with sorafenib in a different manner from that of GSC11 cells regarding both cell death (Fig. 3b) and signaling molecules (data not shown). In addition, phosphorylated Akt and ERK in GSC2 cells treated with erlotinib and sorafenib also showed different alteration patterns from those in GSC11 or GSC20 (data not shown). Given that intratumoral heterogeneity contributes to therapy resistance and recurrence of the tumor in a patient, our findings indicate that intertumoral heterogeneity of GICs in sensitivity to cytotoxic stimuli associated with diverse alterations in intracellular signaling pathways could be another reason for the failure of the “one-size-fit-all” designed treatment for malignant glioma patients [5, 8, 50]. In several types of cancers, novel classifications based on biomarkers or profiling are providing actionable and critical information to make clinical decisions . In glioma patients, there is a possibility that a small number of molecular features such as CpG island methylation phenotype (G-CIMP) or IDH1 mutation status could be useful for prediction of the outcome [32, 52]. In addition to these factors, different status of several key molecules could alter dependency of GICs on the key signaling pathways or crosstalk/feedback network [3, 14, 28, 29, 41, 53]. In this study, we found that at least GSC2, GSC11 and GSC20 are neither G-CIMP nor IDH1 mutant subtype (data not shown) and have different gene expression patterns regarding RTKs, PI3K, RAS and Rb pathways that could induce the diversity in the drug sensitivity of the GICs even in the same subtype (Additional file 3). However, our findings as well as previous studies indicate difficulty in predicting sensitivity of glioma cells to specific treatment by alterations in small number of genes or proteins [30, 54]. These findings have been fueling profiling approach using large data set [19, 32, 55]. However, the current subtype classification of gliomas might not be sufficient for selecting the best treatment for patients while certain extent of the effectiveness was found . Therefore, it is necessary to establish a more effective personalized optimization method for glioma therapy. Although hypothesis-driven studies with omic analyses that attempt to predict sensitivity of tumors are intensely performed, unpredictability of cancer cell response to specific treatment has not yet been overcome. Empiric approaches such as treatment sensitivity testing using cultured tumor cells obtained from each patient have also been studied but have not yet significantly improved outcome of glioma patients [56, 57]. In this study, we showed that synergistic combination treatments for different cell lines with diverse drug sensitivity can be efficiently identified by examining theory-driven combinations, which are RTK inhibitors and RAF/MEK or PI3K inhibitors, indicating that integration of rationale-based candidate selection with empiric screening approach is useful to optimize therapy for individual patients against intertumoral heterogeneity.
Role of autophagy in efficient combination therapy
Among the synergistic combination treatments, we focused on the combination of erlotinib and sorafenib because clinical trials for patients with tumors including malignant gliomas to examine efficacy of the combination therapy were ongoing when we started this study [50, 58]. This is the first report showing autophagic cell death induced by the combination treatment. Autophagy is thought to be a cellular process for homeostasis and adaptation to stress and is thus cytoprotective mechanism . However, unregulated excessive autophagy is known to be cytotoxic, inducing autophagic cell death . Results in this study indicate that monotherapy of erlotinib might induce cytoprotective autophagy but the combination of erlotinib and sorafenib causes cell-killing autophagy. Autophagy might be associated with cellular senescence whose role in antitumor therapy has not yet been established because the senescence in tumor cells can contribute to not only therapy resistance but also growth suppression or sensitization to therapy [60–62]. While these findings suggest the complex role of autophagy in tumor therapy, our observation might indicate the potential of combination therapies that induce autophagy-related remarkable cytotoxicity in tumor-initiating cells.
We found that growth of GICs in agarose gel correlates with the Gompertz function. Combination treatments that affect RTKs and the proximal downstream molecules RAF/MEK or PI3K induce synergistic cytotoxicity in different GICs that show diversity in drug sensitivity. Among these combination treatments, combination of erlotinib with sorafenib induces apoptotic and autophagic cell death associated with attenuated Akt and ERK signaling, decreased nuclear PKM2 and β-catenin and affected Bcl-2 family proteins in GSC11 cells in vitro, and tends to improve survival of the brain tumor xenografts. Although sensitivity of GIC lines to anticancer agents is diverse, the rationale for efficient combination therapy is thought to be inhibition of both RTK and RAF/MEK or PI3K. This rationale and the method in this study could be useful for personalized drug screening using GICs from each patient to optimize targeted therapy for glioma patients.
TS and OB contributed to the concept and the design of the study and drafted the manuscript. TS, LH, VH, QW and LAG acquired and analyzed the data. KRH performed the statistical analysis. KL, AEG, TD, GR, and NPF contributed to the study supervision. WP, YAW and JH contributed to the critical review and revision of the draft. LAG, FFL, LY, EPS, NPF and WP contributed to technical and material support. All authors read and approved the final manuscript.
This work was supported by the National Institute of Health/National Cancer Institute (R01 CA128991 to O.B. and N.P.F, CA16672 to RPPA Core Facility-Functional Proteomics). We thank Drs. N Mizushima, T Yoshimori and ND Tribolet for the generous gifts of experimental materials, and Ms. K Sherman (The University of Texas MD Anderson Cancer Center) for technical support.
Dr. Nicholas P Farrell and Dr. Waldemar Priebe hold patents and have financial interests in the development of compound BBR3610 and WP1066, respectively. The other authors declare that no conflict of interests exists.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Wen PY, Kesari S. Malignant gliomas in adults. New Engl J Med. 2008;359:492–507.View ArticlePubMedGoogle Scholar
- Purow B, Schiff D. Advances in the genetics of glioblastoma: are we reaching critical mass? Nat Rev Neurol. 2009;5:419–26.PubMed CentralView ArticlePubMedGoogle Scholar
- Brennan CW, Verhaak RG, McKenna A, Campos B, Noushmehr H, Salama SR, et al. The somatic genomic landscape of glioblastoma. Cell. 2013;155:462–77.PubMed CentralView ArticlePubMedGoogle Scholar
- Wick W, Weller M, Weiler M, Batchelor T, Yung AW, Platten M. Pathway inhibition: emerging molecular targets for treating glioblastoma. Neuro-oncology. 2011;13:566.PubMed CentralView ArticlePubMedGoogle Scholar
- Sottoriva A, Spiteri I, Piccirillo SG, Touloumis A, Collins VP, Marioni JC, et al. Intratumor heterogeneity in human glioblastoma reflects cancer evolutionary dynamics. Proc Natl Acad Sci USA. 2013;110:4009–14.PubMed CentralView ArticlePubMedGoogle Scholar
- Stommel JM, Kimmelman AC, Ying H, Nabioullin R, Ponugoti AH, Wiedemeyer R, et al. Coactivation of receptor tyrosine kinases affects the response of tumor cells to targeted therapies. Science. 2007;318:287–90.View ArticlePubMedGoogle Scholar
- Iadevaia S, Lu Y, Morales FC, Mills GB, Ram PT. Identification of optimal drug combinations targeting cellular networks: integrating phospho-proteomics and computational network analysis. Cancer Res. 2010;70:6704–14.PubMed CentralView ArticlePubMedGoogle Scholar
- Chang S, Kuhn J, Lamborn K, Cloughesy T, Robins I, Lieberman F, et al. Phase I/II study of erlotinib and temsirolimus for patients with recurrent malignant gliomas (MG)(NABTC 04-02). In: ASCO Annual Meeting Proceedings 2009. p. 2004.
- Singh SK, Clarke ID, Terasaki M, Bonn VE, Hawkins C, Squire J, et al. Identification of a cancer stem cell in human brain tumors. Cancer Res. 2003;63:5821–8.PubMedGoogle Scholar
- Franken NA, Rodermond HM, Stap J, Haveman J, Van Bree C. Clonogenic assay of cells in vitro. Nat Protoc. 2006;1:2315–9.View ArticlePubMedGoogle Scholar
- Friedrich J, Seidel C, Ebner R, Kunz-Schughart LA. Spheroid-based drug screen: considerations and practical approach. Nat Protoc. 2009;4:309–24.View ArticlePubMedGoogle Scholar
- Louis SA, Rietze RL, Deleyrolle L, Wagey RE, Thomas TE, Eaves AC, et al. Enumeration of Neural Stem and Progenitor Cells in the Neural Colony-Forming Cell Assay. Stem Cells. 2008;26:988–96.View ArticlePubMedGoogle Scholar
- Ying M, Wang S, Sang Y, Sun P, Lal B, Goodwin C, et al. Regulation of glioblastoma stem cells by retinoic acid: role for Notch pathway inhibition. Oncogene. 2011;30:3454–67.PubMed CentralView ArticlePubMedGoogle Scholar
- Vivanco I, Robins HI, Rohle D, Campos C, Grommes C, Nghiemphu PL, et al. Differential sensitivity of glioma-versus lung cancer–specific EGFR mutations to EGFR kinase inhibitors. Cancer Discov. 2012;2:458–71.PubMed CentralView ArticlePubMedGoogle Scholar
- Hu J, Ho AL, Yuan L, Hu B, Hua S, Hwang SS, et al. Neutralization of terminal differentiation in gliomagenesis. Proc Natl Acad Sci USA. 2013;110:14520–7.PubMed CentralView ArticlePubMedGoogle Scholar
- Barrett LE, Granot Z, Coker C, Iavarone A, Hambardzumyan D, Holland EC, et al. Self-renewal does not predict tumor growth potential in mouse models of high-grade glioma. Cancer Cell. 2012;21:11–24.View ArticlePubMedGoogle Scholar
- Billecke C, Finniss S, Tahash L, Miller C, Mikkelsen T, Farrell NP, et al. Polynuclear platinum anticancer drugs are more potent than cisplatin and induce cell cycle arrest in glioma. Neuro-oncology. 2006;8:215–26.PubMed CentralView ArticlePubMedGoogle Scholar
- Iwamaru A, Szymanski S, Iwado E, Aoki H, Yokoyama T, Fokt I, et al. A novel inhibitor of the STAT3 pathway induces apoptosis in malignant glioma cells both in vitro and in vivo. Oncogene. 2007;26:2435–44.View ArticlePubMedGoogle Scholar
- Bhat KP, Balasubramaniyan V, Vaillant B, Ezhilarasan R, Hummelink K, Hollingsworth F, et al. Mesenchymal differentiation mediated by NF-κB promotes radiation resistance in glioblastoma. Cancer Cell. 2013;24:331–46.View ArticlePubMedGoogle Scholar
- Saito N, Fu J, Zheng S, Yao J, Wang S, Liu DD, et al. A High Notch Pathway Activation Predicts Response to γ Secretase Inhibitors in Proneural Subtype of Glioma Tumor-Initiating Cells. Stem Cells. 2014;32:301–12.PubMed CentralView ArticlePubMedGoogle Scholar
- Kajiwara Y, Panchabhai S, Levin VA. A new preclinical 3-dimensional agarose colony formation assay. Technol Cancer Res Treat. 2008;7:329–34.View ArticlePubMedGoogle Scholar
- Laird AK. Dynamics of tumour growth. Br J Cancer. 1964;18:490.PubMed CentralView ArticleGoogle Scholar
- Kaufmann SH, Peereboom D, Buckwalter CA, Svingen PA, Grochow LB, Donehower RC, et al. Cytotoxic effects of topotecan combined with various anticancer agents in human cancer cell lines. J Natl Cancer Inst. 1996;88:734–41.View ArticlePubMedGoogle Scholar
- Jorgensen TJ, Tian H, Joseph IB, Menon K, Frost D. Chemosensitization and radiosensitization of human lung and colon cancers by antimitotic agent, ABT-751, in athymic murine xenograft models of subcutaneous tumor growth. Cancer Chemother Pharmacol. 2007;59:725–32.View ArticlePubMedGoogle Scholar
- Lal S, Lacroix M, Tofilon P, Fuller GN, Sawaya R, Lang FF. An implantable guide-screw system for brain tumor studies in small animals. J Neurosurg. 2000;92:326–33.View ArticlePubMedGoogle Scholar
- Mizushima N, Yamamoto A, Matsui M, Yoshimori T, Ohsumi Y. In vivo analysis of autophagy in response to nutrient starvation using transgenic mice expressing a fluorescent autophagosome marker. Mol Biol Cell. 2004;15:1101–11.PubMed CentralView ArticlePubMedGoogle Scholar
- Tjuvajev J, Muraki A, Ginos J, Berk J, Koutcher J, Ballon D, et al. Iododeoxyuridine uptake and retention as a measure of tumor growth. J Nucl Med. 1993;34:1152–62.PubMedGoogle Scholar
- Smith JS, Tachibana I, Passe SM, Huntley BK, Borell TJ, Iturria N, et al. PTEN mutation, EGFR amplification, and outcome in patients with anaplastic astrocytoma and glioblastoma multiforme. J Natl Cancer Inst. 2001;93:1246–56.View ArticlePubMedGoogle Scholar
- Agarwal S, Zerillo C, Kolmakova J, Christensen J, Harris L, Rimm D, et al. Association of constitutively activated hepatocyte growth factor receptor (Met) with resistance to a dual EGFR/Her2 inhibitor in non-small-cell lung cancer cells. Br J Cancer. 2009;100:941–9.PubMed CentralView ArticlePubMedGoogle Scholar
- Griffero F, Daga A, Marubbi D, Capra MC, Melotti A, Pattarozzi A, et al. Different response of human glioma tumor-initiating cells to epidermal growth factor receptor kinase inhibitors. J Biol Chem. 2009;284:7138–48.PubMed CentralView ArticlePubMedGoogle Scholar
- Byers LA, Diao L, Wang J, Saintigny P, Girard L, Peyton M, et al. An epithelial–mesenchymal transition gene signature predicts resistance to EGFR and PI3 K inhibitors and identifies Axl as a therapeutic target for overcoming EGFR inhibitor resistance. Clin Cancer Res. 2013;19:279–90.PubMed CentralView ArticlePubMedGoogle Scholar
- Verhaak RG, Hoadley KA, Purdom E, Wang V, Qi Y, Wilkerson MD, et al. Integrated genomic analysis identifies clinically relevant subtypes of glioblastoma characterized by abnormalities in PDGFRA, IDH1, EGFR, and NF1. Cancer Cell. 2010;17:98–110.PubMed CentralView ArticlePubMedGoogle Scholar
- Nakasu S, Nakasu Y, Fukami T, Jito J, Nozaki K. Growth curve analysis of asymptomatic and symptomatic meningiomas. J Neurooncol. 2011;102:303–10.View ArticlePubMedGoogle Scholar
- Collins VP. Observation on growth rates of human tumors. Am J Roentgenol. 1956;76:988–1000.Google Scholar
- Hart D, Shochat E, Agur Z. The growth law of primary breast cancer as inferred from mammography screening trials data. Br J Cancer. 1998;78:382.PubMed CentralView ArticlePubMedGoogle Scholar
- Guiot C, Degiorgis PG, Delsanto PP, Gabriele P, Deisboeck TS. Does tumor growth follow a “universal law”? J Theor Biol. 2003;225:147–51.View ArticlePubMedGoogle Scholar
- Martinelli E, Troiani T, Morgillo F, Rodolico G, Vitagliano D, Morelli MP, et al. Synergistic antitumor activity of sorafenib in combination with epidermal growth factor receptor inhibitors in colorectal and lung cancer cells. Clin Cancer Res. 2010;16:4990–5001.View ArticlePubMedGoogle Scholar
- Giovannetti E, Labots M, Dekker H, Galvani E, Lind J, Sciarrillo R, et al. Molecular mechanisms and modulation of key pathways underlying the synergistic interaction of sorafenib with erlotinib in non-small-cell-lung cancer (NSCLC) cells. Curr Pharm Des. 2013;19:927–39.View ArticlePubMedGoogle Scholar
- Zhang Z, Zhou X, Shen H, Wang D, Wang Y. Phosphorylated ERK is a potential predictor of sensitivity to sorafenib when treating hepatocellular carcinoma: evidence from an in vitro study. BMC Med. 2009;7:41.PubMed CentralView ArticlePubMedGoogle Scholar
- Ou D-L, Shen Y-C, Liang J-D, Liou J-Y, Yu S-L, Fan H-H, et al. Induction of Bim expression contributes to the antitumor synergy between sorafenib and mitogen-activated protein kinase/extracellular signal-regulated kinase kinase inhibitor CI-1040 in hepatocellular carcinoma. Clin Cancer Res. 2009;15:5820–8.View ArticlePubMedGoogle Scholar
- Ercan D, Xu C, Yanagita M, Monast CS, Pratilas CA, Montero J, et al. Reactivation of ERK signaling causes resistance to EGFR kinase inhibitors. Cancer Discov. 2012;2:934–47.PubMed CentralView ArticlePubMedGoogle Scholar
- Festuccia C, Gravina GL, Muzi P, Millimaggi D, Dolo V, Vicentini C, et al. Akt down-modulation induces apoptosis of human prostate cancer cells and synergizes with EGFR tyrosine kinase inhibitors. Prostate. 2008;68:965–74.View ArticlePubMedGoogle Scholar
- Zito CR, Jilaveanu LB, Anagnostou V, Rimm D, Bepler G, Maira S-M, et al. Multi-level targeting of the phosphatidylinositol-3-kinase pathway in non-small cell lung cancer cells. PLoS One. 2012;7:e31331.PubMed CentralView ArticlePubMedGoogle Scholar
- Glaysher S, Bolton L, Johnson P, Atkey N, Dyson M, Torrance C, et al. Targeting EGFR and PI3 K pathways in ovarian cancer. Br J Cancer. 2013;109:1786–94.PubMed CentralView ArticlePubMedGoogle Scholar
- Mazurek S, Boschek CB, Hugo F, Eigenbrodt E. Pyruvate kinase type M2 and its role in tumor growth and spreading. In Seminars in cancer biology, Elsevier; 2005. p. 300–308.
- Yang W, Xia Y, Ji H, Zheng Y, Liang J, Huang W, et al. Nuclear PKM2 regulates [bgr]-catenin transactivation upon EGFR activation. Nature. 2011;480:118–22.PubMed CentralView ArticlePubMedGoogle Scholar
- Yang W, Xia Y, Hawke D, Li X, Liang J, Xing D, et al. PKM2 phosphorylates histone H3 and promotes gene transcription and tumorigenesis. Cell. 2012;150:685–96.PubMed CentralView ArticlePubMedGoogle Scholar
- Jiang Y, Li X, Yang W, Hawke DH, Zheng Y, Xia Y, et al. PKM2 regulates chromosome segregation and mitosis progression of tumor cells. Mol Cell. 2014;53:75–87.PubMed CentralView ArticlePubMedGoogle Scholar
- Li C, Zhou C, Wang S, Feng Y, Lin W, Lin S, et al. Sensitization of glioma cells to tamoxifen-induced apoptosis by Pl3-kinase inhibitor through the GSK-3beta/beta-catenin signaling pathway. PLoS One. 2011;6.
- Peereboom DM, Ahluwalia MS, Ye X, Supko JG, Hilderbrand SL, Phuphanich S, et al. NABTT 0502: a phase II and pharmacokinetic study of erlotinib and sorafenib for patients with progressive or recurrent glioblastoma multiforme. Neuro-oncology. 2013;15:490–6.PubMed CentralView ArticlePubMedGoogle Scholar
- Perou CM, Sørlie T, Eisen MB, van de Rijn M, Jeffrey SS, Rees CA, et al. Molecular portraits of human breast tumours. Nature. 2000;406:747–52.View ArticlePubMedGoogle Scholar
- Noushmehr H, Weisenberger DJ, Diefes K, Phillips HS, Pujara K, Berman BP, et al. Identification of a CpG island methylator phenotype that defines a distinct subgroup of glioma. Cancer Cell. 2010;17:510–22.PubMed CentralView ArticlePubMedGoogle Scholar
- Akhavan D, Pourzia AL, Nourian AA, Williams KJ, Nathanson D, Babic I, et al. De-repression of PDGFRβ transcription promotes acquired resistance to EGFR tyrosine kinase inhibitors in glioblastoma patients. Cancer Discov. 2013;3:534–47.PubMed CentralView ArticlePubMedGoogle Scholar
- Shingu T, Yamada K, Hara N, Moritake K, Osago H, Terashima M, et al. Growth inhibition of human malignant glioma cells induced by the PI3-K-specific inhibitor. J Neurosurg. 2003;98:154–61.View ArticlePubMedGoogle Scholar
- Phillips HS, Kharbanda S, Chen R, Forrest WF, Soriano RH, Wu TD, et al. Molecular subclasses of high-grade glioma predict prognosis, delineate a pattern of disease progression, and resemble stages in neurogenesis. Cancer Cell. 2006;9:157–73.View ArticlePubMedGoogle Scholar
- Gwak H-S, Park HJ, Yoo H, Youn SM, Rhee CH, Lee SH. Chemotherapy for malignant gliomas based on histoculture drug response assay: a pilot study. J Korean Neurosurg Soc. 2011;50:426–33.PubMed CentralView ArticlePubMedGoogle Scholar
- Yung W-KA, Shapiro JR, Shapiro WR, et al. Heterogeneous chemosensitivities of subpopulations of human glioma cells in culture. Cancer Res. 1982;42:992–8.PubMedGoogle Scholar
- Spigel DR, Burris HA, Greco FA, Shipley DL, Friedman EK, Waterhouse DM, et al. Randomized, double-blind, placebo-controlled, phase II trial of sorafenib and erlotinib or erlotinib alone in previously treated advanced non–small-cell lung cancer. J Clin Oncol. 2011;29:2582–9.View ArticlePubMedGoogle Scholar
- Kondo Y, Kanzawa T, Sawaya R, Kondo S. The role of autophagy in cancer development and response to therapy. Nat Rev Cancer. 2005;5:726–34.View ArticlePubMedGoogle Scholar
- Ewald JA, Desotelle JA, Wilding G, Jarrard DF. Therapy-induced senescence in cancer. J Natl Cancer Inst. 2010;102:1536–46.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang M, Morsbach F, Sander D, Gheorghiu L, Nanda A, Benes C, et al. EGF receptor inhibition radiosensitizes NSCLC cells by inducing senescence in cells sustaining DNA double-strand breaks. Cancer Res. 2011;71:6261–9.PubMed CentralView ArticlePubMedGoogle Scholar
- Achuthan S, Santhoshkumar TR, Prabhakar J, Nair SA, Pillai MR. Drug-induced senescence generates chemoresistant stemlike cells with low reactive oxygen species. J Biol Chem. 2011;286:37813–29.PubMed CentralView ArticlePubMedGoogle Scholar