红树林沉积物中微生物对PBDEs-微塑料复合污染的多聚合物依赖性响应

《Ecotoxicology and Environmental Safety》：Polymer-dependent microbial responses to PBDEs-microplastic co-contamination in mangrove sediments

　　

红树林是介导碳、氮、硫循环的生物地球化学热点，但正日益受到持久性有机污染物与微塑料共现的威胁。尽管多溴二苯醚(PBDEs)与微塑料(MPs)的单独影响已有记载，不同微塑料聚合物如何在PBDE污染下调节微生物响应仍知之甚少。研究人员开展了为期六个月的微宇宙实验，探究红树林沉积物暴露于PBDEs单独处理及PBDEs联合聚丙烯(PP)、聚乙烯(PE)、聚对苯二甲酸乙二醇酯(PET)和聚苯乙烯(PS)时的微生物群落响应。在PBDE污染下，微塑料添加改变了沉积物理化条件，表现为氧化还原电位(Eh)降低及无机氮浓度的聚合物特异性变化。PBDEs显著重塑微生物群落组成(PERMANOVA, R2 = 0.685, p < 0.001)，降低变形菌门(Proteobacteria)相对丰度同时增加绿弯菌门(Chloroflexi)，但未显著影响α多样性。微塑料添加进一步诱导聚合物依赖性偏移。PE和PP提高香农多样性并富集氮循环类群（包括Rhodanobacter和Paenibacillus），伴随硝化与反硝化相关功能预测丰度升高。相反，PS和PET诱导 distinct beta多样性偏移并富集Mariprofundus、Acinetobacter和Pirellula等类群。预测功能谱显示PS/PET暴露下芳香烃降解、硫酸盐呼吸和C1化合物代谢相关功能丰度较高，提示向特化碳利用和硫相关途径的潜在转变。网络分析表明微塑料暴露降低网络复杂性，边数由7736降至2412–3168，模块度由0.26升至0.43–0.49，表明共现模式更加区室化。总体而言，PBDEs作为群落偏移的主要驱动因子，而微塑料以聚合物依赖方式进一步调节微生物组成、网络结构与预测功能潜力。

研究背景与立项依据：红树林位于热带亚热带海岸潮间带，是生产力和生物地球化学活性极高的森林生态系统，在碳固存、营养盐循环及海岸防护中具有关键作用。因其处于陆海界面，红树林沉积物易同时接收来自陆源与海源的污染物，且发达根系可增强污染物滞留，使持久性有机污染物与微塑料在沉积层中长期累积。多溴二苯醚(PBDEs)作为常用溴系阻燃剂，具有高环境持久性与强沉积物亲和性；微塑料(MPs)则通过改变孔隙度、氧渗透及提供定殖界面参与沉积物过程。现场调查已证实PBDEs与MPs在河口及红树林共存，但既往研究多聚焦单一污染物或单一聚合物类型，不同微塑料聚合物在PBDE背景胁迫下如何调节微生物群落、共现网络及功能潜力尚不清楚。研究人员假设含芳环的微塑料（PET、PS）因更易风化，比聚烯烃类（PE、PP）引发更强的微生物与功能偏移，据此在《Ecotoxicology and Environmental Safety》报道了该项系统比较研究。

主要技术方法概览：研究人员自深圳福田红树林自然保护区采集表层沉积物（0–5 cm），均质化后建立6组微宇宙（对照、PBDEs单独、PBDEs分别联合PP/PE/PET/PS），每处理三重复，模拟半日潮盐度15 psu、25℃培养六个月。沉积物理化测定涵盖pH、氧化还原电位(Eh)、总有机碳(TOC)、总凯氏氮(TKN)、铵态氮(NH₄⁺)与硝态氮(NO₃^?)。微生物分析采用16S rDNA V3–V4区扩增子测序（Illumina HiSeq 2500），按97%相似度聚类OTU并比对Greengenes数据库；α多样性、基于Bray–Curtis的PCoA与PERMANOVA、LEfSe差异类群、Mantel检验、Spearman共现网络（|r|≥0.7，BH校正p<0.05）及FAPROTAX功能预测构成核心统计框架；网络按无MP、脂族MP（PP/PE）、芳族MP（PET/PS）归组。

研究结果

3.1 沉积物理化性质变化：六个月末，pH与TKN组间无显著差异；Eh在PBDEs+PE/PET/PS中较对照显著下降。PBDEs+PP与PBDEs+PE的NH₄⁺显著低于对照与PBDEs单独处理，PBDEs+PS的NO₃^?显著最低；PBDEs+PET与PBDEs+PS的TOC显著高于对照与PBDEs单独处理。表明微塑料以聚合物特异方式改变化学计量与还原条件。

3.2 微生物群落组成：优势门为变形菌门、绿弯菌门、浮霉菌门(Planctomycetes)，合计超60%。PBDEs降低变形菌门、升高绿弯菌门；MPs添加未在门水平超越PBDEs效应，但在属水平呈聚合物依赖分化。LEfSe显示对照富集Clostridium等13属；PBDEs单独富集Candidatus Solibacter等7属；PBDEs+PS富集最多（16属，含Sphingobacterium、Mariprofundus、Pirellula等）。

3.3 多样性与群落结构：α多样性在PBDEs单独与对照间无差异，但PBDEs+PE与PBDEs+PP的Shannon、Sobs等指数显著升高。PCoA-PERMANOVA整体R²=0.685（p<0.001）；脂族MP较PBDEs单独未达显著额外偏移（p=0.052/0.066），芳族MP（PET p=0.040，PS p=0.012）则显著再偏移。Mantel检验表明群落与TOC（r=0.47）、Eh（r=0.35）显著关联。

3.4 共现网络响应：无MP网络含161节点7736边、模块度0.26；脂族MP网络边降至3168、模块度0.43；芳族MP网络边2412、模块度0.49，平均度由96.10降至37.49/29.78。芳族MP网络连接器节点最多（27个），整体鲁棒性相近（~0.50）。表明微塑料使共现结构简化且区室化，芳族聚合物关联更强模块分隔。

3.5 预测功能潜力：无MP组硫氧化与甲烷氧化相关功能较高；脂族MP组硝化、好氧氨氧化、反硝化等氮循环功能上调；芳族MP组芳香烃降解、硫酸盐呼吸、C₁代谢（甲基营养、甲醇氧化）功能富集。预测功能与实测理化总体Mantel不显著（r=0.13, p=0.26），但与TKN有弱正关联。

讨论浓缩：PBDEs是群落重构主驱因子，通过氧化应激与脱溴适应筛选耐污类群（如绿弯菌门），不改变α丰富度但驱动β置换。脂族MP因C–C/C–H骨架稳定，仅轻微调制；芳族MP含苯环/酯基，易光氧化释出复杂碳，诱发附加β偏移与TOC累积。网络简化与模块度上升反映微塑料作为离散基质促成微生境分区，连接器增多利于模块间维持关联，芳族MP区室化最强。功能预测显示脂族MP偏向氮循环电子流，芳族MP偏向特化碳利用与硫酸盐呼吸，与NO₃^?降低及还原条件一致；但FAPROTAX仅为潜力推断，需功能基因或速率验证。研究局限含高浓度暴露场景、缺MPs单独处理及未测老化表面化学。

结论翻译：本研究证明PBDEs-微塑料共污染以聚合物依赖方式改变红树林沉积物微生物群落。PBDEs作为主要驱动因子塑造群落组成，富集潜在胁迫适应与污染物转化类群，对α多样性影响有限。微塑料依聚合物类型进一步调节响应：脂族微塑料（PP、PE）关联氮循环预测功能富集，芳族微塑料（PET、PS）与群落变更及芳香化合物降解、硫酸盐呼吸功能更强关联。微塑料暴露还降低网络连通性并提高模块度。综上，PBDEs–微塑料共污染下微生物响应强烈依赖微塑料聚合物类型。所述功能解读基于FAPROTAX预测潜力而非直接过程测量，机制阐明仍需验证。