青海湖高寒湿地优势植物群落间土壤nosZ型反硝化菌的群落分异与装配

《Frontiers in Microbiology》：Community differentiation and assembly of soil nosZ-type denitrifiers across dominant plant communities in an alpine wetland of Qinghai Lake

　　

背景：由携带nosZ的微生物介导的N₂O还原是反硝化的末端步骤，也是湿地土壤中N₂O的重要生物汇，但自然存在的优势植物群落如何与高寒湿地中nosZ型群落的多样性、分类分异和装配相关联仍知之甚少。方法：研究人员调查了青藏高原东北部青海湖周边湿地中由Blysmus sinocompressus（BS）、Carex regescens（CR）、Elymus nutans（EN）和Leymus secalinus（LS）主导的四个优势植物群落，并将靶向nosZ的扩增子测序与植被生物量、土壤理化性质、环境关联分析以及基于系统发育的零模型相结合，以评估群落多样性、分类组成、环境相关性和装配过程的变化。结果：四个植物群落表现出显著的植被—土壤条件差异，并支持不同的nosZ型群落。LS相关土壤显示出最高的丰富度和最大的群落特异性OTUs（operational taxonomic units，操作分类单元）库，而Proteobacteria（变形菌门）和Bradyrhizobium（慢生根瘤菌属）主导总体分类组成。环境关联在群落多样性和分类结构之间不同：多样性更强烈地与植物生物量、土壤温度以及碳氮状况相关，而组成变异更紧密地与土壤C:N比（CNR）、pH、总碳（TC）、水分和硝态氮（NO₃^?-N）相关。基于系统发育的零模型分析进一步表明，群落装配主要与变量选择（variable selection）相关，尽管扩散限制（dispersal limitation）和其他随机过程（stochastic processes）的相对贡献在优势植物群落间存在差异。结论：这些发现表明，青海湖高寒湿地中的优势植物群落与植被—土壤环境、nosZ型群落组成和装配过程的协同分异相关，从而表征了引物定义的土壤nosZ群落在对比植物群落生境间的分类分异和推断装配；然而，对实际N₂O还原的意义需要使用基因丰度、转录、酶活性和气体通量测量进行直接验证。

研究背景：高寒湿地水文与氧化还原梯度强烈，微生物氮转化耦合紧密。N₂O是长寿命温室气体，反硝化末端N₂O还原为N₂由nosZ编码的N₂O还原酶催化，是唯一已知的酶促生物汇；nosZ不仅存在于经典反硝化菌中，也存在于缺乏完整反硝化途径的多样微生物中，非反硝化N₂O还原菌可降低土壤N₂O产生。植被可通过根系生物量、根际沉积、凋落物和根际条件改变有效碳氮底物、土壤水分等，从而影响反硝化菌。然而，自然高寒湿地优势植物群落如何与nosZ型反硝化菌（携带nosZ基因、催化N₂O还原的微生物）群落的多样性、分类分异和装配相关联仍不清楚；现有证据多来自农业管理试验、控制植物—土壤系统或单因子操控，自然优势植物群落是否对应不同装配机制尚不明确。因此，研究人员开展本研究。研究在青海湖周边湿地选择BS、CR、EN、LS四个优势植物群落，整合nosZ靶向扩增子测序、植被生物量、土壤理化、环境关联分析和基于系统发育的零模型，评估群落多样性、分类组成、环境相关性和装配过程。由于植物群落类型与站点条件共变，该观察性设计用于识别生境关联而非植物身份的因果效应。结论表明，四个优势植物群落与显著不同的植被—土壤条件和nosZ型群落相关；LS相关土壤nosZ丰富度最高、独特OTUs（操作分类单元）池最大；Proteobacteria（变形菌门）和Bradyrhizobium（慢生根瘤菌属）主导；多样性与分类组成关联于部分不同的环境维度；群落装配以变量选择（variable selection）为主，扩散限制（dispersal limitation）和其他随机过程（stochastic processes）贡献不同。其意义在于整合植被—土壤异质性、分类分异和群落装配来理解高寒湿地nosZ型反硝化菌生态，但实际N₂O还原后果仍需基因丰度、转录、酶活和气体通量验证。

关键技术方法：样本队列来源于2020年8月青海湖周边湿地10个空间分离站点、50个样方，包括2个BS站点（10样方）、2个CR站点（10样方）、3个EN站点（15样方）、3个LS站点（15样方）；每站点5个1×1 m样方，采集0—10 cm表层土混合。测定地上生物量（AGB）、地下生物量（BGB）、pH、土壤温度、水分、总碳（TC）、总氮（TN）、铵态氮（NH₄⁺-N）、硝态氮（NO₃^?-N）和土壤C:N比（CNR）。提取土壤DNA，以nosZ-F/1622R引物扩增nosZ，Illumina HiSeq 2500测序；PANDAseq合并、PRINSEQ质控、USEARCH去嵌合、UCLUST按97%相似度聚类OTUs；QIIME 1.9.1和R 4.0.5计算α多样性，NMDS/PERMANOVA/PERMDISP/ANOSIM、PCA分析群落分异，相关和Mantel分析环境关联，βNTI（β最近分类群指数）和RCbray零模型推断装配过程。

研究结果：
3.1 植被与土壤性质在优势植物群落间不同：通过植被生物量和土壤理化比较，AGB在CR最高，EN和LS较低，BS居中；BGB以CR最高、BS次之，EN和LS显著较低；CNR和pH在BS、CR较高，CR pH最高；TC和TN在LS最低；NH₄⁺-N在LS低于CR和EN。表明四类优势植物群落存在明显植被—土壤环境异质性。
3.2 nosZ型反硝化菌多样性与OTU分布：稀释曲线随测序深度增加趋于饱和，LS样本丰富度估计较高。ACE和Chao1以LS最高、CR最低，LS显著高于其他三类；Shannon在EN和LS较高；Simpson变化较窄但部分群落对差异显著。共获得40,225个nosZ OTUs，仅859个为四类共享；BS、CR、EN、LS分别含2,411、2,816、6,787、16,644个独特OTUs，LS独特OTU池最大。
3.3 nosZ型反硝化菌群落分异与分类重组：PCA中PC1和PC2分别解释23.84%和13.15%，合计36.99%；NMDS拟合可接受（Stress=0.1147），Adonis R²=0.28、p<0.001，ANOSIM R=0.52、p<0.001，BS与LS分离明显。门水平上Proteobacteria主导，BS和CR较高，EN和LS降低，LS未分类类群增加，Planctomycetes较少；属水平上Bradyrhizobium主导但比例变化，Mesorhizobium、Herbaspirillum、Azospirillum、Paracoccus、Polaromonas、Shinella和未分类类群亦有贡献；Azospirillum、Bradyrhizobium、Herbaspirillum、Paracoccus、Polaromonas和Shinella的测序丰度存在显著组间差异。表明总体组成和多个优势类群相对表示不同。
3.4 nosZ型反硝化菌多样性与群落结构的环境关联：多样性变异中BGB贡献最大（21.1%），其次AGB（9.7%）、土壤温度（9.2%）、TC（8.7%）、TN（7.4%）；群落结构变异中CNR贡献最大（22.5%），其次pH（16.5%）、TC（13.4%）、水分（7.8%）、NO₃^?-N（7.6%）、TN（6.9%）、温度（5.4%）、BGB（5.1%）。Mantel分析显示门水平和属水平群落矩阵与部分植被—土壤变量显著相关，环境变量间存在明显共变。表明多样性和分类组成关联于不同环境维度。
3.5 nosZ型反硝化菌群落装配在优势植物群落间存在差异：βNTI在四类群落中均为正且显著超过+2；BS与CR无显著差异，EN和LS通常较高，LS显著高于EN。变量选择为主要装配成分，在CR和EN尤其突出；BS中变量选择与扩散限制贡献相近，CR、EN、LS中扩散限制较小；漂移在BS、CR、EN中较小，在LS中较大；同质选择和同质扩散贡献很小。表明变量选择主导，但选择、扩散相关过程和漂移的相对平衡在植物群落间不同。

讨论部分总结：4.1 四类优势植物群落的nosZ多样性和组成差异明显，LS丰富度和独特OTU池最高，PCA和NMDS显示群落分异；这与植物群落通过根系生物量、根际沉积、凋落物、养分获取和土壤理化塑造地下生态模板并选择不同微生物 assemblage 的观点一致，Proteobacteria和Bradyrhizobium主导也与既往nosZ调查一致；但植被与土壤共变，应解释为植物群落相关分异而非植物身份直接因果。4.2 多样性与分类组成关联于不同环境维度：BGB对多样性贡献最大，组成更受CNR、pH和TC影响；不同群落属性不必同步响应环境变化，根际可综合影响碳有效性、氧状态、养分供应和微位点，pH可在基因表达、酶成熟和N₂O还原动力学等多层面影响反硝化；但环境变量强共变，贡献值应视为统计重要性而非独立因果，BGB、CNR、pH和TC为候选环境过滤因子。4.3 βNTI一致为正且超过+2，说明系统发育周转强烈偏离零期望，在所用零模型框架内主要与变量选择相关；CR和EN的变量选择较强，BS中扩散限制贡献较大，LS中漂移比例较大，这不否定整体选择信号，而表明局部装配同时含确定性和随机成分；零模型类别是统计推断的残差过程而非直接观测。4.4 生态意义与局限：植物生物量、根源输入、凋落物和植被介导的土壤微环境可生成异质地下生境，多样性和组成代表不同微生物组织维度；由于植物群落位于不同站点，应解释为植物群落相关分异；零模型推断非直接测量。局限包括单生长季采样、站点级重复有限、季节动态和空间依赖不能完全解析；nosZ-F/1622R常规引物可能低估非Proteobacteria类群如Firmicutes和Bacteroidetes；DNA扩增子测序不量化绝对基因丰度、转录、酶活性和实际N₂O→N₂通量。未来需更广或分支解析的nosZ检测、qPCR或digital PCR、转录水平测定、N₂O和N₂通量直接测量、重复季节采样和更多站点重复。

研究结论翻译：青海湖高寒湿地不同优势植物群落与nosZ型反硝化菌群落在多样性、分类组成和装配上的显著分异相关。Leymus secalinus相关土壤支持较高的nosZ丰富度和较大的独特OTUs池，而Proteobacteria和Bradyrhizobium在各类植物群落nosZ群落中占主导。nosZ型反硝化菌的多样性和分类组成与部分不同的环境维度相关，表明群落多样性变化和组成分选可能对植被—土壤环境异质性作出不同响应。零模型分析进一步表明，群落装配主要与变量选择相关，而扩散限制和其他随机过程的推断贡献在优势植物群落生境间不同。总体而言，这些发现强调整合植被—土壤异质性、分类分异和群落装配以理解高寒湿地nosZ型反硝化菌生态的重要性，但实际N₂O还原后果仍需使用绝对丰度、基因表达和气体通量测量进一步验证。