高寒湿地植被演替对nirK型反硝化菌群落组装与功能潜力的影响机制

《Frontiers in Microbiology》：Mechanisms underlying the effects of alpine wetland vegetation succession on the community assembly and functional potential of nirK-type denitrifier

　　

高寒湿地中的反硝化作用（denitrification）严格调控着区域氮收支与N₂O排放，nirK型反硝化菌是参与反硝化过程中亚硝酸盐还原的重要功能群。不同植物群落间的空间异质性可调节土壤微环境，从而塑造N₂O的产生、消耗与排放动态。本研究采用空间代替时间（space-for-time substitution, SFT）方法，调查了青海湖高寒湿地的四个演替阶段——早期（BS）、过渡期（CR）、中期（EN）和成熟期（LS）。结合高通量测序（high-throughput sequencing）与零模型（null model）框架，研究人员表征了微生物群落模式与其功能潜力之间的耦合关系。植被演替通过总氮（total nitrogen, TN）、总碳（total carbon, TC）和碳氮比（carbon-to-nitrogen ratio, CNR）的变化刺激了微生境分化。这种环境过滤（environmental filtering）诱导了定向微生物演替，最终在LS阶段达到最高多样性。CNR成为解释群落方差的主要驱动因子，介导了微生境生态位分离，选择性富集了关键异养反硝化菌（如中华根瘤菌属Sinorhizobium和副球菌属Paracoccus）。βNTI（β-nearest taxon index）零模型揭示确定性过程主导群落组装；值得注意的是，随着演替推进，潜在的组装规则发生根本性转变，在LS阶段以异质性选择（heterogeneous selection）为主导。共生网络（co-occurrence network）在EN阶段复杂性达到峰值，但在LS阶段重组为更加整合、模块化且高度竞争的拓扑结构。此外，功能分析表明，BS阶段在功能上偏向硝酸盐呼吸和固氮作用，而LS阶段在预测的反硝化相关功能谱方面表现出显著优势。总之，植被演替通过环境过滤重塑了nirK型反硝化菌群落，将组装机制从随机性转向确定性；同时，多样化群落结构的"保险效应"（insurance effect）在养分波动中维持了稳定的功能输出。本研究阐明了植被演替是nirK型反硝化菌亚硝酸盐还原潜力的核心驱动力，为气候变化下退化高寒湿地恢复、植被重建及区域氮稳态维持提供了微观尺度机制见解。

**高寒湿地nirK型反硝化菌群落对植被演替的响应机制**

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

一氧化二氮（N₂O）是一种强效非CO₂温室气体，其增温潜势约为二氧化碳（CO₂）的273倍，且可在平流层分解释放一氧化氮，催化破坏臭氧层。高寒湿地因长期低温抑制微生物分解而保有巨大土壤氮库，是N₂O排放的重要潜在源。反硝化作用（NO₃^?→N₂O→N₂）是调控气态氮损失与净氮去除的核心途径，其中亚硝酸盐还原（NO₂^?→NO）由nirK或nirS基因编码的亚硝酸盐还原酶催化，构成整条反硝化通路的限速步骤。青海湖盆地高寒湿地正经历气候变暖与水位波动驱动的强烈植被演替，然而关于植被演替如何驱动反硝化菌群落及其气态氮排放，目前存在冲突性范式。该研究旨在阐明青海湖湿地植被演替序列对nirK型反硝化菌群落结构更替与组装机制的影响，并探究其氮去除功能如何通过"保险效应"（insurance effect）得以维持，以加深对高寒内流湿地氮去除机制的微观尺度理解。

**二、研究内容与总体结论**

研究人员采用空间代替时间（space-for-time substitution, SFT）方法，选取青海湖流域四种代表性演替阶段——早期（Blysmus sinocompressus群落, BS）、过渡期（Carex regescens群落, CR）、中期（Elymus nutans群落, EN）和成熟期（Leymus secalinus群落, LS），共10个典型样地、50个独立样本（BS=10, CR=10, EN=15, LS=15）。研究揭示：植被演替通过环境过滤重塑了nirK型反硝化菌群落，驱动组装机制从随机性向确定性（异质性选择）根本转变，同时多样化群落的"保险效应"在养分波动中维持了稳定的功能输出。相关研究发表于《Frontiers in Microbiology》。

**三、关键技术方法**

研究依托青海湖流域高寒湿地自然演替序列土壤样本，采用Illumina NovaSeq 6000平台对nirK功能基因（引物copper583F/copper909R）进行高通量测序；基于97%相似度聚类获得操作分类单元（OTU）；利用LEfSe（线性判别分析效应大小）筛选差异标志类群；通过BLASTX比对NCBI非冗余（nr）蛋白数据库进行属级分类注释，构建宿主-功能映射框架以预测代谢功能谱；采用βNTI零模型量化划分群落组装的确定性与随机性过程；构建共生网络分析拓扑属性；结合冗余分析（redundancy analysis, RDA）、Mantel检验及方差贡献分解识别核心环境驱动因子。

**四、主要研究结果**

**1. 土壤理化性质与植物生物量沿演替梯度的变化**：地下生物量（belowground biomass, BB）在CR阶段显著达峰，随后在EN和LS阶段显著下降；土壤总碳（TC）与总氮（TN）在BS、CR、EN阶段维持较高水平，在LS阶段显著降低；碳氮比（CNR）在BS与CR阶段显著高于EN与LS阶段；土壤pH沿演替梯度显著下降，至LS阶段呈现明显酸化趋势。

**2. nirK型反硝化菌α多样性及分类学特征**：稀释曲线表明测序深度充分。LS阶段的ACE、Chao1和Shannon指数均显著升高，特有OTU数在LS阶段激增至10,691个，表明成熟阶段物种丰富度与整体多样性结构显著提升。

**3. β多样性空间分异及多层级分类学响应**：非度量多维尺度分析（non-metric multidimensional scaling, NMDS）显示演替阶段间群落差异显著（R²=0.25, p<0.001）。变形菌门（Proteobacteria）和慢生根瘤菌属（Bradyrhizobium）为绝对优势类群；LS阶段显著富集中华根瘤菌属（Sinorhizobium）和副球菌属（Paracoccus），而假单胞菌属（Pseudomonas）、Nitratireductor等早期优势属显著下降，核心属发生剧烈重组。

**4. 共生网络及拓扑属性**：网络复杂性与总边数在EN阶段达峰值；至LS阶段总边数与平均度显著收缩，但节点数与模块度指数达整个演替序列最高值，负相关边比例持续阶梯式升高至LS阶段最高，反映群落向高度模块化、竞争性拓扑结构主动重组。

**5. 功能谱演化与"保险效应"**：CR阶段硝酸盐呼吸、硝酸盐还原及亚硝酸盐呼吸功能预测丰度相对增强；LS阶段总反硝化、亚硝酸盐反硝化与氧化亚氮反硝化功能预测丰度呈阶梯式显著上升，表明成熟阶段在核心反硝化功能谱上具备显著的功能补偿优势。

**6. 核心环境驱动因子与群落组装机制**：Mantel检验表明植物生物量与土壤碳氮养分是关键正向驱动因子；方差贡献分解显示CNR（+17.3%）和pH（+13.6%）是结构与功能差异的首要驱动因子。βNTI分析显示其值随演替呈连续阶梯式显著上升；BS阶段随机过程贡献较高，至LS阶段异质性选择成为绝对主导力量，而随机过程贡献大幅降低。

**五、讨论总结与研究结论**

讨论部分指出，成熟阶段"资源耗竭—多样性峰值"的非线性解耦格局，反映高寒湿地群落通过生态位扩张、网络模块化巩固和功能特化属定向富集，在养分低谷期实现了适应性缓冲。LS阶段网络连接度下降属于高度有序的"主动结构响应"，而非退化性瓦解。研究结论表明：植被演替与nirK型反硝化菌群落组装及预测功能潜力密切关联；成熟期植物定殖诱导的生态位扩张显著提升系统物种承载力；功能预测揭示代谢潜力呈阶段性更替，LS阶段在核心反硝化相关功能谱中表现出显著补偿优势；TC与TN正向驱动群落α多样性，而群落更替与代谢功能受CNR、局部pH（酸化趋势）及BB协同调控；随演替推进，群落组装机制从随机漂变主导彻底转向确定性"异质性选择"控制，赋予系统在养分低谷期独特的适应性缓冲弹性，为退化高寒湿地恢复、植被重建及区域氮稳态维持提供了微观尺度机制依据。