《Journal of Xenobiotics》：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
Divya Mishra
编辑推荐：
草甘膦（Glyphosate）是全球使用最广泛的除草剂之一，其潜在的肝毒性及与肝癌生物学相关的分子变化引发了广泛关注。本研究采用基于RNA测序（RNA-Seq）的计算机系统生物学方法，阐明草甘膦诱导的肝癌细胞差异转录谱。分析揭示了显著的转录谱特征，上调的枢纽基
草甘膦（Glyphosate）是全球使用最广泛的除草剂之一，其潜在的肝毒性及与肝癌生物学相关的分子变化引发了广泛关注。本研究采用基于RNA测序（RNA-Seq）的计算机系统生物学方法，阐明草甘膦诱导的肝癌细胞差异转录谱。分析揭示了显著的转录谱特征，上调的枢纽基因（hub genes）包括ATF3、JUNB、ALDOA、FOSB、PFKFB3、G6PD、ENO2、HK2、FOS和PGK1，这些基因主要与葡萄糖代谢、TNF-α/NF-κB信号、上皮-间质转化（EMT）及细胞应激反应相关。相反，多个关键基因显著下调，包括PIK3R1、FYN、CEBPA、MLXIPL、PPARA、CD36、PCK2、PNPLA3、NR1H4和MGLL，它们参与脂质代谢、免疫调节和非酒精性脂肪肝病（NAFLD）通路。值得注意的是，所有枢纽基因均表现出较强的诊断性能，凸显了其作为草甘膦暴露敏感生物标志物的潜力。该研究全面揭示了草甘膦暴露相关的基因表达变化，将其与肝脏代谢失调和免疫调节联系起来，并提出了一组具有潜在诊断和治疗意义的枢纽基因面板。
论文解读：草甘膦诱导的肝癌细胞分子机制与生物标志物鉴定
研究背景与动机
草甘膦作为全球广泛使用的广谱除草剂，其安全性评估长期处于争议之中。尽管传统观点认为其作用于植物特有的莽草酸途径，对人类相对安全，但国际癌症研究机构（IARC）将其列为可能的人类致癌物，挑战了这一假设。肝脏作为异生物质代谢和解毒的核心器官，是草甘膦毒性作用的潜在靶标。既往研究表明，草甘膦暴露可诱发动物模型出现肝功能障碍、脂代谢紊乱及氧化应激，且在慢性肝病和肝癌患者体内检测到较高水平的草甘膦残留。然而，草甘膦在肝癌细胞中诱导的具体分子改变及其作为生物标志物的潜力仍缺乏系统性探索。为此，研究人员开展了此项研究，旨在通过计算机系统生物学方法，解析草甘膦暴露下肝癌细胞的转录组特征，识别关键调控基因，并评估其诊断与治疗价值。该研究成果已发表于《Journal of Xenobiotics》。
关键技术方法概述
研究人员从NCBI Gene Expression Omnibus (GEO) 检索并下载了草甘膦处理的肝细胞瘤C3A细胞系RNA-Seq原始数据（登录号GSE292868），包含3个处理组和6个对照组样本。利用DESeq2软件包进行差异表达基因（DEGs）分析，筛选标准为调整后p值（padj）<0.05且|log2倍数变化（FC）|>1。随后，利用STRING数据库和Cytoscape软件构建蛋白质-蛋白质相互作用（PPI）网络，并通过MCODE和CytoHubba插件识别关键模块和枢纽基因。利用基因集富集分析（GSEA）和DAVID数据库进行功能与通路富集分析。此外，通过GEPIA2平台利用TCGA-LIHC数据集验证枢纽基因在肿瘤与正常组织中的表达差异，并从比较毒理基因组学数据库（CTD）获取草甘膦相关基因，最后利用受试者工作特征（ROC）曲线评估枢纽基因的诊断效能。
研究结果
3.1. 质量检查与预处理分析
通过对原始RNA-Seq数据进行质量评估，所有样本的测序深度充足，碱基质量高（≥Q30），GC含量分布一致，且未检测到接头污染，表明数据质量可靠，适用于下游转录组分析和系统生物学研究。
3.2. 差异表达基因（DEGs）分析
对比草甘膦处理与未处理的肝癌细胞，共鉴定出1569个显著差异表达基因，其中包含841个上调基因和728个下调基因。热图显示处理组与对照组基因表达谱明显分离。其中，蛋白编码基因LDHAP7上调最为显著（log2FC=6.0），而miRNA编码基因FP671120.6下调最为显著（log2FC=-5.7）。
3.3. 蛋白质-蛋白质相互作用（PPI）网络分析
构建PPI网络后识别出关键枢纽基因。上调枢纽基因为FOSB、PFKFB3、ENO2、HK2、JUNB、FOS、G6PD、ATF3、ALDOA和PGK1；下调枢纽基因为FYN、CD36、MGLL、PIK3R1、CEBPA、PCK2、NR1H4、MLXIPL、PNPLA3和PPARA。这些基因大多富集于网络中识别出的关键模块内。
3.4. 枢纽基因在肿瘤与正常肝组织中的表达分析
利用TCGA-LIHC数据验证发现，PFKFB3和HK2在肿瘤组织中也呈高表达，与草甘膦处理组数据一致。然而，ATF3、ENO2和FOS在肿瘤组织中表达较低，与草甘膦诱导的上调趋势相反。大多数下调枢纽基因在肿瘤与正常肝组织间未显示出显著差异，表明草甘膦诱导的转录改变主要反映了体外暴露驱动的分子变化，而非一致的肿瘤特征性表达。
3.5. 基因集富集分析（GSEA）结果
GSEA分析显示，上调的标志性通路包括上皮-间质转化（EMT）、KRAS信号、肌生成、TNF-α/NF-κB信号和缺氧通路。KEGG分析表明核糖体、心肌收缩等通路上调，而甘油酯代谢、丁酸代谢等通路下调。枢纽基因如FOS、ATF3、PFKFB3等被证实与这些通路密切相关。
3.6. 生物功能与疾病富集分析
GO和KEGG分析进一步揭示，上调基因主要富集于葡萄糖代谢（糖酵解、糖异生）、HIF-1信号通路及碳代谢；下调基因则显著富集于免疫反应、细胞黏附、紧密连接以及非酒精性脂肪肝病（NAFLD）和胰岛素抵抗通路。疾病关联分析显示，这些基因与自身免疫病、心血管疾病及肝脏代谢紊乱密切相关。
3.7. 草甘膦相关基因的检索
通过与CTD数据库比对，研究发现FOS、PGK1和PIK3R1不仅是枢纽基因，同时也是已知的草甘膦关联基因，这进一步证实了研究结果的可靠性。
3.8. 免疫相关基因集的检索
分析表明，多数上调及下调枢纽基因均属于免疫相关基因，提示草甘膦暴露可能通过干扰固有免疫和适应性免疫相关基因的表达，导致免疫微环境失衡。
3.9. 枢纽基因的诊断性能
ROC曲线分析显示，所有10个上调枢纽基因（如FOSB、PFKFB3等）的曲线下面积（AUC）均为1.000，表现出极佳的诊断性能。下调枢纽基因中，除MGLL的AUC为0.944外，其余（如CEBPA、FYN等）AUC亦均为1.000，表明这些基因具有作为草甘膦暴露诊断生物标志物的巨大潜力。
结论与讨论总结
本研究通过整合转录组学与系统生物学方法，揭示了草甘膦暴露诱导肝癌细胞发生显著代谢重编程与免疫调节的分子机制。研究发现，草甘膦处理导致以HK2、PFKFB3、ENO2等为核心的糖酵解和HIF-1信号通路异常激活，同时抑制了以PPARA、CEBPA、CD36为代表的脂质代谢和PPAR信号通路。这种“代谢转换”特征与肿瘤的Warburg效应及缺氧适应机制高度相似。尽管部分枢纽基因在TCGA临床样本中未表现出一致的差异，但其在草甘膦处理模型中的特异性改变及其完美的诊断效能（AUC值接近1），使其有望成为区分草甘膦暴露与肝癌发生发展的特异性生物标志物。研究不仅为草甘膦的肝毒性机制提供了新的视角，也为环境化学物暴露相关疾病的早期诊断提供了潜在的分子靶标。
