ACLY and ACC1 Regulate Hypoxia-Induced Apoptosis by Modulating ETV4 via α-ketoglutarate

During the development of most solid tumors, there are characteristic physiological differences in the tumor that result from tumor cells outgrowing their local blood supply. Two of these physiological differences, or “stresses,” that occur in the tumor are low oxygen levels (hypoxia) and an accumulation of lactic acidic (lactic acidosis). Cancer cells experiencing hypoxia and lactic acidosis tend to be more resistant to chemo- and radio-therapy and metastasize more readily. Therefore, it is important to understand how tumor cells adapt to and survive these stresses. We used a large scale screening experiment in order to find which genes and proteins are involved in tumor cell adaptation and survival under hypoxia or lactic acidosis. We found that inhibiting either of two genes involved in lipid synthesis allowed tumor cells to survive hypoxia. This occurred because silencing these genes led to an increase in the metabolite α-ketoglutarate, which repressed a transcription factor that contributed to cell death under hypoxia. This research specifically advances our understanding of how tumor cells survive hypoxia and lactic acidosis and more broadly enhances our understanding of the cellular biology of solid tumors.

Published in the journal: . PLoS Genet 11(10): e32767. doi:10.1371/journal.pgen.1005599
Category: Research Article
doi: 10.1371/journal.pgen.1005599


During the development of most solid tumors, there are characteristic physiological differences in the tumor that result from tumor cells outgrowing their local blood supply. Two of these physiological differences, or “stresses,” that occur in the tumor are low oxygen levels (hypoxia) and an accumulation of lactic acidic (lactic acidosis). Cancer cells experiencing hypoxia and lactic acidosis tend to be more resistant to chemo- and radio-therapy and metastasize more readily. Therefore, it is important to understand how tumor cells adapt to and survive these stresses. We used a large scale screening experiment in order to find which genes and proteins are involved in tumor cell adaptation and survival under hypoxia or lactic acidosis. We found that inhibiting either of two genes involved in lipid synthesis allowed tumor cells to survive hypoxia. This occurred because silencing these genes led to an increase in the metabolite α-ketoglutarate, which repressed a transcription factor that contributed to cell death under hypoxia. This research specifically advances our understanding of how tumor cells survive hypoxia and lactic acidosis and more broadly enhances our understanding of the cellular biology of solid tumors.


Most solid tumors have substantial physiological deviations from normal tissue, which manifest as tumor microenvironment (TME) stresses [1,2]. These TME “stresses” include, among others, the limited availability of oxygen (hypoxia), glucose or amino acids, and an accumulation of lactic acid (lactic acidosis). In order to grow and propagate a solid tumor, tumor cells must adapt to and survive under these TME stresses. Additionally, tumor cells in regions of LA or hypoxia are more radio- and chemo-resistant and are more likely to metastasize [3]. Since TME stresses are found in most solid tumors, targeting stress-adaptation mechanisms of tumor cells may offer a significant therapeutic window to selectively eradicate tumor cells and improve patient outcomes [4]. Yet, current therapies targeting cells specifically under stress have significant limitations. For example, angiogenesis is a well-established, valuable therapeutic target with agents developed to block it at various stages of tumor development. However, many anti-angiogenic therapies fail over time, through acquired or inherited resistance that may involve the presence of tumor hypoxia [5]. Therefore, a number of other strategies are being developed to directly target hypoxic cells, such as blocking lactate transporters [6,7], or pro-drugs that are activated only in the presence of low oxygen [4]. There remains significant room for improvement to target cells under stress, and thus there remains a need to better understand the genes that impact cellular survival under TME stresses.

Tumor cells employ at least two kinds of adaptive strategies to cope with TME stresses, transcriptional and metabolic, and these are often interconnected. Transcriptional changes are mediated by stress-activated transcription factors. For example, the most significant factors known to regulate a cell’s hypoxia response are the hypoxia-inducible factors (HIFs), which are stabilized under low oxygen and mediate complex transcriptional programs that increase glucose uptake and enhance glycolysis [8,9]. HIF–1α also regulates glutamine metabolism by affecting the ubiquitination of its oxidizing enzyme, AKGDH, to promote reductive carboxylation of glutamine under hypoxic conditions [1012]. In addition to the HIFs, many other transcription factors regulate cellular responses to stresses in the TME. For example, the MondoA:Mlx complex senses and initiates transcriptional changes under both glucose deprivation and lactic acidosis to induce TXNIP and restrict glucose uptake [1315]. Importantly, multiple co-factors and modulators add to the complexity of stress mediated transcriptional responses [16]. While these adaptive mechanisms can be successful to sustain cell growth under stress, in multiple contexts hypoxia induces apoptosis [1719].

Transcriptional responses mediate many metabolic reprogramming events, but, recently, it is also becoming evident that metabolic events can regulate gene expression. For example, the metabolic enzyme ATP citrate lyase (ACLY) generates glucose-derived acetyl-CoA from citrate to alter histone acetylation and, therefore, transcriptional activation [20]. Importantly, the metabolically sensitive mTOR signaling cascade can activate HIF, even under normoxia [21,22]. Mutations in succinate dehydrogenase and fumarate hydratase both lead to increased levels of their substrates (succinate and fumarate, respectively), causing increased HIF–1α stability and alterations of genome-wide histone and DNA methylation [23,24]. Metabolites such as NAD+, NOS and α-ketoglutarate (α-KG) can also affect HIF function and histone and DNA modifications [21,22,2528]. While cellular metabolism and transcriptional changes can provide flexibility for adaptation, cancer cells can also become reliant on, and thus vulnerable to the inhibition of, specific metabolic pathways or gene products [2931]. Therefore, a better understanding of the genes necessary for modulating cancer cell survival under TME stresses will improve the development of targeted therapies that selectively eradicate cancer cells under stresses [4,32].

Functional genetic screens provide an unbiased and powerful means of identifying genes responsible for any phenotype that can be measured experimentally. Unbiased RNAi screens have identified genes that influence the survival of organisms and cells under various stresses [3336]. These studies provide a foundation for using RNAi screens to uncover genes involved in TME-relevant stress responses, but so far have not been applied genome-wide to identify genes that modulate the survival of cancer cells under hypoxia or lactic acidosis. To better understand the genes involved in the adaptations of cancer cells under TME stresses, we performed genome-wide pooled shRNA screens of lung cancer cells under hypoxia and lactic acidosis. Completing these screens revealed that the inhibition of ACC1 or ACLY, two key enzymes of de novo lipogenesis, protected cancer cells from hypoxia-induced apoptosis. ACC1 or ACLY inhibition protected cells by elevating levels of α-ketoglutarate under hypoxia to reduce the activity of the oncogenic transcription factor ETV4. Together, these data provide evidence to support a molecular connection between cellular metabolic and transcriptional hypoxia adaptation via the ACLY-ACC1-ETV4 axis through α-ketoglutarate.


Genome-wide pooled shRNA screens in hypoxia and lactic acidosis

To identify genes that modulate cell survival under lactic acidosis and hypoxia, we conducted genome-wide, shRNA-based, contextual pooled screens in the lung cancer cell line H1975 under hypoxia or lactic acidosis (Fig 1a). To preferentially discover genes important for survival rather than proliferation, the screen was done in low proliferative conditions (see Methods). Cells were transduced with a genome-wide MSCV-based shRNA library [37], selected for successful transduction with puromycin and then grown under control (21% O2, pH 7.4), hypoxia (2% O2, pH 7.4), or lactic acidosis (21% O2, 25mM lactic acid, pH 6.7) for 4 days, as performed previously [15,38]. At these stress treatments, there was a ~50% reduction in cell number, which allowed us to uncover both genes whose suppression reduced or improved survival under stresses. Genomic DNA was isolated from cells under each condition and the genome-incorporated shRNA sequences were amplified by PCR. The amplified PCR products were labeled and competitively hybridized to a custom microarray to identify those shRNA sequences that were either enriched or depleted relative to the control treatment (Fig 1a). The custom array was modified slightly from similar arrays used in other shRNA screens [39,40]. When different ratios of differentially labeled PCR products were hybridized on the arrays, we noted distinguishable differences in the signals, demonstrating the specificity and sensitivity of the array (S1a Fig). Biological triplicates of each condition had highly reproducible signals (S1b Fig). The abundance of each shRNA sequence reflected the effect of its target gene on cell survival under stresses: if the shRNA was depleted in the stress treatment, the gene it targeted had a “synthetic sick/lethal” phenotype; if the shRNA was enriched in the stress treatment, the gene it targeted had a “synthetic survival/protective” phenotype under stress. In order to analyze the effect of each shRNA in stress, we calculated an “R/G” ratio (see Methods). R/G ratios were distributed on a scale of +/- 4.0 that was highly consistent between replicates and stresses (S1c Fig).

Genome-wide shRNA screen identifies the depletion of ACC1 as protective for cells under hypoxia.
Fig. 1. Genome-wide shRNA screen identifies the depletion of ACC1 as protective for cells under hypoxia.
(A) The overall design of genome-wide pooled shRNA screen protocol. Each treatment was performed in n = 3. (B) Western blot showing decreased protein level expression of EPAS1 (HIF–2α) by shRNA. (C) Viable cell number counts by trypan blue exclusion after indicated cells in 4 days of hypoxia (n = 9). (D) Analysis of screen data by RIGER with GENE-E program with the probability of a gene being a hit indicated by 1/NES (normalized enrichment score). With the second best hairpin method, ACC1 is 13th best performing gene overall (indicated by red circle). (E) Venn diagrams showing the overlap of genes in the multiple hairpin analysis (see Methods) across both TME stress treatments, separated by shRNA enrichment (top) or depletion (bottom) in the stress condition. (F) GENE-E heat map results of ACC1 hairpins in hypoxia screen (top 3 rows) and lactic acidosis screen (bottom 3 rows). Each row is one biological replicate (n = 3 per treatment). Each column is a different shRNA targeting ACC1. Green box highlights the hypoxia result where the gene is predicted as a “protective” hit.

To minimize false positives due to off-target effects of individual shRNAs, we focused only on the genes that had at least two distinct shRNA sequences that were enriched or depleted (S1 and S2 Tables, see Methods). Importantly, this “multiple hairpin analysis” identified EPAS1 (hypoxia-inducible factor 2α, HIF–2α) as a synthetic lethal gene under hypoxia. We further validated this by showing that silencing EPAS1 by shRNA reduced cell survival under hypoxia (Fig 1b and 1c). This result was consistent with the critical role of EPAS1 in cellular adaptation to hypoxia [8]. The “re-discovery” of EPAS1 provided confidence in our screen and analysis methods. However, no pathways or gene ontology groups were significantly enriched within the different categories of “multiple hairpin hits”. We then performed a RIGER analysis using a log-fold change and the second best shRNA for each gene criteria (Fig 1d, S3 Table) [41]. This RIGER analysis revealed an enrichment for genes affecting mRNA regulation and binding, as well as membrane dynamics and nuclear localization [42]. Additionally, there was little overlap between the genes targeted by multiple hairpins, either enriched or depleted, in the two stresses (Fig 1e). This was consistent with past reports of distinct responses and adaptations to hypoxia and lactic acidosis [15,43].

Genome-wide screen identified ACC1 depletion as protective to cells under a hypoxic stress

Next, we identified high confidence “hits” in both the multiple hairpin and RIGER analyses for further investigation. From these considerations, we chose ACC1 (acetyl-CoA carboxylase 1 or ACACA) as the top candidate. ACC1 encodes the cytosolic isoform of acetyl-CoA carboxylase, which converts acetyl-CoA to malonyl-CoA in the rate-limiting step of de novo fatty acid synthesis. There was an enrichment of shRNAs targeting ACC1 in the hypoxic versus the control condition, suggesting that ACC1 knockdown allowed for improved survival under hypoxia. ACC1 had 4 hairpins enriched under hypoxia (Fig 1f) and scored as the 13th best gene in the RIGER analysis using the second best shRNA metric (Fig 1d). Additionally, the down-regulation of ACC1 was previously shown to protect cancer cells from glucose deprivation and matrix detachment stresses [44]. Together, these data prompted us to validate and investigate the role of ACC1 under hypoxia.

To validate the shRNA screen result, we silenced ACC1 expression through lentiviral infection of multiple shRNAs that targeted different sequences from those shRNAs used in the screen (see Methods). We confirmed the successful reduction of ACC1 protein by these shRNAs (Fig 2a). In the control cells transduced with a scramble shRNA, hypoxia significantly decreased cell viability and induced apoptosis (Fig 2). However, silencing ACC1 by multiple shRNAs inhibited the hypoxia-induced apoptosis as shown by crystal violet staining (Fig 2b), cell counting (Fig 2c), propidium-iodide staining (flow cytometry) (Fig 2d) and PARP cleavage (Fig 2e). This hypoxic protection associated with ACC1 silencing was also reproduced in additional hypoxia-sensitive cell lines, including MDA-MB–231 (breast cancer; S2a and S2b Fig), and PANC–1 (pancreatic cancer; S2e and S2g Fig). Furthermore, chemical inhibition of ACC1 through the AMPK agonist metformin also protected H1975 cells from hypoxia-induced apoptosis (S2h and S2i Fig). Collectively, these data successfully validated the screen results and showed that the depletion of ACC1 enhanced cell survival under hypoxia in multiple cancer cells from different tissues of origin.

Depletion of ACC1 or ACLY protects cells from hypoxia-induced apoptosis.
Fig. 2. Depletion of ACC1 or ACLY protects cells from hypoxia-induced apoptosis.
(A) Western blot showing efficiency of 4 independent shRNAs targeting ACC1 in H1975 cells. (B) Crystal violet staining of scramble and ACC1 shRNA cells after 4 days of hypoxia. (C) Viable cell numbers of indicated cells after 4 days of hypoxia as determined by counting nuclei (n = 9). (D) Percent of cells in the sub-G1 phase as determined by PI staining and FACS analysis of indicated cell line after 4 days of hypoxia (n = 9). (E) Western blot showing cleaved PARP levels with scramble and shACC1 cells after hypoxia for 48 hours. (F) Crystal violet staining of scramble and shACLY cell lines after 4 days of hypoxia. (G) Viable cell number determined by trypan blue exclusion after 4 days of hypoxia (n = 3). (H) Western blot showing ACLY protein levels and PARP cleavage in the indicated cells after 48 hours of hypoxia. Data are represented as mean values +/- SEM. Data are from the H1975 cell line.

Inhibition of ACC1 or ACLY inhibits hypoxia-induced apoptosis

We next investigated the specificity of the hypoxia protection by ACC1 depletion. Besides ACACA (ACC1), ACACB (ACC2) encodes another isoform of acetyl-CoA carboxylase, located in the outer mitochondrial membrane [45]. ACC2 was not a hit in our screen and its silencing by shRNA did not offer a similar hypoxia protection as seen with ACC1 depletion (S2j Fig). Next, we determined whether depletion of ACC1 protected against other TME stresses (lactic acidosis, glutamine deprivation and glucose deprivation). We found that the protective effect of shACC1 was seen only under hypoxia (S2k and S2l Fig). Since the HIFs are the major transcriptional responders to hypoxia, we examined how loss of ACC1 affected HIF–1α levels. Interestingly, we found decreased levels of HIF–1α under hypoxia with ACC1 depletion across multiple cell lines (S3a–S3c Fig). These data suggested that the protective phenotype was not due to upregulation of the HIF response. Overall, these results indicated that only the cytosolic isoform of acetyl-CoA carboxylase (encoded by ACC1) was essential for apoptosis, specifically under a hypoxic stress.

When we examined the effect of blocking enzymes up- or downstream of ACC1, we found that silencing ATP citrate lyase (ACLY) also enhanced survival under hypoxia (Fig 2f and 2g). ACLY encodes the enzyme immediately upstream of ACC1 in lipogenesis, catalyzing the formation of acetyl-CoA and oxaloacetate from citrate. Similar to ACC1, this protection results from the inhibition of hypoxia-induced apoptosis (Fig 2h). The protective effect of ACLY depletion was also reproduced in MDA-MB–231 cells (S2c and S2d Fig) and PANC–1 cells (S2f and S2g Fig). In H1975 and MDA-MB–231 shACLY cells there was also decreased HIF–1α expression under hypoxia as was seen in the shACC1 cells (S3d and S3e Fig). These data showed that blocking lipogenesis at the points of either ACLY or ACC1 inhibited apoptosis and permitted cell survival under hypoxia in cells of multiple tissue types.

Loss of ACC1 or ACLY did not protect cells from hypoxia-induced apoptosis through NADPH conservation to relieve oxidative stress

Lipogenesis is a highly anabolic process that requires significant amounts of NADPH and ATP. Previously, silencing ACC1 protected cells from death caused by glucose deprivation and matrix detachment by preserving NADPH and ATP to counteract the ensuing oxidative stresses [44,46]. We tested the relevance of these factors in our system. In H1975 cells, silencing ACC1 trended toward increasing the NADP+/NADPH ratio, suggesting a decrease in available NADPH (S4a Fig). We reasoned that if the NADPH were being used to combat elevated reactive oxygen species under hypoxia, then supplementation with antioxidants should protect control cells from hypoxia-induced death similar to [44]. However, neither the addition of N-acetyl cysteine nor glutathione antioxidants rescued hypoxia-induced death in control cells (S4b and S4c Fig). While ATP levels were higher with ACC1 silenced, the change in ATP levels from normoxia to hypoxia was consistent in control and knockdown cells and thus could not readily explain the hypoxia protection (S4d Fig). Therefore, in these cells with ACC1 or ACLY depleted, changes in NADPH and ATP levels may not be the primary mechanism for cell survival under hypoxia. Therefore, we sought to identify another mechanistic explanation for this hypoxia protection phenotype.

Depletion of the transcription factor ETV4 led to similar hypoxia-protection and gene expression phenotypes as with the loss of ACC1 or ACLY

In our “multiple hairpin analysis” of the hypoxia genome-wide screen, there was an enrichment of shRNAs targeting a PEA3 transcription factor family member, ETV4. These results suggested that silencing ETV4 may be protective under hypoxia. A link between lipogenesis and ETV4 was previously established when levels of malonyl-CoA were associated with ETV4 activity [47]. This prompted us to investigate a potential regulatory relationship between ACLY, ACC1 and ETV4. Real-time PCR analysis showed that hypoxia led to a reduction of ETV4 mRNA in the ACC1-depleted, but not control cells (Fig 3a). Additionally, we noted correspondingly reduced ETV4 protein in the shACC1 cells as compared to the scramble cells (Fig 3b). Reduced ETV4 mRNA and protein levels were also noted in the ACLY-depleted cells (S4e and S4f Fig). While ETV4 protein levels were somewhat decreased under normoxia, the down-regulation was stronger under hypoxia. This regulation was mostly specific to ETV4; the other PEA3 subfamily members, ETV1 and ETV5, were not consistently altered by ACC1 depletion (S4g Fig). Additionally, neither ETV1 nor ETV5 were identified as “multiple hairpin hits” in the shRNA screen. We validated the hypoxia-protective phenotype of ETV4 loss with two different shRNAs targeting ETV4 (Fig 3c and 3d). Similar to ACC1/ACLY depletion, the depletion of ETV4 also reduced the percentage of cells in the sub-G1 phase (Fig 3e) and decreased PARP cleavage (Fig 3f) under hypoxia. These data showed that the loss of ETV4 decreased hypoxia-induced apoptosis, similar to the phenotype of reduced ACC1 or ACLY.

Depletion of ETV4 protects cells from hypoxia-induced apoptosis.
Fig. 3. Depletion of ETV4 protects cells from hypoxia-induced apoptosis.
(A) qPCR analysis of ETV4 mRNA levels in shACC1 cells in normoxia and hypoxia (n = 6). (B) Western blot of ETV4 protein levels in shACC1 cell lines under normoxia and hypoxia. (C) Western blot of ETV4 protein levels in shETV4 cell lines under normoxia and hypoxia. (D) Crystal violet staining of two distinct shETV4 cell lines after 4 days of hypoxia. (E) Percent of cells in the sub-G1 phase as determined by PI staining of indicated cell line after 4 days of hypoxia (n = 6). (F) Western blot of intact and cleaved PARP in shETV4 cell lines under normoxia or hypoxia. Data are represented as mean values +/- SEM. Data are from the H1975 cell line.

Since ETV4 is a transcription factor, we investigated the contribution of reduced ETV4 activity to the transcriptional response of ACLY or ACC1 depletion. We used microarrays to analyze the global transcriptional response to the silencing of each ACC1, ACLY or ETV4 by two independent shRNAs under hypoxia (each shRNA was done in triplicate). The transcriptional responses were determined by zero-transformation against the shScramble cells [48]. Next, the data were filtered with a 1.7-fold change in at least six arrays and the selected 641 probesets were grouped by hierarchical clustering (Fig 4a). This analysis revealed a remarkable similarity between the transcriptional responses to the depletion of ETV4, ACC1 or ACLY with the induction and repression of common sets of genes (Fig 4a). Using the GATHER [49] algorithm, we noticed an “anti-apoptotic expression program” that included both the induction of negative regulators of apoptosis signaling such as NQO1, CYP1B1 and SERPINE1 [50,51] and the repression of the apoptosis-promoting genes BIK, TNFRSF9, TNFAIP3, GLIPR1, DDIT and TRIB3 (Fig 4a). The induction and repression of multiple genes in the shACC1, shACLY and shETV4 cells were confirmed by real-time qPCR (Fig 4b and 4c). Using GSEA, these gene expression changes were highly overlapping in all pairwise comparisons with both up- and down-regulated genes (S4h and S4i Fig).

Silencing ETV4 triggers similar gene expression changes as loss of ACC1 or ACLY.
Fig. 4. Silencing ETV4 triggers similar gene expression changes as loss of ACC1 or ACLY.
(A) Heat map of transcriptional responses to the depletion of indicated genes (ETV4, ACC1 or ACLY) by two shRNAs each (n = 3 per shRNA) under hypoxia. Data are log2 values, zero transformed to the scramble shRNA cells. Filtering criteria of at least six occurrences with values greater than 1.7 resulted in 641 probesets. (B,C) qPCR validation of mRNA changes of similarly (B) down-regulated and (C) up-regulated genes with depletion of ETV4, ACLY or ACC1 (n = 6). (D) Binary regressions of shACC1 (left panel) or shACLY (right panel) and shETV4 gene signatures when compared across gene expression patterns of 130 breast cancer tumor samples (see text and Methods).

Next, we evaluated whether the ACC1-affected genes were transcriptional targets of ETV4 by two different approaches. First, we compared publicly available ETV4 ChIP-seq data (Cistrome Finder [52]) from PC3 cells [53] with the genes that were changed in our microarray analysis of H1975 cells with loss of ETV4, ACC1 or ACLY (Fig 4a). While performed in a different cell (PC3), these analyses still identified at least two potential direct ETV4 target genes, PLEC (S5a Fig) and DUSP6 (S5b Fig). For both genes, there were peaks indicating direct ETV4 binding that overlapped with histone H3 lysine 27 acetylation (a mark of actively transcribed gene bodies) and DNase hypersensitive regions of open chromatin (S5a and S5b Fig). DUSP6 has been previously described as an ETS transcription factor family target [54,55]. While PLEC was reported to interact with vimentin, an ETV4 direct target [56], this analysis suggested that PLEC itself may represent a novel ETV4 target. In the second approach, we used qPCR to determine if the ACC1-affected genes could be “rescued” by ETV4 over-expression. We found that CTSS, COL13A1, DUSP6 and SERPINE1 could be reversed by ETV4 re-expression (S6a and S6b Fig). In contrast, the ACC1-altered expression of other genes was either partially or not restored upon ETV4 overexpression (S6c Fig). This analysis suggested that some of the gene expression changes discovered by our microarray analysis may represent direct effects of changed ETV4 transcription, while others likely represent more indirect changes with ETV4 loss. Together, these data indicated that the repression of ETV4 played an important role in a subset of the transcriptional response to ACC1 depletion.

In order to better understand if these changes reflected an in vivo biological regulation between these genes, we developed “gene signatures” associated with the silencing of ETV4, ACC1 or ACLY using the CreateSignature algorithm [57]. These gene expression signatures represent “quantitative phenotypes” that reflect the loss of these genes. Comparing their similarity in different expression datasets allowed us to recognize similar quantitative changes in these genes in both in vitro experimental perturbation and human tumors. Similar “gene signature” approaches have been used to define the influences of oncogenic signaling and TME stresses in multiple cancer types [5860]. Gene expression patterns from human tumor samples [61] were then separated by their similarity to our developed gene signatures associated with loss of ETV4 (shETV4), ACC1 (shACC1) and ACLY (shACLY). Binary regression from this analysis in human tumors showed highly statistically significant correlations between the shACC1 or shACLY and shETV4 signatures (Fig 4d). In other words, patient tumors with gene expression patterns more similar to the ACC1-depletion (shACC1) signature had expression patterns that were also more similar to the ETV4-depletion (shETV4) signature; likewise, patient tumors with gene expression patterns similar to the ACLY-depletion (shACLY) signature also had similar gene expression patterns with the ETV4-depletion (shETV4) signature. Importantly, this showed that the regulation between ACC1/ACLY and ETV4 was relevant in tumor expression datasets. Overall, these analyses demonstrated the similarity of the transcriptional responses to the depletion of ACC1, ACLY or ETV4 and suggested that ETV4 mediated a portion of the transcriptional effect downstream of ACLY or ACC1 both in vitro and in vivo.

Global metabolomics revealed that hypoxia-induced elevated α-ketoglutarate levels protected cells from apoptosis

Considering that ACC1 and ACLY are critical lipogenic enzymes, we performed a metabolomics experiment to interrogate the metabolic effects of ACC1 or ACLY depletion under normoxia or hypoxia. Five cell lines were evaluated in triplicate: 1 control “hypoxia-sensitive” cell line (shScramble line) and four “hypoxia-survival” cell lines (2 shACC1 lines, 2 shACLY lines). After 36 hours of treatment, cells were lysed on ice and collected to measure the intracellular levels of 15 amino acids and 45 acyl-carnitines by tandem mass spectrometry (MS/MS) and levels of 7 organic acids by gas chromatography and mass spectrometry (GC/MS). All measurements were normalized by total protein content per sample.

We observed several expected metabolic changes to validate our approach. Silencing ACC1 depleted basal and hypoxia-induced palmitate levels, reflecting reduced de novo lipogenesis in these cells (S7a Fig). Consistent with hypoxia-induced inhibition of pyruvate utilization in the TCA cycle in favor of anaerobic glycolysis, hypoxia modestly increased the levels of pyruvate and lactate in control cells (S7b and S7c Fig). In addition, as previously noticed [62], hypoxia generally reduced the levels of TCA metabolites succinate, fumarate, malate, and citrate in most cells (Fig 5a). These results indicated that our metabolomics assay accurately detected the expected metabolic changes associated with inhibited lipogenesis and hypoxia exposure.

Loss of ACC1 or ACLY results in elevated levels of α-KG under hypoxia.
Fig. 5. Loss of ACC1 or ACLY results in elevated levels of α-KG under hypoxia.
(A) Heat map showing percent change under hypoxia of various organic acids in indicated cell lines. Yellow = increase in hypoxia; Blue = decrease in hypoxia. Each shRNA cell line (2 per gene) was done in triplicate (n = 3). (B) Protein-normalized fold change of intracellular α-KG levels in hypoxia in indicated cell lines. Dashed line indicates no change in hypoxia from normoxia (n = 3). (C) Normalized intracellular levels of α-KG after supplementation with indicated concentration dimethyl-α-KG for 24 hours (n = 3). (D) Western blot of intact and cleaved PARP after α-KG supplementation in control cells under normoxia or hypoxia for 48 hours. Data are represented as mean values +/- SEM, unless otherwise noted. Data are from the H1975 cell line.

The pattern of α-ketoglutarate (α-KG) levels in this experiment suggested that it may be an interesting candidate for offering protection under hypoxia. In the control cells, hypoxia reduced the levels of α-KG. However, in the hypoxia-resistant cells with depleted ACC1 and ACLY, hypoxia increased the α-KG levels (Fig 5b), as has been seen before in hypoxia-resistant cells [62]. We could test the possibility that levels of α-KG contributed to survival by adding cell-permeable dimethyl-α-KG to H1975 cells. We determined the level of α-KG achieved intracellularly after extracellular supplementation to choose a supplementation treatment that would be relevant to the levels of α-KG seen with ACC1 or ACLY depletion. Mass spectrometry analysis showed increasing amounts of intracellular α-KG after supplementation in a dose-dependent manner (Fig 5c) and that dimethyl-α-KG supplementation at 1mM achieved levels of α-KG comparable to the hypoxia-induced increase found in the ACC1 and ACLY depleted cells under hypoxia (Fig 5b). Supplementation of relevant levels α-KG also inhibited PARP cleavage under hypoxia (Fig 5d). These data indicated that the increased α-KG under hypoxia in the ACC1 and ACLY depleted cells recapitulated the hypoxia-protective phenotype of these cells.

Hypoxia-induced increased α-ketoglutarate levels regulated ETV4, possibly through 2-oxoglutarate/Fe(II)-dependent dioxygenases

With both α-KG and ETV4 acting downstream of ACC1 or ACLY, we next determined if α-KG was mediating the effects of ACC1 or ACLY on ETV4 expression. α-KG supplementation reduced mRNA and protein levels of ETV4 at both lower (Fig 6a and 6b) and higher concentrations (S7d and S7e Fig). Additionally, α-KG supplementation in control cells caused similar changes in repressed and induced genes as was caused by ACC1, ACLY or ETV4 silencing (Fig 6c and 6d), and some of these mRNA effects were dose-dependent with α-KG supplementation (S7f and S7g Fig). Unexpectedly, diM-α-KG supplementation increased HIF–1α protein levels (S7h Fig). Since the regulation of HIF–1α was different with depletion of ACC1/ACLY or α-KG supplementation, the changes in HIF–1α protein did not likely explain the improved hypoxic cell survival of the ACC1/ACLY depleted cells. Combined, these data indicated that, in the ACLY or ACC1 depleted cells, the α-KG increase was a hypoxic trigger that reduced ETV4 levels and activity to mediate an anti-apoptotic gene expression response.

α-KG regulates ETV4 mRNA levels and activity through epigenetic changes.
Fig. 6. α-KG regulates ETV4 mRNA levels and activity through epigenetic changes.
(A) qPCR showing ETV4 mRNA levels after α-KG supplementation (1mM) over time (n = 6). (B) Western blot showing ETV4 protein levels after α-KG supplementation (1mM) over time. (C,D) qPCR of gene expression changes after α-KG supplementation (3.5mM, 24h) for both (C) down-regulated genes (similar to Fig 4b) and (D) up-regulated genes (similar to Fig 4b and 4c) (n = 6). (E) Ratio of α-KG/succinate in normoxia and hypoxia in the indicated cell lines from the metabolomics experiment (n = 3). (F) Western blot showing effect of α-KG supplementation (1mM) on histone methylation marks after indicated time of treatment. Numbers indicate densitometry analysis by Image J relative to that sample’s loading of total H3. Data are represented as mean values +/- SEM. Data are from the H1975 cell line.

In addition to being a TCA cycle intermediate, α-KG is a substrate for the abundant 2-oxoglutarate/Fe(II)-dependent dioxygenases (2-OGDDs) [26]. 2-OGDDs use α-KG and molecular oxygen as substrates to perform a number of different protein modification reactions. These enzymes include families of histone demethylases that recognize and remove methylation marks from histones, as well as the TET family of proteins that facilitate DNA demethylation [26]. Thus, α-KG levels can affect gene expression through the activities of these 2-OGDDs [28]. An elevated ratio of α-KG/succinate (substrate/product ratio) has been proposed as a potential indicator of increased 2-OGDD activity [25,28]. We found that the α-KG/succinate ratio was significantly elevated in all of the shACC1 and shACLY cells under hypoxia (Fig 6e). To better understand if the level of α-KG or the α-KG/succinate ratio determined our hypoxia-survival phenotypes, we supplemented ACLY or ACC1 depleted cells with cell-permeable dimethyl-succinate to theoretically drive the α-KG/succinate ratio in the opposite direction from when α-KG was added. Succinate supplementation did not affect the survival of either shACC1 or shACLY cells under hypoxia (S8a and S8b Fig) and also did not affect the regulation of ETV4 by ACC1 or ACLY (S8c and S8d Fig). Therefore, in our experimental system, we concluded that the levels of α-KG, rather than the α-KG/succinate ratio, were driving the hypoxia survival phenotypes.

While the succinate supplementation did not affect our phenotype, the elevated α-KG in the ACC1 or ACLY depleted cells under hypoxia still suggested that 2-OGDDs may be relevant in our system. Therefore, we hypothesized that α-KG affected ETV4 mRNA abundance by altering the activity of 2-OGDDs and subsequent histone methylations. As a control, we tested the DNA methylation status of the two shore regions and the center region of the ETV4 promoter CpG islands by bisulfite pyro-sequencing and saw no significant change upon α-KG supplementation (S8e Fig). However, α-KG caused a global reduction in two (H3K4me2 and H3K4me3) “active” histone methylation marks and also globally reduced two “repressive” marks, H3K27me3 and, to a lesser extent, H3K9me3 (Fig 6f). When we compared the epigenetic changes associated with α-KG supplementation with Carey et al.[28], we found that H3K27me3 was consistently reduced with the addition of α-KG in both studies. However, there were also differences in the histone methylation changes caused by α-KG across the studies: 1) H4K20me3 was reduced with α-KG previously and we saw an increase in this mark (S8f Fig); 2) whereas we saw a decrease in the levels of H3K4me3, no changes were seen previously. These differences could be due to different cellular contexts (embryonic stem cells vs. cancer cells) or the presence of glutamine deprivation during the previous examination of α-KG effects [28]. Similar changes were observed using either water or DMSO as control (S8g Fig). Collectively, the reduced levels of multiple methylation marks were consistent with our hypothesis that predicted more active histone demethylases as a result of increased levels of α-KG.

To extend evidence in support of our model, we also examined various histone methylation modifications associated with the depletion of ACC1 or ACLY. Overall, we saw similar global histone methylation changes in both the shACC1 (Fig 7a) and shACLY (Fig 7b) cells as compared to α-KG supplementation. Among all the tested epigenetic markers, the H3K4me3 mark was the most pronouncedly decreased across both gene depletions and the α-KG treatment. The H3K4me2 mark was decreased in all three conditions (shACC1, shACLY, α-KG supplementation) to a modest degree. Similarly, H3K9me3 was decreased somewhat with the 8 hour α-KG treatment and in the shACC1 cells, while it was more strongly decreased in the shACLY cells. H3K27me3 was lowered by both α-KG and ACC1 depletion. Levels of H4K20me3 were unchanged in the shACC1 and shACLY cells while they were increased with α-KG treatment, and so this suggested that the changes in this methyl mark were likely not due to the changes we explain in our model.

Histone methylation changes in shACC1 and shACLY cells under hypoxia.
Fig. 7. Histone methylation changes in shACC1 and shACLY cells under hypoxia.
(A) Western blot showing levels of the indicated histone methylation marks in shACC1 cells under normoxia or hypoxia (24h). (B) Western blot showing levels of the indicated histone methylation marks in shACLY cells under normoxia or hypoxia (8h). (C) ChIP-qPCR analysis of the abundance of the H3K4me3 mark at two locations in the promoter region of ETV4 after α-KG treatment (n = 6). (D, E) ChIP-qPCR analysis of the abundance of the H3K4me3 mark at two locations in the promoter region of ETV4 after ACC1 or ACLY depletion under normoxia and hypoxia (n = 9). n.s. = not significant; * = p<0.05, ** = p<0.01, *** = p<0.005, ***** = p<0.0001. Data are represented as mean values +/- SEM. Data are from the H1975 cell line.

Besides global epigenetic changes, we also determined if histone methylation at the ETV4 locus was changed by chromatin immunoprecipitation of the “active” histone H3 lysine 4 tri-methylation (H3K4me3) mark. This mark was chosen because it showed the most consistent and strongest changes across either genetic depletion or with α-KG supplementation. In addition, a loss of this “active” mark H3K4me3 would be consistent with decreased ETV4 expression in these conditions. ChIP experiments showed that α-KG treatment decreased the abundance of the H3K4me3 modification at the ETV4 locus by ~30% (Fig 7c). Additionally, there was a decreased abundance of H3K4me3 in both ACC1 and ACLY depleted cells under hypoxia as compared to control cells (Fig 7d and 7e). Together, these data showed that elevated levels of α-KG affected the histone, but not DNA, methylation status of the ETV4 locus and this pattern was similar to the histone methylation changes seen under hypoxia in the shACC1 and shACLY cells.

Collectively, our data was consistent with a model in which, under hypoxia, the inhibition of ACC1 or ACLY increases levels of α-ketoglutarate to block hypoxia-induced apoptosis by reducing the levels and activity of ETV4, possibly through altered histone methylation patterns (Fig 8). These data offer a novel molecular connection showing that the transcriptional output of altered lipogenic metabolism can modulate the cellular response to and survival under hypoxia.

Model for how loss of ACLY, ACC1, or ETV4 protects cells from hypoxia-induced apoptosis.
Fig. 8. Model for how loss of ACLY, ACC1, or ETV4 protects cells from hypoxia-induced apoptosis.


Here we have described a pooled shRNA screen that successfully identified genes that influence cancer cell survival under hypoxia and lactic acidosis. Specifically, we showed that blocking de novo lipogenesis through the genetic depletion of ACLY or ACC1 protected multiple cancer cells from hypoxia-induced apoptosis through increased levels of α-ketoglutarate and the inhibition of ETV4 and its transcriptional activities. Therefore, inhibition of ACLY or ACC1 affected both metabolism and transcription to protect cells from hypoxia-induced apoptosis. These results also suggest that one important mechanism of hypoxia-induced apoptosis is through the reduction of α-KG, which potentially elevates the levels and activity of ETV4 through histone modifications to promote oncogenesis and trigger apoptosis.

Inhibitors that target ACLY and ACC1 are proposed cancer therapeutics [63,64]. ACC1 has been targeted as it is the rate-limiting enzyme in lipogenesis, a process critical for cancer cells’ rapid proliferation [63]. For ACLY, not only is it a critical enzyme for active lipogenesis, but it also regulates epigenetic states by generating acetyl-CoA, and so it offers two mechanisms to target clinically [20]. Our data provides several insights on the biological effects of ACC1 and ACLY in hypoxic cancer cells that should be considered when targeting these enzymes. First, while we expect ACLY and ACC1 inhibition to have opposite effects on levels of acetyl-CoA, their effects on gene expression were highly similar. This suggests that levels of acetyl-CoA (and subsequent histone acetylation or changes in epigenetic states) may not readily explain the majority of gene expression responses to the inhibition of ACC1 or ACLY under hypoxia. Instead, the ACC1/ACLY-induced reduction in ETV4 levels and activity seemed to account for a significant portion of the hypoxic transcriptional changes. These results suggest that lipogenic inhibitors that block ACC1 or ACLY may be particularly effective for tumors driven by ETV4. Likewise, this reduction of the oncogenic driver ETV4 may account for portions of the therapeutic potential of lipogenic inhibitors. However, our data also showed that these treatments may allow for the survival of slowly proliferating cancer cells in regions of hypoxic tumors. Therefore, targeting lipogenesis in cancer may need to be combined with other therapeutic approaches that target hypoxic regions (such as hypoxia-activated pro-drugs), to eliminate the cancer cells that may be dormant and protected in the hypoxia-niche.

Elevated α-KG levels in the shACLY and shACC1 cells under hypoxia could be due to increased generation or decreased consumption of this metabolite. A portion of α-KG is consumed during lipogenesis, which is further promoted under hypoxia by reductive carboxylation [12,65]. Therefore, blocking lipogenesis under hypoxia may lead to a “build-up” of upstream metabolites, including α-KG. It is well appreciated that α-KG plays a critical role in supporting cell survival by replenishing the metabolic intermediates of the TCA cycle. Here, our data indicated than an additional manner by which α-KG can affect cellular survival under hypoxia was by regulating ETV4 expression, possibly through epigenetic mechanisms.

Transcriptional adaptation to hypoxia is most often orchestrated by the HIFs; however, here we showed that cancer cells’ hypoxic survival can be mediated by a different transcription factor, ETV4. ETV4 was proposed as an essential co-activator of HIF–1α and to have a hypoxic transcriptional program [66]. Our data revealed that ETV4 was critical for hypoxia-induced apoptosis. Interestingly, the levels of ACC1 or ACLY influenced α-KG and ETV4 levels under hypoxia, thereby providing a link between lipogenesis, a TCA cycle intermediate and transcription. Therefore, ETV4 mediated the transcriptional output of varying degrees of active lipogenesis caused by changing ACC1 and ACLY levels. Our gene signature analyses also suggested that this regulation was preserved in human tumors in vivo. While the global and local ETV4 epigenetic changes we describe herein were consistent with decreased activity of ETV4 to protect cells from hypoxia-induced apoptosis, our evidence provides only correlative support with reduced promoter activity and does not provide a causal explanation for the ETV4 repression. Additionally, although this regulation seemed mostly specific to ETV4, we do not fully exclude the possibility that other ETS transcription factor family members contribute to an apoptotic, hypoxic transcriptional program. Overall, our data establishes ETV4 as one critical factor that influences hypoxic cell survival and transcriptional responses downstream from ACC1 and ACLY, and reveals the metabolite α-ketoglutarate as a molecular link between metabolic and transcriptional adaptation to hypoxia.

In addition to its importance in the TCA cycle, our experiments showed that elevated α-KG could alter histone methylation patterns, likely via α-ketoglutarate-dependent dioxygenases, to potentially regulate ETV4. While the causal relationship of several regulatory steps of this hypothesis was not rigorously tested by genetic manipulations here, our data was consistent with such a model. Previous reports suggested that α-KG/succinate ratios determined the direction and activities of dioxygenases [28], yet our data indicated that it was increased α-KG, not succinate or α-KG/succinate, that drove the epigenetic changes, ETV4 repression and hypoxia survival phenotypes. As a substrate, α-KG levels likely affect the enzymatic activities of of the 2-OGDDs histone demethylases, each of which have different sets of lysine residues upon which they act [26]. The number and combination of global histone methylation events that were changed with depletion of ACC1, ACLY or α-KG treatment suggested that many different demethylases could be affected by α-KG or the depletion of ACC1/ACLY in these cellular states. While there are no histone methyltransferases currently known to use α-KG as a substrate, they are regulated by hypoxia [67,68], and so we do not exclude the potential for changes in histone methylation patterns to be due to changes in both methyltransferase and demethylase enzyme activities. Additionally, the changes in global histone methylation patterns are very likely affecting a number of other genes’ expression patterns, in addition to the changes we see at the ETV4 locus.

Since α-ketoglutarate supplementation and loss of ACC1 or ACLY reduced the abundance of H3K4me3 at the ETV4 promoter, we speculate that the JARID1 (KDM5) family of histone demethylases, which specifically demethylate H3K4me2/3 residues, could be affecting the abundance of this mark at the ETV4 locus [69,70]. These enzymes are known to range in their expression and activity by cell type and are differentially influenced by oxygen levels [7173]. A full investigation in to the mRNA, protein and activity levels of many of these family members would be necessary to determine the extent to which each plays a role in regulating ETV4 under hypoxia. Effort has been made to investigate the JARID1 family of H3K4 demethylases as potential cancer drug targets [70,74,75] as the importance of a mis-regulated “histone code” for tumorigenesis is well recognized. While some reports suggest a tumor suppressive role of these enzymes, more suggest an oncogenic function [69]. This information and the data we presented here suggest that various 2-OGDDs in distinctive contexts differentially affect tumorigenesis and tumor cell survival. It will be important to understand the proper context of treatment if these drugs continue into the clinic.

Multiple previous reports show that the activity of ACLY [64,76,77], ACC1 [7881] or ETV4 [8285] is associated with increased tumorigenicity and/or poor patient outcome, or that inhibiting these genes’ activities reduces tumorigenicity and improves patient outcome. We show that the inhibition of ACLY, ACC1 or ETV4 paradoxically allows tumor cells to survive better under hypoxia. To address this apparent conundrum, we propose a conceptual model in which there are two cellular states: one of activating oncogenesis and the other of stress survival. This model includes a trade-off between the two states, such that the promotion of one comes at the expense of the other. Stated another way, the activation of oncogenic programs may also render cells susceptible to apoptosis, especially under stress; likewise, reduced oncogenic programs slow cellular proliferation to a “dormant” state that could allow for better stress survival. A similar model has been proposed for several oncogenes. For example, oncogenesis driven by MYC rendered non-transformed cells vulnerable to hypoxia-induced apoptosis [86,87]; the degradation or cleavage of c-MYC under hypoxia allowed tumor cells to evade hypoxia-induced apoptosis [88,89]. E2F can promote proliferation (oncogenesis) or apoptosis in different contexts, such as with differing PI3K activity [90]. A recent paper also indicated that HIF–1α repressed the ATF4 stress response pathway to allow for the expansion of fetal cardiomyocytes [91]. According to this proposed model and our data, cancer cells treated with ACLY or ACC1 inhibitors (including metformin), may die due to blocked lipogenesis, but may also survive in hypoxic regions. As these cells resist death under stress, they may become the “dormant” cells that recur after treatment regimens end, but they also may become more targetable as they persist under hypoxia. Consistent with past literature and clinical attempts to target these genes, our model advises that treatment regimens be carefully considered. While data to prove such a model remains necessary, this study suggests that the ACLY-ACC1-ETV4 axis might mediate the balance between oncogenesis and stress survival under hypoxia. It is very likely that similar mechanisms of balance between proliferation and stress survival exist in a wide variety of biological contexts.

Materials and Methods

shRNA screen and analysis

H1975 cells were spin-infected with the pMSCV-based retroviral genome-wide library, at an MOI of 0.3, divided into six sub-pools, achieving a final library representation of 1000 cells per shRNA after selection with 1 ug/ml puromycin [37]. After three days of puromycin selection, cells were split into control and stress conditions, maintaining 1000-fold representation of each shRNA per triplicate. Cells were serum starved to 0.1% FBS 24 hours after plating. 24 hours after serum starvation, media was changed to treatment media (0.1% FBS, 25mM Hepes); control and hypoxia media pH = 7.4; lactic acidosis treatment had 25mM lactic acid (Sigma cat. no L6402) adjusted to pH = 6.7 and filter sterilized. After 4 days of treatment, cells were harvested, centrifuged, and frozen at -80°C. Genomic DNA was extracted with the QIAamp DNA Blood Maxi kit (QIAGEN, 51194) then shRNA sequences were PCR amplified. The amplified products from the control and each stress were labeled (Cy3 and Cy5, respectively) and then interrogated by a custom Agilent microarray, which contained probes against the library’s shRNA sequences [37]. We validated the sensitivity and specificity of the array to different ratios of labeled PCR product (S1a Fig). Therefore, relative hybridization of the Cy5/Cy3 labeled shRNA populations determined the abundance of each shRNA under control, hypoxia or LA. The Cy3 and Cy5 signals across the three biological replicates were highly reproducible (S1b Fig). Probes with signal intensities of less than 2-fold above background were discarded. Cy5/Cy3 ratios, also called “R/G” ratios, for remaining probes were calculated, log2 transformed and quantile normalized across pools. The R/G ratios ranged from +/- 4.0, although many fell in the “unchanged” range of +/-0.5 (S1c Fig).

For the “multiple hairpin analysis,” genes were considered a hit when they had 1) at least 2 different shRNAs with (absolute value R/G) > 0.7 in at least 2 of the three biological replicates (2) the (stdev/ave) of the biological replicates was <0.5 (S1 Table).

Cell culture, TME stress treatments and generation of stable shRNA cell lines

H1975 cells were cultured in RPMI media (GIBCO cat. no 11875) supplemented with 10% Fetal bovine serum (heat-inactivated), 1% glucose, 10mM HEPES, 1mM sodium pyruvate, and 1x antibiotics (penicillin, 10,000 UI/ml; streptomycin, 10,000 UI/ml), as directed by the Duke Cell Culture Facility. MDA-MB–231 and PANC–1 cells were cultured in DMEM (GIBCO 11995) supplemented with 10% Fetal bovine serum (heat-inactivated) and 1x antibiotics (penicillin, 10,000 UI/ml; streptomycin, 10,000 UI/ml). Cell lines, obtained from and initially validated by the Duke Cell Culture Facility (Durham, NC, USA), were maintained for fewer than 6 months and validated by microscopy every 1 to 2 days.

Lactic acidosis was generated via addition of lactic acid (Sigma-Aldrich, St. Louis, MO, USA, cat. no L6402) and media pH adjustment to pH 6.7 by HCl immediately before use. Hypoxia was generated with a cell culture incubator with 93–94% N2, 5% CO and 1–2% O2. For the α-KG rescue experiments, media was supplemented with 0.875-4mM dimethyl α-KG as indicated in figure legends (Sigma, cat. no. 349631). For all stress experiments, cells were serum starved (0.5% FBS) for 24 hours before treated with stress under 0.5% FBS. All survival/viability measurements were made after 4 days of stress treatment.

Stable cell line generation

Stable cell lines were created with the pLKO.1 shRNA constructs purchased from the Duke RNAi Core Facility. Virus was generated by transfecting HEK-293T cells with a 1: 0.1: 1 ratio of pMDG2: pVSVG: pLKO.1 with Lipofectamine 2000 in the evening. Media was changed the following morning and virus collected 48 hours after transfection. Stable cell lines were generated by adding 200ul virus to a 60mm dish of parental cells with polybrene (final concentration 8ug/ml). Complete death in blank infection dishes was used to determine success of infection and puromycin selection. The efficiency of silencing or overexpression was determined by western blots. Concentrations of puromycin needed for selection: H1975 cells = 1ug/ml, MDA-MB–231 cells = 1ug/ml, PANC–1 cells = 2ug/ml. For stable overexpression, concentration of blasticidin used was 2.5ug/ml in H1975 cells.

Crystal violet staining

Cells were fixed either in 4% paraformaldehyde (PFA) overnight at 4°C or at room temperature for 30 min. PFA was removed and crystal violet staining solution (0.2% crystal violet, 25% methanol, 75% water) gently shaken on cells for 30+ minutes at room temperature. Staining solution was removed and plates rinsed with tap water 2–3 times. For quantitation, completely dried stain was dissolved by adding 10% acetic acid and shaking gently at room temperature for 30+ min before reading absorbance at 570 nm.

Determination of cell number

Cell number was evaluated by either direct cell counting (trypan blue exclusion) or high-throughput microscopic counting (HTC) of fixed and stained nuclei. For direct cell counting, at designated time after treatment, media was removed, cells were not rinsed for fear of losing loosely-attached cells, trypsinized, diluted 1:1 with trypan blue and immediate counted on a hemocytometer. For HTC experiments, after designated time period, cells were fixed in 4% PFA either overnight at 4°C or for 30 min at RT. Cells were washed 2x, permeabilized with 0.1% Triton-X in PBS, wash 2x, stained with 50ug/ml Hoescht dye (Sigma cat. no B2261) for 30 min, RT in the dark, then washed 2x and PBST added to each well and scanned by the Cellomics high-throughput microscope at the Duke RNAi Core Facility.

Flow cytometry

For cell cycle analysis, after 4 days of stress treatment, media was collected, cells trypsinized and pooled with the media. Cells were centrifuged then fixed by resuspension in ice cold 70% ethanol while gently vortexing. Fixed cells were placed at -20°C until prepared for FACS analysis. Immediately before FACS analysis, cells were centrifuged for 5 min at room temperature, washed twice in PBS then resuspended in 25ug/ml Propidium iodide (Sigma cat. no P4864) and 10ug/ml RNAse A in PBS. Cells were stained for 30+ min in the dark then 8000 events measured on a Canto II Flow cytometer.

Protein lysate collection and Western blots

Cell lysis: Cells were washed once with ice cold PBS, lysed by RIPA buffer with protease and phosphatase inhibitors added fresh, scraped into a microcentrifuge tube, allowed to swell on ice for 15–20 min, vortexed briefly, then spun down at top speed for 15 min at 4°C. Supernatant was transferred to pre-cooled new tube and protein concentration assayed with the Pierce BCA kit (ThermoScientific, cat. no. 23225). Effort was made to immediately rinse and lyse cells coming from a hypoxia condition as cells very quickly re-equilibrate to normoxic conditions. Western blots: Between 15-30ug of lysate was loaded on SDS-PAGE gels, wet-transferred to PDVF membrane, blocked with 5% milk in 1xTBST (0.1% Tween–20), then primary antibodies were incubated overnight at 4°C. For analysis of histones, protein was extracted with the EpiQuick Total Histone Extraction Kit (Epigentek, OP–0006) and 2ug of protein were resolved on 15% SDS-PAGE gels or nuclear fractions were collected by the REAP fractionation method [92] and 7.5–30 ul of lysate were run in each lane. Please contact for details on antibody usage.

Quantitative real-time PCR

RNA was extracted using the RNeasy Kit (QIAGEN). A total of 1 μg of total RNA was reverse transcribed by SuperScript II (Invitrogen) for real-time PCR with Power SYBRGreen Mix (Applied Biosystems/Life Technologies (Grand Island, NY, USA)). Primers were designed across exons whenever possible, verified for specificity by regular PCR prior to use in real-time PCR. Please contact for the sequences of primers used.

Microarrays and analysis

Samples were collected on ice and RNA was isolated with QIAGEN’s RNeasy Mini Kit (cat. no 74104) according to manufacturer’s instructions. After quality control assessment with the Agilent BioAnalyzer, cDNA was amplified from 200ng RNA with the Ambion MessageAmp Premier RNA Amplification (Life Technologies, Grand Island NY, USA). The gene expression pattern of the RNA samples were interrogated with Affymetrix U133A genechips and normalized by the RMA (Robust Multi-Array) algorithm. cDNA synthesis and microarray interrogation was performed by the Duke Microarray Core. The influence of the silencing of ACC1 or ACLY on gene expression was derived by a zero transformation process, in which we compared transcript level for each gene in cells with stably integrated shRNAs targeting ACC1 or ACLY to the average transcript levels in control scramble shRNA cell line samples. Data was then filtered as described with Cluster 3.0 software and heat maps were generated with TreeView. To generate gene signatures of knockdown of ACLY, ACC1 or ETV4, the CreateSignature module in GenePattern ( was used with scramble cells expression pattern as the train0 set, knockdown cells’ gene expression pattern as train1 set and the Gray dataset [61] used as the test set. Default parameters were used for the analysis similar to [58]. The resulting probabilities of gene signature expression in each patient for each knockdown signature were analyzed by simple linear regression in the JMP Pro 11 software.

Statistical analysis

Data, unless otherwise noted, represent the mean +/- the standard error of the mean and n indicates number of replicates used to generate the SEM. P-values were determined either by a two-tailed Student’s t-test in Excel or by a two-way ANOVA with StatView.

Metabolomics profiling and analysis

The measurement of amino acids and acyl carnitines was performed using stable isotope dilution techniques and flow injection tandem mass spectrometry and mass spec sample preparation methods described previously [9395]. Derivatized organic acids were analyzed by capillary gas chromatography/mass spectrometry (GC/MS) using a TRACE DSQ instrument (Thermo Electron Corporation; Austin, TX) [9395]. All MS analyses employed stable-isotope-dilution. The standards serve both to help identify each of the analyte peaks and provide the reference for quantifying their levels. Quantification was facilitated by addition of mixtures of known quantities of stable-isotope internal standards from Isotec (St. Louis, MO), Cambridge Isotope Laboratories (Andover, MA), and CDN Isotopes (Pointe-Claire, Quebec, CN) to samples. Sample Preparation: Biological triplicates of 15cm plates under each treatment condition were placed on ice and washed twice with ice-cold PBS before as much PBS as possible was removed. Cells were lysed in 620ul ml of 0.78% Formic Acid in water and scraped to collect. 30ul were removed for protein quantification. 1x volume of the collected pellet (~800ul) of acetonitrile was added and sample was vortexed vigorously. Aliquots were separated for mass spectrometry measurements (300 ul for organic acids, 100 ul for amino acids/acyl-carnitines) and were immediately frozen on dry ice and transferred to -80°C. Analysis: Protein concentration per replicate was determined by the Pierce BCA Kit (ThermoScientific, Waltham, MA, USA) and used to normalize all metabolite levels.

NADP+/NADPH measurements

A ratio of NADP+/NADPH was calculated after measuring each molecule separately with the Amplite Fluorimetric NADP/NADPH Ratio Assay Kit from AAT Bioquest, Inc (Sunnyvale, CA). Protocol was conducted as the manufacturer suggested and all values were normalized to protein content, as measured by the Pierce BCA kit, on similarly plated and treated samples done in parallel.

Chromatin Immunoprecipitation

3 million H1975 cells were plated in 15 cm dishes, after 24 hours they were serum starved to 0.5% FBS, 24 hours later they were treated with either a 4-hour treatment of α-KG or an 8 hour treatment of hypoxia before collection for a native ChIP. Protocol was carried out as manufacturer suggested with the SimpleChIP Plus Enzymatic Chromatin IP Kit (Agarose Beads) (Cell Signaling Tech., cat. no. 9004). Sonication and digestion were performed to obtain chromatin 1–4 nucleosomes in size, which was verified by gel electrophoresis. Each IP used chromatin prepared from 4–5 million cells and was performed overnight.

DNA Methylation

Genomic DNA was extracted with the DNeasy Blood & Tissue Kit according to the protocol provided by the manufacturer (Qiagen). The genomic DNAs (800 ng) were modified by treatment with sodium bisulfite using the Zymo EZ DNA Methylation kit (Zymo Research, Irvine, CA). Bisulfite treatment of denatured DNA converts all unmethylated cytosines to uracils, leaving methylated cytosines unchanged, allowing for quantitative measurement of cytosine methylation status. Pyrosequencing was performed using a Pyromark Q96 MD pyrosequencer (Qiagen). The bisulfite pyrosequencing assays were used to quantitatively measure the level of methylation at CpG sites contained. Assays were designed to query CpG islands using the Pyromark Assay Design Software (Qiagen). Pyrosequencing was performed using the sequencing primer. PCR conditions were optimized to produce a single, robust amplification product. Defined mixtures of fully methylated and unmethylated control DNAs were used to show a linear increase in detection of methylation values as the level of input DNA methylation increased (Pearson r > 0.98 for all regions). Once optimal conditions were defined, each assay was analyzed using the same amount of input DNA from each specimen (40 ng, assuming complete recovery after bisulfite modification). Percent methylation for each CpG cytosine was determined using Pyro Q-CpG Software (Qiagen).

Supporting Information

Attachment 1

Attachment 2

Attachment 3

Attachment 4

Attachment 5

Attachment 6

Attachment 7

Attachment 8

Attachment 9

Attachment 10

Attachment 11


1. Gatenby RA, Smallbone K, Maini PK, Rose F, Averill J, et al. (2007) Cellular adaptations to hypoxia and acidosis during somatic evolution of breast cancer. Br J Cancer 97: 646–653. 17687336

2. Gatenby RA, Gillies R.J. (2008) A microenvironmental model of carcinogenesis. Nature Reviews Cancer 8: 56–61. 18059462

3. Cardone RA, Casavola V, Reshkin SJ (2005) The role of disturbed pH dynamics and the Na+/H+ exchanger in metastasis. Nat Rev Cancer 5: 786–795. 16175178

4. Wilson WR, Hay MP (2011) Targeting hypoxia in cancer therapy. Nat Rev Cancer 11: 393–410. doi: 10.1038/nrc3064 21606941

5. Abdollahi A, Folkman J (2010) Evading tumor evasion: Current concepts and perspectives of anti-angiogenic cancer therapy. Drug Resistance Updates 13: 16–28. doi: 10.1016/j.drup.2009.12.001 20061178

6. Neri D, Supuran CT (2011) Interfering with pH regulation in tumours as a therapeutic strategy. Nat Rev Drug Discov 10: 767–777. doi: 10.1038/nrd3554 21921921

7. Sonveaux P, Vegran F, Schroeder T, Wergin MC, Verrax J, et al. (2008) Targeting lactate-fueled respiration selectively kills hypoxic tumor cells in mice. J Clin Invest 118: 3930–3942. doi: 10.1172/JCI36843 19033663

8. Majmundar AJ, Wong WJ, Simon MC (2010) Hypoxia-Inducible Factors and the Response to Hypoxic Stress. Molecular Cell 40: 294–309. doi: 10.1016/j.molcel.2010.09.022 20965423

9. Bertout JA, Patel SA, Simon MC (2008) The impact of O2 availability on human cancer. Nat Rev Cancer 8: 967–975. doi: 10.1038/nrc2540 18987634

10. Sun RC, Denko NC (2014) Hypoxic regulation of glutamine metabolism through HIF1 and SIAH2 supports lipid synthesis that is necessary for tumor growth. Cell Metab 19: 285–292. doi: 10.1016/j.cmet.2013.11.022 24506869

11. Metallo CM, Gameiro PA, Bell EL, Mattaini KR, Yang J, et al. (2012) Reductive glutamine metabolism by IDH1 mediates lipogenesis under hypoxia. Nature 481: 380–384.

12. Wise DR, Ward PS, Shay JE, Cross JR, Gruber JJ, et al. (2011) Hypoxia promotes isocitrate dehydrogenase-dependent carboxylation of alpha-ketoglutarate to citrate to support cell growth and viability. Proc Natl Acad Sci U S A 108: 19611–19616. doi: 10.1073/pnas.1117773108 22106302

13. Stoltzman CA, Peterson CW, Breen KT, Muoio DM, Billin AN, et al. (2008) Glucose sensing by MondoA:Mlx complexes: A role for hexokinases and direct regulation of thioredoxin-interacting protein expression. Proceedings of the National Academy of Sciences 105: 6912–6917.

14. Peterson CW, Stoltzman CA, Sighinolfi MP, Han K-S, Ayer DE (2010) Glucose Controls Nuclear Accumulation, Promoter Binding, and Transcriptional Activity of the MondoA-Mlx Heterodimer. Molecular and Cellular Biology 30: 2887–2895. doi: 10.1128/MCB.01613-09 20385767

15. Chen JL-Y, Merl D, Peterson CW, Wu J, Liu PY, et al. (2010) Lactic Acidosis Triggers Starvation Response with Paradoxical Induction of TXNIP through MondoA. PLoS Genet 6: e1001093. doi: 10.1371/journal.pgen.1001093 20844768

16. Pawlus MR, Hu C-J (2013) Enhanceosomes as integrators of hypoxia inducible factor (HIF) and other transcription factors in the hypoxic transcriptional response. Cellular Signalling 25: 1895–1903. doi: 10.1016/j.cellsig.2013.05.018 23707522

17. Hammond EM, Giaccia AJ (2005) The role of p53 in hypoxia-induced apoptosis. Biochem Biophys Res Commun 331: 718–725. 15865928

18. Bacon AL, Harris AL (2004) Hypoxia-inducible factors and hypoxic cell death in tumour physiology. Ann Med 36: 530–539. 15513303

19. Chinnadurai G, Vijayalingam S, Gibson SB (2008) BNIP3 subfamily BH3-only proteins: mitochondrial stress sensors in normal and pathological functions. Oncogene 27 Suppl 1: S114–127. doi: 10.1038/onc.2009.49 19641497

20. Wellen KE, Hatzivassiliou G, Sachdeva UM, Bui TV, Cross JR, et al. (2009) ATP-Citrate Lyase Links Cellular Metabolism to Histone Acetylation. Science 324: 1076–1080. doi: 10.1126/science.1164097 19461003

21. Semenza GL (2010) HIF–1: upstream and downstream of cancer metabolism. Curr Opin Genet Dev 20: 51–56. doi: 10.1016/j.gde.2009.10.009 19942427

22. Cairns RA, Harris IS, Mak TW (2011) Regulation of cancer cell metabolism. Nat Rev Cancer 11: 85–95. doi: 10.1038/nrc2981 21258394

23. King A, Selak MA, Gottlieb E (2006) Succinate dehydrogenase and fumarate hydratase: linking mitochondrial dysfunction and cancer. Oncogene 25: 4675–4682. 16892081

24. Pollard PJ, Briere JJ, Alam NA, Barwell J, Barclay E, et al. (2005) Accumulation of Krebs cycle intermediates and over-expression of HIF1alpha in tumours which result from germline FH and SDH mutations. Hum Mol Genet 14: 2231–2239. 15987702

25. Kaelin WG Jr (2011) Cancer and altered metabolism: potential importance of hypoxia-inducible factor and 2-oxoglutarate-dependent dioxygenases. Cold Spring Harb Symp Quant Biol 76: 335–345. doi: 10.1101/sqb.2011.76.010975 22089927

26. Loenarz C, Schofield CJ (2008) Expanding chemical biology of 2-oxoglutarate oxygenases. Nat Chem Biol 4: 152–156. doi: 10.1038/nchembio0308-152 18277970

27. Quintero M, Brennan PA, Thomas GJ, Moncada S (2006) Nitric Oxide Is a Factor in the Stabilization of Hypoxia-Inducible Factor–1α in Cancer: Role of Free Radical Formation. Cancer Research 66: 770–774. 16424008

28. Carey BW, Finley LW, Cross JR, Allis CD, Thompson CB (2015) Intracellular alpha-ketoglutarate maintains the pluripotency of embryonic stem cells. Nature 518: 413–416. doi: 10.1038/nature13981 25487152

29. Van der Heiden MG (2011) Targeting cancer metabolism: a therapeutic window opens. Nat Rev Drug Discov 10: 671–684. doi: 10.1038/nrd3504 21878982

30. Wise DR, DeBerardinis RJ, Mancuso A, Sayed N, Zhang X-Y, et al. (2008) Myc regulates a transcriptional program that stimulates mitochondrial glutaminolysis and leads to glutamine addiction. Proceedings of the National Academy of Sciences 105: 18782–18787.

31. Yun J, Rago C, Cheong I, Pagliarini R, Angenendt P, et al. (2009) Glucose Deprivation Contributes to the Development of KRAS Pathway Mutations in Tumor Cells. Science 325: 1555–1559. doi: 10.1126/science.1174229 19661383

32. Parks SK, Chiche J, Pouyssegur J (2013) Disrupting proton dynamics and energy metabolism for cancer therapy. Nat Rev Cancer 13: 611–623. doi: 10.1038/nrc3579 23969692

33. Saito S, Furuno A, Sakurai J, Sakamoto A, Park H-R, et al. (2009) Chemical Genomics Identifies the Unfolded Protein Response as a Target for Selective Cancer Cell Killing during Glucose Deprivation. Cancer Research 69: 4225–4234. doi: 10.1158/0008-5472.CAN-08-2689 19435925

34. Mabon ME MX, Jiao Y, Scott BA, Crowder CM. (2009) Systematic identification of gene activities promoting hypoxic death. Genetics 181: 13.

35. Dekanty A, Romero NM, Bertolin AP, Thomas MG, Leishman CC, et al. (2010) Drosophila genome-wide RNAi screen identifies multiple regulators of HIF-dependent transcription in hypoxia. PLoS Genet 6: e1000994. doi: 10.1371/journal.pgen.1000994 20585616

36. Pan J, Zhang J, Hill A, Lapan P, Berasi S, et al. (2013) A kinome-wide siRNA screen identifies multiple roles for protein kinases in hypoxic stress adaptation, including roles for IRAK4 and GAK in protection against apoptosis in VHL-/- renal carcinoma cells, despite activation of the NF-kappaB pathway. J Biomol Screen 18: 782–796. 23591012

37. Schlabach MR, Luo J, Solimini NL, Hu G, Xu Q, et al. (2008) Cancer Proliferation Gene Discovery Through Functional Genomics. Science 319: 620–624. doi: 10.1126/science.1149200 18239126

38. Chen JL-Y, Lucas JE, Schroeder T, Mori S, Wu J, et al. (2008) The Genomic Analysis of Lactic Acidosis and Acidosis Response in Human Cancers. PLoS Genet 4: e1000293. doi: 10.1371/journal.pgen.1000293 19057672

39. Schlabach MR, Luo J, Solimini NL, Hu G, Xu Q, et al. (2008) Cancer proliferation gene discovery through functional genomics. Science 319: 620–624. doi: 10.1126/science.1149200 18239126

40. Luo J, Emanuele MJ, Li D, Creighton CJ, Schlabach MR, et al. (2009) A genome-wide RNAi screen identifies multiple synthetic lethal interactions with the Ras oncogene. Cell 137: 835–848. doi: 10.1016/j.cell.2009.05.006 19490893

41. Gould J GENE-E. Broad Institute.

42. Eden E, Navon R, Steinfeld I, Lipson D, Yakhini Z (2009) GOrilla: a tool for discovery and visualization of enriched GO terms in ranked gene lists. BMC Bioinformatics 10: 48. doi: 10.1186/1471-2105-10-48 19192299

43. Tang X, Lucas JE, Chen JL-Y, LaMonte G, Wu J, et al. (2012) Functional Interaction between Responses to Lactic Acidosis and Hypoxia Regulates Genomic Transcriptional Outputs. Cancer Research 72: 491–502. doi: 10.1158/0008-5472.CAN-11-2076 22135092

44. Jeon S-M, Chandel NS, Hay N (2012) AMPK regulates NADPH homeostasis to promote tumour cell survival during energy stress. Nature 485: 661–665. doi: 10.1038/nature11066 22660331

45. Abu-Elheiga L, Brinkley WR, Zhong L, Chirala SS, Woldegiorgis G, et al. (2000) The subcellular localization of acetyl-CoA carboxylase 2. Proceedings of the National Academy of Sciences 97: 1444–1449.

46. Schafer ZT, Grassian AR, Song L, Jiang Z, Gerhart-Hines Z, et al. (2009) Antioxidant and oncogene rescue of metabolic defects caused by loss of matrix attachment. Nature 461: 109–113. doi: 10.1038/nature08268 19693011

47. Menendez JA, Vellon L, Mehmi I, Oza BP, Ropero S, et al. (2004) Inhibition of Fatty Acid Synthase (FAS) Suppresses HER2/Neu (ErbB–2) Oncogene Overexpression in Cancer Cells. Proceedings of the National Academy of Sciences of the United States of America 101: 10715–10720. 15235125

48. Tang X, Keenan MM, Wu J, Lin CA, Dubois L, et al. (2015) Comprehensive profiling of amino Acid response uncovers unique methionine-deprived response dependent on intact creatine biosynthesis. PLoS Genet 11: e1005158. doi: 10.1371/journal.pgen.1005158 25849282

49. Chang JT, Nevins JR (2006) GATHER: a systems approach to interpreting genomic signatures. Bioinformatics 22: 2926–2933. 17000751

50. Zhu E-D, Li N, Li B-S, Li W, Zhang W-J, et al. (2014) miR-30b, Down-Regulated in Gastric Cancer, Promotes Apoptosis and Suppresses Tumor Growth by Targeting Plasminogen Activator Inhibitor–1. PLoS ONE 9: e106049. doi: 10.1371/journal.pone.0106049 25170877

51. Zhang YP, Wang WL, Liu J, Li WB, Bai LL, et al. (2013) Plasminogen activator inhibitor–1 promotes the proliferation and inhibits the apoptosis of pulmonary fibroblasts by Ca(2+) signaling. Thromb Res 131: 64–71. doi: 10.1016/j.thromres.2012.09.003 23021499

52. Sun H, Qin B, Liu T, Wang Q, Liu J, et al. (2013) CistromeFinder for ChIP-seq and DNase-seq data reuse. Bioinformatics 29: 1352–1354. doi: 10.1093/bioinformatics/btt135 23508969

53. Hollenhorst PC, Ferris MW, Hull MA, Chae H, Kim S, et al. (2011) Oncogenic ETS proteins mimic activated RAS/MAPK signaling in prostate cells. Genes Dev 25: 2147–2157. doi: 10.1101/gad.17546311 22012618

54. Znosko WA, Yu S, Thomas K, Molina GA, Li C, et al. (2010) Overlapping functions of Pea3 ETS transcription factors in FGF signaling during zebrafish development. Dev Biol 342: 11–25. doi: 10.1016/j.ydbio.2010.03.011 20346941

55. Seifert AW, Yamaguchi T, Cohn MJ (2009) Functional and phylogenetic analysis shows that Fgf8 is a marker of genital induction in mammals but is not required for external genital development. Development 136: 2643–2651. doi: 10.1242/dev.036830 19592577

56. Chen JH, Vercamer C, Li Z, Paulin D, Vandenbunder B, et al. (1996) PEA3 transactivates vimentin promoter in mammary epithelial and tumor cells. Oncogene 13: 1667–1675. 8895512

57. Chang JT, Gatza ML, Lucas JE, Barry WT, Vaughn P, et al. (2011) SIGNATURE: a workbench for gene expression signature analysis. BMC Bioinformatics 12: 443. doi: 10.1186/1471-2105-12-443 22078435

58. Gatza ML, Kung HN, Blackwell KL, Dewhirst MW, Marks JR, et al. (2011) Analysis of tumor environmental response and oncogenic pathway activation identifies distinct basal and luminal features in HER2-related breast tumor subtypes. Breast Cancer Res 13: R62. doi: 10.1186/bcr2899 21672245

59. Chen JL, Merl D, Peterson CW, Wu J, Liu PY, et al. (2010) Lactic acidosis triggers starvation response with paradoxical induction of TXNIP through MondoA. PLoS Genet 6.

60. Lucas JE, Kung HN, Chi JT (2010) Latent factor analysis to discover pathway-associated putative segmental aneuploidies in human cancers. PLoS Comput Biol 6: e1000920. doi: 10.1371/journal.pcbi.1000920 20824128

61. Chin K, DeVries S, Fridlyand J, Spellman PT, Roydasgupta R, et al. (2006) Genomic and transcriptional aberrations linked to breast cancer pathophysiologies. Cancer Cell 10: 529–541. 17157792

62. Tennant DA, Frezza C, MacKenzie ED, Nguyen QD, Zheng L, et al. (2009) Reactivating HIF prolyl hydroxylases under hypoxia results in metabolic catastrophe and cell death. Oncogene 28: 4009–4021. doi: 10.1038/onc.2009.250 19718054

63. Currie E, Schulze A, Zechner R, Walther TC, Farese RV Jr (2013) Cellular Fatty Acid Metabolism and Cancer. Cell Metabolism 18: 153–161. doi: 10.1016/j.cmet.2013.05.017 23791484

64. Migita T, Narita T, Nomura K, Miyagi E, Inazuka F, et al. (2008) ATP citrate lyase: activation and therapeutic implications in non-small cell lung cancer. Cancer Res 68: 8547–8554. doi: 10.1158/0008-5472.CAN-08-1235 18922930

65. Metallo CM, Gameiro PA, Bell EL, Mattaini KR, Yang J, et al. (2012) Reductive glutamine metabolism by IDH1 mediates lipogenesis under hypoxia. Nature 481: 380–384.

66. Wollenick K, Hu J, Kristiansen G, Schraml P, Rehrauer H, et al. (2012) Synthetic transactivation screening reveals ETV4 as broad coactivator of hypoxia-inducible factor signaling. Nucleic Acids Research 40: 1928–1943. doi: 10.1093/nar/gkr978 22075993

67. Wu M-Z, Tsai Y-P, Yang M-H, Huang C-H, Chang S-Y, et al. Interplay between HDAC3 and WDR5 Is Essential for Hypoxia-Induced Epithelial-Mesenchymal Transition. Molecular Cell 43: 811–822. doi: 10.1016/j.molcel.2011.07.012 21884981

68. Beyer S, Kristensen MM, Jensen KS, Johansen JV, Staller P (2008) The histone demethylases JMJD1A and JMJD2B are transcriptional targets of hypoxia-inducible factor HIF. J Biol Chem 283: 36542–36552. doi: 10.1074/jbc.M804578200 18984585

69. Johansson C, Tumber A, Che K, Cain P, Nowak R, et al. (2014) The roles of Jumonji-type oxygenases in human disease. Epigenomics 6: 89–120. doi: 10.2217/epi.13.79 24579949

70. Blair LP, Cao J, Zou MR, Sayegh J, Yan Q (2011) Epigenetic Regulation by Lysine Demethylase 5 (KDM5) Enzymes in Cancer. Cancers (Basel) 3: 1383–1404.

71. Yang J, Ledaki I, Turley H, Gatter KC, Montero JC, et al. (2009) Role of hypoxia-inducible factors in epigenetic regulation via histone demethylases. Ann N Y Acad Sci 1177: 185–197. doi: 10.1111/j.1749-6632.2009.05027.x 19845621

72. Shmakova A, Batie M, Druker J, Rocha S (2014) Chromatin and oxygen sensing in the context of JmjC histone demethylases. Biochem J 462: 385–395. doi: 10.1042/BJ20140754 25145438

73. Mimura I, Tanaka T, Wada Y, Kodama T, Nangaku M (2011) Pathophysiological response to hypoxia—from the molecular mechanisms of malady to drug discovery: epigenetic regulation of the hypoxic response via hypoxia-inducible factor and histone modifying enzymes. J Pharmacol Sci 115: 453–458. 21422728

74. Mannironi C, Proietto M, Bufalieri F, Cundari E, Alagia A, et al. (2014) An high-throughput in vivo screening system to select H3K4-specific histone demethylase inhibitors. PLoS One 9: e86002. doi: 10.1371/journal.pone.0086002 24489688

75. Sayegh J, Cao J, Zou MR, Morales A, Blair LP, et al. (2013) Identification of small molecule inhibitors of Jumonji AT-rich interactive domain 1B (JARID1B) histone demethylase by a sensitive high throughput screen. J Biol Chem 288: 9408–9417. doi: 10.1074/jbc.M112.419861 23408432

76. Zaidi N, Royaux I, Swinnen JV, Smans K (2012) ATP Citrate Lyase Knockdown Induces Growth Arrest and Apoptosis through Different Cell- and Environment-Dependent Mechanisms. Molecular Cancer Therapeutics 11: 1925–1935. doi: 10.1158/1535-7163.MCT-12-0095 22718913

77. Bhalla K, Hwang BJ, Dewi RE, Ou L, Twaddel W, et al. (2011) PGC1alpha promotes tumor growth by inducing gene expression programs supporting lipogenesis. Cancer Res 71: 6888–6898. doi: 10.1158/0008-5472.CAN-11-1011 21914785

78. Mason P, Liang B, Li L, Fremgen T, Murphy E, et al. (2012) SCD1 Inhibition Causes Cancer Cell Death by Depleting Mono-Unsaturated Fatty Acids. PLoS ONE 7: e33823. doi: 10.1371/journal.pone.0033823 22457791

79. Zhan Y, Ginanni N, Tota MR, Wu M, Bays NW, et al. (2008) Control of Cell Growth and Survival by Enzymes of the Fatty Acid Synthesis Pathway in HCT–116 Colon Cancer Cells. Clinical Cancer Research 14: 5735–5742. doi: 10.1158/1078-0432.CCR-07-5074 18794082

80. Brusselmans K, De Schrijver E, Verhoeven G, Swinnen JV (2005) RNA interference-mediated silencing of the acetyl-CoA-carboxylase-alpha gene induces growth inhibition and apoptosis of prostate cancer cells. Cancer Res 65: 6719–6725. 16061653

81. Chajes V, Cambot M, Moreau K, Lenoir GM, Joulin V (2006) Acetyl-CoA carboxylase alpha is essential to breast cancer cell survival. Cancer Res 66: 5287–5294. 16707454

82. Aytes A, Mitrofanova A, Kinkade CW, Lefebvre C, Lei M, et al. (2013) ETV4 promotes metastasis in response to activation of PI3-kinase and Ras signaling in a mouse model of advanced prostate cancer. Proceedings of the National Academy of Sciences 110: E3506–E3515.

83. Pellecchia A, Pescucci C, De Lorenzo E, Luceri C, Passaro N, et al. (2012) Overexpression of ETV4 is oncogenic in prostate cells through promotion of both cell proliferation and epithelial to mesenchymal transition. Oncogenesis 1: e20. doi: 10.1038/oncsis.2012.20 23552736

84. Keld R, Guo B, Downey P, Cummins R, Gulmann C, et al. (2011) PEA3/ETV4-related transcription factors coupled with active ERK signalling are associated with poor prognosis in gastric adenocarcinoma. Br J Cancer 105: 124–130. doi: 10.1038/bjc.2011.187 21673681

85. Yuen H-F, Chan Y-K, Grills C, McCrudden CM, Gunasekharan V, et al. (2011) Polyomavirus enhancer activator 3 protein promotes breast cancer metastatic progression through Snail-induced epithelial—mesenchymal transition. The Journal of Pathology 224: 78–89. doi: 10.1002/path.2859 21404275

86. Schmaltz C, Hardenbergh PH, Wells A, Fisher DE (1998) Regulation of Proliferation-Survival Decisions during Tumor Cell Hypoxia. Molecular and Cellular Biology 18: 2845–2854. 9566903

87. Brunelle JK, Santore MT, Budinger GRS, Tang Y, Barrett TA, et al. (2004) c-Myc Sensitization to Oxygen Deprivation-induced Cell Death Is Dependent on Bax/Bak, but Is Independent of p53 and Hypoxia-inducible Factor–1. Journal of Biological Chemistry 279: 4305–4312. 14627695

88. Conacci-Sorrell M, Ngouenet C, Anderson S, Brabletz T, Eisenman RN (2014) Stress-induced cleavage of Myc promotes cancer cell survival. Genes & Development 28: 689–707.

89. Wong WJ, Qiu B, Nakazawa MS, Qing G, Simon MC (2013) MYC Degradation under Low O2 Tension Promotes Survival by Evading Hypoxia-Induced Cell Death. Molecular and Cellular Biology 33: 3494–3504. doi: 10.1128/MCB.00853-12 23816886

90. Guimaraes-Camboa N, Stowe J, Aneas I, Sakabe N, Cattaneo P, et al. (2015) HIF1alpha Represses Cell Stress Pathways to Allow Proliferation of Hypoxic Fetal Cardiomyocytes. Dev Cell 33: 507–521. doi: 10.1016/j.devcel.2015.04.021 26028220

91. Guimarães-Camboa N, Stowe J, Aneas I, Sakabe N, Cattaneo P, et al. HIF1α Represses Cell Stress Pathways to Allow Proliferation of Hypoxic Fetal Cardiomyocytes. Developmental Cell 33: 507–521. doi: 10.1016/j.devcel.2015.04.021 26028220

92. Suzuki K, Bose P, Leong-Quong RY, Fujita DJ, Riabowol K (2010) REAP: A two minute cell fractionation method. BMC Res Notes 3: 294. doi: 10.1186/1756-0500-3-294 21067583

93. Ferrara CT, Wang P, Neto EC, Stevens RD, Bain JR, et al. (2008) Genetic Networks of Liver Metabolism Revealed by Integration of Metabolic and Transcriptional Profiling. PLoS Genet 4: e1000034. doi: 10.1371/journal.pgen.1000034 18369453

94. An J, Muoio DM, Shiota M, Fujimoto Y, Cline GW, et al. (2004) Hepatic expression of malonyl-CoA decarboxylase reverses muscle, liver and whole-animal insulin resistance. Nat Med 10: 268–274. 14770177

95. Lamonte G, Tang X, Chen JL, Wu J, Ding CK, et al. (2013) Acidosis induces reprogramming of cellular metabolism to mitigate oxidative stress. Cancer Metab 1: 23. doi: 10.1186/2049-3002-1-23 24359630

Genetika Reprodukční medicína

Článek vyšel v časopise

PLOS Genetics

2015 Číslo 10

Nejčtenější v tomto čísle
Kurzy Podcasty Doporučená témata Časopisy
Zapomenuté heslo

Nemáte účet?  Registrujte se

Zapomenuté heslo

Zadejte e-mailovou adresu, se kterou jste vytvářel(a) účet, budou Vám na ni zaslány informace k nastavení nového hesla.


Nemáte účet?  Registrujte se