The mTORC1-mediated activation of ATF4 promotes protein and glutathione synthesis downstream of growth signals
Abstract
The mechanistic target of rapamycin complex 1 (mTORC1) stimulates a coordinated anabolic program in response to growth-promoting signals. Paradoxically, recent studies indicate that mTORC1 can activate the transcription factor ATF4 through mechanisms distinct from its canonical induction by the integrated stress response (ISR). However, its broader roles as a downstream target of mTORC1 are unknown. Therefore, we directly compared ATF4-dependent transcriptional changes induced upon insulin-stimulated mTORC1 signaling to those activated by the ISR. In multiple mouse embryo fibroblast and human cancer cell lines, the mTORC1-ATF4 pathway stimulated expression of only a subset of the ATF4 target genes induced by the ISR, including genes involved in amino acid uptake, synthesis, and tRNA charging. We demonstrate that ATF4 is a metabolic effector of mTORC1 involved in both its established role in promoting protein synthesis and in a previously unappreciated function for mTORC1 in stimulating cellular cystine uptake and glutathione synthesis.
eLife digest
When building healthy tissue, the human body must carefully control the growth of new cells to prevent them from becoming cancerous. A core component of this regulation is the protein mTORC1, which responds to various growth-stimulating factors and nutrients, and activates the chemical reactions cells need to grow. Part of this process involves controlling ‘nutrient-sensing transcription factors’ – proteins that regulate the activity of specific genes based on the availability of different nutrients.
One of these nutrient-sensing transcription factors, ATF4, has recently been shown to be involved in some of the processes triggered by mTORC1. The role this factor plays in how cells respond to stress – such as when specific nutrients are depleted, protein folding is disrupted or toxins are present – is well-studied. But how it reacts to the activation of mTORC1 is less clear. To bridge this gap, Torrence et al. studied mouse embryonic cells and human prostate cancer cells grown in the laboratory, to see whether mTORC1 influenced the behavior of ATF4 differently than cellular stress.
Cells were treated either with insulin, which activates mTORC1, or an antibiotic that sparks the stress response. The cells were then analyzed using a molecular tool to see which genes were switched on by ATF4 following treatment. This revealed that less than 10% of the genes activated by ATF4 during cellular stress are also activated in response to mTORC1-driven growth.
Many of the genes activated in both scenarios were involved in synthesizing and preparing the building blocks that make up proteins. This was consistent with the discovery that ATF4 helps mTORC1 stimulate growth by promoting protein synthesis. Torrence et al. also found that mTORC1’s regulation of ATF4 stimulated the synthesis of glutathione, the most abundant antioxidant in cells.
The central role mTORC1 plays in controlling cell growth means it is important to understand how it works and how it can lead to uncontrolled growth in human diseases. mTORC1 is activated in many overgrowth syndromes and the majority of human cancers. These new findings could provide insight into how tumors coordinate their drive for growth while adapting to cellular stress, and reveal new drug targets for cancer treatment.
Introduction
Pro-growth signals in the form of growth factors, hormones, and nutrients impinge on cellular metabolic programs in a coordinated fashion involving both acute, post-translational regulation and transcriptional control of nutrient transporters and metabolic enzymes. The mechanistic target of rapamycin complex 1 (mTORC1) acts as a central point of integration for these signals and propagates a downstream metabolic response that increases anabolic processes while decreasing specific catabolic processes. Through a variety of downstream effectors, mTORC1 stimulates the synthesis of the major macromolecules comprising cellular biomass, including protein, lipid, and nucleic acids, along with metabolic and adaptive pathways that support this anabolic program (Valvezan and Manning, 2019).
A particularly interesting feature of this coordinated metabolic program downstream of mTORC1 is the co-opting of key nutrient-sensing transcription factors that are established to be activated, independent of mTORC1, in response to depletion of specific nutrients (Torrence and Manning, 2018). In their canonical roles, these transcription factors serve to mount an adaptive response by upregulating genes that allow cells to overcome the specific nutrient deficiency. Perhaps the best characterized of these transcription factors with dual regulation is the hypoxia-inducible factor 1 (HIF1) comprising the labile HIF1α protein heterodimerized with the aryl hydrocarbon receptor nuclear translocator (ARNT or HIF1β). Oxygen depletion (i.e., hypoxia) results in the rapid stabilization of HIF1α and allows the HIF1 heterodimer to induce genes involved in glucose uptake, glycolysis, and angiogenesis to adapt to hypoxia and decrease mitochondrial respiration (Nakazawa et al., 2016). On the other hand, in response to upstream growth factor signaling pathways, activation of mTORC1 stimulates an increase in HIF1α protein synthesis, leading to elevated expression of HIF1 gene targets (Brugarolas et al., 2003; Düvel et al., 2010; Hudson et al., 2002; Laughner et al., 2001; Thomas et al., 2006; Zhong et al., 2000). The result is an mTORC1-mediated increase in glucose uptake and glycolysis even when oxygen is not limiting (e.g., normoxia), a process referred to as aerobic glycolysis, which can support the production of biosynthetic precursors in the form of glycolytic intermediates. Similarly, the sterol regulatory element (SRE)-binding protein (SREBP) family of transcription factors are independently regulated by both adaptive nutrient signals and growth signals controlling mTORC1. The SREBPs are canonically activated upon sterol depletion and induce expression of the enzymes required for de novo synthesis of fatty acid and sterol lipids (Horton et al., 2002). However, insulin and growth factor signaling can also induce lipid synthesis via mTORC1-stimulated activation of SREBP and its lipogenic gene targets (Düvel et al., 2010; Owen et al., 2012; Peterson et al., 2011; Porstmann et al., 2008). Recent studies have suggested that the regulation of nutrient-sensing transcription factors by mTORC1 signaling extends to the activating transcription factor 4 (ATF4) (Adams, 2007; Ben-Sahra et al., 2016; Torrence and Manning, 2018).
ATF4 is a basic leucine zipper (bZIP) transcription factor that is selectively translated in response to specific forms of cellular stress to induce the expression of genes involved in adaptation to stress (Walter and Ron, 2011). This adaptive program is referred to as the integrated stress response (ISR) and is initiated by stress-activated protein kinases, including general control nonderepressible 2 (GCN2) activated upon amino acid deprivation and protein kinase RNA-like endoplasmic reticulum kinase (PERK) activated by ER stress, among others, which phosphorylate eIF2α on Ser51 (Harding et al., 2000). Phosphorylation of eIF2α serves to globally attenuate mRNA translation to conserve amino acids and energy and decrease the cellular protein load as one adaptive measure to overcome these stresses (Clemens, 1996). Importantly, a small number of mRNAs, including that encoding ATF4, exhibit increased translation upon eIF2α-Ser51 phosphorylation (Vattem and Wek, 2004). The stress-induced increase in ATF4 leads to the expression of a canonical set of ATF4 target genes, including those involved in nonessential amino acid (NEAA) biosynthesis and amino acid transport, as part of the adaptive cellular response specific to stresses such as amino acid depletion (Harding et al., 2003).
ATF4 functions in heterodimers with other bZIP transcription factors and also co-regulates many of its target genes with additional transcription factors as part of the cellular stress response (Newman and Keating, 2003; Wortel et al., 2017). However, whether and how the many distinct upstream stresses that activate ATF4 influence its heterodimerization partners and the induction of specific sets of genes is not well understood.
While ATF4 is a major downstream effector of the ISR, evidence has emerged that ATF4 can also be activated by pro-growth signals that stimulate mTORC1 signaling (Adams, 2007; Ben-Sahra et al., 2016), and cis-regulatory elements for ATF4 binding are enriched in the promoters of mTORC1-induced genes (Düvel et al., 2010). Importantly, the mTORC1-mediated activation of ATF4 involves its increased translation in a manner that is independent of the ISR and phosphorylation of eIF2α (Ben-Sahra et al., 2016; Park et al., 2017). These findings suggest that, similar to HIF1 and SREBP, ATF4 induction may be mobilized as part of the broader anabolic program downstream of mTORC1. Indeed, our previous findings indicate that mTORC1 promotes de novo purine synthesis, in part, through induction of mitochondrial one-carbon metabolism via ATF4 activation and expression of its gene target MTHFD2 (Ben-Sahra et al., 2016).
How the ATF4-dependent gene program compares between its adaptive role in the ISR and its activation as a downstream effector of mTORC1 signaling and whether ATF4 contributes to established or new functions of mTORC1 are unknown (Figure 1A). Here, we find that the mTORC1-ATF4 program represents a small subset of ATF4-dependent genes induced by ER stress and includes genes encoding the enzymes required for tRNA charging, NEAA synthesis, and amino acid uptake. Consistent with regulation of these enzymes by mTORC1 through ATF4, ATF4 contributes to the induction of protein synthesis downstream of mTORC1. We also find that mTORC1 signaling promotes glutathione synthesis through ATF4 and its specific regulation of the cystine transporter SLC7A11. Thus, ATF4 is an anabolic effector of mTORC1 signaling, necessary for both its canonical regulation of protein synthesis and its induced synthesis of glutathione, the most abundant antioxidant in cells.
Results
mTORC1 signaling activates a subset of ATF4-dependent genes also activated by the ISR
To identify the ATF4-dependent gene targets downstream of mTORC1, we compared the insulin-induced, rapamycin-sensitive transcripts between wild-type MEFs and those with biallelic loss of ATF4 via CRISPR/Cas9 gene deletion (see Materials and methods). Consistent with our previous studies (Ben-Sahra et al., 2016), insulin stimulated an increase in ATF4 protein in MEFs, which was decreased with rapamycin (Figure 1B, Figure 1—figure supplement 1). In parallel, these cells were treated with a time course of tunicamycin, an inhibitor of N-glycosylation that potently induces ER stress and an increase in ATF4, to identify ATF4 gene targets downstream of the ISR. RNA-seq analysis revealed that 20% of transcripts (253 total) significantly upregulated upon insulin stimulation were significantly blocked in their induction with rapamycin treatment. Approximately 30% of these mTORC1-regulated genes (77 total) lost their insulin responsiveness with ATF4 deletion. In comparison, 36% of transcripts significantly induced with tunicamycin treatment at 4 hr were dependent on ATF4 (774 total). Importantly, the expression of just 61 genes was found to overlap between these two modes of ATF4 regulation, being ATF4 dependent in response to both mTORC1 activation and the ISR (Figure 1C, Figure 1—source data 1).
The RNA-seq analysis demonstrated that only 8% of ATF4 gene targets induced by ER stress were also significantly stimulated by mTORC1 signaling (e.g., insulin induced and rapamycin sensitive). Interestingly, these 61 shared genes showed significant KEGG pathway enrichment for aminoacyl tRNA biosynthesis, amino acid metabolism, and one-carbon metabolism (Figure 1D). The genes shared between the ISR and mTORC1 signaling were greatly enriched among those exhibiting the most significant increase upon tunicamycin treatment, with 75% (46 genes) lying within the top 100 of the 774 tunicamycin-induced genes (Figure 1E). It is worth noting that many of the top ATF4-dependent genes that scored as being induced by the ISR alone showed some degree of rapamycin-sensitive induction with insulin but did not reach statistical significance in the RNA-seq analyses. These data indicate that the subset of ATF4-dependent genes induced by mTORC1 signaling largely comprised those that are also most sensitive to ATF4 induction by the ISR. Among the 61 shared ATF4-dependent transcripts, those involved in amino acid synthesis and transport, one-carbon metabolism, and aminoacyl tRNA charging often displayed comparable fold changes between insulin stimulation and tunicamycin treatment, while canonical genes of the ER stress response, such as Herpud1 and Ddit3/Chop, showed much greater induction with tunicamycin (Figure 1F). In addition to aminoacyl tRNA synthetase genes, expression of the Xpot gene encoding Exportin-T, which is the major Ran GTPase family member for nuclear to cytosolic export of mature tRNAs (Arts et al., 1998; Kutay et al., 1998), was found to be similarly regulated by ATF4 in response to mTORC1 and ISR activation. Transcripts encoding known negative regulators of mTORC1 signaling are also among these 61 shared ATF4-induced genes, including Ddit4/Redd1 and Sesn2 (Brugarolas et al., 2004; Condon et al., 2021; Lee et al., 2010; Reiling and Hafen, 2004; Wolfson et al., 2016). These targets likely contribute to the ATF4-dependent inhibition of mTORC1 signaling (S6K1 phosphorylation) observed upon tunicamycin treatment (Figure 1B), while in the context of mTORC1 signaling, these ATF4 targets might play a role in negative feedback regulation of mTORC1. These data suggest that a specific subset of ATF4-dependent ISR-induced genes are likewise regulated by growth factor signaling through mTORC1 and are enriched for specific processes including aminoacyl tRNA synthesis, amino acid synthesis and uptake, and one-carbon metabolism.
ATF4 is known to form heterodimers with other bZIP transcription factors to engage its gene targets, while also co-regulating genes with transcription factors that bind additional promoter elements (Kilberg et al., 2009; Newman and Keating, 2003; Wortel et al., 2017). Thus, we used bioinformatic tools to determine whether the promoters of ATF4 gene targets shared between the ISR and mTORC1 signaling might be distinct from those induced by the ISR alone. Indeed, CiiiDER analysis (Gearing et al., 2019) revealed that there are predicted promoter-binding sequences that distinguish the 61 shared target genes from the top 200 ATF4-dependent genes induced by the ISR alone (Figure 2A, Figure 2—source data 1). Regulatory elements for the C/EBP family of transcription factors, which are well established to heterodimerize with ATF4 to induce its canonical downstream targets (Cohen et al., 2015; Ebert et al., 2020; Huggins et al., 2015), were the most enriched in the promoters of ATF4 gene targets with shared regulation. On the other hand, binding elements for the TEAD family of transcription factors, which function with YAP/TAZ in the Hippo signaling pathway, were enriched in the promoters of ATF4-dependent gene targets significantly induced only by tunicamycin, consistent with published work indicating a functional connection between the unfolded protein response and YAP-TEAD activation (Wu et al., 2015; Xiao et al., 2019). We next analyzed the 61 ATF4-dependent genes with shared regulation for physical evidence of promoter binding of specific transcription factors using the Cistrome Data Browser, a portal for mining existing chromatin immunoprecipitation-DNA sequencing (ChIP-seq) data (Mei et al., 2017; Zheng et al., 2019). Importantly, this second unbiased analysis also revealed C/EBP isoforms as most commonly binding to the promoters of these genes (Figure 2B).
As all members of the C/EBP family have the potential to heterodimerize with ATF4 and contribute to the induction of these gene targets (Newman and Keating, 2003), we first determined the effects of siRNA-mediated knockdown of individual isoforms, relative to ATF4 knockdown, on expression of three representative genes in Tsc2-/- MEFs, which exhibit growth factor-independent activation of mTORC1 signaling. This analysis revealed that knockdown of ATF4, C/EBPβ, C/EBPδ, or C/EBPγ each led to decreased transcript levels of the shared mTORC1 and ISR gene targets Slc7a5, Mthfd2, and Aars (Figure 2C). However, this analysis was complicated by the finding of substantial co-dependence for expression among these bZIP transcription factors, with knockdown of any one of the C/EBP family members or ATF4 significantly changing expression of at least one other family member. C/EBPδ knockdown, for instance, decreased expression of all genes measured, including ATF4, which was also reflected in loss of ATF4 protein (Figure 2C, D). It is worth noting that we were unable to identify reliable antibodies to specific C/EBP family members for use in MEFs. Among C/EBP family members, C/EBPγ has been found in other settings to regulate many of the genes revealed in our analysis to be induced through shared regulation of ATF4 (Huggins et al., 2015), and its knockdown significantly decreased expression of the three ATF4 target genes tested without effects on ATF4 protein levels (Figure 2C, D). Based on this finding, we knocked down ATF4 or C/EBPγ in wild-type MEFs and stimulated the cells with insulin in the presence or absence of rapamycin to determine whether C/EBPγ impacted the mTORC1 and ATF4-dependent regulation of these genes. Indeed, knockdown of C/EBPγ attenuated the insulin-induced expression of these genes, albeit to a lesser extent than ATF4 knockdown (Figure 2E). C/EBPγ knockdown also blocked the ability of insulin to increase ATF4 transcript levels, suggesting that following the induction of ATF4 mRNA translation downstream of mTORC1 (Ben-Sahra et al., 2016; Park et al., 2017), it stimulates its own expression via ATF4-C/EBPγ heterodimers. Thus, C/EBPγ and likely other ATF4-binding partners of the C/EBP family contribute to the induction of ATF4 gene targets following mTORC1-mediated activation of ATF4.
mTORC1 signaling induces genes involved in amino acid synthesis, transport, and tRNA charging through ATF4 activation
To validate and expand the findings from the RNA-seq analysis, a NanoString codeset was designed to simultaneously quantify transcripts of genes involved in the enriched processes above (see Materials and methods). As positive and negative controls, respectively, we included the glycolytic targets of HIF1, established previously to be regulated downstream of mTORC1 (Düvel et al., 2010), and the mitochondrial tRNA synthetases, not believed to be regulated by ATF4. Using this codeset, we analyzed gene expression in three settings of mTORC1 activation: (1) wild-type MEFs stimulated with insulin in the presence or absence of rapamycin, (2) growth factor-independent activation of mTORC1 via genetic loss of the TSC protein complex in Tsc2-/- MEFs, and (3) Tsc2-/- MEFs with siRNA-mediated knockdown of ATF4 (Figure 3A, Figure 3—source data 1). ATF4 protein levels were robustly upregulated with either genetic or insulin-stimulated mTORC1 activation in these settings, with both rapamycin and ATF4-targeting siRNAs blocking this induction (Figure 3—figure supplement 1A). Interestingly, the majority of transcripts analyzed in the functional categories that encode enzymes of NEAA synthesis, one-carbon metabolism, amino acid transporters, and cytosolic aminoacyl tRNA synthetases (and Xpot) were increased with mTORC1 activation in a manner sensitive to both rapamycin and siRNA-mediated knockdown of ATF4. However, the HIF1 targets of glycolysis were mTORC1-regulated but independent of ATF4, and transcripts encoding the mitochondrial tRNA synthetases were reproducibly regulated by neither mTORC1- nor ATF4. We further confirmed the mTORC1- and ATF4-mediated regulation of a representative subset of these transcripts via qPCR (Figure 3B, C, Figure 3—figure supplement 1B). Consistent with previous studies (Ben-Sahra et al., 2016; Park et al., 2017), both ATF4 transcript and protein levels were induced by mTORC1 signaling in these settings (Figure 3B–E). Transcriptional changes in ATF4 gene targets were reflected in corresponding changes in the abundance of representative protein products, with varying degrees of rapamycin sensitivity, likely reflecting inherent differences in the turnover rates of these proteins (Figure 3D, E, Figure 3—figure supplement 1C, D). We next wanted to confirm that these mTORC1- and ATF4-induced changes were independent of the ISR. While we have shown previously that mTORC1 regulates ATF4 in a manner that is independent of eIF2α-S51 phosphorylation (Ben-Sahra et al., 2016), chronic activation of mTORC1 upon loss of TSC2 is known to cause a basal increase in ER stress and activation of the ISR (Ozcan et al., 2008). Therefore, we utilized MEFs with endogenous, homozygous knock-in of the Eif2a-S51A mutation (Eif2aA/A), which fail to induce ATF4 downstream of cellular stress (Figure 3—figure supplement 1E; Scheuner et al., 2001). Consistent with mTORC1-dependent, ISR-independent regulation, insulin increased ATF4 protein levels and expression of its gene targets involved in amino acid synthesis and one-carbon metabolism, amino acid transporters, and aminoacyl tRNA synthetases in a rapamycin-sensitive manner in these cells, as shown by NanoString analysis and confirmed for a subset of genes by qPCR (Figure 3F, G, Figure 3—figure supplement 1F). Notably, like insulin stimulation, genetic activation of mTORC1 via siRNA-mediated knockdown of Tsc2 in the eIF2αA/A MEFs also increased ATF4 protein levels, further confirming that this regulation can occur independent of the ISR (Figure 3H). Growth factor-independent activation of mTORC1 also occurs upon loss of the PTEN tumor suppressor, and rapamycin was found to decrease ATF4 protein levels and ATF4-dependent expression of representative gene targets in established PTEN-deficient prostate cancer cells, LNCaP and PC3 (Figure 3I–L, Figure 3—figure supplement 1G–J).
The transcription factor c-MYC can be activated downstream of mTORC1 in some settings and has previously been shown to regulate many of the same ATF4 target genes encoding the enzymes of amino acid synthesis and transport, as well as aminoacyl tRNA synthetases (Csibi et al., 2014; Stine et al., 2015; Wall et al., 2008; Zirin et al., 2019). A survey of expression for 11 of such genes in Tsc2-/- cells finds that the majority are not significantly affected by siRNA-mediated knockdown of c-MYC, whereas a few others are modestly but significantly decreased (e.g., Gars, Slc1a5, Psat1), albeit to a lesser degree than with siRNA-mediated knockdown of ATF4 (Figure 3—figure supplement 1K, L).
To determine whether ATF4 activation is both necessary and sufficient for mTORC1-mediated regulation of these gene targets related to amino acid acquisition and utilization, we knocked out Atf4 using CRISPR/Cas9 in Tsc2-/- MEFs and confirmed biallelic disruption (Figure 4A). Protein levels of identified ATF4 targets were decreased in Tsc2-/- Atf4-/- MEFs and fully rescued with expression of wild-type ATF4 but not a DNA-binding domain mutant (ATF4DBD) (Figure 4B). As mTORC1 regulates ATF4 translation through a mechanism requiring its 5′-UTR (Ben-Sahra et al., 2016; Park et al., 2017), the stably rescued cell lines, which express an ATF4 cDNA lacking the 5′-UTR, exhibit protein expression of ATF4 and its encoded gene targets that are resistant to rapamycin (Figure 4C, Figure 4—figure supplement 1). The expression of select ATF4 target transcripts was markedly decreased in Tsc2-/- Atf4-/- MEFs to a similar extent to that measured in Atf4 wild-type cells treated with rapamycin (Figure 4D). These transcript levels were rescued with the expression of wild-type ATF4, but not ATF4DBD, in a manner that was completely or partially rapamycin resistant (Figure 4D). These collective data show that mTORC1 signaling drives the expression of genes involved in tRNA export and charging, amino acid uptake, and NEAA synthesis through its downstream regulation of the ATF4 transcription factor.
Activation of ATF4 contributes to the induction of protein synthesis downstream of mTORC1
mTORC1 induces protein synthesis through multiple downstream targets (Valvezan and Manning, 2019). Given that the major mTORC1-regulated ATF4 target genes identified above are involved in amino acid uptake, synthesis, and tRNA charging, we hypothesized that ATF4 induction through mTORC1 signaling might contribute to the canonical increase in protein synthesis upon mTORC1 activation. Relative rates of protein synthesis were measured via [35S]-methionine incorporation into newly synthesized proteins in Tsc2+/+ and Tsc2-/- cells treated with control siRNAs or Tsc2-/- cells treated with siRNAs targeting ATF4 or Rheb, the small GTPase target of TSC2 that is an essential upstream activator of mTORC1. siRNA-mediated knockdown of either ATF4 or Rheb substantially decreased ATF4 protein levels in Tsc2-/- cells (Figure 5A). Importantly, the elevated rate of protein synthesis in Tsc2-/- MEFs was decreased with ATF4 knockdown to a similar extent to that observed with Rheb knockdown (Figure 5A, B). No change in mTORC1 signaling or phosphorylation of eIF2α was observed with ATF4 knockdown in this setting (Figure 5—figure supplement 1A). Protein synthesis was also measured in the Tsc2-/- Atf4-/- cell lines described above. Notably, the cells lacking ATF4 exhibited increased uptake of [35S]-methionine relative to parental cells and those reconstituted with ATF4 (Figure 5—figure supplement 1B). Despite this unexplained difference in methionine uptake, Atf4 knockout cells exhibited a reduced rate of protein synthesis that was similar to the parental lines treated with rapamycin (Figure 5C, D). However, rapamycin treatment further reduced protein synthesis in the Atf4 knockout cells. Importantly, cells reconstituted with the rapamycin-resistant ATF4 cDNA exhibited rescued protein synthesis, surpassing that observed in cells with endogenous ATF4, but this enhanced protein synthesis was still significantly reduced with rapamycin. Like Tsc2-/- cells, wild-type cells cultured in the presence of growth factors exhibited reduced protein synthesis upon deletion of Atf4, again with a reduction similar to that from rapamycin treatment (Figure 5—figure supplement 1C, D). Together, these data suggest that ATF4 induction downstream of mTORC1 is necessary but not sufficient for mTORC1-regulated protein synthesis, consistent with the multiple mechanisms through which mTORC1 controls this key anabolic process.
mTORC1 regulates cystine uptake through ATF4
Among the 61 shared mTORC1- and ISR-induced ATF4 gene targets identified, the cystine-glutamate antiporter Slc7a11 was the gene with the highest fold induction by RNA-seq analysis upon insulin treatment (Figure 1F). SLC7A11 (also known as xCT) associates with SLC3A2 (CD98) at the plasma membrane and serves as the primary transporter of cystine, the oxidized form of cysteine and predominant cysteine species in both plasma and cell culture media, whereas reduced cysteine is transported through neutral amino acid systems (Figure 6A; Bannai and Kitamura, 1980; Conrad and Sato, 2012). Transcript levels of Slc7a11 were sensitive to rapamycin in Tsc2-/- MEFs and greatly decreased with ATF4 knockout (Figure 6B). Upon reconstitution with rapamycin-resistant ATF4, expression of Slc7a11 was rescued and no longer sensitive to rapamycin, while the ATF4DBD mutant was unable to restore Slc7a11 transcript levels. A similar pattern of mTORC1- and ATF4-regulated expression was measured for Slc3a2 (Figure 6C). SLC7A11 protein, detected using an antibody validated with siRNA knockdown (Figure 6—figure supplement 1A), decreased in Tsc2-/- MEFs treated with mTOR inhibitors and were increased in wild-type MEFs stimulated with insulin in an mTOR-dependent manner (Figure 6D, E, Figure 6—figure supplement 1B, C). SLC7A11 transcript levels were also decreased with both ATF4 knockdown and rapamycin in the PTEN-deficient cancer cell lines LNCaP and PC3, although SLC7A11 expression was relatively more resistant to rapamycin in PC3 cells (Figure 6—figure supplement 1D, E). SLC7A11 protein levels likewise decreased in LNCaP and PC3 cells treated with mTOR inhibitors, without significant changes to the SLC3A2 gene product CD98 (Figure 6F, Figure 6—figure supplement 1F–H). These data confirm and extend the findings from RNA-seq and NanoString analyses (Figure 1F, Figure 3A, F) and demonstrate that ATF4 is both necessary and sufficient for the mTORC1-mediated induction of Slc7a11 expression.
ATF4 is known to be important for the uptake and synthesis of NEAAs (Harding et al., 2003). In agreement with this, we observed that Tsc2-/- Atf4-/- cells fail to proliferate in Dulbecco’s Modified Eagle’s Medium (DMEM), which only contains a subset of NEAAs, while addback of wild-type ATF4, but not ATF4DBD, restored proliferation (Figure 6G). Supplementation of DMEM with a mixture of all NEAAs, including cysteine, allowed the Tsc2-/- Atf4-/- MEFs to proliferate at the same rate as the ATF4-reconstituted cells, while NEAAs lacking cysteine completely failed to support proliferation of these cells (Figure 6H). Furthermore, supplementation with excess reduced cysteine alone, but not equimolar concentrations of oxidized cysteine in the form of cystine, was able to significantly increase proliferation of the Tsc2-/- Atf4-/- MEFs, albeit to a lesser extent than NEAAs plus cysteine. The majority of these cells die after 72 hr in DMEM, and exogenous expression of either ATF4 or SLC7A11 restores their survival (Figure 6I). Taken together, these data indicate that a defect in the acquisition of cysteine, which normally occurs through SLC7A11-mediated uptake of cystine, underlies the inability of Tsc2-/- Atf4-/- cells to proliferate or survive in DMEM and suggest a key role for mTORC1 signaling in controlling cystine uptake through ATF4.
To directly test whether mTORC1 influences cystine uptake, we employed both genetic (Tsc2 loss) and physiological (insulin stimulation) activation of mTORC1, measuring [14C]-cystine uptake in the presence or absence of rapamycin or the xCT inhibitor erastin (Figure 6A; Yang and Stockwell, 2008). Both rapamycin and erastin significantly decreased [14C]-cystine uptake into Tsc2-/- MEFs (Figure 6J). In wild-type MEFs, insulin stimulated an increase in cystine uptake that was inhibited with rapamycin, and this mTORC1-regulated cystine transport was completely lost with ATF4 knockout, reduced to levels of erastin-treated cells (Figure 6K). Additionally, Tsc2-/- Atf4-/- MEFs showed a decrease in cystine uptake when compared to parental Tsc2-/- MEFs, which could be rescued with re-expression of ATF4 (Figure 6L). However, cystine consumption in cells reconstituted with rapamycin-resistant ATF4 was still significantly sensitive to rapamycin treatment, suggesting the existence of additional, ATF4-independent mechanisms influencing the transport or cellular incorporation of cystine downstream of mTORC1. As mTORC2 has been previously suggested to directly regulate xCT (Gu et al., 2017), we utilized Rictor-/- MEFs, which lack mTORC2 activity, to determine whether mTORC2 was contributing to the decreased cystine uptake observed upon treatment with mTOR inhibitors. While Rictor-/- MEFs displayed increased uptake of cystine relative to their wild-type counterparts, cystine uptake, as well as ATF4 protein levels, was sensitive to rapamycin and Torin1 in both cell lines (Figure 6—figure supplement 1I, J). Thus, mTORC1 promotes cellular cystine uptake, at least in part, through the activation of ATF4 and induction of Slc7a11 expression, which supports cell proliferation and survival.
mTORC1 regulates glutathione levels through ATF4-mediated induction of Slc7a11
Cysteine, generally acquired through cystine uptake and reduction, is an essential component of the tripeptide glutathione (Figure 6A), the most abundant antioxidant in cells (Meister, 1983). We hypothesized that the regulation of cystine uptake through mTORC1 and ATF4 might influence cellular glutathione content. Indeed, mTORC1 inhibition with rapamycin or Torin1 significantly decreased total glutathione levels in Tsc2-/- MEFs, albeit less than buthionine-sulfoximine (BSO), a direct inhibitor of glutathione synthesis (Figure 7A; Griffith and Meister, 1979). Similar to BSO treatment, mTOR inhibitors decreased both reduced (GSH) and oxidized (GSSG) forms of glutathione to the same degree, indicating effects on total glutathione abundance rather than its redox state (Figure 7—figure supplement 1A). Stable reconstitution of Tsc2-/- MEFs with TSC2 also decreased total glutathione levels (Figure 7B). To examine this response in vivo, we employed a mouse model of tuberous sclerosis complex involving xenograft tumors derived from the rat TSC2-/- tumor cell line ELT3 (Hodges et al., 2002). To avoid major differences in tumor size from the treatments, we treated tumor-bearing mice for just 5 days with either vehicle or rapamycin, prior to harvesting tumors for immunoblot analysis and metabolite profiling. Importantly, we found that rapamycin treatment strongly decreased ATF4 protein levels in these tumors with a concomitant decrease in glutathione levels, measured by LC-MS in tumor metabolite extracts (Figure 7C, D). An analysis of published metabolomics data (Tang et al., 2019) also revealed that rapamycin treatment significantly decreased glutathione levels in human TSC2-deficient angiomyolipoma cells (Figure 7—figure supplement 1B). Likewise, inhibition of mTORC1 signaling with rapamycin or Torin1 in LNCaP and PC3 cells resulted in a significant decrease in total glutathione levels, although the degree of decrease varied between the two cell lines (Figure 7E), perhaps reflecting the above finding that SLC7A11 expression is more resistant to mTOR inhibitors in PC3 cells (Figure 6—figure supplement 1D–H).
To determine the role of ATF4 and SLC7A11-dependent cystine uptake in glutathione synthesis downstream of mTORC1 signaling, we compared Tsc2-/- MEFs with or without Atf4 knockout. Total glutathione levels were greatly decreased in Tsc2-/- Atf4-/- MEFs, and exogenous expression of ATF4 or SLC7A11, but not ATF4DBD, was able to restore glutathione levels to these cells (Figure 7F). Supplementation with all NEAAs or just cysteine, transported through neutral amino acid systems, but not equimolar concentrations of cystine, transported through SLC7A11, also rescued total glutathione levels, as measured by either enzymatic assay or LC-MS (Figure 7F, Figure 7—figure supplement 1C). Furthermore, insulin stimulated an increase in glutathione levels in wild-type MEFs in a manner completely sensitive to rapamycin, an effect ablated in Atf4 knockout cells, which had a very low abundance of glutathione (Figure 7G). Likewise, the rapamycin-sensitive nature of glutathione in Tsc2-/- MEFs was completely lost in Tsc2-/- Atf4-/- MEFs (Figure 7H). Glutathione levels were restored to these cells upon exogenous expression of ATF4, but not the ATF4DBD mutant. However, glutathione was still significantly reduced with rapamycin treatment in cells expressing the rapamycin-resistant ATF4, suggesting possible ATF4-independent mechanisms also contributing to this regulation (Lam et al., 2017). As one possible contributing factor, we found that the transcript encoding both the catalytic (GCLC) and regulatory (GCLM) subunits of glutamate-cysteine ligase, the first enzyme of glutathione synthesis, was sensitive rapamycin, in a manner unaffected by ATF4 knockdown (Figure 7—figure supplement 1D, E). We also found that GCLC and GCLM protein levels could be modestly induced by insulin through mTORC1 signaling in both wild-type and ATF4 knockout cells, but their protein abundance was unaffected by rapamycin in Tsc2-/- MEFs (Figure 7—figure supplement 1F, G). Thus, the mechanism underlying the apparent ATF4-independent effects of mTORC1 signaling on glutathione levels remains unknown. These collective data show that mTORC1 signaling induces glutathione synthesis, at least in part, through the activation of ATF4 and SLC7A11-dependent cystine uptake.
Discussion
Our findings expand the functional repertoire of mTORC1 signaling as it relates to the control of anabolic processes and cellular metabolism through its noncanonical activation of ATF4. Importantly, less than 10% of stress-responsive, ATF4-dependent targets were found to be significantly stimulated through the mTORC1-mediated activation of ATF4 in response to insulin. Among others, we found that genes involved in amino acid biosynthesis, transport, and tRNA charging were induced by mTORC1-ATF4 signaling, many to a comparable level to that of ER-stress induction with tunicamycin. While the molecular nature of this selective induction remains unknown, our data suggest that the 61 ATF4-dependent genes shared in their induction between mTORC1 signaling and the ISR represent targets most highly responsive to increases in ATF4 levels. Since mTORC1 signaling leads to a more modest increase in ATF4 protein levels than does the ISR, the selective induction of these genes might be reminiscent of the dose-dependent activation of MYC target genes documented in other studies (Sabò et al., 2014; Schuhmacher and Eick, 2013; Walz et al., 2014). Our bioinformatic analyses and functional data also indicate involvement of the C/EBP family of transcription factors as heterodimerization partners of ATF4 for the regulation of these gene targets shared between mTORC1 signaling and the ISR.
Consistent with the specific ATF4 target genes induced by mTORC1 signaling, including those involved in amino acid acquisition and tRNA charging, we found that ATF4 activation contributes to both canonical (e.g., protein synthesis) and new (e.g., glutathione synthesis) functions of mTORC1. As mTORC1 stimulates protein synthesis through multiple downstream targets (Holz et al., 2005; Jefferies et al., 1994; Ma and Blenis, 2009; Raught et al., 2004; Thoreen et al., 2012), it was not surprising to find that ATF4 was necessary but not sufficient for the increased rate of protein synthesis accompanying mTORC1 activation. We also demonstrate that mTORC1 signaling regulates the abundance of total cellular glutathione, both reduced and oxidized, at least in part through the ATF4-dependent induction of the cystine transporter SLC7A11, a major source of the cysteine that is limiting for glutathione synthesis. Importantly, this mTORC1-ATF4-mediated transcriptional upregulation of SLC7A11 leading to increased cystine uptake would temporally follow the inhibition of SLC7A11 recently found to be mediated through mTORC2 and Akt-dependent transient phosphorylation of the transporter (Gu et al., 2017; Lien et al., 2017). Our findings are consistent with a recent study indicating rapamycin-sensitive expression of xCT in TSC models, which the authors attribute to the OCT1 transcription factor (Li et al., 2019). However, our study indicates that the mTORC1-mediated activation of ATF4 is both necessary and sufficient for this regulation. The transcription factor NRF2 (also known as NFE2L2) is activated by oxidative stress and is a master regulator of the enzymes required for glutathione synthesis, as well as SLC7A11 to increase cystine uptake (Habib et al., 2015; Sasaki et al., 2002; Ye et al., 2014). While NRF2 depletion has been described to decrease the viability of cells with TSC gene loss (Zarei et al., 2019), we have no evidence from this or previous studies that mTORC1 signaling influences the levels or activity of NRF2 (Zhang et al., 2014). mTORC1 serves to couple growth signals to the coordinated control of anabolic processes, including the biosynthesis of protein, lipids, and nucleotides, as well as metabolic pathways that support this anabolic state (Valvezan and Manning, 2019). This metabolic program is orchestrated to provide biosynthetic precursors and directly promote the synthesis of macromolecules while also maintaining cellular homeostasis and preventing nutrient or metabolic stress. For example, mTORC1 signaling promotes metabolic flux through the NADPH-producing oxidative branch of the pentose phosphate pathway, thereby providing the reducing power essential to support an mTORC1-stimulated increase in de novo lipid synthesis (Düvel et al., 2010). Importantly, NADPH is also essential to reduce cystine, taken up through SLC7A11, into two molecules of cysteine for use in glutathione synthesis, in addition to being required to regenerate reduced glutathione following its oxidation. Supporting this logic of a coordinated metabolic program downstream of mTORC1, pro-growth signaling through mTORC1 likely promotes glutathione synthesis to help buffer against the oxidative stress that accompanies anabolic metabolism and increased rates of protein synthesis (Han et al., 2013; Harding et al., 2003; Kong and Chandel, 2018).
Our findings further support the addition of ATF4 to SREBP and HIF1, as nutrient- and stress-sensing transcription factors that are independently co-opted by mTORC1 signaling to drive the expression of metabolic enzymes and nutrient transporters. Unlike adaptive signals stemming from the depletion of individual nutrients, such as amino acids, sterols, or oxygen, which generally attenuate mTORC1 signaling as part of the adaptive response, pro-growth signals that activate mTORC1 can stimulate these transcription factors in concert to support a broader anabolic program. It will be important in future studies to understand the dual regulation of these transcription factors by both pro-growth and adaptive mechanisms as it relates to settings of physiological (fasting and feeding) and pathological (tumor development) nutrient fluctuations.
Materials and methods
Cell culture
Request a detailed protocolMEFs and PC3 cells were maintained in DMEM (Corning/Cellgro, 10-017-CV) with 10% fetal bovine serum (FBS, Corning/Gibco). LNCaP cells were maintained in RPMI-1640 (Corning/Cellgro 10-040-CV) with 10% FBS. Tsc2-/- (p53-/-) MEFs and littermate-derived wild-type counterparts were provided by David Kwiatkowski (Brigham and Women’s Hospital, Boston, MA). eIF2αS/S (WT) and eIF2αA/A (S51A knock-in mutant) MEFs were provided by Randal Kaufman (Sanford-Burnham-Prebys Medical Discovery Institute, La Jolla, CA) and were not used above passage 3 (after received). Rictor+/+ and Rictor-/- MEFs were provided by D.A. Guertin and D.M. Sabatini (Whitehead Institute, Massachusetts Institute of Technology, Cambridge, MA). Cancer cell lines were obtained from ATCC. Atf4-/- MEF lines generated in this study were maintained in DMEM with 10% FBS, supplemented with 55 µM 2-mercaptoethanol (Thermo, 21985023), and 1X MEM NEAA mix (NEAA, final concentrations: 100 µM each of alanine, aspartate, asparagine, glutamate, glycine, proline, and serine; Thermo 11140050). In experiments with supplementation of excess cysteine, cells were plated in DMEM with 10% FBS, 1 mM cysteine, and, where indicated, 1X MEM NEAA mix.
siRNA knockdowns
Request a detailed protocolCells were transfected with 20 nM of the indicated siRNAs using Opti-MEM (Thermo, 31985062) and RNAimax (Thermo, 13778150) according to the manufacturer’s protocol. siRNAs were from GE Life Sciences/Dharmacon: non-targeting pool (D-001810-10-50), Myc (L-040813-00-0010), Atf4 (mouse, L-042737-01-0020), Rheb (L-057044-00-0020), RhebL1 (L-056074-01-0020), Tsc2 (L-047050-00-0020), C/ebpα (L-040561-00-0005), C/ebpβ (L-043110-00-0005), C/ebpδ (L-060294-01-0005), C/ebpγ (L-065627-00-0005), and ATF4 (human, L-005125-00-0020). Forty-eight hours post transfection, cells were treated as indicated prior to lysis for immunoblotting, RNA extraction, or protein synthesis assays. For C/EBP isoform knockdown experiments, transfection of siRNAs was performed a second time, 24 hr after the first transfection. For protein synthesis assays involving siRNAs, transfection was also performed a second time, 24 hr after the first transfection, which was necessary to achieve sufficient knockdown of Rheb to reduce mTORC1 signaling.
Generation and validation of Atf4-/- and reconstituted cell lines
Request a detailed protocolTsc2+/+ (WT) and Tsc2-/- MEFs lacking Atf4 were generated by CRISPR-Cas9-mediated deletion using pSpCas9n(BB)−2A-GFP (PX461) vector (Addgene, 48140) according to the previously described protocol (Ran et al., 2013). The paired nickase guides were designed using E-CRISP (Heigwer et al., 2014) and targeted the sequences AGCATAGCCCCTCCACCTCC and GACAATCTGCCTTCTCCAGG in exon 2 of ATF4. Forty-eight hours post transfection, single GFP-positive cells were sorted into 96-well plates. Cells were cultured in DMEM with 10% FBS supplemented with 1X MEM NEAA and 55 µM 2-mercaptoethanol. Single cell clones were grown for immunoblot analysis, and those showing loss of ATF4 protein were selected for sequence analysis involving the isolation of genomic DNA (Qiagen, 69504), PCR amplification using KOD Xtreme Hot Start DNA Polymerase (Millipore, 71975), and the primers TCGATGCTCTGTTTCGAATG and CTTCTTCCCCCTTGCCTTAC flanking the targeted deletion site, with sequencing on an ABI3730xl DNA analyzer at the DNA Resource Core of Dana-Farber/Harvard Cancer Center (funded in part by NCI Cancer Center support grant 2P30CA006516-48). The mutations were identified using CRISP-ID software (Dehairs et al., 2016). For the final clones selected, the mutations generated in the WT MEFs include an out-of-frame 17-bp deletion starting at the codon encoding T237 and a large out-of-frame 73-bp deletion starting after the codon encoding G219, both resulting in premature STOP codons. The mutations generated in the Tsc2-/- MEFs include an out-of-frame 245-bp deletion starting at the codon encoding G190, resulting in a premature STOP codon after the D191 codon, and an in-frame deletion removing the sequences encoded between E210 and E284.
For generation of Atf4 expression vectors, the murine Atf4 cDNA was amplified from the plasmid 21845 from Addgene (Harding et al., 2000). Restriction enzyme cloning with AgeI and ClaI was used to insert the Atf4-coding sequence (lacking the 5′ and 3′-UTR) into the pTRIPZ plasmid for doxycycline-inducible expression. The ATF4DBD mutant, in which amino acids 292–298 of the DNA-binding domain are changed from RYRQKKR to GYLEAAA (He et al., 2001; Lange et al., 2008), was generated by DpnI-mediated site-directed mutagenesis using KOD Xtreme Hot Start DNA Polymerase. The dox-inducible pTRIPZ and SLC7A11 plasmids were a gift from Alex Toker (Beth Israel Deaconess Medical Center, Boston, MA). cDNA expression was induced with 1 μg/mL of doxycycline (Sigma-Aldrich, D3447) for 12–24 hr before assays were conducted. GFP was inserted into pTRIPZ to produce the control vector. Lentivirus was generated in HEK293T cells transfected with pMD2.G and psPAX2 (Addgene, 12259 and 12260) and the given pTRIPZ constructs. Forty-eight hours post transfection, the virus-containing medium was used to infect the Atf4 knockout cells, which were selected with 8 μg/mL puromycin. TSC2 addback cell lines were generated by retroviral infection following transfection of PT67 cells with pBabe hygro IRES-EV or pBabe hygro IRES-TSC2. Cells were selected with 400 μg/mL hygromycin B (Thermo, 10687010).
RNA-sequencing
Request a detailed protocolWild-type and Atf4-/- MEFs were grown to 70% confluence in 6 cm plates and were serum-starved in the presence of 2-mercaptoethanol and 1X MEM NEAA mixture and treated with vehicle (DMSO) or 20 nM rapamycin (LC Laboratories, R5000) for 30 min prior to stimulation with vehicle (water) or 500 nM insulin (Alpha Diagnostic, INSL 16 N-5) for 16 hr or treated with vehicle (DMSO) or 2 μg/mL tunicamycin (Sigma-Aldrich, T7765) for 4, 8, or 16 hr. RNA was harvested with TRIzol according to the manufacturer’s protocol (Thermo, 15596018). All samples passed RNA quality control measured by NanoDrop 1000 Spectrophotometer (NanoDrop Technologies) and 2100 Bioanalyzer (Agilent Technologies). cDNA libraries were generated to produce 150-bp paired-end reads on an Illumina NovaSeq with read depth of 20 million paired reads per sample. Reads were aligned and annotated to the Ensembl Mus musculus GRCm38.p6 genome assembly using the align and featureCounts functions from the Rsubread package (2.0.1) in R (3.6.3) (Liao et al., 2019). Differential gene expression analysis was performed using the voom and eBayes functions from the EgdeR (3.28.1) and Limma (3.42.2) packages, respectively (Ritchie et al., 2015; Robinson et al., 2010). Transcripts found to be significantly induced by tunicamycin were further limited to those with a greater than 1.2-fold increase. The enrichKEGG function from the clusterProfiler package (3.14.3) was used to perform KEGG pathway over-representation tests (Yu et al., 2012). Gene set enrichment analysis was evaluated using GSEA software from the Broad Institute (Subramanian et al., 2005). Computations were run on the FASRC Cannon cluster, supported by the Faculty of Arts and Sciences Division of Science Research Computing Group at Harvard University. Pseudogenes and unannotated genes were excluded from Figure 1C and Figure 1—source data 1. The complete RNA-seq data can be found at GEO under the accession number GSE158605.
CiiiDER analysis
Request a detailed protocolCiiiDER software was downloaded from CiiiDER.org with the M. musculus GRCm38.94 genome files. Searches were run against the JASPAR transcription factor-binding profile database. Searches were run on promoter regions spanning +1500 to −500 bp from the predicted transcriptional start site using a site identification deficit threshold of 0.1. The background gene list (ISR only) comprised the 200 ATF4-dependent genes most significantly increased in expression upon tunicamycin treatment that were not in the list of 61 genes shared in their regulation by mTORC1 signaling and the ISR. Results of this analysis are included in Figure 2—source data 1.
Cistrome analysis
Request a detailed protocolEach of the 61 shared mTORC1 and ISR genes was analyzed using the CistromeDB Toolkit (http://dbtoolkit.cistrome.org/) of existing genome-wide ChIP-seq data. A half-decay distance of 1 kb to the transcription start site was used. The top 20 transcription factors or chromatin regulators found in ChIP-seq experiments to bind to each gene were compiled, and the number of genes each factor bound to within the list of 61 was determined, with the top 9 regulators graphed, excluding the general factors EP300 and POL2RA.
Immunoblotting
Request a detailed protocolCells were lysed in ice-cold Triton lysis buffer (40 mM HEPES pH 7.4, 120 mM NaCl, 1 mM EDTA, 1% Triton X-100, 10 mM sodium pyrophosphate, 10 mM glycerol 2-phosphate, 50 mM NaF, 0.5 mM sodium orthovanadate, 1 μM Microcystin-LR, and Sigma protease inhibitor cocktail P8340). For immunoblots on SLC7A11, cells were lysed in 1% SDS lysis buffer (10 mM Tris pH 7.4, 1 mM EDTA, 10 mM sodium pyrophosphate, 10 mM glycerol 2-phosphate, 50 mM NaF, 0.5 mM sodium orthovanadate, and Sigma protease inhibitor cocktail P8340). Samples were centrifuged at 20,000 × g for 10 min at 4°C, and protein concentration in the supernatant was determined by Bradford assay (Bio-Rad, 5000202) and normalized across samples. Proteins were separated by SDS-PAGE, transferred to nitrocellulose membranes, and immunoblotted with indicated antibody. Primary antibodies used were MTHFD2 (Protein Tech, 12270-1-AP), PSAT1 (Protein Tech 20180-1-AP), phospho (P)-S6K1 T389 (Cell Signaling Technologies [CST], 9234), ATF4 (CST, 11815), eIF2α (CST, 9722), P-eIF2α S51 (CST, 9721), S6K1 (CST, 2708), CD98 (CST, 13180), AARS (Bethyl Antibodies, A303-475A-M), GARS (Bethyl Antibodies, A304-746A-M), LARS (Bethyl Antibodies, A304-316A-M), XPOT (Aviva Biotechnologies, ARP40711_P050), TSC2 (CST, 4308), SLC7A11 (mouse, Abcam, ab175186), SLC7A11 (human, CST, 12691), GCLC (Abcam, ab190685), GCLM (Abcam, ab126704), RICTOR (CST, 9476), RHEB (CST, 13879), 4EBP1 (CST 9644), c-MYC (CST, 9402), and β-actin (Sigma-Aldrich, A5316); secondary antibodies used were IRDye 800CW Donkey anti-Mouse IgG (H + L) (LI-COR, 926-32212) and Donkey anti-Rabbit IgG (H + L) (LI-COR, 925–32213), and HRP-conjugated anti-mouse and anti-rabbit secondary antibodies from (CST, 7074 and 7076). Immunoblots of MTHFD2, PSAT1, AARS, GARS, LARS, and XPOT were imaged using Odyssey CLx Imaging System (LI-COR Biosciences). β-Actin was developed with Odyssey CLx Imaging System or enhanced chemiluminescence assay (ECL). All remaining immunoblots were developed using ECL. Immunoblots were quantified using the Odyssey CLx Imaging System and were normalized to β-actin. The ATF4 immunoblot corresponding to Figure 1B (Figure 1—figure supplement 1A) and the SLC7A11 immunoblots in Figure 6D, E (Figure 6—figure supplement 1B, C) was quantified using ImageJ and normalized to β-actin.
NanoString analysis
Request a detailed protocolRNA was harvested using TRIzol from cells at 70% confluence. All samples passed RNA quality control measured by NanoDrop 1000 Spectrophotometer. Parallel plates were lysed for immunoblots. The isolated RNA was analyzed using a custom NanoString probe library according to the manufacturer’s instructions (NanoString Technologies). Briefly, sample RNA was hybridized to RNA-probes at 65°C for 16 hr, excess probe was washed away, and an nCounter SPRINT was used to quantify specific mRNA molecules present in each sample. Direct mRNA counts were normalized to internal control genes, and mRNA expression was analyzed using nSOLVER software. Heatmaps were generated with Morpheus software from the Broad Institute (https://software.broadinstitute.org/morpheus). Transcripts with 100 counts or fewer, a value based off of negative control samples, were not included in our analyses. Results of this analysis are included in Figure 3—source data 1.
qPCR
Request a detailed protocolFor gene expression analysis, RNA was isolated with TRIzol according to the manufacturer’s protocol. All samples passed RNA quality control measured by NanoDrop 1000 Spectrophotometer. cDNA was generated with Superscript III First Strand Synthesis System for RT-PCR (Thermo, 18080051). Quantitative RT-PCR was performed with a CFX Connect Realtime PCR Detection System (Bio-Rad) using iTaq Universal SYBR Green Supermix (Bio-Rad, 1725125). Samples were quantified by the ∆∆CT method, normalized to β-actin (mouse samples) or RPLP0 (human samples) to quantify relative mRNA expression levels. qPCR primers:
Atf4 | GATGGGTTCTCCAGCGACAAG | CCGGAAAAGGCATCCTCCTTC |
Psat1 | GCTGTCGCCTTAGCACCA | TGGATCTCCAACAATACCGAGTG |
Mthfd2 | TCCTTGTTGTCTGCGTTGGC | CTTCATTTCGCACTGCCGCC |
Slc7a5 | GGACAAGGTGATGCGTCCAA | GCCAACACAATGTTCCCCAC |
Aars | TTGCTATTCCCTCGGAGCAC | CTCCTCGGGAACCTTAGCTC |
Gars | GGCAGAGGTCTCTGAGCTG | GCACGATGGTCATAAGCTGC |
Cars | GAGCAGGCTGCCGACTACA | TATAGCTACGCGTGCTGAGG |
Nars | GAGCCGGCCTGTGTAAAGAT | GACCCAGCCAAACACCTTCA |
Iars | TATTGCATCACCTCCAGACGC | TGAACCATTCTGTTGCTGGGA |
Gclc | TGTACTCCACCCTCGTCACCC | CTGCTGTCCCAAGGCTCG |
Gclm | TGGGCACAGGTAAAACCCAA | CTGGGCTTCAATGTCAGGGA |
Slc7a11 | ATCTCCCCCAAGGGCATACT | GCATAGGACAGGGCTCCAAA |
Slc3a2 | TGATGAATGCACCCTTGTACTTG | GCTCCCCAGTGAAAGTGGA |
Slc1a5 | GTAAAATACCGCAATCCTGTATCC | CGATAGCGAAGACCACCAGG |
Xpot | GCTTCAGGCTCAGATGCAGA | AAAGCAAGGCGAACACTTGG |
c-Myc | AGAGCTCCTCGAGCTGTTTG | TTCTCTTCCTCGTCGCAGAT |
C/ebpα | CAAGAACAGCAACGAGTACCG | GTCACTGGTCAACTCCAGCAC |
C/ebpβ | CGCCTTATAAACCTCCCGCT | TGGCCACTTCCATGGGTCTA |
C/ebpδ | CGACTTCAGCGCCTACATTGA | CTAGCGACAGACCCCACAC |
C/ebpγ | TCGGATCACATTGCTCTGATTTC | TGTGCCTGAGTATGAATGACACT |
Actin | CACTGTCGAGTCGCGTCC | TCATCCATGGCGAACTGGTG |
PSAT1 | AAAAACAATGGAGGTGCCGC | GGCTCCACTGGACAAACGTA |
ASNS | TGGCTGCCTTTTATCAGGGG | TCTGCCACCTTTCTAGCAGC |
MTHFD2 | GGCAGTTCGAAATGAAGCTGTT | GCCAACCAGGATCACACTCA |
SLC7A5 | GACTACGCCTACATGCTGGA | AGCAGCAGCACGCAGAG |
AARS | CCATTCAGAAGGGCACAGGT | TATCCACGCCCTGTGTTGTC |
GARS | GCCAGCAGGGAGATCTTGTG | CCAGCTCCTTTGCTTCCAGA |
XPOT | GACGCAGAGCGACTAGAGG | TAAACATCTTCCCTATCACTCCATC |
SLC7A11 | AAGGTGCCACTGTTCATCCC | ATGTTCTGGTTATTTTCTCCGACA |
RPLP0 | CCTCGTGGAAGTGACATCGT | ATCTGCTTGGAGCCCACATT |
Protein synthesis assay
Request a detailed protocolCells were cultured as indicated, washed twice with PBS, and changed to methionine/cystine/glutamine-free DMEM (Thermo, 21013024) supplemented with 0.5 mM L-cystine (Sigma, 57579) and L-glutamine (4 mM) with 50 μCi/mL 35S-methionine (PerkinElmer, NEG009L005MC) for 20 min. Cells were washed twice in ice-cold PBS and lysed in ice-cold Triton lysis buffer. Total protein concentrations were normalized following a Bradford assay, and normalized samples were separated by SDS-PAGE and transferred to nitrocellulose. 35S-methionine incorporation into protein was analyzed by autoradiography, and relative rates of protein synthesis were quantified using ImageJ Software (NIH) to quantify radiolabeled protein per lane for each sample. For isogenic Atf4-/- cell lines, exogenous cDNA expression was induced for 16 hr with 150 ng/mL doxycycline, and cells were treated as indicated in the presence of doxycycline (150 ng/mL) and 1X MEM NEAA mixture plus 1 mM cysteine (Sigma, C7477). Protein synthesis was assayed with 20 min labeling in the absence of NEAA and cysteine to avoid competition of 35S-methionine uptake with the supplemented amino acids.
Methionine and cystine uptake
Request a detailed protocolFor methionine uptake assays, cells were cultured and labeled as described above for the protein synthesis assay, were washed three times in cold PBS, and lysed in Triton lysis buffer. For cystine uptake assays, cells were treated the same but labeled for the final 10 min with medium containing 0.1 µCi L-[1, 2, 1', 2'-14C]-Cystine (PerkinElmer, NEC854010UC) and washed three times in ice-cold PBS containing cold cystine (1 mM), prior to lysis in Triton lysis buffer. Whole-cell radiolabel incorporation was quantified with a Beckman LS6500 scintillation counter. Cells from identically treated parallel plates were counted using a Beckman Z1-Coulter Counter to normalize uptake measurements to cell number.
Proliferation assay
Request a detailed protocolTo quantify cell proliferation, cells were plated in DMEM in 6-well plates in triplicate in the presence of 2-mercaptoethanol, 1X MEM NEAA mixture, 10% FBS, and doxycycline (1 μg/mL). Twenty-four hours after plating, cells were washed twice with PBS and media was changed to DMEM with 10% dialyzed FBS, doxycycline (1 μg/mL), and the amino acid supplements indicated for each experiment, with media refreshed daily. Starting on day 0, viable cells from triplicate wells corresponding to each condition or cell line were counted using a hemocytometer, excluding dead cells detected by trypan blue stain (Sigma-Aldrich, T8154).
Analysis of cell death
Request a detailed protocolCells were plated in DMEM in 6-well plates in triplicate in the presence of 2-mercaptoethanol, 1X MEM NEAA mixture, 10% FBS, and doxycycline (1 μg/mL). Twenty-four hours after plating, cells were washed twice with PBS and media was changed to DMEM with 10% dialyzed FBS and doxycycline (1 μg/mL). Seventy-two hours later, cells were detached with Accumax (Sigma-Aldrich, A7089) and washed twice with cold PBS on ice. Cells were stained with annexin V and propidium iodide (PI) according to the manufacturer’s instruction (BD, 556547). Samples were analyzed using an LSRFortessa (BD) flow cytometer, and the fractions of stained cells were quantified using FloJo 10.6, with Staurosporine (4 hr, 5 μM) (Tocris, 1285) used as a positive control for cell death and to help establish gating of the sorted cells.
Measurements of cellular and tumor glutathione
Request a detailed protocolCells were plated in 96-well plates at 5000 cells/well. Total glutathione, GSH, and GSSG levels were measured using the GSH/GSSG-Glo Assay (Promega, V6611) according to the manufacturer’s protocol. Total glutathione levels were normalized to cell number determined from parallel plates. BSO (Sigma, B2515) was used as a positive control to inhibit glutathione synthesis.
For measurements via LC-MS/MS, metabolites were extracted from cells on dry ice using 80% methanol, and extracts were dried under nitrogen gas for metabolite profiling via selected reaction monitoring with polarity switching using a 5500 QTRAP mass spectrometer. Data were analyzed using MultiQuant 2.1.1 software (AB/SCIEX) to calculate the Q3 peak area. Normalized peak area of glutathione from human TSC2-/- angiomyolipoma (621-101) cells was determined from previously published data (Tang et al., 2019). For xenograft tumor studies, experimental details were provided previously (Valvezan et al., 2017). Briefly, mice bearing TSC2-/- ELT3 xenograft tumors were treated every other day for 5 days with vehicle or rapamycin (1 mg/kg on days 1, 3, and 5) and tumors were harvested for metabolite extraction, as above, 3 hr after the final treatment.
Statistics
For RNA-sequencing analysis, Benjamini–Hochberg false discovery rate (FDR)-adjusted p values were determined from empirical Bayes moderated t-statistics using the voom and eBayes functions from the limma package. Comparisons with FDR-adjusted p<0.05 were considered significant for the gene groups denoted compared to vehicle-treated controls. For KEGG enrichment, p values were FDR corrected. For CiiiDER transcription factor over-representation analysis, Fisher's exact test p values were used. Transcription factor binding elements with p<0.01 and test statistic >0 were considered over-represented in genes of interest. Unpaired two-tailed t-tests were used for NanoString analyses to calculate p values for rank ordering. All remaining statistical analyses were performed with Prism 8 software (GraphPad Software, La Jolla, CA). Statistical analyses for qPCR data with two treatment groups were determined by unpaired two-tailed t-test, while those with greater than two treatment groups were determined by one-way analysis of variance (ANOVA) with Holm–Sidak method for multiple comparisons. Statistical analyses for protein synthesis assays were determined by one-way ANOVA with Holm–Sidak method for multiple comparisons from values quantified with ImageJ software (US National Institutes of Health, Bethesda, MD). Statistical analyses for immunoblot quantification data with two treatment groups were determined by unpaired two-tailed t-test, while those with greater than two treatment groups were determined by one-way ANOVA with Holm–Sidak method for multiple comparisons. For proliferation assays, unpaired two-tailed t-test was used for comparisons to GFP-expressing cells. For cell death analysis, one-way ANOVA with Holm–Sidak method for multiple comparisons, summing the annexin V+/PI-, annexin V-/PI+, and annexin V+/PI+ populations for each conditions. For glutathione quantification of experiments with two conditions, an unpaired two-tailed t-test was performed. For remaining glutathione assays and all cystine and methionine uptake experiments, one-way ANOVA with Holm–Sidak method for multiple comparisons was used.
Source data
Request a detailed protocolThe source data for the RNA-sequencing experiment can be found at GEO under the accession number GSE158605. The source data for Figure 1C can be found in Figure 1—source data 1. The source data for Figure 2A can be found in Figure 2—source data 1. The source data for NanoString heatmaps shown in Figure 3 can be found in Figure 3—source data 1.
Data availability
All data generated or analyzed during this study are included in the manuscript and supporting files. Source data files have been provided for Figures 1, 2, and 3. RNA-Seq data have been deposited in GEO under accession code GSE158605.
-
NCBI Gene Expression OmnibusID GSE158605. The mTORC1-mediated activation of ATF4 promotes protein and glutathione synthesis (Tunicamycin).
References
-
Role of the transcription factor ATF4 in the anabolic actions of insulin and the anti-anabolic actions of glucocorticoidsJournal of Biological Chemistry 282:16744–16753.https://doi.org/10.1074/jbc.M610510200
-
The role of exportin-t in selective nuclear export of mature tRNAsThe EMBO Journal 17:7430–7441.https://doi.org/10.1093/emboj/17.24.7430
-
Transport interaction of L-cystine and L-glutamate in human diploid fibroblasts in cultureJournal of Biological Chemistry 255:2372–2376.https://doi.org/10.1016/S0021-9258(19)85901-X
-
CRISP-ID: decoding CRISPR mediated indels by Sanger sequencingScientific Reports 6:28973.https://doi.org/10.1038/srep28973
-
Activating transcription factor 4 (ATF4) promotes skeletal muscle atrophy by forming a heterodimer with the transcriptional regulator C/EBPβJournal of Biological Chemistry 295:2787–2803.https://doi.org/10.1074/jbc.RA119.012095
-
Potent and specific inhibition of glutathione synthesis by buthionine sulfoximine (S-n-butyl homocysteine sulfoximine)Journal of Biological Chemistry 254:7558–7560.https://doi.org/10.1016/S0021-9258(18)35980-5
-
Identification of activating transcription factor 4 (ATF4) as an Nrf2-interacting protein. implication for heme oxygenase-1 gene regulationThe Journal of Biological Chemistry 276:20858–20865.https://doi.org/10.1074/jbc.M101198200
-
E-CRISP: fast CRISPR target site identificationNature Methods 11:122–123.https://doi.org/10.1038/nmeth.2812
-
Transdominant suppression of estrogen receptor signaling by progesterone receptor ligands in uterine leiomyoma cellsMolecular and Cellular Endocrinology 196:11–20.https://doi.org/10.1016/S0303-7207(02)00230-7
-
SREBPs: activators of the complete program of cholesterol and fatty acid synthesis in the liverJournal of Clinical Investigation 109:1125–1131.https://doi.org/10.1172/JCI0215593
-
Regulation of hypoxia-inducible factor 1alpha expression and function by the mammalian target of rapamycinMolecular and Cellular Biology 22:7004–7014.https://doi.org/10.1128/MCB.22.20.7004-7014.2002
-
C/EBPγ is a critical regulator of cellular stress response networks through heterodimerization with ATF4Molecular and Cellular Biology 36:693–713.https://doi.org/10.1128/MCB.00911-15
-
ATF4-dependent transcription mediates signaling of amino acid limitationTrends in Endocrinology & Metabolism 20:436–443.https://doi.org/10.1016/j.tem.2009.05.008
-
Regulation of redox balance in Cancer and T cellsJournal of Biological Chemistry 293:7499–7507.https://doi.org/10.1074/jbc.TM117.000257
-
ATF4 is an oxidative stress-inducible, prodeath transcription factor in neurons in vitro and in vivoJournal of Experimental Medicine 205:1227–1242.https://doi.org/10.1084/jem.20071460
-
mTOR-dependent upregulation of xCT blocks melanin synthesis and promotes tumorigenesisCell Death & Differentiation 26:2015–2028.https://doi.org/10.1038/s41418-019-0274-0
-
Molecular mechanisms of mTOR-mediated translational controlNature Reviews. Molecular Cell Biology 10:307–318.https://doi.org/10.1038/nrm2672
-
Cistrome data browser: a data portal for ChIP-Seq and chromatin accessibility data in human and mouseNucleic Acids Research 45:D658–D662.https://doi.org/10.1093/nar/gkw983
-
Oxygen availability and metabolic adaptationsNature Reviews Cancer 16:663–673.https://doi.org/10.1038/nrc.2016.84
-
Genome engineering using the CRISPR-Cas9 systemNature Protocols 8:2281–2308.https://doi.org/10.1038/nprot.2013.143
-
Electrophile response element-mediated induction of the cystine/glutamate exchange transporter gene expressionJournal of Biological Chemistry 277:44765–44771.https://doi.org/10.1074/jbc.M208704200
-
MYC, metabolism, and cancerCancer Discovery 5:1024–1039.https://doi.org/10.1158/2159-8290.CD-15-0507
-
Nutrient sensing in cancerAnnual Review of Cancer Biology 2:251–269.https://doi.org/10.1146/annurev-cancerbio-030617-050329
-
Molecular logic of mTORC1 signalling as a metabolic rheostatNature Metabolism 1:321–333.https://doi.org/10.1038/s42255-019-0038-7
-
Surviving stress: modulation of ATF4-Mediated stress responses in normal and malignant cellsTrends in Endocrinology & Metabolism 28:794–806.https://doi.org/10.1016/j.tem.2017.07.003
-
Nrf2- and ATF4-dependent upregulation of xCT modulates the sensitivity of T24 bladder carcinoma cells to proteasome inhibitionMolecular and Cellular Biology 34:3421–3434.https://doi.org/10.1128/MCB.00221-14
-
clusterProfiler: an R package for comparing biological themes among gene clustersOMICS: A Journal of Integrative Biology 16:284–287.https://doi.org/10.1089/omi.2011.0118
-
Tumors with TSC mutations are sensitive to CDK7 inhibition through NRF2 and glutathione depletionJournal of Experimental Medicine 216:2635–2652.https://doi.org/10.1084/jem.20190251
-
Cistrome data browser: expanded datasets and new tools for gene regulatory analysisNucleic Acids Research 47:D729–D735.https://doi.org/10.1093/nar/gky1094
-
Modulation of hypoxia-inducible factor 1alpha expression by the epidermal growth factor/phosphatidylinositol 3-kinase/PTEN/AKT/FRAP pathway in human prostate Cancer cells: implications for tumor angiogenesis and therapeuticsCancer Research 60:1541–1545.
Article and author information
Author details
Funding
National Institutes of Health (R35-CA197459)
- Brendan D Manning
National Institutes of Health (P01-CA120964)
- John M Asara
- Brendan D Manning
U.S. Department of Defense (W81XWH-18-1- 0659)
- Brendan D Manning
National Institutes of Health (T32-ES016645)
- Margaret E Torrence
National Institutes of Health (F31-CA228332)
- Margaret E Torrence
Larimar Therapeutics
- Brendan D Manning
- James R Mitchell
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank David J Kwiatkowski, Elizabeth P Henske, Randal J Kaufman, Alex Toker, and David M Sabatini for cell lines and materials and Min Yuan, Gerta Hoxhaj, and Issam Ben-Sahra for technical assistance and critical discussions. This research was supported by grants from the NIH: T32-ES016645 (MET), F31-CA228332 (MET), R35-CA197459 (BDM), and P01-CA120964 (BDM and JA); the Department of Defense’s Congressionally Directed Medical Research Program on Tuberous Sclerosis Complex, award no. W81XWH-18-1-0659 (BDM); and a research grant from Zafgen (BDM and JRM). BDM is a shareholder and scientific advisory board member of LAM Therapeutics and Navitor Pharmaceuticals. All other authors declare no competing financial interests.
Ethics
Animal experimentation: All animal procedures were conducted under strict adherence to recommendations in the Guide for the Care and Use of Laboratory Animals of the National Institutes of Health and were approved by the Harvard Institutional Animal Care and Use Committee (#IS00000780).
Copyright
© 2021, Torrence 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
-
- 17,700
- views
-
- 2,091
- downloads
-
- 145
- 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
The prevalence and mortality rates of colorectal cancer (CRC) are increasing worldwide. Radiation resistance hinders radiotherapy, a standard treatment for advanced CRC, leading to local recurrence and metastasis. Elucidating the molecular mechanisms underlying radioresistance in CRC is critical to enhance therapeutic efficacy and patient outcomes. Bioinformatic analysis and tumour tissue examination were conducted to investigate the CPT1A mRNA and protein levels in CRC and their correlation with radiotherapy efficacy. Furthermore, lentiviral overexpression and CRISPR/Cas9 lentiviral vectors, along with in vitro and in vivo radiation experiments, were used to explore the effect of CPT1A on radiosensitivity. Additionally, transcriptomic sequencing, molecular biology experiments, and bioinformatic analyses were employed to elucidate the molecular mechanisms by which CPT1A regulates radiosensitivity. CPT1A was significantly downregulated in CRC and negatively correlated with responsiveness to neoadjuvant radiotherapy. Functional studies suggested that CPT1A mediates radiosensitivity, influencing reactive oxygen species (ROS) scavenging and DNA damage response. Transcriptomic and molecular analyses highlighted the involvement of the peroxisomal pathway. Mechanistic exploration revealed that CPT1A downregulates the FOXM1-SOD1/SOD2/CAT axis, moderating cellular ROS levels after irradiation and enhancing radiosensitivity. CPT1A downregulation contributes to radioresistance in CRC by augmenting the FOXM1-mediated antioxidant response. Thus, CPT1A is a potential biomarker of radiosensitivity and a novel target for overcoming radioresistance, offering a future direction to enhance CRC radiotherapy.
-
- Cancer Biology
- Evolutionary Biology
In growing cell populations such as tumours, mutations can serve as markers that allow tracking the past evolution from current samples. The genomic analyses of bulk samples and samples from multiple regions have shed light on the evolutionary forces acting on tumours. However, little is known empirically on the spatio-temporal dynamics of tumour evolution. Here, we leverage published data from resected hepatocellular carcinomas, each with several hundred samples taken in two and three dimensions. Using spatial metrics of evolution, we find that tumour cells grow predominantly uniformly within the tumour volume instead of at the surface. We determine how mutations and cells are dispersed throughout the tumour and how cell death contributes to the overall tumour growth. Our methods shed light on the early evolution of tumours in vivo and can be applied to high-resolution data in the emerging field of spatial biology.