IF50.5｜农药残留重塑土壤生物多样性与功能基因

欧洲373个位点揭示农药是仅次于土壤性质的第二大驱动因子。

这篇按生信解读号的读法拆：先看背景和问题，再看作者怎么做、做出了哪几条硬发现，最后把数据和代码收走。

研究背景

土壤生物多样性对维持生态系统功能至关重要，但农药在土壤中广泛残留，其对地下生物群落的影响尚缺乏大尺度评估。已有研究表明农药对鸟类、蜜蜂等地上生物有负面影响，而对土壤生物的研究多局限于特定区域、少数物种或单一农药，难以反映真实田间多种农药共存下的复杂效应。LUCAS Soil 2018调查显示欧洲87%的土壤样本至少检出一种农药残留，凸显了评估农药对土壤生物多样性影响的紧迫性。

科学问题

卡在哪里： 此前缺乏在洲际尺度上同时涵盖多种农药、多个生物类群和功能基因的系统研究，且现有风险评估依赖少数模式物种和单一化合物测试，无法捕捉田间复杂群落对农药的响应。

本文要回答： 在欧洲大陆尺度上评估多种农药残留对土壤古菌、细菌、真菌、原生生物、线虫和节肢动物及关键功能基因多样性的影响，并量化农药相对于环境因子的重要性。

技术路线

作者主要用了：DNA宏条形码、宏基因组学、广义线性模型、变异分解、VIP分析。

研究利用LUCAS Soil 2018调查中373个位点的土壤样本，这些位点同时测定了农药残留和土壤生物多样性。通过DNA宏条形码技术获取古菌、细菌、真菌、原生生物、线虫和节肢动物的群落组成，并计算物种丰富度和Shannon多样性；通过宏基因组学分析48个与碳、氮、磷循环相关的功能基因组的多样性。随后构建广义线性模型，纳入农药浓度、土壤性质、气候和生态系统类型作为预测变量，通过逐步回归筛选重要变量，并利用变异分解量化各驱动因子的相对贡献。

图 1　Fig. 1。图1展示了研究框架：利用广义线性模型结合农药浓度、土壤性质、气候和生态系统类型，分析农药对土壤生物多样性（分类多样性、功能群和功能基因）的影响，并量化农药相对于环境变量的贡献，同时比较农田与所有生态系统类型中的效应。

核心亮点

亮点 1｜农药残留普遍存在且类型多样

在373个位点中，70%的位点至少检出一种农药残留，共检测到63种不同的农药。其中杀菌剂占54%，除草剂占34.9%，杀虫剂占11.1%。草甘膦及其代谢物AMPA是最常检出的农药，其次是啶酰菌胺、二甲戊灵和氟环唑。农药残留数量和累积浓度在农田中最高，其次是草地和林地。

图 2　Fig. 2。图2a显示了373个采样位点的分布，涵盖5种生态系统类型；图2b显示检测到的农药类型中杀菌剂占54%，除草剂占34.9%，杀虫剂占11.1%；图2c箱线图显示农田中检测到的农药残留数量最高，其次是永久农田、近期转化草地、粗放草地和林地。

亮点 2｜农药是土壤生物多样性的第二大驱动因子

变异分解显示，在农田中，农药浓度解释了土壤生物多样性变异的最高达29.5%，仅次于土壤性质，高于气候和生态系统类型。对于真菌丰富度，农药单独解释了12.3%的变异，是首要驱动因子。在48个功能基因组中，有20个基因组的多样性受农药影响程度与土壤性质相当或更高。

图 3　Fig. 3。图3展示了农药浓度与土壤生物多样性指标之间的VIP网络。节点代表不同的生物类群和功能基因群，连线颜色表示关联方向。真菌丰富度与多种杀菌剂和除草剂呈负相关，而某些功能基因群（如细菌降解基因）与农药浓度呈正相关。

亮点 3｜农药对非靶标生物产生复杂影响

农药对非靶标生物的影响具有生物类群特异性和农药种类特异性。例如，杀菌剂联苯吡菌胺浓度与原生生物、线虫和节肢动物丰富度下降以及古菌、细菌和节肢动物多样性降低相关。除草剂二甲戊灵降低了有益菌根真菌和食细菌线虫的相对丰度，却促进了植物寄生原生生物。草甘膦浓度升高与真菌植物病原体和植食性线虫相对丰度增加相关，同时降低了原生生物和线虫丰富度。

图 4　Fig. 4。图4展示了变异分解结果，比较了农药残留、土壤性质、气候和生态系统类型对土壤生物多样性指标的解释方差百分比。在农田中，农药对真菌丰富度、细菌化学异养菌、丛枝菌根真菌、原生生物植物寄生生物、植食性线虫和食细菌线虫的相对丰度有显著贡献，最高达29.5%。

亮点 4｜农药改变碳氮磷循环相关功能基因多样性

宏基因组分析显示，农药对参与反硝化和几丁质降解的细菌基因影响最为显著。在农田中，70%的细菌基因群和84%的真菌基因群多样性对农药呈正响应，而古菌和动物基因群中约50%的关联为负。例如，草甘膦代谢物AMPA浓度与细菌膦酸盐降解基因多样性正相关，表明细菌可能利用草甘膦作为磷源。

亮点 5｜农药效应在农田中更强且存在溢出效应

当仅分析农田时，农药对土壤生物多样性的解释力高于包含所有生态系统类型的分析（功能基因解释力提高5.5倍）。然而，在草地和林地中也检测到农药残留，且某些关联（如二甲戊灵与真菌植物病原体和植食性线虫的正相关）仅在纳入非农田生态系统时出现，表明农药可能从农田漂移到周边生态系统并产生影响。

生信可带走

这一块是给做分析的人用的，不是科普点缀。

• 方法栈： DNA宏条形码、宏基因组学、广义线性模型、变异分解、VIP分析
• 公开数据： 原文未给出公开组学登录号
• 代码： 原文未给出公开 GitHub/GitLab 仓库
• 谁该点开原文： 从事土壤生态学、微生物组学、农药风险评估或环境生物信息学的研究者。

带走一句

做土壤微生物组分析时，别只盯着土壤理化性质，农药残留是解释群落变异的关键变量，应纳入模型。