石油烃污染地下水羽流中微生物群落的构建：功能组织优于空间隔离

《Frontiers in Microbiology》：Functional organization outweighs spatial separation in structuring microbial communities across a petroleum-hydrocarbon groundwater plume

　　

摘要 目的： 石油烃污染羽流（petroleum-hydrocarbon plume）受水文地质运移、氧化还原条件及原位微生物活动共同塑造。然而，地下水微生物组之间的空间差异究竟主要反映井间物理隔离，还是反映羽流内与不同衰减状态相关的生态重组，目前仍不明确。 方法： 研究人员整合了某石油烃污染场地的水文地质背景、石油烃浓度、水化学变量及来自15口地下水监测井的鸟枪法宏基因组（shotgun metagenomic）图谱。评估了群落β多样性、空间与功能距离关系、基于置换的组间检验，以及自然衰减潜力与分类-功能多样性的探索性复合指数。 结果： 属水平群落相异度随估算空间距离的增加而增大（Mantel r = 0.240，p = 0.034），但与功能模块相异度的关联更强（r = 0.670，p = 0.001）。在控制空间距离后，生物降解潜力与群落相异度仍保持关联（偏Mantel r = 0.428，p = 0.001）。运行性衰减状态解释了属水平组成变异的35.4%（PERMANOVA p = 0.008）。高自然衰减潜力与高分类-功能多样性出现在部分不同的井位中。 结论： 结果表明，羽流内地下水微生物组的更替具有空间结构，但与功能重组的耦合程度高于与单纯距离的关联。联合考虑水文地质背景、微生物降解潜力与功能多样性，有助于改进地下水自然衰减评估。

石油烃污染地下水羽流中微生物群落构建与自然衰减潜力：基于宏基因组学的多维度整合研究解读

一、研究背景与科学问题

石油烃释放对地下水质量构成持续威胁，溶解态芳香烃及残余石油组分可在含水层中迁移，并同时经历稀释、吸附、相分配及生物降解等过程。其环境归趋不仅受污染物负荷控制，还受地下水流动、含水层结构、氧化还原条件及微生物代谢的共同影响。因此，监测自然衰减（monitored natural attenuation）成为石油烃污染地下水的重要管理选项。然而，可靠的衰减评估需要证明浓度变化反映的是衰减过程而非单纯的稀释或混合。

微生物转化是含水层中石油烃衰减的核心环节。烃类输入可改变电子受体可利用性，形成空间分异的氧化还原格局，进而塑造好氧降解菌、硝酸盐还原菌、铁还原菌、硫酸盐还原菌、发酵菌及产甲烷菌等不同类群的组合。但浓度数据本身不足以区分某口井代表高污染物压力、活跃响应的微生物群落，还是具有遗留功能容量的残余状态。

另一个挑战涉及地下水微生物组的空间组织方式。地下水微生物群落受含水层岩性、流场、补给、理化选择及扩散限制的共同影响。当前仍不明确石油烃污染羽流中地下水微生物群落的差异主要反映物理隔离与水文连通性，还是与烃类转化及衰减状态相关的生态重组。鸟枪法宏基因组学可同时表征分类组成与编码功能潜力，为弥合这一认知缺口提供了契机。

二、研究设计与核心发现

研究人员整合了位于中国山东省曲阜市某油库场地15口监测井（W1–W15）的水文地质背景、石油烃浓度、水化学变量及鸟枪法宏基因组图谱，系统评估了属水平微生物群落更替与空间隔离、功能组织及生物降解潜力之间的关系，并构建了运行性衰减状态框架以综合刻画污染物压力、生物响应与宏基因组功能信息。

研究得出以下核心结论：属水平群落相异度随估算空间距离增加而显著增大，但功能模块距离与群落相异度的关联远强于空间距离；生物降解潜力距离在控制空间距离后仍与群落相异度保持显著关联；运行性衰减状态解释了属水平组成变异的35.4%，高于污染分组（29.5%）的解释度；高自然衰减潜力井（W1、W2、W3、W7）与高分类-功能多样性井（W13、W15、W12、W9、W14）在空间上仅部分重合。

三、主要关键技术方法

研究采用以下关键技术方法：场地水文地质调查与地下水位插值；石油烃及水化学指标测定，包括苯系物（BTEX）的吹扫捕集-气相色谱-质谱联用分析及总石油烃（TPH C10-C40）的气相色谱分析；鸟枪法宏基因组测序，基于Illumina NovaSeq平台完成；生物信息学分析流程包括Kraken2/Bracken分类学注释、MEGAHIT拼接、Prokka基因预测、MMseqs2非冗余基因目录构建及KofamScan功能注释；统计分析方法包括Mantel检验、偏Mantel检验、置换多元方差分析（PERMANOVA）、Spearman置换相关分析及Benjamini-Hochberg错误发现率校正。

四、研究结果

4.1 水文地质背景与烃类空间异质性

15口监测井的地下水位高程介于61.890至62.382 m之间，插值水位面显示局部高值区叠加于整体较小差异之上，表明监测网络尺度上存在水文地质异质性。总BTEX浓度介于2.77至31.21 mg·L^?1，最高值出现在W3；TPH C10-C40浓度介于1.00至5.88 mg·L^?1，最高值出现在W2。BTEX与TPH C10-C40的空间分布并不完全重合。

4.2 群落更替的空间与功能关联

属水平Bray-Curtis相异度随估算井间距离增加而显著增大（Mantel r = 0.240，p = 0.034），表明群落更替具有空间结构。功能模块距离与群落相异度的关联最强（Mantel r = 0.670，p = 0.001），显著超过空间距离的关联。生物降解潜力距离亦与群落相异度显著关联（r = 0.455，p = 0.003），且在控制空间距离后仍保持显著（偏Mantel r = 0.428，p = 0.001）。环境距离与群落相异度的关联不显著。pH与聚合降解通路丰度（ρ = 0.760）及降解相关功能多样性（ρ = 0.687）呈显著正相关，乙苯、萘及氟苯甲酸降解通路与pH的正相关在校正后仍显著。

4.3 运行性衰减状态下的群落变异

15口井被划分为四种运行性衰减状态：高压力下活跃衰减（n = 5）、过渡衰减带（n = 5）、残余/恢复伴遗留响应（n = 3）及低影响/背景（n = 2）。PERMANOVA确认衰减状态与属水平群落组成显著关联（R² = 0.354，p = 0.008），为所有测试分类框架中解释度最高者。活跃衰减组中Acidovorax平均相对丰度最高（30.4%），过渡带中Sulfurospirillum相对富集（9.8%），低影响/背景组中Rhodoferax相对丰度较高（10.2%）。

4.4 衰减潜力与分类-功能多样性的井间分布

W7、W3、W2、W1在所有替代权重方案中均位列NAPI前五（top-five frequency = 1.00），表明其高衰减潜力判定对权重方案不敏感。分类-功能多样性最高值出现在W13（0.920）和W15（0.766），其后为W9、W12、W14。高NAPI井与高分类-功能多样性井在监测网络中占据部分不同的空间位置，二者并非同一组井。

4.5 功能模块图谱与分类群-模块关联

芳香烃响应模块在活跃衰减组中平均z标准化得分最高（0.978），羽流边缘氧化还原过渡模块在过渡组（0.205）及残余/恢复组（0.243）中较高，PAH/萘残余模块在残余/恢复组中得分最高（1.047）。Acidovorax与芳香烃响应模块呈名义显著正相关（ρ = 0.746，p = 0.002），但所有40项属-模块关联在校正后均不显著。

五、讨论与结论

讨论部分指出，空间距离对群落更替构成约束但非充分解释；功能组织与群落变异的关联远强于空间隔离。pH与部分降解通路丰度的正相关应视为位点特异性共变而非直接地球化学控制的证据。NAPI与分类-功能多样性刻画了自然衰减的不同维度：前者强调当代压力-响应关系的强度，后者反映降解相关分类与功能特征的丰富度和均匀度。监测策略应整合水文地质位置、污染物化学及宏基因组微生物指标，同时纳入高响应位点与高功能多样性位点。

研究结论为：属水平群落更替与空间距离存在显著但适度的关联，表明水文地质对微生物分布构成约束；但群落相异度与功能模块组织及生物降解潜力的关联更强，且在控制空间距离后仍成立。运行性衰减状态捕捉了群落组成与功能模块图谱的显著变异。高自然衰减潜力在高压力井子集中被稳健识别，而高分类-功能多样性主要出现在不同的井位中，表明当前衰减潜力与降解相关多样性是互补而非可互换的指标。研究基于单次采样，描述的是空间共变而非时间动态；宏基因组图谱代表编码功能潜力而非实际原位代谢活性。在此解释边界内，整合水文地质、化学、分类与功能证据可改进石油烃污染地下水中异质性衰减状态的识别与解释。