Epithelial to mesenchymal transition (EMT) has implications in tumor progression and metastasis. Metabolic alterations have been described in cancer development but studies focused on the metabolic re-wiring that takes place during EMT are still limited. We performed metabolomics profiling of a breast epithelial cell line and its EMT derived mesenchymal phenotype to create genome-scale metabolic models descriptive of both cell lines. Glycolysis and OXPHOS were higher in the epithelial phenotype while amino acid anaplerosis and fatty acid oxidation fueled the mesenchymal phenotype. Through comparative bioinformatics analysis, PPAR-γ1, PPAR- γ2 and AP-1 were found to be the most influential transcription factors associated with metabolic re-wiring. In silico gene essentiality analysis predicts that the LAT1 neutral amino acid transporter is essential for mesenchymal cell survival. Our results define metabolic traits that distinguish an EMT derived mesenchymal cell line from its epithelial progenitor and may have implications in cancer progression and metastasis. Furthermore, the tools presented here can aid in identifying critical metabolic nodes that may serve as therapeutic targets aiming to prevent EMT and inhibit metastatic dissemination.

Signaling between cancer and nonmalignant (stromal) cells in the tumor microenvironment (TME) is a key to tumor progression. Here, we deconvoluted bulk tumor transcriptomes to infer cross-talk between ligands and receptors on cancer and stromal cells in the TME of 20 solid tumor types. This approach recovered known transcriptional hallmarks of cancer and stromal cells and was concordant with single-cell, in situ hybridization and IHC data. Inferred autocrine cancer cell interactions varied between tissues but often converged on Ephrin, BMP, and FGFR-signaling pathways. Analysis of immune checkpoints nominated interactions with high levels of cancer-to-immune cross-talk across distinct tumor types. Strikingly, PD-L1 was found to be highly expressed in stromal rather than cancer cells. Overall, our study presents a new resource for hypothesis generation and exploration of cross-talk in the TME. Significance: This study provides deconvoluted bulk tumor transcriptomes across multiple cancer types to infer cross-talk in the tumor microenvironment.

The metabolism of gluconate is well characterized in prokaryotes where it is known to be degraded following phosphorylation by gluconokinase. Less is known of gluconate metabolism in humans. Human gluconokinase activity was recently identified proposing questions about the metabolic role of gluconate in humans. Here we report the recombinant expression, purification and biochemical characterization of isoform I of human gluconokinase alongside substrate specificity and kinetic assays of the enzyme catalyzed reaction. The enzyme, shown to be a dimer, had ATP dependent phosphorylation activity and strict specificity towards gluconate out of 122 substrates tested. In order to evaluate the metabolic impact of gluconate in humans we modeled gluconate metabolism using steady state metabolic network analysis. The results indicate that significant metabolic flux changes in anabolic pathways linked to the hexose monophosphate shunt (HMS) are induced through a small increase in gluconate concentration. We argue that the enzyme takes part in a context specific carbon flux route into the HMS that, in humans, remains incompletely explored. Apart from the biochemical description of human gluconokinase, the results highlight that little is known of the mechanism of gluconate metabolism in humans despite its widespread use in medicine and consumer products.

Epithelial to mesenchymal transition (EMT) is an important event during development and cancer metastasis. There is limited understanding of the metabolic alterations that give rise to and take place during EMT. Dysregulation of signalling pathways that impact metabolism, including epidermal growth factor receptor (EGFR), are however a hallmark of EMT and metastasis. In this study, we report the investigation into EGFR signalling and metabolic crosstalk of EMT through constraint-based modelling and analysis of the breast epithelial EMT cell model D492 and its mesenchymal counterpart D492M. We built an EGFR signalling network for EMT based on stoichiometric coefficients and constrained the network with gene expression data to build epithelial (EGFR_E) and mesenchymal (EGFR_M) networks. Metabolic alterations arising from differential expression of EGFR genes was derived from a literature review of AKT regulated metabolic genes. Signaling flux differences between EGFR_E and EGFR_M models subsequently allowed metabolism in D492 and D492M cells to be assessed. Higher flux within AKT pathway in the D492 cells compared to D492M suggested higher glycolytic activity in D492 that we confirmed experimentally through measurements of glucose uptake and lactate secretion rates. The signaling genes from the AKT, RAS/MAPK and CaM pathways were predicted to revert D492M to D492 phenotype. Follow-up analysis of EGFR signaling metabolic crosstalk in three additional breast epithelial cell lines highlighted variability in in vitro cell models of EMT. This study shows that the metabolic phenotype may be predicted by in silico analyses of gene expression data of EGFR signaling genes, but this phenomenon is cell-specific and does not follow a simple trend.