- Open Access
A post-transcriptional program of chemoresistance by AU-rich elements and TTP in quiescent leukemic cells
Genome Biology volume 21, Article number: 33 (2020)
Quiescence (G0) is a transient, cell cycle-arrested state. By entering G0, cancer cells survive unfavorable conditions such as chemotherapy and cause relapse. While G0 cells have been studied at the transcriptome level, how post-transcriptional regulation contributes to their chemoresistance remains unknown.
We induce chemoresistant and G0 leukemic cells by serum starvation or chemotherapy treatment. To study post-transcriptional regulation in G0 leukemic cells, we systematically analyzed their transcriptome, translatome, and proteome. We find that our resistant G0 cells recapitulate gene expression profiles of in vivo chemoresistant leukemic and G0 models. In G0 cells, canonical translation initiation is inhibited; yet we find that inflammatory genes are highly translated, indicating alternative post-transcriptional regulation. Importantly, AU-rich elements (AREs) are significantly enriched in the upregulated G0 translatome and transcriptome. Mechanistically, we find the stress-responsive p38 MAPK-MK2 signaling pathway stabilizes ARE mRNAs by phosphorylation and inactivation of mRNA decay factor, Tristetraprolin (TTP) in G0. This permits expression of ARE mRNAs that promote chemoresistance. Conversely, inhibition of TTP phosphorylation by p38 MAPK inhibitors and non-phosphorylatable TTP mutant decreases ARE-bearing TNFα and DUSP1 mRNAs and sensitizes leukemic cells to chemotherapy. Furthermore, co-inhibiting p38 MAPK and TNFα prior to or along with chemotherapy substantially reduces chemoresistance in primary leukemic cells ex vivo and in vivo.
These studies uncover post-transcriptional regulation underlying chemoresistance in leukemia. Our data reveal the p38 MAPK-MK2-TTP axis as a key regulator of expression of ARE-bearing mRNAs that promote chemoresistance. By disrupting this pathway, we develop an effective combination therapy against chemosurvival.
Quiescent (G0) cells are an assortment of reversibly arrested cells, including dormant stem cells, which are found as a clinically relevant subpopulation in cancers [1,2,3,4]. Such cells are anti-proliferative, anti-differentiation, and anti-apoptotic and show distinct properties including resistance to harsh conditions [1, 2, 5,6,7,8,9,10]. G0 cells show specific gene expression that may underlie their resistance and other properties [1, 2, 8,9,10]. Analyses from multiple groups revealed some genes upregulated at the transcriptional level [1, 8, 11]. Altered polyadenylation site selection on mRNAs produces longer 3′-untranslated regions (3′UTRs) in G0 compared to proliferating cells—which increases 3′UTR elements that can mediate post-transcriptional gene expression regulation . Our previous data demonstrated that translation mechanisms are distinct in G0 leukemic cells, with decreased canonical translation mechanisms and increase in alternative mechanisms that involve non-canonical translation initiation factors  and 3′UTR-mediated specific mRNA translation . These data suggest that alternate post-transcriptional mechanisms in G0 cancer cells may regulate a distinct translatome to mediate their resistance. Translated genes in G0, the post-transcriptional mechanisms involved, and outcomes on cancer persistence remain to be investigated.
We analyzed the translatome and proteome of chemotherapy-surviving G0 cancer cells, focusing on acute monocytic leukemia (AML), to provide comprehensive information that complement and expand previous transcriptome analyses [1, 2, 8, 11, 15, 16], by uncovering critical genes that are post-transcriptionally regulated for chemosurvival. G0 can be induced by growth factor deprivation or serum starvation and other conditions that isolate dormant cancer stem cells in distinct cell types [1, 6, 7]. Our data demonstrate that serum starvation-induced G0 AML cells are chemoresistant—similar to surviving AML cells, isolated after chemotherapy. Chemoresistant cells isolated via serum starvation, or as surviving cells post-chemotherapy, show inhibition of canonical translation mechanisms, indicating that non-canonical mechanisms express specific mRNAs when these cells are chemoresistant. Consistently, the translatomes and proteomes of serum-starved G0 and chemosurviving cells show greater similarity than the transcriptomes alone. Our data reveal that DNA damage and stress signaling cause post-transcriptional alterations to produce a specialized gene expression program of pro-inflammatory, immune effectors that elicit chemosurvival.
Serum starvation or AraC treatment induces a quiescent and chemoresistant state of leukemic cells
To study clinical resistance in cancer, THP1 human AML cells were used as they show significant resistance to AraC  (cytosine arabinoside, Additional file 1: Figure S1A), a standard anti-leukemic chemotherapeutic that targets DNA replication and thus proliferating cells (referred to as S+). Our data and others find that serum starvation of THP1  and other cell lines [1, 8, 11, 18] induces a transient G0 state with known G0 and cell cycle arrest markers expressed (Fig. 1a and Additional file 1: Figure S1B-C). Such serum starvation-induced G0 cells (referred to as SS) can be returned to the cell cycle upon serum addition (Fig. 1b), verifying that they are quiescent and transiently arrested, unlike senescence or differentiation that are not easily reversed . We find that serum starvation-induced G0 SS cells show resistance to AraC chemotherapy. Serum-grown S+ cells show a dose-dependent decrease in cell viability with AraC as expected, while SS cells persist, indicating their chemoresistance (Fig. 1c). Chemoresistant cancer cells include cancer stem cells and are a subpopulation that can be isolated from cancers after treatment with chemotherapy [2, 6,7,8,9,10] that targets and eliminates S+ cells. We find that AraC-surviving THP1 (referred to as AraCS) cells are transiently arrested, like SS cells (Fig. 1b and Additional file 1: Figure S1B); both AraCS and SS cells survive chemotherapy (Fig. 1c). AraCS cells recover from their transient arrest upon AraC removal and proliferate (Fig. 1b), affirming the reversible G0 arrest state of chemoresistant cells, similar to SS cells [1, 2, 6,7,8,9,10].
G0 cells induced by SS or AraC have similar translatomes and proteome features that recapitulate gene expression profiles of in vivo chemoresistant leukemic and G0 models
To study post-transcriptionally regulated genes in G0, we profiled S+ cells, SS cells, and AraCS cells at the proteome, translatome, and transcriptome levels using multiplexed quantitative proteomics , microarray analysis of heavy polysome-associated mRNAs [13, 14, 19], and total RNAs respectively (Fig. 1d and Additional file 1: Figure S1D-E). Notably, we find that AraCS and SS cells show more similar gene expression profiles at the proteome and translatome levels, compared to transcriptome levels (Fig. 1e). These data suggest that although these chemoresistant G0 cells are isolated via two different methods, they exhibit a common set of translatome and proteome, which could underlie their common characteristic of chemoresistance. These data indicate the relevance of examining both the translatome and transcriptome. Time-course translatome analysis revealed that SS G0 cells that were serum-starved for short periods (4 h and 1 day) are distinct from SS G0 cells that were serum-starved for long periods (2 days and 4 days) (Additional file 1: Figure S1F). This is consistent with G0 as a continuum of assorted, arrested states , with temporal differences in the underlying gene expression in early G0, compared to more homogeneity at late G0. SS and AraCS cells provide sufficient material to perform concurrent translatome, proteome, and transcriptome profiling, compared to limited cells from in vivo resistance models where only transcriptomes were profiled. To test whether our G0 leukemic cells are relevant models to study chemoresistance and G0, gene expression profiles of AraCS and SS cells were compared to published transcriptome profiles of leukemia stem cells (LSC) from AML , dormant leukemic cells (LRC), and minimal residual disease (MRD) from chemotherapy surviving patient samples with acute lymphocytic leukemia (ALL) , as well as SS G0 fibroblasts (G0 HFF) . Importantly, we find that these published transcriptome signatures for in vivo chemoresistance and G0 models were significantly upregulated in our SS and AraCS cells (referred to as resistant G0 leukemic cells), compared to S+ cells (Fig. 1f and Additional file 1: Figure S1G). These data indicate that our resistant G0 leukemic cells are relevant models to study post-transcriptional regulation in chemoresistance as they have similar gene expression profiles to known transcriptional profiles from in vivo chemoresistance models.
Inhibition of canonical translation initiation in resistant G0 leukemic cells
Mechanistically, both rate-limiting steps in canonical translation initiation: recruitment of initiator tRNA, and mRNA cap recognition to recruit mRNAs to ribosomes, are inhibited in G0 leukemic cells (Fig. 2a–d) [13, 14]. We find overall protein synthesis is reduced at least twofold in AraCS, compared to S+ cells (Fig. 2b and Additional file 1: Figure S1D). Recruitment of initiator tRNA by eIF2 can be blocked by eIF2α phosphorylation as a stress response [13, 20,21,22,23,24,25]. We find that two eIF2 kinases, PKR and PERK, are activated and significantly increase eIF2α phosphorylation in SS and AraCS G0 leukemic cells (Fig. 2c, 5.3-fold (3.9/0.74) in SS and 4.2-fold (3.2/0.75) in AraCS cells of increase of phospho-eIF2α over total eIF2α, based on quantitation below the blots in Fig. 2c), which inhibits canonical translation initiation at one of the two rate-limiting steps. Consistent with our previous study , we observed moderate dephosphorylation in SS and AraCS G0 leukemic cells of eIF4EBP (4EBP, Fig. 2d) that can inhibit canonical translation initiation at the other rate-limiting step [26,27,28]. Decreased canonical translation by the above mechanistic changes can enable post-transcriptional regulation of specific genes, as observed previously [13, 14], and lead to survival of G0 leukemic cells.
Although 4EBP, a translation inhibitor that is a downstream target of mTOR signaling, is dephosphorylated in SS and AraCS cells (Fig. 2d), the effect is more pronounced for SS cells than AraCS cells (50–59% decreased phosphorylation of 4EBP in SS cells versus 31–33% in AraCS cells at two key sites; Thr37/40 and Ser65). This difference is reflected in the translation regulation of terminal oligopyrimidine tract (TOP) mRNAs, such as ribosomal protein mRNAs, that are regulated by mTOR [28,29,30]: TOP mRNAs are translationally decreased in SS cells (Additional file 1: Figure S1H) but not in AraCS cells. This indicates differences in SS and AraCS cells, in mTOR and other signals, which could differentially affect 4EBP [19, 28, 31] and other downstream effectors that are implicated in regulating TOP mRNAs and vary in distinct conditions [32,33,34,35]. Phosphorylation of the downstream mTOR target, RPS6 (S6), is moderately decreased in SS and AraCS leukemic cells (Additional file 1: Figure S1I, 36% in SS and 27% in AraCS cells). While the translation inhibitor, 4EBP that is downstream of the mTOR pathway, is dephosphorylated in these conditions and can partially reduce canonical translation (Fig. 2d), coordinate dephosphorylation of mTOR is not significantly observed (Additional file 1: Figure S1I) in SS (15% reduction at one site) and AraCS cells. This indicates that the mTOR pathway shows differences at different levels in SS and AraCS cells. This is likely due to feedback on mTOR from downstream S6K and other kinases, as observed in other systems [26, 36,37,38,39,40,41], while 4EBP is also known to be regulated by other kinases, independent of mTOR [42,43,44,45]. Thus, while the eIF2 pathway is strongly inhibited in SS and AraCS cells at similar levels (Fig. 2c), reducing canonical translation—and 4EBP is also moderately regulated with more dephosphorylation in SS than in AraCS cells, which could partially affect translation (Fig. 2d)—other levels of the mTOR pathway are not coordinately modulated or similarly altered in SS and AraCS cells.
Global translatome analysis shows that inflammatory response genes are selectively translated in resistant G0 cancer cells
We measured the number of genes upregulated at the transcriptome, translatome, and proteome levels in resistant G0 leukemic cells, compared to S+ cells. A significant number of genes were upregulated in the transcriptome (465 genes), the translatome (490 genes, Additional file 2: Table S1), and proteome (716 genes) as shown in Fig. 2e. Importantly, 67% of upregulated genes were upregulated only at the translatome level (Fig. 2f) but not in the transcriptome, indicating post-transcriptional regulation. To investigate the biological function of these differentially expressed genes, gene ontology (GO) analysis was performed. Gene categories upregulated in G0 translatomes include inflammatory response (pro-inflammatory cytokines and inflammation regulators), immune response genes (immune modulators that are not pro-inflammatory cytokines or inflammation regulators, including interferon-stimulated genes, immune receptors, antigen presentation, and processing genes), cell adhesion, cell migration, lipid biosynthesis, and cholesterol pathway genes (Fig. 2g and Additional file 1: Figure S2A-B). Downregulated genes include RNA processing and ribosome genes (Fig. 2g). To identify translationally upregulated genes, we measured the change in ribosome occupancy (RO) which is the ratio of polysome-associated mRNA levels to total mRNA levels of each gene (Fig. 2f heatmap and Additional file 2: Table S1). RO values are increased for some genes, indicating translational upregulation. These genes include antigen processing and presentation genes (HLA-G)  and immune receptors (CD47, Fig. 2g and Additional file 1: Figure S2C) [47,48,49] that regulate antitumor immune response and are associated with leukemic stem cells and resistance [50, 51].
We asked if this specific gene expression profile in resistant G0 leukemic cells is conserved in G0 cells of other tumors and cell types. Therefore, global translatome profiling was conducted in G0 cells from four different cell lines: breast cancer (MCF7), liver cancer (HEPG2), and osteosarcoma (U2OS) as well as non-cancerous fibroblasts (HFF) (Additional file 1: Figure S2D-G). Their translatome profiles were compared with resistant G0 leukemic cells, using GSEA and DAVID tools (Fig. 2h, i and Additional file 1: Figure S2A). We find that the 490 signature genes (upregulated translatome) of resistant G0 leukemic cells (Additional file 2: Table S1) were highly upregulated at the translatome level in G0 cells of these other cell types (Fig. 2h). As expected for these arrested cells, genes related to cell cycle, ribosome biogenesis, and DNA replication were commonly downregulated (Fig. 2i and Additional file 1: Figure S2A). We focused on inflammatory response genes as these were commonly upregulated in G0 cells from cancer cell lines and do not significantly overlap with the senescence-associated secretory pathway (SASP) (Fig. 2i and Additional file 1: Figure S2H) [52, 53].
Stabilization of ARE-bearing mRNAs is mediated by phosphorylation of TTP in resistant G0 leukemic cells
To identify cis-acting elements that mediate post-transcriptional regulation, the untranslated regions (UTRs) of differentially expressed genes were examined. We find that a GC-rich motif was enriched on 5′UTRs of translationally upregulated genes and an AU-rich motif, on 5′UTRs of downregulated genes, indicating that mRNAs with structured 5′UTRs are highly translated in G0 cells (Additional file 1: Figure S3A-B). Importantly, 3′UTR AU-rich elements (AREs) are significantly enriched in the upregulated translatome as well as transcriptome (Fig. 3a). Furthermore, 25% of the translatome signature of G0 leukemic cells bear AREs (Additional file 3: Table S2), including pro-inflammatory cytokines such as TNFα and chemokines (Fig. 3b, c) according to the ARE database . AREs are important post-transcriptional regulatory elements that mediate rapid degradation and repression of mRNAs . To understand how ARE mRNAs are highly expressed in G0 cells, we assessed the expression level of RNA-binding proteins. As expected, most ARE-binding proteins known to cause mRNA decay or translation repression [55, 56] are significantly reduced in G0 cells (Additional file 1: Figure S3C-D). Additionally, the exosome and proteasome complexes that are implicated in ARE mRNA decay [57, 58] are reduced (Additional file 1: Figure S3E-F). However, a key ARE mRNA decay factor, Tristetraprolin (TTP), was surprisingly increased in AraCS from multiple AML cell lines (Fig. 3d, e). However, we find that TTP is phosphorylated in SS and AraCS cells (Fig. 3e, right blot). TTP phosphorylation is established to increase its levels  and blocks its ability to destabilize ARE mRNAs, thus enabling ARE mRNA translation upon lipopolysaccharide (LPS) treatment in immune cells [60, 61]. To test whether phosphorylation of TTP was required for the increased expression of ARE mRNAs in G0 leukemic cells, we generated non-phosphorylatable mutant TTP with key phosphorylation sites (Ser 52, 178) replaced by alanine (TTP-AA). TTP-AA has been shown to mediate ARE mRNA decay activity and reduce pro-inflammatory cytokines like TNFα in immune cells, as it cannot be phosphorylated and inactivated [59,60,61]. Expression of myc-tagged TTP-AA significantly reduced TNFα mRNA in both THP1 and K562 AraCS cells (Fig. 3f), as it restored a form of TTP that cannot be phosphorylated and can pursue its decay function, unlike the endogenous TTP that gets phosphorylated and inactivated for its decay function. To determine the effect of TTP phosphorylation on the stability of ARE mRNAs, we measured the half-life of TNFα mRNA. Expression of TTP-AA mutant reduced the half-life of TNFα mRNA more significantly than TTP wild-type (TTP-WT) expressed in AraC-treated TTP-deficient cells (Fig. 3g). We compared the G0 translatome and transcriptome with TTP-CLIP datasets  to identify how many G0-expressed genes are TTP targets. The upregulated G0 translatome and RNA profiles (166 out of 490 translatome genes with p < 2.71e−21; 174 out of 465 RNA profile genes with p < 1.322e−26), and those with AREs (49–53%; 59 out of 121 translatome genes with p value < 7.302e−16; 75 out of 142 RNA profile genes with p < 1.535e−22) include known TTP targets (Additional file 1: Figure S3G, using hypergeometric probability test). Furthermore, immunoprecipitation demonstrated that TTP-AA was associated with TNFα mRNA in AraCS cells (Fig. 3h, GFP-tagged TTP-AA). In addition, to determine how many genes are regulated by TTP phosphorylation, we profiled AraC-treated cells that lack endogenous TTP but stably express TTP-AA or TTP-WT at the RNA level . Expression of TTP-AA mutant downregulated 58 genes at the mRNA level compared to TTP-deficient cells. TTP-CLIP data suggests that 40% of mRNAs affected by TTP-AA (23 out of 58, p < 3.585e−05, Additional file 1: Figure S3Hi) are associated with TTP. The other RNAs may be indirect targets, not directly associated with TTP. Furthermore, these mRNAs are stabilized by phosphorylation of TTP; expression of TTP-WT that allows TTP phosphorylation upregulates 53% of genes downregulated by TTP-AA (31 out of 58, p value < 0.05, fold change > 1.5, Additional file 1: Figure S3Hii and S3Hiii). Of these 58, at least 18 genes have AREs that are recorded in the ARE database  and are also stabilized by phosphorylation of TTP (S3Hiv). These data indicate that inactivation of the ARE mRNA decay function of TTP by TTP phosphorylation [59, 61, 62] is a key regulator of expression of a pro-inflammatory gene, TNFα, in chemoresistant G0 cells. These results are consistent with our findings of increased levels and translation of ARE-bearing mRNAs due to decreased ARE mRNA decay activity in G0 cells (Fig. 3a–c and Additional file 1: Figure S3C-F).
The p38 MAPK-MK2 pathway phosphorylates TTP to promote expression of ARE-bearing mRNAs in resistant G0 leukemic cells
To investigate how TTP is phosphorylated in resistant G0 leukemic cells, we examined key signaling molecules involved in DNA-damage response (DDR) (Fig. 4a) that is induced by chemotherapies like AraC [63,64,65,66]. As expected, AraC treatment induced rapid phosphorylation and activation of ATM (Fig. 4b and Additional file 1: Figure S4A). Importantly, we find that these conditions lead to phosphorylation and activation of p38 MAPK and its downstream effector, MAPKAPK2 (MK2) [67, 68] (Fig. 4b). MK2 has been shown to phosphorylate TTP in macrophages treated with LPS [59, 61, 62]. To examine whether the p38 MAPK-MK2 pathway phosphorylates TTP in resistant G0 leukemic cells, two different inhibitors of p38 MAPK were tested. Treatment with p38 MAPKα/β inhibitor, LY2228820 (LY) [68, 69], or a pan-p38 MAPK inhibitor that targets all isoforms, BIRB796 (BIRB) , blocked phosphorylation of MK2 and prevented MK2-mediated TTP phosphorylation and reduces TNFα in AraCS cells (Fig. 4c). These results suggest that p38 MAPK-MK2 phosphorylates and inactivates TTP, resulting in enhanced expression of ARE mRNAs such as TNFα upon AraC treatment (Fig. 4a). To test if the p38 MAPK-MK2-TTP pathway regulates TNFα expression via its ARE, a firefly luciferase reporter bearing the 3′UTR ARE of TNFα, and as control, Renilla luciferase, were co-transfected. Luciferase activity of the ARE reporter increased by twofold in AraCS cells compared to S+ cells but not when p38 MAPK was inhibited (Fig. 4d). These data suggest that the p38 MAPK-MK2-TTP axis upregulates expression of specific genes via AREs in G0 leukemic cells.
Phosphorylation of TTP induced by p38 MAPK-MK2 promotes chemoresistance
We noted that the p38 MAPK-MK2 pathway was rapidly activated to phosphorylate TTP within 1 day of SS or AraC treatment (Fig. 4b and Additional file 1: Figure S4A-B). To test the effect of inhibition of TTP phosphorylation on chemoresistance, p38 MAPK was inhibited before (or along with) as well as after treatment with AraC—and then, chemosurvival was measured using multiple assays, including cell death and two cell viability assays (Fig. 4e–g). Inhibition of p38 MAPK with BIRB or LY, 1 day after AraC treatment, when TTP was already phosphorylated, did not show any significant reduction in survival of AraC-resistant cells (Fig. 4f, g). Conversely, inhibition of p38 MAPK at earlier time points prior to AraC treatment, when TTP was not phosphorylated, increased apoptosis, and reduced survival of AraC-resistant cells (Fig. 4f, g). As a control, p38 MAPK inhibition alone does not affect viability of S+ cells that are not treated with AraC (Fig. 4f, g). These results suggest that p38 MAPK is rapidly activated upon AraC treatment to turn on downstream survival pathways such as phosphorylation of TTP. Thus, to inhibit phosphorylation of TTP and hence overcome AraC resistance effectively, p38 MAPK needs to be targeted at early time points.
To confirm that phosphorylation of TTP induces chemoresistance, we overexpressed TTP mutant (TTP-AA) that cannot be phosphorylated by p38 MAPK-MK2, followed by AraC treatment. Importantly, we find that TTP-AA mutant expression reduces survival of AraC-resistant cells in THP1 and K562 leukemic cell lines (Fig. 4h). Furthermore, TTP-AA mutant, expressed in TTP-knockout macrophages, induced apoptosis of AraC-surviving cells more significantly, compared to TTP wild-type (Fig. 4i). Consistently, in multiple AML cell lines, early inhibition of p38 MAPK showed dramatically reduced chemosurvival but not in non-cancerous CD34+ cells (Fig. 4j). When treated with p38 MAPK inhibitor alone, viability of S+ cells in multiple AML cell lines remained unchanged, indicating the synergism of AraC and p38 MAPK inhibitors (Fig. 4j). Interestingly, p38 MAPK inhibition eliminated resistant cells more significantly at increasing concentrations of AraC (Fig. 4k). This indicates that treatment with high concentrations of AraC would increase the number of cells induced into the resistant G0 state with strong phosphorylation of p38 MAPK-MK2-TTP. Conversely, even low concentrations of BIRB were sufficient to reduce chemoresistance (Additional file 1: Figure S4C). Unlike in solid tumors, where activation of p38 MAPK-MK2 induces resistance by arresting the cell cycle [30, 67, 68], p38 MAPK inhibition did not affect the cell cycle in AML cells (Additional file 1: Figure S4D). These data uncover rapid activation of a p38 MAPK-MK2 pathway that enables chemosurvival of G0 leukemic cells via inhibition of TTP activity.
TNFα, induced by phosphorylation of TTP, promotes chemoresistance
We demonstrated that TTP inactivation in SS and AraCS cells regulates the stability of ARE mRNAs such as TNFα in AraCS cells (Figs. 3g and 5a). This allowed such resistant G0 leukemic cells to show elevated TNFα translatome and protein levels (Fig. 5b, c). To assess the effect of TNFα on chemoresistance, we altered TNFα levels genetically and phamacologically in G0 cells. Induction of TNFα depletion prior to AraC effectively reduced AraC resistance, compared to depleting TNFα after AraC treatment, while no effect was observed with TNFα depletion alone without AraC (Fig. 5d). In contrast, addition of recombinant TNFα enhanced survival of AraCS cells (Fig. 5d). TNFα-mediated chemoresistance is not due to arrested cell cycle as TNFα treatment without subsequent AraC does not alter the cell cycle (Additional file 1: Figure S5A). These data suggest that phosphorylation of TTP and subsequent expression of TNFα, which are induced by p38 MAPK-MK2, are responsible for survival of G0 leukemic cells.
TNFα can also be inhibited pharmacologically with the drug pirfenidone (PFD) that can block TNFα translation in RAW264.7 cells and is used to treat idiopathic pulmonary fibrosis [68, 71, 72]. In G0 leukemic cells, PFD reduced TNFα translatome and protein levels but not mRNA levels (Fig. 5e and Additional file 1: Figure S5B). PFD treatment at least 18 h prior to or along with AraC or SS significantly reduced viability of G0 leukemic cells but failed to reduce resistance when added after AraC treatment (Fig. 5f and Additional file 1: Figure S5C). As observed with p38 MAPK-MK2 activation (Fig. 4a, b), TNFα translatome level also is rapidly and dramatically increased upon SS treatment (Fig. 5b). These data indicate that activation of TNFα is an early event in G0 induction, which leads to resistance, and needs to be inhibited early to preclude downstream survival regulators. PFD treatment alone does not affect the viability of untreated S+ cells, indicating that the cytotoxic effect of PFD is specific to G0 leukemic cells (Fig. 5f). PFD treatment reduced chemotherapy survival in multiple AML cell lines (Fig. 5g). Similar results were observed in MCF7 cells, where PFD reduced doxorubicin resistance (Additional file 1: Figure S5D).
TNFα activates the NFκB pathway that increases anti-apoptotic gene expression to promote cell survival [73,74,75]. Our observation of early activation of p38 MAPK-MK2 (Fig. 4a, b) suggested that TNFα could be rapidly upregulated upon G0 induction. Time-course translatome analysis affirmed that TNFα is highly increased (16-fold) at the earliest time point of 4 h after serum starvation or AraC treatment (Fig. 5b) along with its receptors, leading to rapid elevation of downstream NFκB target genes including antiapoptotic BCL family members [75,76,77] (Fig. 5b and Additional file 1: Figure S5E-F). Similar to our observations with TNFα inhibitor PFD (Fig. 5f), NFκB inhibitor, BAY11-7082 , prior to or along with AraC or SS decreases the viability of G0 cells, while treatment after AraC or SS had no effect (Fig. 5h). TNFα shRNA (Fig. 5d) or inhibition (Fig. 5f), and NFkB inhibition (Fig. 5h), effectively reduce survival of resistant cells as noted by viability assays. Apoptosis or caspase 3/7 activity mediated by TNFα shRNA or TNFα inhibitor, PFD, and NFkB inhibitor, BAY 11-7082, vary, although they all cause significant decrease in chemoresistant cell viability. While other pathways downstream of TNFα could affect apoptosis [79, 80], the differences in caspase activity can be due to differences in inhibition by shRNA depletion versus drug effects, as BAY 11-7082 can mediate NFκB-independent pathways and non-apoptotic cell death mechanisms [81, 82]. These data suggest that the TNFα-NFκB inflammatory pathway is upregulated as an early survival pathway in G0 cells.
TTP regulates a pro-apoptotic JNK pathway via targeting DUSP1
We asked what other ARE mRNAs are targeted by TTP and affect cell survival. DUSP1 mRNA contains AREs in its 3′ UTR. TTP has been shown to target DUSP1 mRNA for degradation upon LPS treatment of macrophages or dendritic cells [60, 61, 83]. Consistently, DUSP1 in AraC treatment in both THP1 and MOLM13 cells is decreased upon treatment with BIRB (Additional file 1: Figure S5G), indicating its regulation by p38 MAPK. To determine if TTP phosphorylation regulates DUSP1 in AraCS, we expressed TTP-AA mutant that is not phosphorylated in cells that lack TTP (Fig. 6a). Expression of TTP-AA mutant more significantly reduced DUSP1 mRNA and protein levels compared to cells expressing TTP wild-type (Fig. 6b, c). Immunoprecipitation showed that TTP-AA associated with DUSP1 mRNA in AraCS cells (Fig. 6d), similar to its association with TNFα mRNA (Fig. 3h). Furthermore, inhibition of phosphorylation of TTP by p38 MAPK inhibitor decreased DUSP1 protein level (Fig. 6e). DUSP1 is a MAPK phosphatase which dephosphorylates JNK . In AraCS cells, DUSP1 protein level is negatively correlated with phosphorylated JNK (Fig. 6e), consistent with DUSP1-mediated suppression of JNK . To determine the effect of JNK on survival of leukemic cells, JNK inhibitor, JNK-IN-8 was used (Fig. 6a). Importantly, JNK inhibition reversed apoptosis of leukemic cells treated with AraC, LY and PFD, but did not affect the viability of untreated cells (Fig. 6f, graph), indicating that inhibition of JNK pathway contributes to chemoresistance. Together, these results suggest that TTP-DUSP1 axis promotes chemoresistance via suppressing JNK-mediated apoptosis (Fig. 6a).
Co-inhibition of p38 MAPK and TNFα sensitizes resistant leukemic cells to AraC treatment
Although chemoresistant cells are sensitive to individual inhibition of either TNFα or p38 MAPK by PFD or LY respectively, a substantial proportion of cells still survived (Figs. 4g and 5f). Therefore, we asked if co-inhibition of p38 MAPK and TNFα with LY and PFD respectively, could eliminate the remaining resistant cells (Fig. 7a). We find that individual treatment with either of LY or PFD (at half the dosages used in Figs. 4g and 5f) prior to or along with AraC, reduces approximately 50% of surviving leukemic cells (Fig. 7b). Importantly, this combination of PFD and LY2228820 prior to AraC treatment—called PLA therapy—eliminates about 90% of chemoresistant cells in multiple AML cell lines and not just THP1 cells (Fig. 7a–c). Furthermore, PLA therapy decreased colony formation of leukemic cells on methylcellulose by 10-fold, compared to AraC treatment alone (Fig. 7d). These data indicate decreased survival of leukemic cells treated with PLA therapy. In contrast, in the absence of AraC treatment, the combination of PFD and LY2228820 did not affect cell viability, apoptosis and colony formation, indicating the synergistic effect between AraC and anti-inflammatory drugs (Fig. 7b–d). Despite the fact that stromal niche cells have been shown to protect leukemic cells from chemotherapy , we find that AML cells co-cultured with stromal cells remained sensitive to PLA therapy (Additional file 1: Figure S5H). We examined the molecular mechanism by which PLA therapy enhanced chemosensitivity. We find that LY treatment destabilizes TNFα mRNAs by TTP dephosphorylation  (Figs. 3g and 4c), while PFD suppresses translation of TNFα mRNA  (Fig. 5e and Additional file 1: Figure S5B). Therefore, in PLA therapy, TNFα remains more effectively blocked, compared to individual drug treatments (Figs. 6e and 7b). Furthermore, a pro-apoptotic JNK pathway was more significantly activated in cells treated with PLA therapy than single-drug treatments (Fig. 6e). Together, these results suggest that PLA therapy reduces TNFα and promotes a pro-apoptotic JNK pathway, leading to apoptosis of chemoresistant cells.
PLA therapy reduces chemoresistance in primary AML cells ex vivo and in vivo
To test the anti-leukemic activity of PLA therapy in primary AML , primary cells from AML patients (Fig. 8a and Additional file 1: Figure S5I) as well as two murine AML models driven by Hoxa9/Meis1 or MLL-AF9 (Fig. 8b), were used. When either p38 MAPK or TNFα was inhibited prior to AraC treatment, moderate apoptosis of chemoresistant cells was observed in primary AML cells (Fig. 8a, b and Additional file 1: Figure S5I). Importantly, co-inhibition of p38 MAPK and TNFα by PLA therapy (pre-treatment before AraC) significantly reduced AraC resistance in AML patient samples (Fig. 8a and Additional file 1: Figure S5I) as well as in primary cells from two AML mouse models ex vivo (Fig. 8b). In contrast, the viability of normal CD34+ cells from healthy donors was not affected by treatment with LY or PFD (Fig. 4j and Additional file 1: Figure S5I), consistent with clinical studies that have shown that PFD and LY have acceptable safety and tolerance [69, 71]. To further investigate the therapeutic potential of PLA therapy in vivo, human AML cells expressing luciferase (MOLM13-Luc) were intravenously or subcutaneously injected into NSG mice. After confirmation of engraftment by measuring tumor volume or bioluminescent imaging (BLI), the mice were treated with PLA therapy or AraC for two weeks. Consistent with ex vivo results (Fig. 7b), PLA therapy significantly decreased the leukemic burden and tumor volume by 6-fold, compared to AraC treatment alone (Fig. 8c, d). Next, primary Hoxa9/Meis1 or MLL-AF9 leukemia cells were generated as described previously , and transplanted to second recipient mice. These mice were treated with PLA therapy or AraC, with two different dosage schedules. Consistently, BLI shows that PLA therapy eliminated 78% or 96% of chemoresistant cells in a dosage-dependent manner (Fig. 8e, f; 8 times with drugs over 4 days versus 6 times with drugs over 2 weeks). In the absence of AraC treatment, the combination of PFD and LY2228820 did not affect leukemic burden, suggesting that cytotoxic effects of this combination are limited to AraC-resistant cells, rather than proliferating cells (Fig. 8g). Correspondingly, PLA therapy extended mice survival (Fig. 8h and Additional file 1: Figure S5J). Together, these results suggest PLA therapy has potential for improving AraC-mediated apoptosis in AML.
G0 cells are a transiently arrested, clinically relevant subpopulation in cancers [1, 2, 5,6,7,8,9,10]. Our previous data and others, revealed altered gene expression mechanisms in G0 leukemic cells, at the post-transcriptional [8, 12] and translational levels [13, 14, 18]. This would lead to a distinct gene expression profile to enable G0 cell survival in harsh conditions. G0 cells are resistant to stress conditions like serum starvation, with transient inhibition of apoptosis, and proliferation [1, 11, 18]. Importantly, we find that serum-starved leukemic SS G0 cells exhibit chemoresistance (Fig. 1c); consistently, true chemosurviving AraCS cells are transiently arrested and chemoresistant (Fig. 1a, b and Additional file 1: Figure S1B). In accord, we find that SS cells are similar in translatome and proteome to AraCS cells (Fig. 1e), indicating that consistent with their common features of G0 arrest and chemosurvival, they show similar post-transcription gene expression. Published transcriptional signatures of in vivo chemoresistance leukemic models [1, 2, 8, 11, 15, 16], are also highly expressed in SS and AraCS cells (Fig. 1f and Additional file 1: Figure S1G). Thus, the common G0 resistance gene expression profile observed in AraCS and SS G0 cells likely comprises genes that control survival and resistance. These data revealed that in addition to known transcriptional profiles, altered post-transcriptional mechanisms in G0 resistant cells contribute to their unique gene expression profile that underlies their chemoresistance.
Our findings reveal the importance of DNA damage and stress signaling that can initiate a pro-inflammatory response that causes survival (Fig. 4). Differential genomic instability in cancers would lead to subpopulations within a tumor with disparate DDR and stress signaling [63,64,65] that we find, enables their chemotherapy survival via pro-inflammatory cytokines. Cytokines upregulated in SS and AraCS cells include some SASP factors but also other unique cytokines [52, 53] (Additional file 1: Figure S2H). This is consistent with similarities and differences between G0 and senescence : both show inhibition of the cell cycle but—unlike in senescence—G0 shows reversible cell cycle arrest (Fig. 1a, b and Additional file 1: Figure S1B-C), increased stem cell markers (Fig. 1f and Additional file 1: Figure S1G), markers of maintenance of G0 such as HES1 (Additional file 1: Figure S1C, S2D)  that are not expressed in senescence , low p53 , and lack of common senescence markers (Additional file 1: Figure S2H) [18, 52, 53]. These data indicate that a quiescence- and resistance-specific set of pro-inflammatory and signaling genes are expressed in these resistant cells (Fig. 2g). These include inflammatory cytokine, TNFα, & its receptors that promote downstream NFκB activated pro-survival target genes [73,74,75] including BCL family members of antiapoptotic genes [75,76,77] (Fig. 5a–c and Additional file 1: Figure S5E-F). Treatment with anti-inflammatory reagents after chemotherapy is not very effective as the downstream survival effectors have already been induced (Additional file 1: Figure S5F); thus, targeting their upstream cytokine regulators would not be effective at this later time (Figs. 4f, g and 5f–h and Additional file 1: Figure S5C). Therefore, treatment with reagents that block these resistance pathways prior to (and continued with) or along with chemotherapy, enables the most effective reduction of resistance, as they prevent further enrichment of such resistant cells by blocking induction of pro-survival signaling.
Increasing AraC, a nucleotide analog that inhibits replication , would activate DDR and downstream p38 MAPK signaling [63,64,65] and should lead to more cells expressing this inflammatory pathway that enables resistance. Consistently, increased AraC treatment leads to more cells in the inflammatory phase that can be targeted by LY to curb resistance (Fig. 4k). Non-cancerous cells are not affected by these inhibitors (Fig. 4j and Additional file 1: Figure S5I). These data suggest that certain chemotherapies and stresses like serum starvation induce stress signaling (Fig. 4a–c and Additional file 1: Figure S4A-B) and enrich for resistant G0 cells—in addition to pre-existing subpopulations with genomic instability that trigger DDR and stress [63,64,65]. Importantly, this resistance mechanism can be blocked, not only in different AML cell lines (Figs. 4j, 5g, and 7b) but also in vivo (Fig. 8c–h) and in multiple patient-derived primary AML—without affecting normal cells (Fig. 8a and Additional file 1: Figure S5I)—supporting their potential applicability as a therapeutic against chemoresistance in AML.
We find key signaling pathways induced by AraCS and SS treatments, which alter post-transcriptional and translational gene expression to enable resistance. These include: 1. DNA damage ATM [63,64,65] and stress activated p38 MAPK that in turn promotes MK2 [67, 68], to post-transcriptionally upregulate ARE-bearing mRNAs [59, 61, 62]. The expressed mRNAs include ARE-bearing pro-inflammatory cytokine TNFα [73, 74] & its receptors that activates downstream anti-apoptosis signals (Fig. 4a–d and Additional file 1: Figure S4A-B, Fig. 5a–c and Additional file 1: Figure S5E-F) [75,76,77], and ARE-bearing signaling regulator DUSP1 [83, 84] that blocks JNK-mediated apoptosis (Fig. 6), to promote resistance. 2. UPR and PKR stress signaling are induced downstream of p38 MAPK  and DNA damage [89, 90], and inhibit canonical translation via PERK and PKR phosphorylation of eIF2α (Fig. 2a–c). This enables non-canonical translation of specific mRNAs when this rate-limiting step of canonical translation initiation is reduced [89, 90]. 3. In addition, DNA damage signaling can also cause suppression of the other rate-limiting step of canonical translation initiation, by dephosphorylation of 4EBP [63, 64]. Consistently, 4EBP dephosphorylation is observed here, although more moderately in AraCS cells compared to SS cells (Fig. 2d), and moderately in both conditions for S6, a second canonical translation regulator (Additional file 1: Figure S1I). These changes in post-transcriptional and translational mechanisms allow specific translation of pro-inflammatory cytokines  (Fig. 3a–c) and immune modulators  (HLA-G, CD47, Fig. 2f–g and Additional file 1: Figure S2C) [47,48,49] that regulate antitumor immune response and resistance [50, 51]. While the translation inhibitor, 4EBP that is downstream of the mTOR pathway, is dephosphorylated in these conditions and can reduce canonical translation (Fig. 2d), mTOR phosphorylation was not coordinately altered (Additional file 1: Figure S1I), indicating that the mTOR pathway shows differences at different levels in SS and AraCS cells. This is likely due to feedback regulation from S6K and other downstream kinases that can affect mTOR, as observed in other systems [26, 36,37,38,39,40,41], while 4EBP is also known to be regulated by other kinases, independent of mTOR [42,43,44,45].
Blocking the p38 MAPKα/β pathway with LY [68, 69] (Fig. 4c), in combination with the anti-inflammatory PFD [68, 71, 72] that precludes downstream TNFα expression [71, 72] (Fig. 5e)—prior to (and continued with) AraC chemotherapy—lead to effective loss of chemoresistance in multiple AML cell lines (Fig. 7b), in tumors in vivo in AML mouse models (Fig. 8c–h), and in patient samples (Fig. 8a and Additional file 1: Figure S5I), validating their ability to reduce resistance and tumors in vitro and in vivo. LY destabilizes TNFα mRNA by TTP dephosphorylation (Fig. 3g and 4c) , while PFD suppresses TNFα selectively at the translation level  (Fig. 5e and Additional file 1: Figure S5B) and thus enables PLA combination therapy to more effectively curb resistance than the individual drugs (Figs. 7b and 8a, b). Apart from its effect on TNFα translation, PFD blocks inflammation regulator p38 MAPKγ [91, 92] that can be increased upon p38MAPKα/β inhibition, preventing feedback reactivation of inflammation, and enabling PLA combination therapy to remain more efficacious than the individual drugs. Therefore, the combination of PFD and LY suppresses the inflammatory and stress response more effectively in vitro and in vivo (Figs. 7 and 8). Upon inhibition of p38 MAPK, in addition to reduction of TNFα and its downstream antiapoptotic signals, we find the ARE-bearing DUSP1 is reduced, leading to activation [83, 84] of the JNK pathway  to promote apoptosis (Fig. 6e, f). These data indicate that blocking pro-inflammatory effectors—that are induced by chemotherapy mediated DNA damage and stress signaling—leads to increased chemosensitivity and decreased resistant cell survival.
Our findings revealed that these pro-inflammatory and signaling genes upregulated in G0, have AREs and other UTR sequences that regulate mRNA levels and translation (Fig. 3a–c and Additional file 1: Figure S3A-B). The ATM-p38 MAPK-MK2 axis stabilizes these ARE-bearing pro-inflammatory cytokine and signaling mRNAs by phosphorylating ARE-binding mRNA decay factor, TTP, to prevent its mRNA decay activity on pro-inflammatory cytokine TNFα (Figs. 3e–h and 4c, d, h, i) and signaling regulator, DUSP1 (Fig. 6a–e and Additional file 1: Figure S5G). There may be additional contributors to TTP besides MK2: p38 MAPK also directly phosphorylates TTP in macrophages [94, 95] while MEKK1 can act as a TTP kinase with TNF receptor-associated factor 2 (TRAF2) . These could be involved in prolonged serum starvation, as phosphorylation of p38 MAPK is reduced after 1 day of serum starvation. In addition, protein phosphatase (PP2A) is known to dephosphorylate TTP. If PP2A activity is reduced in prolonged serum starvation, TTP can remain phosphorylated.
In support of the critical role of TTP regulation in chemoresistance, overexpression of TTP-AA—that cannot be phosphorylated and is a dominant active form that restores ARE mRNA decay [59,60,61]—decreases TNFα and DUSP1 expression (Figs. 3f–h and 6a–d), and thereby reduces chemoresistance (Figs. 4h, i and 6e, f). This is consistent with previous studies on AREs in cancers [14, 30, 59, 97,98,99,100]. These data suggest that phospho-TTP level or TTP activity is an important regulator of inflammatory response-mediated chemoresistance, which can be harnessed as a marker and target against AML resistance. Consistently, published in vivo leukemia resistance models show increased expression of TTP and ARE-bearing genes [15, 101], similar to our studies (Fig. 3a–e). Our studies on TTP and ARE regulated immune and signaling modulators that promote chemoresistance, are consistent with recent findings of TTP regulation of PDL1 to mediate immuno-resistance in solid tumors . Importantly, inhibition of these pathways curtails chemoresistance and tumor survival in vivo in primary AML patients and tumor models (Fig. 8 and Additional file 1: Figure S5I-J). Together, these pathways that are upregulated in resistant cells (Figs. 4a, 5a, and 6a) via chemotherapy and stress-induced signaling—decrease canonical translation and permit non-canonical post-transcriptional regulation of specific genes (Additional file 1: Figure S6)—to promote chemosurvival of G0 cancer cells.
Our studies reveal that G0 leukemic cells are chemoresistant, indicating their clinical importance in cancer persistence. We find a specific proteomic and translation profile that is induced commonly between G0 cells and chemosurviving leukemic cells. We uncovered critical genes that are upregulated post-transcriptionally for cell survival in these conditions by key, survival signaling pathways. These studies reveal the significance of post-transcriptional regulation of pro-inflammatory genes and signaling modulators in chemoresistance in leukemia. Our data enabled the development of a new combination therapy to effectively reduce resistance in cancer cell lines, in tumors in vivo, and in patient tumor samples.
Overview, aim, design, and setting
Therapeutic targeting of minimal residual disease or chemoresistant, leukemic stem cells in leukemias, particularly acute myeloid leukemia, has been ineffective thus far and refractory leukemia is fatal. The mechanisms of translation and post-transcriptional control, and the critical translation profile that control the ultimate, specific protein profile, and thereby survival of such clinically resistant cells, are largely undiscovered. Therefore, we globally analyzed gene expression at every level—RNA levels, translatome, and proteome—in chemotherapy-surviving G0 cancer cells in acute monocytic leukemia and other cancers, the specialized post-transcriptional and translational mechanistic changes, and their key signaling regulatory pathways, as well as developed a new, resistance-gene expression targeting therapy to understand and reduce chemoresistance.
Detailed description of characteristics, materials used, methods, and statistical analyses including cell culture, patient samples, tumor models, profiling, plasmids, cell viability assays, flow cytometry, protein analysis, drugs, and motif analysis is described in detail below.
THP1 cells were cultured in Roswell Park Memorial Institute (RPMI)1640 media supplemented with 10% fetal bovine serum (FBS), 2 mM l-glutamine, 100 μg/mL streptomycin, and 100 U/ml penicillin at 37 °C in 5% CO2. SS THP1 cells were prepared by washing with PBS followed by serum starvation at a density of 2 × 105 cells/mL and AraCS cells, by treatment with 5 μM AraC for 3 days or 9 days. MCF7, HFF, HEPG2, and U2OS cells were cultured in Dulbecco’s modified Eagle’s medium (DMEM) media with 10% FBS, 2 mM l-glutamine, 100 μg/mL streptomycin, and 100 U/ml penicillin, as done previously [13, 14]. MCF7 cells were serum-starved or treated with 150 μM doxorubicin. THP1 (TIB-202), MV4:11 (CRL-9591), K562 (CCL243), HFF (SCRC-1041), MCF7 (HTB-22), U2OS (HTB-96), and HEPG2 (HB-8065) were obtained from ATCC. MOLM13 (ACC554), NOMO1 (ACC542) and MONOMAC6 (ACC124) were obtained from DSMZ. Cell lines kindly provided by David Scadden  and MOLM13-GFP-Luc by Monica Guzman . As previously described [60, 61], we used bone marrow-derived macrophages (BMDMs) transduced with plasmids coding for doxycycline-inducible GFP-TTP, GFP-TTP-AA, or GFP. Cell lines were tested for Mycoplasma (Promega) and authenticated by the ATCC Cell Authentication Testing Service .
Primary AML patient samples and human monocytes
All human samples (de-identified) were handled in accordance with IRB protocols to SV (2015P000998/MGH), approved by the Partners Human Research Committee Institutional Review Board (IRB)/MGH IRB, to DAS and to TG (DF/HCC 13-583), approved by DF/HCC Office for Human Research Studies. AML samples used in this study were obtained by DAS including the following: MGH15—bone marrow 60% blasts, karyotype 46, XX, t(9;11)(p22;q23)[20/20]; MGH22—peripheral blood, 60% blasts, karyotype 46,XX,t(3;21)(q26;q22),t(9;22)(q34;q11.2) /46,XX ; and MGH25—bone marrow, 90% blasts, karyotype 46,XX  and by JL-S and TG including bone marrow samples: EQ1899, CI2095, PO2038, LA2053, NC1866, GO1122, CM2164, MV2192, VD2160, XD2101, VL2317, and OA2500. Bone marrow or peripheral blood mononuclear cells were isolated from de novo AML patients by ficoll density gradient centrifugation and cryopreserved with DMSO in a liquid nitrogen tank. Thawed cells were maintained in RPMI media with 10% FBS for several days before drug treatment and analyses. Human CD34+ monocytes (2 M-101) were obtained from Lonza. Primary cells from MLL-AF9, HoxA9/Meis1 mouse models were provided by DS . Mouse primary cells were maintained in RPMI media with 10% FBS, 2 mM l-glutamine, 100 μg/mL streptomycin, 100 U/ml penicillin, 5 ng/ml murine IL-3, and 25 ng/ml murine stem cell factor (SCF).
In vivo AML mouse models
AML mouse models have been shown to predict therapy response accurately . C57BL/6 and NSG were obtained from MGH Cox-7 Gnotobiotic animal facility of the AAALAC-accredited Center for Comparative Medicine and Services at MGH. C57BL/6 or NSG mice were injected intravenously or subcutaneously with MOLM13 cells expressing luciferase or intravenously with HoxA9/Meis1 or MLL-AF9 [103, 105]. IVIS imaging system (Perkin Elmer) were used to confirm engraftment of AML cells. Mice were intraperitoneally injected with 200 μl of luciferase substrate D-Luciferin (15 mg/ml) and anesthetized. Images were taken 5 or 10 min after D-Luciferin injection. After confirmation of engraftment by IVIS imaging, mice were randomly assigned to two groups and treated with pirfenidone (100 mg/kg, intraperitoneally), LY2228820 (20 mg/kg, intraperitoneally), AraC (30 mg/kg, intraperitoneally), or saline according at indicated combinations and dosages. Tumor volumes were measured by IVIS imaging at indicated time points.
Polysome profiling with microarray
Sucrose was dissolved in lysis buffer containing 100 mM KCl, 5 mM MgCl2, 100 μg/ml cycloheximide, 2 mM DTT, and 10 mM Tris-HCl (pH 7.4). Sucrose gradients from 15 to 50% were prepared in ultracentrifuge tubes (Beckman) as previously described [13, 106,107,108]. Cells were treated with 100 μg/mL cycloheximide at 37 °C for 5 min before collecting them. Harvested cell were rinsed with ice-cold PBS having 100 μg/mL cycloheximide and then were resuspended in lysis buffer with 1% Triton X-100 and 40 U/mL murine (New England Biolabs) for 20 min. After centrifugation of cell lysates at 12,000×g for 20 min, supernatants were loaded onto sucrose gradients followed by ultracentrifugation (Beckman Coulter Optima L90) at 34,000 × rpm at 4 °C for 2 h in the SW40 rotor. Samples were separated by density gradient fractionation system (Teledyne Isco). RNAs were purified by using TRIzol (Invitrogen) from heavy polysome fractions and whole cell lysates. The synthesized cDNA probes from WT Expression Kit (Ambion) were hybridized to Gene Chip Human Transcriptome Array 2.0 (Affymetrix) and analyzed by the Partners Healthcare Center for Personalized Genetic Medicine Microarray and BUMC facilities. Gene ontology analysis for differentially expressed translatome or proteome was conducted by DAVID 6.7 tools [109, 110]. Molecular signatures enriched in AraCS or SS were identified by Gene Set Enrichment Analysis (GSEA) .
TRIPZ plasmids expressing shRNA against human TNFα (V2THS_111606), and miR30a primiR sequences used as control (RHS4750), were obtained from Open Biosystems and MGH cancer center, respectively. Stable cell lines were constructed as described by Open Biosystems. The stable cells expressing shRNA against TNFα were induced with 1 μg/mL doxycycline at indicated time points to knockdown TNFα. Luciferase reporters to test ARE expression were previously described . Cells were treated with 10 ng/ml recombinant TNFα (R&D Systems) to activate the NFκB pathway. Myc-tagged TTP-AA [112, 113] was a gift from Nancy Kedersha and Shawn Lyons from Paul Anderson’s lab.
MTS assay, a colorimetric quantification of viable cells, was conducted as described by the manufacturer, Promega. A volume of 100 μl cells was placed in a 96-well plate after drug treatment. A volume of 20 μl MTS reagent (CellTiter 96® Aqueous Non-Radioactive Cell Proliferation Assay) was added to each well followed by incubation at 37 °C for 1 h. Absorbance was measured at 490 nm by using a microplate reader.
Caspase 3/7 assay
After drug treatment, cell death was measured by using caspase-glo® 3/7 assay kit (Promega) according to the protocol provided by the manufacturer. The equal volume of caspase-glo reagent was added to cells, and samples were gently mixed with pipetting. The plates were incubated at room temperature in the dark for 2 h. The luminescence of each sample was measured in a luminometer (Turner BioSystems).
Flow cytometry and cell cycle analysis
Cell proliferation was determined by flow cytometry of cells labeled with propidium iodide and bromodeoxyuridine (BrdU). The cells were incubated with 10 μM BrdU for 90 min at 37 °C in 5% CO2 before harvesting. Collected cells were fixed in ice-cold 70% ethanol overnight. Cells were washed in PBS and treated with 2 M HCl for 30 min. Cells were incubated for 1 h with anti-BrdU antibody conjugated to FITC (eBioscience) in the dark, washed, and stained with propidium iodide. Samples were filtered through a nylon mesh filter and cell cycle analysis performed on the flow cytometry .
Western blot analysis
Cells were collected and resuspended in lysis buffer containing 40 mM Tris-HCl (pH 7.4), 6 mM MgCl2, 150 mM NaCl, 0.1% NP-40, 1 mM DTT, and protease inhibitors (Roche). Samples containing 80 μg of protein were loaded onto 10% or 12% SDS-PAGE (Bio-Rad), transferred to Nitrocellulose membranes and processed for immunoblotting. Antibodies against p27 (#06–445) and tubulin (#05–829) were obtained from Millipore. Antibodies against HES1 (#sc-25,392), eIF2α (#sc-11,386), and GFP (#sc-9996) were from Santa Cruz. Antibodies against phospho-ATM (#ab81292), phospho-PKR (#ab32036), DUSP1 (#ab138265), and phospho-IRE1 (#ab124945) were from Abcam. Antibody against RPS6 (#66886-1-Ig) was from Proteintech. Antibody against phospho-PERK (#649401) was from Biolegend. Antibodies against phospho-mTOR (Ser2448, #2971), phospho-mTOR (Ser2481, #2974), mTOR (#2983), phospho-4E-BP1 (Thr37/46, #2855), phospho-4E-BP1 (Ser65, #9451), phospho-S6 ribosomal protein (Ser235/236, #2211),TNFα (#3707), phospho-p38 MAPK (#4511), phospho-MK2 (#3007), phospho-eIF2α (#9721), TTP (#71632), JNK (#9252), phospho-JNK (#9251), and 4EBP1 (#9452) were from Cell Signaling Technology.
Total RNA was isolated using TRIzol (Invitrogen) according to the manufacturer’s instructions. The cDNA was synthesized from 1 μg of RNA using M-MuLV Reverse Transcriptase (NEB) and random hexamer primer (Promega). qPCRs were run on LightCycler® 480 Instrument II (Roche) using 2X SYBR green mix (Bio-rad). The primers used in the qPCR were as follows: mouse TNF-α sense 5′-GCCTCTTCTCATTCCTGCTTG-3′, antisense 5′-CTGATGAGAGGGAGGCCATT-3′; mouse Gapdh sense 5′-CATGGCCTTCCGTGTTCCT-3′, antisense 5′-TGATGTCATCATACTTGGCAGGTT-3′; Dusp1 sense 5′-GGCCAGCTGCTGCAGTTTGAG-3′, antisense 5′-AGGTGCCCCGGTCAAGGACA-3′.
Leukemic cells were treated with indicated drug combinations. Annexin V FITC/PI staining was performed with FITC Annexin V Apoptosis Detection Kit I (BD Pharmingen). Flow cytometry analysis and FlowJo software were used to quantify the percentages of apoptotic cells.
After treatment with indicated drug combinations, the same number of cells was plated in methylcellulose-based media with human recombinant cytokines (stem cell technology, MethoCult™ H4435). Number of colonies was quantified in each plate after 10 days.
Multiplex quantitative proteomics analysis was conducted, as performed previously , in S+-, SS-, and AraC-treated THP1 leukemic cells.
Expression of GFP-tagged TTP-AA mutant was induced with 1 μg/ml doxycycline prior to 1 μM AraC treatment in TTP-deficient BMDM cells. The cells were cross-linked with UV 254 nm. Cells were lysed in lysis buffer (20 mM Tris, pH 7.5, 150 mM NaCl, 1 mM EDTA, 1 mM EGTA, 1% Triton X-100, 2.5 mM sodium pyrophosphate, 1 mM β-glycerophosphate, 1 mM Na3VO4, protease inhibitor, RNase inhibitor). Cell lysates were incubated overnight at 4 °C with either IgG control or GFP antibody. Protein G agarose (Santa Cruz) was used to pull down antibody bound RNA-protein complexes.
Pirfenidone (10 to 300 μg/ml [71, 116,117,118]) was obtained from Chemietek. AraC (1 to 10 μM [119, 120]), LY2228820 (0.03 to 2 μM [68, 69, 121, 122]), BIRB796 (BIRB, 5 μM [70, 123,124,125,126]), and JNK-IN-8 (1 μM  were from Selleckchem. KU55933 (10 μM , tested but was toxic for the cells tested), BAY 11-7082 (10 μM ), and D-luciferin were from Cayman Chemical and doxorubicin (10 to 500 nM ) was from Tocris Bioscience.
Motif, AREs, RNA-binding proteins, ribosome occupancy, and GSEA analysis
The Multiple Em for Motif Elicitation (MEME) software was used to search for cis-elements enriched in 5′ UTR of translationally regulated genes . Human 5′ UTR sequences were retrieved from UCSC table browser . In a discriminative mode, 5′ UTR sequences of translationally up- or downregulated genes were used as the primary sequences and 5′ UTR sequences of translationally unchanged genes, the control sequences. Motifs were found in the given strand with 6–30 nt motif width. We compared polysome-associated mRNAs with their total RNA levels in serum-starved and AraCS cells to generate the change in ribosome occupancy (RO) [132,133,134]—which is the ratio of the level of mRNA that is associated with heavy polysomes compared to the total mRNA level of each gene (Fig. 2f, heat map, Additional file 2: Table S1). ARE Score algorithm  was used to assess scores of AU-rich elements quantitatively. The list of RNA-binding protein genes were obtained from RBPDB database . Gene Set Enrichment Analysis (GSEA) was performed using all 50 gene sets of the Hallmarks, and gene sets from KEGG, reactome, and GO pathways from the Molecular Signatures Database (MSigDB) [111, 137] with our transcriptome, translatome, and proteome datasets.
Statistical tests and differential gene expression analyses
All experiments were performed with at least three replicates except for experiments with AML patients. Sample sizes were estimated on the basis of availability and previous experiments [13, 14]. No samples were excluded from analyses. Statistical methods were not used to pre-determine sample size. Two-tailed t-test and Wilcoxon rank sum test were performed for statistical tests. SEM (standard error of mean) values are shown as error bars in all figures. Means were used as center values in box plots. p values less than 0.05 were indicated with an asterisk. E-values were used for the statistical significance in the motif analysis. Affymetrix microarray data were normalized and summarized using the RMA method implemented in the affy R package . Genes with small variation or a consistently low signal across samples were filtered by the varFilter function in the genefilter package. A robust linear regression model was then used to fit to the probe intensities using the lmFit function, followed by the detection of differentially expressed genes using the eBayes function in the limma R package . Differentially expressed genes were identified using p < 0.05 and log2 fold change of ± 0.585 (1.5-fold change). The statistical significance of overlaps between two different groups of genes was assessed using hypergeometric probability test (http://nemates.org/MA/progs/overlap_stats.cgi), (Additional file 1: Figure S3G), with the total number of proteins being 26,809, based on our arrays.
Availability of data and materials
Raw datasets are available on the public repository, GEO, with series accession numbers GSE141075 and GSE141329 which are included in Super Series GSE141332 .
All materials will be made available publicly on publication and on request.
Coller HA, Sang L, Roberts JM. A new description of cellular quiescence. PLoS Biol. 2006;4(3):e83.
Ng SW, Mitchell A, Kennedy JA, Chen WC, McLeod J, Ibrahimova N, et al. A 17-gene stemness score for rapid determination of risk in acute leukaemia. Nature. 2016;540(7633):433–7.
Meacham CE, Morrison SJ. Tumor heterogeneity and cancer cell plasticity. Nature. 2013;501(7467):328–37.
Crews LA, Jamieson CH. Selective elimination of leukemia stem cells: hitting a moving target. Cancer Lett. 2013;338(1):15–22.
Bhola PD, Mar BG, Lindsley RC, Ryan JA, Hogdal LJ, Vo TT, et al. Functionally identifiable apoptosis-insensitive subpopulations determine chemoresistance in acute myeloid leukemia. J Clin Invest. 2016;126(10):3827–36.
Tavaluc RT, Hart LS, Dicker DT, El-Deiry WS. Effects of low confluency, serum starvation and hypoxia on the side population of cancer cell lines. Cell Cycle. 2007;6(20):2554–62.
Gupta PB, Onder TT, Jiang G, Tao K, Kuperwasser C, Weinberg RA, et al. Identification of selective inhibitors of cancer stem cells by high-throughput screening. Cell. 2009;138(4):645–59.
Salony, Sole X, Alves CP, Dey-Guha I, Ritsma L, Boukhali M, et al. AKT inhibition promotes nonautonomous cancer cell survival. Mol Cancer Ther. 2016;15(1):142–53.
Li L, Bhatia R. Stem cell quiescence. Clin Cancer Res. 2011;17(15):4936–41.
Giles FJ, DeAngelo DJ, Baccarani M, Deininger M, Guilhot F, Hughes T, et al. Optimizing outcomes for patients with advanced disease in chronic myelogenous leukemia. Semin Oncol. 2008;35(1 Suppl 1):S1–17.
Liu H, Adler AS, Segal E, Chang HY. A transcriptional program mediating entry into cellular quiescence. PLoS Genet. 2007;3(6):e91.
Sandberg R, Neilson JR, Sarma A, Sharp PA, Burge CB. Proliferating cells express mRNAs with shortened 3′ untranslated regions and fewer microRNA target sites. Science. 2008;320(5883):1643–7.
Lee S, Truesdell SS, Bukhari SI, Lee JH, Letonqueze O, Vasudevan S. Upregulation of eIF5B controls cell-cycle arrest and specific developmental stages. Proc Natl Acad Sci U S A. 2014;111(41):E4315–E22.
Bukhari SI, Truesdell SS, Lee S, Kollu S, Classon A, Boukhali M, et al. A specialized mechanism of translation mediated by FXR1a-associated MicroRNP in cellular quiescence. Mol Cell. 2016;61(5):760–73.
Ebinger S, Ozdemir EZ, Ziegenhain C, Tiedt S, Castro AC, Grunert M, et al. Characterization of rare, dormant, and therapy-resistant cells in acute lymphoblastic leukemia. Cancer Cell. 2016;30(6):849–62.
Saito Y, Kitamura H, Hijikata A, Tomizawa-Murasawa M, Tanaka S, Takagi S, et al. Identification of therapeutic targets for quiescent, chemotherapy-resistant human leukemia stem cells. Sci Transl Med. 2010;2(17):17ra9.
Forbes SA, Beare D, Gunasekaran P, Leung K, Bindal N, Boutselakis H, et al. COSMIC: exploring the world's knowledge of somatic mutations in human cancer. Nucleic Acids Res. 2015;43(Database issue):D805–D11.
Loayza-Puch F, Drost J, Rooijers K, Lopes R, Elkon R, Agami R. p53 induces transcriptional and translational programs to suppress cell proliferation and growth. Genome Biol. 2013;14(4):R32.
Hsieh AC, Liu Y, Edlind MP, Ingolia NT, Janes MR, Sher A, et al. The translational landscape of mTOR signalling steers cancer initiation and metastasis. Nature. 2012;485(7396):55–61.
Sendoel A, Dunn JG, Rodriguez EH, Naik S, Gomez NC, Hurwitz B, et al. Translation from unconventional 5′ start sites drives tumour initiation. Nature. 2017;541(7638):494–9.
Holcik M, Sonenberg N. Translational control in stress and apoptosis. Nat Rev Mol Cell Biol. 2005;6(4):318–27.
Zeenko VV, Wang C, Majumder M, Komar AA, Snider MD, Merrick WC, et al. An efficient in vitro translation system from mammalian cells lacking the translational inhibition caused by eIF2 phosphorylation. RNA. 2008;14(3):593–602.
Lorsch JR, Dever TE. Molecular view of 43 S complex formation and start site selection in eukaryotic translation initiation. J Biol Chem. 2010;285(28):21203–7.
Ron D, Walter P. Signal integration in the endoplasmic reticulum unfolded protein response. Nat Rev Mol Cell Biol. 2007;8(7):519–29.
Zismanov V, Chichkov V, Colangelo V, Jamet S, Wang S, Syme A, et al. Phosphorylation of eIF2alpha; is a translational control mechanism regulating muscle stem cell quiescence and self-renewal. Cell Stem Cell. 2016;18(1):79–90.
Sonenberg N, Hinnebusch AG. Regulation of translation initiation in eukaryotes: mechanisms and biological targets. Cell. 2009;136(4):731–45.
Culjkovic B, Topisirovic I, Borden KL. Controlling gene expression through RNA regulons: the role of the eukaryotic translation initiation factor eIF4E. Cell Cycle. 2007;6(1):65–9.
Thoreen CC, Chantranupong L, Keys HR, Wang T, Gray NS, Sabatini DM. A unifying model for mTORC1-mediated regulation of mRNA translation. Nature. 2012;485(7396):109–13.
Han K, Jaimovich A, Dey G, Ruggero D, Meyuhas O, Sonenberg N, et al. Parallel measurement of dynamic changes in translation rates in single cells. Nat Methods. 2014;11(1):86–93.
Damgaard CK, Lykke-Andersen J. Regulation of ARE-mRNA stability by cellular signaling: implications for human cancer. Cancer Treat Res. 2013;158:153–80.
Alain T, Morita M, Fonseca BD, Yanagiya A, Siddiqui N, Bhat M, et al. eIF4E/4E-BP ratio predicts the efficacy of mTOR targeted therapies. Cancer Res. 2012;72(24):6468–76.
Damgaard CK, Lykke-Andersen J. Translational coregulation of 5'TOP mRNAs by TIA-1 and TIAR. Genes Dev. 2011;25(19):2057–68.
Tcherkezian J, Cargnello M, Romeo Y, Huttlin EL, Lavoie G, Gygi SP, et al. Proteomic analysis of cap-dependent translation identifies LARP1 as a key regulator of 5'TOP mRNA translation. Genes Dev. 2014;28(4):357–71.
Miloslavski R, Cohen E, Avraham A, Iluz Y, Hayouka Z, Kasir J, et al. Oxygen sufficiency controls TOP mRNA translation via the TSC-Rheb-mTOR pathway in a 4E-BP-independent manner. J Mol Cell Biol. 2014;6(3):255–66.
Fonseca BD, Zakaria C, Jia J-J, Graber TE, Svitkin Y, Tahmasebi S, et al. La-related protein 1 (LARP1) represses terminal oligopyrimidine (TOP) mRNA translation downstream of mTOR complex 1 (mTORC1). J Biol Chem. 2015;290(26):15996–6020.
Koritzinsky M, Magagnin MG, van den Beucken T, Seigneuric R, Savelkouls K, Dostie J, et al. Gene expression during acute and prolonged hypoxia is regulated by distinct mechanisms of translational control. EMBO J. 2006;25(5):1114–25.
Bodur C, Kazyken D, Huang K, Ekim Ustunel B, Siroky KA, Tooley AS, et al. The IKK-related kinase TBK1 activates mTORC1 directly in response to growth factors and innate immune agonists. EMBO J. 2018;37(1):19–38.
Caron E, Ghosh S, Matsuoka Y, Ashton-Beaucage D, Therrien M, Lemieux S, et al. A comprehensive map of the mTOR signaling network. Mol Syst Biol. 2010;6:453.
Foster KG, Fingar DC. Mammalian target of rapamycin (mTOR): conducting the cellular signaling symphony. J Biol Chem. 2010;285(19):14071–7.
Levin VA, Panchabhai SC, Shen L, Kornblau SM, Qiu Y, Baggerly KA. Different changes in protein and phosphoprotein levels result from serum starvation of high-grade glioma and adenocarcinoma cell lines. J Proteome Res. 2010;9(1):179–91.
Pirkmajer S, Chibalin AV. Serum starvation: caveat emptor. Am J Physiol Cell Physiol. 2011;301(2):C272–9.
Hay N, Sonenberg N. Upstream and downstream of mTOR. Genes Dev. 2004;18(16):1926–45.
Villalonga P, Fernandez de Mattos S, Ridley AJ. RhoE inhibits 4E-BP1 phosphorylation and eIF4E function impairing cap-dependent translation. J Biol Chem. 2009;284(51):35287–96.
Nawroth R, Stellwagen F, Schulz WA, Stoehr R, Hartmann A, Krause BJ, et al. S6K1 and 4E-BP1 are independent regulated and control cellular growth in bladder cancer. PLoS One. 2011;6(11):e27509.
Batool A, Aashaq S, Andrabi KI. Reappraisal to the study of 4E-BP1 as an mTOR substrate - a normative critique. Eur J Cell Biol. 2017;96(4):325–36.
Shukla SA, Rooney MS, Rajasagi M, Tiao G, Dixon PM, Lawrence MS, et al. Comprehensive analysis of cancer-associated somatic mutations in class I HLA genes. Nat Biotechnol. 2015;33(11):1152–8.
Sosale NG, Spinler KR, Alvey C, Discher DE. Macrophage engulfment of a cell or nanoparticle is regulated by unavoidable opsonization, a species-specific ‘Marker of Self’ CD47, and target physical properties. Curr Opin Immunol. 2015;35:107–12. https://doi.org/10.1016/j.coi.2015.06.013. Epub;%2015 Jul 13.:107–12
Soto-Pantoja DR, Kaur S, Roberts DD. CD47 signaling pathways controlling cellular differentiation and responses to stress. Crit Rev Biochem Mol Biol. 2015;50(3):212–30.
Zhang H, Lu H, Xiang L, Bullen JW, Zhang C, Samanta D, et al. HIF-1 regulates CD47 expression in breast cancer cells to promote evasion of phagocytosis and maintenance of cancer stem cells. Proc Natl Acad Sci U S A. 2015;112(45):E6215–E23.
de Kruijf EM, Sajet A, van Nes JG, Natanov R, Putter H, Smit VT, et al. HLA-E and HLA-G expression in classical HLA class I-negative tumors is of prognostic value for clinical outcome of early breast cancer patients. J Immunol. 2010;185(12):7452–9.
Majeti R, Chao MP, Alizadeh AA, Pang WW, Jaiswal S, Gibbs KD, et al. CD47 is an adverse prognostic factor and therapeutic antibody target on human acute myeloid leukemia stem cells. Cell. 2009;138(2):286–99.
Serrano M, Lin AW, McCurrach ME, Beach D, Lowe SW. Oncogenic ras provokes premature cell senescence associated with accumulation of p53 and p16INK4a. Cell. 1997;88(5):593–602.
Coppe JP, Desprez PY, Krtolica A, Campisi J. The senescence-associated secretory phenotype: the dark side of tumor suppression. Annu Rev Pathol. 2010;5:99–118. https://doi.org/10.1146/annurev-pathol-121808-102144.:99-118.
Bakheet T, Hitti E, Khabar Khalid SA. ARED-Plus: an updated and expanded database of AU-rich element-containing mRNAs and pre-mRNAs. Nucleic Acids Res. 2017;46(D1):D218–D20.
Moore AE, Chenette DM, Larkin LC, Schneider RJ. Physiological networks and disease functions of RNA-binding protein AUF1. Wiley Interdiscip Rev RNA. 2014;5(4):549–64.
White EJ, Brewer G, Wilson GM. Post-transcriptional control of gene expression by AUF1: mechanisms, physiological targets, and regulation. Biochim Biophys Acta. 2013;1829(6–7):680–8.
Zinder JC, Lima CD. Targeting RNA for processing or destruction by the eukaryotic RNA exosome and its cofactors. Genes Dev. 2017;31(2):88–100.
Laroia G, Cuesta R, Brewer G, Schneider RJ. Control of mRNA decay by heat shock-ubiquitin-proteasome pathway. Science. 1999;284(5413):499–502.
Hitti E, Iakovleva T, Brook M, Deppenmeier S, Gruber AD, Radzioch D, et al. Mitogen-activated protein kinase-activated protein kinase 2 regulates tumor necrosis factor mRNA stability and translation mainly by altering tristetraprolin expression, stability, and binding to adenine/uridine-rich element. Mol Cell Biol. 2006;26(6):2399–407.
Clement SL, Scheckel C, Stoecklin G, Lykke-Andersen J. Phosphorylation of tristetraprolin by MK2 impairs AU-rich element mRNA decay by preventing deadenylase recruitment. Mol Cell Biol. 2011;31(2):256–66.
Tiedje C, Diaz-Munoz MD, Trulley P, Ahlfors H, Laaß K, Blackshear PJ, et al. The RNA-binding protein TTP is a global post-transcriptional regulator of feedback control in inflammation. Nucleic Acids Res. 2016;44(15):7418–40.
Brooks SA, Blackshear PJ. Tristetraprolin (TTP): interactions with mRNA and proteins, and current thoughts on mechanisms of action. Biochim Biophys Acta. 2013;1829(6–7):666–79.
Shiloh Y, Ziv Y. The ATM protein kinase: regulating the cellular response to genotoxic stress, and more. Nat Rev Mol Cell Biol. 2013;14(4):197–210.
Tee AR, Proud CG. DNA-damaging agents cause inactivation of translational regulators linked to mTOR signalling. Oncogene. 2000;19(26):3021–31.
Jackson SP, Bartek J. The DNA-damage response in human biology and disease. Nature. 2009;461(7267):1071–8.
Holohan C, Van Schaeybroeck S, Longley DB, Johnston PG. Cancer drug resistance: an evolving paradigm. Nat Rev Cancer. 2013;13(10):714–26.
Cannell IG, Merrick KA, Morandell S, Zhu CQ, Braun CJ, Grant RA, et al. A pleiotropic RNA-binding protein controls distinct cell cycle checkpoints to drive resistance of p53-defective tumors to chemotherapy. Cancer Cell. 2015;28(5):623–37.
Lalaoui N, Hanggi K, Brumatti G, Chau D, Nguyen NYN, Vasilikos L, et al. Targeting p38 or MK2 enhances the anti-leukemic activity of Smac-Mimetics. Cancer Cell. 2016;29(2):145–58.
Patnaik A, Haluska P, Tolcher AW, Erlichman C, Papadopoulos KP, Lensing JL, et al. A first-in-human phase I study of the oral p38 MAPK inhibitor, Ralimetinib (LY2228820 dimesylate), in patients with advanced cancer. Clin Cancer Res. 2016;22(5):1095.
Kuma Y, Sabio G, Bain J, Shpiro N, Marquez R, Cuenda A. BIRB796 inhibits all p38 MAPK isoforms in vitro and in vivo. J Biol Chem. 2005;280(20):19472–9.
Grattendick KJ, Nakashima JM, Feng L, Giri SN, Margolin SB. Effects of three anti-TNF-alpha drugs: etanercept, infliximab and pirfenidone on release of TNF-alpha in medium and TNF-alpha associated with the cell in vitro. Int Immunopharmacol. 2008;8(5):679–87.
Nakazato H, Oku H, Yamane S, Tsuruta Y, Suzuki R. A novel anti-fibrotic agent pirfenidone suppresses tumor necrosis factor-α at the translational level. Eur J Pharmacol. 2002;446(1–3):177–85.
Kagoya Y, Yoshimi A, Kataoka K, Nakagawa M, Kumano K, Arai S, et al. Positive feedback between NF-kappaB and TNF-alpha promotes leukemia-initiating cell capacity. J Clin Invest. 2014;124(2):528–42.
Frelin C, Imbert V, Griessinger E, Peyron AC, Rochet N, Philip P, et al. Targeting NF-κB activation via pharmacologic inhibition of IKK2-induced apoptosis of human acute myeloid leukemia cells. Blood. 2005;105(2):804.
Chang TP, Vancurova I. Bcl3 regulates pro-survival and pro-inflammatory gene expression in cutaneous T-cell lymphoma. Biochim Biophys Acta. 2014;1843(11):2620–30.
Haq R, Yokoyama S, Hawryluk EB, Jönsson GB, Frederick DT, McHenry K, et al. BCL2A1 is a lineage-specific antiapoptotic melanoma oncogene that confers resistance to BRAF inhibition. Proc Natl Acad Sci. 2013;110(11):4321–6.
Kurosu T, Fukuda T, Miki T, Miura O. BCL6 overexpression prevents increase in reactive oxygen species and inhibits apoptosis induced by chemotherapeutic reagents in B-cell lymphoma cells. Oncogene. 2003;22(29):4459–68.
Rushworth SA, Bowles KM, Raninga P, MacEwan DJ. NF-kappaB-inhibited acute myeloid leukemia cells are rescued from apoptosis by heme oxygenase-1 induction. Cancer Res. 2010;70(7):2973–83.
Hsu H, Shu HB, Pan MG, Goeddel DV. TRADD-TRAF2 and TRADD-FADD interactions define two distinct TNF receptor 1 signal transduction pathways. Cell. 1996;84(2):299–308.
Annibaldi A, Meier P. Checkpoints in TNF-induced cell death: implications in inflammation and Cancer. Trends Mol Med. 2018;24(1):49–65.
Rauert-Wunderlich H, Siegmund D, Maier E, Giner T, Bargou RC, Wajant H, et al. The IKK inhibitor Bay 11-7082 induces cell death independent from inhibition of activation of NFkappaB transcription factors. PLoS One. 2013;8(3):e59292.
Strickson S, Campbell D, Emmerich C, Knebel A, Plater L, Ritorto MS, et al. The anti-inflammatory drug BAY 11–7082 suppresses the MyD88-dependent signaling network by targeting the ubiquitin system. Biochem J. 2013;451:427–37.
Emmons J, Townley-Tilson WHD, Deleault KM, Skinner SJ, Gross RH, Whitfield ML, et al. Identification of TTP mRNA targets in human dendritic cells reveals TTP as a critical regulator of dendritic cell maturation. RNA. 2008;14(5):888–902.
Kondoh K, Nishida E. Regulation of MAP kinases by MAP kinase phosphatases. Biochim Biophys Acta. 2007;1773(8):1227–37.
Li ZW, Dalton WS. Tumor microenvironment and drug resistance in hematologic malignancies. Blood Rev. 2006;20(6):333–42.
Townsend E-C, Murakami M-A, Christodoulou A, Christie A-Á, Kοster J, Desouza T-A, et al. The public repository of xenografts enables discovery and randomized phase II-like trials in mice. Cancer Cell. 2016;30(1):183.
Sykes DB, Kfoury YS, Mercier FE, Wawer MJ, Law JM, Haynes MK, et al. Inhibition of dihydroorotate dehydrogenase overcomes differentiation blockade in acute myeloid leukemia. Cell. 2016;167(1):171–86.
Sosa MS, Avivar-Valderas A, Bragado P, Wen HC, Aguirre-Ghiso JA. ERK1/2 and p38α/β signaling in tumor cell quiescence: opportunities to control dormant residual disease. Clin Cancer Res. 2011;17(18):5850–7.
Claudio N, Dalet A, Gatti E, Pierre P. Mapping the crossroads of immune activation and cellular stress response pathways. EMBO J. 2013;32(9):1214–24.
Senft D, Ronai ZA. UPR, autophagy, and mitochondria crosstalk underlies the ER stress response. Trends Biochem Sci. 2015;40(3):141–8.
Cuenda A, Rousseau S. p38 MAP-kinases pathway regulation, function and role in human diseases. Biochim Biophys Acta. 2007;1773(8):1358–75.
Korb A, Tohidast-Akrad M, Cetin E, Axmann R, Smolen J, Schett G. Differential tissue expression and activation of p38 MAPK alpha, beta, gamma, and delta isoforms in rheumatoid arthritis. Arthritis Rheum. 2006;54(9):2745–56.
Liu J, Lin A. Role of JNK activation in apoptosis: a double-edged sword. Cell Res. 2005;15:36.
Carballo E, Cao H, Lai WS, Kennington EA, Campbell D, Blackshear PJ. Decreased sensitivity of tristetraprolin-deficient cells to p38 inhibitors suggests the involvement of tristetraprolin in the p38 signaling pathway. J Biol Chem. 2001;276(45):42580–7.
Chrestensen CA, Schroeder MJ, Shabanowitz J, Hunt DF, Pelo JW, Worthington MT, et al. MAPKAP kinase 2 phosphorylates tristetraprolin on in vivo sites including Ser178, a site required for 14-3-3 binding. J Biol Chem. 2004;279(11):10176–84.
Schichl YM, Resch U, Lemberger CE, Stichlberger D, de Martin R. Novel phosphorylation-dependent ubiquitination of tristetraprolin by mitogen-activated protein kinase/extracellular signal-regulated kinase kinase kinase 1 (MEKK1) and tumor necrosis factor receptor-associated factor 2 (TRAF2). J Biol Chem. 2011;286(44):38466–77.
Nichols RC, Botson J, Wang XW, Hamilton BJ, Collins JE, Uribe V, et al. A flexible approach to studying post-transcriptional gene regulation in stably transfected mammalian cells. Mol Biotechnol. 2011;48(3):210–7.
Zhang T, Kruys V, Huez G, Gueydan C. AU-rich element-mediated translational control: complexity and multiple activities of trans-activating factors. Biochem Soc Trans. 2001;30(6):952–8.
Mukherjee N, Corcoran DL, Nusbaum JD, Reid DW, Georgiev S, Hafner M, et al. Integrative regulatory mapping indicates that the RNA-binding protein HuR couples pre-mRNA processing and mRNA stability. Mol Cell. 2011;43(3):327–39.
Lal A, Mazan-Mamczarz K, Kawai T, Yang X, Martindale JL, Gorospe M. Concurrent versus individual binding of HuR and AUF1 to common labile target mRNAs. EMBO J. 2004;23(15):3092–102.
Kesarwani M, Kincaid Z, Gomaa A, Huber E, Rohrabaugh S, Siddiqui Z, et al. Targeting c-FOS and DUSP1 abrogates intrinsic resistance to tyrosine-kinase inhibitor therapy in BCR-ABL-induced leukemia. Nat Med. 2017;23:472.
Coelho MA, de Carne TS, Rana S, Zecchin D, Moore C, Molina-Arcas M, et al. Oncogenic RAS signaling promotes tumor Immunoresistance by stabilizing PD-L1 mRNA. Immunity. 2017;47(6):1083–99.
Zong H, Gozman A, Caldas-Lopes E, Taldone T, Sturgill E, Brennan S, et al. A hyperactive signalosome in acute myeloid leukemia drives addiction to a tumor-specific Hsp90 species. Cell Rep. 2015;13(10):2159–73.
Zuber J, Radtke I, Pardee TS, Zhao Z, Rappaport AR, Luo W, et al. Mouse models of human AML accurately predict chemotherapy response. Genes Dev. 2009;23(7):877–89.
Chen Y, Jacamo R, Konopleva M, Garzon R, Croce C, Andreeff M. CXCR4 downregulation of let-7a drives chemoresistance in acute myeloid leukemia. J Clin Invest. 2013;123(6):2395–407.
Vasudevan S, Steitz JA. AU-rich-element-mediated upregulation of translation by FXR1 and Argonaute 2. Cell. 2007;128(6):1105–18.
Gandin V, Sikstrom K, Alain T, Morita M, McLaughlan S, Larsson O, et al. Polysome fractionation and analysis of mammalian translatomes on a genome-wide scale. J Vis Exp. 2014;87:10.
Truitt ML, Conn CS, Shi Z, Pang X, Tokuyasu T, Coady AM, et al. Differential requirements for eIF4E dose in normal development and cancer. Cell. 2015;162(1):59–71.
Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protocols. 2008;4(1):44–57.
Huang DW, Sherman BT, Lempicki RA. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009;37(1):1–13.
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci. 2005;102(43):15545–50.
Stoecklin G, Stubbs T, Kedersha N, Wax S, Rigby WF, Blackwell TK, et al. MK2-induced tristetraprolin:14-3-3 complexes prevent stress granule association and ARE-mRNA decay. EMBO J. 2004;23(6):1313–24.
Sun L, Stoecklin G, Van WS, Hinkovska-Galcheva V, Guo RF, Anderson P, et al. Tristetraprolin (TTP)-14-3-3 complex formation protects TTP from dephosphorylation by protein phosphatase 2a and stabilizes tumor necrosis factor-alpha mRNA. J Biol Chem. 2007;282(6):3766–77.
Hoy CA, Seamer LC, Schimke RT. Thermal denaturation of DNA for immunochemical staining of incorporated bromodeoxyuridine (BrdUrd): critical factors that affect the amount of fluorescence and the shape of BrdUrd/DNA histogram. Cytometry. 1989;10(6):718–25.
Ting L, Rad R, Gygi SP, Haas W. MS3 eliminates ratio distortion in isobaric labeling-based multiplexed quantitative proteomics. Nat Methods. 2011;8(11):937–40.
Ozes O, Blatt LM, Seiwert SD. Use of pirfenidone in therapeutic regimens. United States Patent-US 7,407,973; 2008. p. 1–46.
Schaefer CJ, Ruhrmund DW, Pan L, Seiwert SD, Kossen K. Antifibrotic activities of pirfenidone in animal models. Eur Respir Rev. 2011;20(120):85.
King TE, Bradford WZ, Castro-Bernardini S, Fagan EA, Glaspole I, Glassberg MK, et al. A phase 3 trial of pirfenidone in patients with idiopathic pulmonary fibrosis. N Engl J Med. 2014;370(22):2083–92.
Kojima K, Konopleva M, Samudio IJ, Shikami M, Cabreira-Hansen M, McQueen T, et al. MDM2 antagonists induce p53-dependent apoptosis in AML: implications for leukemia therapy. Blood. 2005;106(9):3150–9.
Grant S. Ara-C: cellular and molecular pharmacology. Adv Cancer Res. 1998;72:197-233.
Tate C, Blosser W, Wyss L, Evans G, Xue Q, Pan Y, et al. LY2228820 dimesylate, a selective inhibitor of p38 mitogen-activated protein kinase, reduces angiogenic endothelial cord formation in vitro and in vivo. J Biol Chem. 2013;288(9):6743–53.
Campbell RM, Anderson BD, Brooks NA, Brooks HB, Chan EM, De DA, et al. Characterization of LY2228820 dimesylate, a potent and selective inhibitor of p38 MAPK with antitumor activity. Mol Cancer Ther. 2014;13(2):364–74.
Pargellis C, Tong L, Churchill L, Cirillo PF, Gilmore T, Graham AG, et al. Inhibition of p38 MAP kinase by utilizing a novel allosteric binding site. Nat Struct Biol. 2002;9(4):268–72.
Yasui H, Hideshima T, Ikeda H, Jin J, Ocio EM, Kiziltepe T, et al. BIRB 796 enhances cytotoxicity triggered by bortezomib, heat shock protein (Hsp) 90 inhibitor, and dexamethasone via inhibition of p38 mitogen-activated protein kinase/Hsp27 pathway in multiple myeloma cell lines and inhibits paracrine tumour growth. Br J Haematol. 2007;136(3):414–23.
Regan J, Breitfelder S, Cirillo P, Gilmore T, Graham AG, Hickey E, et al. Pyrazole urea-based inhibitors of p38 MAP kinase: from lead compound to clinical candidate. J Med Chem. 2002;45(14):2994–3008.
Regan J, Capolino A, Cirillo PF, Gilmore T, Graham AG, Hickey E, et al. Structure-activity relationships of the p38alpha MAP kinase inhibitor 1-(5-tert-butyl-2-p-tolyl-2H-pyrazol-3-yl)-3-[4-(2-morpholin-4-yl-ethoxy)naph- thalen-1-yl]urea (BIRB 796). J Med Chem. 2003;46(22):4676–86.
Zhang T, Inesta-Vaquera F, Niepel M, Zhang J, Ficarro SB, Machleidt T, et al. Discovery of potent and selective covalent inhibitors of JNK. Chem Biol. 2012;19(1):140–54.
Batey MA, Zhao Y, Kyle S, Richardson C, Slade A, Martin NM, et al. Preclinical evaluation of a novel ATM inhibitor, KU59403, in vitro and in vivo in p53 functional and dysfunctional models of human cancer. Mol Cancer Ther. 2013;12(6):959–67.
Xu F, Wang F, Yang T, Sheng Y, Zhong T, Chen Y. Differential drug resistance acquisition to doxorubicin and paclitaxel in breast cancer cells. Cancer Cell Int. 2014;14(1):538.
Bailey TL, Elkan C. Fitting a mixture model by expectation maximization to discover motifs in biopolymers. Proc Int Conf Intell Syst Mol Biol. 1994;2(28–36):28–36.
Karolchik D, Hinrichs AS, Furey TS, Roskin KM, Sugnet CW, Haussler D, et al. The UCSC table browser data retrieval tool. Nucleic Acids Res. 2004;32(Database issue):D493–D6.
Arava Y, Wang Y, Storey JD, Liu CL, Brown PO, Herschlag D. Genome-wide analysis of mRNA translation profiles in Saccharomyces cerevisiae. Proc Natl Acad Sci. 2003;100(7):3889–94.
Piques M, Schulze WX, Höhne M, Usadel BR, Gibon Y, Rohwer J, et al. Ribosome and transcript copy numbers, polysome occupancy and enzyme dynamics in Arabidopsis. Mol Syst Biol. 2009;5:314.
Liu MJ, Wu SH, Chen HM, Wu SH. Widespread translational control contributes to the regulation of Arabidopsis photomorphogenesis. Mol Syst Biol. 2012;8:566.
Spasic M, Friedel CC, Schott J, Kreth J, Leppek K, Hofmann S, et al. Genome-wide assessment of AU-rich elements by the AREScore algorithm. PLoS Genet. 2012;8(1):e1002433.
Cook KB, Kazan H, Zuberi K, Morris Q, Hughes TR. RBPDB: a database of RNA-binding specificities. Nucleic Acids Res. 2011;39(Database issue):D301–D8.
Mootha VK, Lindgren CM, Eriksson K-F, Subramanian A, Sihag S, Lehar J, et al. PGC-1α-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes. Nat Genet. 2003;34(3):267–73.
Gautier L, Cope L, Bolstad BM, Irizarry RA. affy--analysis of Affymetrix GeneChip data at the probe level. Bioinformatics (Oxford, England). 2004;20(3):307–15.
Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43(7):e47.
Lee S, Micalizzi D, Truesdell SS, Bukhari SI, Boukhali M, Lombardi-Story J, Kato Y, Choo M, Dey-Guha I, Ji F, Nicholson BT, Myers DT, Lee D, Mazzola MA, Raheja R, Langenbucher A, Haradhvala NJ, Lawrence MS, Gandhi R, Tiedje C, Diaz-Munoz M, Sweetser DA, Sadreyev R, Sykes D, Haas W, Haber DA, Maheswaran S, Vasudevan S. A post-transcriptional program of chemoresistance by AU-rich elements and TTP in quiescent leukemic cells. Datasets. Gene Expression Omnibus. (2019). https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE141332.
We thank Partners Healthcare Center for Personalized Genetic Medicine, & Microarray & Sequencing Resource, BUMC facilities for microarray data; N. Kedersha, S. Lyons, and P. Anderson for plasmids and antibody; T. Graubert for patient samples; D. Bloch, S. Wu, A. Naar, M. Gaestel, M. Guzman, N. Bardeesy, D. Scadden, S. Ramaswamy, and C. Benes for reagents.
Peer review information
Anahita Bishop was the primary editor on this article and managed its editorial process and peer review in collaboration with the rest of the editorial team.
The study is funded by Cancer Research Institute, D. & M-E Ryder, Leukemia & Lymphoma Society, Smith Family Foundation, GM100202 grants to SV & CA185086. SL was funded by Fund for Medical Discovery fellowship and by a postdoctoral fellowship from the Basic Science Research Program through the National Research Foundation of Korea (2015017218). DAS is funded by NCI CA115772. DAH is funded by Howard Hughes Medical Institute & 2R01 CA129933.
Ethics approval and consent to participate
Statement on human data: All human samples (de-identified) were handled in accordance with Institutional Review Board (IRB) protocols to SV (2015P000998/MGH), approved by the Partners Human Research Committee Institutional Review Board/MGH IRB. Samples were obtained by David A. Sweetser, and Tim Graubert and Jennifer Lombardi-Story of the Hematologic Malignancies tissue bank oversight committee (protocol DF/HCC 13-583), approved by DF/HCC Office for Human Research Studies, and include written consent for publication and use of patient samples from all patients. All experimental methods used comply with the Helsinki Declaration.
Statement on animal data: nod-scid-gamma (NSG), C57Black/6 mice, 10–12 weeks, are obtained from MGH Cox-7 Gnotobiotic animal facility of the AAALAC-accredited Center for Comparative Medicine and Services at MGH. These facilities are supervised by veterinarians in the Center for Comparative Medicine and the MGH Subcommittee for Animal Research (SRAC) and maintained according to the protocol approved by SRAC, and provide services for breeding, regular health checks, histopathology, and macropathology.
Consent for publication
All samples are de-identified, described under the above IRB protocols. The protocols include written consent for publication and use of patient samples from all patients. All experimental methods used comply with the Helsinki Declaration.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Uncropped Western blots.
About this article
Cite this article
Lee, S., Micalizzi, D., Truesdell, S.S. et al. A post-transcriptional program of chemoresistance by AU-rich elements and TTP in quiescent leukemic cells. Genome Biol 21, 33 (2020). https://doi.org/10.1186/s13059-020-1936-4
- Post-transcriptional regulation
- AU-rich elements