MAPK activity dynamics regulate non-cell autonomous effects of oncogene expression
Abstract
A large fraction of human cancers contain genetic alterations within the Mitogen Activated Protein Kinase (MAPK) signaling network that promote unpredictable phenotypes. Previous studies have shown that the temporal patterns of MAPK activity (i.e. signaling dynamics) differentially regulate cell behavior. However, the role of signaling dynamics in mediating the effects of cancer driving mutations has not been systematically explored. Here, we show that oncogene expression leads to either pulsatile or sustained ERK activity that correlate with opposing cellular behaviors (i.e. proliferation vs. cell cycle arrest, respectively). Moreover, sustained–but not pulsatile–ERK activity triggers ERK activity waves in unperturbed neighboring cells that depend on the membrane metalloprotease ADAM17 and EGFR activity. Interestingly, the ADAM17-EGFR signaling axis coordinates neighboring cell migration toward oncogenic cells and is required for oncogenic cell extrusion. Overall, our data suggests that the temporal patterns of MAPK activity differentially regulate cell autonomous and non-cell autonomous effects of oncogene expression.
eLife digest
In animals, the MAPK pathway is a network of genes that helps a cell to detect and then respond to an external signal by switching on or off a specific genetic program. In particular, cells use this pathway to communicate with each other. In an individual cell, the MAPK pathway shows fluctuations in activity over time.
Mutations in the genes belonging to the MAPK pathway are often one of the first events that lead to the emergence of cancers. However, different mutations in the genes of the pathway can have diverse effects on a cell’s behavior: some mutations cause the cell to divide while others make it migrate. Recent research has suggested that these effects may be caused by changes in the pattern of MAPK signaling activity over time.
Here, Aikin et al. used fluorescent markers to document how different MAPK mutations influence the behavior of a human breast cell and its healthy neighbors. The experiments showed that cells with different MAPK mutations behaved in one of two ways: the signaling quickly pulsed between high and low levels of activity, or it remained at a sustained high level. In turn, these two signaling patterns altered cell behavior in different ways. Pulsed signaling led to more cell division, while sustained signaling stopped division and increased migration.
Aikin et al. then examined the effect of the MAPK mutations on neighboring healthy cells. Sustained signaling from the cancerous cell caused a wave of signaling activity in the surrounding cells. This led the healthy cells to divide and migrate toward the cancerous cell, pushing it out of the tissue layer. It is not clear if these changes protect against or promote cancer progression in living tissue. However, these results explain why specific cancer mutations cause different behaviors in cells.
Introduction
The Receptor-Tyrosine Kinase (RTK)/RAS/ERK signaling axis (Figure 1A) is mutated in most human cancers (Sanchez-Vega et al., 2018). In normal conditions, the ERK pathway promotes proliferation, differentiation, survival and cell migration (Johnson and Lapadat, 2002). During oncogenesis, mutations or amplification of ERK pathway components can also promote oncogene-induced senescence (Hahn and Weinberg, 2002) (OIS) or oncogenic cell extrusion from epithelial monolayers in the so-called Epithelial Defense Against Cancer response (EDAC) (Hogan et al., 2009; Kajita et al., 2010). The mechanisms underlying dose-dependent effects of ERK signaling have been intensely studied using bulk cell population assays. However, the advent of single-cell analysis has shown that single cells often behave qualitatively different than bulk populations. In fact, in vivo and in vitro studies have now shown that pulsatile or sustained ERK activity have different effects on cell behavior (Albeck et al., 2013; Aoki et al., 2013; de la Cova et al., 2017; Johnson and Toettcher, 2019; Santos et al., 2007; Bugaj et al., 2018; Aoki et al., 2017). Whether different oncogenic perturbations also have different functional outcomes depending on downstream signaling dynamics remains unknown. To address this question, an isogenic single-cell approach with temporal control of oncogene expression is needed.
Recent in vivo studies revealed that oncogene expression can trigger tissue level responses involving normal neighboring cells (Brown et al., 2017; Ellis et al., 2019; Clavería et al., 2013; Sancho et al., 2013). In specific cases, mosaic oncogene expression leads to either basal extrusion or apical extrusion (Hogan et al., 2009; Kajita et al., 2010); however, the signaling mechanism responsible for recognition between normal and diseased cells is poorly understood (Kajita and Fujita, 2015; Clavería and Torres, 2016; Maruyama and Fujita, 2017). Coincidentally, propagating ERK signaling waves requiring the sheddase ADAM17 have been observed in mouse epidermis and intestinal organoids, but the physiological role of these signaling events remains unclear (Hiratsuka et al., 2015; Muta et al., 2018). Observation of interactions between oncogenic and neighboring epithelium with live-cell biosensors could provide insights into the collective signaling preceding oncogenic extrusion. In fact, a recent study using live imaging of calcium biosensors during EDAC of HRASG12V cells showed a calcium signaling wave which propagated through neighboring epithelium to coordinate actin rearrangements and polarized movements during apical extrusion (Takeuchi et al., 2020). The mechanistic basis underlying EDAC calcium waves remains unknown.
Here, we combine live cell imaging of MAPK activity biosensors with inducible expression of oncogenes to study the effects of oncogene expression on signaling dynamics and how altered MAPK dynamics impact both cell autonomous and non-cell autonomous behaviors in epithelial tissues. Our data shows that pulsatile or sustained ERK signaling resulting from oncogenic perturbations triggers different dynamics-dependent cell behaviors including oncogene-induced paracrine signaling via the ADAM17-AREG-EGFR signaling axis. The resulting signaling gradients are required to coordinate neighboring cell migration and active oncogenic cell extrusion (EDAC). Our study highlights the role of MAPK signaling dynamics in coordinating individual and collective cell behaviors.
Results
To study the effects of oncogene expression on the temporal patterns of MAPK signaling we generated a reporter cell line derived from the chromosomally-normal human breast epithelial line, MCF10A, expressing the ERK Kinase Translocation Reporter (Regot et al., 2014) (ERK KTR) and a fluorescently tagged ERK kinase (ERK-mRuby2). This combination of biosensors allowed independent measurement of ERK activity and ERK localization in live single cells at high temporal resolution. Then, we introduced 12 different doxycycline-inducible oncogenic perturbations via lentiviral infection and measured ERK signaling dynamics during overexpression (Figure 1B). Our results revealed two qualitatively different responses to oncogene induction: (i) increased frequency of ERK activity pulses with no change in ERK kinase localization (i.e. EGFR, B-RafWT), and (ii) sustained ERK activity with subsequent nuclear translocation of ERK kinase (i.e. B-RafV600E, MEK2DD) (Figure 1, and Figure 1—Video 1). We refer to these distinct dynamics as pulsatile or sustained ERK, respectively. Of note, MEK1/2WT expression is capable of exporting ERK into the cytoplasm without changing kinase activity (Adachi et al., 2000; Figure 1—figure supplement 1) and ERK nuclear accumulation occurs only when activity is sustained, suggesting that ERK activity and ERK localization are not always correlated. Interestingly, expression of B-RafWT or B-RafV600E elicit qualitatively different downstream dynamics even though they differ in a single amino acid and show similar expression levels by immunoblotting (Figure 1—figure supplement 2). Given that the B-RafV600E is insensitive to negative feedback regulation by ERK (Yao et al., 2015), this result suggests that ERK inhibition to B-RafWT is mechanistically involved in the characteristic pulsatile dynamics.
Next, we assessed how ERK dynamics affect cell behaviors by measuring cell migration and proliferation. While pulsatile ERK dynamics (i.e. EGFR or B-RafWT) consistently correlated with increased cell cycle progression, sustained ERK activity (i.e. B-RafV600E or MEK2DD) caused cell cycle arrest and increased migration (Figure 1D). Importantly, observed differences in cell behavior correlated with dynamics independently of the point in the cascade that perturbations were introduced (EGFR, Raf or MEK), suggesting that ERK is responsible for differences in cell behaviors rather than alternate downstream pathways. Moreover, expression of B-RafWT or B-RafV600E, which activate the cascade at the same point, caused different ERK activity dynamics (i.e. pulsatile or sustained respectively) and triggered opposing cellular behaviors (Figure 1D). Taken together, these data suggest that ERK activity dynamics can either promote or inhibit proliferation cell autonomously.
The sudden increase in migration and the loss of cell-cell contacts observed in cases where ERK activity is sustained (B-RafV600E and MEK2DD, Figure 1—Video 1) are reminiscent of phenotypes described for cells undergoing Epithelial-to-Mesenchymal Transition, or EMT (Hao et al., 2019). We sought out to address the role of EMT in oncogene-dependent cell behaviors by immunofluorescent staining of an epithelial marker E-Cadherin (E-Cad) and the mesenchymal marker N-Cadherin (N-Cad). While cell migration was clearly increased at 24 hr post-oncogene expression, cells retain E-Cad expression with no clear increase in N-Cad, as was observed in TGFβ-induced EMT (Figure 1—figure supplement 3). These results indicate that at the time points studied here, altered cell behaviors are either distinct from or precede those resulting from EMT.
To examine the non-cell autonomous effects of oncogene expression in epithelial monolayers, we cocultured ‘inducible’ cells (expressing doxycycline-inducible oncogenes, a constitutively expressed H2B-mClover, and the ERK biosensors) with ‘neighboring’ reporter cells (expressing ERK biosensors without inducible oncogenes) and monitored signaling dynamics upon induction (Figure 2A). Interestingly, expression of B-RafV600E, but not B-RafWT, resulted in waves of ERK activation of neighboring cells (Figure 2B–C and Figure 2—Video 1). This comparison suggests that oncogenic perturbations that elicit sustained ERK activity propagate ERK activity pulses to neighboring cells. In agreement, other oncogenes that triggered sustained, but not pulsatile, ERK activity also promoted ERK activity waves in the neighboring cells (Figure 2—figure supplement 1). By using KTRs for p38 and JNK, we observed that neighboring epithelia did not activate other the MAPK pathways (Figure 2—figure supplement 1). Notably, spontaneous cell death events were also followed by similar ERK signaling waves (Figure 2—Video 2), indicating that oncogene expression and cell death may be similarly perceived by neighboring cells.
We then addressed the mechanistic basis of oncogene-dependent paracrine signaling. Previous studies demonstrated that ERK waves in epithelial monolayers depend on the membrane-tethered sheddase ADAM17, which releases membrane-bound growth factors that activate EGFR signaling in adjacent cells (Aoki et al., 2013; Aoki et al., 2017; Hiratsuka et al., 2015). Thus, we hypothesized that oncogenic cell ADAM17 may be decoding ERK signaling dynamics to trigger growth factor release. To test this hypothesis, we generated an ADAM17 knockout (ADAM17KO) cell line (Figure 3A) and used it as either ‘inducible’ or ‘neighboring’ cells in our coculture assay. Live imaging of WT and ADAM17KO cocultures indicated that ADAM17 is necessary in inducible, but not neighboring cells, to trigger ERK waves in the monolayer (Figure 3B–C and Figure 2—Video 1). Therefore, ADAM17 decodes ERK activity dynamics in inducible cells to transmit ERK signaling to neighboring cells. Previous work has shown that ADAM17 is weakly phosphorylated compared to other ERK substrates (Díaz-Rodríguez et al., 2002), thus the phosphorylation-dephosphorylation kinetics of ADAM17 and the temporal patterns of ERK activity may explain dynamics-specific ADAM17 activation.
ADAM17-released growth factors include HB-EGF, TGF-α, Epiregulin, and Amphiregulin (Zunke and Rose-John, 2017; Rios-Doria et al., 2015). In order to identify the factors mediating oncogene-induced paracrine signaling we used Tandem-Mass-Tag Mass Spectrometry of supernatant proteins following induction of sustained ERK activity in WT and ADAM17KO cells. A variety of known and unknown ADAM17 substrates were present in the induced cell supernatants, including immune surveillance (HLA-A/B/C), Delta-Notch (JAG1), and Wnt (SFRP) signaling proteins (Figure 3D and Supplementary file 1). Of note, the EGFR ligand Amphiregulin (AREG) was the most upregulated, ADAM17-dependent protein in the supernatant, suggesting that AREG released from inducible cells could act as an oncogene-dependent paracrine signaling molecule. Accordingly, cocultures pre-incubated with AREG function-blocking antibodies or EGFR inhibitors prevented neighboring cell ERK activation without affecting ERK signaling in inducible cells (Figure 3E–F). These results indicate that oncogene-dependent ERK waves are mediated by ADAM17 (in inducible cells), AREG release, and EGFR signaling (Figure 3G).
Given that cells surrounding B-RafV600E expressing cells showed pulsatile ERK activity (Figure 2C), we hypothesized that oncogene expression may promote cell proliferation in a non-cell autonomous manner. Accordingly, sustained ERK signaling in inducible cells increased proliferation of neighboring cells up to 10-fold (Figure 4A–C). Together, these data indicate that, depending on ERK dynamics, oncogenic cells can have either cell autonomous or non-cell autonomous contributions to tissue growth.
In addition to proliferation, ERK waves have been shown to orient collective cell migration during wound healing (Aoki et al., 2017). In cocultures, sustained ERK activity in B-RafV600E-inducible cells correlated with neighboring cell migration towards inducible cells in an ADAM17 and EGFR-dependent manner (Figure 5A). We hypothesized that coordinated migration of neighboring cells could physically contribute to oncogenic cell extrusion (Hogan et al., 2009; Leung and Brugge, 2012; Slattum et al., 2014). To address this hypothesis, we used confocal Z stacks to quantify extrusion of oncogene-expressing cells from monolayers (Figure 5B and Figure 5—Video 1). Interestingly, while pulsatile ERK activity (i.e. EGFR and B-Raf) was not sufficient to extrude cells, sustained ERK activity (i.e. B-RafV600E and MEK2DD) led to efficient epithelial cell extrusion apically (Figure 5C). KRASG12V induction did not result in apical extrusion to the extent observed for HRASG12V (Hogan et al., 2009, Figure 5—figure supplement 1). However, since sustained ERK activation in KRASG12V occurs later than B-RafV600E (Figure 1—figure supplement 1) further apical extrusion may also occur at a later time. Taken together, our data suggests that apical extrusion occurs when oncogenic perturbations trigger sustained ERK activity.
In mammalian epithelia, apical extrusion eliminates apoptotic cells or crowded cells to maintain homeostasis (Rosenblatt et al., 2001; Eisenhoffer et al., 2012). Similar, but mechanistically different, apical extrusion has been observed for some oncogenic cells during EDAC (Kajita and Fujita, 2015). We wanted to know whether the extrusion of inducible cells with sustained ERK activity resembled EDAC. To differentiate between pseudostratified or de-laminated (Grieve and Rabouille, 2014) epithelium and extrusion we analyzed confocal images of E-cad membrane staining in induced cocultures. Cells with sustained ERK activity were fully extruded, sitting above WT cells in the plane of the monolayer (Figure 5—figure supplement 2). These images also demonstrate maintenance of E-Cad at the junctions between WT cells below extruded cells. Quantification of fully-extruded cells at several timepoints showed the majority of oncogenic cells being extruded from 4 to 8 hr after induction, but continuing until 24 hr, when 91% of inducible cells are fully extruded.
Both the oncogenic and apoptotic extrusion models involve cytoskeletal rearrangements at the site of extrusion (Rosenblatt et al., 2001; Kajita and Fujita, 2015). To observe live actin dynamics in cocultures, we made cell lines stably expressing Utrophin-261-EGFP (Belin et al., 2014). Using this tool, we observed transient accumulation of actin at the basal interface of B-RafV600E expressing and neighboring cells that first closed off the basal attachments of inducible cells before they were pushed apically out of the monolayer (Figure 5—figure supplement 2 and Figure 5—Video 2). These polarized, actin-containing basal protrusions were dependent on EGFR activity as they could be inhibited by EGFR inhibitor. Apoptotic extrusion relies on Sphingosine-1-phosphate (S1P) signaling through intrinsic S1P production and juxtracrine activation of the GPCR S1PR2 (Gu et al., 2011), yet inhibition of S1P production had only moderate effects on apical extrusion of MEK2DD cells (Figure 5—figure supplement 2). Together, our results suggest that the apical extrusion of oncogenic cells observed in our experiments are similar to the EDAC mechanism previously described for HRASG12V, V-Src, and other cells (Hogan et al., 2009; Kajita et al., 2010; Kajita and Fujita, 2015).
The requirement for paracrine signals in collective migration led to the question of whether paracrine signals were also required for extrusion. To test the role of ADAM17-mediated AREG-EGFR paracrine signals in promoting extrusion, we performed extrusion assays using ADAM17KO cells or in the presence of EGFR inhibitor or AREG function-blocking antibodies. Extrusion of inducible cells was abolished in these conditions (Figure 5D–G), suggesting that ERK signaling waves are required for extrusion. Of note, since ADAM17KO and EGFR inhibition affect ERK activation of neighboring cells without altering ERK dynamics in inducible cells, we hypothesized that that activation of inducible cells alone is not sufficient for extrusion, but that neighboring cell ERK activation may be required. To address this question, ERK-independent ADAM17 activation is needed.
Previous studies have shown that the stress MAPK p38 phosphorylates and activates ADAM17 (Xu and Derynck, 2010). Thus, we used our doxycycline-inducible system to drive the expression of MKK3DD, a constitutively-active MAP2K specific for p38 (Figure 6—figure supplement 1), to activate ADAM17 independently of ERK. As expected, we found that p38 activation leads to ERK signaling waves (Figure 6A–B), proliferation (Xu and Derynck, 2010; Figure 6—figure supplement 2), directed migration (Figure 6C) and extrusion in an ADAM17 and EGFR dependent manner (Figure 6D and Figure 6—Video 1). However, B-RafV600E oncogenic signaling, extrusion and proliferation were unaffected by p38 inhibition (Figure 6—figure supplement 3), suggesting that sustained ERK or p38 activity are each capable of activating ADAM17 paracrine signaling. Using this ERK-independent MKK3DD extrusion system, we found that MEK inhibition decreased directed migration and prevented extrusion, confirming that extrusion requires ERK activity in the neighboring cells (Figure 6C–D). Thus, ERK activity is required for extrusion in both oncogenic and neighboring cells with qualitatively different temporal dynamics. Notably, inhibition of S1P signaling in MKK3DD cocultures also reduced extrusion efficiency despite having unaffected signaling (Figure 6—figure supplement 4). This result suggests that cell-autonomous ERK or p38 activation in extruded cells may underlie the fundamental differences reported between extrusion of oncogenic and stressed cells.
Finally, we asked whether the spatiotemporal properties of paracrine ERK signaling waves are important to coordinate extrusion. We first tested the efficiency of extrusion with altered proportions of B-RafV600E cells in the coculture, as higher proportions will have de-centralized and overlapping signaling events. The proportion of inducible cells was inversely correlated with extrusion efficiency (Figure 7A). Moreover, exogenous addition of AREG, which triggers widespread ERK activation preventing any spatially defined waves, eliminated directed migration of neighboring cells and extrusion (Figure 7B–C). The observation that in cocultures, polarized actin enrichment in neighboring cell basal protrusions is absent with EGFR inhibition, also indicates that growth factor signaling provides directional information during extrusion (Figure 5—figure supplement 2). Together this data suggests that locally generated paracrine signaling coordinates directed migration of neighboring epithelia to promote extrusion of oncogenic cells (Figure 8).
Discussion
A wide variety of ERK pathway alterations occur across human tumors, often resulting in different cancer phenotypes (Sanchez-Vega et al., 2018; Bailey et al., 2018; Hoadley et al., 2018). To better understand the signaling effects of oncogenic mutations, we used live cell imaging of signaling biosensors upon induction of oncogenes and found that different ERK pathway oncogenes trigger two distinct temporal patterns: pulsatile or sustained ERK activity. While our approach is admittedly different than acquisition of point mutations in vivo, ERK dynamics resulting from oncogene overexpression robustly correlated with the same cellular phenotypes: (i) pulsatile ERK activity correlates with increased proliferation and, (ii) sustained ERK activity leads to cell cycle arrest similar to OIS (Hahn and Weinberg, 2002; Courtois-Cox et al., 2006; Kuilman et al., 2008). Moreover, we showed that sustained ERK activity in oncogenic cells triggers ERK signaling waves through unperturbed neighboring cells. These signaling waves depend on the ADAM17-EGFR paracrine signaling axis and lead to different non cell-autonomous behaviors such as (i) proliferation, (ii) directed migration of neighboring cells toward oncogenic cells, and (iii) oncogenic cell extrusion (Figure 8).
Our data indicates that cancer mutations can have non-cell autonomous contributions to tissue growth (Figure 4). Interestingly, studies in mouse epidermis have shown that mosaic oncogene expression promotes proliferation of wild type surrounding cells, which is required to expel mutant outgrowths from the tissue (Brown et al., 2017). While the role of ADAM17 in this phenomenon and during early tumorigenesis in vivo is yet unknown, it is tempting to speculate that polypous outgrowths may occur in the presence of non-proliferative oncogenic cells that release growth factors via ADAM17. The mechanisms that mediate tissue expelling in vivo remain unknown.
Previous work in described the process of oncogenic cell extrusion as part of the so called Epithelial Defense Against Cancer (i.e. EDAC) (Hogan et al., 2009). However, the signals involved in recognition of oncogenic cells, and why only specific oncogenes trigger oncogenic cell extrusion was unclear (Kajita and Fujita, 2015; Clavería and Torres, 2016; Maruyama and Fujita, 2017). Our data suggests that perturbations that elicit sustained ERK activity (eg. B-RafV600E, MEK2DD), activate ADAM17, which in turn releases EGFR ligands (Figure 3). This paracrine signal is critical for oncogenic extrusion (Figure 5). However, we acknowledge that different cellular states such as apoptosis or overcrowding lead to extrusion by different mechanisms. Of note, our data showed that ERK activation drives extrusion to a higher extent than p38 activation (Figure 6), which may result from a difference in overall cell autonomous migration in these two cases (Figure 6—figure supplement 1). Moreover, sphingosine kinase inhibition caused greater defects in extrusion of p38-active cells than ERK-active cells (Figure 5—figure supplement 2 and Figure 6—figure supplement 4). This finding agrees with work showing that EDAC of transformed HRASG12V cells is less dependent on sphingosine-1-phosphate production than extrusion of crowded or apoptotic cells (Yamamoto et al., 2016), where stress signaling may be involved.
We and others have identified AREG as one of the key EGFR ligands in mammary epithelial cells (Sternlicht et al., 2005; Figure 3); however, different ligands may be required in other tissues. These ligands, released by ADAM17, coordinate the migration of neighboring cells by mechanisms that remain unclear. Cultured monolayers are fundamentally different than in vivo tissues; however, the chemo-attractive properties of growth factors for directed migration have been modeled and studied in cell culture (Devreotes et al., 2017; Tranquillo et al., 1988). We propose that local signaling gradients are created by oncogenic cells to coordinate directed migration of neighbors. To support this idea, we show that addition of exogenous AREG or increased fractions of oncogenic cells both prevent directed migration of neighboring cells (Figure 7), and that during extrusion, polarized actin-containing basal protrusions require growth factor signaling (Figure 5—figure supplement 2). Localized ERK signaling gradients have also been observed during morphogenesis of Drosophila, avian, and mammalian embryos (Yang et al., 2002; Ogura et al., 2018; Corson et al., 2003), and in preserving homeostasis of mammalian epidermis and intestinal organoids (Hiratsuka et al., 2015; Muta et al., 2018; Liang et al., 2017). Thus, in addition to roles in oncogenesis, the ADAM17-EGFR paracrine signaling axis may direct collective behaviors during development.
Overall, our results highlight the importance of quantitative live-cell approaches to understand the effects of genetic perturbations and cell-cell communication in tissues. We propose a critical role for ERK signaling dynamics and the ADAM17-EGFR signaling axis in coordinating cell behaviors at the tissue level.
Materials and methods
Cell lines & reagents
Request a detailed protocolMCF10A human mammary epithelial cells (ATCC) were grown at 37° and 5% CO2 in DMEM/F12 (Gibco) with 5% horse serum (HS) (Sigma), 10 µg/ml Insulin (Sigma), 20 ng/ml EGF (Peprotech), 1x Penicillin-Streptomycin (P/S) (Gibco), 0.5 mg/ml Hydrocortisone (Sigma), 100 ng/ml Cholera Toxin (Sigma). Cells were passaged every 3 days with 0.25% Trypsin-EDTA (Gibco), are mycoplasma free, and were verified by STR-profiling (ATCC).
Cell lines were generated with lentivirus produced in HEK293-FTs (Thermo) with third-generation packaging plasmids and Lipofectamine 2000 (Thermo). Viral supernatants were collected 48 hr after transfection and incubated in MCF10As with polybrene (10 µg/ml, EMD Millipore). To create dual-sensor cells, MCF10As were infected with a lentiviral H2B-iRFP vector (Addgene) and sorted. We used gateway cloning (Campeau et al., 2009) to introduce ERK-KTR-mCer3 and ERK1-mRuby2 into PGK pLenti DEST vectors (Addgene), infected and selected the H2B-iRFP MCF10As (Blasticidin 3 µg/ml and Hygromycin 10 µg/ml Corning). We isolated moderately expressing clones using cloning cylinders (EMD Milipore). For inducible cells, a gateway-ready reverse TET trans-activator (rtTA) plasmid was created by adding the rtTA with a 2A peptide to the Puromycin resistance gene in a CMV Puro DEST plasmid (Addgene) by gibson cloning (Gibson et al., 2009). Human coding sequences were acquired from either Addgene or the Thermo Ultimate ORF Collection, sequence verified, and introduced in the rtTA CMV Puro DEST plasmid by gateway cloning (Campeau et al., 2009). These plasmids were used for lentivirus, and infected cells were selected with Puromycin (1 µg/ml, Sigma). Utrophin-261-EGFP cell lines were made by cloning the coding region from pEGFP-C1 Utr261-EGFP (Addgene) into a pENTR backbone by Gibson cloning, and then introduced into the pLenti PGK Puro DEST plasmid by gateway cloning. These plasmids were used to generate lentivirus, and infected cells were selected with Puromycin.
For inhibitor experiments, small molecules or antibodies and doxycycline were dissolved to a 10X working concentration in imaging media before addition. Final DMSO concentration did not exceed 0.15%. Inhibitors used include the MEK inhibitor PD-0325901, the MMP/ADAM inhibitor Batimastat, the EGFR inhibitor Gefitinib, the p38 inhibitor BIRB-796, the Sphingosine Kinase inhibitor SKII, and the S1PR2, inhibitor JTE-013 all from Selleck Chemicals. The p38 inhibitor SB-203580 was obtained from Sigma. Amphiregulin was ordered from Peprotech. Amphiregulin function-blocking antibody is from R and D systems.
The ADAM17KO cell lines were created using the CRISPR V2 Neo system (a gift from Dr. Andrew Holland) and gRNA oligos targeting R241 of exon 6. Dual sensor cells were infected with lentivirus carrying this plasmid, selected with Neomycin (500 µg/ml, Sigma) and clonally expanded before western blot validation (Figure 2B).
Live imaging
Request a detailed protocolCells were plated at 3*105 cells/well in fibronectin-treated (EMD Millipore) 96-well glass-bottom plates (Thermo Scientific) 48 hr before imaging. The following day, monolayers were serum-starved with 0.5% HS, phenol-red-free DMEM/F12 containing P/S with 1 mM Na Pyruvate and 10 mM HEPES. For signaling experiments in Figure 1 and Figure 1—figure supplement 1, media was switched to 0% HS several hours before imaging to limit basal signaling. Monolayers were imaged using a Metamorph-controlled Nikon Eclipse Ti-E epifluorescence microscope with a 20x air objective and a Hamamatsu sCMOS camera. The multi LED light source SpectraX (Lumencor) and the multiband dichroic mirrors DAPI/FITC/Cy3/Cy5 and CFP/YPF/mCherry (Chroma) where used for illumination and imaging without any spectral overlap. For extrusion and live-actin experiments, a Metamorph-controlled Nikon Eclipse Ti-E spinning-disc confocal (Yokogawa W1) with a 20x or 40X objective, Prime 95-B sCMOS camera (Photometrics) and a Multiline laser launch (Cairn Research) was used to capture H2B-iRFP and H2B-mClover or Utrophin-261-EGFP images every 1 µm of a 25–30 µm range through monolayers. Temperature (37°C), humidity and CO2 (5%) were maintained throughout all imaging using OKO Labs control units. Sample sizes were selected by attempting to capture at least 100 cells from each population, with several hundred cells preferred. Key conditions from imaging experiments were performed at least twice, with one replicate presented in figures.
Image analysis and quantification
Request a detailed protocolPrimary time-lapse images were subjected to flat-fielding and registration (custom software Aikin et al., 2020) before object segmentation and measurements in Cell Profiler. Nuclear positions were used to track individual cells through time-series (custom software Aikin et al., 2020) and intensity ratios were calculated as previously described (Regot et al., 2014). Minimal cleaning of traces excluded cells where tracks switched between two objects, where the KTR ratios were affected by segmentation errors, or where traces represent less than two thirds of the entire time-course. In conditions where cells move rapidly, such as B-RafV600E and MEK2DD, and traces are shorter due tracking errors, track-length restraints were relaxed to include more cells for analysis. Single-cell traces were chosen by random plotting of distinct cells and selection of those that were tracked throughout the whole experiment. Peak counting was performed with software based on findPeaks (O'Haver, 2014; Mathworks.com) and modified to detect peaks based on the rate of change between gaussian-fitted minima and maxima from single-cell traces.
For directed migration, positions were selected where distinct groups of inducible cells were present in the center of the field of view. Migration was quantified by positional changes over 20 min intervals for specified time windows, from all WT neighboring cells within a 200 µM X 200 µM area centered on the group of inducible cells. The migration angles of neighboring cells are plotted as radial histograms where 0° indicates migration directly towards, and 180° directly away from the center of isolated inducible cell groups. Migration datasets contain many sampled angles from large populations of cells. To overcome issues with high power, we applied subsampling techniques using 1000 iterations of 1000 randomly-selected migration angles each, and presented the median Two-Sample Kolmogorov-Smirnov (KS) Test P-values from these iterations (‘ns’, not significant, *p<0.05, **p<0.01, ***p<0.001).
For extrusion experiments, histograms of mClover and mRuby pixel intensities across each z-stack were fit to gaussian curves using Matlab. The difference in gaussian fitted maxima of inducible cells and neighboring cells for each observation are plotted. Extrusion experiment sample size represents all non-overlapping positions from 2 to 3 independent wells excluding outliers resulting from imaging artifacts. Two-sample T-test significance values compare indicated conditions (‘ns’, not significant, *p<0.05, **p<0.01, ***p<0.001).
For live-actin imaging experiments, ERK-KTR, H2B-iRFP were infected with the TRE3G::BRAFV600E. These cells were plated in 1% coculture with neighboring ERK-KTR, H2B-iRFP cells containing the Utrophin-261-EGFP construct. Oncogenic cells were identified by lack of green fluorescence and confirmed by images of the KTR, showing activation after induction with doxycycline (2 µg/ml). Actin enrichment was quantified by manually tracing the border of adjacent Utrophin-261-EGFP cells at the leading edge in contact with oncogenic cells on Fiji. All cells from a single extrusion event are plotted along with their mean.
Immunoblotting
Request a detailed protocolFor assessment of protein expression in Figure 1—figure supplement 2 , parental, TRE3G::BRAFWT and TRE3G::BRAFV600E cells were plated in 6-well plastic culture plates, and starved with 0.5% HS, DMEM/F12 containing P/S with 1 mM Na Pyruvate and 10 mM HEPES overnight before treatment with media or doxycycline (2 µg/ml) for 24 hr. Samples were lysed with RIPA buffer (CST) containing HALT protease and phosphatase inhibitors (Thermo), and reduced in Laemelli SDS buffer (BioRad) with BME (Sigma). Samples underwent electrophoresis on 4–15% gradient polyacrylamide gels (BioRad) and were immunoblotted with Rabbit anti-BRAF (CST) and mouse anti-HSC70 (Santa Cruz Biotechnology), and IRDye donkey anti-rabbit 800 and goat anti-mouse 680 secondary antibodies (LiCor) before imaging. For validation of ADAM17 CRISPR-KOs in Figure 3 , suspected clones were grown, lysed, and run on a gel as described above, before immunoblotting with Rabbit anti-ADAM17 (CST) and mouse anti-HSC70 (Santa Cruz Biotechnology) primary and IRDye donkey anti-rabbit 800 and goat anti-mouse 680 secondary antibodies (LiCor). All images were acquired on an Odyssey Infrared Scanner (LiCor).
Proteomics
Request a detailed protocolFor mass spectrometry, cells were grown to 90% confluency in T175 flasks and serum starved 24 hr (see live imaging) before switching to 15 mL growth factor/serum-free DMEM/F12 +/- Dox for 4 hr. The supernatant was collected and concentrated using 3 kDa cut-off centrifugal filters (Millipore-Sigma). Triplicate samples were quantified by the Pierce Assay (Thermo Scientific), reduced, alkylated, and trypsin digested before labeling with Tandem Mass Tag labels. Peptide fractions were analyzed by LC/MSMS using an Easy-LC 1200 HPLC system interfaced with an Orbitrap Fusion Lumos Tribrid Mass Spectrometer (Thermo Fisher Scientific). Isotopically resolved masses in precursor and fragmentation spectra were processed in Proteome Discoverer software (v2.3, Thermo Scientific). Data were searched using Mascot (2.6.2, Matrix Science) against the 2017_Refseq 83 Human database and filtered at a 1% FDR confidence threshold.
Cell proliferation assay
Request a detailed protocolMonolayers were plated and starved as described above and treated with doxycycline (Dox, 2 µg/ml) in the presence of indicated inhibitors for 24 hr. During the final 4 hr, EdU (10 µM, Thermo Fischer Scientific) was added into cultures to label S phase cells then fixed with methanol and washed before Alexa-Fluor Azide 488 click labelling (Thermo Fischer Scientific) and DAPI staining (Thermo Scientific). Monolayers were imaged by epifluorescence. Because methanol fixation eliminates fluorescence from fluorescent proteins, cocultures were imaged just before fixation and registered with DAPI and EdU images to determine positions of inducible and neighboring cells. Sample size for population EdU experiments represents all non-overlapping positions from 2 to 3 independent wells, excluding outliers resulting from imaging artifacts. Key conditions were replicated at least twice.
Immunofluorescence
Request a detailed protocolMonolayers were plated and starved as described above, and treated with media or doxycycline (Dox, 2 µg/ml) in the presence of any indicated inhibitors for 24 hr or timepoints as marked. To induce EMT, parental cells were maintained in full serum supplemented with TGFβ (5 ng/ml or 50 ng/ml, R and D Systems) through splittings over 8 days to induce EMT (Hao et al., 2019), then cells were plated and starved as described with consistent TGFβ. Cells were fixed 15 min with 4% PFA in PBS, washed with PBS before incubating 1.5 hr in blocking buffer (PBS + 0.3% Triton X-100 + 5% BSA), followed by PBS washing and incubation overnight in blocking buffer with added primary antibodies (Rabbit anti-E-Cadherin, or Rabbit anti-N-Cadherin, both CST). The following day, cells were washed in PBS before incubating 2 hr in blocking buffer with secondary antibody (Donkey anti-Rabbit IgG Alexa Fluor 405, Abcam). Cells were then washed with PBS and stored at 4°C until imaging via spinning disk confocal as described above. All incubations occurred at room temp in the dark, except the overnight primary, which was incubated at 4°C.
Appendix 1
Data availability
All data generated or analysed during this study are included in the manuscript and supporting files.
References
-
Nuclear export of MAP kinase (ERK) involves a MAP kinase kinase (MEK)-dependent active transport mechanismThe Journal of Cell Biology 148:849–856.https://doi.org/10.1083/jcb.148.5.849
-
Propagating wave of ERK activation orients collective cell migrationDevelopmental Cell 43:305–317.https://doi.org/10.1016/j.devcel.2017.10.016
-
Cell competition: mechanisms and physiological rolesAnnual Review of Cell and Developmental Biology 32:411–439.https://doi.org/10.1146/annurev-cellbio-111315-125142
-
Excitable signal transduction networks in directed cell migrationAnnual Review of Cell and Developmental Biology 33:103–125.https://doi.org/10.1146/annurev-cellbio-100616-060739
-
Extracellular cleavage of E-cadherin promotes epithelial cell extrusionJournal of Cell Science 127:3331–3346.https://doi.org/10.1242/jcs.147926
-
Epithelial cell extrusion requires the sphingosine-1-phosphate receptor 2 pathwayThe Journal of Cell Biology 193:667–676.https://doi.org/10.1083/jcb.201010075
-
Modelling the molecular circuitry of CancerNature Reviews Cancer 2:331–341.https://doi.org/10.1038/nrc795
-
TGF-β-Mediated Epithelial-Mesenchymal transition and Cancer metastasisInternational Journal of Molecular Sciences 20:2767.https://doi.org/10.3390/ijms20112767
-
Characterization of the interface between normal and transformed epithelial cellsNature Cell Biology 11:460–467.https://doi.org/10.1038/ncb1853
-
Signaling dynamics control cell fate in the early Drosophila embryoDevelopmental Cell 48:361–370.https://doi.org/10.1016/j.devcel.2019.01.009
-
Cell competition in mammals - novel homeostatic machinery for embryonic development and Cancer preventionCurrent Opinion in Cell Biology 48:106–112.https://doi.org/10.1016/j.ceb.2017.06.007
-
SoftwarefindPeaksMathWorks.
-
A monoclonal antibody to ADAM17 inhibits tumor growth by inhibiting EGFR and Non-EGFR-Mediated pathwaysMolecular Cancer Therapeutics 14:1637–1649.https://doi.org/10.1158/1535-7163.MCT-14-1040
-
Calcium wave promotes cell extrusionCurrent Biology 30:670–681.https://doi.org/10.1016/j.cub.2019.11.089
-
A stochastic model for leukocyte random motility and chemotaxis based on receptor binding fluctuationsThe Journal of Cell Biology 106:303–309.https://doi.org/10.1083/jcb.106.2.303
-
A role of the sphingosine-1-phosphate (S1P)-S1P receptor 2 pathway in epithelial defense against Cancer (EDAC)Molecular Biology of the Cell 27:491–499.https://doi.org/10.1091/mbc.e15-03-0161
-
The shedding protease ADAM17: physiology and pathophysiologyBiochimica Et Biophysica Acta (BBA) - Molecular Cell Research 1864:2059–2070.https://doi.org/10.1016/j.bbamcr.2017.07.001
Article and author information
Author details
Funding
National Institute of General Medical Sciences (1R35GM133499)
- Sergi Regot
National Science Foundation (MCB-1844994)
- Sergi Regot
American Cancer Society (133537-RSG-19-005-01-CCG)
- Sergi Regot
National Institute of General Medical Sciences (T32-GM007445)
- Timothy J Aikin
- Amy F Peterson
- Michael J Pokrass
- Helen R Clark
National Science Foundation (DGE-1746891)
- Timothy J Aikin
- Michael J Pokrass
- Helen R Clark
The Jerome L Greene Foundation (Jerome L Greene Foundation Discovery Award)
- Sergi Regot
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank all members of the Regot lab, and members of the Holland lab (A Holland, B Lambrus, and C Gleitch) and Wantanabe lab (T Ogunmowo) for helpful discussions and technical advice, and the Mass Spectrometry and Proteomics Core at JHSOM for assistance processing samples. We acknowledge our funding sources: NIH T32 pre-doctoral training grants to TA, AP, MP, and HC (GM007445). NSF Graduate Research Fellowships to TA, MP and HC (DGE-1746891). An NSF CAREER award (MCB-1844994), NIGMS R35 (1R35GM133499), American Cancer Society Research Scholar Grant (133537-RSG-19-005-01-CCG), and Jerome L Greene Foundation Discovery Award to SR.
Copyright
© 2020, Aikin 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,787
- views
-
- 504
- downloads
-
- 58
- 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
- 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.
-
- Cancer Biology
- Computational and Systems Biology
Effects from aging in single cells are heterogenous, whereas at the organ- and tissue-levels aging phenotypes tend to appear as stereotypical changes. The mammary epithelium is a bilayer of two major phenotypically and functionally distinct cell lineages: luminal epithelial and myoepithelial cells. Mammary luminal epithelia exhibit substantial stereotypical changes with age that merit attention because these cells are the putative cells-of-origin for breast cancers. We hypothesize that effects from aging that impinge upon maintenance of lineage fidelity increase susceptibility to cancer initiation. We generated and analyzed transcriptomes from primary luminal epithelial and myoepithelial cells from younger <30 (y)ears old and older >55y women. In addition to age-dependent directional changes in gene expression, we observed increased transcriptional variance with age that contributed to genome-wide loss of lineage fidelity. Age-dependent variant responses were common to both lineages, whereas directional changes were almost exclusively detected in luminal epithelia and involved altered regulation of chromatin and genome organizers such as SATB1. Epithelial expression of gap junction protein GJB6 increased with age, and modulation of GJB6 expression in heterochronous co-cultures revealed that it provided a communication conduit from myoepithelial cells that drove directional change in luminal cells. Age-dependent luminal transcriptomes comprised a prominent signal that could be detected in bulk tissue during aging and transition into cancers. A machine learning classifier based on luminal-specific aging distinguished normal from cancer tissue and was highly predictive of breast cancer subtype. We speculate that luminal epithelia are the ultimate site of integration of the variant responses to aging in their surrounding tissue, and that their emergent phenotype both endows cells with the ability to become cancer-cells-of-origin and represents a biosensor that presages cancer susceptibility.