《Emerging Contaminants》：Manure inputs are strongly associated with soil microbial community structure in pastoral farmland: disentangling shared variation with antibiotics and microplastics
编辑推荐：
随着畜牧生产的集约化，畜禽粪便已成为抗生素和微塑料进入农业土壤的主要途径。然而，粪便有机肥、抗生素和微塑料对土壤微生物群落的相对贡献仍不明确。本研究调查了中国西北地区某大型养牛场周边的十块牧区农田。利用粪甾醇分子标志物表征牛粪投入的相对强度，同时量化了区域普遍
随着畜牧生产的集约化，畜禽粪便已成为抗生素和微塑料进入农业土壤的主要途径。然而，粪便有机肥、抗生素和微塑料对土壤微生物群落的相对贡献仍不明确。本研究调查了中国西北地区某大型养牛场周边的十块牧区农田。利用粪甾醇分子标志物表征牛粪投入的相对强度，同时量化了区域普遍存在的抗生素和土壤微塑料，并将其与微生物群落结构相关联。结果表明，豆甾烷醇浓度与土壤电导率（EC）、全氮（TN）、有机质（OM）和有效磷（AP）呈显著正相关，证实了豆甾烷醇作为牛粪投入替代指标的可靠性。12种目标抗生素的总和浓度（Σ12抗生素）与粪便投入梯度无显著关联；仅培氟沙星（Pef）、环丙沙星（Cip）和恩诺沙星（Enr）与豆甾烷醇呈显著正相关，而四环素类则与粪便替代指标无显著关联。微塑料主要以细颗粒（<50 μm）为主，聚合物组成的变化表明其来源既有粪便相关也有非粪便来源。主坐标分析（PCoA）显示，粪便投入较高的样点沿PCo1轴（解释45.6%变异）独立聚类。单变量基于距离的冗余分析（db-RDA）中，豆甾烷醇是微生物群落结构最强的实测相关因子（校正R2 = 0.394, p = 0.003）。方差分解分析（VPA）进一步显示，与粪便投入、所选抗生素替代指标和微塑料丰度相关的解释部分存在大量重叠：培氟沙星的边际关联大部分与粪便投入替代指标共享，而与微塑料相关的独特部分仍处于探索性阶段，因其单变量模型不显著。本研究确立了粪便投入在塑造牧区农田土壤微生物群落结构分化中的主导作用，排除了当前残留浓度下抗生素独立的群落水平效应，并为后续研究粪便-抗生素-微塑料相互作用的控制实验提供了清晰的研究假设和优先方向。
研究背景与科学问题
随着畜牧生产集约化程度的不断提高，大量畜禽粪便的产生与处置已成为全球性的环境挑战。畜禽粪便富含有机质及氮、磷等养分，施用于农田作为有机肥料是一种广泛采用的资源化利用策略，对提升土壤肥力和改善土壤结构具有重要作用。然而，集约化养殖中抗生素的大量使用以及饲料和养殖设施中塑料材料的普遍存在，意味着畜禽粪便在作为养分载体进入农田的同时，也构成了抗生素和微塑料等新兴污染物进入土壤环境的主要途径。抗生素进入土壤后，不仅可能被作物吸收并通过食物链传递，还会对土壤微生物群落产生选择压力，诱导抗生素抗性基因的富集与扩散。微塑料则可能作为其他污染物的载体，并通过改变土壤物理结构和养分有效性间接影响土壤微生物组成。
在此背景下，一个核心的科学问题凸显出来：在粪便还田的实际农田场景中，粪便有机肥本身带来的营养效应、共输入的抗生素的生态毒理效应，以及微塑料的物理化学效应，这三者对土壤微生物群落结构的相对贡献究竟如何？由于在实际田间条件下，粪便投入、抗生素和微塑料往往高度共现且相互关联，传统的观测研究很难将它们的效应进行有效分离。因此，亟需一种能够定量表征粪便投入强度的可靠方法，并结合合理的统计分析框架，以区分这些因素对微生物群落的独立与共享影响。
研究内容与主要结论
为解决上述问题，本研究聚焦于中国西北地区某大型养牛场周边的十块牧区农田。研究人员创新性地引入了粪甾醇分子标志物，特别是豆甾烷醇（Stigmastanol, Sti），来定量表征各田块间牛粪投入的相对强度。同时，系统分析了土壤中多种抗生素和微塑料的赋存特征，并利用宏基因组测序技术解析了土壤微生物群落结构。通过一系列相关性分析和方差分解分析（Variance Partitioning Analysis, VPA），研究人员旨在定量分离粪便投入、抗生素和微塑料三者对微生物群落结构变异的独立与共享贡献。
研究得出的核心结论是：粪便投入是驱动牧区农田土壤微生物群落结构分化的主导因素。具体而言，豆甾烷醇浓度与土壤养分和理化性质紧密相关，并且是微生物群落结构最强的单一相关因子。相比之下，在当前研究的残留浓度水平下，未发现抗生素存在独立的群落水平效应；培氟沙星（Pefloxacin, Pef）与微生物群落的关联几乎完全源于其与粪便投入梯度的强烈共变。微塑料的独立贡献虽在方差分解中有所体现，但由于其单变量模型不显著，仍需审慎对待。
研究的重要意义
本研究的科学意义在于，它首次在一个实际田间尺度上，利用创新的分子标志物方法，定量揭示了粪便投入相对于其携带的抗生素和微塑料，在塑造土壤微生物群落结构中的主导地位。这一发现为评估畜禽粪便农用的综合生态效应提供了关键的定量依据。更重要的是，研究结果清晰地指出了后续研究的方向：未来的控制实验必须将粪便基质本身的效应与其中特定污染物的效应进行独立操纵，才能验证本研究中观察到的统计模式是否对应着真实的因果生态效应。从管理角度看，减少畜牧生产和粪便处理过程中的不必要塑料输入，结合改进的粪便处理和土壤有机质管理措施，可能有助于限制粪便伴生污染物的迁移和环境持久性。
主要关键技术方法
- 1.
粪甾醇分子标志物分析法：采用超高效液相色谱-串联质谱法（UPLC-MS/MS）测定土壤中粪甾醇（如豆甾烷醇Sti、粪甾醇Cop、β-谷甾醇b-Sit）含量，用以定量表征不同田块间牛粪投入的相对强度。
- 2.
抗生素检测方法：采用液相色谱-串联质谱法（LC-MS/MS）定量分析土壤中12种目标抗生素（主要为氟喹诺酮类和四环素类）的浓度。
- 3.
微塑料分析方法：采用密度分离（碘化钠溶液）、过氧化氢消解结合激光直接红外化学成像系统（LDIR）鉴定土壤中微塑料的丰度、尺寸分布和聚合物组成。
- 4.
宏基因组测序与生物信息学分析：提取土壤微生物总DNA，构建宏基因组鸟枪法文库并在Illumina NovaSeq平台测序，利用Kraken2软件对序列进行物种分类注释。
- 5.
多元统计分析方法：运用Spearman秩相关、偏相关分析、主坐标分析（PCoA）、单变量距离冗余分析（db-RDA）和方差分解分析（VPA）等方法，解析粪便投入、抗生素、微塑料与微生物群落结构之间的关联及其相对贡献。所有样品均采集自中国西北地区某大型养牛场周边的十块农田。
研究结果
3.1 粪甾醇浓度与基本土壤理化性质
3.1.1 粪甾醇浓度与污染源识别：通过对十个样点土壤中粪甾醇浓度的测定发现，S1和S2样点的粪甾醇（Cop）、豆甾烷醇（Sti）和β-谷甾醇（b-Sit）浓度显著高于其他样点，形成高浓度簇。所有样点的Cop/(Cop+Sti)比值均低于0.3，证实了草食动物（牛）粪便信号的主导地位。S1和S2样点较高的S/(S+b)比值（Sti/(Sti+b-Sit)）进一步表明这两个田块接受了更高的牛粪投入量。即使在低投入组（S3-S10）内部，S/(S+b)比值也存在可分辨的梯度，反映了施肥频率、用量或时间的差异。
3.1.2 土壤理化性质及其与粪甾醇的关系：豆甾烷醇（Sti）浓度与土壤电导率（EC）、全氮（TN）、有机质（OM）和有效磷（AP）呈显著正相关，与pH呈显著负相关。这证实了Sti作为牛粪投入替代指标的可靠性，并表明粪便投入改变了土壤化学环境，进而可能影响共存污染物的环境行为及微生物群落。
3.2 土壤中的抗生素
3.2.1 抗生素的赋存特征：所有12种目标抗生素均在S1和S2样点检出。四环素类抗生素浓度远高于氟喹诺酮类。值得注意的是，Σ12抗生素总浓度的空间分布并未遵循粪便投入梯度，高投入样点S1和S2的总浓度仅处于中等水平，而最高总浓度出现在低投入样点（S9、S5、S3）。这表明牛粪投入并非控制该区域抗生素残留空间分布的唯—因素，可能存在废水灌溉、历史残留等其他来源。
3.2.2 甾醇与抗生素的相关性：在控制土壤有机质和pH后进行偏相关分析，结果显示三种氟喹诺酮类药物——培氟沙星（Pef）、环丙沙星（Cip）和恩诺沙星（Enr）与豆甾烷醇（Sti）呈显著正相关。而所有四环素类抗生素与Sti均无显著正相关，部分甚至呈负相关。这揭示了两种抗生素类别与牛粪投入截然不同的关联模式：氟喹诺酮类与粪便投入的空间耦合更强，而四环素类的分布则受到更复杂因素的影响。
3.3 土壤中的微塑料
3.3.1 微塑料的丰度、组成与粒径分布：所有样点均检出微塑料，总丰度范围为140至755 items/g。S1和S2样点总丰度最高。微塑料粒径分布以<50 μm的细小颗粒占绝对优势（超过80%）。聚合物组成方面，共检出9种类型，其中聚酰胺（PA）分布最广，聚乙烯（PE）丰度最高。S1和S2样点检出的聚合物类型最多，表明粪便投入可能同时引入多种来源的微塑料。
3.3.2 微塑料与粪甾醇的关联：豆甾烷醇（Sti）与<50 μm微塑料丰度呈显著正相关，但与50-100 μm和>100 μm的较大颗粒无显著相关。在聚合物层面，Sti与聚碳酸酯（PC）、聚四氟乙烯（PTFE）和聚氯乙烯（PVC）的丰度呈显著正相关。这种选择性关联表明，牛粪投入是<50 μm细小颗粒和某些特定聚合物的一个重要来源，而聚乙烯（PE）等聚合物则更多来自地膜残留等其他非粪便途径。
3.4 抗生素与微塑料复合污染的耦合特征：通过交叉分析抗生素与微塑料的相关性，并结合它们与粪便标志物的关系，可以识别出三种耦合模式：一是粪便关联共现模式（如培氟沙星Pef与聚对苯二甲酸乙二醇酯PET共同与Sti相关）；二是非粪便关联模式（如聚乙烯PE与多数抗生素呈负相关，且其高值出现在低投入样点）；三是输入后行为分歧模式（四环素类与微塑料的相关性复杂，表明即使共输入，其后继环境归趋也可能截然不同）。
3.5 土壤微生物群落的响应
3.5.1 微生物群落组成与多样性：α多样性指数在各样点间重叠，表明长期牛粪施用并未显著增加或降低物种丰富度和均匀度。然而，主坐标分析（PCoA）显示，高投入样点S1和S2在PCo1轴上与其他样点明显分离，表明粪便投入梯度与微生物群落结构的空间分化存在强关联。例如，放线菌（Actinobacteria）在高投入样点的相对丰度较低，可能与新鲜有机质的持续补充抑制了寡营养型分解者的竞争优势有关。
3.5.2 粪便投入、抗生素和微塑料对微生物群落变异的贡献划分：单变量距离冗余分析（db-RDA）显示，豆甾烷醇（Sti）单独解释了39.4%的群落变异（最强），培氟沙星（Pef）解释了27.6%，而总微塑料丰度（TMP）的解释度仅为6.9%且不显著。方差分解分析（VPA）表明，三者共同解释了35.7%的变异。重要的是，Pef的独特贡献几乎为零（0.5%），说明其显著的边际关联完全源于与Sti的共变，而非直接的抗生素毒性效应。TMP的独特贡献（8.1%）虽最大，但其单变量模型不显著，因此只能视为初步探索性结果。总体而言，粪便投入相关的变异与抗生素和微塑料的变异存在大量重叠，表明在观测研究中，明确考虑粪便源协方差对于正确解释新兴污染物与微生物群落的关联至关重要。
讨论与结论总结
讨论部分强调，粪便投入可能通过多种途径影响微生物群落，包括提供有机碳和养分、改变土壤理化条件，而这些变化又会影响抗生素的吸附和持久性以及微塑料的行为。微塑料可通过改变土壤孔隙度、持水性和团聚体稳定性，并提供生物膜附着表面，从而间接影响微生物。这些并发过程解释了为何粪便投入替代指标表现出最强的边际关联，以及三者的变异存在大量重叠。研究同时指出，由于样本量有限（n=10），统计效力较低，尤其对弱关联的检测能力不足。因此，需要未来通过独立操纵粪便、抗生素和微塑料的控制实验，来验证这些统计模式是否对应真实的因果生态效应。
研究结论翻译：本研究利用粪甾醇分子标志物表征了牧区农田牛粪投入的相对梯度，并考察了其与抗生素、微塑料和土壤微生物群落结构的关联。豆甾烷醇与土壤养分和理化性质密切相关。几种氟喹诺酮类药物与豆甾烷醇共变，而四环素类与粪便梯度无显著关联。微塑料丰度和聚合物组成进一步表明其来源既有粪便相关，也有非粪便农业源。在所选的预测因子中，粪便投入替代指标与微生物群落结构具有最强的边际关联，而方差分解显示，粪便投入、培氟沙星和总微塑料丰度之间存在大量的共享解释变异。与培氟沙星相关的群落变异无法在统计上与粪便投入梯度所共享的变异分开，而与微塑料相关的独特部分仍处于探索性阶段，因为其边际关联不显著。在类似的区域农田系统中，未能找到真正的零粪便对照田块，因为牛粪历来被广泛用于土地开垦和土壤肥力管理。尽管如此，使用来源特异性的粪甾醇替代指标提供了一个实用框架，用于在粪便改良土壤中区分与粪便源协方差共享的污染物关联和预测因子特异性关联。从管理角度来看，减少畜牧生产和粪便处理过程中不必要的塑料输入，结合改进的粪便处理、土壤有机质管理及其他增强污染物滞留或降解的土壤基措施，可能有助于限制粪便伴生污染物的转移和环境持久性。未来的工作需要将更广泛的空间重复与受控的粪便-抗生素-微塑料实验相结合，以检验这些田间模式所暗示的因果机制。
