Multiomic Metabolic Enrichment Network Analysis Reveals Metabolite–Protein Physical Interaction Subnetworks Altered in Cancer

Metabolism is recognized as an important driver of cancer progression and other complex diseases, but global metabolite profiling remains a challenge. Protein expression profiling is often a poor proxy since existing pathway enrichment models provide an incomplete mapping between the proteome and metabolism. To overcome these gaps, we introduce multiomic metabolic enrichment network analysis (MOMENTA), an integrative multiomic data analysis framework for more accurately deducing metabolic pathway changes from proteomics data alone in a gene set analysis context by leveraging protein interaction networks to extend annotated metabolic models. We apply MOMENTA to proteomic data from diverse cancer cell lines and human tumors to demonstrate its utility at revealing variation in metabolic pathway activity across cancer types, which we verify using independent metabolomics measurements. The novel metabolic networks we uncover in breast cancer and other tumors are linked to clinical outcomes, underscoring the pathophysiological relevance of the findings.     

Four-carbon dicarboxylic acid production through the reductive branch of the open cyanobacterial tricarboxylic acid cycle in Synechocystis sp. PCC 6803

Succinate, fumarate, and malate are valuable four-carbon (C4) dicarboxylic acids used for producing plastics and food additives. C4 dicarboxylic acid is biologically produced by heterotrophic organisms. However, current biological production requires organic carbon sources that compete with food uses. Herein, we report C4 dicarboxylic acid production from CO₂ using metabolically engineered Synechocystis sp. PCC 6803. Overexpression of citH, encoding malate dehydrogenase (MDH), resulted in the enhanced production of succinate, fumarate, and malate. citH overexpression increased the reductive branch of the open cyanobacterial tricarboxylic acid (TCA) cycle flux. Furthermore, product stripping by medium exchanges increased the C4 dicarboxylic acid levels; product inhibition and acidification of the media were the limiting factors for succinate production. Our results demonstrate that MDH is a key regulator that activates the reductive branch of the open cyanobacterial TCA cycle. The study findings suggest that cyanobacteria can act as a biocatalyst for converting CO₂ to carboxylic acids.     

Sorting through the extensive and confusing roles of sortilin in metabolic disease

Sortilin is a post-Golgi trafficking receptor homologous to the yeast vacuolar protein sorting receptor 10 (VPS10). The VPS10 motif on sortilin is a 10-bladed β-propeller structure capable of binding more than 50 proteins, covering a wide range of biological functions including lipid and lipoprotein metabolism, neuronal growth and death, inflammation, and lysosomal degradation. Sortilin has a complex cellular trafficking itinerary, where it functions as a receptor in the trans-Golgi network, endosomes, secretory vesicles, multivesicular bodies, and at the cell surface. In addition, sortilin is associated with hypercholesterolemia, Alzheimer’s disease, prion diseases, Parkinson’s disease, and inflammation syndromes. The 1p13.3 locus containing SORT1, the gene encoding sortilin, carries the strongest association with LDL-C of all loci in human genome-wide association studies. However, the mechanism by which sortilin influences LDL-C is unclear. Here, we review the role sortilin plays in cardiovascular and metabolic diseases and describe in detail the large and often contradictory literature on the role of sortilin in the regulation of LDL-C levels.     

Enhanced targeted DNA methylation of the CMV and endogenous promoters with dCas9-DNMT3A3L entails distinct subsequent histone modification changes in CHO cells

With the emergence of new CRISPR/dCas9 tools that enable site specific modulation of DNA methylation and histone modifications, more detailed investigations of the contribution of epigenetic regulation to the precise phenotype of cells in culture, including recombinant production subclones, is now possible. These also allow a wide range of applications in metabolic engineering once the impact of such epigenetic modifications on the chromatin state is available.In this study, enhanced DNA methylation tools were targeted to a recombinant viral promoter (CMV), an endogenous promoter that is silenced in its native state in CHO cells, but had been reactivated previously (β-galactoside α-2,6-sialyltransferase 1) and an active endogenous promoter (α-1,6-fucosyltransferase), respectively. Comparative ChIP-analysis of histone modifications revealed a general loss of active promoter histone marks and the acquisition of distinct repressive heterochromatin marks after targeted methylation. On the other hand, targeted demethylation resulted in autologous acquisition of active promoter histone marks and loss of repressive heterochromatin marks. These data suggest that DNA methylation directs the removal or deposition of specific histone marks associated with either active, poised or silenced chromatin. Moreover, we show that de novo methylation of the CMV promoter results in reduced transgene expression in CHO cells. Although targeted DNA methylation is not efficient, the transgene is repressed, thus offering an explanation for seemingly conflicting reports about the source of CMV promoter instability in CHO cells.Importantly, modulation of epigenetic marks enables to nudge the cell into a specific gene expression pattern or phenotype, which is stabilized in the cell by autologous addition of further epigenetic marks. Such engineering strategies have the added advantage of being reversible and potentially tunable to not only turn on or off a targeted gene, but also to achieve the setting of a desirable expression level.     

Predictive macroscopic modeling of cell growth,metabolism and monoclonal antibody production: Case study of a CHO fed-batch production

We describe a systematic approach to establish predictive models of CHO cell growth, cell metabolism and monoclonal antibody (mAb) formation during biopharmaceutical production. The prediction is based on a combination of an empirical metabolic model connecting extracellular metabolic fluxes with cellular growth and product formation with mixed Monod-inhibition type kinetics that we generalized to every possible external metabolite. We describe the maximum specific growth rate as a function of the integral viable cell density (IVCD). Moreover, we also take into account the accumulation of metabolites in intracellular pools that can influence cell growth. This is possible even without identification and quantification of these metabolites as illustrated with fed-batch cultures of Chinese Hamster Ovary (CHO) cells producing a mAb. The impact of cysteine and tryptophan on cell growth and cell productivity was assessed, and the resulting macroscopic model was successfully used to predict the impact of new, untested feeding strategies on cell growth and mAb production. This model combining piecewise linear relationships between metabolic rates, growth rate and production rate together with Monod-inhibition type models for cell growth did well in predicting cell culture performance in fed-batch cultures even outside the range of experimental data used for establishing the model. It could therefore also successfully be applied for in silico prediction of optimal operating conditions.     

ABHD4 regulates adipocyte differentiation in vitro but does not affect adipose tissue lipid metabolism in mice

Alpha/beta hydrolase domain-containing protein 4 (ABHD4) catalyzes the deacylation of N-acyl phosphatidyl-ethanolamine (NAPE) and lyso-NAPE to produce glycerophospho-N-acyl ethanolamine (GP-NAE). Through a variety of metabolic enzymes, NAPE, lyso-NAPE, and GP-NAE are ultimately converted into NAE, a group of bioactive lipids that control many physiological processes including inflammation, cognition, food intake, and lipolysis (i.e., oleoylethanolamide or OEA). In a diet-induced obese mouse model, adipose tissue Abhd4 gene expression positively correlated with adiposity. However, it is unknown whether Abhd4 is a causal or a reactive gene to obesity. To fill this knowledge gap, we generated an Abhd4 knockout (KO) 3T3-L1 pre-adipocyte. During adipogenic stimulation, Abhd4 KO pre-adipocytes had increased adipogenesis and lipid accumulation, suggesting Abhd4 is responding to (a reactive gene), not contributing to (not a causal gene), adiposity, and may serve as a mechanism for protecting against obesity. However, we did not observe any differences in adiposity and metabolic outcomes between whole-body Abhd4 KO or adipocyte-specific Abhd4 KO mice and their littermate control mice (both male and female) on chow or a high-fat diet. This might be because we found that deletion of Abhd4 did not affect NAE such as OEA production, even though Abhd4 was highly expressed in adipose tissue and correlated with fasting adipose OEA levels and lipolysis. These data suggest that ABHD4 regulates adipocyte differentiation in vitro but does not affect adipose tissue lipid metabolism in mice despite nutrient overload, possibly due to compensation from other NAPE and NAE metabolic enzymes.     

Ionizing Radiation Drives Key Regulators of Antigen Presentation and a Global Expansion of the Immunopeptidome

Little is known about the pathways regulating MHC antigen presentation and the identity of treatment-specific T cell antigens induced by ionizing radiation. For this reason, we investigated the radiation-specific changes in the colorectal tumor cell proteome. We found an increase in DDX58 and ZBP1 protein expression, two nucleic acid sensing molecules likely involved in induction of the dominant interferon response signature observed after genotoxic insult. We further observed treatment-induced changes in key regulators and effector proteins of the antigen processing and presentation machinery. Differential regulation of MHC allele expression was further driving the presentation of a significantly broader MHC-associated peptidome postirradiation, defining a radiation-specific peptide repertoire. Interestingly, treatment-induced peptides originated predominantly from proteins involved in catecholamine synthesis and metabolic pathways. A nuanced relationship between protein expression and antigen presentation was observed where radiation-induced changes in proteins do not correlate with increased presentation of associated peptides. Finally, we detected an increase in the presentation of a tumor-specific neoantigen derived from Mtch1. This study provides new insights into how radiation enhances antigen processing and presentation that could be suitable for the development of combinatorial therapies. Data are available via ProteomeXchange with identifier PXD032003.     

Deconvoluting the biology and druggability of protein lipidation using chemical proteomics

Lipids are indispensable cellular building blocks, and their post-translational attachment to proteins makes them important regulators of many biological processes. Dysfunction of protein lipidation is also implicated in many pathological states, yet its systematic analysis presents significant challenges. Thanks to innovations in chemical proteomics, lipidation can now be readily studied by metabolic tagging using functionalized lipid analogs, enabling global profiling of lipidated substrates using mass spectrometry. This has spearheaded the first deconvolution of their full scope in a range of contexts, from cells to pathogens and multicellular organisms. Protein N-myristoylation, S-acylation, and S-prenylation are the most well-studied lipid post-translational modifications because of their extensive contribution to the regulation of diverse cellular processes. In this review, we focus on recent advances in the study of these post-translational modifications, with an emphasis on how novel mass spectrometry methods have elucidated their roles in fundamental biological processes.     

Identification of two lipid phosphatases that regulate sphingosine-1-phosphate cellular uptake and recycling

Sphingosine-1-phosphate (S1P) is a sphingolipid metabolite that serves as a potent extracellular signaling molecule. Metabolic regulation of extracellular S1P levels impacts key cellular activities through altered S1P receptor signaling. Although the pathway through which S1P is degraded within the cell and thereby eliminated from reuse has been previously described, the mechanism used for S1P cellular uptake and the subsequent recycling of its sphingoid base into the sphingolipid synthesis pathway is not completely understood. To identify the genes within this S1P uptake and recycling pathway, we performed a genome-wide CRISPR/Cas9 KO screen using a positive-selection scheme with Shiga toxin, which binds a cell-surface glycosphingolipid receptor, globotriaosylceramide (Gb3), and causes lethality upon internalization. The screen was performed in HeLa cells with their sphingolipid de novo pathway disabled so that Gb3 cell-surface expression was dependent on salvage of the sphingoid base of S1P taken up from the medium. The screen identified a suite of genes necessary for S1P uptake and the recycling of its sphingoid base to synthesize Gb3, including two lipid phosphatases, PLPP3 (phospholipid phosphatase 3) and SGPP1 (S1P phosphatase 1). The results delineate a pathway in which plasma membrane–bound PLPP3 dephosphorylates extracellular S1P to sphingosine, which then enters cells and is rephosphorylated to S1P by the sphingosine kinases. This rephosphorylation step is important to regenerate intracellular S1P as a branch-point substrate that can be routed either for dephosphorylation to salvage sphingosine for recycling into complex sphingolipid synthesis or for degradation to remove it from the sphingolipid synthesis pathway.     

Transgelin: a new gene involved in LDL endocytosis identified by a genome-wide CRISPR-Cas9 screen

A significant proportion of patients with elevated LDL and a clinical presentation of familial hypercholesterolemia do not carry known genetic mutations associated with hypercholesterolemia, such as defects in the LDL receptor. To identify new genes involved in the cellular uptake of LDL, we developed a novel whole-genome clustered regularly interspaced short palindromic repeat-Cas9 KO screen in HepG2 cells. We identified transgelin (TAGLN), an actin-binding protein, as a potentially new gene involved in LDL endocytosis. In silico validation demonstrated that genetically predicted differences in expression of TAGLN in human populations were significantly associated with elevated plasma lipids (triglycerides, total cholesterol, and LDL-C) in the Global Lipids Genetics Consortium and lipid-related phenotypes in the UK Biobank. In biochemical studies, TAGLN-KO HepG2 cells showed a reduction in cellular LDL uptake, as measured by flow cytometry. In confocal microscopy imaging, TAGLN-KO cells had disrupted actin filaments as well as an accumulation of LDL receptor on their surface because of decreased receptor internalization. Furthermore, TAGLN-KO cells exhibited a reduction in total and free cholesterol content, activation of SREBP2, and a compensatory increase in cholesterol biosynthesis. TAGLN deficiency also disrupted the uptake of VLDL and transferrin, other known cargoes for receptors that depend upon clathrin-mediated endocytosis. Our data suggest that TAGLN is a novel factor involved in the actin-dependent phase of clathrin-mediated endocytosis of LDL. The identification of novel genes involved in the endocytic uptake of LDL may improve the diagnosis of hypercholesterolemia and provide future therapeutic targets for the prevention of cardiovascular disease.     

New Global Insights on the Regulation of the Biphasic Life Cycle and Virulence Via ClpP-Dependent Proteolysis in Legionella pneumophila

Legionella pneumophila, an environmental bacterium that parasitizes protozoa, causes Legionnaires’ disease in humans that is characterized by severe pneumonia. This bacterium adopts a distinct biphasic life cycle consisting of a nonvirulent replicative phase and a virulent transmissive phase in response to different environmental conditions. Hence, the timely and fine-tuned expression of growth and virulence factors in a life cycle–dependent manner is crucial for survival and replication. Here, we report that the completion of the biphasic life cycle and bacterial pathogenesis is greatly dependent on the protein homeostasis regulated by caseinolytic protease P (ClpP)-dependent proteolysis. We characterized the ClpP-dependent dynamic profiles of the regulatory and substrate proteins during the biphasic life cycle of L. pneumophila using proteomic approaches and discovered that ClpP-dependent proteolysis specifically and conditionally degraded the substrate proteins, thereby directly playing a regulatory role or indirectly controlling cellular events via the regulatory proteins. We further observed that ClpP-dependent proteolysis is required to monitor the abundance of fatty acid biosynthesis–related protein Lpg0102/Lpg0361/Lpg0362 and SpoT for the normal regulation of L. pneumophila differentiation. We also found that the control of the biphasic life cycle and bacterial virulence is independent. Furthermore, the ClpP-dependent proteolysis of Dot/Icm (defect in organelle trafficking/intracellular multiplication) type IVB secretion system and effector proteins at a specific phase of the life cycle is essential for bacterial pathogenesis. Therefore, our findings provide novel insights on ClpP-dependent proteolysis, which spans a broad physiological spectrum involving key metabolic pathways that regulate the transition of the biphasic life cycle and bacterial virulence of L. pneumophila, facilitating adaptation to aquatic and intracellular niches.     

Interactome Analysis of Human Phospholipase D and Phosphatidic Acid-Associated Protein Network

Mammalian phospholipase D (PLD) enzyme family consists of six members. Among them, PLD1/2/6 catalyzes phosphatidic acid (PA) production, while PLD3/4/5 has no catalytic activities. Deregulation of the PLD-PA lipid signaling has been associated with various human diseases including cancer. However, a comprehensive analysis of the regulators and effectors for this crucial lipid metabolic pathway has not been fully achieved. Using a proteomic approach, we defined the protein interaction network for the human PLD family of enzymes and PA and revealed diverse cellular signaling events involving them. Through it, we identified PJA2 as a novel E3 ubiquitin ligase for PLD1 involved in control of the PLD1-mediated mammalian target of rapamycin signaling. Additionally, we showed that PA interacted with and positively regulated sphingosine kinase 1. Taken together, our study not only generates a rich interactome resource for further characterizing the human PLD-PA lipid signaling but also connects this important metabolic pathway with numerous biological processes.     

In-planta production of the biodegradable polyester precursor 2-pyrone-4,6-dicarboxylic acid (PDC): Stacking reduced biomass recalcitrance with value-added co-product

2-Pyrone-4,6-dicarboxylic acid (PDC), a chemically stable intermediate that naturally occurs during microbial degradation of lignin by bacteria, represents a promising building block for diverse biomaterials and polyesters such as biodegradable plastics. The lack of a chemical synthesis method has hindered large-scale utilization of PDC and metabolic engineering approaches for its biosynthesis have recently emerged. In this study, we demonstrate a strategy for the production of PDC via manipulation of the shikimate pathway using plants as green factories. In tobacco leaves, we first showed that transient expression of bacterial feedback-resistant 3-deoxy-D-arabinoheptulosonate 7-phosphate synthase (AroG) and 3-dehydroshikimate dehydratase (QsuB) produced high titers of protocatechuate (PCA), which was in turn efficiently converted into PDC upon co-expression of PCA 4,5-dioxygenase (PmdAB) and 4-carboxy-2-hydroxymuconate-6-semialdehyde dehydrogenase (PmdC) derived from Comamonas testosteroni. We validated that stable expression of AroG in Arabidopsis in a genetic background containing the QsuB gene enhanced PCA content in plant biomass, presumably via an increase of the carbon flux through the shikimate pathway. Further, introducing AroG and the PDC biosynthetic genes (PmdA, PmdB, and PmdC) into the Arabidopsis QsuB background, or introducing the five genes (AroG, QsuB, PmdA, PmdB, and PmdC) stacked on a single construct into wild-type plants, resulted in PDC titers of ~1% and ~3% dry weight in plant biomass, respectively. Consistent with previous studies of plants expressing QsuB, all PDC producing lines showed strong reduction in lignin content in stems. This low lignin trait was accompanied with improvements of biomass saccharification efficiency due to reduced cell wall recalcitrance to enzymatic degradation. Importantly, most transgenic lines showed no reduction in biomass yields. Therefore, we conclude that engineering plants with the proposed de-novo PDC pathway provides an avenue to enrich biomass with a value-added co-product while simultaneously improving biomass quality for the supply of fermentable sugars. Implementing this strategy into bioenergy crops has the potential to support existing microbial fermentation approaches that exploit lignocellulosic biomass feedstocks for PDC production.     

Low Cell Number Proteomic Analysis Using In-Cell Protease Digests Reveals a Robust Signature for Cell Cycle State Classification

Comprehensive proteome analysis of rare cell phenotypes remains a significant challenge. We report a method for low cell number MS-based proteomics using protease digestion of mildly formaldehyde-fixed cells in cellulo, which we call the “in-cell digest.” We combined this with averaged MS1 precursor library matching to quantitatively characterize proteomes from low cell numbers of human lymphoblasts. About 4500 proteins were detected from 2000 cells, and 2500 proteins were quantitated from 200 lymphoblasts. The ease of sample processing and high sensitivity makes this method exceptionally suited for the proteomic analysis of rare cell states, including immune cell subsets and cell cycle subphases. To demonstrate the method, we characterized the proteome changes across 16 cell cycle states (CCSs) isolated from an asynchronous TK6 cells, avoiding synchronization. States included late mitotic cells present at extremely low frequency. We identified 119 pseudoperiodic proteins that vary across the cell cycle. Clustering of the pseudoperiodic proteins showed abundance patterns consistent with “waves” of protein degradation in late S, at the G2&M border, midmitosis, and at mitotic exit. These clusters were distinguished by significant differences in predicted nuclear localization and interaction with the anaphase-promoting complex/cyclosome. The dataset also identifies putative anaphase-promoting complex/cyclosome substrates in mitosis and the temporal order in which they are targeted for degradation. We demonstrate that a protein signature made of these 119 high-confidence cell cycle–regulated proteins can be used to perform unbiased classification of proteomes into CCSs. We applied this signature to 296 proteomes that encompass a range of quantitation methods, cell types, and experimental conditions. The analysis confidently assigns a CCS for 49 proteomes, including correct classification for proteomes from synchronized cells. We anticipate that this robust cell cycle protein signature will be crucial for classifying cell states in single-cell proteomes.     

Quantification and Proteomic Characterization of β-Hydroxybutyrylation Modification in the Hearts of AMPKα2 Knockout Mice

AMP-activated protein kinase alpha 2 (AMPKα2) regulates energy metabolism, protein synthesis, and glucolipid metabolism myocardial cells. Ketone bodies produced by fatty acid β-oxidation, especially β-hydroxybutyrate, are fatty energy–supplying substances for the heart, brain, and other organs during fasting and long-term exercise. They also regulate metabolic signaling for multiple cellular functions. Lysine β-hydroxybutyrylation (Kbhb) is a β-hydroxybutyrate–mediated protein posttranslational modification. Histone Kbhb has been identified in yeast, mouse, and human cells. However, whether AMPK regulates protein Kbhb is yet unclear. Hence, the present study explored the changes in proteomics and Kbhb modification omics in the hearts of AMPKα2 knockout mice using a comprehensive quantitative proteomic analysis. Based on mass spectrometry (LC-MS/MS) analysis, the number of 1181 Kbhb modified sites in 455 proteins were quantified between AMPKα2 knockout mice and wildtype mice; 244 Kbhb sites in 142 proteins decreased or increased after AMPKα2 knockout (fold change >1.5 or <1/1.5, p < 0.05). The regulation of Kbhb sites in 26 key enzymes of fatty acid degradation and tricarboxylic acid cycle was noted in AMPKα2 knockout mouse cardiomyocytes. These findings, for the first time, identified proteomic features and Kbhb modification of cardiomyocytes after AMPKα2 knockout, suggesting that AMPKα2 regulates energy metabolism by modifying protein Kbhb.     

Comparative Proteomic Analysis Identifies Key Metabolic Regulators of Gemcitabine Resistance in Pancreatic Cancer

Pancreatic adenocarcinoma (PDAC) is highly refractory to treatment. Standard-of-care gemcitabine (Gem) provides only modest survival benefits, and development of Gem resistance (GemR) compromises its efficacy. Highly GemR clones of Gem-sensitive MIAPaCa-2 cells were developed to investigate the molecular mechanisms of GemR and implemented global quantitative differential proteomics analysis with a comprehensive, reproducible ion-current–based MS1 workflow to quantify ～6000 proteins in all samples. In GemR clone MIA-GR8, cellular metabolism, proliferation, migration, and ‘drug response’ mechanisms were the predominant biological processes altered, consistent with cell phenotypic alterations in cell cycle and motility. S100 calcium binding protein A4 was the most downregulated protein, as were proteins associated with glycolytic and oxidative energy production. Both responses would reduce tumor proliferation. Upregulation of mesenchymal markers was prominent, and cellular invasiveness increased. Key enzymes in Gem metabolism pathways were altered such that intracellular utilization of Gem would decrease. Ribonucleoside-diphosphate reductase large subunit was the most elevated Gem metabolizing protein, supporting its critical role in GemR. Lower Ribonucleoside-diphosphate reductase large subunit expression is associated with better clinical outcomes in PDAC, and its downregulation paralleled reduced MIAPaCa-2 proliferation and migration and increased Gem sensitivity. Temporal protein-level Gem responses of MIAPaCa-2 versus GemR cell lines (intrinsically GemR PANC-1 and acquired GemR MIA-GR8) implicate adaptive changes in cellular response systems for cell proliferation and drug transport and metabolism, which reduce cytotoxic Gem metabolites, in DNA repair, and additional responses, as key contributors to the complexity of GemR in PDAC. These findings additionally suggest targetable therapeutic vulnerabilities for GemR PDAC patients.     

Glutamine-Fructose-6-Phosphate Transaminase 2 (GFPT2) Is Upregulated in Breast Epithelial–Mesenchymal Transition and Responds to Oxidative Stress

Breast cancer cells that have undergone partial epithelial–mesenchymal transition (EMT) are believed to be more invasive than cells that have completed EMT. To study metabolic reprogramming in different mesenchymal states, we analyzed protein expression following EMT in the breast epithelial cell model D492 with single-shot LFQ supported by a SILAC proteomics approach. The D492 EMT cell model contains three cell lines: the epithelial D492 cells, the mesenchymal D492M cells, and a partial mesenchymal, tumorigenic variant of D492 that overexpresses the oncogene HER2. The analysis classified the D492 and D492M cells as basal-like and D492HER2 as claudin-low. Comparative analysis of D492 and D492M to tumorigenic D492HER2 differentiated metabolic markers of migration from those of invasion. Glutamine-fructose-6-phosphate transaminase 2 (GFPT2) was one of the top dysregulated enzymes in D492HER2. Gene expression analysis of the cancer genome atlas showed that GFPT2 expression was a characteristic of claudin-low breast cancer. siRNA-mediated knockdown of GFPT2 influenced the EMT marker vimentin and both cell growth and invasion in vitro and was accompanied by lowered metabolic flux through the hexosamine biosynthesis pathway (HBP). Knockdown of GFPT2 decreased cystathionine and sulfide:quinone oxidoreductase (SQOR) in the transsulfuration pathway that regulates H₂S production and mitochondrial homeostasis. Moreover, GFPT2 was within the regulation network of insulin and EGF, and its expression was regulated by reduced glutathione (GSH) and suppressed by the oxidative stress regulator GSK3-β. Our results demonstrate that GFPT2 controls growth and invasion in the D492 EMT model, is a marker for oxidative stress, and associated with poor prognosis in claudin-low breast cancer.