关键K策略菌塑造四川盆地页岩气井场土壤微生物共现网络的生态机制

《Biology》：Keystone K-Strategists Shape Microbial Co-Occurrence Networks in Shale-Gas Wellpad Soils in the Field, Sichuan Basin, China

　　

页岩气开采通过返排水盐渍化、钻屑沉积及中度石油烃污染对陆地土壤生态系统施加持续、低强度、多因子扰动。然而，此类复合慢性扰动下微生物群落的组装机制及其相应生态策略仍认识不足。本研究采用16S rRNA基因测序、定量PCR（quantitative PCR, qPCR）和分子生态网络（molecular ecological network, MEN）分析，对中国四川盆地长宁页岩气田的土壤微生物组进行了表征。结果表明，页岩气井场土壤微生物群落形成了高度模块化（M = 0.786）、非随机的“小世界”共现网络，指示其具有内在稳定的微生物系统。通过Zi-Pi拓扑分类，研究人员鉴定出6个低丰度（<0.63%）的关键物种：3个模块枢纽（OTU11572、OTU14908、OTU11976）和3个连接者（OTU8845、OTU13035、OTU8962），系统发育归属于酸杆菌门（Acidobacteriota）、绿弯菌门（Chloroflexi）和放线菌门（Actinobacteriota）。这些关键微生物表现出典型的K策略特征，包括高网络连通度（度=21–24）、代谢稳定性及对种间协作互作的依赖。研究人员检测到这些K策略关键物种与负责反硝化和硫酸盐还原的厌氧功能基因（nirS、nirK、dsrAB）之间存在显著正相关，凸显其对钻屑堆积所产生的波动性氧化还原环境的适应能力。基于距离的冗余分析（distance-based redundancy analysis, Db-RDA）进一步证明，电导率（代表返排水侵入产生的盐渍化胁迫）和土壤有机质是塑造微生物群落组成的主导环境过滤器，其作用超过易利用的无机养分。总体而言，本研究为页岩气开采相关的慢性土壤扰动下K选择主导微生物群落组装提供了新证据。关键K策略菌通过代谢可塑性与厌氧呼吸适应的结合维持模块化网络稳定性。上述发现将r/K选择理论扩展到人为扰动的土壤生境，并为页岩气开采区域的土壤生态恢复和风险评估提供了理论指导。

**研究背景与科学问题**

页岩气开采已成为全球能源供应的重要组成部分，中国页岩气储量主要集中于四川盆地，其中长宁—威远国家级页岩气示范区自2009年以来已钻探超过2000口井，年产能超过146亿立方米。然而，页岩气开发通过钻井液泄漏、油基钻屑堆积和降水径流等途径对土壤生态系统引入复杂扰动。与急性溢油污染不同，这些扰动具有低浓度、持续性和多源复合特征，对土壤微生物群落的结构与功能构成持续累积压力。现有r/K选择理论框架多源于溢油扰动生态系统，关于慢性物理土地扰动下微生物网络拓扑、关键物种动态与r/K策略响应的实证证据仍然不足。分子生态网络（molecular ecological network, MEN）分析虽已应用于重金属污染土壤、农业生态系统和废水处理反应器，但尚未见诸于页岩气井场周边土壤。因此，本研究旨在填补这一空白，揭示慢性人为扰动下微生物群落的组装机制与稳定性维持策略，为页岩气开发区域的生态风险评估和土壤管理提供理论依据。

**研究内容与主要技术方法**

研究人员以四川盆地南部长宁—威远国家级页岩气示范区为研究区域，在珙县和兴文县选取3个运行超过两年的代表性井场，设置9个采样点，采集0–20 cm表层土壤。采用Illumina NovaSeq PE250平台对16S rRNA基因V3-V4高变区进行双端测序（2 × 250 bp），序列经QIIME 2（DADA2）降噪去嵌合后按97%相似度聚类为OTU，并基于SILVA 138数据库进行物种注释；厌氧功能基因（nirS、nirK、dsrAB）丰度采用ABI 7500实时荧光定量PCR系统测定；共现网络基于Spearman秩相关（|ρ| > 0.7, p < 0.01）通过随机矩阵理论（random matrix theory, RMT）框架构建，利用Cytoscape实现可视化；关键物种通过Zi-Pi拓扑分类鉴定；环境驱动因子采用基于距离的冗余分析（distance-based redundancy analysis, Db-RDA）和方差分解解析。

**3.1 微生物群落结构、多样性与分类特征**

高通量测序获得1,247,583条高质量细菌序列，聚类为12,456个OTU。NMDS排序显示采样点间微生物群落存在显著空间分异（ANOSIM: R2 = 0.326, p = 0.0399）。核心井场扰动区的Shannon多样性指数（5.23 ± 0.41）显著低于未扰动缓冲区（6.08 ± 0.35）（p < 0.05），而Simpson优势度指数呈相反趋势，表明长期扰动简化了群落结构并增强了耐逆类群的竞争优势。门水平优势菌门为Chloroflexi（26.0%）、Acidobacteriota（23.3%）、Proteobacteria（14.2%）和Actinobacteriota（11.7%），其中Chloroflexi和Acidobacteriota在核心扰动土壤中显著富集（p < 0.01），两者均为典型寡营养K策略菌门，其富集证实慢性扰动选择性淘汰了富营养r策略菌。属水平上，Sphingomonas（4.8%）、Mycobacterium（3.2%）和Alkanindiges（2.1%）为核心指示属，其共富集反映了井场土壤盐碱胁迫与慢性工业扰动的双重选择压力。Welch's t检验鉴定出4个显著差异OTU（OTU1204、OTU120、OTU322、OTU201），涉及硫代谢与耐盐功能。

**3.2 共现网络构建与关键物种鉴定**

MEN分析构建了包含892个节点和10,744条边的共现网络（|ρ| > 0.7, p < 0.01），平均度=24.0、平均路径长度=5.29、聚类系数=0.44、模块度（modularity, M）=0.786，显著高于随机网络（M=0.27），呈现非随机“小世界”拓扑特征；节点度分布拟合幂律模型（R2 = 0.77），体现“核心—外围”组织结构。Zi-Pi拓扑分类鉴定出6个关键物种：3个模块枢纽（OTU11572、OTU14908、OTU11976）和3个连接者（OTU8845、OTU13035、OTU8962），分属Acidobacteriota、Chloroflexi和Actinobacteriota。这些关键物种丰度极低（<0.63%，均值0.18%），但具有高连通度（度=21–24）、代谢稳定性、种间协作依赖和环境耐受性等典型K策略特征，支持“低丰度—高影响力”生态范式。

**3.3 K策略物种与厌氧适应性的关联**

qPCR定量显示，nirS、nirK和dsrAB基因拷贝数范围分别为2.1×10⁶至8.7×10⁷、1.5×10⁶至5.3×10⁷和8.2×10⁵至3.6×10⁷拷贝/克干土，各基因丰度在不同采样区间存在显著空间异质性（p < 0.05）。Mantel检验显示模块枢纽OTU11572（Acidobacteriota）与nirS丰度显著正相关（r = 0.43, p = 0.02），连接者OTU8845（Chloroflexi）与dsrAB丰度强正相关（ρ = 0.65），表明K策略核心类群通过反硝化和硫酸盐还原等厌氧呼吸代谢适应井场土壤频繁的氧化还原波动，厌氧功能通路的整合增强了其代谢稳定性和环境适应能力。

**3.4 群落与网络格局的环境驱动因子**

Db-RDA识别出电导率（electrical conductivity, EC，解释22.4%变异）、Ca²⁺（18.7%）、Mg²⁺（16.6%）和NO₃^?（15.3%）为塑造微生物群落结构的主要环境因子。方差分解显示Ca²⁺对共享细菌群落解释力最强（31.5%, p = 0.004），Mg²⁺主要解释位点特异性群落的变异（20.4%, p = 0.044）。关键物种分布与土壤有机质显著正相关（r = 0.55），而与NH₄⁺、NO₃^?等易利用养分相关性较弱；EC与K策略微生物丰度正相关（r = 0.58, p < 0.05），证实返排水侵入引起的盐渍化胁迫对耐盐K策略菌产生强定向选择，离子强度与碳源可用性共同驱动了胁迫耐受型紧密关联微生物网络的形成。

**3.5 理论意义与管理启示**

研究人员提出“扰动—稳定性”假说，将r/K选择理论扩展至人为扰动陆地生态系统。与急性溢油污染选择快速生长的r策略烃降解菌不同，页岩气井场土壤以持续离子盐渍化为主导环境过滤器，偏好稳定低质量有机碳资源，形成K策略主导的微生物共现网络，赋予群落更高的功能稳定性和抗离子胁迫能力，但突发碳输入下代谢恢复较慢。管理上，保护本土K策略关键类群、增加土壤有机质并缓解离子胁迫优于外源易利用养分添加。

**结论**

本研究得出四点核心结论：（1）井场土壤形成高度模块化（M = 0.786）“小世界”微生物网络，慢性低强度扰动可培育稳定且功能冗余的微生物群落；（2）6个低丰度（<0.63%）高连通度关键物种（OTU11572、OTU14908、OTU11976、OTU8845、OTU13035、OTU8962）被鉴定为K策略菌，具有代谢稳定、种间协作和厌氧适应特征；（3）K策略菌与厌氧功能基因（nirS、nirK、dsrAB）显著正相关，反硝化和硫酸盐还原代谢支持其对波动氧化还原条件的适应并维持群落功能完整性；（4）EC和土壤有机质是塑造K策略菌空间分布的主要环境过滤器。研究同时指出样本量有限和基于相关性的网络推断等局限性，建议未来结合宏基因组分箱、代谢重构、稳定同位素探针（stable isotope probing, SIP）和定向培养等手段验证K策略关键物种的代谢协作机制。上述发现为页岩气开采区域的原位土壤生态恢复与风险评估提供了微生物学理论支撑和决策参考。