Lactic Acidosis Triggers Starvation Response with Paradoxical Induction of TXNIP through MondoA
Although lactic acidosis is a prominent feature of solid tumors, we still have limited understanding of the mechanisms by which lactic acidosis influences metabolic phenotypes of cancer cells. We compared global transcriptional responses of breast cancer cells in response to three distinct tumor microenvironmental stresses: lactic acidosis, glucose deprivation, and hypoxia. We found that lactic acidosis and glucose deprivation trigger highly similar transcriptional responses, each inducing features of starvation response. In contrast to their comparable effects on gene expression, lactic acidosis and glucose deprivation have opposing effects on glucose uptake. This divergence of metabolic responses in the context of highly similar transcriptional responses allows the identification of a small subset of genes that are regulated in opposite directions by these two conditions. Among these selected genes, TXNIP and its paralogue ARRDC4 are both induced under lactic acidosis and repressed with glucose deprivation. This induction of TXNIP under lactic acidosis is caused by the activation of the glucose-sensing helix-loop-helix transcriptional complex MondoA:Mlx, which is usually triggered upon glucose exposure. Therefore, the upregulation of TXNIP significantly contributes to inhibition of tumor glycolytic phenotypes under lactic acidosis. Expression levels of TXNIP and ARRDC4 in human cancers are also highly correlated with predicted lactic acidosis pathway activities and associated with favorable clinical outcomes. Lactic acidosis triggers features of starvation response while activating the glucose-sensing MondoA-TXNIP pathways and contributing to the “anti-Warburg” metabolic effects and anti-tumor properties of cancer cells. These results stem from integrative analysis of transcriptome and metabolic response data under various tumor microenvironmental stresses and open new paths to explore how these stresses influence phenotypic and metabolic adaptations in human cancers.
Published in the journal:
. PLoS Genet 6(9): e32767. doi:10.1371/journal.pgen.1001093
Category:
Research Article
doi:
https://doi.org/10.1371/journal.pgen.1001093
Summary
Although lactic acidosis is a prominent feature of solid tumors, we still have limited understanding of the mechanisms by which lactic acidosis influences metabolic phenotypes of cancer cells. We compared global transcriptional responses of breast cancer cells in response to three distinct tumor microenvironmental stresses: lactic acidosis, glucose deprivation, and hypoxia. We found that lactic acidosis and glucose deprivation trigger highly similar transcriptional responses, each inducing features of starvation response. In contrast to their comparable effects on gene expression, lactic acidosis and glucose deprivation have opposing effects on glucose uptake. This divergence of metabolic responses in the context of highly similar transcriptional responses allows the identification of a small subset of genes that are regulated in opposite directions by these two conditions. Among these selected genes, TXNIP and its paralogue ARRDC4 are both induced under lactic acidosis and repressed with glucose deprivation. This induction of TXNIP under lactic acidosis is caused by the activation of the glucose-sensing helix-loop-helix transcriptional complex MondoA:Mlx, which is usually triggered upon glucose exposure. Therefore, the upregulation of TXNIP significantly contributes to inhibition of tumor glycolytic phenotypes under lactic acidosis. Expression levels of TXNIP and ARRDC4 in human cancers are also highly correlated with predicted lactic acidosis pathway activities and associated with favorable clinical outcomes. Lactic acidosis triggers features of starvation response while activating the glucose-sensing MondoA-TXNIP pathways and contributing to the “anti-Warburg” metabolic effects and anti-tumor properties of cancer cells. These results stem from integrative analysis of transcriptome and metabolic response data under various tumor microenvironmental stresses and open new paths to explore how these stresses influence phenotypic and metabolic adaptations in human cancers.
Introduction
Human cancers are extremely heterogeneous due to multiple mutations in oncogenes and tumor suppressor genes, a range of inherited germline variations and varying degrees of microenvironmental stresses. These tumor microenvironmental stresses include tumor hypoxia, accumulation of lactic acid (lactic acidosis) and depletion of glucose, glutamine and other nutrients [1]. These stresses are largely caused by a combination of poor tissue perfusion, abnormal tumor vasculature, uncontrolled proliferation and dysregulated energy metabolism of cancer cells during tumor development and progression. Importantly, these microenvironmental stresses also directly modulate physiological and metabolic phenotypes of cancer cells and ultimately affect the clinical outcomes of patients. With major variations known to exist among different tumors, advances in the pretreatment assessment of the influences of these stresses will aid in improved selection of appropriate therapeutic strategies for individual patients. These stresses and their downstream effects are also the targets of cancer therapeutics, including anti-angiogenesis and hyperthermia treatments. It is therefore important to fully understand the impact and mechanism of how these stresses affect various tumor and non-tumor cells in human cancers.
It is well known that cells resort to glycolysis instead of oxidative phosphorylation to utilize glucose as energy source during hypoxia. In addition, cancer cells have a preferential use of glycolysis pathways for energy generation even in the presence of oxygen – so called “aerobic glycolysis” as first proposed by Dr. Otto Warburg [2]. These factors all likely contribute to high glucose flux and form the basis of using glucose analog 18F-FDG to detect tumor cells. Such dysregulated metabolisms in cancer cells also lead to the accumulation of the metabolic product of glycolysis – lactic acids in solid tumors. Many measurements have been performed to determine the level of tumor lactate and significant variations were found, with the medium range of 7–10 mM/g and up to 25.9 mM/g [3]–[5]. These studies show that high tumor lactate levels are typically associated with more aggressive tumors and resistance to treatment [3]–[5]. How lactic acidosis affects tumor and non-tumor cells in human cancers has been the focus of many elegant studies as summarized in several reviews [6]–[8]. The exposure of cultured cells to lactic acidosis in vitro has been shown to trigger calcium signaling [9], induce angiogensis gene expression (e.g. VEGF, IL8) [10]–[12], HIF-1α stabilization [13] and cell death [14]. Recent genomic analyses identify the transcriptional responses of different cell types to acidosis and high lactate [15]–[17]. These in vitro studies have clearly shown the significant impact of lactic acidosis on the gene expression and phenotypes of cancer cells.
However, it is often challenging to relate the influences of these microenvironmental stresses in vitro to the complex cancer phenotypes in vivo. We have previously overcome this challenge by defining “gene signatures” based on sets of genes whose expression levels are altered by lactic acidosis in vitro as quantitative “common phenotypes”, and projecting them to the in vivo microarray expression data of patients' tumors [18]. We have used this approach to successfully investigate the pathways of hypoxia [19], vascular injury [20] and lactic acidosis [18] in human cancers. This reciprocal flow of information between the in vitro and in vivo systems demonstrate that lactic acidosis triggers significant metabolic reprogramming by forcing cells to rely more on oxidative phosphorylation as an energy source with the suppression of glycolytic phenotypes [18]. Therefore, cellular metabolism is critical in determining the impact of lactic acidosis to tumor phenotypes and clinical outcomes.
In the use of microarrays to connect cultured cells with human cancers, primary epithelial cells are often used to provide the common and shared responses to the defined perturbations due to their intact genetic materials and signaling circuitry [18]–[22]. However, given our intention to relate gene signatures to human cancers, it is probably more relevant to assess signatures in cancer cells. Cancer cell lines are also often easier to be transfected, thus allowing genetic manipulations for mechanistic studies. In this study, we used breast cancer cell line, MCF-7, to conduct a detailed temporal analysis of lactic acidosis response and compare the transcriptomic and metabolic responses of cancer cells to hypoxia, lactic acidosis and glucose deprivation in parallel to gain a further understanding of the underlying molecular mechanisms.
One important goal in understanding the cellular responses to tumor microenvironmental stresses is to identify the key regulator(s) responsible for these observed gene expression response since such understanding will lead to important insights into the development and progression of human cancers. For example, the identification of hypoxia-inducible transcription factors (HIFs) as central regulators in the hypoxia response and its regulation at the level of protein stability has become crucial in our understanding of tumor hypoxia [23]–[26]. HIF-1α protein stabilization can also be seen in the development of multiple neoplasms in patients with von Hippel-Lindau disease [27] or mutations in several enzymes of the tricarboxylic acid (TCA) cycle, such as succinate dehydrogenase (SDH) and fumarate hydratase (FH) [28]–[30]. The glucose deprivation is another feature of tumor microenvironmental stress caused by the imbalance between supply and consumption [31]. Glucose deprivation triggers “starvation-like” signaling through the activation of AMPK and LKB-1, which in turn activates TSC1/TSC2 and inhibits central energy sensor mTOR activities to inhibit ribosomal biogenesis, translation activities and proliferation [32]. The mutations in LKB-1 and TSCs in the glucose sensing pathway lead to cancer development in the Peutz-Jeghers Cancer Syndrome [33] and Tuberous Sclerosis [34], [35], respectively. In contrast, very little is known about the transcriptional regulation of the lactic acidosis gene response program. To uncover the molecular mechanisms by which tumor cells respond to their microenvironment, we conducted a detailed temporal transcriptional analysis of the MCF7 breast cancer cell line under lactic acidosis and compared this response to those elicited by glucose deprivation and hypoxia. We identified a novel growth suppressive pathway (MondoA-TXNIP) that portends better prognosis in breast cancer and contribute to the anti-Warburg effects of lactic acidosis.
Results
The temporal transcriptional responses of MCF-7 under lactic acidosis
We examined the temporal gene expression patterns of a breast cancer cell line MCF-7 exposed to lactic acidosis (25 mM lactic acid with pH 6.7) at various time points during the first 24 hours. Cells were first brought to replicative arrest by serum withdrawal for 24 hours before exposure to either control or lactic acidosis conditions composed of pH 6.7 and 25mM lactic acid in the presence of high levels of glucose (4.5 g/L) in triplicate of different time points at 1, 4, 12 and 24 hrs to characterize the temporal changes of gene expression patterns. The RNA samples harvested from these MCF-7 cells were interrogated with Affymetrix GeneChip U133 plus 2.0 arrays (∼54,000 probe sets on ∼47,000 transcripts and variants) with results deposited in Gene Expression Omnibus (GSE19123).
Gene expression profiles of all MCF-7 cells were normalized by RMA, zero-transformed against the average expression levels of the same probe sets of the time-matched control samples as performed previously [19], [20]. 1761 probes sets showing with at least two fold changes in at least two samples were selected and arranged by hierarchical clustering according to similarities in expression patterns (Figure 1A). This analysis showed that lactic acidosis induced a dramatic change in the gene expression with significant temporal patterns (Figure 1A and 1B).
Among the genes induced in MCF-7 by lactic acidosis, we found PLAU, major histocompatility complex (MHC) type I and CD44, and REDD1 (Figure 1B), a p53 transcriptional target following DNA damage [36]. The expression of CD44 has been reported to be induced by lactosis in cancer cells [37]. Tumor necrosis factor (TNF), Fas and connective tissue growth factors (CTGF) were only induced in earlier time points and returned to baseline in the later time points (Figure 1B). Among clusters of genes repressed by lactic acidosis in the later time points is a large group whose expression is closely linked to cell proliferation (cyclin D, PCNA, CCRK, E2F3, and E2F6) (Figure 1C). These lactic acidosis-repressed genes may reflect a physiological alteration that halts cell proliferation as the cells try to preserve energy consumption under metabolic stress [31], [32].
To define the lactic acidosis response in MCF-7 at a pathway level, we performed Gene Set Enrichment Analysis (GSEA) [38] to compare the pathway composition in the gene expression of all control vs. lactic acidosis samples of MCF-7 cells. We found that samples exposed to lactic acidosis were enriched in gene sets representing nutrient deprivation [39], the treatment of histone deacetylase inhibitor trichostatin A (TSA) in cancer cells, breast cancer good prognosis [40] and exposure of DNA damaging agent (bleomycin) (Table S1). The gene expression of MCF-7 which has been exposed to lactic acidosis was depleted in gene sets representing E2F1 target genes, DNA replication, breast cancer poor prognosis [40], mitotic cycles and RNA processing (Table S1).
We previously showed that the hypoxia and lactic acidosis response signatures elicited in cultured primary non-cancerous cells HMECs could provide a molecular gauge of hypoxia and lactic acidosis response for cancerous human tissues in vivo, and also predict clinical outcomes. To test robustness of this gene signature approach based now on cancerous cells, the lactic acidosis gene signatures generated in MCF-7 were projected to the Miller breast cancer data to evaluate a tumor-specific numerical score representing the predicted signature of lactic acidosis – i.e., quantifying lactic acidosis pathway activity across tumors [41]; see Materials and Methods section below and the detailed statistical supplemental materials from a previous study [18]. Statistical survival analysis then indicates that patients with tumors showing higher levels of lactic acidosis pathway activity – defined by projected signatures of 12 and 24 hours exposure in MCF-7 – had significantly better clinical outcomes (Figure 1D), consistent with our previous studies using HMEC [18]. The lactic acidosis gene signature at 12 hour time point had the most consistent prognostic significance across different tumor datasets, including three other breast cancer expression studies with different stages of diseases (Figure S1). These datasets include a study of 286 lymph node negative early breast cancers from NKI (Wang) [42], and two studies of invasive breast carcinomas (Sotiriou, Pawitan) [43], [44].
We also measured the reproducibility and consistency of the predicted lactic acidosis pathway activities in each tumor using gene signatures generated in both MCF-7 and HMECs. We found that the two variants of the lactic acidosis signature resulted in highly similar predicted lactic acidosis activity, evidenced by strongly correlated signature scores in Miller and other tumor datasets (Figure 1E and Figure S2). This result suggests that although there are differences between the lactic acidosis responses at the gene level, the pathway level predictions based on the gene signatures are similar and reproducible in both normal primary and breast cancer cells.
Lactic acidosis triggers “starvation” response
In the GSEA analysis for the MCF-7 cells exposed to lactic acidosis, we found enrichment in the starvation pathways caused by the deprivation of glucose and glutamine obtained in an independent study [39]. Such association between lactic acidosis and nutrient deprivation was consistent with the reduced ATP production under lactic acidosis [18]. To test for such association, we directly compared the gene expression of MCF-7 cells of lactic acidosis with cellular starvation stress caused by glucose deprivation. We performed parallel global transcriptional analysis on the MCF7 cells in quintuplicate which had been exposed to lactic acidosis and 1% oxygen (hypoxia) in the media with high glucose (4.5g/L) and glucose deprivation (using media 0g/L glucose) for 4 hours. Similar parallel analysis was performed for subsets of genes involved in metastasis and invasive capacity [45], [46]. RNAs from these MCF-7 cells were extracted and hybridized to Affymetrix GeneChip U133 plus 2.0 arrays, normalized by RMA and zero-transformed by deducting the mean expression values for samples under control condition. 3903 probe sets were selected by the criteria of at least two observations with at least two fold changes and arranged by hierarchically clustering (Figure 2A and 2B). Unexpectedly, we noted a highly similar transcriptional response to lactic acidosis and glucose deprivation which is distinct from the hypoxia response (Figure 2A). The similarity of the lactic acidosis and glucose deprivation response was observed in both induction and repression of a large number of common sets of genes, including major histocompatility complex (MHC) type I & II, DNA repair gene, alkB (alkyation repair homolog 7), and CTGF, Jun and TNF (Figure 2B), which were also seen in our previous time course experiment (Figure 1B). In contrast, hypoxia elicited a very distinct transcriptional response with the upregulation of many known hypoxia-inducible genes such as EGLN3, CA9, stanniocalcin1 (STC1), BNIP3 and genes involved in glycolysis, including pyruvate kinase (PKM2) (Figure 2B) [19]. In a previous microarray study of the stress response in yeast Saccharomyces cerevisiae, the induction of a share set of “common stress genes” was a prominent feature [47]. In contrast, there was only a small cluster of genes which were induced by lactic acidosis, glucose deprivation and hypoxia, which may represent “common stress genes” in MCF-7 cells (Figure 2B). These genes include HIG-2 (hypoxia-inducible gene) and REDD, both genes reported to be induced by hypoxia (Figure 2B).
Since glucose deprivation is known to trigger cellular starvation response, this similarity in the transcriptional responses suggests that lactic acidosis may also trigger a “starvation” response seen for nutrient deprivation. To test this possibility, we measured the effects of lactic acidosis on several biochemical markers of starvation response in MCF-7 (Figure 2C). AMP-activated protein kinase (AMPK) is a highly conserved energy-sensing heterotrimeric complex that plays a key role in the regulation of energy homeostasis; it becomes phosphorylated at Thr172 of AMPK α by an elevated AMP/ATP ratio due to various stress indicating energy stress and starvation. We find that lactic acidosis significantly increased Thr172 phosphorylation of AMPK, even in the presence of high levels of glucose and amino acids in the media (Figure 2C). The Mammalian Target of Rapamcyin (mTOR) is another crucial cellular sensor for energy status and a crucial downstream target of AMPK [48], [49]. We examined how lactic acidosis affected mTORC1 activities through the phosphorylation of its downstream target S6 kinase (S6K) and found it dramatically reduced the S6K phophorylation (Figure 2C). Taken together, these data supported the notion that lactic acidosis triggers a cellular starvation response similar to that of glucose starvation, with AMPK activation and mTOR inhibition. These biochemical changes are likely to contribute to reduced cell growth and proliferation under lactic acidosis as evidenced by gene expression (Figure 1C).
The distinct metabolic effects of lactic acidosis and glucose deprivation
Knowing that lactic acidosis can inhibit tumor glycolysis [18], we evaluated how lactic acidosis and glucose deprivation affect glucose uptake using 2-deoxy-D-[2,6-3H]-glucose. While glucose deprivation increased the uptake of the glucose by 84%, lactic acidosis significantly reduced the glucose uptake by 67% (Figure 3A). The inhibition of glucose uptake by lactic acidosis was consistent with the reduced glucose consumption and decreased lactate production we have reported previously. However, this result was not consistent with the expected increased uptake of glucose to increase energy generation associated with AMPK activation and cellular starvation, as seen in glucose deprivation [50]. Therefore, in spite of the shared AMPK activation, mTORC1 inhibition and similar transcriptional responses, lactic acidosis and glucose starvation triggered opposite effects on glucose uptake.
To identify transcriptionally regulated genes that may contribute to these opposite metabolic responses, we used the discordance to select genes with significant but opposite expression changes under lactic acidosis and glucose deprivation. This analysis identified 115 probe sets (Table S2); 49 were induced by lactic acidosis and repressed by glucose deprivation while 66 were repressed by lactic acidosis and induced by glucose deprivation (Figure 3B). The three probe sets with the largest induction levels under lactic acidosis were all associated with thioredoxin interacting protein (TXNIP or Vitamin D3-upregulated protein 1, VDUP1) (Figure 3B). TXNIP was strongly induced by lactic acidosis, and suppressed by glucose deprivation. While hypoxia was reported to trigger TXNIP expression [51], we only noted a modest increase (∼40%) in our experiment. We tested for the prognostic significance of this differentially expressed signature by projecting the 115 gene signature into the Miller breast tumor data set. Due to different arrays used in that data set, this results in 109 of the 115 signature probes present in the Miller data set. We find that higher levels of this signature also indicate favorable outcomes (Figure S3A). The three TXNIP probe sets alone in fact have prognostic significance comparable with that of the signature (Figure S3B); on re-evaluating the signature after removal of these three TXNIP probe sets, we find that the reduced signature maintains prognostic significance albeit at a slightly reduced level (Figure S3C).
When we examined expression of MCF-7 and HMEC under lactic acidosis at different time points, we also noted significant induction of TXNIP and its parralogue alpha-arrestin domain containing 4 (ARRDC4) (Figure S4). Real-time RT-PCR further confirmed the induction of TXNIP (up to 4–5 fold) and ARRDC4 transcripts (up to 3.4 fold) in MCF-7 (Figure 3C and Figure S5). In addition, lactic acidosis also induced TXNIP expression in other cancer cell lines, including WiDr (colon cancer cell) and SiHa (cervical cancer cell) (Figure S6).
TXNIP expression is known to be induced by glucose exposure [52], [53]. Since there were high levels of glucose (4.5g/dL) in both control and lactic acidosis media, we further clarified the role of glucose vs. lactic acidosis in the TXNIP induction using real-time PCR to determine the level of TXNIP in MCF-7 cultured under different glucose and lactic acidosis conditions (Figure 3D). Glucose (4.5g/dL) or lactic acidosis alone increased TXNIP levels by approximately 4–5 fold (Figure 3D). In the presence of both high glucose levels and lactic acidosis, TXNIP was induced up to 18 fold (Figure 3D). Thus lactic acidosis and glucose exposure are both potent inducers of TXNIP with synergistic induction potential, suggesting distinct mechanisms of TXNIP induction by these two stimuli.
To define the individual contribution of lactosis vs. acidosis, we tested the level of TXNIP induction under different degrees of lactosis (12.5 or 25 mM lactate) and acidosis (pH 6.7) (Figure 3E and 3F). We found that acidosis alone, but not lactosis alone, led to the TXNIP induction, consistent with our previous study [18] (Figure 3E). Although lactosis alone did not induce TXNIP, the addition of lactate to the acidosis conditions led to a dose-dependent augmentation of the TXNIP induction (Figure 3E). Importantly, this effect on TXNIP induction was present even under 10 or 12.5 mM lactate, a level seen in many human tumors (Figure 3E and 3F). We also measured how lactosis and acidosis affected glucose uptake, and found that acidosis, but not lactosis, led to a reduction in glucose uptake, a pattern consistent with TXNIP induction (Figure 3G).
Although TXNIP was initially identified as a protein interacting with thioredoxin and modulating cellular responses to oxidative stresses, it has the ability to inhibit glucose uptake and is being recognized as an important regulator of dysregulated metabolism in diabetes [54], [55]. Therefore, the upregulation of TXNIP under lactic acidosis make it an attractive candidate contributing to the anti-Warburg effects and inhibition of tumor glycolysis under lactic acidosis [18]. To investigate this, we silenced the TXNIP transcripts with two independent siRNAs and confirmed the successful reduction of the TXNIP protein (Figure 4A). Consistent with the previously known role of TXNIP to inhibit glucose uptake [55], the silencing of TXNIP by siRNAs led to significant decrease in the reduction of glucose uptake under lactic acidosis culture condition (Figure 4B). Lactic acidosis caused 52% repression of glucose uptake in MCF7 cells transfected with non-targeting siRNAs (-). In cells transfected with two different siTXNIPs (T1, T2), glucose uptake was increased while the repressing effect of lactic acidosis was decreased to 39% and 44% respectively (Figure 4B). In addition, we found the silencing of TXNIP in MCF-7 increased both the glucose consumption and lactate productions under normal media, but lactic acidosis significantly reduced both the glucose consumption and lactate production in all treated cells (Figure S7).
Although the TXNIP gene silencing reduced the glycolysis inhibition under lactic acidosis, the effects were modest, which may be due to the remaining level of TXNIP. We further tested the effect of TXNIP knocking out in the TXNIP deficient MEF cells [56]. Lactic acidosis caused 68% reduction in glucose uptake in the wild-type (WT) MEF cells. In contrasts, lactic acidosis only reduced the glucose uptake by 28% in TXNIP knockout (TKO) MEF cells (Figure 4C). To further test for the role of TXNIP upregulation in the lactic acidosis response, we exposed the wild type and TXNIP deficient MEF cells to the lactic acidosis conditions (10 mM lactate, pH 6.7) and performed gene expression analysis using the Affymetrix mouse 430A2 GeneChip and normalized the expression data by RMA. We first examined the effects of the disruption of TXNIP on the gene expression of MEF cells by performing zero transformation of TXNIP deficient against the average expression in the wild type MEF cells (Figure S8). We found 798 probe sets whose expression was altered at least 1.7 fold in at least two samples (Figure S8). Among the genes repressed in the TXNIP KO MEF cells were TXNIP, COX2, collagens and many other genes in the HOX genes involved in pattern specification (Figure S8). Among the induced in the TXNIP KO MEF cells were genes in the complement activation, AP2, and Notch3 (Figure S8).
To define the role of TXNIP in the lactic acidosis gene expression, we compared the respective lactic acidosis response of the wild type and TXNIP deficient MEF cells by zero-transformation against the corresponding cells cultured in the control conditions. 1327 probe sets showing with at least 1.7 fold changes in at least two samples were selected and arranged by hierarchical clustering according to similarities in expression patterns (Figure 4D). This analysis showed that the loss of TXNIP reduced the induction and repression of a significant number of genes as shown in the three gene clusters (Figure 4D), including the repression of many cell-cycle related genes under lactic acidosis. In addition, the loss of TXNIP led to a statistically significant reduction in the overall degrees of induction and repression (p<0.0001) for the average expression of the 1049 repressed and 278 induced genes (Figure S9A and Figure 9B). We also tested the ability of lactic acidosis to repress the proliferation of MEF cells and found that the proliferation of the wild type, but not TXNIP deficient MEF cells, were repressed by lactic acidosis (Figure 4E). This observation is consistent with the reduced repression of cell cycle genes under lactic acidosis seen for the TXNIP deficient MEF cells (Figure 4D). Taken together, these data showed that TXNIP induction contributes to the inhibition of glycolysis phenotypes, metabolic reprogramming, and cell cycle arrest and gene expression under lactic acidosis.
The role of the MondoA:Mlx complex in the activation of TXNIP by lactic acidosis
Glucose-induced TXNIP transcription depends on a short proximal region of the TXNIP promoter. Specifically, this region includes a well conserved Carbohydrate Response Elements (ChoRE) consisting of two E-boxes [53]. We first tested the influence of lactic acidosis on the reporter constructs driven by the promoters of TXNIP and ARRDC4 and found that lactic acidosis could induce the reporter activities of both constructs by more than 6 folds (Figure 5A). Importantly this induction of reporter activity under lactic acidosis was reduced by 61% with mutation in the ChoRE of the TXNIP promoter (Figure 5A), indicating the importance of the ChoRE regions to the transcriptional activation of TXNIP under lactic acidosis.
The transcriptional activation through two E-boxes in the promoters of TXNIP upon glucose exposure has been reported to be caused by the binding MondoA:Mlx [57] or Carbohydrate Response Elements-binding protein (ChREBP) [58]. MondoA is likely to be more relevant in the observed lactic acidosis response of MCF-7 given its higher expression levels in MCF-7 and the simultaneous upregulation of both TXNIP and ARRDC4 under lactic acidosis [57]. MondoA:Mlx complexes are held latently at the outer mitochondrial membrane (OMM), yet shuttle between the OMM and the nucleus suggesting that they facilitate communication between these two essential organelles. MondoA:Mlx complexes are sensors of intracellular glucose concentration and accumulate in the nucleus following increases in glucose-6 phosphate to occupy the E-box-containing promoters of targets such as TXNIP and ARRDC4 [57]. To directly test the inducible physical binding of MondoA to the promoter regions of TXNIP and ARRDC4 under lactic acidosis, chromatin-immunoprecipitation (CHIP) was performed with the antibody against MondoA in MCF-7 cells which have been exposed to lactic acidosis conditions of different acidity. We detected increased specific binding of MondoA to the promoter regions of TXNIP and ARRDC4 with more acidic lactic acidosis environments (Figure 5B). These results indicated that the MonoA became activated to occupy the promoters of TXNIP and ARRDC4 under lactic acidosis and these binding sites were important for their transcriptional inductions.
To further determine the role of MondoA in the induction of TXNIP under lactic acidosis, we used two different sets of siRNAs to knock down MondoA by gene silencing (Figure 5C). During lactic acidosis, the level of MondoA protein did not change, while TXNIP was induced significantly (Figure 5C). This induction of TXNIP was significantly reduced at both levels when MondoA was reduced by both sets of siRNAs targeting MondoA (Figure 5C). The silencing of MondoA also led to increased glucose uptake under both control and lactic acidosis condition (Figure 5D), a result similar to the silencing of TXNIP (Figure 3F and 3H). While lactic acidosis caused 57% repression in MCF7 cells transfected with non-targeting siRNA, such repression was reduced to 44% and 40% with MCF7 cells transfected with two different siRNAs targeting MondoA (M1 and M2). The degree of the lactic acidosis-induced repression of the glucose uptake was lower after the silencing of MondoA by two independent siRNAs (Figure S10).
Although gene silencing led to significant reduction of MondoA, the remaining MondoA level may still account for the slight induction of TXNIP under lactic acidosis. To further examine the effects in the absence of MondoA, we tested the MondoA knockout MEF (mouse embryonic fibroblasts) cells created by the cre-loxP system (CWP and DEA, manuscript submitted). In the complete absence of MondoA, both the basal and inducible level of TXNIP under 2-DG or lactic acidosis was completed abolished (Figure 5E). Importantly, the re-introduction of full length MondoA back to the MEF cells fully restored the induction of TXNIP under lactic acidosis (Figure 5E). Taken together, these results demonstrate the critical role of MondoA in regulating TXNIP under both glucose exposure and lactic acidosis.
TXNIP expression and lactic acidosis pathways in human cancers
To determine the prognostic significance of TXNIP expression in human cancers, we performed survival analyses using TXNIP expression as the sole predictor of survival time. Breast cancers in the Miller dataset stratified based on TXNIP expression were found to have significant differences in clinical phenotypes; tumors with high TXNIP have better survival and clinical outcomes (Figure 6A, p = 0.00115). Similar results were also obtained in the three other breast cancer datasets. To test whether the in vitro correlation of TXNIP induction with lactic acidosis pathway activity persists in vivo, we compared the predicted lactic acidosis pathway activity with the TXNIP expression levels in the breast cancer data sets and found a positive correlation in all four data sets; although the relationship is of limited predictive value with low R values reflecting a noisy relationship, the positive relationship is statistically significant and consistent across tumor data sets (Figure 6B). Similar results were found for ARRDC4, another paralouge of TXNIP and downstream target of MondoA. When breast cancers were stratified by the level of ARRDC4, tumors with higher levels of ARRDC4 had better prognosis and clinical outcome in the both breast cancer data (Miller and Pawitan) in which ARRDC4 was measured (Figure 6C). The expression of ARRDC4 was also positively associated with the predicted lactic acidosis pathway activity in these two breast cancer datasets (Figure 6D). These observations demonstrate the high degree of in vivo correlation between the expression of TXNIP and ARRDC4 and the lactic acidosis pathway activity in human cancers. It is also consistent with the known role of TXNIP as a probable tumor suppressor in several cancer types where it has been shown to suppress oncogenic phenotypes and its loss of function linked with cancer development [59]–[61]. Therefore, the expression of both TXNIP and ARRDC4 are potentially important mediators of the lactic acidosis response and suggest the important prognostic significance of molecular pathways driven by Mondo-Mlx in human cancers.
Discussion
Hypoxia, glucose deprivation and lactic acidosis are all well-recognized tumor microenvironmental stresses. Each of these stresses has distinct features – either limited availability of energy fuel (glucose), depletion of cofactor (oxygen) or accumulation of metabolic byproducts (lactic acidosis). Previous studies in various settings have shown that glucose deprivation causes metabolic stress and induces adaptive responses —so-called “starvation response” – manifest by AMPK activation, mTOR inhibition and multiple other biochemical and metabolic changes in cancer cells. Tumor hypoxia has also been reported to induce such metabolic stresses [62]. Here, we have shown that lactic acidosis triggers a “starvation” response similar to glucose deprivation as evidenced by the biochemical and transcriptional responses. Since this occurs even in the presence of abundant oxygen and nutrients, this response may represent “pseudo-starvation”. This is especially unexpected given the presence of high levels of lactate as an additional potential energy source for cancer cells [63].
Adaptive changes during the starvation response often reflect a need to preserve energy homeostasis during energy depletion by switching off cell growth and proliferation and many anabolic processes, such as protein, carbohydrate and lipid biosynthesis. In addition, there is also a simultaneous increase in cellular catabolism and breakdown of various energy sources to increase energy production. Even with this shared energy need, different tumor microenvironmental stresses employ different means to accomplish such adaptations to preserve energy homesostasis. During hypoxia, the essential cofactor (oxygen) required for oxidative phosphorylation becomes limited. As a result, the transcriptional factor HIF-1α activates genes encoding glucose transporters and enzymes involved in the glycolytic process shift to the use of glycolysis as the main source of energy generation. In the setting of glucose deprivation, there is a lack of fuel required for glycolysis as energy source. During these energy shortages, AMPK become activated to modify many proteins controlling energy flow to increase the energy generation in mitochondria from other source of fuels, such as oxidative of fatty acid and amino acids [64]. In addition, the inhibition of mTOR turns off many energy-consumption processes (e.g., translation, cell growth and proliferations) in an effort to restore energy homeostasis. The lactic acidosis-induced cell cycle arrest may be caused by many changes, including the induction of AMPK, mTOR inhibition, TXNIP and cell cycle arrest caused by the induction p57, p21 and other inhibitors of cellular proliferation. Moreover, the expression of many genes involved in glycolysis is affected by changes in histone acetylation of chromatin [65]. These responses triggered by glucose deprivation may be also important for energy maintenance and survival in acidosis-induced apoptosis [66], [67]. For example, AMPK may help in re-directing cells to utilize non-glucose energy sources (e.g., fatty acids and amino acids) and increased mitochondria activities under lactic acidosis [18], [68]. In addition, since the AMPK activator AICAR reduces the expression of lactate importer monocarboxylate transporter (MCT)-1 but increases the expression of the lactate exporter in MCT4 [69], AMPK activation under lactic acidosis may also inhibit the uptake of excessive cellular lactate. These stress environments are also likely to select for tumor cells that have developed strategies to survive energy deprivation with better energetic balance. The genetic mutations of many genes involved in these processes of adaptation to hypoxia and glucose deprivation are known to associate with tumor development. These findings highlight the crucial role of biochemical, metabolic control of energy homeostasis in tumor development.
There are many possible explanations as to how lactic acidosis leads to AMPK activation, mTOR inhibition and other features of the starvation response. For example, lactic acidosis can reduce ATP generation in cells [18] and thus lead to a high AMP/ATP ratio to activate the AMPK. Extracellular acidosis also triggers an increase in cytosolic Ca++ [9] which may activate CaMKK – an alternative AMPK upstream kinase [70]–[72] – to phophorylate and activate AMPK. Lactic acidosis may also cause intracellular nutrient depletion with reduced uptake of glucose (this study) and glutamine through the inhibition of acidosis-sensitive glutamine pumps [73], [74]. This depletion of intracellular pools of nutrients may in turn repress the energy sensor mTOR, the translation activities and ribosomal biogenesis required for cellular proliferation. In the future, it will be important to further dissect the contribution of each factor and signaling components to the induction of these starvation responses under lactic acidosis to further our detailed understanding of their impact on the metabolisms and phenotypes of cancer cells.
In both hypoxia and glucose deprivation, there is an increase in glucose uptake and glycolysis to provide essential fuel through induction (by HIFs), modification of glucose transporters (by AMPK) [75], [76] and histone acetylation via ATP-citrate lyase [65]. In contrast, lactic acidosis presents a different instance of starvation with inhibition of glucose uptake and other glycolysis activities in cancer cells [77]. Through these comparisons, we have also dissected two distinct molecular pathways (AMPK-mTORC1, MondoA-TXNIP) by which various microenvironmental stresses influence cancer metabolic phenotypes. While the AMPK-mTOR response is similar under lactic acidosis and glucose deprivation, the MondoA-TXNIP is affected in opposite directions by these two stresses. Given the continuous need for energy generation with reduced ATP and metabolic substrate from the glycolysis pathways under lactic acidosis, cells are likely to undergo extensive metabolic reprogramming to utilize other nutrients as energy sources. This idea is supported by the increased reliance on mitochondria for ATP generation. Our analysis has highlighted the potential roles of TXNIP and AMPK in this metabolic reprogramming. The induction of TXNIP under lactic acidosis and its ability to inhibit glucose uptake and reduce lactic acidosis production from glycolysis form a negative feedback loop. In addition, the loss of TXNIP leads to many features of Warburg effects and glycolytic phenotypes of cancer cells, opposite to the influences of lactic acidosis [18]. High TXNIP expression is associated with favorable outcomes, consistent with its postulated role as a tumor suppressor gene based on its growth-suppressing activity and the increased occurrences of tumors with deficiency of TXNIP [59], [60], [78], [79]. These findings may provide important insights into the regulatory mechanisms of TXNIP as well as the phenotypic alterations under lactic acidosis. TXNIP may affect the glucose uptake through the suppression of two important regulators of glycolysis, Akt [56] and HIF-1α [80]. In addition, TXNIP is known for other properties, which may explain its effect on the gene expression during lactic acidosis. TXNIP is a negative regulator of cellular oxidative tolerance by binding thioredoxin [81]–[83] and as a feedback regulator of S-nitrosylation [84] relevant in the cellular adaptive response to tumor microenvironmental stresses. For example, repression of TXNIP enhancing the anti-oxidative capacity during glucose deprivation may be required to cope with increased mitochondria oxidative stresses [85]. 2-Deoxy-D-glucose (2-DG) can enhance the cytotoxicity of cisplatin through mechanisms involving increased oxidative stress [86]. Since 2-DG is a strong inducer of TXNIP [57], a lower anti-oxidative capacity caused by 2-DG exposure may also involve neutralization of the anti-oxidative capacity of thioredoxin by TXNIP induction.
TXNIP and ARRDC4 are both transcriptionally regulated by the MondoA:Mlx complex to coordinate the fuel status of cellular metabolism and proliferation [57], [87]. Our findings identify lactic acidosis as a novel stimulus for the activation of MondoA and induction of TXNIP/ARRDC4. Even though we have identified the importance of MondoA-TXNIP of one component of the lactic acidosis response, there are still many aspects of lactic acidosis which remained unexplained, such as similarities to the glucose deprivation response. At least two factors are likely to be relevant for the shared gene expression response of glucose deprivation and lactic acidosis. First, these changes may lead to significant changes in the level of acetyl-CoA, which in turn impacts on histone acetylation of the chromatin of many target genes [65]. In addition, AMPK activation under both stresses may also impact on several transcription factors and co-activators to affect gene expression [88], [89]. It is likely that the anti-tumor activities of AMPK activators (e.g., metformin or AICAR) [90]–[93] may involve similar pathways triggered by lactic acidosis.
Many studies on the lactate levels in human cancers have found that tumors with high lactate levels are associated with poorer clinical outcomes, tumor aggression and treatment failure [3]–[5]. How do our finding of the association of the lactic acidosis response with favorable outcomes reconcile with these studies? There are at least four different ways how our results can be reconciled with the clinical observation. First, our studies mainly focus on the cellular response of cancer cells to short term exposure of lactic acidosis. This result may be different from the long term selection in the high lactic acidosis environments of human tumors. For example, glucose deprivation triggers starvation response with AMPK activation, mTOR inhibition and cell cycle arrest during short term exposure [48]. During the long term exposure to these stresses in human cancers, these undesirable conditions may select for tumor cells with the somatic mutations (such as K-Ras mutation [94]) that confers ability to adapt these stress conditions and strong metastasis potential. Similar selection pressure has also been suggested for the lactic acidosis [1], [15], [95]. In addition, the high level of lactate in the tumors is the downstream effects of the preferential use of glycolysis pathways due to tumor hypoxia, which may lead to more aggressive tumor behaviors and worse clinical outcomes. While lactic acidosis itself may exert some anti-tumor influences, this may not be enough to counter the influence of these somatic mutations and hypoxia in driving the tumor aggressiveness. It is also important to point out that the experimental evidence suggests that the lactic acidosis response is mainly due to acidosis instead of lactosis. Although we have found that lactosis may further augment the acidosis response, lactosis by itself have relatively little effects on the gene expression. Although the degree of the cellular response to these stresses (such as CA9) can reflect the levels of stresses (low tumor pO2) and serve as “endogenous markers” of such stresses in human cancers, such connection are not absolute and direct [96]. Therefore, the degree of lactic acidosis response in human cancers may or may not correlate directly with the tumor lactate levels.
Through the dissection of individual stresses in vitro, we have shown here that lactic acidosis simultaneously triggers two anti-tumor pathways (AMPK-mTOR and MondoA-TXNIP). Therefore, triggering such responses in cancer cells using small molecule compounds may have therapeutic potentials. It is of interest to investigate the mechanisms by which the lactic acidosis response is sensed and triggered in cancer cells. Extracellular lactate enters the cells through the lactate transporter proteins of monocarboxylate transporter (MCT), which may be also important for the lactic acidosis response. Lactate has also been recognized for its role as an energy source [63] and a signaling molecule to affect tumor cell phenotypes and target of cancer therapeutics [97]. Previous studies have shown that lowering the extracellular pH from 7.4 to ∼6.7 will lead to a slight lowering of intracellular pH (pHi) from 7.4 to 6.9–7.0 [98], [99]. Since many surface or cytosolic molecules exhibit high sensitivity to pH, this drop in pHe and the corresponding slight decrease in pHi may induce conformational changes to trigger signaling events [100], [101]. For example, it has been postulated that inhibition of the acidosis-sensing glutamine pumps leads to amino acid depletion [73], [74] and increases in the proteosome activity [102] in muscle cells, also contributing to muscle loss during metabolic acidosis. Low pHi has also been shown to directly enhance both DNA-binding and the TBP-interacting capacity of the transcriptional factor Sp1 [103]. Similarly, it is possible that this slight drop in pHi may destabilize the hydrogen bonds in several histone residues in the basic regions of MondoA to causes its transcriptional activation of target genes [57]. Homeostasis of pHi is mainly regulated by the sodium/hydrogen exchanger (NHE) family of proteins [104], [105] and such regulation is dysregulated in cancer cells [106]. TXNIP is known to be transcriptionally regulated by a variety of stresses and stimuli [84], [87], [107] and it is possible that some of these stimuli may act through modulating pHi to affect MondoA. TXNIP has been proposed as an inhibitor of cell growth [107] and the inhibition of the growth-inhibiting TXNIP expression under alkaline pHi may help to explain its permissiveness for cellular proliferation in response to growth factors [108], [109]. In addition, evidence is accumulating for the role of membrane acid-sensing receptors in either GPR4 family of G-protein coupled receptors (GPCR) and Acid-Sensing Ion Channels (ASICs) in many cell types. These acid-sensing receptors may also be critical in the observed response in both the AMPK-mTOR and MondoA-TXNIP under lactic acidosis. Importantly, these acid-sensing receptors can be modulated by small compounds pharmacologically [110], [111] and thus have potential as cancer therapeutics [112], [113]. For instance, several glycolipids are known natural ligands of these acid-sensing GPCRs and can activate or inhibit receptor responses to acidosis. Small compounds blocking the ASICs have been used to improve the symptoms and severity of stroke after vascular blockage [114]. Similarly, the use of compounds modulating the lactic acidosis response may mimic the tumor suppressor activities of lactic acidosis in novel cancer therapeutics. The definition of key genes and pathways will allow the use of genetic and chemical means to identify how to modulate cellular lactic acidosis for therapeutic purposes and thus potentially improve outcomes for cancer patients.
Materials and Methods
Cell culture and conditions modeling different tumor microenvironmental stresses
MCF7 breast cancer cell lines were cultured in DMEM (GIBCO11995) with 4.5 g/L glucose, supplemented with 10% fetal bovine serum, 1× non-essential amino acid and 1× antibiotics (penicillin, 10000UI/ml; streptomycin, 10000UI/ml). MEF cells were cultured in DMEM media supplemented with 15% FBS. Lactic acidosis conditions were created with the addition of 25mM or 10mM lactic acid (Sigma) to respective media and adjusted to the desired pHs with HCL or NaOH. Glucose deprivation was created by using 0 glucose/L media (GIBCO11966). Hypoxia was created by lowering the oxygen level to 1% oxygen. To stabilize the pH of DMEM media better, 25mM HEPES was also added into the media.
RNA isolation and microarray analysis
RNAs from MCF7 cells exposed to control or lactic acidosis culture conditions for 1, 4, 12, 24 hr were collected at respective time points and extracted with miRVana kits (Ambion), followed by hybridization to Affymetrix Hu133 plus2 gene chips with standard protocols. RNAs from MCF7 cells exposed to lactic acidosis, glucose deprivation and hypoxia were extracted after 4 hour exposure with miRVana kits (Ambion) and hybridized to Affymetrix Hu133 plus2 gene chips in a similar fashion. CEL files were normalized by RMA using Expression console (Affymetrix), filtered by indicated criteria, clustered with cluster 3.0, and displayed with treeview. All the microarray results have been submitted into GEO with accession number GSE19123. The RNA from the treated MEF cells were interrogated by hybridization to Affymetrix mouse 430A2 GeneChip with standard protocols and processed in a similar fashion.
Statistical analyses
Gene expression signatures associated with 1,4,12, and 24 hour lactic acidosis treatments were derived from RMA expression values using standard Bayesian sparse multivariate regression techniques, full details of which appear in previous publications [18], [115]–[117] For each treatment group (i.e. 1-,4-,12-, and 24-hr lactic acidosis treatment), this estimated the probability of differential expression for each Affymetrix probe as well as the fold change of the differentially expressed transcripts, as compared to the control. These values define the in-vitro signature of the treatment. Projection of these signatures into the primary tumor data sets is derived from the weighted inner product of the vector of tumor gene expression values and the vector of regression coefficients associated with the signature, as described in [117]: the score associated with signature k in sample i is defined as , where is the estimated probability of differential expression of probe g in treatment k, is the estimated fold change of the transcript, is the estimated residual variance of probe g, and is the RMA normalized expression value of probe g in sample i. The signature score thus provides a relative measure of the extent to which the pattern of expression described by the signature is present or reversed in a given vector of expression values.
P-values associated with the significance of Kaplan Meier curves were derived from survival analyses conducted using a Cox proportional hazards model in MATLAB. Kaplan Meier curves show the optimal stratification of high-risk and low-risk groups as determined by choosing the partition threshold that maximizes the area between curves. Measures of correlation and the associated P-values for the signature score scatter plots were calculated using the Pearson coefficient.
Full details of statistical analysis, including data and computer code to replicate the analyses, are available as Dataset S1.
Realtime RT–PCR
RNAs were reverse-transcribed to cDNAs with SuperScript II reverse transcription kit following the manufacturer's protocol (Invitrogen). cDNAs were then used as the substrate for gene expression level measured by qPCR with Power SYBRGreen PCR Mix (Applied Biosystem) and primers specific for TXNIP (Forward: CTGGCGTAAGCTTTTCAAGG, Reverse: AGTGCACAAAGGGGAAACAC), ARRDC4 (Forward: CCCCCTCCCACATGGTCACA, Reverse: TCCCTGGCTCCCTTCCATGTGT), Actin-beta (Forward: CTCTTCCAGCCTTCCTTCCT, Reverse: AGCACTGTGTTGGCGTACAG) following the manufacturer's protocol.
Glucose uptake assays
MCF7 cells were plated in 6-well/12-well plates at the density of 800,000/200,000cells per well. Once cells were more than 75% confluent, they were washed with 1× PBS twice, followed by application of serum starvation media (0.1%FBS) for three hours. Cells were then treated with respective conditions for the desired time. For MEF cells, they were plated in 12-well plates at the density of 100,000 cells per well. Once they reached more than 70% confluence, applying respective indicated conditions for four hours. For glucose uptake measurement, cells would then be washed with 37°C KRH buffer twice, followed by adding in 500ul/200ul KRH buffer containing 0.5uCi/0.2uCi 2-deoxy-D-glucose (GE Healthcare) for one hour in 37°C incubator. 20uM cytochalasin B was added for negative controls. After incubation, cells would be washed three times with 1ml/400ul of ice-cold KRH buffer containing 20mM glucose and 0.5mM phloretin to quench the glucose uptake. Finally, cells were lysed with 1ml/400ul RIPA buffer and the lysates were subjected to liquid scintillation counting. Protein concentrations were measured with Bradford assay. To measure the glucose uptake of genetically-manipulated cells, glucose uptake was measured 24 or 48 hours after the transfections.
RNA interference
MCF7 cells were plated in 12-well plates for the density of 200,000 cells per well. Once the cells reached 60% confluence, 100nM siRNAs were transfected by using lipofectamine. To verify the successful knocking down of the intended transcripts, RNAs were collected with miRVana kit 24 hours after the transfection, whereas 48 hours after the transfection, proteins were collected by lysing the cells with RIPA buffer.
Western blot analysis
Proteins were collected with RIPA buffer and their concentrations were measured with Bradford assay. Equal amounts of proteins were loaded for the protein analyses. Primary antibodies of AMPK, S6K (cell signaling), TXNIP (MBL) and anti-MondoA antibodies were applied following the manufacturers' protocols or as described previously [57].
Cell proliferation assay
Wild type mouse embryo fibroblast (MEF) and TXNIP null MEF cells were plated at the density of 25,000 cells per ml. The multi-channel pipette was used to plate 100ul of cell suspension evenly into the 96-well cell culture plate. Respective conditions of 10mM or 25mM lactic acidosis were applied the next day and cell number was estimation by using a standard MTS (3-(4,5-dimethylthiazol-2-yl)-5-(3-carboxymethoxyphenyl)-2-(4-sulphophenyl)-2H-tetrazolium)colorimetric assay (Promega) 48 hours after respective treatments.
Chromatin Immunoprecipitation
ChIP studies were performed as described previously [57] using an off-target region located on chromosome 10 to calculate fold enrichment. Primer sequences for the promoters of TXNIP, ARRDC4, and the off target control of the CHIP analysis are available upon request to D.E.A.
Glucose and lactate measurement in media
MCF cells were plated in six-well plates at a density of 500,000 cells per well. The next day, non-targeting siRNAs and siRNAs against TXNIP were transfected as described previously. After 24 hours, fresh media for the respective conditions, including control and 25mM lactic acidosis were applied to cells for 48 hrs when media were collected for glucose (ACCU-CHECK Aviva, Roche) and lactate (ARKRAY) measurements with respective meters. The results were normalized by cell numbers to obtain the glucose consumption and the lactate production amount per million cells.
Supporting Information
Zdroje
1. VaupelP
2004 Tumor microenvironmental physiology and its implications for radiation oncology. Semin Radiat Oncol 14 198 206
2. WarburgO
1956 On respiratory impairment in cancer cells. Science 124 269 270
3. QuennetV
YarominaA
ZipsD
RosnerA
WalentaS
2006 Tumor lactate content predicts for response to fractionated irradiation of human squamous cell carcinomas in nude mice. Radiother Oncol 81 130 135
4. BrizelDM
SchroederT
ScherRL
WalentaS
CloughRW
2001 Elevated tumor lactate concentrations predict for an increased risk of metastases in head-and-neck cancer. Int J Radiat Oncol Biol Phys 51 349 353
5. WalentaS
WetterlingM
LehrkeM
SchwickertG
SundforK
2000 High lactate levels predict likelihood of metastases, tumor recurrence, and restricted patient survival in human cervical cancers. Cancer Res 60 916 921
6. FangJS
GilliesRD
GatenbyRA
2008 Adaptation to hypoxia and acidosis in carcinogenesis and tumor progression. Semin Cancer Biol 18 330 337
7. GatenbyRA
GilliesRJ
2008 A microenvironmental model of carcinogenesis. Nat Rev Cancer 8 56 61
8. LaconiE
2007 The evolving concept of tumor microenvironments. Bioessays 29 738 744
9. HuangWC
SwietachP
Vaughan-JonesRD
AnsorgeO
GlitschMD
2008 Extracellular acidification elicits spatially and temporally distinct Ca2+ signals. Curr Biol 18 781 785
10. XuL
FidlerIJ
2000 Acidic pH-induced elevation in interleukin 8 expression by human ovarian carcinoma cells. Cancer Res 60 4610 4616
11. FukumuraD
XuL
ChenY
GohongiT
SeedB
2001 Hypoxia and acidosis independently up-regulate vascular endothelial growth factor transcription in brain tumors in vivo. Cancer Res 61 6020 6024
12. ShiQ
LeX
WangB
AbbruzzeseJL
XiongQ
2001 Regulation of vascular endothelial growth factor expression by acidosis in human cancer cells. Oncogene 20 3751 3756
13. MekhailK
GunaratnamL
BonicalziME
LeeS
2004 HIF activation by pH-dependent nucleolar sequestration of VHL. Nat Cell Biol 6 642 647
14. GrahamRM
FrazierDP
ThompsonJW
HalikoS
LiH
2004 A unique pathway of cardiac myocyte death caused by hypoxia-acidosis. J Exp Biol 207 3189 3200
15. MoelleringRE
BlackKC
KrishnamurtyC
BaggettBK
StaffordP
2008 Acid treatment of melanoma cells selects for invasive phenotypes. Clin Exp Metastasis 25 411 425
16. ZiekerD
SchaferR
GlatzleJ
NieseltK
CoerperS
2008 Lactate modulates gene expression in human mesenchymal stem cells. Langenbecks Arch Surg 393 297 301
17. NowikM
LeccaMR
VelicA
RehrauerH
BrandliAW
2008 Genome-wide gene expression profiling reveals renal genes regulated during metabolic acidosis. Physiol Genomics 32 322 334
18. ChenJL
LucasJE
SchroederT
MoriS
WuJ
2008 The genomic analysis of lactic acidosis and acidosis response in human cancers. PLoS Genet 4 e1000293 doi:10.1371/journal.pgen.1000293
19. ChiJT
WangZ
NuytenDS
RodriguezEH
SchanerME
2006 Gene Expression Programs in Response to Hypoxia: Cell Type Specificity and Prognostic Significance in Human Cancers. PLoS Med 3 e47 doi:10.1371/journal.pmed.0030370
20. ChiJ-T
RodriguezEH
WangZ
NuytenDSA
MukherjeeS
2007 Gene Expression Programs of Human Smooth Muscle Cells: Tissue-Specific Differentiation and Prognostic Significance in Breast Cancers. PLoS Genetics 3 e164 doi:10.1371/journal.pgen.0030164
21. BildAH
YaoG
ChangJT
WangQ
PottiA
2006 Oncogenic pathway signatures in human cancers as a guide to targeted therapies. Nature 439 353 357
22. ChangHY
SneddonJB
AlizadehAA
SoodR
WestRB
2004 Gene Expression Signature of Fibroblast Serum Response Predicts Human Cancer Progression: Similarities between Tumors and Wounds. PLoS Biol 2 e7 doi:10.1371/journal.pbio.0020007
23. GiacciaAJ
SimonMC
JohnsonR
2004 The biology of hypoxia: the role of oxygen sensing in development, normal function, and disease. Genes Dev 18 2183 2194
24. HarrisAL
2002 Hypoxia–a key regulatory factor in tumour growth. Nat Rev Cancer 2 38 47
25. SchofieldCJ
RatcliffePJ
2004 Oxygen sensing by HIF hydroxylases. Nat Rev Mol Cell Biol 5 343 354
26. SemenzaGL
2002 HIF-1 and tumor progression: pathophysiology and therapeutics. Trends Mol Med 8 S62 67
27. MotzerRJ
BanderNH
NanusDM
1996 Renal-cell carcinoma. N Engl J Med 335 865 875
28. IsaacsJS
JungYJ
MoleDR
LeeS
Torres-CabalaC
2005 HIF overexpression correlates with biallelic loss of fumarate hydratase in renal cancer: novel role of fumarate in regulation of HIF stability. Cancer Cell 8 143 153
29. PollardPJ
Spencer-DeneB
ShuklaD
HowarthK
NyeE
2007 Targeted inactivation of fh1 causes proliferative renal cyst development and activation of the hypoxia pathway. Cancer Cell 11 311 319
30. SelakMA
ArmourSM
MacKenzieED
BoulahbelH
WatsonDG
2005 Succinate links TCA cycle dysfunction to oncogenesis by inhibiting HIF-alpha prolyl hydroxylase. Cancer Cell 7 77 85
31. SchroederT
YuanH
VigliantiBL
PeltzC
AsopaS
2005 Spatial heterogeneity and oxygen dependence of glucose consumption in R3230Ac and fibrosarcomas of the Fischer 344 rat. Cancer Res 65 5163 5171
32. ShawRJ
KosmatkaM
BardeesyN
HurleyRL
WittersLA
2004 The tumor suppressor LKB1 kinase directly activates AMP-activated kinase and regulates apoptosis in response to energy stress. Proc Natl Acad Sci U S A 101 3329 3335
33. YooLI
ChungDC
YuanJ
2002 LKB1–a master tumour suppressor of the small intestine and beyond. Nat Rev Cancer 2 529 535
34. MaL
ChenZ
Erdjument-BromageH
TempstP
PandolfiPP
2005 Phosphorylation and functional inactivation of TSC2 by Erk implications for tuberous sclerosis and cancer pathogenesis. Cell 121 179 193
35. PanD
DongJ
ZhangY
GaoX
2004 Tuberous sclerosis complex: from Drosophila to human disease. Trends Cell Biol 14 78 85
36. EllisenLW
RamsayerKD
JohannessenCM
YangA
BeppuH
2002 REDD1, a developmentally regulated transcriptional target of p63 and p53, links p63 to regulation of reactive oxygen species. Mol Cell 10 995 1005
37. SternR
ShusterS
NeudeckerBA
FormbyB
2002 Lactate stimulates fibroblast expression of hyaluronan and CD44: the Warburg effect revisited. Exp Cell Res 276 24 31
38. SubramanianA
TamayoP
MoothaVK
MukherjeeS
EbertBL
2005 Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A 102 15545 15550
39. PengT
GolubTR
SabatiniDM
2002 The immunosuppressant rapamycin mimics a starvation-like signal distinct from amino acid and glucose deprivation. Mol Cell Biol 22 5575 5584
40. van 't VeerLJ
DaiH
van de VijverMJ
HeYD
HartAA
2002 Gene expression profiling predicts clinical outcome of breast cancer. Nature 415 530 536
41. MillerLD
SmedsJ
GeorgeJ
VegaVB
VergaraL
2005 An expression signature for p53 status in human breast cancer predicts mutation status, transcriptional effects, and patient survival. Proc Natl Acad Sci U S A 102 13550 13555
42. WangY
KlijnJG
ZhangY
SieuwertsAM
LookMP
2005 Gene-expression profiles to predict distant metastasis of lymph-node-negative primary breast cancer. Lancet 365 671 679
43. SotiriouC
WirapatiP
LoiS
HarrisA
FoxS
2006 Gene expression profiling in breast cancer: understanding the molecular basis of histologic grade to improve prognosis. J Natl Cancer Inst 98 262 272
44. PawitanY
BjohleJ
AmlerL
BorgAL
EgyhaziS
2005 Gene expression profiling spares early breast cancer patients from adjuvant therapy: derived and validated in two population-based cohorts. Breast Cancer Res 7 R953 964
45. JangA
HillRP
1997 An examination of the effects of hypoxia, acidosis, and glucose starvation on the expression of metastasis-associated genes in murine tumor cells. Clin Exp Metastasis 15 469 483
46. CuvierC
JangA
HillRP
1997 Exposure to hypoxia, glucose starvation and acidosis: effect on invasive capacity of murine tumor cells and correlation with cathepsin (L+B) secretion. Clin Exp Metastasis 15 19 25
47. GaschAP
SpellmanPT
KaoCM
Carmel-HarelO
EisenMB
2000 Genomic expression programs in the response of yeast cells to environmental changes. Mol Biol Cell 11 4241 4257
48. ShackelfordDB
ShawRJ
2009 The LKB1-AMPK pathway: metabolism and growth control in tumour suppression. Nat Rev Cancer 9 563 575
49. HayN
SonenbergN
2004 Upstream and downstream of mTOR. Genes Dev 18 1926 1945
50. HayashiT
HirshmanMF
KurthEJ
WinderWW
GoodyearLJ
1998 Evidence for 5′ AMP-activated protein kinase mediation of the effect of muscle contraction on glucose transport. Diabetes 47 1369 1373
51. BakerAF
KohMY
WilliamsRR
JamesB
WangH
2008 Identification of thioredoxin-interacting protein 1 as a hypoxia-inducible factor 1alpha-induced gene in pancreatic cancer. Pancreas 36 178 186
52. SchulzePC
YoshiokaJ
TakahashiT
HeZ
KingGL
2004 Hyperglycemia promotes oxidative stress through inhibition of thioredoxin function by thioredoxin-interacting protein. J Biol Chem 279 30369 30374
53. MinnAH
HafeleC
ShalevA
2005 Thioredoxin-interacting protein is stimulated by glucose through a carbohydrate response element and induces beta-cell apoptosis. Endocrinology 146 2397 2405
54. MuoioDM
2007 TXNIP links redox circuitry to glucose control. Cell Metab 5 412 414
55. ParikhH
CarlssonE
ChutkowWA
JohanssonLE
StorgaardH
2007 TXNIP regulates peripheral glucose metabolism in humans. PLoS Med 4 e158 doi:10.1371/journal.pmed.0040158
56. HuiST
AndresAM
MillerAK
SpannNJ
PotterDW
2008 Txnip balances metabolic and growth signaling via PTEN disulfide reduction. Proc Natl Acad Sci U S A 105 3921 3926
57. StoltzmanCA
PetersonCW
BreenKT
MuoioDM
BillinAN
2008 Glucose sensing by MondoA:Mlx complexes: a role for hexokinases and direct regulation of thioredoxin-interacting protein expression. Proc Natl Acad Sci U S A 105 6912 6917
58. Cha-MolstadH
SaxenaG
ChenJ
ShalevA
2009 Glucose-stimulated expression of Txnip is mediated by carbohydrate response element-binding protein, p300, and histone H4 acetylation in pancreatic beta cells. J Biol Chem 284 16898 16905
59. GoldbergSF
MieleME
HattaN
TakataM
Paquette-StraubC
2003 Melanoma metastasis suppression by chromosome 6: evidence for a pathway regulated by CRSP3 and TXNIP. Cancer Res 63 432 440
60. ShethSS
BodnarJS
GhazalpourA
ThipphavongCK
TsutsumiS
2006 Hepatocellular carcinoma in Txnip-deficient mice. Oncogene 25 3528 3536
61. HanSH
JeonJH
JuHR
JungU
KimKY
2003 VDUP1 upregulated by TGF-beta1 and 1,25-dihydorxyvitamin D3 inhibits tumor cell growth by blocking cell-cycle progression. Oncogene 22 4035 4046
62. LiuL
CashTP
JonesRG
KeithB
ThompsonCB
2006 Hypoxia-induced energy stress regulates mRNA translation and cell growth. Mol Cell 21 521 531
63. SonveauxP
VegranF
SchroederT
WerginMC
VerraxJ
2008 Targeting lactate-fueled respiration selectively kills hypoxic tumor cells in mice. J Clin Invest 118 3930 3942
64. HardieDG
2007 AMP-activated/SNF1 protein kinases: conserved guardians of cellular energy. Nat Rev Mol Cell Biol 8 774 785
65. WellenKE
HatzivassiliouG
SachdevaUM
BuiTV
CrossJR
2009 ATP-citrate lyase links cellular metabolism to histone acetylation. Science 324 1076 1080
66. Harhaji-TrajkovicL
VilimanovichU
Kravic-StevovicT
BumbasirevicV
TrajkovicV
2009 AMPK-mediated autophagy inhibits apoptosis in cisplatin-treated tumor cells. J Cell Mol Med
67. KubasiakLA
HernandezOM
BishopricNH
WebsterKA
2002 Hypoxia and acidosis activate cardiac myocyte death through the Bcl-2 family protein BNIP3. Proc Natl Acad Sci U S A 99 12825 12830
68. PutmanCT
KiricsiM
PearceyJ
MacLeanIM
BamfordJA
2003 AMPK activation increases uncoupling protein-3 expression and mitochondrial enzyme activities in rat muscle without fibre type transitions. J Physiol 551 169 178
69. GalardoMN
RieraMF
PellizzariEH
CigorragaSB
MeroniSB
2007 The AMP-activated protein kinase activator, 5-aminoimidazole-4-carboxamide-1-b-D-ribonucleoside, regulates lactate production in rat Sertoli cells. J Mol Endocrinol 39 279 288
70. HawleySA
PanDA
MustardKJ
RossL
BainJ
2005 Calmodulin-dependent protein kinase kinase-beta is an alternative upstream kinase for AMP-activated protein kinase. Cell Metab 2 9 19
71. WoodsA
DickersonK
HeathR
HongSP
MomcilovicM
2005 Ca2+/calmodulin-dependent protein kinase kinase-beta acts upstream of AMP-activated protein kinase in mammalian cells. Cell Metab 2 21 33
72. HurleyRL
AndersonKA
FranzoneJM
KempBE
MeansAR
2005 The Ca2+/calmodulin-dependent protein kinase kinases are AMP-activated protein kinase kinases. J Biol Chem 280 29060 29066
73. EvansK
NasimZ
BrownJ
ClappE
AminA
2008 Inhibition of SNAT2 by metabolic acidosis enhances proteolysis in skeletal muscle. J Am Soc Nephrol 19 2119 2129
74. EvansK
NasimZ
BrownJ
ButlerH
KauserS
2007 Acidosis-sensing glutamine pump SNAT2 determines amino acid levels and mammalian target of rapamycin signalling to protein synthesis in L6 muscle cells. J Am Soc Nephrol 18 1426 1436
75. RussellRR3rd
BergeronR
ShulmanGI
YoungLH
1999 Translocation of myocardial GLUT-4 and increased glucose uptake through activation of AMPK by AICAR. Am J Physiol 277 H643 649
76. BergeronR
RussellRR3rd
YoungLH
RenJM
MarcucciM
1999 Effect of AMPK activation on muscle glucose metabolism in conscious rats. Am J Physiol 276 E938 944
77. ChenJL
LucasJE
SchroederT
MoriS
NevinsJR
2008 Genomic analysis of response to lactic acidosis and acidosis in human cancers. PLoS Genetics 4 e1000293 doi:10.1371/journal.pgen.1000293
78. ErkelandSJ
PalandeKK
ValkhofM
GitsJ
Danen-van OorschotA
2009 The gene encoding thioredoxin-interacting protein (TXNIP) is a frequent virus integration site in virus-induced mouse leukemia and is overexpressed in a subset of AML patients. Leuk Res 33 1367 1371
79. NishinakaY
NishiyamaA
MasutaniH
OkaS
AhsanKM
2004 Loss of thioredoxin-binding protein-2/vitamin D3 up-regulated protein 1 in human T-cell leukemia virus type I-dependent T-cell transformation: implications for adult T-cell leukemia leukemogenesis. Cancer Res 64 1287 1292
80. ShinD
JeonJH
JeongM
SuhHW
KimS
2008 VDUP1 mediates nuclear export of HIF1alpha via CRM1-dependent pathway. Biochim Biophys Acta 1783 838 848
81. JunnE
HanSH
ImJY
YangY
ChoEW
2000 Vitamin D3 up-regulated protein 1 mediates oxidative stress via suppressing the thioredoxin function. J Immunol 164 6287 6295
82. YamanakaH
MaehiraF
OshiroM
AsatoT
YanagawaY
2000 A possible interaction of thioredoxin with VDUP1 in HeLa cells detected in a yeast two-hybrid system. Biochem Biophys Res Commun 271 796 800
83. NishiyamaA
MatsuiM
IwataS
HirotaK
MasutaniH
1999 Identification of thioredoxin-binding protein-2/vitamin D(3) up-regulated protein 1 as a negative regulator of thioredoxin function and expression. J Biol Chem 274 21645 21650
84. ForresterMT
SethD
HausladenA
EylerCE
FosterMW
2009 Thioredoxin interacting protein (Txnip) is a feedback regulator of S-nitrosylation. J Biol Chem
85. AhmadIM
Aykin-BurnsN
SimJE
WalshSA
HigashikuboR
2005 Mitochondrial O2*- and H2O2 mediate glucose deprivation-induced stress in human cancer cells. J Biol Chem 280 4254 4263
86. SimonsAL
AhmadIM
MattsonDM
DornfeldKJ
SpitzDR
2007 2-Deoxy-D-glucose combined with cisplatin enhances cytotoxicity via metabolic oxidative stress in human head and neck cancer cells. Cancer Res 67 3364 3370
87. KaadigeMR
LooperRE
KamalanaadhanS
AyerDE
2009 Glutamine-dependent anapleurosis dictates glucose uptake and cell growth by regulating MondoA transcriptional activity. Proc Natl Acad Sci U S A 106 14878 14883
88. YangW
HongYH
ShenXQ
FrankowskiC
CampHS
2001 Regulation of transcription by AMP-activated protein kinase: phosphorylation of p300 blocks its interaction with nuclear receptors. J Biol Chem 276 38341 38344
89. KawaguchiT
OsatomiK
YamashitaH
KabashimaT
UyedaK
2002 Mechanism for fatty acid “sparing” effect on glucose-induced transcription: regulation of carbohydrate-responsive element-binding protein by AMP-activated protein kinase. J Biol Chem 277 3829 3835
90. SwinnenJV
BeckersA
BrusselmansK
OrganeS
SegersJ
2005 Mimicry of a cellular low energy status blocks tumor cell anabolism and suppresses the malignant phenotype. Cancer Res 65 2441 2448
91. RattanR
GiriS
SinghAK
SinghI
2005 5-Aminoimidazole-4-carboxamide-1-beta-D-ribofuranoside inhibits cancer cell proliferation in vitro and in vivo via AMP-activated protein kinase. J Biol Chem 280 39582 39593
92. BuzzaiM
JonesRG
AmaravadiRK
LumJJ
DeBerardinisRJ
2007 Systemic treatment with the antidiabetic drug metformin selectively impairs p53-deficient tumor cell growth. Cancer Res 67 6745 6752
93. ZakikhaniM
DowlingR
FantusIG
SonenbergN
PollakM
2006 Metformin is an AMP kinase-dependent growth inhibitor for breast cancer cells. Cancer Res 66 10269 10273
94. YunJ
RagoC
CheongI
PagliariniR
AngenendtP
2009 Glucose deprivation contributes to the development of KRAS pathway mutations in tumor cells. Science 325 1555 1559
95. GatenbyRA
GilliesRJ
2004 Why do cancers have high aerobic glycolysis? Nat Rev Cancer 4 891 899
96. MayerA
HockerM
VaupelP
2008 Endogenous Hypoxia Markers: Case Not Proven;
VaupelP
Springer US
97. WalentaS
Mueller-KlieserWF
2004 Lactate: mirror and motor of tumor malignancy. Semin Radiat Oncol 14 267 274
98. AdamsDJ
WahlML
FlowersJL
SenB
ColvinM
2006 Camptothecin analogs with enhanced activity against human breast cancer cells. II. Impact of the tumor pH gradient. Cancer Chemother Pharmacol 57 145 154
99. WahlML
PoolerPM
BriandP
LeeperDB
OwenCS
2000 Intracellular pH regulation in a nonmalignant and a derived malignant human breast cell line. J Cell Physiol 183 373 380
100. BusaWB
NuccitelliR
1984 Metabolic regulation via intracellular pH. Am J Physiol 246 R409 438
101. ShrodeLD
TapperH
GrinsteinS
1997 Role of intracellular pH in proliferation, transformation, and apoptosis. J Bioenerg Biomembr 29 393 399
102. MitchWE
MedinaR
GrieberS
MayRC
EnglandBK
1994 Metabolic acidosis stimulates muscle protein degradation by activating the adenosine triphosphate-dependent pathway involving ubiquitin and proteasomes. J Clin Invest 93 2127 2133
103. TorigoeT
IzumiH
YoshidaY
IshiguchiH
OkamotoT
2003 Low pH enhances Sp1 DNA binding activity and interaction with TBP. Nucleic Acids Res 31 4523 4530
104. GrinsteinS
WoodsideM
WaddellTK
DowneyGP
OrlowskiJ
1993 Focal localization of the NHE-1 isoform of the Na+/H+ antiport: assessment of effects on intracellular pH. EMBO J 12 5209 5218
105. CounillonL
PouyssegurJ
2000 The expanding family of eucaryotic Na(+)/H(+) exchangers. J Biol Chem 275 1 4
106. CardoneRA
CasavolaV
ReshkinSJ
2005 The role of disturbed pH dynamics and the Na+/H+ exchanger in metastasis. Nat Rev Cancer 5 786 795
107. KaadigeMR
LooperRE
KamalanaadhanS
AyerDE
2009 Glutamine-dependent anapleurosis dictates glucose uptake and cell growth by regulating MondoA transcriptional activity. Proc Natl Acad Sci U S A
108. PouyssegurJ
FranchiA
L'AllemainG
ParisS
1985 Cytoplasmic pH, a key determinant of growth factor-induced DNA synthesis in quiescent fibroblasts. FEBS Lett 190 115 119
109. PutneyLK
BarberDL
2003 Na-H exchange-dependent increase in intracellular pH times G2/M entry and transition. J Biol Chem 278 44645 44649
110. TomuraH
MogiC
SatoK
OkajimaF
2005 Proton-sensing and lysolipid-sensitive G-protein-coupled receptors: a novel type of multi-functional receptors. Cell Signal 17 1466 1476
111. RaduCG
NijagalA
McLaughlinJ
WangL
WitteON
2005 Differential proton sensitivity of related G protein-coupled receptors T cell death-associated gene 8 and G2A expressed in immune cells. Proc Natl Acad Sci U S A 102 1632 1637
112. ImDS
2005 Two ligands for a GPCR, proton vs lysolipid. Acta Pharmacol Sin 26 1435 1441
113. MogiC
TomuraH
ToboM
WangJQ
DamirinA
2005 Sphingosylphosphorylcholine antagonizes proton-sensing ovarian cancer G-protein-coupled receptor 1 (OGR1)-mediated inositol phosphate production and cAMP accumulation. J Pharmacol Sci 99 160 167
114. XiongZG
ZhuXM
ChuXP
MinamiM
HeyJ
2004 Neuroprotection in ischemia: blocking calcium-permeable acid-sensing ion channels. Cell 118 687 698
115. CarvalhoC
ChangJ
LucasJ
WangQ
NevinsJR
2008 High-dimensional sparse factor modeling: Applications in gene expression genomics. Journal of American Statistical Society 103 1438 1456
116. LucasJE
CarvalhoCM
ChenJL
ChiJT
WestM
2009 Cross-study projections of genomic biomarkers: an evaluation in cancer genomics. PLoS ONE 4 e4523 doi:10.1371/journal.pone.0004523
117. MerlD
ChenJLY
ChiJT
WestM
2009 An integrative analysis of cancer gene expression studies using Bayesian latent factor modeling. Annals of Applied Statistics 3 1675 1694
Štítky
Genetika Reprodukční medicínaČlánek vyšel v časopise
PLOS Genetics
2010 Číslo 9
- Primární hyperoxalurie – aktuální možnosti diagnostiky a léčby
- Srdeční frekvence embrya může být faktorem užitečným v předpovídání výsledku IVF
- Akutní intermitentní porfyrie
- Vztah užívání alkoholu a mužské fertility
- Vliv kvality morfologie spermií na úspěšnost intrauterinní inseminace
Nejčtenější v tomto čísle
- Synthesizing and Salvaging NAD: Lessons Learned from
- Optimal Strategy for Competence Differentiation in Bacteria
- Long- and Short-Term Selective Forces on Malaria Parasite Genomes
- Identifying Signatures of Natural Selection in Tibetan and Andean Populations Using Dense Genome Scan Data