青海高原蚕豆根际细菌群落多样性、稳定性及生态位宽度沿海拔梯度的分布格局

《Microorganisms》：Altitudinal Patterns of Diversity, Stability, and Niche Breadth in Faba Bean Rhizosphere Bacterial Communities on the Qinghai Plateau

　　

本研究通过对青海高原蚕豆种植区开展大范围调查，表征了蚕豆（Vicia faba L.）根际细菌群落的区域变异特征，并检验其与土壤及空间因子的关联。研究在青海省海东市、西宁市及海南藏族自治州15个蚕豆田间位点采集根际土壤样本，海拔跨度1842.9至3206.3 m，每个田块设9个嵌套子样本，以田块为独立重复单元。测定土壤理化性质，并对细菌16S rRNA基因进行高通量测序。位点水平上，三组海拔区间的细菌α多样性（alpha-diversity）指数呈数值波动但无显著差异；基于Bray–Curtis相异度的位点水平PERMANOVA显示群落组成在组间差异显著（pseudo-F2,12 = 1.789，R2 = 0.230，p = 0.031），但位点间重叠仍较大，保守位点水平敏感性分析不显著（p = 0.060）。地理距离与细菌群落相异度显著相关，而控制地理距离后海拔差无独立可检关联，表明空间结构对群落分化贡献显著。优势菌门组成跨位点广泛保守，属水平相对丰度呈描述性空间变异。基于RAVD（Reciprocal of Average Variation Degree，平均变异度倒数）的组分一致性与群落水平生态位宽度（niche breadth）组间均无显著差异。所测土壤性质中仅pH差异显著，多数养分变量无显著差异。总体表明，青海高原蚕豆根际细菌群落呈空间结构化区域分化，观测格局可能反映地理结构、海拔、背景土壤条件、农业管理、植物物候及其他未测环境因素的联合作用，本次调查中无法独立解耦。

研究背景方面，青藏高原作为生态脆弱区对全球气候变化高度敏感，高海拔农业系统面临低温、强辐射与短生长季等不同于低海拔的约束。蚕豆是青海重要粮饲豆科作物，依托与根瘤菌共生固氮，其根际微生物参与有机质分解、养分转化与植物适应过程。海拔通过气候、土壤与植物因素联合调控微生物多样性与组成，但已有报道在山地生态系统中结论不一；同时农业生态系统叠加作物身份与管理扰动，使豆科作物根际细菌沿海拔的多样性、稳定性及生态位宽度研究尤为不足。鉴于此，研究人员在青海高原蚕豆主产区布设15个田块位点，系统刻画细菌群落沿海拔（低海拔1800–2300 m、中海拔2300–2800 m、高海拔2800–3300 m）的分化，并发表于《Microorganisms》。

关键技术方法上，研究人员于2024年7月25日在青海海东、西宁、海南州15个蚕豆田块各设9个1 m×1 m样方，采集0–15 cm土柱测理化，并以根表0–2 mm附着土为根际样本；以田块为独立重复（每海拔组5位点，共135子样本）。土壤有机碳（SOC）、pH、全氮（TN）、碱解氮（AN）、全磷（TP）、速效磷（AP）、全钾（TK）、速效钾（AK）按常规法测定。细菌群落通过FastDNA? SPIN Kit提取DNA，扩增16S rRNA V4区（引物515F/806R），Illumina MiSeq PE300测序，UPARSE聚OTU（97%相似度），SILVA 138注释。统计上采用嵌套设计，位点水平聚合后做ANOVA、Tukey HSD、NMDS、PERMANOVA、Mantel及偏Mantel；稳定性用RAVD，生态位宽度用Levins指数（B）及群落加权Bcom。

研究结果部分，3.1节土壤理化变异显示SOC与TN在中高海拔显著升高，AN在中海拔达峰，TP与AK随海拔下降，pH组间差异显著（F_2,12 = 7.781，p = 0.0068），高海拔pH显著高于中海拔。3.2节α多样性Richness、Shannon、Pielou均匀度数值上随海拔降低（LF最高），但位点水平ANOVA不显著；NMDS与PERMANOVA表明群落组成组间显著分化（pseudo-F_2,12 = 1.789，R² = 0.230，p = 0.031），地理距离偏Mantel检验显示海拔无独立效应。3.3节组成分析中，菌门级Pseudomonadota、Actinomycetota、Acidobacteriota、Bacteroidota跨组保守，Pseudomonadota占比随海拔微增；属级Allorhizobium–Neorhizobium–Pararhizobium–Rhizobium组在高海拔均值高但因SAF、SBF两位点集中（占88.37%）且OTU21444（Rhizobium phaseoli）占90.19%，属描述性空间异质。3.4节Mantel显示AN与TP在三组均与群落相异度关联，LF另关TP，MF广关联SOC/AN/AP/TK/TP/TN，HF关联SOC/AN/AK/AP/TP/TN/pH，但属探索性模式。3.5节RAVD在LF、MF、HF均值为2.249、2.547、2.466，ANOVA不显著（F_2,12 = 1.851，p = 0.199）。3.6节Bcom均值61.77、62.54、54.27，高海拔数值最低但ANOVA不显著（F_2,12 = 1.754，p = 0.215）。

讨论部分总结指出，海拔并非孤立驱动因子，而是与温度、湿度、pH、养分矿化及植物输入协变；α多样性与稳定性、生态位宽度未呈显著海拔递变，提示高海拔强过滤假设在本农业系统中未获支持。群落组成分化虽显著但R值低、重叠大，反映宿主选择、农业管理与局部空间异质性的共享作用。根瘤菌高丰度聚集于个别高海拔位点，不能推广为全高海拔富集。 bulk soil（本体土）理化与根际菌群存在空间错配，关系应作位点水平关联解读而非根际直接因果。

结论部分翻译：本研究揭示了青海高原15个1842.9–3206.3 m蚕豆田块根际细菌群落的明确区域与空间分化。α多样性组间数值波动但位点水平不显著；群落组成具可检区域分化但位点重叠大；地理距离与相异度显著相关，控制后海拔无独立关联，表明格局受强空间结构化。优势菌门广泛保守，属级呈描述性空间变异；RAVD一致性及群落生态位宽度组间不显著。总体应将变异解释为跨大海拔幅度的空间结构化区域分化，而非海拔独立效应。地理结构、海拔、背景土壤、农业管理、植物物候及其他未测因素协变且难解耦。未来需空间均衡采样、季节重复、根际专一化学、气候观测与功能解析以厘清贡献。