Inhibition of mutant EGFR in lung cancer cells triggers SOX2-FOXO6-dependent survival pathways
Abstract
Treatment of EGFR-mutant lung cancer with erlotinib results in dramatic tumor regression but it is invariably followed by drug resistance. In characterizing early transcriptional changes following drug treatment of mutant EGFR-addicted cells, we identified the stem cell transcriptional regulator SOX2 as being rapidly and specifically induced, both in vitro and in vivo. Suppression of SOX2 sensitizes cells to erlotinib-mediated apoptosis, ultimately decreasing the emergence of acquired resistance, whereas its ectopic expression reduces drug-induced cell death. We show that erlotinib relieves EGFR-dependent suppression of FOXO6, leading to its induction of SOX2, which in turn represses the pro-apoptotic BH3-only genes BIM and BMF. Together, these observations point to a physiological feedback mechanism that attenuates oncogene addiction-mediated cell death associated with the withdrawal of growth factor signaling and may therefore contribute to the development of resistance.
https://doi.org/10.7554/eLife.06132.001eLife digest
Tumors can form when cells gain mutations in genes that enable them to grow and divide rapidly. In some human lung cancers, genetic mutations are found in a gene that makes a protein called EGFR. This protein encourages cells to divide and the mutations can lead to the cancer cells producing more EGFR, or producing a form of the protein that is more active.
Treating these cancers with a drug called erlotinib inhibits EGFR and makes the tumors shrink dramatically, but the tumors will usually re-grow because any tumor cells that survive often become resistant to the drug. There are several ways that the tumor cells can become resistant, which makes the task of developing a solution to this problem more difficult. It has been suggested that the tumor cells may enter a temporary ‘drug-tolerant’ state that helps them to survive and makes it more likely that they will develop resistance to the drug. However, it is not clear how this drug-tolerant state might work.
To address this question, Rothenberg et al. examined which genes are switched on (or ‘expressed’) in tumor cells with a mutant version of EGFR after they were treated with the erlotinib drug. The experiments show that a gene called SOX2 is expressed in these cells. Cells that had lower levels of SOX2 expression were more sensitive to the effects of the drug and fewer cells developed resistance. On the other hand, cells that had higher levels of SOX2 expression were less sensitive to the drug and resistance was more likely to develop.
A protein called FOXO6—which is usually suppressed by EGFR—activates the SOX2 gene in these cells. Therefore, using erlotinib to inhibit EGFR to kill the cancer cells increases the activity of FOXO6, which in turn promotes the survival of some of the cells by activating the SOX2 gene. A better understanding of the ways in which cancer cells adapt to erlotinib and other drugs may help us to design more effective treatments with better outcomes for patients.
https://doi.org/10.7554/eLife.06132.002Introduction
The invariable development of drug resistance presents a critical challenge to the success of targeted cancer therapies (Jänne et al., 2005; O'Hare et al., 2006; Poulikakos and Rosen, 2011). Several mechanisms leading to such acquired resistance have been identified in patients with EGFR-mutant non-small cell lung cancer (NSCLC) treated with small molecule EGFR inhibitors such as erlotinib. Following dramatic initial tumor shrinkage, tumor regrowth is most frequently associated with the emergence of a secondary genetic change, the T790M ‘gatekeeper’ mutation within the EGFR kinase domain, which restores ATP binding in the presence of drug (Pao et al., 2004; Kobayashi et al., 2005; Yun et al., 2007). In other cases, amplification of related receptor tyrosine kinases (e.g., MET) or mutational activation of downstream kinases (e.g., BRAF, PIK3CA) may bypass the effect of EGFR inhibition (Engelman et al., 2007; Sequist et al., 2011; Ohashi et al., 2012). Different drug resistance mechanisms may coexist within different metastatic lesions of individual patients. Recently, clinical trials involving rebiopsy of tumor lesions at the earliest sign of drug resistance have also revealed phenotypic conversions that may contribute to drug resistance, including activation of epithelial-to-mesenchymal transition (EMT) and the remarkable trans-differentiation of lung cancers from adenocarcinoma to small cell histologies (Thomson et al., 2005; Witta et al., 2006; Sequist et al., 2011). While some well-defined resistance mechanisms, such as the T790M-EGFR gatekeeper mutation and MET amplification, may be addressed using second line targeted drugs, the plasticity of cancer cell adaptation to disrupted oncogenic signaling poses a major challenge to the long-term success of these promising therapies.
In vitro modeling of acquired resistance to EGFR inhibitors has raised the possibility that a transient so-called ‘drug-tolerant’ state may precede the development of mutationally defined, heritable drug resistance (Sharma et al., 2010). By analogy with bacterial models of antibiotic resistance, such an intermediate state may be unstable, but enable treated cells to survive in the presence of drug long enough to acquire mutations that ultimately confer sustained drug resistance (Balaban et al., 2004). In PC9 mutant, EGFR-addicted lung cancer cells, EGFR inhibition triggers apoptosis in the vast majority of cells in vitro, uncovering approximately 0.3% that are drug tolerant, quiescent, and expressing the stem cell marker CD133 and the histone H3K4 demethylase KDM5A (Sharma et al., 2010). These drug tolerant cells readily revert to a drug-sensitive state following removal of the EGFR inhibitor, and their emergence in vitro is suppressed by treatment with an EGFR inhibitor combined with inhibitors of either histone deacetylases (HDACs) or the IGF-1 receptor. While this intermediate resistance mechanism remains to be validated in the clinical setting, it raises the possibility of suppressing pre-conditions that favor the acquisition of drug resistance, in order to circumvent the challenge of treating multiple established drug-resistant pathways.
Beyond the selection of cancer cell populations with transient drug-resistant phenotypes, recent studies of targeted cancer drugs have defined more rapid signaling feedback loops that modulate the cellular response to growth factor inhibition. For instance, acute loss of ERK signaling triggered by RAF or MEK inhibitors in BRAF mutant melanoma cells relieves ERK-dependent inhibition of RAS and CRAF, whose activation through ErbB receptor signaling may lead to paradoxical proliferative signals (Pratilas et al., 2009; Paraiso et al., 2010; Lito et al., 2012). Similarly, in BRAF mutant colorectal cancers, feedback activation of EGFR-dependent signaling attenuates the consequences of mutant BRAF inhibition, suppressing the apoptotic effect (Corcoran et al., 2012; Prahallad et al., 2012). In addition to signaling feedback loops, transcriptional outputs that generally limit cell proliferation have also been implicated following disruption of EGFR activity, including the expression of transcriptional repressors, regulators of mRNA stability and microRNAs (Kobayashi et al., 2006; Amit et al., 2007; Avraham et al., 2010).
Here, we screened for early, unique transcriptional changes following erlotinib treatment in mutant EGFR-addicted cells, identifying highly specific induction of SOX2, a master transcriptional regulator required for embryonic stem cell maintenance. SOX2 represses the expression of pro-apoptotic molecules that mediate death following oncogene withdrawal in these cells. The induction of SOX2 results from the activation of FOXO6, a forkhead family transcription factor, following EGFR inhibition. Knockdown or ectopic expression of SOX2 modulates the degree of apoptosis observed following oncogene withdrawal and promotes drug resistance, pointing to a novel homeostatic mechanism that may contribute to cellular adaptation to the withdrawal of growth factor signaling, which underlies most approaches to targeted cancer therapy.
Results
SOX2 is specifically induced in EGFR-mutated lung cancer cells following treatment with the EGFR inhibitor erlotinib
To interrogate the transcriptional response to EGFR inhibition, we used HCC827 lung cancer cells, harboring an amplified mutated EGFR allele (in-frame deletion of 15 nucleotides in exon 19) and displaying exquisite sensitivity to the EGFR inhibitor erlotinib. Cell cultures were treated in triplicate with 1 µM erlotinib for 6 hr, followed by mRNA isolation and whole transcriptome analysis (Affymetrix U133 Plus 2.0 expression arrays) (Rothenberg, 2015). A total of 35 genes showed >fourfold change in expression (FDR <0.05), including 22 downregulated and 13 upregulated transcripts (represented by 48 unique probe sets; Figure 1—figure supplement 1A). Among induced transcripts, SOX2 was unique in the specificity and rapidity of its induction following EGFR inhibition (Figure 1, Figure 1—figure supplement 1B). Thus, SOX2 was strongly induced in three mutant EGFR-addicted lung cancer cell lines (HCC827, PC9, H3255) following treatment with physiologically relevant concentrations of erlotinib (0.1 µM), but not when these cells were treated with comparably effective doses of cytotoxic chemotherapy (Figure 1A,B and Figure 1—figure supplement 2A). SOX2 was also not induced in other oncogene-dependent models, such as ALK-translocated lung cancer cells treated with crizotinib, HER2-amplified breast cancer cells exposed to lapatinib or BRAF-mutant melanoma cells treated with AZD6244 (Figure 1A and Figure 1—figure supplement 2B). Consistent with its dependence on suppression of mutant EGFR signaling in the context of EGFR ‘addiction’, SOX2 was not induced following erlotinib treatment of H1975 cells, which harbor both an EGFR activating mutation and the T790M gatekeeper mutation that confers resistant to erlotinib; or in H1650 cells with mutated EGFR that are relatively resistant to the effects of EGFR inhibition in part through genetic loss of PTEN (Figure 1—figure supplement 2B) (Sos et al., 2009). However, treatment of H1975 cells with the L858R/T790M mutation-selective inhibitor WZ4002 resulted in SOX2 induction (Figure 1—figure supplement 2B, right) (Zhou et al., 2009). In cells that show erlotinib-mediated induction of SOX2, siRNA-mediated knockdown of EGFR also led to strong induction of SOX2 (in the absence of erlotinib), confirming the specificity of the drug effect (Figure 1C). Simultaneous treatment of cells with actinomycin D and erlotinib suppressed the induction of SOX2, consistent with a primary effect of EGFR inhibition in increasing SOX2 transcript levels (Figure 1—figure supplement 2C).
Other transcripts induced or repressed following erlotinib treatment of mutant EGFR-addicted cells were not selective to EGFR signaling. Downregulated genes included known direct transcriptional targets of ERK signaling (CCND1, FOSL1, EGR1, IER3, IL-8) and shared feedback inhibitors of receptor tyrosine kinase (RTK) signaling (DUSP6) (Amit et al., 2007). This gene set overlaps with genes known to be differentially expressed following treatment of BRAF-mutant melanoma cells with a MEK inhibitor and exposure of EGFR-mutant lung cancer cells to an irreversible EGFR inhibitor (Figure 1—figure supplement 3) (Kobayashi et al., 2006; Pratilas et al., 2009). In addition to SOX2, the 12 other transcripts induced by EGFR inhibition included genes encoding metabolizing enzymes (CYP1B1, CYP1A1) that are normally induced by treatment with a variety of small chemical entities, genes that we also found to be equally well induced by inhibition of downstream signaling pathways (MEK/MAPK inhibition using AZD6244 and PI3K/mTOR inhibition using BEZ235), one transcript (CCNG2) previously described in another EGFR model (Kobayashi et al., 2006) and a long noncoding RNA (NEAT1). Taken all together, these results indicate that suppression of EGFR signaling in mutant EGFR-addicted lung cancer cells is highly specific in triggering transcriptional induction of SOX2.
Induction of SOX2 in erlotinib-treated cells
SOX2 encodes a master transcriptional regulator, implicated in stem cell maintenance and iPS cell generation. It is also required for upper aerodigestive tract development, and is known to be amplified in a subset of esophageal and squamous lung cancers, although it has not been previously implicated in lung adenocarcinomas, including the subset driven by mutant EGFR (Ellis et al., 2004; Gontan et al., 2008; Bass et al., 2009). Remarkably, SOX2 expression following exposure of HCC827 cells to erlotinib was transient, peaking at 24 hr after exposure to therapeutic levels of the drug (Figure 1D). Thereafter, SOX2 expression returned to basal levels despite continued erlotinib treatment in surviving cells (Figure 1D, Figure 1—figure supplement 4).
The level of SOX2 induction in cultured cells exposed to erlotinib showed considerable heterogeneity, with a subset of cells (∼20%, with some experimental variability) expressing high levels (Figure 2A). The SOX2+ fraction was not increased by higher drug dosage, beyond that required for full inhibition of EGFR (Figure 2—figure supplement 1). Given the link between SOX2 expression and cellular reprogramming, we first asked whether cells with the high SOX2 expression represent a subset with stem cell markers. However, SOX2 expression did not correlate with expression of the putative stem cell markers CD133, CD44, CD24, OCT-4, or KLF-4 (Figure 2—figure supplement 2) nor did microarray-based expression profiling of high SOX2-sorted cells identify a stem-like signature (data not shown). Nonetheless, SOX2-expressing cells had a very low proliferative index, as measured by Ki67 staining (0.5% of Ki67+/SOX2+ vs 51% Ki67+/SOX2− HCC827 cells at baseline [p = 0.015]; and 0.15% Ki67+/SOX2+ vs 6.4% Ki67+/SOX2− cells following erlotinib [p < 0.0001]) (Figure 2C, Figure 2—figure supplement 3).
-
Figure 2—source data 1
- https://doi.org/10.7554/eLife.06132.009
-
Figure 2—source data 2
- https://doi.org/10.7554/eLife.06132.010
-
Figure 2—source data 3
- https://doi.org/10.7554/eLife.06132.011
-
Figure 2—source data 4
- https://doi.org/10.7554/eLife.06132.012
-
Figure 2—source data 5
- https://doi.org/10.7554/eLife.06132.013
-
Figure 2—source data 6
- https://doi.org/10.7554/eLife.06132.014
To test whether the heterogeneous induction of SOX2 following EGFR inhibition represents a stochastic event or a heritable property shared by a subset of the parental population, we treated cells sequentially with pulses of erlotinib. Retreatment of cells after a period of recovery produced similar heterogeneity of SOX2+ cells as the initial treatment, pointing to the absence of enrichment for highly inducible cells (Figure 2B and Figure 2—figure supplement 4A). In addition, SOX2 could still be induced in cells made resistant to erlotinib, but only at the much higher doses of drug required to fully inhibit EGFR in resistant cells (Figure 2—figure supplement 4B). Finally, cloning of 5–6 individual HCC827 and PC9 cells consistently generated mixed populations, including high level SOX2 inducers together with non-expressing cells, demonstrating that this heterogeneity is likely stochastic, rather than heritable (Figure 2C and Figure 2—figure supplement 5). Thus, erlotinib treatment of EGFR-mutant cells results in transient and heterogeneous induction of SOX2, with a stochastic distribution, integrally tied to inhibition of EGFR, in which the cells with the highest expression have a low proliferative index.
Induction of SOX2 in EGFR-mutant tumors following erlotinib treatment
To test the physiological significance of SOX2 induction following withdrawal of mutant EGFR signaling, we first made use of mouse tumor models. The effectiveness of EGFR inhibitors in treating patients with EGFR-mutant NSCLC is well modeled in mouse xenograft assays, where oral administration of erlotinib for a few days is sufficient to cause massive regression of established tumors. We generated PC9 cell-derived subcutaneous tumors in nude mice and treated these with a single oral dose of 100 mg/kg erlotinib when the tumors had reached approximately 500 mm3, harvesting tumors 24 hr after treatment. Immunohistochemical (IHC) analysis revealed minimal SOX2 expression in mock-treated xenografts (mean 0.2 SOX2+ nuclei per field), but clearly increased (and heterogeneous) SOX2 positive cells after a single dose of erlotinib (mean: 7.4 SOX2+ nuclei per field, N = 147–151 fields, p < 0.0001) (Figure 3). Similar studies in HCC827 cell-derived xenografts revealed low (but detectable) levels of SOX2 expression in mock-treated tumors; again, a single oral dose of erlotinib increased both the number of SOX2-positive cells and the level of SOX2 expression per nucleus (Figure 3—figure supplement 1). Thus, in a physiological setting that mimics the initial therapeutic response to EGFR inhibitors in EGFR-mutant NSCLC, treated cancer cells rapidly induce SOX2.
-
Figure 3—source data 1
- https://doi.org/10.7554/eLife.06132.022
-
Figure 3—source data 2
- https://doi.org/10.7554/eLife.06132.023
Induction of SOX2 in patient-derived EGFR-mutant tumor cells
While established cancer cell lines, such as PC9 and HCC827, recapitulate the phenomenon of oncogene addiction, patient-derived cell lines directly cultured from biopsies may be more representative of heterogeneous primary cultures (Crystal et al., 2014). Such biopsies are typically obtained at the time of disease progression, where defining drug resistance mechanisms may shape further therapy. We therefore analyzed short-term cultures of EGFR-mutant cells derived by re-biopsy of two patients who had initially responded to erlotinib therapy but subsequently developed progressive disease due to the acquisition of a T790M gatekeeper mutation. SOX2 induction was absent following treatment with erlotinib, which was ineffective in inhibiting EGFR in these resistant patient-derived cells (Figure 4). However, the novel ‘third line’ irreversible, EGFR-mutant-specific inhibitor WZ4002 demonstrated potent EGFR inhibition in these cells, along with induction of SOX2 (Figure 4). Thus, SOX2 induction is consistently observed following acute withdrawal of EGFR signals in cancer models as well as in patient-derived cells that exhibit oncogene dependence on the EGFR pathway.
-
Figure 4—source data 1
- https://doi.org/10.7554/eLife.06132.026
Regulation of anti-apoptotic signals by SOX2
To examine the consequence of ectopic SOX2 expression, we generated HCC827 cells with a lentiviral-driven, doxycycline-inducible construct. Careful titration of doxycycline-permitted induction of ectopic SOX2 to physiologic levels in the absence of erlotinib, comparable at the single cell level to endogenous SOX2 induction in the presence of erlotinib, though in a larger fraction of cells (Figure 5A and Figure 5—figure supplement 1A). As expected, treatment of control HCC827 cells with erlotinib led to the inhibition of EGFR signaling, as measured by decreased phosphorylation of EGFR, AKT, and ERK, and caused dramatic apoptosis, as determined by increased PARP and caspase-3 cleavage (Figure 5A, lanes 1–2). In contrast, expression of ectopic SOX2 significantly decreased erlotinib-mediated apoptosis, resulting in decreased PARP and caspase-3 cleavage (Figure 5A, lane 2 vs 4). Exogenous SOX2 itself did not alter basal or erlotinib-inhibited phosphorylation of EGFR, AKT, or ERK (Figure 5A, lanes 3–4).
-
Figure 5—source data 1
- https://doi.org/10.7554/eLife.06132.028
-
Figure 5—source data 2
- https://doi.org/10.7554/eLife.06132.029
To determine the functional significance of endogenous SOX2 induction, we next used siRNAs to block its induction in erlotinib-treated cells. While both HCC827 and PC9 cells are highly sensitive to EGFR inhibition at baseline, SOX2 knockdown further increased erlotinib-induced apoptosis, as determined by PARP and caspase-3 cleavage assays and by cell enumeration (Figure 5B,D and Figure 5—figure supplement 2). The apoptotic effect of the most potent siRNA was rescued by expression of an ectopic siRNA-resistant SOX2 construct (Figure 5C) and individual siRNAs-induced apoptosis in proportion to the degree of knockdown (Figure 5—figure supplement 3), confirming the specificity of the effect. As with expression of exogenous SOX2, endogenous SOX2 suppression itself did not have a consistent effect on EGFR signaling, as measured by phosphorylation of EGFR, AKT, or ERK (Figure 5B).
Given the effect of transient SOX2 induction in acutely modulating cell survival following erlotinib treatment of EGFR-mutant cells, we tested whether this translates into a longer term impact on acquired drug resistance. We used the well-characterized PC9 cell model of EGFR-mutant NSCLC to transfect siRNAs against SOX2 (the same duplex validated by rescue in Figure 5C) or control, followed by continuous exposure to 1.0 µM erlotinib (Figure 5E). As expected, SOX2 knockdown alone had no significant toxicity, while erlotinib treatment led to massive cell death at 3 days (Figure 5E, left and middle images). At high magnification, many more individual surviving cells remained in control transfected cells than after transfection with siRNA targeting SOX2 (Figure 5E, middle insets). With continued incubation in erlotinib, drug-resistant colonies emerged in control transfected cultures at higher frequency than in those treated with siRNA targeting SOX2 (Figure 5E, right). Consistent with the short-term duration of siRNA effectiveness, these results suggest that preventing the short-term induction of SOX2 allows fewer cells to survive the initial exposure to erlotinib, preventing an adaptation response that ultimately delays the emergence of erlotinib-resistant colonies.
Erlotinib-induced SOX2 directly regulates expression of BIM and BMF
The mechanisms underlying the cell death response of cancer cells following withdrawal of oncogene-addicting signals are complex but appear to result in part from activation of the pro-apoptotic proteins BIM and PUMA (Costa et al., 2007; Gong et al., 2007; Bean et al., 2013). To search for additional targets that mediate SOX2's anti-apoptotic effect, we screened for altered expression of multiple pro-apoptotic and anti-apoptotic proteins in cells with overexpression of ectopic SOX2, following treatment with erlotinib. In both HCC827 and PC9 EGFR-mutant lung cancer cells, three of 14 BH3-domain containing proteins tested showed strongly reduced mRNA induction by erlotinib in the presence of ectopic SOX2: BIM, BMF, and HRK (but not PUMA) (Figure 6A and Figure 6—figure supplement 1). All three of these are known to induce apoptosis to different degrees, depending on the type of apoptotic stimulus. The effect of ectopic SOX2 on BIM protein levels was confirmed using two different expression strategies (Figure 5—figure supplement 1B,C). Ectopic expression of SOX2, both at high and physiological levels, attenuated the erlotinib-induced induction of SOX2. In contrast, SOX2 knockdown coincident with erlotinib treatment led to an overall increase in BIM levels, although this effect was attenuated when averaged across a cell population with heterogeneous induction of endogenous SOX2 (Figure 5—figure supplement 4).
To determine whether BIM and BMF are endogenous SOX2 target genes, we carried out chromatin immunoprecipitation sequencing (ChIP seq) experiments, using erlotinib-treated HCC827 cells. Strong SOX2-binding peaks were present within ∼2 kB of the transcriptional start sites (TSS) of both BIM and BMF genes (for BMF, the peak spans the TSS; for BIM, it is located within the first intron) (Figure 6B, upper). Importantly, ChIP qPCR demonstrated significant enhancement of SOX2 binding with erlotinib treatment at both peaks compared to untreated cells, suggesting that binding is functionally significant (Figure 6B, lower).
Consistent with their functional roles, knockdown of either BIM or BMF (but not HRK) decreased erlotinib-induced apoptosis, with combined BIM and BMF knockdown displaying an additive effect (Figure 6C). The role of BMF could be even more important, since BMF knockdown increases BIM levels, which may blunt the anti-apoptotic effect of BMF loss (Figure 6C, left). Together, these results suggest that SOX2 induction following erlotinib exposure in EGFR-mutant cells suppresses transcriptional induction of the BH3-only BIM and BMF genes, which contribute to apoptosis following oncogene withdrawal.
Induction of SOX2 following EGFR inhibition is regulated by FOXO6
To search for mediators of SOX2 induction, we explored the Molecular Signatures and TRANSFAC databases for transcription factor target sequences within the promoters of the 12 highest erlotinib-induced genes (Wingender et al., 2000; Subramanian et al., 2005). Several binding motifs for FOXO proteins were highly significantly enriched (q-value = 0.003 or less): for SOX2, multiple sites were present within 2 kb of the transcriptional start site (Figure 7A and Figure 7—figure supplement 1). Expression of all of the FOXO family members was detectable at baseline in HCC827 cells and erlotinib treatment (8 hr) was associated with a 1.6–4.4-fold induction (Figure 7B), as well as with loss of the AKT-mediated inhibitory N-terminal threonine phosphorylation of the FOXO proteins (Figure 7—figure supplement 2A).
-
Figure 7—source data 1
- https://doi.org/10.7554/eLife.06132.037
-
Figure 7—source data 2
- https://doi.org/10.7554/eLife.06132.038
Given the evidence of erlotinib-mediated FOXO activation and its potential regulation of SOX2, we tested the consequence of siRNA-mediated knockdown of each gene family member, alone and in combination. Knockdown of FOXO6 using pooled siRNA constructs, but not the other FOXO proteins (individually or simultaneously), dramatically reduced erlotinib-mediated induction of SOX2 (Figure 7B,C and Figure 7—figure supplement 2B). The effect of FOXO6 knockdown on SOX2 was evident using multiple individual siRNAs targeting FOXO6 in both HCC827 and PC9 cells (Figure 7D; effect on other FOXO isoforms is shown in Figure 7—figure supplement 2C). Although some individual siRNAs targeting FOXO6 had off target effects on other FOXOs, direct targeting of FOXOs 1, 3a, and 4 had no effect on SOX2 expression, further supporting the specificity of the FOXO6 effect (Figure 7B). The effect of FOXO6 knockdown on SOX2 was not associated with any consistent effect on other aspects of EGFR signaling, although a moderate decrease in phospho-EGFR was observed in some experiments without significant differences in phospho-AKT or phospho-ERK (Figure 7E). Notably, FOXO6 expression was also heterogeneous and partially colocalized with SOX2 expression among populations of both untreated and treated cells (Figure 7—figure supplement 3). FOXO6 differs from other FOXO proteins in that even the inactive protein is localized in the nucleus (Jacobs et al., 2003; van der Heide et al., 2005). However, it shares the FOXO protein AKT-dependent inhibitory phosphorylation, whose suppression following repression of mutant EGFR signaling may in part explain the erlotinib-mediated FOXO6 activation.
Erlotinib resistance in SOX2-expressing EGFR-mutant cells
Virtually all EGFR-mutant lung cancer cell lines established from patients who have not been treated with erlotinib are highly sensitive to this drug, although a few cells lines appear to be intrinsically resistant. The HCC2935 human lung cancer cell line is remarkable for harboring a characteristic oncogene-addicting EGFR mutation (exon 19 deletion) yet having unexplained resistance to erlotinib (Figure 8A), including absence of the common T790M gatekeeper mutation within EGFR and no amplification of the MET bypass signaling pathway (Zheng et al., 2011). Notably, SOX2 expression at baseline is detectable in 90% of HCC2935 cells (compared to 3% of HCC827 and <1% PC9 cells), and its expression per cell is further increased upon EGFR inhibition (Figure 8B). To test whether increased SOX2 contributes to decreased erlotinib sensitivity in HCC2935 cells, we knocked down SOX2 using siRNA. A striking increase in erlotinib cytotoxicity was evident following SOX2 suppression (IC50 0.8 µM for siCTRL cells with no further toxicity up to 10 µM, IC50 0.1 µM for siSOX2 cells), associated with higher levels of BIM and increased PARP and caspase-3 cleavage (Figure 8A,C). Thus, increased baseline SOX2 contributes to erlotinib resistance in these EGFR-mutant cells.
-
Figure 8—source data 1
- https://doi.org/10.7554/eLife.06132.044
Discussion
The dramatic responsiveness of EGFR-mutated NSCLC to small molecule inhibitors of the EGFR kinase has provided a paradigm for the targeted therapy of epithelial cancers and has established a new standard of care for a genetically defined subset of patients with lung cancer. However, the invariable development of drug resistance greatly limits the effectiveness of this therapy, despite efforts to circumvent acquired genetic abnormalities, such as the recurrent T790M-EGFR gatekeeper mutation or MET amplification. Understanding how immediate signaling feedback loops modulate the cellular response to targeted inhibitors and how cellular heterogeneity may lead to transient drug tolerant states may thus provide important therapeutic opportunities.
FOXO6 regulates SOX2 expression
The regulation of SOX2 expression by FOXO6, whose activation is normally repressed by EGFR signaling, is consistent with the critical role played by FOXO proteins as integrators of cellular signaling pathways. The best studied isoform, FOXO1, is highly expressed in ES cells and has been implicated in the maintenance of pluripotency through activation of SOX2 transcription (Zhang et al., 2011). All FOXO proteins bind to similar DNA sequences, with isoform-specific activity presumably conferred by cellular and promoter context (Furuyama et al., 2000). Indeed, all four FOXO isoforms are expressed in EGFR-mutant lung cancer cells, transcriptionally induced following EGFR inhibition and phosphorylated at homologous Serine residues, yet only FOXO6 regulates SOX2 in these cells. Given variability in knockdown efficacy (Figure 7B, column 6), we cannot exclude some contribution from the other FOXO family members on SOX2 expression, but FOXO6 has the dominant effect in the cells tested. Activation of FOXO6 may occur through both AKT dependent and independent pathways (Jacobs et al., 2003; van der Heide et al., 2005), and indeed we observed that treatment with PI3K inhibitors alone is insufficient for induction of SOX2 (Figure 1B). Other pathways that have been implicated in SOX2 regulation in the developing lung, including FGF10, WNT/beta-Catenin signaling, and TTF1 (Que et al., 2007; Gontan et al., 2008; Hashimoto et al., 2012), had relatively modest effects on its induction following erlotinib treatment in EGFR-mutant cancer cells (Figure 7—figure supplement 4), pointing to FOXO6 as the dominant pathway in this model of oncogene-dependent signaling.
SOX2 regulates apoptosis through BIM and BMF
Our study extends the pro-apoptotic signals implicated in withdrawal of mutant EGFR signaling to include BMF, in addition to BIM and PUMA (Gong et al., 2007; Bean et al., 2013). In contrast to the latter, BMF binds with significant affinity to a subset of BCL-2 family members (BCL-2, BCL-xL, and BCL-w) (Chen et al., 2005; Kuwana et al., 2005), yet it clearly contributes to the apoptotic response to erlotinib (Figure 6). The direct binding by SOX2 of the BIM and BMF genes is consistent with the ability of the four reprogramming factors (SOX2, OCT-4, KLF-4, and c-MYC) to bind to the promoters of several anti-apoptotic genes (including BMF) which are induced early during the reprogramming process (Kim et al., 2008; Soufi et al., 2012). Interestingly, control of SOX2 expression and its modulation of apoptosis may differ in EGFR-mutant lung cancer, compared with other forms of NSCLC. In a recent study of lung cancers with wild-type EGFR, the activation of EGFR upregulated SOX2, thereby decreasing apoptosis through BCL2L1, a phenomenon that was notably absent in EGFR-mutant lung cancers (Chou et al., 2013). A similar pathway was reported in a prostate cancer cell line, also with wild-type EGFR (Rybak and Tang, 2013). Thus, an EGFR-SOX2-BCL2L1 pathway may be implicated in cancer cells with wild-type EGFR, whereas the EGFR-FOXO6-SOX2-BIM/BMF pathway we describe is specific for cells that are dependent on mutant EGFR signals for their survival.
Heterogeneity of SOX2 induction and implications for targeted therapy
The heterogeneous induction of SOX2 within a clonally derived cell population could reflect the existence of an intrinsic subpopulation with heritable traits, as recently proposed in medulloblastoma and skin carcinoma (Boumahdi et al., 2014; Vanner et al., 2014). Alternatively, it could result from stochastic variation between cells in the activity of cellular signaling pathways. We favor the latter model, given the absent coexpression of SOX2 with putative stem cell markers, the failure to enrich for SOX2 positive cells following repeated erlotinib treatments, and the regeneration of heterogeneous SOX2 inducibility following single cell cloning experiments. A stochastic cell killing model for TRAIL-induced apoptosis of HeLa cells was recently described, in which naturally occurring cell-to-cell variation in the levels or activity of upstream signaling proteins leads to differential induction of mitochondrial membrane permeability and apoptosis only in some cells (Spencer et al., 2009).
Critically, the ability of siRNA targeting SOX2 to substantially decrease the number and rate at which resistant subclones of EGFR-mutant cells emerge following continuous erlotinib treatment suggests that, despite its heterogeneous, transient, and stochastic expression, SOX2 contributes to the emergence of stably acquired resistance (Figure 5E). In these respects, our observations are reminiscent of transient drug-tolerant persister cells (DTPs), also observed following erlotinib treatment (Sharma et al., 2010). However, there are also significant differences between the mechanisms underlying DTPs and those observed here. DTPs emerge following longer erlotinib exposure at much higher concentrations, and they constitute a lower percentage of cells within the population. Their expression of the stem cell marker CD133 and the chromatin remodeling protein KDM5A was not reproduced by SOX2+ cells (Figure 2—figure supplements 2, 6). Furthermore, DTPs are detectable in many cancer cell lines following treatment with multiple cytotoxic and targeted agents, whereas induction of SOX2 appears to be strictly limited to targeted EGFR inhibition in cells addicted to mutant EGFR signaling. SOX2 induction may thus be a specific and early signaling response to EGFR withdrawal, which contributes to increased cell survival, thus enhancing the likelihood of epigenetic events leading to DTPs and ultimately to stable genetic mechanisms of acquired drug resistance. Our results with HCC2935 cells (Figure 8) further suggest that in a subset of cases, high basal SOX2 may even blunt the initial signaling response to EGFR inhibitors.
In summary, the rewiring of cellular signaling pathways driven by a dominant mutationally activated kinase underlies oncogene addiction in cancer, providing powerful opportunities for targeted therapy. At the same time, feedback pathways that attenuate these signaling readouts may play a major role in enhancing the development of drug resistance. The EGFR-FOXO6-SOX2 signaling pathway regulating expression of the BIM and BMF apoptotic factors thus identifies a feedback loop that may attenuate the effectiveness of anti-EGFR therapy in cancer and contribute toward the ultimate development of drug resistance (Figure 9). The contribution of key embryonic regulators such as SOX2 points to the conservation of critical developmental pathways in cancer cells, which modulate their response to the disruption of oncogenic signals.
Materials and methods
Cell lines and reagents
Request a detailed protocolAll cell lines were grown in RPMI (GIBCO) with 10% FBS and were obtained from the ATCC or the Massachusetts General Hospital Center for Molecular Therapeutics, which performs routine cell line authentication testing by SNP and STR analysis. Erlotinib, AZD6244, BEZ235, crizotonib, lapatinib, WZ4002 (Selleckchem, Houston, TX), and BIO (Sigma, St. Louis, MO) were dissolved in DMSO.
Microarray analysis
Request a detailed protocolHCC827 cells were plated in triplicates and treated the following day for 6 hr with DMSO, erlotinib, AZD6244, and BEZ235 (each 1 µM) or for 0, 3, 6, 12, and 24 hr with erlotinib (single samples for each time point). Total mRNA was isolated using the RNeasy Mini Kit (Qiagen). Generation of cRNA, hybridization to GeneChip Human Genome U133 Plus 2.0 mRNA expression arrays and array scanning were done according to the manufacturer's standard protocols (Affymetrix, Inc.). Raw Affymetrix CEL files were converted to a single value for each probe set using Robust Multi-array Average (RMA) and normalized using quantile normalization. Quality control was performed using the distribution analysis, correlation and principle variance components analysis functions in Jmp Genomics (SAS Institute). Individual genes with statistically significantly altered expression after treatment (compared to untreated cells) were identified using ANOVA after model adjustment for multiple hypothesis testing across LSMeans differences using the False Discovery Method of Benjamini and Hochberg (FDR) with Alpha set to 0.05. The complete data set is available at the NCBI Gene Expression Omnibus (http://www.ncbi.nlm.nih.gov/geo/), Accession GSE51212.
Quantitative PCR
Request a detailed protocolFor RT-qPCR, 1 µg mRNA was converted to cDNA using the First-strand cDNA Synthesis Kit (GE Healthcare, Pittsburgh, PA). cDNA was analyzed on a 7500 Real Time PCR System (Applied Biosystems) using TaqMan Gene Expression Master Mix and TaqMan gene expression assays (with GAPDH or ACTB as control, individual assays listed in Supplementary file 1—Thermo Fisher Scientific, Grand Island, NY). For ChIP-qPCR, immunoprecipitated chromatin (see below) was analyzed on the same system, using paired DNA PCR primers (Supplementary file 1) and Power SYBR Green PCR Master Mix (Applied Biosystems).
Immunoblot analysis and antibodies
Request a detailed protocolImmunoblotting was performed using standard methods. After treatment with the indicated drugs, cells were washed with cold PBS and lysed in buffer containing 20 mM Tris pH 7.5, 150 mM NaCl, 100 mM MgCl2, 1% Nonidet P-40 and 10% glycerol supplemented with HALT protease and phosphatase inhibitor cocktail (Thermo Fisher Scientific) using a Q800R sonicator (Qsonica, Newtown, CT). Lysates were centrifuged at 16,000×g for 5 min at 4°C. Protein concentrations were determined by BCA assay (Thermo Fisher Scientific). Proteins were resolved by SDS-PAGE and transferred to a polyvinylidenes difluoride membranes (Biorad) using the Transblot Turbo Transfer System (Biorad, Hercules, CA). Immunoblotting was performed per each antibody manufacturer's specifications. Antibodies used were pEGFR (Y1068), EGFR, pAKT (S473), AKT, pERK (T202/Y204), ERK, cleaved PARP, cleaved caspase-3, SOX2, BIM, phospho-FOXO1 (T24)/FoxO3a (T32), FOXO1, FOXO3 (Cell Signaling Technology, Beverly, MA), GAPDH (EMD Millipore, Billerica, MA), Ki67 (Epitomics, Burlingame, CA), phospho-FOXO6 (S184) (Abcam, Cambridge, MA), FOXO6 (Proteintech, Chicago, IL), CD133, GKLF, OCT4, MYC (Santa Cruz, Dallas, TX), CD44 and CD24 (BD Biosciences, San Jose, CA), and TY1 (Diagenode, Denville, NJ).
Chromatin immunoprecipitation
Request a detailed protocolChIP assays were carried out using approximately 5–10 × 106 HCC827 cells, following the procedures described previously (Mikkelsen et al., 2007; Ku et al., 2008). In brief, chromatin from formaldehyde-fixed cells was fragmented to a size range of 200–700 bases with a Branson 250 sonifier. Solubilized chromatin was immunoprecipitated overnight with goat anti-SOX2 antibody or goat IgG as a negative control (both antibodies are from R&D Systems, Minneapolis, MN). Antibody–chromatin complexes were pulled down with protein G-Dynabeads (Thermo Fisher Scientific), washed, and then eluted. After crosslink reversal, RNase A, and proteinase K treatment, immunoprecipitated DNA was extracted with the Agencourt AMPure XP PCR Purification Kit (Beckman Coulter, Brea, CA). ChIP DNA was quantified with Qubit (Thermo Fisher Scientific).
ChIP seq
Request a detailed protocol5 ng purified DNA (immunoprecipiated chromatin and input controls) were used to prepare Illumina compatible sequencing libraries for sequencing using the MiSeq Desktop Sequencer (Illumina, San Diego, CA). ChIP seq reads were aligned to the hg19 reference genome using BWA (Li and Durbin, 2009). Aligned reads were extended to 200 bp to approximate fragment sizes, and then 25-bp resolution density maps were derived by counting the number of fragments overlapping each position, using IGV tools (Robinson et al., 2011). The density maps were normalized to 5 million reads, and IGV was used to visualize ChIP seq coverage maps (Thorvaldsdottir et al., 2013).
Quantitative immunofluorescence analysis/immunohistochemistry
Request a detailed protocolFor immunofluorescence (IF) analysis, cells plated in chamber slides were fixed with 4% formaldehyde, permeabilized/blocked with 5% normal goat serum/0.3% Triton X-100 and then incubated overnight at 4°C with antibody to SOX2 (either rabbit—Cell Signaling Technology or goat—R&D systems). SOX2 staining was visualized with appropriate Alexa Fluor conjugated secondary antibodies (Jackson Immunoresearch, West Grove, PA). Ki-67 costaining was performed with antibody to Ki-67 conjugated to Alexa Fluor 488 (Epitomics). Nuclei were visualized with DAPI. Immunohistochemistry (IHC) of SOX2 on formalin-fixed, paraffin embedded tumor tissue was performed according to the antibody manufacturer's suggested protocol by the Specialized Histopathology Laboratory of the Massachusetts General Hospital using the SignalStain Boost IHC detection reagent (Cell Signaling Technology). All IF/IHC samples were imaged and quantitated using the Vectra Automated Multispectral Imaging System (PerkinElmer, Waltham, MA). Images were initially scanned at 4× magnification and then multiple high-powered fields were automatically acquired. The emission spectra of each fluorophore/IHC stain was computationally unmixed by preparing matched single stained samples. Unmixed images were segmented to identify individual cell nuclei based on DAPI or hematoxylin signal, and the mean nuclear signal for each fluorophore/IHC stain was calculated using inForm Advanced Image Analysis Software (PerkinElmer). For some samples, in addition to analyzing the mean signal for each stain on every cell analyzed, individual nuclei were scored as positive or negative for SOX2 using the scoring function in inForm to set a fixed threshold for SOX2 signal based on the presence or absence of any degree of SOX2 staining. For fluorescence images, when samples to be compared were acquired using different exposure times, data were normalized to exposure time using the normalized counts setting/function in inForm.
Toxicity assays
Request a detailed protocolCells were plated in multiple replicate wells at 2500 cells per well in 96-well format, treated 24 hr after plating and analyzed 72 hr later. Wells were fixed with 4% formaldehyde and stained with Syto60 red fluorescence nucleic acid stain (Thermo Fisher Scientific) for 1 hr at room temperature. After washing each well three times with water, the fluorescence of each well was analyzed using a Spectramax M5 plate reader (Molecular Devices, Sunnyvale, CA) with excitation 630 nm, emission 697 nm, and cutoff 695 nm, background corrected by subtracting the mean signal from empty wells and normalized to the mean value of untreated wells.
siRNA
Request a detailed protocolAll Dharmacon ON-TARGET plus siRNA pools were purchased from Thermo Fisher Scientific. Sequences are included in Supplementary file 1. Cells were plated at 12,500–25,000 per ml in 96-well plates (0.2 ml), 8-chamber slides (0.5 ml), 6-well plates (3 ml), or 60 mm plates (6 ml) in antibiotic-free medium and transfected the following day with each siRNA (12.5 nM final concentration) with the Dharmafect I transfection reagent (2 µl/200 μl for PC9 and 4 µl/200 μl for HCC827 and HCC2935) according to the manufacturer's standard protocol. Media were changed 6–24 hr after transfection, and isolation of mRNA and protein lysates and immunofluorescence analysis was carried out 48–96 hr after transfection (and treatment with erlotinib).
Lentiviral-inducible expression
Request a detailed protocolLentiviral-inducible expression constructs containing SOX2 (wild type and epitope tagged) or WNT pathway mutants under the control of a doxycycline-inducible promoter or control vector containing the GUS gene were constructed by subcloning each ORF into the pInducer 20 lentiviral vector (Meerbrey et al., 2011) using the Gateway cloning system (Thermo Fisher Scientific). This vector also contains a neomycin resistance cassette. Lentiviral expression constructs were cotransfected into 293T cells with the HIV-1 packaging construct pCMVdeltaR8.91 and the VSV-G envelope construct pMD.G using TransIT-LT-1 transfection reagent (Thermo Fisher Scientific). Viral supernatants were collected at 48 and 72 hr post transfection in DMEM with 30% FCS, filtered through 0.45 μm syringe filters to remove cell debris and stored in aliquots at −80°C. Transductions were carried out in the presence of 8 μg/ml Polybrene (Sigma) by spinoculation at 1200×g and at 32°C for 60 min in a Sorvall Legend RT table-top centrifuge. Viral supernatant was exchanged for fresh media 24 hr after spinoculation. To generate stable cDNA expressing cell lines, cells were selected in G418. Tetracycline-reduced FBS (Clontech, Mountain View, CA) was substituted for all media for cells transduced with the pInducer 20 vectors. To induce expression of SOX2 or WNT pathway constructs, doxycycline (Sigma) at 0.1 µg/ml was added to all cultures at least 3 hr prior to addition of erlotinib.
Luciferase reporter assay
Request a detailed protocolStable cell lines expressing the TOP-FLASH reporter were generated by transducing cells with 7TFP recombinant lentiviruses (Fuerer and Nusse, 2010), and luciferase assay was performed as previously described (Singh et al., 2012). Briefly, Rediject D-Luciferin Ultra (Perkin Elmer) was added in 0.2 ml fresh media (1–200 dilution) to each well of cells in a 96-well plate and incubated for 15 min at 37°C. Luciferase activity was imaged with the IVIS Lumina II In Vivo Imaging System (Perkin Elmer). The radiance of each well was determined using Living Image 4.2 software (Perkin Elmer), background corrected by subtracting the mean signal from empty wells and normalized both to the relative cell number of each well as determined by Syto60 assay and the resulting normalized mean value of untreated wells.
Xenograft tumor studies
Request a detailed protocolAll mouse studies were carried out according to Institutional IUCAC guidelines. Mice were treated by oral gavage with a single 100 mg/kg dose of erlotinib when subcutaneous tumors had reached ∼500 mm3 in sizes (∼21–28 days). PC9 and HCC827 xenograft tumors were harvested 21 hr after erlotinib treatment.
Statistical analyses
Request a detailed protocolAll statistical analyses, including number of replicates and statistical method used, are included in the relevant figure legends.
Data availability
-
Whole transcriptome analysis of erlotinib treatment in EGFR-mutant cellsID GSE51212. Publically available at the NCBI Gene Expression Omnibus (http://www.ncbi.nlm.nih.gov/geo/).
References
-
A module of negative feedback regulators defines growth factor signalingNature Genetics 39:503–512.https://doi.org/10.1038/ng1987
-
Bacterial persistence as a phenotypic switchScience 305:1622–1625.https://doi.org/10.1126/science.1099390
-
β-Catenin-SOX2 signaling regulates the fate of developing airway epitheliumJournal of Cell Science 125:932–942.https://doi.org/10.1242/jcs.092734
-
FoxO6, a novel member of the FoxO class of transcription factors with distinct shuttling dynamicsThe Journal of Biological Chemistry 278:35959–35967.https://doi.org/10.1074/jbc.M302804200
-
Epidermal growth factor receptor mutations in non-small-cell lung cancer: implications for treatment and tumor biologyJournal of Clinical Oncology 23:3227–3234.https://doi.org/10.1200/JCO.2005.09.985
-
EGFR mutation and resistance of non-small-cell lung cancer to gefitinibThe New England Journal of Medicine 352:786–792.https://doi.org/10.1056/NEJMoa044238
-
The pINDUCER lentiviral toolkit for inducible RNA interference in vitro and in vivoProceedings of the National Academy of Sciences of USA 108:3665–3670.https://doi.org/10.1073/pnas.1019736108
-
Targeted CML therapy: controlling drug resistance, seeking cureCurrent Opinion in Genetics & Development 16:92–99.https://doi.org/10.1016/j.gde.2005.11.002
-
Lung cancers with acquired resistance to EGFR inhibitors occasionally harbor BRAF gene mutations but lack mutations in KRAS, NRAS, or MEK1Proceedings of the National Academy of Sciences of USA 109:E2127–E2133.https://doi.org/10.1073/pnas.1203530109
-
EGF receptor gene mutations are common in lung cancers from ‘never smokers’ and are associated with sensitivity of tumors to gefitinib and erlotinibProceedings of the National Academy of Sciences of USA 101:13306–13311.https://doi.org/10.1073/pnas.0405220101
-
Recovery of phospho-ERK activity allows melanoma cells to escape from BRAF inhibitor therapyBritish Journal of Cancer 102:1724–1730.https://doi.org/10.1038/sj.bjc.6605714
-
(V600E)BRAF is associated with disabled feedback inhibition of RAF-MEK signaling and elevated transcriptional output of the pathwayProceedings of the National Academy of Sciences of USA 106:4519–4524.https://doi.org/10.1073/pnas.0900780106
-
Whole transcriptome analysis of erlotinib treatment in EGFR-mutant cellsGene Expression Omnibus, Accession GSE51212.
-
Genotypic and histological evolution of lung cancers acquiring resistance to EGFR inhibitorsScience Translational Medicine 3:ra26.https://doi.org/10.1126/scitranslmed.3002003
-
Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profilesProceedings of the National Academy of Sciences of USA 102:15545–15550.https://doi.org/10.1073/pnas.0506580102
-
Integrative Genomics Viewer (IGV): high-performance genomics data visualization and explorationBriefings in Bioinformatics 14:178–192.https://doi.org/10.1093/bib/bbs017
-
FoxO6 transcriptional activity is regulated by Thr26 and Ser184, independent of nucleo-cytoplasmic shuttlingThe Biochemical Journal 391:623–629.https://doi.org/10.1042/BJ20050525
-
TRANSFAC: an integrated system for gene expression regulationNucleic Acids Research 28:316–319.https://doi.org/10.1093/nar/28.1.316
-
FOXO1 is an essential regulator of pluripotency in human embryonic stem cellsNature Cell Biology 13:1092–1099.https://doi.org/10.1038/ncb2293
-
Establishment and characterization of primary lung cancer cell lines from Chinese populationActa Pharmacologica Sinica 32:385–392.https://doi.org/10.1038/aps.2010.214
Article and author information
Author details
Funding
Howard Hughes Medical Institute
- Daniel A Haber
National Institutes of Health (RO1 CA207186)
- Daniel A Haber
National Institute of Dental and Craniofacial Research (NIDCR) (K08DE020139)
- S Michael Rothenberg
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
The authors wish to thank Doug Robinson (JMP Life Sciences) for help with array analysis, Charles Vanderburg (MGH Department of Neurology) for technical help and Matthew J Niederst for providing patient-derived cell lines. This work was supported by NIH RO1 CA207186 (DAH), HHMI (DAH, MNR), National Institute of Health/National Institute of Dental & Craniofacial Research K08DE020139 (SMR) and the Burroughs Wellcome Fund (MNR).
Ethics
Animal experimentation: All animal studies were conducted through Institutional Animal Care and Use Committee (IUCAC)-approved animal protocol 2010N000006 from the Massachusetts General Hospital. Mice were maintained in laminar flow units in aseptic condition and the care and treatment of all mice was in in accordance with institutional guidelines.
Copyright
© 2015, Rothenberg et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 4,184
- views
-
- 823
- downloads
-
- 32
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Cancer Biology
- Evolutionary Biology
In asexual populations that don’t undergo recombination, such as cancer, deleterious mutations are expected to accrue readily due to genome-wide linkage between mutations. Despite this mutational load of often thousands of deleterious mutations, many tumors thrive. How tumors survive the damaging consequences of this mutational load is not well understood. Here, we investigate the functional consequences of mutational load in 10,295 human tumors by quantifying their phenotypic response through changes in gene expression. Using a generalized linear mixed model (GLMM), we find that high mutational load tumors up-regulate proteostasis machinery related to the mitigation and prevention of protein misfolding. We replicate these expression responses in cancer cell lines and show that the viability in high mutational load cancer cells is strongly dependent on complexes that degrade and refold proteins. This indicates that the upregulation of proteostasis machinery is causally important for high mutational burden tumors and uncovers new therapeutic vulnerabilities.
-
- Cancer Biology
- Cell Biology
Understanding the cell cycle at the single-cell level is crucial for cellular biology and cancer research. While current methods using fluorescent markers have improved the study of adherent cells, non-adherent cells remain challenging. In this study, we addressed this gap by combining a specialized surface to enhance cell attachment, the FUCCI(CA)2 sensor, an automated image analysis pipeline, and a custom machine learning algorithm. This approach enabled precise measurement of cell cycle phase durations in non-adherent cells. This method was validated in acute myeloid leukemia cell lines NB4 and Kasumi-1, which have unique cell cycle characteristics, and we tested the impact of cell cycle-modulating drugs on NB4 cells. Our cell cycle analysis system, which is also compatible with adherent cells, is fully automated and freely available, providing detailed insights from hundreds of cells under various conditions. This report presents a valuable tool for advancing cancer research and drug development by enabling comprehensive, automated cell cycle analysis in both adherent and non-adherent cells.