GENI: A web server to identify gene set enrichments in tumor samples.

The Cancer Genome Atlas (TCGA) and analogous projects have yielded invaluable tumor-associated genomic data. Despite several web-based platforms designed to enhance accessibility, certain analyses require prior bioinformatic expertise. To address this need, we developed Gene ENrichment Identifier (GENI, https://www.shaullab.com/geni), which is designed to promptly compute correlations for genes of interest against the entire transcriptome and rank them against well-established biological gene sets. Additionally, it generates comprehensive tables containing genes of interest and their corresponding correlation coefficients, presented in publication-quality graphs. Furthermore, GENI has the capability to analyze multiple genes simultaneously within a given gene set, elucidating their significance within a specific biological context. Overall, GENI's user-friendly interface simplifies the biological interpretation and analysis of cancer patient-associated data, advancing the understanding of cancer biology and accelerating scientific discoveries.

Depletion of Fumarate Hydratase, an Essential TCA Cycle Enzyme, Drives Proliferation in a Two-Step Model.

Fumarate hydratase (FH) is an evolutionary conserved TCA cycle enzyme that reversibly catalyzes the hydration of fumarate to L-malate and has a moonlight function in the DNA damage response (DDR). Interestingly, FH has a contradictory cellular function, as it is pro-survival through its role in the TCA cycle, yet its loss can drive tumorigenesis. Here, we found that in both non-cancerous (HEK-293T) and cancerous cell lines (HepG2), the cell response to FH loss is separated into two distinct time frames based on cell proliferation and DNA damage repair. During the early stages of FH loss, cell proliferation rate and DNA damage repair are inhibited. However, over time the cells overcome the FH loss and form knockout clones, indistinguishable from WT cells with respect to their proliferation rate. Due to the FH loss effect on DNA damage repair, we assumed that the recovered cells bear adaptive mutations. Therefore, we applied whole-exome sequencing to identify such mutated genes systematically. Indeed, we identified recurring mutations in genes belonging to central oncogenic signaling pathways, such as JAK/STAT3, which we validated in impaired FH-KO clones. Intriguingly, we demonstrate that these adaptive mutations are responsible for FH-KO cell proliferation under TCA cycle malfunction.

Overexpression of UBA5 in Cells Mimics the Phenotype of Cells Lacking UBA5.

Ufmylation is a posttranslational modification in which the modifier UFM1 is attached to target proteins. This conjugation requires the concerted work of three enzymes named UBA5, UFC1, and UFL1. Initially, UBA5 activates UFM1 in a process that ends with UFM1 attached to UBA5's active site Cys. Then, in a trans-thiolation reaction, UFM1 is transferred from UBA5 to UFC1, forming a thioester bond with the latter. Finally, with the help of UFL1, UFM1 is transferred to the final destination-a lysine residue on a target protein. Therefore, not surprisingly, deletion of one of these enzymes abrogates the conjugation process. However, how overexpression of these enzymes affects this process is not yet clear. Here we found, unexpectedly, that overexpression of UBA5, but not UFC1, damages the ability of cells to migrate, in a similar way to cells lacking UBA5 or UFC1. At the mechanistic level, we found that overexpression of UBA5 reverses the trans-thiolation reaction, thereby leading to a back transfer of UFM1 from UFC1 to UBA5. This, as seen in cells lacking UBA5, reduces the level of charged UFC1 and therefore harms the conjugation process. In contrast, co-expression of UBA5 with UFM1 abolishes this effect, suggesting that the reverse transfer of UFM1 from UFC1 to UBA5 depends on the level of free UFM1. Overall, our results propose that the cellular expression level of the UFM1 conjugation enzymes has to be tightly regulated to ensure the proper directionality of UFM1 transfer.

DPYSL2 interacts with JAK1 to mediate breast cancer cell migration.

The intricate neuronal wiring during development requires cytoskeletal reorganization orchestrated by signaling cues. Because cytoskeletal remodeling is a hallmark of cell migration, we investigated whether metastatic cancer cells exploit axon guidance proteins to migrate. Indeed, in breast cancer patients, we found a significant correlation between mesenchymal markers and the expression of dihydropyrimidinase-like 2 (DPYSL2), a regulator of cytoskeletal dynamics in growing axons. Strikingly, DPYSL2 knockout in mesenchymal-like breast cancer cells profoundly inhibited cell migration, invasion, stemness features, tumor growth rate, and metastasis. Next, we decoded the molecular mechanism underlying this phenomenon and revealed an interaction between DPYSL2 and Janus kinase 1 (JAK1). This binding is crucial for activating signal transducer and activator of transcription 3 (STAT3) and the subsequent expression of vimentin, the promigratory intermediate filament. These findings identify DPYSL2 as a molecular link between oncogenic signaling pathways and cytoskeletal reorganization in migrating breast cancer cells.

Monitoring Breast Cancer Growth and Metastatic Colony Formation in Mice using Bioluminescence.

Breast cancer is a frequent heterogeneous malignancy and the second leading cause of mortality in women, mainly due to distant organ metastasis. Several animal models have been generated, including the widely used orthotopic mouse models, where cancer cells are injected into the mammary fat pad. However, these models cannot help monitor tumor growth kinetics and metastatic colonization. Cutting-edge tools to monitor cancer cells in real time in mice will significantly advance the understanding of tumor biology. Here, breast cancer cell lines stably expressing luciferase and green fluorescent protein (GFP) were established. Specifically, this technique contains two sequential steps initiated by measuring the luciferase activity in vitro and followed by the implantation of the cancer cells into mammary fat pads of nonobese diabetic-severe combined immunodeficiency (NOD-SCID) mice. After the injection, both the tumor growth and metastatic colonization are monitored in real time by the noninvasive bioluminescence imaging system. Then, the quantification of GFP-expressing metastases in the lungs will be examined by fluorescence microscopy to validate the observed bioluminescence results. This sophisticated system combining luciferase and fluorescence-based detection tools evaluates cancer metastasis in vivo, which has great potential for use in breast cancer therapeutics and disease management.

The glutathione peroxidase 8 (GPX8)/IL-6/STAT3 axis is essential in maintaining an aggressive breast cancer phenotype.

One of the emerging hallmarks of cancer illustrates the importance of metabolic reprogramming, necessary to synthesize the building blocks required to fulfill the high demands of rapidly proliferating cells. However, the proliferation-independent instructive role of metabolic enzymes in tumor plasticity is still unclear. Here, we provide evidence that glutathione peroxidase 8 (GPX8), a poorly characterized enzyme that resides in the endoplasmic reticulum, is an essential regulator of tumor aggressiveness. We found that GPX8 expression was induced by the epithelial-mesenchymal transition (EMT) program. Moreover, in breast cancer patients, GPX8 expression significantly correlated with known mesenchymal markers and poor prognosis. Strikingly, GPX8 knockout in mesenchymal-like cells (MDA-MB-231) resulted in an epithelial-like morphology, down-regulation of EMT characteristics, and loss of cancer stemness features. In addition, GPX8 knockout significantly delayed tumor initiation and decreased its growth rate in mice. We found that these GPX8 loss-dependent phenotypes were accompanied by the repression of crucial autocrine factors, in particular, interleukin-6 (IL-6). In these cells, IL-6 bound to the soluble receptor (sIL6R), stimulating the JAK/STAT3 signaling pathway by IL-6 trans-signaling mechanisms, so promoting cancer aggressiveness. We observed that in GPX8 knockout cells, this signaling mechanism was impaired as sIL6R failed to activate the JAK/STAT3 signaling pathway. Altogether, we present the GPX8/IL-6/STAT3 axis as a metabolic-inflammatory pathway that acts as a robust regulator of cancer cell aggressiveness.

Large-Scale Differential Gene Expression Transcriptomic Analysis Identifies a Metabolic Signature Shared by All Cancer Cells.

Cancer-dependent metabolic rewiring is often manifested by selective expression of enzymes essential for the transformed cells' viability. However, the metabolic variations between normal and transformed cells are not fully characterized, and therefore, a systematic analysis will result in the identification of unknown cellular mechanisms crucial for tumorigenesis. Here, we applied differential gene expression transcriptome analysis to examine the changes in metabolic gene profiles between a wide range of normal tissues and cancer samples. We found that, in contrast to normal tissues which exhibit a tissue-specific expression profile, cancer samples are more homogenous despite their diverse origins. This similarity is due to a "proliferation metabolic signature" (PMS), composed of 158 genes (87 upregulated and 71 downregulated gene sets), where 143 are common to all proliferative cells but 15 are cancer specific. Intriguingly, the PMS gene set is enriched for genes encoding rate-limiting enzymes, and its upregulated set with genes associated with poor patient outcome and essential genes. Among these essential genes is ribulose-5-phosphate-3-epimerase (RPE), which encodes a pentose phosphate pathway enzyme and whose role in cancer is still unclear. Collectively, we identified a set of metabolic genes that can serve as novel cancer biomarkers and potential targets for drug development.

2-APB and CBD-Mediated Targeting of Charged Cytotoxic Compounds Into Tumor Cells Suggests the Involvement of TRPV2 Channels.

Targeted delivery of therapeutic compounds to particular cell types such that they only affect the target cells is of great clinical importance since it can minimize undesired side effects. For example, typical chemotherapeutic treatments used in the treatment of neoplastic disorders are cytotoxic not only to cancer cells but also to most normal cells when exposed to a critical concentration of the compound. As such, many chemotherapeutics exhibit severe side effects, often prohibiting their effective use in the treatment of cancer. Here, we describe a new means for facilitated delivery of a clinically used chemotherapy compound' doxorubicin, into hepatocellular carcinoma cell line (BNL1 ME). We demonstrate that these cells express a large pore, cation non-selective transient receptor potential (TRP) channel V2. We utilized this channel to shuttle doxorubicin into BNL1 ME cells. We show that co-application of either cannabidiol (CBD) or 2-APB, the activators of TRPV2 channels, together with doxorubicin leads to significantly higher accumulation of doxorubicin in BNL1 ME cells than in BNL1 ME cells that were exposed to doxorubicin alone. Moreover, we demonstrate that sub-effective doses of doxorubicin when co-applied with either 2-APB or CBD lead to a significant decrease in the number of living BNL1 ME cell and BNL1 ME cell colonies in comparison to application of doxorubicin alone. Finally, we demonstrate that the doxorubicin-mediated cell death is significantly more potent, requiring an order of magnitude lower dose, when co-applied with CBD than with 2-APB. We suggest that CBD may have a dual effect in promoting doxorubicin-mediated cell death by facilitating the entry of doxorubicin via TRPV2 channels and preventing its clearance from the cells by inhibiting P-glycoprotein ATPase transporter. Collectively, these results provide a foundation for the use of large pore cation-non selective channels as "natural" drug delivery systems for targeting specific cell types.

Microenvironmental Cues Determine Tumor Cell Susceptibility to Neutrophil Cytotoxicity.

We have recently shown that neutrophil antitumor cytotoxicity is Ca2+ dependent and is mediated by TRPM2, an H2O2-dependent Ca2+ channel. However, neutrophil antitumor activity is dependent on context and is manifested in the premetastatic niche, but not at the primary site. We therefore hypothesized that expression of TRPM2 and the consequent susceptibility to neutrophil cytotoxicity may be associated with the epithelial/mesenchymal cellular state. We found that TRPM2 expression was upregulated during epithelial-to-mesenchymal transition (EMT), and mesenchymal cells were more susceptible to neutrophil cytotoxicity. Conversely, cells undergoing mesenchymal-to-epithelial transition (MET) expressed reduced levels of TRPM2, rendering them resistant to neutrophil cytotoxicity. Cells expressing reduced levels of TRPM2 were protected from neutrophil cytotoxicity and seeded more efficiently in the premetastatic lung. These data identify TRPM2 as the link between environmental cues at the primary tumor site, tumor cell susceptibility to neutrophil cytotoxicity, and disease progression. Furthermore, these data identify EMT as a process enhancing tumor-cell immune susceptibility and, by contrast, MET as a novel mode of immune evasion.Significance: EMT is required for metastatic spread and concomitantly enhances tumor cell susceptibility to neutrophil cytotoxicity. Cancer Res; 78(17); 5050-9. ©2018 AACR.

TRPM2 Mediates Neutrophil Killing of Disseminated Tumor Cells.

Neutrophils play a critical role in cancer, with both protumor and antitumor neutrophil subpopulations reported. The antitumor neutrophil subpopulation has the capacity to kill tumor cells and limit metastatic spread, yet not all tumor cells are equally susceptible to neutrophil cytotoxicity. Because cells that evade neutrophils have greater chances of forming metastases, we explored the mechanism neutrophils use to kill tumor cells. Neutrophil cytotoxicity was previously shown to be mediated by secretion of H2O2 We report here that neutrophil cytotoxicity is Ca2+ dependent and is mediated by TRPM2, a ubiquitously expressed H2O2-dependent Ca2+ channel. Perturbing TRPM2 expression limited tumor cell proliferation, leading to attenuated tumor growth. Concomitantly, cells expressing reduced levels of TRPM2 were protected from neutrophil cytotoxicity and seeded more efficiently in the premetastatic lung.Significance: These findings identify the mechanism utilized by neutrophils to kill disseminated tumor cells and to limit metastatic spread. Cancer Res; 78(10); 2680-90. ©2018 AACR.

A PHGDH inhibitor reveals coordination of serine synthesis and one-carbon unit fate.

Serine is both a proteinogenic amino acid and the source of one-carbon units essential for de novo purine and deoxythymidine synthesis. In the canonical pathway of glucose-derived serine synthesis, Homo sapiens phosphoglycerate dehydrogenase (PHGDH) catalyzes the first, rate-limiting step. Genetic loss of PHGDH is toxic toward PHGDH-overexpressing breast cancer cell lines even in the presence of exogenous serine. Here, we used a quantitative high-throughput screen to identify small-molecule PHGDH inhibitors. These compounds reduce the production of glucose-derived serine in cells and suppress the growth of PHGDH-dependent cancer cells in culture and in orthotopic xenograft tumors. Surprisingly, PHGDH inhibition reduced the incorporation into nucleotides of one-carbon units from glucose-derived and exogenous serine. We conclude that glycolytic serine synthesis coordinates the use of one-carbon units from endogenous and exogenous serine in nucleotide synthesis, and we suggest that one-carbon unit wasting thus may contribute to the efficacy of PHGDH inhibitors in vitro and in vivo.

MERAV: a tool for comparing gene expression across human tissues and cell types.

The oncogenic transformation of normal cells into malignant, rapidly proliferating cells requires major alterations in cell physiology. For example, the transformed cells remodel their metabolic processes to supply the additional demand for cellular building blocks. We have recently demonstrated essential metabolic processes in tumor progression through the development of a methodological analysis of gene expression. Here, we present the Metabolic gEne RApid Visualizer (MERAV, http://merav.wi.mit.edu), a web-based tool that can query a database comprising ∼4300 microarrays, representing human gene expression in normal tissues, cancer cell lines and primary tumors. MERAV has been designed as a powerful tool for whole genome analysis which offers multiple advantages: one can search many genes in parallel; compare gene expression among different tissue types as well as between normal and cancer cells; download raw data; and generate heatmaps; and finally, use its internal statistical tool. Most importantly, MERAV has been designed as a unique tool for analyzing metabolic processes as it includes matrixes specifically focused on metabolic genes and is linked to the Kyoto Encyclopedia of Genes and Genomes pathway search.

A respiratory chain controlled signal transduction cascade in the mitochondrial intermembrane space mediates hydrogen peroxide signaling.

Reactive oxygen species (ROS) such as hydrogen peroxide (H2O2) govern cellular homeostasis by inducing signaling. H2O2 modulates the activity of phosphatases and many other signaling molecules through oxidation of critical cysteine residues, which led to the notion that initiation of ROS signaling is broad and nonspecific, and thus fundamentally distinct from other signaling pathways. Here, we report that H2O2 signaling bears hallmarks of a regular signal transduction cascade. It is controlled by hierarchical signaling events resulting in a focused response as the results place the mitochondrial respiratory chain upstream of tyrosine-protein kinase Lyn, Lyn upstream of tyrosine-protein kinase SYK (Syk), and Syk upstream of numerous targets involved in signaling, transcription, translation, metabolism, and cell cycle regulation. The active mediators of H2O2 signaling colocalize as H2O2 induces mitochondria-associated Lyn and Syk phosphorylation, and a pool of Lyn and Syk reside in the mitochondrial intermembrane space. Finally, the same intermediaries control the signaling response in tissues and species responsive to H2O2 as the respiratory chain, Lyn, and Syk were similarly required for H2O2 signaling in mouse B cells, fibroblasts, and chicken DT40 B cells. Consistent with a broad role, the Syk pathway is coexpressed across tissues, is of early metazoan origin, and displays evidence of evolutionary constraint in the human. These results suggest that H2O2 signaling is under control of a signal transduction pathway that links the respiratory chain to the mitochondrial intermembrane space-localized, ubiquitous, and ancient Syk pathway in hematopoietic and nonhematopoietic cells.

Dihydropyrimidine accumulation is required for the epithelial-mesenchymal transition.

It is increasingly appreciated that oncogenic transformation alters cellular metabolism to facilitate cell proliferation, but less is known about the metabolic changes that promote cancer cell aggressiveness. Here, we analyzed metabolic gene expression in cancer cell lines and found that a set of high-grade carcinoma lines expressing mesenchymal markers share a unique 44 gene signature, designated the "mesenchymal metabolic signature" (MMS). A FACS-based shRNA screen identified several MMS genes as essential for the epithelial-mesenchymal transition (EMT), but not for cell proliferation. Dihydropyrimidine dehydrogenase (DPYD), a pyrimidine-degrading enzyme, was highly expressed upon EMT induction and was necessary for cells to acquire mesenchymal characteristics in vitro and for tumorigenic cells to extravasate into the mouse lung. This role of DPYD was mediated through its catalytic activity and enzymatic products, the dihydropyrimidines. Thus, we identify metabolic processes essential for the EMT, a program associated with the acquisition of metastatic and aggressive cancer cell traits.

Critical role for IL-1ß in DNA damage-induced mucositis.

ß-TrCP, the substrate recognition subunit of SCF-type ubiquitin ligases, is ubiquitously expressed from two distinct paralogs, targeting for degradation many regulatory proteins, among which is the NF-κB inhibitor IκB. To appreciate tissue-specific roles of ß-TrCP, we studied the consequences of inducible ablation of three or all four alleles of the E3 in the mouse gut. The ablation resulted in mucositis, a destructive gut mucosal inflammation, which is a common complication of different cancer therapies and represents a major obstacle to successful chemoradiation therapy. We identified epithelial-derived IL-1ß as the culprit of mucositis onset, inducing mucosal barrier breach. Surprisingly, epithelial IL-1ß is induced by DNA damage via an NF-κB-independent mechanism. Tissue damage caused by gut barrier disruption is exacerbated in the absence of NF-κB, with failure to express the endogenous IL-1ß receptor antagonist IL-1Ra upon four-allele loss. Antibody neutralization of IL-1ß prevents epithelial tight junction dysfunction and alleviates mucositis in ß-TrCP-deficient mice. IL-1ß antagonists should thus be considered for prevention and treatment of severe morbidity associated with mucositis.

Functional genomics reveal that the serine synthesis pathway is essential in breast cancer.

Cancer cells adapt their metabolic processes to drive macromolecular biosynthesis for rapid cell growth and proliferation. RNA interference (RNAi)-based loss-of-function screening has proven powerful for the identification of new and interesting cancer targets, and recent studies have used this technology in vivo to identify novel tumour suppressor genes. Here we developed a method for identifying novel cancer targets via negative-selection RNAi screening using a human breast cancer xenograft model at an orthotopic site in the mouse. Using this method, we screened a set of metabolic genes associated with aggressive breast cancer and stemness to identify those required for in vivo tumorigenesis. Among the genes identified, phosphoglycerate dehydrogenase (PHGDH) is in a genomic region of recurrent copy number gain in breast cancer and PHGDH protein levels are elevated in 70% of oestrogen receptor (ER)-negative breast cancers. PHGDH catalyses the first step in the serine biosynthesis pathway, and breast cancer cells with high PHGDH expression have increased serine synthesis flux. Suppression of PHGDH in cell lines with elevated PHGDH expression, but not in those without, causes a strong decrease in cell proliferation and a reduction in serine synthesis. We find that PHGDH suppression does not affect intracellular serine levels, but causes a drop in the levels of α-ketoglutarate, another output of the pathway and a tricarboxylic acid (TCA) cycle intermediate. In cells with high PHGDH expression, the serine synthesis pathway contributes approximately 50% of the total anaplerotic flux of glutamine into the TCA cycle. These results reveal that certain breast cancers are dependent upon increased serine pathway flux caused by PHGDH overexpression and demonstrate the utility of in vivo negative-selection RNAi screens for finding potential anticancer targets.

Genome-scale RNAi on living-cell microarrays identifies novel regulators of Drosophila melanogaster TORC1-S6K pathway signaling.

The evolutionarily conserved target of rapamycin complex 1 (TORC1) controls cell growth in response to nutrient availability and growth factors. TORC1 signaling is hyperactive in cancer, and regulators of TORC1 signaling may represent therapeutic targets for human diseases. To identify novel regulators of TORC1 signaling, we performed a genome-scale RNA interference screen on microarrays of Drosophila melanogaster cells expressing human RPS6, a TORC1 effector whose phosphorylated form we detected by immunofluorescence. Our screen revealed that the TORC1-S6K-RPS6 signaling axis is regulated by many subcellular components, including the Class I vesicle coat (COPI), the spliceosome, the proteasome, the nuclear pore, and the translation initiation machinery. Using additional RNAi reagents, we confirmed 70 novel genes as significant on-target regulators of RPS6 phosphorylation, and we characterized them with extensive secondary assays probing various arms of the TORC1 pathways, identifying functional relationships among those genes. We conclude that cell-based microarrays are a useful platform for genome-scale and secondary screening in Drosophila, revealing regulators that may represent drug targets for cancers and other diseases of deregulated TORC1 signaling.

Specific phosphorylation and activation of ERK1c by MEK1b: a unique route in the ERK cascade.

Extracellular signal-regulated kinases (ERKs) are key signaling molecules that regulate a large number of cellular processes, including mitosis. We showed previously that ERK1c, an alternatively spliced form of ERK1, facilitates mitotic Golgi fragmentation without the involvement of ERK1 and ERK2. Here we demonstrate that activation of ERK1c is mainly mediated by mitogen-activated protein kinase (MAPK)/ERK kinase 1b (MEK1b), which is an alternatively spliced form of MEK1 that was previously considered an inactive kinase. MEK1b phosphorylation and activity are preferentially stimulated by nocodazole, to induce its specific activity toward ERK1c. MEK1/2, on the other hand, preferentially target ERK1/2 in response to growth factors, such as EGF. As previously demonstrated for ERK1c, also MEK1b expression and activity are elevated during mitosis, and thereby enhance Golgi fragmentation and mitotic rate. MEK1 activity is also increased during mitosis, but this isoform facilitates mitotic progression without affecting the Golgi architecture. These results illustrate that the ERK cascade is divided into two routes: the classic MEK1/2-ERK1/2 and the splice-variant MEK1b-ERK1c, each of which regulates distinct cellular processes and thus extends the cascade specificity.