**论文解读：土壤真菌α多样性与草地生产力在自然干旱梯度上的关联**

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

草地覆盖地球陆地表面约四分之一，支撑着碳固存、气候调节、生物多样性维持和畜牧业生产等关键生态系统功能。这类生态系统对气候变率高度敏感，尤其是对全球气候变化加剧下降水格局的变动。干旱频率和强度的增加已引发草地稳定性和生产力的广泛退化，对陆地碳平衡和区域粮食安全造成严重后果。净初级生产力（NPP）作为生态系统功能的核心指标，对水分可用性响应强烈，但其变化幅度和方向因草地类型和环境背景而异。因此，理解干旱化背景下NPP的生物和环境决定因素，对于预测未来气候情景下草地的韧性至关重要。

土壤微生物群落是养分循环、植物生长和生态系统稳定的重要组成部分，构成了草地响应干旱的关键媒介。微生物多样性和群落组成影响植物养分获取、竞争动态、分解过程和土壤碳周转。然而，细菌和真菌在耐旱策略、生态功能和对环境过滤的响应上存在显著差异：细菌通常表现出快速的更替和生理可塑性，而真菌凭借广泛的菌丝网络和孢子传播机制，可能对水分胁迫具有更强的抵抗能力。这些功能差异暗示细菌和真菌群落可能在干旱化过程中对生态系统功能做出不同的贡献。尽管已有研究记录了干旱引起的微生物群落组成和组装机制变化，但这些变化如何级联影响自然干旱梯度上的草地生产力仍不清楚。

内蒙古草地样带横跨自然的干旱梯度（多年平均降水量175–362 mm；多年平均温度?2.4至5.6°C），草地类型从草甸草原过渡到荒漠草原，为研究土壤微生物对干旱化的响应及其对生态系统功能的影响提供了理想的天然实验室。基于这一空间梯度，本研究整合了高通量测序、共现网络分析、零模型推断和结构方程模型，旨在评估干旱条件下微生物组的重组及其对草地生产力的影响。具体研究目标包括：（1）细菌和真菌多样性、组成及生物标志物类群的变迁；（2）微生物互作网络、关键类群和组装过程的变化；（3）干旱、土壤性质和微生物组调控地上和地下NPP的直接与间接路径。

**二、研究设计与主要结论**

研究在内蒙古四个横跨东西干旱梯度的样地开展，包括草甸草原（ET：119°58′E, 49°19′N；HS：119°30′E, 50°12′N）、典型草原（XY：116°33′E, 43°32′N）和荒漠草原（WL：106°58′E, 41°25′N）。四个样地的干旱指数（AI）分别为1.50（ET，半湿润）、1.41（HS，半湿润）、2.21（XY，半干旱）和6.53（WL，干旱），涵盖了从半湿润到干旱的条件。每个样地于2020年8月建立6个6 m×6 m样方，2021年8月采集植物和土壤样品。研究测定土壤理化性质，通过收获法估算地上净初级生产力（ANPP），采用根内生长环法测定地下净初级生产力（BNPP）。土壤DNA提取后，针对细菌16S rRNA基因（338F/806R）和真菌ITS区（ITS1F/ITS2R）进行Illumina MiSeq高通量测序。

研究得出以下主要结论：真菌系统发育多样性沿干旱梯度显著下降（从200.09降至86.03，降幅57.0%），而细菌多样性相对稳定；真菌网络复杂性和稳定性随干旱加剧而降低，细菌网络虽节点和边数增加但网络稳定性下降；细菌群落组装从确定性过程转向随机过程，而真菌群落的确定性过滤随干旱加剧而增强；干旱条件下耐旱类群（如壶菌门Chytridiomycota和球囊菌门Glomeromycota）富集；真菌多样性而非细菌多样性与ANPP和BNPP均表现出最强的正相关关系。

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

研究采用以下核心分析方法：（1）基于Illumina MiSeq平台的高通量测序技术，对细菌16S rRNA基因和真菌ITS区域进行扩增测序；（2）共现网络分析（Spearman相关性，|r|>0.6, p<0.05），结合igraph包和Gephi软件进行网络拓扑结构分析和可视化；（3）零模型推断，通过β-NTI和RC_bray量化确定性与随机性组装过程的相对贡献；（4）偏最小二乘路径模型（PLS-PM）评估干旱、土壤性质和微生物因素对草地生产力组分的直接与间接效应；（5）随机森林模型（randomForest包）、LEfSe分析和Wilcoxon秩和检验鉴定差异丰度类群和生物标志物。所有测序数据已提交至NCBI数据库（登录号PRJNA1226389）。

**四、研究结果**

**4.1 微生物群落组成与多样性沿干旱梯度的变化**

在细菌群落中，放线菌门（Actinobacteriota）、变形菌门（Proteobacteria）、酸杆菌门（Acidobacteriota）和绿弯菌门（Chloroflexi）占总序列的81.15–85.04%。随干旱加剧，绿弯菌门、芽单胞菌门（Gemmatimonadota）、粘球菌门（Myxococcota）、厚壁菌门（Firmicutes）和蓝藻门（Cyanobacteria）的相对丰度增加，而变形菌门、疣微菌门（Verrucomicrobiota）和甲基暖菌门（Methylomirabilota）下降。真菌群落中，子囊菌门（Ascomycota）、担子菌门（Basidiomycota）和被孢霉门（Mortierellomycota）占85.29–96.59%，壶菌门（Chytridiomycota）和球囊菌门（Glomeromycota）在高干旱条件下上升而担子菌门下降。细菌Shannon多样性在WL和XY样地高于ET和HS样地，而真菌系统发育多样性沿干旱梯度显著下降。主坐标分析（PCoA）显示四个区域的细菌和真菌群落均呈现明显分离。

**4.2 微生物组中显著富集的生物标志物**

LEfSe分析鉴定出四个干旱区域从门到属水平显著富集的细菌和真菌类群。细菌在ET、HS、XY和WL中显著富集的类群数分别为18、8、11和21，真菌对应的数量为12、13、10和26。Wilcoxon秩和检验显示细菌属Rubrobacter在四个区域间存在显著差异；真菌属Mortierella在ET和WL之间、Knufia在HS和WL之间存在显著差异。随机森林分析确定了预测区域干旱度的关键生物标志物：细菌的前三位为g_unclassified_o_Burkholderiales、g_unclassified_o_Ardenticatenales和g_unclassified_c_AD3；真菌的最重要生物标志物为Pleotrichocladium、Ilyonectria和Geoglossum。共鉴定出72个细菌OTU和21个真菌OTU作为关键判别类群。五个细菌类群（g_unclassified_f_67-14、g_unclassified_f_Vicinamibacteraceae、g_unclassified_o_Rokubacteriales、g_unclassified_f_JG30-KF-CM45和g_unclassified_o_0319-7L14）和三个真菌类群（Clavaria、Microdochium和Monosporascus）被确定为跨区域共享类群。

**4.3 微生物互作网络的复杂性与稳定性**

共现网络分析显示细菌和真菌群落对干旱的响应存在显著差异。细菌网络的节点数随干旱加剧而增加（从ET的197增至WL的222），边数从ET的2,144增至HS的2,789；平均度在HS最高（13.74）而在ET最低（10.88）。真菌网络方面，ET具有最高的节点和边数（173和1,096），XY最低（143和786）；平均度沿干旱梯度从ET的6.34降至WL的5.43。这表明真菌网络对干旱比细菌网络更敏感，干旱显著降低了真菌网络复杂性。网络枢纽节点鉴定结果表明，ET、HS、XY和WL中分别有14、27、13和38个枢纽节点。细菌群落稳定性（以平均变异度AVD衡量）随干旱加剧显著下降，而真菌群落稳定性在四个区域间无显著变化。

**4.4 关键类群的分布与功能特征**

在门水平上，细菌关键类群中绿弯菌门（Chloroflexi）和芽单胞菌门（Gemmatimonadota）的相对丰度随干旱加剧而增加，而疣微菌门（Verrucomicrobiota）和甲基暖菌门（Methylomirabilota）逐渐减少。蓝藻门（Cyanobacteria）仅在最高干旱区域（WL）被鉴定为关键门，丰度为1.03%；甲基暖菌门主要存在于较低干旱区域（ET为3.38%，HS为1.84%）。真菌方面，壶菌门（Chytridiomycota，6.52%）和球囊菌门（Glomeromycota，5.26%）仅在WL被鉴定为关键门，且真菌关键类群随干旱加剧而变得更加多样化。

真菌关键类群的功能预测表明，其前三位功能为参与三磷酸腺苷酶、L-阿拉伯糖异构酶和葡聚糖1,4-α-葡萄糖苷酶的合成。高干旱区域的真菌关键类群功能丰度显著高于低干旱区域。细菌关键类群的主要预测功能为化能异养、好氧化能异养和硝酸盐还原，而光合营养功能（包括光异养、光合营养和光自养）的丰度随干旱加剧而增加。

**4.5 细菌和真菌群落的组装过程与生物地理模式**

β-NTI与RC_bray的整合分析揭示了沿干旱梯度细菌和真菌群落组装机制的差异。细菌β-NTI从ET的3.52和HS的3.30逐渐降至XY的2.82和WL的0.96，而RC_bray在所有样地均为负值（?0.82至?0.52），表明细菌群落从草甸草原的确定性过程（异质选择）转变为荒漠草原的随机过程（同质扩散），XY代表中间过渡状态。真菌β-NTI在所有样地均接近零（?0.11至0.24），而RC_bray从ET的0.87和HS的0.85降至XY的0.77和WL的0.06，表明随机过程在整个干旱梯度上持续主导真菌群落组装，但随着干旱加剧，确定性过程的影响逐渐增强。距离衰减分析显示，随着地理距离增加，细菌的|ΔAVD|和β-NTI显著增加而|Δ系统发育多样性|显著降低，真菌的|ΔAVD|显著下降而|Δ系统发育多样性|和β-NTI均明显增加。

**4.6 净初级生产力的生物与非生物驱动因素**

偏最小二乘路径模型（PLS-PM）结果显示，干旱对ANPP和BNPP同时存在直接和间接影响。干旱与土壤物理变异性及细菌群落变异呈正相关，与土壤养分可用性呈负相关。土壤养分与细菌和真菌多样性正相关，与真菌关键类群丰度负相关。真菌多样性是ANPP和BNPP最强的正相关因素，而真菌关键类群与BNPP存在特异性关联。干旱通过正相关于真菌关键类群丰度对BNPP产生间接正效应，通过负相关于真菌多样性对两个生产力组分产生间接负效应。模型将干旱确定为ANPP和BNPP的最强负相关因素，将真菌多样性确定为最强正相关因素。

三项补充分析验证了这些发现：偏相关分析在控制土壤理化性质（pH、容重、总磷、速效磷、总氮、土壤有机碳、NO₃^?-N和NH₄⁺-N）后，真菌Shannon多样性与BNPP仍显著相关（偏r=?0.5299, p=0.0348），但ANPP相关性不显著；逐一剔除样地（leave-one-site-out）敏感性分析表明"真菌多样性→生产力"的核心路径在大多数子样本中方向一致，但剔除最湿润（ET）或最干旱（WL）样地时方向反转，表明该关联具有环境依赖性；随机森林分析确定了pH和关键养分（碳、氮、磷）为微生物多样性的主要预测因子。

**五、讨论与结论**

讨论部分指出，干旱引起的养分限制显著改变了土壤微生物群落的结构和组成。真菌系统发育多样性的急剧下降可能归因于干旱条件下土壤有机碳可用性的减少，这对真菌的影响强于细菌。细菌和真菌对干旱的差异响应源于其固有的微生物策略差异：细菌依赖快速的生理可塑性和抗逆性状（如内生孢子形成），而真菌受益于菌丝觅食网络。

关于网络特征，细菌网络在干旱条件下表现出节点和边数增加但网络稳定性下降，真菌网络则表现出节点数和平均度降低，表明真菌网络对干旱胁迫更为敏感。关键类群的分布变化反映了不同类群的生理适应策略：绿弯菌门凭借生理韧性和代谢灵活性在干旱中增加，壶菌门仅在最高干旱位点出现，体现其耐旱性、低养分需求和低生物压力下的竞争优势。

关于组装过程，细菌群落从确定性组装转向随机组装可能与极端pH和土壤性质变化有关，而真菌群落中确定性过滤的增强表明气候改变的资源可用性驱动了更强的环境筛选。细菌作为r-策略者能快速响应环境变化，而真菌作为K-策略者生长较慢但资源利用效率更高，可能强化了相对中性的过程。

关于生产力关联，研究发现真菌多样性而非细菌多样性与草地生产力表现出最强的正相关，这与真菌通过提高资源利用效率和植物互作调节生物多样性-生产力关系的角色一致。PLS-PM与偏相关分析结果方向的差异可能反映了分析框架的不同：PLS-PM捕捉变量交互网络中的总效应，而偏相关分离了控制协变量后的直接残差关联。负的直接关联可能提示资源竞争，而正的总效应可能由真菌关键类群介导的间接路径所主导。leave-one-site-out分析揭示这种关系具有环境依赖性，且由于缺乏植物群落数据，观察到的正相关关系可能部分反映了植物群落沿梯度的更替。

研究局限性包括单一年份单一季节的采样设计、仅分析浅层土壤、以及缺乏植物群落数据。未来的研究需纳入多年、多季节、深剖面采样和全面的植被调查。

**研究结论：** 本研究表明，日益加剧的干旱与土壤微生物群落的根本性转变及其与草地生产力的关系密切相关。真菌群落，尤其是其多样性和关键类群，对干旱的敏感性远高于细菌群落，并与ANPP和BNPP均表现出强烈的正相关。相比之下，细菌多样性保持相对稳定，尽管其网络稳定性和组装过程发生了变化。结构方程分析表明，真菌多样性的丧失与干旱引起的生态系统生产力限制密切相关。总体而言，这项工作突出了土壤真菌群落与旱地生态系统功能和韧性的紧密联系，为理解草地应对持续气候变化的响应以及考虑基于微生物组的减缓策略提供了基础。