西南喀斯特石漠化梯度对豌豆根瘤内生固氮细菌多样性的影响

《Horticulturae》：Endophytic and Diazotrophic Bacterial Diversity in Pisum sativum Root Nodules Across Southwest China’s Rocky Desertification Gradients Qiuli Yan, Chengyi Wu, Wuxian Zhang, Yating He and Jinhua Wang

本研究发现，喀斯特石漠化(RD)梯度通过改变土壤理化性质，影响了豌豆（Pisum sativum）根瘤内生及固氮细菌群落结构。其中，Allorhizobium–Neorhizobium–Pararhizobium–Rhizobium和Pseudomonas为优势属，Rhizobium是主要固氮菌。轻度石漠化(SRD)区域根瘤内生菌多样性更高，且碳水化合物代谢潜能预测更强。土壤pH、总磷(TP)、全氮(TN)和土壤含水量(SWC)是驱动群落变化的关键环境因子。该研究揭示了微生物对石漠化胁迫的响应机制，为喀斯特退化地区利用微生物辅助豆科植物进行生态恢复提供了科学依据。

　　

西南喀斯特石漠化梯度对豌豆根瘤内生及固氮细菌多样性的影响研究

1. 引言

豆科植物与固氮根瘤菌的共生关系是自然界最高效的生物固氮(BNF)系统之一。这种以形成根瘤为特征的共生关系，在根瘤内，根瘤菌能将大气中的氮转化为植物可利用的氨，对维持全球氮循环和可持续农业发展至关重要。然而，在我国西南喀斯特地区，石漠化(RD)作为一种严重的土地退化形式，导致植被丧失、土壤侵蚀加剧、基岩大面积裸露，严重威胁着生态安全和当地生计。

在传统作物难以生存的贫瘠土壤中，引入豆科植物进行生态修复与农业生产具有重要意义。豌豆（Pisum sativum）因其耐旱、根系发达、可进行共生固氮等特性，被视为一个理想的候选物种。然而，在石漠化胁迫下，根瘤这一关键的共生微环境及其内部的细菌群落如何响应，目前尚不清楚。本研究旨在探究中国西南喀斯特石漠化梯度下，豌豆根瘤内生细菌及固氮菌群的多样性与演替规律。我们假设：(i) 石漠化程度加剧会驱动豌豆根瘤微生物群落发生定向变化，其特征是固氮共生菌减少而耐胁迫类群增加；(ii) 这种演替由沿RD梯度的土壤环境因子介导。

2. 材料与方法

2.1. 研究地点

实验材料采自红河州不同RD程度（基岩裸露率30-50%为轻度石漠化(SRD)，50-70%为中度石漠化(MRD)）的区域，包括个旧市、建水县和开远市的四个采样点(S_GJ2, S_JS3, M_GJ7, M_KY1)，每个点设置3个生物学重复，共获得12个样本。

2.2. 样品采集与处理

采集具有良好根瘤的完整豌豆植株。将根系分离，根瘤经过表面灭菌后用于提取内生细菌总DNA。同时采集附着在根表的根际土壤（0-2 mm范围）用于测定土壤理化性质。

2.3. 土壤环境因子测定

测定了土壤pH、土壤含水量(SWC)、全氮(TN)、碱解氮(AN)、全磷(TP)、有效磷(AP)等理化指标。

2.4. DNA提取与测序

使用HiPure Soil DNA提取试剂盒提取豌豆根瘤内生细菌总DNA。以基因组DNA为模板，分别扩增16S rRNA基因的V5-V7区（用于分析总内生菌群落）和固氮酶nifH基因（用于分析固氮菌群落）。扩增产物在Illumina MiSeq平台上进行双端测序。

2.5. 生物信息学与统计分析

对原始测序数据进行质控、合并、去除嵌合体，在97%相似度水平下聚类为操作分类单元(OTU)。计算了Chao1、Shannon等α多样性指数，并基于Bray-Curtis距离进行主坐标分析(PCoA)。通过PICRUSt2预测群落的功能潜力。利用冗余分析(RDA)和Mantel检验等方法分析了微生物群落结构与土壤环境因子之间的关系。

3. 结果

3.1. 不同RD程度的土壤理化性质存在差异

土壤理化分析揭示了采样点间显著的空间异质性。中度石漠化点位M_GJ7的TN和AN含量显著高于其他点位。建水县的S_JS3点位则具有较高的pH、AP和SWC。相关分析显示，AP与TP、AN均呈显著负相关，表明在喀斯特土壤中，磷的有效性可能在高总磷和总氮条件下受到限制。

3.2. 豌豆根瘤内生菌群落结构与多样性对喀斯特RD程度的响应

通过对16S rRNA基因测序分析发现，豌豆根瘤内生菌的α多样性（Shannon指数）在轻度石漠化(SRD)点位显著高于中度石漠化(MRD)点位。该指数与土壤含水量(SWC)呈显著正相关，与RD程度呈显著负相关，表明根瘤内生菌多样性随土壤水分可利用性增加而增加，随石漠化加剧而下降。

β多样性分析（PCoA）显示，不同区域的样品群落结构存在明显差异。在属水平上，根瘤内生细菌群落主要由Allorhizobium–Neorhizobium–Pararhizobium–Rhizobium复合群和Pseudomonas属主导。值得注意的是，Pseudomonas的相对丰度在MRD点位显著高于SRD点位，并且与RD程度呈显著正相关，显示出其对环境胁迫的强耐受性。而具有固氮能力的Allorhizobium–Neorhizobium–Pararhizobium–Rhizobium复合群的相对丰度则与土壤pH呈显著负相关。

功能预测(PICRUSt2)显示，碳水化合物代谢是丰度最高的代谢类别之一，其在SRD样品中的预测基因丰度高于MRD样品，这与较高的SRD微生物多样性模式一致。

3.3. 土壤因子塑造RD生境中豌豆根瘤核心内生菌群

核心细菌属的数量从SRD点的47个下降到MRD点的35个，表明RD梯度产生了环境过滤效应。在所有四个样品中共有的29个核心属占总高质量序列的98.71%，表明存在一个跨RD梯度的稳定核心菌群，对共生功能至关重要。

冗余分析(RDA)表明，土壤pH和总磷(TP)是影响整个内生细菌群落沿RD梯度变化的关键环境因子，共同解释了79.51%的群落变异。这暗示土壤酸化和磷有效性是驱动豌豆根瘤内生菌群落演替的胁迫因子。

3.4. 土壤因子驱动豌豆根瘤固氮微生物的多样性格局与构建

通过对nifH基因的测序分析发现，携带nifH基因的内生菌丰富度和多样性在不同样品间存在差异。SRD点位S_JS3的Chao1指数和Shannon指数最高。物种丰富度指数与土壤pH呈显著正相关，Shannon指数与有效磷(AP)呈显著正相关。

在属水平上，Rhizobium是绝对优势的固氮菌，在四个点位的序列占比分别为83.8%、34.5%、13.6%和81.2%。主坐标分析(PCoA)显示，不同样品组的固氮内生菌群落结构存在显著差异。相关分析显示，Rhizobium的丰度与土壤pH呈负相关。Bradyrhizobium的丰度与TN、TP、AN呈负相关。

冗余分析(RDA)进一步揭示，土壤全氮(TN)、总磷(TP)和土壤含水量(SWC)是驱动固氮细菌群落沿RD梯度变化的关键环境因子，其中TN的影响最强。

4. 讨论

4.1. 石漠化程度本身并非土壤养分的唯一决定因素

研究发现，中度石漠化点位M_GJ7的土壤氮含量反而最高，这表明RD程度并非土壤性质的唯一决定因素。土壤养分水平是RD程度、海拔、地形、植被和人为活动等多因素综合作用的结果。M_GJ7点位海拔（2340.59米）显著高于其他点位，这种海拔差异可能通过影响小气候和植被进而影响养分循环。因此，在评估RD的生态效应时，应将海拔作为关键协变量纳入分析框架。

4.2. RD对植物内生菌多样性的影响

研究显示，RD强度与豌豆根瘤内生细菌多样性相关。轻度石漠化地区内生菌多样性更高，预测的碳水化合物代谢潜能也更强，这可能为耗能的固氮过程提供能量。随着RD加剧，多样性下降，这可能反映了环境过滤作用，筛选出耐胁迫的类群。

群落分析发现，Pseudomonas的丰度在MRD地区显著高于SRD地区，并与RD程度呈正相关。我们推测，随着石漠化加剧，相关的土壤胁迫（如pH和磷限制）可能抑制了核心共生固氮菌，而有利于像Pseudomonas这样的耐胁迫类群。这使得植物-微生物互作可能从以固氮为主导，转变为以耐胁迫为主导。这种功能群的演替在盐碱、重金属等其他胁迫生态系统中也有报道。

4.3. RD对植物内生固氮菌群结构的影响

在RD地区，生态系统退化与氮素流失显著。豆科植物与根瘤菌的共生固氮可以改善土壤结构和养分状况。本研究中，固氮菌的丰度主要受TP、TN和SWC影响，其中TN效应最强。Shannon指数与AP显著正相关，表明有效磷显著增强了固氮内生菌的物种多样性。Rhizobium是豌豆根瘤中的优势固氮内生菌。

土壤氮、磷有效性是决定固氮菌群落结构和功能的核心因子。固氮是一个消耗ATP的生物过程，其酶促反应和细胞结构的构建都需要磷的参与。因此，土壤磷和氮的有效性，而非RD严重程度本身，是驱动豌豆根瘤内生固氮菌群落构建的关键化学因子。此外，根瘤菌携带的phoD基因编码的碱性磷酸酶可矿化有机磷，在喀斯特钙质土壤中，释放的磷酸根与丰富的Ca²⁺离子反应沉淀，可能间接促进碳酸盐岩的溶解，为恢复这些脆弱生态系统提供了一个潜在的微生物机制。

4.4. 研究的局限性与普遍性

本研究存在一定局限性。受限于豌豆在RD区域的分布，采样点仅涵盖四个地点（SRD和MRD各两个），每个地点三个重复。采样点间存在海拔等环境条件的差异，这些因素与RD程度在空间上存在共线性。因此，本研究核心发现（如Allorhizobium–Neorhizobium–Pararhizobium–Rhizobium类群相对丰度下降和Pseudomonas富集）应谨慎解释为沿复杂环境梯度观察到的关联，而非RD过程本身的直接因果效应。尽管如此，这项探索性工作为未来的因果机制研究（如跨海拔梯度的分层嵌套采样、同质园实验、合成群落接种试验）提供了生态学假设和微生物资源靶点。

5. 结论

本研究表明，石漠化梯度通过改变土壤理化性质，与豌豆根瘤内生及固氮细菌群落的变化相关。Allorhizobium–Neorhizobium–Pararhizobium–Rhizobium和Pseudomonas占主导地位，其中Rhizobium是主要的推定固氮菌。土壤pH、含水量、氮和磷是与群落组成相关的关键环境因子。与中度石漠化点位相比，轻度石漠化点位的根瘤表现出更高的内生菌多样性和预测的碳水化合物代谢潜力。Pseudomonas和Rhizobium沿胁迫梯度的差异分布，表明微生物接种剂在改善豆科植物对退化喀斯特环境适应性方面具有潜力。然而，考虑到RD强度与海拔在空间上的共线性，观察到的模式代表的是沿复杂梯度的关联，而非因果效应。因此，这些类群的促生长功能及其在RD适应中的具体作用，有待通过靶向分离、受控接种和分层采样等未来研究进行验证。