农业土壤微宇宙中聚乙烯与纸基覆盖物残留物对微生物响应的对比研究

目的（Aims）：农用覆盖材料可在土壤中产生颗粒态残留物。聚乙烯（Polyethylene, PE）覆盖物可能贡献持久性的微塑料（Microplastics, MPs）碎片，而纸基覆盖物则产生更易降解的残留物，其对微生物的影响尚不清楚。微生物对这些残留物的响应难以与培养期间自然的时间驱动群落演替区分。本研究通过在农业土壤微宇宙中比较聚乙烯微塑料与纸基微碎片，考察相对于时间匹配未处理对照组，颗粒类型和浓度如何影响细菌及真菌群落。 方法（Methods）：将聚乙烯微塑料和纸基微碎片以低、高浓度添加到取自葡萄牙Fataca地区农业土壤的无植物土壤微宇宙中。各采样时间点均设置未处理对照以分离改良剂相关效应与孵育驱动微生物演替。培养2个月和6个月后采用扩增子测序对细菌16S rRNA基因和真菌内部转录间隔区（Internal Transcribed Spacer, ITS）群落进行表征。每种孵育处理设6个独立生物学重复微宇宙。 结果（Results）：细菌群落主要受孵育时间塑造。未处理对照显示2个月时丰富度和多样性显著下降，6个月时恢复，表明存在强时间基线。相对于时间匹配对照，聚乙烯微塑料和纸基微碎片引起较轻微的细菌变化，包括高浓度下2个月时酸杆菌门（Acidobacteriota）增加、聚乙烯微塑料处理中慢生根瘤菌属（Bradyrhizobium）相对丰度升高，以及改良土壤6个月时芽孢杆菌纲/芽胞杆菌门（Bacilli/Bacillota）减少。真菌群落显示出更明确的改良剂相关响应。高剂量纸基微碎片使观测真菌丰富度较2个月对照降低约22%，较6个月对照降低约49%。该处理还促进高优势度真菌群落，早期壶菌门（Chytridiomycota）富集，后期子囊菌门（Ascomycota）特别是粪壳菌纲（Sordariomycetes）、球毛壳霉科（Chaetomiaceae）和Parachaetomium占主导。 结论（Conclusions）：在无植物微宇宙条件下，细菌响应受时间演替主导，而真菌群落表现出更强的改良剂相关响应，尤其在高量纸基微碎片添加时。结果表明塑料与纸基农用覆盖物残留可以不同方式影响土壤微生物轨迹，评估覆盖材料残留时应考虑颗粒类型、浓度及孵育时间。

论文解读——《Contrasting microbial responses to polyethylene and paper-based cover residues in agricultural soil microcosms》

该研究发表于《Plant and Soil》。

研究背景与立项依据

微塑料（Microplastics, MPs）因农用覆膜、温室覆盖、灌溉组件及有机肥施用已成为陆地生态系统尤其是农业土壤中的重要污染物。聚乙烯（Polyethylene, PE）覆膜在使用—风化—破碎—清除不全的过程中会在土壤中残留持久性塑料微粒，而纸基覆盖材料作为替代方案虽可减少持久性塑料输入，其降解产生的纤维素类微碎片（paper-based microfragments）作为可降解颗粒及潜在碳源进入土壤后的微生物学效应尚不明确。已有报道中土壤微生物对MPs的响应存在较大变异，且培养实验中微生物群落本身存在时间驱动的演替（temporal succession），若不设时间匹配对照易将背景演替误判为处理效应。因此，有必要在含时间匹配未处理对照的微宇宙体系中，比较PE-MPs与纸基微碎片对土壤原核（细菌）和真核（真菌）微生物群落α/β多样性及分类组成的差异化影响。

研究人员以葡萄牙Fataca地区农业土壤建立无植物土壤微宇宙，设低（0.001% w/w）、高（0.05% w/w）PE-MPs及纸基微碎片添加组，并设同期未处理对照，于孵育2个月和6个月采集样本，通过细菌16S rRNA基因V₃–V₄区及真菌ITS1区扩增子测序分析微生物群落变化，并以时间匹配对照为基准评估改良剂关联效应。

主要技术方法概述

研究材料为葡萄牙Fataca地区农田表土建无植物土壤微宇宙，添加研磨商用牛皮纸制备纸基微碎片及研磨40 μm厚PE膜制备聚乙烯微塑料，设低(0.001% w/w)、高(0.05% w/w)两浓度及相应未处理对照，初始(t=0)、孵育2月(t=2)、6月(t=6)采样，各处理6个生物学重复，维持60%田间持水量温室孵育。总DNA用DNeasy PowerSoil Pro Kit提取，细菌16S rRNA基因V₃–V₄区引物341F/806R、真菌ITS1区引物ITS5-1737F/ITS2-2043R扩增，Illumina双端测序。原始数据按样本拆分后FLASH合并、fastp质控、VSEARCH去嵌合体，QIIME2(v2022.02)中DADA2推断扩增子序列变体(Amplicon Sequence Variant, ASV)，细菌比对SILVA(v138.2)、真菌比对UNITE数据库分类注释。特征表按16S(46317 reads)和ITS(12432 reads)稀疏化(rarefy)后计算α多样性(Observed features/Chao1/Shannon/Simpson/Pielou evenness/dominance)及β多样性(加权/非加权UniFrac, PCoA, ANOSIM/PERMANOVA)。细菌功能预测用PICRUSt2(KO/EC/COG/PFAM)，真菌生态功能位(Guild)用FUNGuild注释，并做归一化随机率比(Normalized Stochasticity Ratio, NST, 含系统发育pNST与分类学tNST)探索性分析。原始数据存NCBI BioProject PRJNA1420942。

研究结果

细菌群落（Bacterial community）

—时间动态主导细菌多样性，为处理效应设定移动基线（Temporal dynamics dominate bacterial diversity, setting a moving baseline for treatment effects）：初始土细菌Observed features=3642、Shannon=10.32；2月未处理对照降至2308和8.40；6月未处理对照恢复至3509和10.09。证明须以同时间未处理对照为参照解读添加效应。

—细菌整体群落结构显示有限的处理分离（Bacterial whole-community structure shows limited treatment separation）：加权UniFrac分析未显示处理组间稳健分离，细菌群落主要由孵育时间而非颗粒添加结构化，添加效应表现为谱系水平（lineage-level）偏离而非整体重组。

—分类谱显示时间驱动演替伴选择性改良剂关联偏离（Taxonomic profiles show time-driven turnover with selective amendment-associated deviations）：未处理对照2月蓝藻菌门（Cyanobacteriota，主要为未分类叶绿体unclassified_Chloroplast）升至23.32%、酸杆菌门（Acidobacteriota）降至2.45%；6月Cyanobacteriota回降至6.61%、Acidobacteriota恢复至13.66%。相对于2月对照，所有改良处理均降低unclassified_Chloroplast相对丰度；高剂量处理(PE-high/Paper-high)使Acidobacteriae由2.17%分别升至6.12%/6.93%；PE处理中Bradyrhizobium轻微升高；6月所有改良土壤Bacilli/Bacillota较对照(6.91%→约1.88%–3.05%)降低。

—预测细菌功能谱及NST为探索性描述（Predicted bacterial functional profiles vary across time and amendments / Exploratory NST analyses suggest group-level differences in bacterial community turnover structure）：PICRUSt2推断功能潜力随时间及处理变化但与直接功能测定不等价；NST组间有差异(pNST:F=5.21,p=0.00103;tNST:F=15.82,p=2.90×10^-7)仅作周转结构补充描述。

真菌群落（Fungal community）

—真菌多样性具强时间信号，高纸基微碎片下持续抑制（Fungal diversity shows a strong time signal, with persistent suppression under high paper-based microfragments）：初始Observed features=1170、Shannon=6.561；2月对照降至594和5.804；6月对照恢复至1170和7.145。高纸基微碎片组2月Observed features=462(较2月对照↓≈22%)、Shannon=4.537；6月Observed features=598(较6月对照↓≈49%)、Shannon=4.014，优势度(dominance)升高，指示简化、少数类群主导。

—真菌群落结构较细菌更具处理敏感性（Fungal community structure is more treatment-sensitive than bacterial community structure）：UniFrac分析显示加纸基微碎片(尤高剂量)及部分PE处理在2月即与同期对照分离，6月高纸基微碎片组仍明显偏离。

—分类谱揭示独特改良剂关联真菌轨迹，高纸基微碎片驱动最强偏离（Taxonomic profiles reveal distinct amendment-linked fungal trajectories, with high micropaper driving the strongest departures）：高纸基微碎片2月壶菌门(Chytridiomycota)升至43.47%(主要为Rhizophlyctidomycetes/Rhizophlyctidales incertae sedis)，6月子囊菌门(Ascomycota)达77.49%且粪壳菌纲(Sordariomycetes)68.67%、球毛壳霉科(Chaetomiaceae)63.21%、Parachaetomium属升高，表明持续多样性抑制伴随窄谱类群优势。

—FUNGuild以未赋值为主（FUNGuild outputs were dominated by unassigned annotations）：高比例未赋值reads限制营养模式解读，高纸基微碎片6月未赋值达86.71%，仅作探索性生态背景参考。

讨论与结论总结

讨论指出孵育时间强烈塑造土壤微生物尤其细菌群落，不设时间匹配对照会高估添加效应；细菌对PE-MPs及纸基微碎片响应温和（谱系水平偏移如Acidobacteriota↑、Bacilli↓），不应过度解读为功能改变；真菌对高纸基微碎片响应显著——纸基微碎片兼具物理颗粒与可降解木质纤维素底体双重属性，可能选择性促进能利用分解残体或耐受微生境改变的狭谱腐生真菌（Sordariomycetes/Chaetomiaceae/Parachaetomium），但该机制需理化参数验证。研究局限含无植物微宇宙、仅两时间点、缺颗粒形貌/FTIR/土壤理化时序测定，结论限于扩增子观测模式。

结论（Conclusions）翻译：

在无植物微宇宙条件下，细菌群落响应受时间演替主导，而真菌群落尤其在高剂量纸基微碎片添加时表现出更明显的改良剂关联响应。这些发现表明塑料与纸基农用覆盖材料残留可以不同方式影响土壤微生物轨迹，评估覆盖材料残留时应考虑颗粒类型、浓度及孵育时间。