铁死亡相关的分子特征连接百草枯毒理基因组学扰动与帕金森病相关的黑质生物学

《Journal of Hazardous Materials Advances》：Ferroptosis-related molecular signatures linking paraquat toxicogenomic perturbation with Parkinson’s disease-relevant substantia nigra biology

　　

背景与目的：百草枯（paraquat）的氧化还原循环促进氧化应激和铁依赖性脂质过氧化，可能将其与帕金森病（Parkinson's disease, PD）中的铁死亡（ferroptosis）联系起来。研究人员优先选择了连接百草枯毒理基因组学证据与PD相关黑质（substantia nigra）生物学的候选基因。方法：对与百草枯、PD和铁死亡相关的基因进行交集分析。优先排序结合了供体分层的黑质差异表达、供体水平共表达和机器学习（machine learning）。所有36个可测量的候选基因直接进入LASSO（Least Absolute Shrinkage and Selection Operator）分析，未根据差异表达或GSE7621结果进行预选。内部性能通过供体分组的重复嵌套交叉验证（nested cross-validation）在GSE8397中评估。在GSE7621中比较了四种分类器；由于稀疏性和可解释性，保留了LASSO，然后在不重新训练的情况下在四个额外的大体组织队列（bulk cohorts）中评估。对单核（single-nucleus）和蛋白质组学（proteomic）数据进行了描述性评估。结果：交集鉴定出37个候选基因。GCH1、DDIT4和ATF4具有最强的供体感知边缘转录组学证据，并包含在固定的八基因LASSO评分中。重复供体分组的嵌套交叉验证在组织标本水平上产生了平均ROC-AUC（receiver operating characteristic-area under the curve）为0.898，在供体聚合后为0.899。固定评分在GSE7621中的ROC-AUC为0.854（95%置信区间，0.694–0.972），在四个额外的大体组织队列中为0.625–0.812。当在特征构建过程中省略先前的PD基因关联时，供体水平平均ROC-AUC为0.871。在错误发现率（false discovery rate, FDR）校正后，单核分析中的任何测试基因和任一优先集合中的任何定量蛋白质均未达到显著。结论：三个边缘候选基因嵌套在一个先验知识约束特征空间中的八基因多变量评分内。这些输出生成可测试的假设，但未建立机制联系或诊断效用。

**论文解读：基于铁死亡相关分子特征的百草枯毒理基因组学与帕金森病黑质生物学关联的计算研究**

**1. 研究背景与问题**

帕金森病（Parkinson's disease, PD）是第二常见的神经退行性疾病，其核心病理特征为黑质（substantia nigra）多巴胺能神经元的进行性退化和α-突触核蛋白（alpha-synuclein）的聚集。尽管遗传易感性贡献了疾病风险，但越来越多的证据表明环境毒物暴露是PD相关分子应激的重要修饰因子。百草枯（paraquat）是一种广泛使用的除草剂，其氧化还原循环机制生成活性氧（reactive oxygen species, ROS），引发线粒体损伤，并在实验模型中诱导多巴胺能神经毒性。流行病学荟萃分析报告了百草枯暴露与PD风险增加之间的关联，但暴露评估、共暴露和观察性设计的局限性要求谨慎解释因果关联。

铁死亡（ferroptosis）是一种铁依赖性的、调控性细胞死亡形式，由磷脂过氧化和细胞过氧化物解毒系统失败驱动。百草枯的氧化还原循环可能通过超越一般氧化应激的途径与铁死亡交汇：持续的过氧化物产生增加氧化负担，消耗谷胱甘肽（glutathione），削弱谷胱甘肽依赖性抗氧化防御，并促进铁催化的脂质自由基传播。这些过程与PD中涉及的线粒体功能障碍、铁失调、自噬和脂质重塑相交。然而，连接百草枯相关毒理基因组学信号与人类PD相关黑质生物学的特异性铁死亡相关基因尚未完全明确。

**2. 研究目的与内容**

研究人员开展本项研究，旨在整合毒理基因组学、疾病基因和铁死亡资源，结合人死后黑质转录组学、加权基因共表达网络分析（WGCNA）、机器学习建模、单核转录组学、动物暴露-反应数据、蛋白质组学和计算机模拟分析，以优先识别连接百草枯毒理基因组学扰动与PD相关黑质生物学的候选基因。论文发表在《Journal of Hazardous Materials Advances》。研究结论是：识别出三个边缘候选基因（GCH1、DDIT4和ATF4）嵌套在一个八基因多变量评分内，该评分在知识约束的特征空间中生成可测试假设，但未建立机制联系或诊断效用。其重要意义在于为环境神经毒理学和PD风险评估提供了可验证的分子假设，并展示了计算整合在直接人类暴露转录组学不可用时的价值与局限。

**3. 主要关键技术方法（不超过250字）**

本研究完全基于公开数据，主要关键技术方法包括：（1）基因集构建：从Comparative Toxicogenomics Database（CTD）获取百草枯相关基因，从GeneCards、DisGeNET和CTD获取PD相关基因，从FerrDb获取铁死亡相关基因，三者交集生成37个知识约束候选基因；（2）人黑质转录组学分析：使用GSE8397队列（23个生物学供体，15例PD和8例对照）进行差异表达分析（limma，以供体为分组变量），并进行WGCNA；（3）机器学习建模：将所有36个可测量候选基因直接输入LASSO逻辑回归，采用20次重复供体分组嵌套交叉验证评估内部性能，并在GSE7621中进行模型族比较（LASSO、弹性网、随机森林、线性支持向量机），最终固定八基因及系数；（4）外部验证：在四个额外大体组织队列（GSE49036、GSE114517、GSE20292、GSE168496）中评估固定评分，并在GSE182622激光捕获显微切割队列中进行敏感性分析；（5）多组学上下文分析：使用GSE157783单核数据集、GSE36232小鼠暴露数据、PXD037684蛋白质组数据以及计算机扰动网络分析。

**4. 研究结果**

**4.1 三向交集定义了一个37基因知识约束候选集**
通过整合百草枯相关基因（4,841个）、PD相关基因（1,213个）和铁死亡相关基因（104个），三向交集获得37个候选基因。KEGG富集分析显示铁死亡为最显著通路（9/37基因，FDR校正P=4.44×10^-12），其次是自噬、线粒体自噬、神经退行性疾病和帕金森病等。GO分析包括氧化应激、化学应激、氧水平、金属离子和活性氧代谢等反应。这些注释描述了候选集的内部组成，但PD和铁死亡注释部分由构建过程预期，不构成独立验证。

**4.2 边缘转录组学和共表达证据优先排序了候选基因**
在GSE8397中，基于供体感知的差异表达分析，三个具有最小FDR的候选基因被优先排序：GCH1（log₂FC=-1.797，FDR=0.00151）、DDIT4（log₂FC=1.151，FDR=0.00151）和ATF4（log₂FC=0.640，FDR=0.00313）。WGCNA显示，包含DDIT4和ATF4的棕色模块与PD正相关（r=0.705，FDR=0.00111），而包含GCH1的绿松石模块与PD负相关（r=-0.715，FDR=0.00111）。这三个边缘候选均被后续LASSO保留。

**4.3 八基因LASSO评分在队列间表现出异质性性能**
所有36个可测量候选基因未经过滤直接进入LASSO训练。在23个生物学供体的20次重复供体分组嵌套交叉验证中，平均标本水平ROC-AUC为0.898（SD=0.038），供体聚合后为0.899（SD=0.037）。最终模型保留DDIT4、NFE2L2、BECN1、GSR、ATF4、ATG7、SLC7A11和GCH1。在GSE7621中，LASSO的ROC-AUC为0.854（95% CI, 0.694–0.972），但随机森林最高（0.917），LASSO因稀疏性和可解释性被保留。固定评分在四个额外大体组织队列中的ROC-AUC范围为0.625–0.812，GSE168496中为0.625（95% CI, 0.328–0.891）。在GSE182622细胞富集制备中，评分未区分PD与对照（ROC-AUC=0.467）。敏感性分析（省略PD基因关联）显示类似分离，但未提供独立验证。

**4.4 单核转录组学提供了细胞表达上下文**
在GSE157783中，八基因在注释的多巴胺能神经元中均有检测，但供体水平的PD-对照比较均未通过FDR校正（最低名义P为GCH1：P=0.052，FDR=0.415），仅提供细胞表达上下文。

**4.5 小鼠百草枯暴露数据提供了描述性跨物种上下文**
在GSE36232中，八基因在腹侧中脑和纹状体中有部分检测，暴露后平均log₂FC为-0.051（腹侧中脑）和-0.206（纹状体），描述性结果提示暴露响应性，但未直接复制人PD方向。

**4.6 人黑质蛋白质组学对优先基因集的评估**
在PXD037684发现集（15例PD/15例对照）中，定量了8个LASSO蛋白中的5个（包括GCH1），但FDR校正后无一显著（最小名义P为GSR：P=0.013，q=0.081）。复制集（10例PD/9例对照）中未包含任何边缘候选。更广泛的37个候选基因空间中有4个蛋白（HSPB1、ACO1、KEAP1、IDH1）在发现集中q<0.05，但均不属于优先集。未进行集合级推断。

**4.7 计算机扰动和通路分析提供了探索性上下文**
网络扰动、通路评分和化学扰动匹配为5个LASSO基因（DDIT4、GCH1、ATG7、NFE2L2、BECN1）提供了预测性扰动效应，但结果来自模型，非实验干预。

**5. 讨论总结与结论翻译**

**讨论总结：** 本研究的主要输出是三个具有最强边缘转录组学证据的基因（GCH1、DDIT4、ATF4）和一个八基因条件LASSO评分。两种程序回答不同统计问题，应视为互补优先排序视角。研究设计的关键特征包括：将知识约束构建与下游优先排序分离、通过供体分组处理重复标本、在嵌套交叉验证中分组重复标本、报告固定系数和稳定性。敏感性分析表明，即使省略PD基因关联，分离仍存在，但无法区分泛应激或组织组成信号。DDIT4（RTP801/REDD1）在PD黑质神经元中升高，其生物学背景依赖，可抑制mTORC1，但可能改变自噬通量和抗氧化适应。GCH1调节四氢生物蝶呤生物合成，连接多巴胺产生和氧化还原生物学，但蛋白质组学未显示转录-蛋白一致性。单核分析、小鼠数据和蛋白质组学均未提供FDR显著支持。研究局限性包括：先验知识引入的循环性、缺乏人类百草枯暴露黑质转录组、训练队列仅23个供体、外部队列性能异质性、蛋白质组学覆盖不全。未来方向包括在人类iPSC衍生的多巴胺能神经元中进行前瞻性实验，靶向优先候选基因，并评估脂质过氧化等指标。

**研究结论翻译：** 本研究完全基于公开数据，在一个百草枯-PD-铁死亡知识约束空间中，识别出三个供体感知的边缘候选基因，嵌套在一个八基因多变量评分内。该评分在大体组织队列中表现出异质性性能，在细胞富集的GSE182622敏感性分析中未显示区分能力。结合选择不稳定性和任一优先集缺乏FDR显著蛋白水平证据，这些发现支持假设生成，而非机制性或诊断性主张。需要在人类相关神经元系统中进行受控实验。