HIF-1-alpha links mitochondrial perturbation to the dynamic acquisition of breast cancer tumorigenicity
Metrics: HTML 2244 views | ?
Ching-Ying Kuo1, Chun-Ting Cheng1,2, Peifeng Hou1,3,4, Yi-Pei Lin5, Huimin Ma1, Yiyin Chung1, Kevin Chi1, Yuan Chen6, Wei Li7, Hsing-Jien Kung5,8, David K. Ann1,2
1Diabetes and Metabolism Research Institute, City of Hope, Duarte, CA, USA
2Irell & Manella Graduate School of Biological Sciences, City of Hope, Duarte, CA, USA
3Department of Oncology, Fujian Medical University Union Hospital, Fuzhou, Fujian, China
4Fujian Key Laboratory of Translational Cancer Medicine, Fuzhou, Fujian, China
5Integrated Laboratory, Center of Translational Medicine, Taipei Medical University, Taipei, Taiwan, ROC
6Department of Molecular Medicine, Beckman Research Institute, City of Hope, Duarte, CA, USA
7Department of Dermatology, Keck School of Medicine, University of Southern California, Los Angeles, CA, USA
8Institute of Molecular and Genomic Medicine, National Health Research Institutes, Taiwan, ROC
David K. Ann, e-mail: firstname.lastname@example.org
Keywords: α-ketoglutarate, mitochondria, metabolic reprogramming, hypoxia-inducible factor-1α (HIF-1α), breast cancer tumorigenicity
Received: November 19, 2015 Accepted: March 10, 2016 Published: April 4, 2016
Up-regulation of hypoxia-inducible factor-1α (HIF-1α), even in normoxia, is a common feature of solid malignancies. However, the mechanisms of increased HIF-1α abundance, and its role in regulating breast cancer plasticity are not fully understood. We have previously demonstrated that dimethyl-2-ketoglutarate (DKG), a widely used cell membrane-permeable α-ketoglutarate (α-KG) analogue, transiently stabilizes HIF-1α by inhibiting prolyl hydroxylase 2. Here, we report that breast cancer tumorigenicity can be acquired through prolonged treatment with DKG. Our results indicate that, in response to prolonged DKG treatment, mitochondrial respiration becomes uncoupled, leading to the accumulation of succinate and fumarate in breast cancer cells. Further, we found that an early increase in the oxygen flux rate was accompanied by a delayed enhancement of glycolysis. Together, our results indicate that these events trigger a dynamic enrichment for cells with pluripotent/stem-like cell markers and tumorsphere-forming capacity. Moreover, DKG-mediated metabolic reprogramming results in HIF-1α induction and reductive carboxylation pathway activation. Both HIF-1α accumulation and the tumor-promoting metabolic state are required for DKG-promoted tumor repopulation capacity in vivo. Our data suggest that mitochondrial adaptation to DKG elevates the ratio of succinate or fumarate to α-KG, which in turn stabilizes HIF-1α and reprograms breast cancer cells into a stem-like state. Therefore, our results demonstrate that metabolic regulation, with succinate and/or fumarate accumulation, governs the dynamic transition of breast cancer tumorigenic states and we suggest that HIF-1α is indispensable for breast cancer tumorigenicity.
Altered energy metabolism, a hallmark of cancer, plays a key role in tumor development for rapid proliferation, maintenance of redox homeostasis, and epigenetic reprogramming [1, 2]. Experimental evidence has supported the hypothesis that various metabolic enzymes are reprogrammed by oncogenes and dysfunctional or non-functional tumor suppressor genes [3, 4]. Consequently, the metabolic profiles in proliferating cancer cells are shifted toward glycolysis, even though there is sufficient oxygen available to support oxidative phosphorylation (OXPHOS); this shift to glycolysis is also known as the Warburg effect [5, 6]. Recently, attention has turned to the mitochondria and their role in promoting stem cell self-renewal and tumorigenic cell expansion [7–10]. Implicit in the ability of cancer cells to switch metabolic states is the concept that alternating between differing metabolic states is required to effectively adapt to the rapidly changing tumor nutritional microenvironment.
There is increasing awareness that individual breast cancer (BC) tumors contain a heterogeneous cell population, and that this intratumoral heterogeneity contributes to therapy failure and disease progression . In the past decade, studies have largely focused on the unidirectional hierarchy from tumorigenic cells, so-called cancer stem cells (CSCs) to non-tumorigenic progeny, emphasizing the critical role of CSCs in cancer development, treatment resistance and metastasis . Recently, however, it has been proposed that non-stem-like cancer cells can spontaneously acquire or preserve CSC-like properties [12–14]. Although emerging studies reveal that CSCs are metabolically distinct from cancer cells, no studies have explored the tumor nutritional microenvironment and its effects on the dynamics of CSC-like subpopulations.
Hypoxia-inducible factor-1α (HIF-1α) is a master transcriptional regulator that responds to a hypoxic environment and induces the transcription of genes, including vascular endothelial growth factor (VEGF), glucose transporter 1 (GLUT1) and pyruvate dehydrogenase kinase 1 (PDK1), which promote vascularization and tumor cell survival . For example, overexpression of HIF-1α is associated with BC metastasis and poor survival of BC patients [16–19]. In some cancers, accumulation of HIF-1α can result from loss-of-function mutations in proteins that cause HIF-1α degradation, such as the Von Hippel-Lindau tumor suppressor (pVHL) or the enzymes that produce cofactors for the prolyl hydroxylases (PHDs). For example, sporadic renal cell carcinomas with mutations or deletions in the VHL gene exhibit elevated HIF-1α levels . In addition, mutations in succinate dehydrogenase (SDH) and fumarate hydratase (FH), enzymes that produce competitive metabolites for PHD cofactors, are found in cancers [21–30]. SDH and FH hydrolyze succinate and fumarate, respectively, to fuel the tricarboxylic acid (TCA) cycle. Mutations in SDH or FH cause succinate or fumarate to accumulate and compete with α-ketoglutarate (α-KG) for PHD binding, thereby inhibiting PHD and stabilizing HIF-1α [31, 32]. Mutations have also been identified in isocitrate dehydrogenase 1 (IDH1) that inhibit IDH1 catalytic activity in gliomas, thereby reducing the production of α-KG, inhibiting PHD, increasing HIF-1α, and presumably, promoting tumorigenesis .
Although the mechanism is not totally understood, some evidence suggests that α-KG can increase the stem or stem-like potential of embryonic stem cells (ESCs) . Here, we have addressed this fundamental biological question in the context of BC cell metabolic state. Our laboratory initially identified that dimethyl-2-ketoglutarate (DKG), which has been widely used as an α-KG-supplement [35, 36], transiently stabilizes HIF-1α by inhibiting PHD2-mediated hydroxylation/degradation of HIF-1α under normoxia . HIF-1α, along with its complex signaling network, has been proposed as a key mediator of BC malignancies [16, 38]. Nonetheless, nothing is known about the mechanism of DKG-induced PHD2 inhibition and the consequences of prolonged DKG exposure on BC cells. Here, we studied the CSC-like properties of a panel of established and patient-derived BC cells treated with DKG. The metabolic and transcriptional landscape and the underlying mechanism were analyzed. We found that sustained DKG treatment triggered the accumulation of succinate and fumarate, while reducing the abundance of mRNAs encoding SDH, FH, and subunits of the mitochondrial electron transport chain (ETC) complex I and V. Our data suggest that differential regulation of mitochondrial respiration, glycolysis and fatty acid oxidation (FAO), coupled with accumulated HIF-1α, aggravate tumorigenicity in vivo. To our knowledge, this is the first use of a metabolite precursor, DKG, to perturb mitochondrial function, which in turn regulates BC tumorigenicity. Together, our results provide insight into how the tumor metabolic state is dynamically regulated to promote breast cancer aggressiveness.
DKG causes HIF-1α induction, succinate and fumarate accumulation, and transcriptional reprogramming in BC cells
To extend our previous observations , multiple BC cell lines were treated with a pharmacological dose of DKG (10 mM). Consistent with our initial results, DKG induced HIF-1α in four established BC cell lines (Figure 1A, upper panel) and two primary BC cell lines (Figure 1A, lower panel). The induction of HIF-1α was rapid (within 2 h after treatment [Figure 1B, upper panel]) and transient (Figure 1B, lower panel). In the remaining studies, BC cells were replenished with DKG daily to sustain HIF-1α at a higher level. To understand intracellular action of DKG, we measured several key intermediary metabolites in the TCA cycle using NMR-based metabolomic analysis. The intracellular levels of α-KG, succinate, fumarate, citrate and pyruvate all increased significantly, albeit with different kinetics and to varying degrees, within 24 h after DKG treatment (Figures 1C and S1A). Notably, although we treated BC cells with a pharmacological dose of DKG (10 mM), the intracellular level of α-KG remained in the physiological range (~μM) (Figure S1A), excluding the possibility that the effects we observed were due to excess intracellular α-KG. The disproportionate increase in these metabolites elevated the ratios of succinate to α-KG (from 0.54 to 1.02, 1.14 and 0.99) and fumarate to α-KG (from 0.11 to 0.32, 0.36 and 0.38) in a time-dependent manner (Figure S1B). Because we also observed a time-dependent ATP reduction (Figure 1D), we speculated that the DKG-fueled TCA cycle was uncoupled from ATP generation. Because succinate and fumarate are known to inhibit PHD2 [28, 31, 32], their accumulation in the DKG-treated cells corroborated our previous observations that DKG stabilized HIF-1α through inhibiting PHD2-mediated hydroxylation .
Figure 1: Dimethyl-2-ketoglutarate (DKG) induces HIF-1α and alters the metabolic profile in breast cancer (BC) cells. A. DKG stabilizes HIF-1α in multiple BC cells. HIF-1α levels in the DKG-treated BC cells (10 mM, 2 and 6 h) were assessed by western blot analysis. B. The induction of HIF-1α by DKG is fast and transient. HIF-1α levels in DKG-treated MDA-MB-231 cells were assessed by western blot analysis. C, D. Succinate and fumarate accumulate (C) and ATP is decreased (D) in DKG-treated cells. Metabolites in the MDA-MB-231 cells treated with DKG (10 mM, 0, 2, 8, 24 h) were analyzed by NMR spectroscopy. n = 3, *: p < 0.05; **: p < 0.01; ***: p < 0.005. (A, B) One representative blot from n = 3 is shown. Italic numbers indicate the relative protein level.
Because HIF-1α is known to regulate transcription, we therefore compared the gene expression profiles in MDA-MB-231 cells and two primary BC cells with or without DKG administration by performing RNA-sequencing (RNA-seq) analysis. The top five DKG-affected pathways were HIF-1α signaling, ubiquinol-10 biosynthesis, cell cycle control, chromosomal replication and TGF-β signaling (Figure S1C). We concluded that DKG treatment, in addition to inducing HIF-1α (Figure 1A), creates a pseudohypoxic state under normoxia. From our RNA-seq analysis, we also observed that the message abundance of SDH and FH was down-regulated in the DKG-treated cells (Figure S1D). We further postulated that the increase in both succinate and fumarate, as well as the decrease in SDH and FH mRNA levels, resulted in an imbalance of TCA metabolites. This metabolite imbalance could then impair PHD2 activity, thereby stabilizing HIF-1α and reprogramming the transcriptional landscape in BC cells.
DKG promotes the acquisition of breast cancer stem cell-like properties
HIF-1α signaling has been proposed to be a key mediator of BC malignancies [16, 38]; we therefore investigated the effects of prolonged DKG treatment on the tumorigenic properties of BC cells. Prolonged treatment with DKG (10 days) reduced the clonogenicity of MDA-MB-231 cells (Figure S1E, upper panel). However, after terminating DKG treatment, DKG pre-treated cells resumed proliferation and exhibited comparable colony numbers to untreated cells (Figure S1E, lower panel), ruling out the hypothesis that DKG is toxic and suggesting that its effect is dynamic and reversible. In parallel, the subset of Ki67Low/Hoechst 33342Low cells, representing the G0 quiescent subpopulation, increased in DKG-treated MDA-MB-231 cells (Figure S1F, upper panel, lower-left quadrants). After 4 to 7 days of DKG treatment, the G0 subpopulation increased by approximately 2-fold (12.7% to 24.5% and 31.3%) (Figure S1F, lower panel).
Given that tumorigenic cells can grow as 3-dimensional non-adherent structures when clonally seeded in serum-free media, we pre-treated a panel of established and primary BC cells with DKG (10 mM) for 7 days and subjected the cells to tumorsphere formation assays. Consistently, DKG pre-treatment increased both the size (Figure S2A) and the number (Figure 2A, left panel) of tumorspheres in multiple BC cells. In addition, DKG pre-treatment supported the in vitro propagation of tumorspheres (Figure 2A, right panel), suggesting that DKG enriches the stem-like subpopulation in culture.
Figure 2: DKG promotes tumorigenic properties in BC cells. A. DKG pre-treatment promotes tumorsphere formation. Left panel: Quantification of tumorspheres formed by the untreated and DKG pre-treated BC cells. Right panel: In vitro serial passaging of tumorspheres formed by the untreated and DKG-treated MCF7 cells. *: p < 0.05; **: p < 0.01, n = 3. B. DKG regulates the abundance of cancer stem cell (CSC) surface markers in BC cells. Flow cytometric analyses of surface markers in DKG-treated BC cells (10 mM, 4, 7 days). CD133 was assessed in MDA-MB-231 cells (a). CD44 and CD24 were assessed in MCF7 (b), MDA-MB-468 (c) and primary BC cells (d). The percentage of CD133-positive or CD44HighCD24Low subpopulations in the untreated sample was set as 1. Bar graphs represent the mean ± SD, n = 3. C. DKG converts non-tumorigenic subpopulations to tumorigenic subpopulations. MDA-MB-468 cells were sorted based on CD44 and CD24 expression. Sorted cells were treated with DKG (10 mM, 7 days). CD44 and CD24 expression was assessed. APC: allophycocyanin-conjugated. PE: phycoerythin-conjugated. Representative graphs. n = 3. D. DKG elevates OCT4 abundance. Upper panel: OCT4 levels in multiple BC cell lines treated with DKG (10 mM for 4 or 7 days) were determined by western blot analysis. Lower panel: OCT4 is expressed in the CD44High subpopulations after DKG treatment. The DKG-treated MDA-MB-468 cells (10 mM, 7 days) were sorted into CD44HighCD24Low (HL), CD44HighCD24High (HH) and CD44Low (L) subpopulations for analyzing OCT4 expression by western blot analysis. Representative blots. n = 3. E. DKG-induced OCT4 expression depends on HIF-1α. F. HIF-1α is responsible for inducing a tumorigenic subpopulation in BC cells. **: p < 0.01, n = 3. (D, E) One representative blot from n = 3 is shown. Italic numbers indicate the relative protein level. ND: not detectable.
To further characterize the tumorsphere-promoting effects of DKG, flow cytometric analyses were performed to examine the abundance of CD44, CD24 and CD133. These are cell surface markers frequently used to identify CSCs in BC [39–42]. Because more than 70% of MDA-MB-231 cells exhibit the CD44HighCD24Low phenotype  and CD133 expression is highly correlated with the invasiveness and poor prognosis of triple-negative BC (TNBC) [40, 42], we assessed whether DKG treatment increased the surface expression of CD133 in MDA-MB-231 cells by flow cytometry. A 4-fold increase in the percentage of CD133-positive was observed in DKG-treated MDA-MB-231 cells (4 and 7 days; Figures 2B, panel a, and S2B). In parallel, the CD44HighCD24Low subpopulation, within the bulk of MCF7 cells, was elevated from 1.8% to 6.4% and 9.5% following 4 and 7 days of DKG treatment (Figures 2B, panel b, and S2C). Likewise, the CD44HighCD24Low subpopulation also increased in the DKG-treated MDA-MB-468 cells, by approximately 3-fold (Figures 2B, panel c, and S2D). In addition, DKG treatment caused a significant enrichment of the CD44HighCD24Low subpopulation in two patient-derived primary BC cell lines, BC_01 and BC_03 (Figures 2B, panel d, and S2E). To differentiate whether DKG-mediated modulation in cell surface marker expression was stable or transient, we performed fluorescence-activated cell sorting (FACS) to enrich for four subpopulations, according to the surface expression of CD44 and CD24 in MDA-MB-468 cells, and we then incubated the respective subpopulations with DKG for an additional 7 days. We found that DKG not only expanded the CD44HighCD24Low subpopulation, but also converted the other three non-CD44HighCD24Low subpopulations to CD44HighCD24Low (Figure 2C). These results suggest that DKG treatment enriched for the subpopulation of BC cells with increased tumorsphere formation capacity and increased CSC surface marker abundance. These two traits are frequently used as surrogate markers for tumorigenic CSC-like properties.
Next, we examined the expression of pluripotency factors, which were reported to express at higher levels in breast CSC-like cells and are associated with BC progression [43–47]. For example, the subpopulation of BC cells expressing a high level of octamer-binding transcription factor 4 (OCT4) displays tumorigenic properties [43, 48, 49]. After DKG treatment for 4 and 7 days, the steady-state OCT4 protein abundance was elevated in multiple BC cell lines (Figure 2D, upper panel). Notably, OCT4 was enriched in the sorted CD44High subpopulations (Figure 2D, lower panel). In parallel, the mRNA abundances of pluripotency factors, including OCT4, NANOG, SOX2, MYC and KLF4, as well as the cell surface marker CD44, increased in DKG-treated MDA-MB-231 and MCF7 cells (Figure S2F). Altogether, these results demonstrated that DKG induced the expression of surrogate markers for the CSC-like phenotype in multiple established and primary BC cells.
Lastly, we assessed whether the up-regulation of the surrogate markers for the tumorigenic phenotype was dependent on HIF-1α. Knocking down HIF-1α by a short hairpin (sh)RNA reduced the induction of OCT4 and CSC surface markers in DKG-treated cells (Figures 2E and 2F), suggesting the DKG-induced expression of OCT4 and the CSC surface markers was HIF-1α-dependent. Taken together, these data suggest that DKG treatment promotes the accumulation of succinate and fumarate, resulting in HIF-1α induction and changes in subsequent transcriptional reprogramming and thereby impacting multiple cellular processes in BC cells.
DKG enhances tumorigenic capacity
Mouse xenografts assay was used to confirm the effects of DKG on promoting the tumorigenic properties of MDA-MB-231 cells in vivo. Serial dilutions (104, 103, 102 cells) of untreated and DKG pre-treated MDA-MB-231 cells were injected in the left and right flanks of the same NOD/SCID/IL-2Rγ null (NSG) mouse, respectively. Visible tumors were formed on the both flanks in all groups of mice (n = 8 for each dilution). Consistent with results shown in Figure 2, DKG pre-treated cells gave rise to bigger tumors than the vehicle pre-treated cells (Figure 3A). Next, to determine the contribution of HIF-1α in DKG-mediated tumor propagation, MDA-MB-231/shCON or/shHIF-1α cells (104) were pre-treated with DKG, and then injected into the flanks of NSG mice. For comparison, we also injected 106 vehicle pre-treated MDA-MB-231/shCON cells. Measurable tumors formed in the mice injected with 106 MDA-MB-231/shCON cells by day 30 after engrafting. Strikingly, after 40-days of engraftment, 104 DKG-pretreated MDA-MB-231/shCON cells propagated rapidly to form tumors similar in size to the tumors formed by 106 MDA-MB-231/shCON cells. In contrast, neither 104 untreated MDA-MB-231/shCON cells nor 104 DKG-pretreated MDA-MB-231/shHIF-1α cells formed measurable tumors (Figure 3B). Tumor sizes and weights at day 65 after engraftment are shown in Figure 3C. 104 DKG-pretreated MDA-MB-231/shCON cells produced tumors with sizes comparable to those from 106 untreated MDA-MB-231/shCON cells. Next, we examined a well-established HIF-1α downstream protein, carbonic anhydrase 9 (CAIX), in the xenograft tumors by immunohistochemistry (IHC) and found a significant increase in CAIX signals in the DKG-pretreated tumors (Figure 3D). As expected, the level of CAIX was much lower in the HIF-1α knockdown tumors. This result suggests that DKG increases CAIX abundance via HIF-1α signaling. Consistent with our in vitro study, we observed an increase of CD44 and CD133 in DKG-pretreated tumors and HIF-1α knockdown reduced the intensity of CD44 and CD133. Moreover, we also observed more vascularization in DKG-pretreated tumors in a HIF-1α-dependent manner, consistently indicating an increase in HIF-1α signal in the context of DKG treatment (Figure 3D).
Figure 3: DKG augments the tumorigenic properties of BC cells in vivo. A. DKG enhances the growth of tumors in vivo. The untreated or DKG pre-treated MDA-MB-231 cells were injected subcutaneously in the left and right flank of the same NSG mouse in serial dilutions (104, 103, 102 cells). The plot shows the comparison of corresponding tumor weight (g), measured after harvesting on day 60, in the same mouse. *: p < 0.05; ***: p < 0.005, n = 8. B, C. DKG augments the tumorigenic properties of BC cells through HIF-1α in vivo. (B) Different numbers of MDA-MB-231/shCON or shHIF-1α cells with or without pre-treatment with DKG were injected into NSG mice. Tumor volume was measured twice weekly for 65 days and tumors were harvested on day 65. (C) Images show that pre-treatment with DKG boosted expansion of xenografted MDA-MB-231 cells in a HIF-1α-dependent manner (upper panel). Tumor weight (g) was measured after harvesting on day 65 and is shown as mean ± SD (lower panel). *: p < 0.05. n = 4. D. Pre-treatment with DKG increases CAIX abundance in xenografted tumors. Immunohistochemical (IHC) staining of CAIX, CD44, CD133 and H&E staining of the tumors. Bar: 50 μm. Insets show the enlargement of the indicated areas. E. DKG confers relative resistance of BC cells to irradiation. DKG-pre-treated xenografted MDA-MB-231/shCON and MDA-MB-231/shHIF-1α cells were subjected to a clonogenic survival assay followed by irradiation. Colonies were counted, corrected for plating efficiency and presented as the surviving fraction relative to untreated controls. *: p < 0.05; **: p < 0.01, shCON/DKG compared with shCON; #: p < 0.05, shCON/DKG compared with shHIF-1α/DKG. F. Harvested tumor cells from xenografted tumors are still responsive to DKG. HIF-1α and OCT4 levels in MDA-MB-231 cells harvested from xenografted tumors were assessed by western blot analysis. One representative blot from n = 3 is shown. Italic numbers indicate the relative protein level. ND: not detectable.
Another important breast CSC-like feature is treatment resistance. To address this possibility, we established cell lines from the harvested tumors, treated them with DKG for 7 days, and performed a clonogenic survival assay to evaluate cell survival following ionizing radiation (IR). As shown in Figure 3E, the DKG-treated tumor cells were relatively resistant to IR, and this phenomenon was reversed in the HIF-1α knockdown tumor cells. The observation of a lack of OCT4 induction in DKG-treated cells derived from MDA-MB-231/shHIF-1α-xenografted tumors (Figure 3F) supported the notion that both HIF-1α and DKG-mediated metabolic adaptation are required for the expansion of cells with tumorigenic properties.
DKG increases oxygen consumption rate
To investigate how DKG impacts metabolic pathways, we measured the metabolic flux in the DKG-treated BC cells. Intriguingly, DKG treatment first preferentially increased the oxygen consumption rate (OCR) (day 2), and then mainly elevated the extracellular acidification rate (ECAR) (day 4), with both values reaching a plateau by day 4 (Figures 4A, upper panel, and S3A). The enhancement of both OCR and ECAR was diminished by HIF-1α knockdown (Figures 4A, lower panel, S3B). Interestingly, this increase in basal OCR was not abrogated by oligomycin treatment (Figures S3C and S3D, panel a), suggesting an increase in the oligomycin-insensitive fraction, which could be caused by proton leakage from the intermembrane space . However, additional analysis indicated that the increased proton leakage in DKG-treated cells was very modest (Figure S3D, panel b). Instead, a marked expansion of non-mitochondrial respiration was detected in DKG-treated cells (Figure S3D, panel c). In addition, an increase in the maximal respiration (uncoupled respiration) was observed in DKG-treated cells (Figure S3D, panel d). We therefore concluded that one of the DKG-associated metabolic phenotypes was the elevated respiration in the setting of mitochondrial uncoupling, coexisting with reduced ATP (Figure 1D). Moreover, the increase in uncoupled oxygen utilization was HIF-1α-dependent, as knockdown of HIF-1α significantly dampened the DKG-induced basal, non-mitochondrial and maximal respiration (Figures S3C and S3D, panels a, c, d).
To explore the underlying mechanism(s), we next evaluated the mRNA abundance of glycolytic genes from our RNA-seq data. Consistent with Figure 4A, the message abundance of glycolytic genes was up-regulated in DKG-treated MDA-MB-231 and two primary BC cells (Figure 4B, upper panel). The abundance of the GLUT1 and PDK1 messages was further validated by qRT-PCR in MDA-MB-231 and MCF7 cells treated with DKG (Figure 4B, lower panels). In contrast, relative message abundances for mitochondrial ETC proteins, especially complexes I and V, decreased after DKG treatment (Figure 4C, left panel). Indeed, unlike its role in restoring steady-state ATP level in the cells deprived of glucose (Glc) for 48h (Figure 3C, upper right panel), but consistent with Figure 1D, prolonged DKG treatment reduced ATP levels in the cells cultured in nutrient-rich conditions (Figure 4C, lower right panel). It is also possible that increased α-KG (Figure 1C) binds and inhibits ATP synthase to reduce ATP production . Lastly, given that DKG increased OCR, which could serve as a source of intracellular ROS, we examined the ROS levels in the DKG-treated cells by 2’,7’-dichlorofluorescin diacetate (DCFDA) labeling. After 2 and 4 days of DKG treatment, a modest, albeit significant, N-acetylcysteine (NAC)-sensitive increase in oxidized 2’,7’-dichlorofluorescin (DCF) signal was observed (Figure 4D). Taken together, our data support a model where DKG regulates specific sets of genes to increase both OCR and ECAR in BC cells, while uncoupling the ETC with ATP synthesis, preventing a massive surge in ROS.
Figure 4: DKG enhances metabolic flux in BC cells. A. Energy metabolism is enhanced by DKG-induced HIF-1α in MDA-MB-231 cells. The oxygen-consumption rate (OCR) and extracellular acidification rate (ECAR) were measured during the course of DKG treatment in MDA-MB-231 cells (upper panel) and in MDA-MB-231/shCON versus MDA-MB-231/shHIF-1α cells (lower panel). The PhenoGram profiles show the metabolic switch of OCR and ECAR at the sixth point (22.5 mins after Glc injection during the ECAR measurement) for the ECAR versus OCR plots, as shown in Figures S3A and S3B; n = 3. B. DKG increases the message abundance of genes regulating glucose metabolism. Upper panel: heat map generated from RNA-seq analysis of the untreated and DKG-treated MDA-MB-231 cells and primary BC cells. Lower panel: qRT-PCR analysis revealed that glucose transporter 1 (GLUT1) and pyruvate dehydrogenase kinase 1 (PDK1) mRNA levels are elevated in DKG-treated MDA-MB-231 cells (10 mM, 4, 7 days). *: p < 0.05; **: p < 0.01, n = 3. C. Genes regulating the mitochondrial electron transport chain are reduced by DKG treatment. Left panel: heat map generated from RNA-seq analysis. Right panel: MDA-MB-231 cells were cultured in Glc-depleted medium (upper right panel) or Glc-complete (25 mM) medium (lower right panel) supplemented with or without DKG for 48 h. The ATP levels were measured using a luciferase assay. *: p < 0.05; **: p < 0.01. n = 3. D. The levels of reactive oxygen species (ROS) increase in DKG-treated cells. ROS levels were measured by DCFDA staining. The oxidized DCF signal is shown. **: p < 0.01 (DKG-treated vs. untreated); #: p < 0.05 (DKG+NAC vs. DKG only), n = 3.
DKG increases reductive carboxylation and FAO beyond HIF-1α
From our NMR-based metabolite analysis, we also noticed that glutamine was decreased, whereas glutamate and citrate were increased, after DKG treatment (Figure 1C and 5A, upper panel). We postulated that, together with HIF-1α accumulation, the DKG-treated cells underwent reductive carboxylation, a process that utilizes glutamine to generate fatty acids. Indeed, the increased OCR in DKG-treated cells could be partially dampened by adding etomoxir (ETO), an inhibitor of carnitine palmitoyltransferase (CPT), which inhibits FAO  (Figure 5A, lower panel). We further speculated that lipids accumulated in DKG-treated cells. As expected, we detected elevated lipid droplet accumulation in DKG-treated MDA-MB-231 cells using Oil Red O staining (Figure 5B, upper panel) and a lipophilic fluorescent dye, BODIPY  (Figure 5B, lower panel). Notably, the magnitude of accumulated lipid droplets was partially dependent on the HIF-1α level. HIF-1α knockdown raised the basal lipid level and reduced the magnitude of induction (Figure 5C). Together with the results shown in Figure 5A, it is conceivable that DKG treatment also increases lipid utilization, presumably by FAO. To investigate a specific role of altered metabolism in the acquisition of tumorigenic properties, we pre-treated MDA-MB-231 cells with DKG for 4 days, added inhibitors of glycolysis (dichloroacetate, DCA) and FAO (ETO) for another 3 days, and assessed tumorsphere formation. As shown in Figure 5D, although treatment with DCA or ETO alone did not affect the tumorsphere number or other DKG-mediated phenotypes, including the CD133-positive or BODIPY-stained MDA-MB-231 fraction (Figure S4), combined DCA and ETO treatment notably reduced the tumorsphere number without affecting cell viability (data not shown).
Figure 5: DKG regulates glutamine-dependent reductive carboxylation. A. Upper panel: increased utilization of glutamine in DKG-treated cells. The levels of glutamine and glutamate in DKG-treated MDA-MB-231 cells were measured by NMR spectroscopy. **: p < 0.05; **: p < 0.01, n = 3. Lower panel: The enhanced OCR by DKG can be dampened by adding a fatty acid oxidation (FAO) inhibitor, etomoxir (ETO). OCR was measured 10 mins after ETO (40 nM) was added to MDA-MB-231 cells treated with DKG (10 mM, 4 days). n = 3. B. DKG promotes lipid droplet accumulation. Upper panel: Oil Red O staining was used to detect lipid droplets in DKG-treated MDA-MB-231 cells (10 mM, 4, 7 days). Bar: 50 μm. Representative images. n = 2. Lower panel: BODIPY 493/503 staining, followed by flow cytometric analysis, was performed to quantify the levels of neutral lipids in DKG-treated BC cells (10 mM, 7 days). **: p < 0.01, n = 3. C. The induction of lipid droplet accumulation by DKG is partially dependent on HIF-1α. BODIPY 493/503 staining was performed in DKG-treated BC cells harboring control shRNA or shRNA against HIF-1α. *: p < 0.05; **: p < 0.01; N.S.: not significant (DKG-treated vs. untreated). #: p < 0.05; ##: p < 0.01 (shCON vs. shHIF-1α), n = 3. D. Blocking both glycolysis and FAO reduces DKG-enhanced tumorsphere formation. MDA-MB-231 cells were pre-treated with DKG (10 mM, 7 days). Dichloroacetate (DCA) (10 mM) or ETO (10 μM) or DCA+ETO were added on day 5 and the cells were subjected to tumorsphere formation assays. *: p < 0.05.
To assess whether stabilization of HIF-1α alone was sufficient to mimic DKG and promote acquisition of tumorigenic phenotypes, dimethyloxaloylglycine (DMOG), an established PHD inhibitor, was used to stabilize HIF-1α. Unlike DKG, DMOG exclusively increased ECAR, at the expense of OCR (Figure S5A). Moreover, DMOG was less effective in enriching for the CD133-positive subpopulation (Figure S5B) and promoting lipid droplet accumulation (Figure S5C) than DKG, in MDA-MB-231 cells. However, DMOG and DKG induced HIF-1α and OCT4 abundance to a similar level in MDA-MB-231 and MCF7 cells (Figure S5D). Lastly, DKG rescued cells from death by long-term Glc deprivation (Figure S5E), but DMOG did not (Figure S5E), further distinguishing DKG from DMOG. Altogether, we concluded that the stabilization of HIF-1α is required, but alone is not sufficient, to establish the full suite of tumorigenic properties in BC cells. Therefore, we propose that increasing HIF-1α alone is not sufficient to promote tumorigenic phenotypes.
There is a growing appreciation for the idea that cancer cells undergo comprehensive metabolic reprogramming, leading to a glycolytic metabolic state that permits rapid growth and adaptation to the tumor microenvironment [1, 6, 54, 55]. Nonetheless, whereas altered glycolysis, glutaminolysis, and lipid synthesis have been demonstrated [6, 56], whether the exogenous nutrient availability alters the intrinsic tumor cell properties has not been clearly elucidated. Here, for the first time, we provide a comprehensive analysis of BC cell metabolic reprogramming, using DKG to uncouple mitochondrial respiration and also to increase succinate and fumarate. In response to DKG, the excess succinate/fumarate is then transported out of the mitochondria into the cytoplasm, where it impairs the activity of PHD, leading to the stabilization and activation of HIF-1α, and creating a state known as pseudohypoxia [32, 57]. This HIF-1α-induced pseudohypoxia then increases the glycolysis and oxygen flux rates. We propose that this HIF-1α-induced pseudohypoxia, together with increased FAO and transcriptional re-programming, promotes the acquisition of a tumorigenic CSC-like phenotype (Figure 6).
Figure 6: Model depicting the dual function of DKG to induce tumorigenic properties in BC cells through co-opting a pseudohypoxic pathway and metabolic rewiring. DKG promotes both OXPHOS and HIF-1α-dependent glycolysis to enhance tumorigenic expansion.
It is intriguing to note that the DKG-induced pseudohypoxic phenotype persisted for up to three passages (Figure 2A), which was further corroborated by an increased CAIX expression in tumors derived from DKG-pretreated BC cells (Figure 5D). Supported by the finding that CAIX is predominantly expressed in high-grade tumors , we propose that induction of a sustained pseudohypoxic state, with mitochondrial dysregulation, induces tumorigenicity.
The current studies reveal an important, but counterintuitive, ability of DKG to trigger a combined metabolic and transcriptional response, promoting a tumorigenic phenotype. DKG has been widely used as an analogue of α-KG, which is important for the TCA cycle, fatty acids and glutamine metabolism, and α-KG-dependent enzyme function. Therefore, one might suspect that an α-KG analogue would increase PHD enzyme function and de-stabilize HIF-1α. However, because DKG also increases the succinate and fumarate levels, the overall ratio of succinate/fumarate to α-KG remains high and PHD enzymes are inhibited, thereby increasing HIF-1α levels.
Although our in vitro and in vivo findings have established that increased HIF-1α levels, generated by DKG, are indispensable for the DKG-promoted tumor expansion, we cannot rule out the potential contribution of other α-KG-dependent enzymes. For instance, a group of jumonji domain-containing histone lysine demethylases (KDMs) utilizes α-KG as a cofactor for histone demethylation. α-KG is also important for DNA demethylation by the ten-eleven translocation (TET) enzymes . The potential role of these and other DKG targets will require further investigation.
Our results are in marked contrast to a recent report indicating that a low intracellular ratio of succinate to α-KG is important to maintain the proliferation and the pluripotency of mouse embryonic stem cells in glutamine-free media containing inhibitors of GSK-3β and MAPK/ERK . It has been shown that pluripotent stem cells have compromised mitochondrial function compared to differentiated cells . Therefore, we speculate that there are important differences between the mouse ESCs and differentiated BC cells in intrinsic metabolic preferences and maintenance of the ratio of succinate/fumarate to α-KG. Also, there might be additional effects caused by the inhibition of GSK-3β and MAPK/ERK signaling pathways, as reported by Carey et al. , on the metabolic state. We therefore propose that increasing the ratio of succinate/fumarate to α-KG may result in opposite effects in ESCs and BC cells. Because both succinate and α-KG are crucial for regulating cellular signaling, the effective ratio of succinate to α-KG may be pre-determined by the cellular context and experimental paradigm.
It is well known that tumorigenic cancer cells or CSC-like cells increase tumor malignancy, but how cancer cells retain or acquire tumorigenic properties remains elusive. Metabolic reprogramming, such as the Warburg effect, is emerging as a mechanism for meeting the cancer cell’s need for efficient anabolic metabolism. In the current study, we describe a model whereby cells acquire or maintain tumorigenicity by metabolic reprogramming, without additional mutations or exogenous manipulation of proto-oncogenes or tumor suppressor genes. The importance of tumor tissue niches, such as the hypoxic niche, in tumorigenicity, has been recognized [16, 61]. However, the impact of key nutrients on tumorigenic populations remains largely unknown. The classic CSC model follows a unidirectional, hierarchical structure, whereby CSCs differentiate into non-CSC progeny cells . Recent studies have revealed that cues from the environment perturb the phenotypic equilibrium of cancer cell populations and enable the bidirectional inter-conversion of BC cells [12–14]. Based on our results, we further postulate that altered metabolism, induced by the nutritional microenvironment, could regulate cancer cell plasticity by promoting the conversion of non-CSC to CSC cells.
Moreover, the DKG-induced overall metabolic phenotype is consistent with previous reports on the role of HIF-1α in tumor expansion. It has been shown that CSCs, or tumorigenic cells, rely on glycolysis and using DCA to inhibit glycolysis can eliminate CSC populations and retard tumor growth . Therefore, there is interest in the importance of mitochondrial metabolism in cancer progression and the maintenance of tumorigenicity [7, 8, 63]. Additional studies also highlighted a critical role for reductive carboxylation, mediated by hypoxia or HIF-1α, in promoting the lipogenic phenotype and supporting tumour growth [36, 56, 64–66]. Recently, a report indicated that ovarian cancer spheroid cells with tumorigenic properties demonstrate hypoxia-resistant metabolism, including enhanced glycolysis, pentose phosphate pathway activity and reductive carboxylation . Another study demonstrated that leukemic cells exhibit increased FAO, uncoupled with ATP synthesis, and rely on de novo lipogenesis . It is possible that DKG treatment of BC cells increases HIF-1α-dependent glycolysis and reductive carboxylation, which fuels FAO. Consistent with reports that FAO is a key mechanism of pancreatic tumor relapse  and is implicated in resistance to oxidative stress in glioblastoma cells , we demonstrate that inhibiting both glycolysis and FAO partially reduces DKG-induced tumorsphere formation.
In summary, metabolic reprogramming and the switch to glycolysis in cancer cells have received a substantial amount of interests over the past several decades, but the metabolic state has been poorly linked to the cancer cellular hierarchy. Our results suggest that breast cancer tumorigenicity can increase due to metabolic reprogramming and HIF-1α-dependent signaling. Treatment with the α-KG precursor, DKG, causes a cascade of metabolic changes and the induction of pluripotency markers, contributing to tumor expansion in vivo. Notably, cancer cells have multiple systems to sense energy imbalance and can modify signaling pathways to favor survival and growth. Along these lines, a dramatic increase in the prevalence of obesity and related metabolic disorders over the last few decades is associated with increased mortality from heart disease, stroke, diabetes, and cancer . The nutrient excess associated with the obese state may affect cancer risk and aggressiveness in previously unknown ways, and this is an area of active research. Reminiscent of our report, the key TCA cycle intermediary metabolite, α-KG, is elevated in female ob/ob mice  and morbidly obese patients . Our results are the first, to our knowledge, to provide insight into the potential tumor-promoting effects of excess exogenous α-KG. Therefore, our findings support a more comprehensive approach when devising therapeutic strategies in cancer, which should consider not only the presence of oncogenes and tumor suppressors, but also how the nutritional microenvironment might affect the subpopulation of cancer cells with tumorigenic CSC-like traits.
MATERIALS AND METHODS
Cell lines and reagents
All cell lines were cultured in a humidified 5% CO2 incubator at 37°C. The human BC cell lines MDA-MB-231, MCF7 and MDA-MB-468 were all obtained from American Type Culture Collection (ATCC) and cultured in Dulbecco’s modified Eagle’s medium (DMEM) containing 10% fetal bovine serum (FBS) and penicillin (100 U/ml)-streptomycin (100 μg/ml). T47D cells obtained from ATCC were cultured in RPMI supplemented with 10% FBS and antibiotics. MDA-MB-231 cells transduced with control shRNA (shCON) or shRNA targeting HIF-1α (shHIF-1α) were maintained in full DMEM containing puromycin (1 μg/ml). DKG, DCA, ETO and DMOG were purchased from Sigma-Aldrich.
Human primary BC cells
Human primary BC cells, BC_01 and BC_03, were derived from patients treated at Taipei Medical University according to Institutional Review Board-approved protocols. Both BC_01 and BC_03 are invasive ductal carcinoma, ER-positive, PR-positive and HER2-negative. Isolated BC cells were grown on NIH J2 feeder cells in the presence of the Rho-associated protein kinase (ROCK1) inhibitor Y-27632 (Enzo Life Sciences). The growth medium contained 3:1 F12 nutrient mixture: DMEM with FBS (5%), EGF (10 ng/ml, Invitrogen), cholera toxin (10 ng/ml, Calbiochem), insulin (5 μg/ml, Gibco), Y-27632 (5 μM) and dexamethasone (0.1 μM, Sigma-Aldrich). The combination of the growth media and the usage of a J2 fibroblast feeder layer induces normal and tumor epithelial cells from many tissues to proliferate in vitro without transduction of exogenous viral or cellular genes [73, 74].
Western blots and antibodies
For protein expression studies, whole cell lysates were extracted using SDS lysis buffer and then supplemented with cOmplete protease inhibitor mixture (Roche Applied Science). Equal amounts of whole cell lysates were separated by SDS-PAGE and immunoblotted with antibodies that recognized HIF-1α (BD Bioscience), OCT4 (BioLegend). Anti-GAPDH (Santa Cruz Biotechnology, Inc.) and anti-Actin (Millipore) antibodies were used to assess equal protein loading. Immunoblots were visualized using an enhanced chemiluminescence detection kit (ECL-Plus, Amersham Pharmacia Biotech) and were imaged with a Versadoc 3000 Imaging System (Bio-Rad). Densitometric data were obtained and quantified with Quantity One Software (Bio-Rad). Western analyses shown are representative of two to four independent experiments.
NMR spectrometry-based metabolite measurement
The metabolites in the DKG-treated MDA-MB-231 cells were measured using solution nuclear magnetic resonance (NMR) spectroscopy. Briefly, metabolites from an equal number of cells treated with DKG (10 mM, 0, 2, 8, 24 h, three independent samples for each time point) were extracted as described previously . The lyophilized, water-soluble metabolites were resuspended in 0.5 ml 99.96% D2O containing 50.9 μM 4,4-dimethyl-4-silapentane-1-sulfonic acid (DSS), which served as the 1H chemical-shift reference and an internal concentration standard. Three replicates were prepared. NMR spectra were acquired at 25°C on a Bruker Avance spectrometer equipped with a cryoprobe operating at 600 MHz 1H frequency. 1H NMR spectra were acquired with the Bruker noesygppr1d pulse sequence using a mixing time of 100 ms, recycle delay of 1 s and 2048 scans. Water suppression was achieved by saturation during the recycle delay and during the mixing time. The acquisition time is 4 s with 64 k data points and 13 ppm spectral width. 2D TOCSY spectra were acquired using the Bruker dipsi2phpr pulse sequence on one of each triplicate samples. The number of scans for each free-induction delay (FID) is 64. Complex data points of 4096 and 512 with spectra width of 13.3 and 11 ppm were used for direct and indirect dimensions, respectively. The recycle delay was 1.5 s, and the mixing time was 80 ms. The NMR spectra were processed using the Bruker topspin software or Chenomx, and the metabolite concentrations were calibrated against DSS using Chenomx NMR Suite Processor. The Chenomx NMR Suite Profiler software was used to identify the metabolites. 2D TOCSY spectra were used to confirm the assignments of some metabolites when ambiguity occurred in the 1D spectra. The fold changes were calculated by comparing the concentrations of metabolites (millimolar) in the treated samples (2, 8, 24 h) with the untreated samples (0 h).
RNA sequencing and bioinformatic analysis
Total RNA was isolated from cells using the TRIzol reagent (Ambion) for RNA-sequencing (RNA-seq) according to the manufacturer’s instructions. RNA-seq was conducted by the Integrative Genomics Core at City of Hope. Ingenuity Pathway Analysis was used for pathway analysis and functional annotation among untreated and DKG-treated cells. The genes involved in glycolysis and the electron transport chain were identified and then the heat map was generated using the GENE-E software (Broad Institute). For MDA-MB-231 cells, two biological replicates were subjected to RNA-seq analysis. The results were analyzed separately, as shown in the figures. For primary BC cells, one sample from each cell line was subjected to RNA-seq analysis. The RNA-seq data are available in Gene Expression Omnibus (GEO) database (GSE47031).
Tumorsphere formation assay
The assay was performed according to a previous report . The cells were seeded in triplicate in an ultra-low attachment plate (Corning) at 20,000 cells/ml. The spheres were counted 7 to 10 days after seeding the cells. Spheres with diameters ≥ 200 μm were enumerated and the bar graphs show the quantitative results from three independent wells. For serial passaging, day-3 tumorspheres were collected, dissociated by trypsin, and then seeded at 5,000 cells/ml and 1,000 cells/ml for the secondary and tertiary cultures, respectively.
Immunofluorescent staining, intracellular ROS staining, flow cytometry and cell sorting
To detect cell surface markers, the cells were trypsinized, blocked in ice-cold PBS containing bovine serum albumin (BSA, 2%) and 0.1% NaN3, then incubated at 4°C with an anti-CD133/2 (293C3)-APC (Miltenyi Biotec), anti-CD44-APC or anti-CD24-PE (BD Pharmingen) antibody for 30 mins. ROS detection was done by DCFDA (Sigma) labeling followed by detection of the oxidized DCF signal. Cellular fluorescence was measured using 50,000 cells on an Accuri C6 Flow Cytometer (BD Biosciences). FACS was performed using an AriaSORT cell sorter (BD Biosciences). Data were processed using the FlowJo software package (Tree Star Inc.).
ECAR and OCR
ECAR and OCR were measured as previously described  by the Seahorse Bioscience XFe24 Extracellular Flux Analyzer (Seahorse Bioscience) to track glycolysis and OXPHOS, respectively. The cells (6 × 104) were seeded into each well the night before the measurements. The measurements were performed according to the manufacturer’s instructions. Briefly, the cells were washed twice with XF Base medium (Seahorse Bioscience) plus L-Glutamine (2 mM final concentration), and the pH was adjusted to 7.35. ECAR was assayed by sequential injections of D-glucose (Glc; 10 mM), oligomycin (1 μM) and 2-deoxyglucose (2-DG, inhibitor of hexokinase; 50 mM; Sigma-Aldrich). OCR was also measured independently, as described previously . The mito-stress test was performed by sequential injections of oligomycin (1 μM), FCCP (0.5 μM) and Rotenone (2.5 μM), according to the manufacturer’s instructions. To integrate the ECAR and OCR data together, we chose the sixth point (22.5 mins after Glc injection) of the ECAR measurement for the ECAR vs. OCR plots. After the measurements, the cells were trypsinized for cell number normalization. The data are presented as the mean ± SD, n = 5.
RNA extraction and qRT-PCR
Total RNA was extracted from cells using an RNeasy Mini Kit (Qiagen) according to the manufacturer instructions. The cDNA was synthesized from the total RNA using an iScript cDNA Synthesis Kit (Bio-Rad). Quantitative PCR analyses of targeted sequences were generated using the iTaq SYBR Green Supermix (Bio-Rad), a fraction of each cDNA sample, and gene-specific primers (Table S1). PCR amplification and fluorescence were detected using a MyIQ real-time PCR detection system, and threshold cycles were determined by the iCycler program (default setting). Fold induction was determined using the ΔΔCT method, normalized to GAPDH.
ATP production was measured using the ENLITEN® ATP Assay System (Promega, FF2000) according to the manufacturer’s instructions. The cells were harvested with cold PBS and then 5% trichloroacetic acid (TCA) was added. Tris-acetate-EDTA buffer (pH 7.75) was used to neutralize and dilute the TCA solution to a final concentration of 0.1%. The extracts were further diluted 1:100 and added to an equal volume of rL/L Reagent (Promega, FF2000). Then, the luminescence was measured using a TD-20e luminometer (Turner). The ATP standard (Promega, FF2000) was serially diluted to generate a regression curve for calculation of the exact number of ATP molecules, which was then divided by the cell number to derive the number of ATP molecules per cell.
Lipid droplet staining
To visualize lipid droplets in the cells, Oil Red O staining was performed, as previously described . BODIPY 493/503 (Life Technologies) staining was performed at 2.5 ng/100 μl and analyzed by an Accuri C6 Flow Cytometer (BD Biosciences).
Tumor xenografts and isolation of tumor cells
To evaluate tumorigenesis, 104, 103, or 102 untreated or DKG pre-treated cells were mixed with Matrigel (1:1, Corning) and injected subcutaneously in the left and right flanks of the same NOD.Cg-PrkdcscidIl2rgtm1Wjl/SzJ (NSG) mouse, respectively. Tumor weight was measured 2 months after injection. To further investigate the role of HIF-1α in the xenograft model, 104 or 106 of MDA-MB-231/shCON or shHIF-1α cells pre-treated with DKG were resuspended in 100 μl serum-free DMEM media, and then injected subcutaneously into the right flanks of six-week-old female NSG mice. The tumor volumes were measured twice weekly and the mice were euthanized on day 65, and the tumors were harvested and weighed. The tissue sectioning and the IHC staining were done by the Pathology Core at City of Hope, using an anti-CAIX antibody (Genetex), anti-CD44 antibody (Cell Signaling) and anti-CD133 antibody (Fitzgerald). Isolation and culture of tumor cells from the mouse xenografts were performed as previously described . Briefly, the tumors were minced into small pieces and then digested in collagenase-hyaluronidase-DNase I (Sigma) buffer for 60 minutes at 37 °C. The cell suspension was then passed through a 100 μm cell strainer (BD Biosciences) and the resulting cells were plated on petri dishes and maintained in the growth medium for culturing MDA-MB-231 cells.
Ionizing radiation (IR)
A Shepherd Mark I Cesium-137 γ irradiator was used to irradiate cells at a fixed dose rate of 1 to 2 Gy/min. 1-4 Gy was used in this study.
Clonogenic survival assay
DKG-pretreated cells were seeded in triplicate in six-well plates at 100 cells/well in complete growth medium one day before IR. After 10 days, cells were washed with PBS and fixed in ice-cold methanol for 10 mins, then stained with crystal violet solution (0.5%, 10 mins). Colonies were counted and the surviving fraction was calculated using the plating efficiency.
Data are presented as the mean ± SD. Statistical significance was determined using the two-tailed Student’s t-test. p < 0.05 is considered significant. *: p < 0.05, **: p < 0.01, ***: p < 0.005.
We are sincerely grateful to Drs. Arthur Riggs, Rama Natarajan, Emily Wang and Mei Kong for their helpful suggestions and critical reading of manuscript. We thank Dr. Jun Wu, Ms. Yen Wang and Ms. Meng-Ying Hsieh of the Animal Tumor Core for xenograft efficacy studies, Ms. Sofia Loera of the Pathology Core for IHC analyses, Ms. Lucy Brown of the Analytical Cytometry Core for flow cytometry analyses, Dr. Weidong Hu of the NMR core for metabolomics studies, members of Dr. Ann’s laboratory for helpful discussions and Dr. Nancy Linford for editing.
CONFLICTS OF INTEREST
The authors indicate no potential conflicts of interest.
National Institute of Health Research Grants R01DE10742 and R01DE14183 and The Mary Kay Foundation Research Grant number 005-13 (to D.K.A.) and P30CA33572. MOST 104-2321-B-400-009 MOHW104-TDU-M-221-123017 NIH CA165263 and CA150197 (to H.J.K.).
1. Hanahan D, Weinberg R. Hallmarks of Cancer: The Next Generation. Cell. 2011; 144:646–674.
2. Lu C, Thompson CB. Metabolic Regulation of Epigenetics. Cell Metabolism. 2012; 16:9–17.
3. Ward PS, Thompson CB. Metabolic Reprogramming: A Cancer Hallmark Even Warburg Did Not Anticipate. Cancer Cell. 2012; 21:297–308.
4. Vogelstein B, Kinzler KW. Cancer genes and the pathways they control. Nat Med. 2004; 10:789–799.
5. Warburg O, Wind F, Negelein E. The metabolism of tumors in the body. The Journal of General Physiology. 1927; 8:519–530.
6. Vander Heiden MG, Cantley LC, Thompson CB. Understanding the Warburg Effect: The Metabolic Requirements of Cell Proliferation. Science. 2009; 324:1029–1033.
7. LeBleu VS, O’Connell JT, Gonzalez Herrera KN, Wikman H, Pantel K, Haigis Marcia C, de Carvalho FM, Damascena A, Domingos Chinen LT, Rocha RM, Asara JM, Kalluri R. PGC-1α mediates mitochondrial biogenesis and oxidative phosphorylation in cancer cells to promote metastasis. Nat Cell Biol. 2014; 16:992–1003.
8. Tan An S, Baty James W, Dong L-F, Bezawork-Geleta A, Endaya B, Goodwin J, Bajzikova M, Kovarova J, Peterka M, Yan B, Pesdar Elham A, Sobol M, Filimonenko A et al. Mitochondrial Genome Acquisition Restores Respiratory Function and Tumorigenic Potential of Cancer Cells without Mitochondrial DNA. Cell Metabolism. 2015; 21:81–94.
9. Katajisto P, Döhla J, Chaffer CL, Pentinmikko N, Marjanovic N, Iqbal S, Zoncu R, Chen W, Weinberg RA, Sabatini DM. Asymmetric apportioning of aged mitochondria between daughter cells is required for stemness. Science. 2015; 348:340–343.
10. Xie Q, Wu Q, Horbinski CM, Flavahan WA, Yang K, Zhou W, Dombrowski SM, Huang Z, Fang X, Shi Y, Ferguson AN, Kashatus DF, Bao S et al. Mitochondrial control by DRP1 in brain tumor initiating cells. Nat Neurosci. 2015; 18:501–510.
11. Visvader J, Lindeman G. Cancer Stem Cells: Current Status and Evolving Complexities. Cell Stem Cell. 2012; 10:717–728.
12. Chaffer C, Brueckmann I, Scheel C, Kaestli A, Wiggins P, Rodrigues L, Brooks M, Reinhardt F, Su Y, Polyak K, Arendt L, Kuperwasser C, Bierie B et al. Normal and neoplastic nonstem cells can spontaneously convert to a stem-like state. Proceedings of the National Academy of Sciences. 2011; 108:7950–7955.
13. Chaffer C, Marjanovic N, Lee T, Bell G, Kleer C, Reinhardt F, D’Alessio A, Young R, Weinberg R. Poised Chromatin at the ZEB1 Promoter Enables Breast Cancer Cell Plasticity and Enhances Tumorigenicity. Cell. 2013; 154:61–74.
14. Gupta P, Fillmore C, Jiang G, Shapira S, Tao K, Kuperwasser C, Lander E. Stochastic State Transitions Give Rise to Phenotypic Equilibrium in Populations of Cancer Cells. Cell. 2011; 146:633–644.
15. Semenza G. Hypoxia-Inducible Factors in Physiology and Medicine. Cell. 2012; 148:399–408.
16. Wong CC-L, Gilkes DM, Zhang H, Chen J, Wei H, Chaturvedi P, Fraley SI, Wong C-M, Khoo U-S, Ng IO-L, Wirtz D, Semenza GL. Hypoxia-inducible factor 1 is a master regulator of breast cancer metastatic niche formation. Proceedings of the National Academy of Sciences. 2011; 108:16369–16374.
17. Semenza G. Molecular mechanisms mediating metastasis of hypoxic breast cancer cells. Trends in Molecular Medicine. 2012; 18:534–543.
18. Wong C-L, Zhang H, Gilkes D, Chen J, Wei H, Chaturvedi P, Hubbi M, Semenza G. Inhibitors of hypoxia-inducible factor 1 block breast cancer metastatic niche formation and lung metastasis. Journal of Molecular Medicine. 2012; 90:803–815.
19. Kaya AO, Gunel N, Benekli M, Akyurek N, Buyukberber S, Tatli H, Coskun U, Yildiz R, Yaman E, B O. Hypoxia inducible factor-1 alpha and carbonic anhydrase IX overexpression are associated with poor survival in breast cancer patients. Journal of BUON. 2012; 17:663–668.
20. Pugh CW, Ratcliffe PJ. The von Hippel–Lindau tumor suppressor, hypoxia-inducible factor-1 (HIF-1) degradation, and cancer pathogenesis. Seminars in Cancer Biology. 2003; 13:83–89.
21. Astuti D, Douglas F, Lennard TWJ, Aligianis IA, Woodward ER, Evans DGR, Eng C, Latif F, Maher ER. Germline SDHD mutation in familial phaeochromocytoma. The Lancet. 2001; 357:1181–1182.
22. Castro-Vega LJ, Buffet A, De Cubas AA, Cascón A, Menara M, Khalifa E, Amar L, Azriel S, Bourdeau I, Chabre O, Currás-Freixes M, Franco-Vidal V, Guillaud-Bataille M et al. Germline mutations in FH confer predisposition to malignant pheochromocytomas and paragangliomas. Human Molecular Genetics. 2014; 23:2440–2446.
23. Gimenez-Roqueplo A-P, Favier J, Rustin P, Rieubland C, Crespin M, Nau V, Van Kien PK, Corvol P, Plouin P-F, Jeunemaitre X, Network ftC. Mutations in the SDHB Gene Are Associated with Extra-adrenal and/or Malignant Phaeochromocytomas. Cancer Research. 2003; 63:5615–5621.
24. Hao H-X, Khalimonchuk O, Schraders M, Dephoure N, Bayley J-P, Kunst H, Devilee P, Cremers CWRJ, Schiffman JD, Bentz BG, Gygi SP, Winge DR, Kremer H et al. SDH5, a Gene Required for Flavination of Succinate Dehydrogenase, Is Mutated in Paraganglioma. Science. 2009; 325:1139–1142.
25. Kim S, Kim D, Jung W-H, Koo J. Succinate dehydrogenase expression in breast cancer. SpringerPlus. 2013; 2:299.
26. Niemann S, Muller U. Mutations in SDHC cause autosomal dominant paraganglioma, type 3. Nature Genetics. 2000; 26:268–270.
27. Pantaleo MA, Astolfi A, Indio V, Moore R, Thiessen N, Heinrich MC, Gnocchi C, Santini D, Catena F, Formica S, Martelli PL, Casadio R, Pession A et al. SDHA Loss-of-Function Mutations in KIT–PDGFRA Wild-Type Gastrointestinal Stromal Tumors Identified by Massively Parallel Sequencing. Journal of the National Cancer Institute. 2011; 103:983–987.
28. Pollard PJ, Briere JJ, Alam NA, Barwell J, Barclay E, Wortham NC, Hunt T, Mitchell M, Olpin S, Moat SJ, Hargreaves IP, Heales SJ, Chung YL et al. Accumulation of Krebs cycle intermediates and over-expression of HIF1 alpha in tumours which result from germline FH and SDH mutations. Human Molecular Genetics. 2005; 14:2231–2239.
29. Vanharanta S, Buchta M, McWhinney SR, Virta SK, Peczkowska M, Morrison CD, Lehtonen R, Januszewicz A, Jarvinen H, Juhola M, Mecklin JP, Pukkala E, Herva R et al. Early-onset renal cell carcinoma as a novel extraparaganglial component of SDHB-associated heritable paraganglioma. American Journal of Human Genetics. 2004; 74:153–159.
30. Xiao M, Yang H, Xu W, Ma S, Lin H, Zhu H, Liu L, Liu Y, Yang C, Xu Y, Zhao S, Ye D, Xiong Y et al. Inhibition of α-KG-dependent histone and DNA demethylases by fumarate and succinate that are accumulated in mutations of FH and SDH tumor suppressors. Genes & Development. 2012; 26:1326–1338.
31. Isaacs JS, Jung YJ, Mole DR, Lee S, Torres-Cabala C, Chung Y-L, Merino M, Trepel J, Zbar B, Toro J, Ratcliffe PJ, Linehan WM, Neckers L. HIF overexpression correlates with biallelic loss of fumarate hydratase in renal cancer: Novel role of fumarate in regulation of HIF stability. Cancer Cell. 2005; 8:143–153.
32. Selak MA, Armour SM, MacKenzie ED, Boulahbel H, Watson DG, Mansfield KD, Pan Y, Simon MC, Thompson CB, Gottlieb E. Succinate links TCA cycle dysfunction to oncogenesis by inhibiting HIF-α prolyl hydroxylase. Cancer Cell. 2005; 7:77–85.
33. Zhao S, Lin Y, Xu W, Jiang W, Zha Z, Wang P, Yu W, Li Z, Gong L, Peng Y, Ding J, Lei Q, Guan K-L et al. Glioma-Derived Mutations in IDH1 Dominantly Inhibit IDH1 Catalytic Activity and Induce HIF-1α. Science. 2009; 324:261–265.
34. Carey BW, Finley LWS, Cross JR, Allis CD, Thompson CB. Intracellular [agr]-ketoglutarate maintains the pluripotency of embryonic stem cells. Nature. 2015; 518:413–416.
35. Lin T-C, Chen Y-R, Kensicki E, Li AY-J, Kong M, Li Y, Mohney RP, Shen H-M, Stiles B, Mizushima N, Lin L-I, Ann DK. Autophagy: Resetting glutamine-dependent metabolism and oxygen consumption. Autophagy. 2012; 8:1477–1493.
36. Wise DR, Ward PS, Shay JES, Cross JR, Gruber JJ, Sachdeva UM, Platt JM, DeMatteo RG, Simon MC, Thompson CB. Hypoxia promotes isocitrate dehydrogenase-dependent carboxylation of α-ketoglutarate to citrate to support cell growth and viability. Proceedings of the National Academy of Sciences. 2011; 108:19611–19616.
37. Hou P, Kuo C-Y, Cheng C-T, Liou J-P, Ann DK, Chen Q. Intermediary Metabolite Precursor Dimethyl-2-Ketoglutarate Stabilizes Hypoxia-Inducible Factor-1α by Inhibiting Prolyl-4-Hydroxylase PHD2. PLoS ONE. 2014; 9:e113865.
38. Chaturvedi P, Gilkes DM, Wong CCL, Kshitiz, Luo W, Zhang H, Wei H, Takano N, Schito L, Levchenko A, Semenza GL. Hypoxia-inducible factor–dependent breast cancer–mesenchymal stem cell bidirectional signaling promotes metastasis. The Journal of Clinical Investigation. 2013; 123:189–205.
39. Al-Hajj M, Wicha M, Benito-Hernandez A, Morrison S, Clarke M. Prospective identification of tumorigenic breast cancer cells. Proceedings of the National Academy of Sciences. 2003; 100:3983–3988.
40. Brugnoli F, Grassilli S, Piazzi M, Palomba M, Nika E, Bavelloni A, Capitani S, Bertagnolo V. In triple negative breast tumor cells, PLC-beta2 promotes the conversion of CD133high to CD133low phenotype and reduces the CD133-related invasiveness. Molecular Cancer. 2013; 12:165.
41. Ricardo S, Vieira AF, Gerhard R, Leitão D, Pinto R, Cameselle-Teijeiro JF, Milanezi F, Schmitt F, Paredes J. Breast cancer stem cell markers CD44, CD24 and ALDH1: expression distribution within intrinsic molecular subtype. Journal of Clinical Pathology. 2011; 64:937–946.
42. Zhao P, Lu Y, Jiang X, Li X. Clinicopathological significance and prognostic value of CD133 expression in triple-negative breast carcinoma. Cancer Science. 2011; 102:1107–1111.
43. Patel SA, Ramkissoon SH, Bryan M, Pliner LF, Dontu G, Patel PS, Amiri S, Pine SR, Rameshwar P. Delineation of breast cancer cell hierarchy identifies the subset responsible for dormancy. Sci Rep. 2012; 2.
44. Lu X, Mazur SJ, Lin T, Appella E, Xu Y. The pluripotency factor nanog promotes breast cancer tumorigenesis and metastasis. Oncogene. 2014; 33:2655–2664.
45. Piva M, Domenici G, Iriondo O, Rábano M, Simões BM, Comaills V, Barredo I, López-Ruiz JA, Zabalza I, Kypta R, Vivanco MdM. Sox2 promotes tamoxifen resistance in breast cancer cells. EMBO Mol Med. 2014; 6:66–79.
46. Wolfer A, Wittner BS, Irimia D, Flavin RJ, Lupien M, Gunawardane RN, Meyer CA, Lightcap ES, Tamayo P, Mesirov JP, Liu XS, Shioda T, Toner M et al. MYC regulation of a “poor-prognosis” metastatic cancer cell state. Proceedings of the National Academy of Sciences. 2010; 107:3698–3703.
47. Yu F, Li J, Chen H, Fu J, Ray S, Huang S, Zheng H, Ai W. Kruppel-like factor 4 (KLF4) is required for maintenance of breast cancer stem cells and for cell migration and invasion. Oncogene. 2011; 30:2161–2172.
48. Kim R, Nam J. OCT4 expression enhances features of cancer stem cells in a mouse model of breast cancer. Lab Animal Research. 2011; 27:147–152.
49. Nishi M, Sakai Y, Akutsu H, Nagashima Y, Quinn G, Masui S, Kimura H, Perrem K, Umezawa A, Yamamoto N, Lee SW, Ryo A. Induction of cells with cancer stem cell properties from nontumorigenic human mammary epithelial cells by defined reprogramming factors. Oncogene. 2014; 33:643–652.
50. Brand MD, Pakay JL, Ocloo A, Kokoszka J, Wallace DC, Brookes PS, Cornwall EJ. The basal proton conductance of mitochondria depends on adenine nucleotide translocase content. Biochemical Journal. 2005; 392:353–362.
51. Chin R, Fu X, Pai M, Vergnes L, Hwang H, Deng G, Diep S, Lomenick B, Meli V, Monsalve G, Hu E, Whelan S, Wang J et al. The metabolite alpha-ketoglutarate extends lifespan by inhibiting ATP synthase and TOR. Nature. 2014; 510:397–401.
52. Samudio I, Harmancey R, Fiegl M, Kantarjian H, Konopleva M, Korchin B, Kaluarachchi K, Bornmann W, Duvvuri S, Taegtmeyer H, Andreeff M. Pharmacologic inhibition of fatty acid oxidation sensitizes human leukemia cells to apoptosis induction. The Journal of Clinical Investigation. 2010; 120:142–156.
53. Listenberger L, Brown D. Fluorescent detection of lipid droplets and associated proteins. Current Protocols in Cell Biology. 2007:24.22.21–24.22.11.
54. DeBerardinis RJ, Thompson CB. Cellular Metabolism and Disease: What Do Metabolic Outliers Teach Us? Cell. 2012; 148:1132–1144.
55. Thompson CB. Metabolic Enzymes as Oncogenes or Tumor Suppressors. New England Journal of Medicine. 2009; 360:813–815.
56. Metallo CM, Gameiro PA, Bell EL, Mattaini KR, Yang J, Hiller K, Jewell CM, Johnson ZR, Irvine DJ, Guarente L, Kelleher JK, Vander Heiden MG, Iliopoulos O et al. Reductive glutamine metabolism by IDH1 mediates lipogenesis under hypoxia. Nature. 2012; 481:380–384.
57. King A, Selak MA, Gottlieb E. Succinate dehydrogenase and fumarate hydratase: linking mitochondrial dysfunction and cancer. Oncogene. 0000; 25:4675–4682.
58. Céline Pinheiro, Bárbara Sousa, André Albergaria, Joana Paredes, Rozany Dufloth, Daniella Vieira, Schmitt F, Baltazar F. GLUT1 and CAIX expression profiles in breast cancer correlate with adverse prognostic factors and MCT1 overexpression. Histology and Histopathology. 2011; 26:1279–1286.
59. Loenarz C, Schofield CJ. Expanding chemical biology of 2-oxoglutarate oxygenases. Nat Chem Biol. 2008; 4:152–156.
60. Xu X, Duan S, Yi F, Ocampo A, Liu G-H, Izpisua Belmonte Juan C. Mitochondrial Regulation in Pluripotent Stem Cells. Cell Metabolism. 2013; 18:325–332.
61. Louie E, Nik S, Chen J-s, Schmidt M, Song B, Pacson C, Chen XF, Park S, Ju J, Chen E. Identification of a stem-like cell population by exposing metastatic breast cancer cell lines to repetitive cycles of hypoxia and reoxygenation. Breast Cancer Research. 2010; 12:R94.
62. Feng W, Gentles A, Nair RV, Huang M, Lin Y, Lee CY, Cai S, Scheeren FA, Kuo AH, Diehn M. Targeting Unique Metabolic Properties of Breast Tumor Initiating Cells. Stem cells. 2014; 32:1734–1745.
63. Janiszewska M, Suvà ML, Riggi N, Houtkooper RH, Auwerx J, Clément-Schatlo V, Radovanovic I, Rheinbay E, Provero P, Stamenkovic I. Imp2 controls oxidative phosphorylation and is crucial for preserving glioblastoma cancer stem cells. Genes & Development. 2012; 26:1926–1944.
64. Gameiro P, Yang J, Metelo A, Pérez-Carro R, Baker R, Wang Z, Arreola A, Rathmell WK, Olumi A, López-Larrubia P, Stephanopoulos G, Iliopoulos O. In Vivo HIF-Mediated Reductive Carboxylation Is Regulated by Citrate Levels and Sensitizes VHL-Deficient Cells to Glutamine Deprivation. Cell Metabolism. 2013; 17:372–385.
65. Mullen AR, Hu Z, Shi X, Jiang L, Boroughs Lindsey K, Kovacs Z, Boriack R, Rakheja D, Sullivan Lucas B, Linehan WM, Chandel Navdeep S, DeBerardinis Ralph J. Oxidation of Alpha-Ketoglutarate Is Required for Reductive Carboxylation in Cancer Cells with Mitochondrial Defects. Cell Reports. 2014; 7:1679–1690.
66. Mullen AR, Wheaton WW, Jin ES, Chen P-H, Sullivan LB, Cheng T, Yang Y, Linehan WM, Chandel NS, DeBerardinis RJ. Reductive carboxylation supports growth in tumour cells with defective mitochondria. Nature. 2012; 481:385–388.
67. Liao J, Qian F, Tchabo N, Mhawech-Fauceglia P, Beck A, Qian Z, Wang X, Huss WJ, Lele SB, Morrison CD, Odunsi K. Ovarian Cancer Spheroid Cells with Stem Cell-Like Properties Contribute to Tumor Generation, Metastasis and Chemotherapy Resistance through Hypoxia-Resistant Metabolism. PLoS ONE. 2014; 9:e84941.
68. Viale A, Pettazzoni P, Lyssiotis CA, Ying H, Sanchez N, Marchesini M, Carugo A, Green T, Seth S, Giuliani V, Kost-Alimova M, Muller F, Colla S et al. Oncogene ablation-resistant pancreatic cancer cells depend on mitochondrial function. Nature. 2014; 514:628–632.
69. Pike LS, Smift AL, Croteau NJ, Ferrick DA, Wu M. Inhibition of fatty acid oxidation by etomoxir impairs NADPH production and increases reactive oxygen species resulting in ATP depletion and cell death in human glioblastoma cells. Biochimica et Biophysica Acta. 2011; 1807:726–734.
70. Prospective Studies C. Body-mass index and cause-specific mortality in 900 000 adults: collaborative analyses of 57 prospective studies. The Lancet. 2009; 373:1083–1096.
71. Won E-Y, Yoon M-K, Kim S-W, Jung Y, Bae H-W, Lee D, Park SG, Lee C-H, Hwang G-S, Chi S-W. Gender-Specific Metabolomic Profiling of Obesity in Leptin-Deficient ob/ob Mice by 1H NMR Spectroscopy. PLoS ONE. 2013; 8:e75998.
72. Rodriguez-Gallego E, Guirro M, Riera-Borrull M, Hernandez-Aguilera A, Marine-Casado R, Fernandez-Arroyo S, Beltran-Debon R, Sabench F, Hernandez M, del Castillo D, Menendez JA, Camps J, Ras R et al. Mapping of the circulating metabolome reveals [alpha]-ketoglutarate as a predictor of morbid obesity-associated non-alcoholic fatty liver disease. Int J Obes. 2015; 39:279–287.
73. Gurkar AU, Chu K, Raj L, Bouley R, Lee S-H, Kim Y-B, Dunn SE, Mandinova A, Lee SW. Identification of ROCK1 kinase as a critical regulator of Beclin1-mediated autophagy during metabolic stress. Nat Commun. 2013; 4:2189–2201.
74. Liu X, Ory V, Chapman S, Yuan H, Albanese C, Kallakury B, Timofeeva OA, Nealon C, Dakic A, Simic V, Haddad BR, Rhim JS, Dritschilo A et al. ROCK Inhibitor and Feeder Cells Induce the Conditional Reprogramming of Epithelial Cells. The American Journal of Pathology. 2012; 180:599–607.
75. Cano KE, Li Y-J, Chen Y. NMR Metabolomic Profiling Reveals New Roles of SUMOylation in DNA Damage Response. Journal of Proteome Research. 2010; 9:5382–5388.
76. Dong C, Yuan T, Wu Y, Wang Y, Fan Teresa WM, Miriyala S, Lin Y, Yao J, Shi J, Kang T, Lorkiewicz P, St Clair D, Hung M-C et al. Loss of FBP1 by Snail-Mediated Repression Provides Metabolic Advantages in Basal-like Breast Cancer. Cancer Cell. 2013; 23:316–331.
77. Qiu F, Chen Y-R, Liu X, Chu C-Y, Shen L-J, Xu J, Gaur S, Forman HJ, Zhang H, Zheng S, Yen Y, Huang J, Kung H-J et al. Arginine Starvation Impairs Mitochondrial Respiratory Function in ASS1-Deficient Breast Cancer Cells. Sci Signal. 2014; 7:ra31-.
78. Fu Y, Cruz-Monserrate Z, Helen Lin H, Chung Y, Ji B, Lin S-m, Vonderfecht S, Logsdon CD, Li C-F, Ann DK. Ductal activation of oncogenic KRAS alone induces sarcomatoid phenotype. Scientific Reports. 2015; 5:13347.
All site content, except where otherwise noted, is licensed under a Creative Commons Attribution 3.0 License.