# Glyphosate-Induced Metabolic and Immune Modulation in Hepatoma Cells: Identification of Key Genes as Diagnostic and Therapeutic Targets Using an In Silico Systems Biology Approach

**Authors:** Divya Mishra

PMC · DOI: 10.3390/jox16020051 · 2026-03-19

## TL;DR

This study uses in silico methods to identify key genes affected by glyphosate in liver cancer cells, offering potential biomarkers for diagnosis and treatment.

## Contribution

The study identifies novel hub genes linked to glyphosate-induced metabolic and immune changes in hepatoma cells.

## Key findings

- Glyphosate exposure upregulates genes involved in glucose metabolism and stress responses.
- Downregulated genes are associated with lipid metabolism and immune regulation pathways.
- Hub genes show strong diagnostic potential for glyphosate exposure in liver cancer cells.

## Abstract

Glyphosate, one of the most widely used herbicides worldwide, has raised significant concerns regarding its potential involvement in hepatotoxicity and molecular changes associated with liver cancer biology. These concerns highlight the need to better understand its underlying molecular mechanisms in hepatoma cells. Emerging evidence suggests that glyphosate exposure may increase the risk of liver cancer and chronic liver disease. However, the precise molecular alterations and promising biomarkers associated with glyphosate-induced hepatic toxicity and disease remain largely unexplored. In this study, an RNA-Seq-based in silico systems biology approach was employed to elucidate glyphosate-induced differential transcriptional profiling in hepatoma cells. This analysis revealed significant transcriptional profiling characterized by the upregulated hub genes ATF3, JUNB, ALDOA, FOSB, PFKFB3, G6PD, ENO2, HK2, FOS and PGK1. These genes were primarily associated with glucose metabolism, TNF-α/NF-κB signaling, epithelial–mesenchymal transition (EMT) and cellular stress responses. Conversely, several key genes were significantly downregulated, including PIK3R1, FYN, CEBPA, MLXIPL, PPARA, CD36, PCK2, PNPLA3, NR1H4 and MGLL, which were involved in lipid metabolism, immune regulation and non-alcoholic fatty liver disease (NAFLD) pathways. Notably, all hub genes demonstrated strong diagnostic performance, highlighting their potential as sensitive biomarkers of glyphosate exposure. Collectively, this study provides comprehensive insights into gene expression changes associated with glyphosate exposure in hepatoma cells, linking them to hepatic metabolic dysregulation and immune modulation and suggesting a panel of hub genes with potential diagnostic and therapeutic significance.

## Linked entities

- **Genes:** ATF3 (activating transcription factor 3) [NCBI Gene 467], JUNB (JunB proto-oncogene, AP-1 transcription factor subunit) [NCBI Gene 3726], ALDOA (aldolase, fructose-bisphosphate A) [NCBI Gene 226], FOSB (FosB proto-oncogene, AP-1 transcription factor subunit) [NCBI Gene 2354], PFKFB3 (6-phosphofructo-2-kinase/fructose-2,6-biphosphatase 3) [NCBI Gene 5209], G6PD (glucose-6-phosphate dehydrogenase) [NCBI Gene 2539], ENO2 (enolase 2) [NCBI Gene 2026], HK2 (hexokinase 2) [NCBI Gene 3099], FOS (Fos proto-oncogene, AP-1 transcription factor subunit) [NCBI Gene 2353], PGK1 (phosphoglycerate kinase 1) [NCBI Gene 5230], PIK3R1 (phosphoinositide-3-kinase regulatory subunit 1) [NCBI Gene 5295], FYN (FYN proto-oncogene, Src family tyrosine kinase) [NCBI Gene 2534], CEBPA (CCAAT enhancer binding protein alpha) [NCBI Gene 1050], MLXIPL (MLX interacting protein like) [NCBI Gene 51085], PPARA (peroxisome proliferator activated receptor alpha) [NCBI Gene 5465], CD36 (CD36 molecule (CD36 blood group)) [NCBI Gene 948], PCK2 (phosphoenolpyruvate carboxykinase 2, mitochondrial) [NCBI Gene 5106], PNPLA3 (patatin like domain 3, 1-acylglycerol-3-phosphate O-acyltransferase) [NCBI Gene 80339], NR1H4 (nuclear receptor subfamily 1 group H member 4) [NCBI Gene 9971], MGLL (monoglyceride lipase) [NCBI Gene 11343]
- **Chemicals:** glyphosate (PubChem CID 3496)
- **Diseases:** liver cancer (MONDO:0002691), non-alcoholic fatty liver disease (MONDO:0013209)

## Full-text entities

- **Genes:** Pfkfb3 (6-phosphofructo-2-kinase/fructose-2,6-biphosphatase 3) [NCBI Gene 170768] {aka E330010H22Rik, iPFK-2, uPFK-2}, PGK1 (phosphoglycerate kinase 1) [NCBI Gene 5230] {aka HEL-S-68p, MIG10, PGKA}, PLIN1 (perilipin 1) [NCBI Gene 5346] {aka FPLD4, PERI, PLIN}, Pck2 (phosphoenolpyruvate carboxykinase 2 (mitochondrial)) [NCBI Gene 74551] {aka 1810010O14Rik, 9130022B02Rik, PEPCK-M}, PIK3CB (phosphatidylinositol-4,5-bisphosphate 3-kinase catalytic subunit beta) [NCBI Gene 5291] {aka P110BETA, PI3K, PI3KBETA, PIK3C1}, GADD45B (growth arrest and DNA damage inducible beta) [NCBI Gene 4616] {aka GADD45BETA, MYD118}, Cebpa (CCAAT/enhancer binding protein alpha) [NCBI Gene 12606] {aka C/ebpalpha, CBF-A, Cebp}, PIK3R2 (phosphoinositide-3-kinase regulatory subunit 2) [NCBI Gene 5296] {aka MPPH, MPPH1, P85B, p85, p85-BETA, p85beta}, Fosb (Fos B proto-oncogene, AP-1 transcription factor subunit) [NCBI Gene 14282], PNPLA3 (patatin like domain 3, 1-acylglycerol-3-phosphate O-acyltransferase) [NCBI Gene 80339] {aka ADPN, C22orf20, iPLA(2)epsilon}, INS (insulin) [NCBI Gene 3630] {aka IDDM, IDDM1, IDDM2, ILPR, IRDN, MODY10}, Fyn (Fyn proto-oncogene, Src family tyrosine kinase) [NCBI Gene 14360], Pik3r1 (phosphoinositide-3-kinase regulatory subunit 1) [NCBI Gene 18708] {aka PI3K, p50alpha, p55alpha, p85alpha}, Pgk1 (phosphoglycerate kinase 1) [NCBI Gene 18655] {aka Pgk-1}, Mgll (monoglyceride lipase) [NCBI Gene 23945] {aka Magl, Mgl}, Crp (C-reactive protein, pentraxin-related) [NCBI Gene 12944], Nfkb1 (nuclear factor of kappa light polypeptide gene enhancer in B cells 1, p105) [NCBI Gene 18033] {aka NF-KB1, NF-kappaB, NF-kappaB1, p105, p50, p50/p105}, Hk2 (hexokinase 2) [NCBI Gene 15277] {aka HKII}, Mlxipl (MLX interacting protein-like) [NCBI Gene 58805] {aka ChREBP, Mlx, WS-bHLH, Wbscr14, bHLHd14}, Fos (Fos proto-oncogene, AP-1 transcription factor subunit) [NCBI Gene 14281] {aka D12Rfj1, c-fos, cFos}, G6pd2 (glucose-6-phosphate dehydrogenase 2) [NCBI Gene 14380] {aka G6pdx-ps1, Gpd-2, Gpd2}, PFKFB3 (6-phosphofructo-2-kinase/fructose-2,6-biphosphatase 3) [NCBI Gene 5209] {aka IPFK2, PFK2, iPFK-2}, TRBV20OR9-2 (T cell receptor beta variable 20/OR9-2 (non-functional)) [NCBI Gene 6962] {aka CDR3, TCRBV20S2, TCRBV2O, TCRBV2S2O}, G6PD (glucose-6-phosphate dehydrogenase) [NCBI Gene 2539] {aka CNSHA1, G6PD1}, LDHAP7 (lactate dehydrogenase A pseudogene 7) [NCBI Gene 100190800] {aka LDHAL7}, Nr1h4 (nuclear receptor subfamily 1, group H, member 4) [NCBI Gene 20186] {aka Fxr, HRR1, RIP14, Rxrip14}, VEGFA (vascular endothelial growth factor A) [NCBI Gene 7422] {aka L-VEGF, MVCD1, VEGF, VPF}, Pnpla3 (patatin-like phospholipase domain containing 3) [NCBI Gene 116939] {aka Adpn}, ALDOA (aldolase, fructose-bisphosphate A) [NCBI Gene 226] {aka ALDA, GSD12, HEL-S-87p}, HIF1A (hypoxia inducible factor 1 subunit alpha) [NCBI Gene 3091] {aka HIF-1-alpha, HIF-1A, HIF-1alpha, HIF1, HIF1-ALPHA, MOP1}, DUSP1 (dual specificity phosphatase 1) [NCBI Gene 1843] {aka CL100, HVH1, MKP-1, MKP1, PTPN10}, ENO2 (enolase 2) [NCBI Gene 2026] {aka HEL-S-279, NSE}, HK2 (hexokinase 2) [NCBI Gene 3099] {aka HKII, HXK2}, KRAS (KRAS proto-oncogene, GTPase) [NCBI Gene 3845] {aka 'C-K-RAS, C-K-RAS, CFC2, K-RAS2A, K-RAS2B, K-RAS4A}, Junb (jun B proto-oncogene, AP-1 transcription factor subunit) [NCBI Gene 16477], PFKP (phosphofructokinase, platelet) [NCBI Gene 5214] {aka ATP-PFK, PFK-C, PFK-P, PFKF}, ATF3 (activating transcription factor 3) [NCBI Gene 467], Aldoa (aldolase A, fructose-bisphosphate) [NCBI Gene 11674] {aka Aldo-1, Aldo1}, THBS1 (thrombospondin 1) [NCBI Gene 7057] {aka THBS, THBS-1, TSP, TSP-1, TSP1}, FYN (FYN proto-oncogene, Src family tyrosine kinase) [NCBI Gene 2534] {aka SLK, SYN, p59-FYN}, Tnf (tumor necrosis factor) [NCBI Gene 21926] {aka DIF, TNF-a, TNF-alpha, TNFSF2, TNFalpha, Tnfa}, Alb (albumin) [NCBI Gene 11657] {aka Alb-1, Alb1, BCL001, BCL002, BPL001}, PIK3R1 (phosphoinositide-3-kinase regulatory subunit 1) [NCBI Gene 5295] {aka AGM7, GRB1, IMD36, p85, p85-ALPHA, p85alpha}, FOS (Fos proto-oncogene, AP-1 transcription factor subunit) [NCBI Gene 2353] {aka AP-1, C-FOS, p55}, Atf3 (activating transcription factor 3) [NCBI Gene 11910] {aka LRG-21}, HSPA6 (heat shock protein family A (Hsp70) member 6) [NCBI Gene 3310] {aka HSP70B'}, Il6 (interleukin 6) [NCBI Gene 16193] {aka Il-6}, AKT1 (AKT serine/threonine kinase 1) [NCBI Gene 207] {aka AKT, PKB, PKB-ALPHA, PRKBA, RAC, RAC-ALPHA}, Ppara (peroxisome proliferator activated receptor alpha) [NCBI Gene 19013] {aka 4933429D07Rik, Nr1c1, PPAR-alpha, PPARalpha, Ppar}, MGLL (monoglyceride lipase) [NCBI Gene 11343] {aka HU-K5, HUK5, MAGL, MGL}, BCL2 (BCL2 apoptosis regulator) [NCBI Gene 596] {aka Bcl-2, PPP1R50}, UCP2 (uncoupling protein 2) [NCBI Gene 7351] {aka BMIQ4, SLC25A8, UCPH}, Il1b (interleukin 1 beta) [NCBI Gene 16176] {aka IL-1beta, Il-1b}
- **Diseases:** chronic liver disease (MESH:D008107), hypoxic (MESH:D002534), cytotoxicity (MESH:D064420), Parkinson's disease (MESH:D010300), bladder cancer (MESH:D001749), LIHC (MESH:D017093), autoimmune arthritis (MESH:D001168), autoimmune, cardiovascular, metabolic and cerebrovascular disorders (MESH:D024821), HCC (MESH:D006528), carcinogenic (MESH:D011230), myocardial ischemia (MESH:D017202), hepatic fibrosis (MESH:D008103), arthrogryposis (MESH:D001176), inflammation (MESH:D007249), Huntington's disease (MESH:D006816), fatty liver (MESH:D005234), hepatic metabolic dysregulation (MESH:D021081), metabolic disorder (MESH:D008659), cardiac muscle contraction (MESH:C536214), injury to (MESH:D014947), hepatic toxicity (MESH:D056486), glycogen storage disease (MESH:D006008), tumor necrosis factors (MESH:C536657), Insulin resistance (MESH:D007333), hypoxia (MESH:D000860), liver tumor (MESH:D008113), gastrointestinal cancers (MESH:D005770), primary immunodeficiencies (MESH:D000081207), metastasis (MESH:D009362), fibrosis (MESH:D005355), carcinogenesis (MESH:D063646), Cancer (MESH:D009369), NAFLD (MESH:D065626), hypertrophic cardiomyopathy (MESH:D002312)
- **Chemicals:** pyruvate (MESH:D019289), acetyl-CoA (MESH:D000105), lipopolysaccharide (MESH:D008070), ATP (MESH:D000255), fructose-2,6-bisphosphate (MESH:C027652), free fatty acids (MESH:D005230), PEP (MESH:D010728), glycerol (MESH:D005990), monoacylglycerols (MESH:D050178), triglyceride (MESH:D014280), ADP (MESH:D000244), pentose phosphate (MESH:D010428), D-glucose (MESH:D005947), lipid (MESH:D008055), dihydroxyacetone phosphate (MESH:D004099), mannose (MESH:D008358), fructose-1,6-bisphosphate (MESH:C029063), NADPH (MESH:D009249), calcium (MESH:D002118), shikimate (MESH:C000723335), 2-phosphoglycerate (MESH:C008885), cholesterol (MESH:D002784), ROS (MESH:D017382), amino (-), glyceraldehyde-3-phosphate (MESH:D005986), phosphate (MESH:D010710), fatty acid (MESH:D005227), TCA (MESH:D014233), carbon (MESH:D002244), lactate (MESH:D019344), oxygen (MESH:D010100), amino acid (MESH:D000596), fructose (MESH:D005632), 1,3-bisphosphoglycerate (MESH:C015891), Glyphosate (MESH:C010974), glycogen (MESH:D006003)
- **Species:** Homo sapiens (human, species) [taxon 9606], Mus musculus (house mouse, species) [taxon 10090], Rattus norvegicus (brown rat, species) [taxon 10116]
- **Mutations:** AUC of 0, AUC of 1, I148M
- **Cell lines:** C3A — Homo sapiens (Human), Hepatoblastoma, Cancer cell line (CVCL_1098), HepG2 — Homo sapiens (Human), Hepatoblastoma, Cancer cell line (CVCL_0027)

## Figures

8 figures with captions in the complete paper: https://tomesphere.com/paper/PMC13010713/full.md

---
Source: https://tomesphere.com/paper/PMC13010713