### 中文标题 沿草地退化梯度的土壤微生物共现网络：非线性阈值、细菌-真菌差异响应及环境驱动因素

《Agriculture》：Soil Microbial Co-Occurrence Networks Along a Grassland Degradation Gradient: Nonlinear Thresholds, Divergent Bacterial–Fungal Responses, and Environmental Drivers

　　

草地退化是一场全球性的生态危机，深刻改变了地上植被和土壤性质，然而其对土壤微生物共现网络的影响仍知之甚少。研究人员利用Illumina MiSeq测序和共现网络分析，沿一条明确界定的五阶段退化梯度（从非退化羊草草地到极度退化的裸斑盐斑），调查了中国东北松嫩草甸草原的土壤细菌和真菌群落。研究结果表明，细菌α-多样性（α-diversity）呈现单峰（驼峰状）响应，在中度退化（MD）阶段达到峰值，而真菌多样性沿梯度单调下降，表明真菌对退化胁迫更敏感。细菌和真菌群落组成均随退化方向性变化，主要由土壤碱化（pH）和电导率（EC）驱动。两个界别的网络复杂性均遵循单峰模式，在MD阶段最大化，在重度退化（SD）阶段崩溃，表明中度与重度退化之间存在关键的生态阈值，而重度退化下正相关比例增加暗示了微生物为应对环境胁迫而增强的合作。结构方程模型进一步揭示了不同的调控路径：细菌网络受直接环境过滤和间接多样性介导效应的共同支配，而真菌网络对直接的pH/EC约束和组成变化响应更强。研究人员的研究结果表明，微生物网络对草地退化表现出非线性阈值响应，其中真菌比细菌更敏感，并强调了将网络层面属性整合到退化监测和恢复框架中的重要性。

### 论文解读文章

#### 研究背景与问题
草地退化是全球性生态危机，影响约40%的陆地表面，威胁养分循环、生物多样性和土壤稳定。尽管已有研究关注退化对微生物群落组成的影响，但微生物共现网络（co-occurrence networks）沿完整退化梯度的响应机制仍不清楚，尤其是细菌与真菌的差异响应、非线性阈值以及环境驱动路径（如直接土壤过滤 vs. 间接植被介导效应）尚未系统阐明。为此，研究人员在松嫩草甸草原（位于吉林松嫩草地生态系统国家观测研究站）沿五阶段梯度（非退化ND→轻度退化LD→中度退化MD→重度退化SD→极度退化ED）开展研究，旨在揭示微生物网络变化规律、关键驱动因素及潜在生物指示物。该论文发表在《Agriculture》。

#### 关键技术与方法
研究样本采集自松嫩草地生态系统国家观测研究站（44°45′ N, 123°45′ E），沿五阶段退化梯度设置6个重复样地（25 m × 25 m）。主要技术方法包括：
- **Illumina MiSeq测序**：对土壤细菌16S rRNA基因V3-V4区及真菌ITS1-ITS2区进行扩增子测序，生成操作分类单元（OTU）。
- **共现网络分析**：基于Spearman秩相关（|ρ| > 0.6, p < 0.01），在属水平构建网络，计算节点、链接、平均度、模块性等拓扑参数。
- **偏最小二乘结构方程模型（PLS-SEM）**：整合土壤pH、EC、植被生物量、微生物α-多样性（α-diversity）和β-多样性（β-diversity），评估对网络属性的直接与间接效应。

#### 研究结果
**3.1 植被与土壤性质沿退化梯度的变化**
通过野外调查与实验室分析，发现植物物种丰富度从ND的5.67降至ED的近零（F = 32.10, p < 0.001），总生物量在MD后急剧下降，凋落物生物量减少约84%（F = 105.87, p < 0.001）。土壤pH从ND的8.38升至MD-ED的9.70-10.01（F = 212.50, p < 0.001），Na⁺浓度增加235倍，土壤有机碳（SOC）下降37.7%，而NO₃^--N和有效磷（AP）在ED异常积累，表明养分循环解耦。

**3.2 土壤微生物多样性的变化**
细菌α-多样性（α-diversity）呈单峰模式，Shannon指数在MD最高后于SD骤降；真菌α-多样性则单调下降（F = 13.55, p < 0.001）。NMDS排序与PERMANOVA分析显示，退化阶段解释了细菌和真菌群落变异的61.2%和58.3%（p < 0.001），且真菌群落离散度显著增加（F = 22.839, p = 0.001），表明真菌组成更易受退化影响。

**3.3 退化诱导的微生物共现网络拓扑变化**
细菌和真菌网络的节点数和链接数均呈单峰模式：细菌网络从ND（1220节点，12860链接）增至MD（1311节点，16989链接），后在SD降至768节点、6993链接；平均度在MD最高（25.918），SD最低（18.211）。真菌网络链接从LD（4340）增至MD（6584），后降至SD（2800），平均路径长度在MD最小（3.620）、SD最大（5.764），表明网络断裂。关键阈值出现在MD至SD过渡期，链接数下降约59.5%（细菌）和57.5%（真菌）。

**3.4 微生物枢纽类群的更替**
通过节点度和中介中心性鉴定枢纽类群（hub taxa），发现细菌早期阶段以Hirschia、Taibaiella、Tumebacillus为主，MD阶段出现Actinobacteriota相关枢纽（如Phaselicystis、Glycomyces），SD阶段则由耐逆类群（如Acidobacteria相关类群、Desulfuribacillus）主导。真菌枢纽从Paraphaeosphaeria、Sarocladium更替为Plenodomus、Colletotrichum、Preussia等耐逆类群，表明网络简化伴随功能重组。

**3.5 微生物网络变化的驱动因素**
PLS-SEM模型解释了细菌和真菌网络属性（Net property，即网络拓扑主成分）变异的99.2%和90.6%。对于细菌，土壤pH和EC显著影响α-多样性（α-diversity）和β-多样性（β-diversity），β-多样性（β-diversity）对网络属性有强正效应（β = 0.97, p < 0.001），表明组成更替是主要驱动路径。对于真菌，pH和EC直接调控α-多样性（α-diversity）和β-多样性（β-diversity），β-多样性（β-diversity）同样与网络属性紧密关联（β = 0.73, p < 0.001），而植被生物量影响不显著，暗示环境过滤主导真菌网络重组。

#### 讨论与结论
讨论部分指出，土壤碱化和盐积累（pH、EC）是主要环境过滤器，驱动微生物群落组成与网络结构变化。细菌α-多样性（α-diversity）的单峰模式与中度干扰假说（IDH）一致，但需谨慎解释；真菌多样性单调下降反映其依赖于植物宿主和更弱的渗透调节能力，可作为更敏感的退化指示生物。网络复杂性在MD至SD阶段崩溃，但正相关比例和聚类系数增加，符合胁迫梯度假说（SGH），表明微生物在重度胁迫下增强合作。细菌网络受环境过滤和组成更替共同调控，真菌网络则更依赖直接pH/EC约束和组成变化。
**结论翻译**：总体而言，研究结果支持了三个假设。第一，微生物共现网络沿退化梯度呈非线性响应，网络复杂性在中度退化阶段达到峰值，并在向重度退化过渡时崩溃，证实退化诱导的网络简化在生态阈值后最为显著。第二，真菌群落对退化的敏感性高于细菌，表现为真菌多样性持续下降和更强的组成更替，支持真菌更易受植被损失和环境胁迫的假设。第三，土壤碱化和盐化被确定为调控微生物网络变化的主要驱动因素，而植被介导效应通过微生物多样性变化间接贡献，证实了环境过滤机制。这些结果强调，微生物网络重构和枢纽类群更替是与草地退化相关的敏感响应，经过更广泛时空尺度验证后，可为未来微生物指示物的开发提供候选。