来源：Native Pseudomonas inoculants alleviate microbial resource limitations and enhance plant biomass in degraded alpine meadow soils（Journal of Environmental Management）| 编译：DNA Lab Space | 原文许可：elsevier-subscription
导读
青藏高原高寒草甸退化破坏碳氮磷循环，造成微生物资源限制并制约植物生长。本研究通过三年田间试验，评估由四种本地假单胞菌组成的复合接种剂在退化梯度上的修复效果。结果表明，接种剂持续缓解微生物磷限制、增强碳氮获取酶活性、提升微生物生物量与植物地上地下生物量，且调控路径随退化程度而异——中度退化土壤以缓解碳限制为主，重度退化土壤则以微生物量积累为前提。
研究背景
青藏高原拥有全球面积最大的高寒草地生态系统，覆盖约250万平方公里，在生物多样性保育、碳固持和亚洲主要河流水系调节方面提供关键生态系统服务。然而，过度放牧、气候变暖以及高原鼠兔等小型哺乳动物扰动，已引发大面积高寒草甸退化。退化表现为植被盖度下降、土壤有机碳与氮流失、养分循环紊乱。据报道，青藏高原超过60%的高寒草甸已发生不同程度退化，严重退化区原始植被盖度损失可达80%。这一退化削弱了生态系统碳汇潜力、水文调节功能与牧区支撑能力，对区域生态安全与可持续发展构成重大挑战。全球其他寒冷草地地区，包括安第斯普纳、喜马拉雅等地，也面临类似退化问题。本地微生物接种剂因环境风险低，被视为退化土壤修复的潜在替代方案，但其在退化梯度上的作用机制尚不明确。
研究方法
2.1 研究地点
研究在青藏高原东部甘肃省玛曲县大水种畜场进行（34°0′33″N，102°16′6″E），海拔3420 m。该区域属高原大陆性气候，年均气温1.2°C（范围−28°C至23.6°C），年均降水量594 mm（主要集中在6月至9月），年均蒸发量1355.7 mm（Cao et al., 2012）。优势植被包括嵩草属、羊茅属、早熟禾属、披碱草属、风毛菊属和橐吾属物种，以及其他草本如高山金莲花、钝裂银莲花、细叶葱和花苜蓿。畜牧业经济活动叠加高原鼠兔扰动，导致景观尺度上出现不同程度退化（Niu et al., 2009）。基于植被盖度、物种组成、地上生物量和裸地面积，本研究选取三个退化水平：未退化（ND）、中度退化（MD）和重度退化（SD）（Li et al., 2021, 2024a）。其他样地特征见表S1。根据FAO世界土壤资源参考基础，研究样地土壤分类为亚高山草甸土（相当于暗色土），发育于第四纪冲积和湖积沉积物。0–10 cm土层的详细剖面描述和基本性质见表S2。这些退化土壤代表生物资源应用的目标生境，其独特的理化性质指导了针对不同退化程度的接种策略。
【数据】地点：甘肃省玛曲县大水种畜场；坐标：34°0′33″N，102°16′6″E；海拔：3420 m；年均气温：1.2°C；温度范围：−28°C至23.6°C；年均降水量：594 mm；年均蒸发量：1355.7 mm；退化水平：ND、MD、SD；土壤类型：亚高山草甸土（暗色土）；采样深度：0–10 cm
2.2 微生物接种剂生产、田间试验设计与采样
#### 2.2.1 微生物菌株分离与表征
微生物接种剂（MI）由四株细菌组成，采用标准稀释涂布法从青藏高原高寒草甸优势植物（嵩草属、垂穗披碱草）根际分离，使用选择性培养基（无氮苹果酸盐培养基用于 putative 固氮菌；Pikovskaya琼脂用于解磷菌）。分离物经反复划线纯化，并参照Li et al.（2020）表征植物促生特性。四株选定菌株经16S rRNA基因测序鉴定为 Pseudomonas marginalis PGPR1、Pseudomonas silesiensis PGPR2、Pseudomonas antarctica PGPR3 和 Pseudomonas putida PGPR4（图S1）。其功能特征包括：（1）固氮（乙炔还原活性：PGPR1和PGPR4为229.47–545.23 C2H4 nmol·mg−1 h−1）；（2）无机磷酸盐溶解（131.26–438.35 μg mL−1）；（3）有机磷矿化（13.36–35.46 μg mL−1）；（4）吲哚-3-乙酸（IAA）产生（PGPR1和PGPR4为11.72–41.04 μg mL−1）。详细菌株特征见表S3。
【数据】菌株：Pseudomonas marginalis PGPR1、Pseudomonas silesiensis PGPR2、Pseudomonas antarctica PGPR3、Pseudomonas putida PGPR4；固氮活性：229.47–545.23 C2H4 nmol·mg−1 h−1；无机磷溶解：131.26–438.35 μg mL−1；有机磷矿化：13.36–35.46 μg mL−1；IAA产量：11.72–41.04 μg mL−1；分离培养基：无氮苹果酸盐培养基、Pikovskaya琼脂
#### 2.2.2 接种剂生产与配方
接种剂生产时，各菌株分别在含100 mL LB肉汤（pH 7.0）的250 mL三角瓶中、28°C、180 rpm振荡培养48 h，至对数晚期（OD600 ≈ 1.8，对应约1.0 × 10⁹ CFU mL−1）。扩大培养时，将5 mL各种子培养液转移至含200 mL改良LB培养基的500 mL挡板瓶中，在相同条件下培养。细胞通过8000×g、4°C离心10 min收获，用无菌磷酸盐缓冲液（PBS，pH 7.2）洗涤三次去除残留培养基，重悬于无菌PBS中至终浓度约1.0 × 10⁹ CFU mL−1（LB琼脂平板计数确定）。四株菌按等体积比（1:1:1:1）混合配制成复合液体接种剂。田间施用前立即通过平板计数验证接种剂活力；4°C储存7天内活力保持>95%。每年施用前新鲜制备接种剂以确保质量一致。
【数据】培养基：LB肉汤（pH 7.0）、改良LB培养基；培养温度：28°C；转速：180 rpm；时间：48 h；OD600：约1.8；浓度：约1.0 × 10⁹ CFU mL−1；离心条件：8000×g、4°C、10 min；洗涤液：无菌PBS（pH 7.2）；混合比例：1:1:1:1（体积比）；储存条件：4°C；活力：>95%（7天内）
#### 2.2.3 田间试验设计
田间试验于2021年5月中旬植被返青期启动，采用随机完全区组设计。该时间基于初步监测确定的微生物活性最适土壤温度（10 cm深度≥5°C）和水分条件（SWC约25–30%）选择。设置五个处理：（1）ND，不施MI；（2）MD，不施MI；（3）MD + MI（MD施用75 kg ha−1 MI）；（4）SD，不施MI；（5）SD + MI（SD施用75 kg ha−1 MI）。施用量（75 kg ha−1）基于前期温室试验计算（Li et al., 2020），对应约1.5 × 10¹² CFU ha−1。每个处理设三个重复小区（5 m × 5 m），间隔1 m缓冲带以减少交叉污染。MI于2021、2022和2023年每年5月中旬使用背负式喷雾器以土壤浇灌方式施用，随后轻灌（5 mm）促进渗透。研究期间试验区域围栏禁牧以防止外部扰动。
【数据】试验启动时间：2021年5月中旬；设计：随机完全区组；处理：ND、MD、MD + MI、SD、SD + MI；施用量：75 kg ha−1（约1.5 × 10¹² CFU ha−1）；小区面积：5 m × 5 m；缓冲带：1 m；重复：3；施用时间：2021、2022、2023年5月中旬；灌溉量：5 mm；土壤温度：≥5°C（10 cm）；SWC：约25–30%
#### 2.2.4 样品采集
植物和土壤采样遵循多年计划，旨在捕捉短期响应和年际动态：（1）2021年，在MI施用后30、60和90天采样，分别标记为I-30、I-60和I-90，其中I-90恰逢植物生长高峰（8月中旬）；（2）2022和2023年，仅在MI施用后90天（生长高峰期）采样，标记为II-90和III-90。该设计允许评估初始定殖动态（2021年）和重复施用下的长期处理效应（2022–2023年）。每个小区内设三个随机50 cm × 50 cm样方。因此，n = 3小区 × 每小区1个混合样 = 每处理每采样时间3个独立重复。地上生物量（AGB）在每个样方内齐地面剪取。根系小心挖至30 cm深度，用水冲洗去除土壤，105°C干燥30 min灭活酶，然后65°C烘干至恒重。植被收获后，使用10 cm直径土钻采集0–10 cm土层土壤样品（每小区五钻，混合）。每个混合样分为三个子样：一份风干、研磨、过筛（2 mm，然后0.25 mm）用于理化分析；第二份4°C储存≤2周用于生物学分析（酶活性和微生物生物量）；第三份−80°C储存用于后续DNA提取和16S rRNA基因扩增子测序。所有分析在储存期内完成以尽量减少样品降解。
【数据】采样时间：2021年I-30、I-60、I-90；2022年II-90；2023年III-90；样方：50 cm × 50 cm；重复：n = 3；根系深度：30 cm；干燥条件：105°C 30 min，然后65°C至恒重；土钻直径：10 cm；采样深度：0–10 cm；每小区土钻数：5；筛孔：2 mm、0.25 mm；储存条件：4°C（≤2周）、−80°C
2.3 土壤理化性质测定
土壤pH和电导率（EC）分别在1:2.5和1:5土水比悬浮液中测定，使用多参数pH/电导率仪（PCTestr™ 35，日本）。土壤含水量（SWC）通过105°C烘箱干燥新鲜样品至恒重的重量法测定。SOC采用K2Cr2O7氧化还原法分析（Nelson and Sommers, 1996）。全氮（TN）采用凯氏消解法定量（Keeny and Nelsson, 1982）。全磷（TP）和速效磷（AP）采用钼锑抗比色法测定（Olsen and Sommers, 1982）。溶解性有机碳（DOC）用0.5 M K2SO4提取，使用TOC分析仪（TOC-L CPN，Shimadzu，日本）按Brookes et al.（1985）分析。速效氮（AN）采用连续碱解扩散法测定。铵态氮（NH4+-N）和硝态氮（NO3−-N）用2 M KCl提取（1:5土液比），使用连续流动自动分析仪（Seal Analytical，AA3，德国）定量。所有分析设三次重复并含质量控制标准。
【数据】pH测定：1:2.5土水比；EC测定：1:5土水比；SWC：105°C烘箱干燥；SOC：K2Cr2O7氧化还原法；TN：凯氏消解法；TP/AP：钼锑抗比色法；DOC：0.5 M K2SO4提取，TOC分析仪；AN：连续碱解扩散法；NH4+-N/NO3−-N：2 M KCl提取（1:5），连续流动自动分析仪；重复：3
2.4 土壤微生物生物量和胞外酶活性测定
土壤微生物生物量碳（MBC）、氮（MBN）和磷（MBP）采用氯仿熏蒸提取法测定（Diaz-Raviña et al., 1992; Vance et al.）。
【数据】方法：氯仿熏蒸提取法；指标：MBC、MBN、MBP；原文未详述具体用量与培养条件
研究结果
3.1 植物生物量对假单胞菌接种剂的响应及伴随的土壤性质变化
连续三年每年施用MI，导致中度（MD）和重度（SD）退化高寒草甸地上（AGB）和地下（BGB）生物量发生显著变化（图1）。
两因素线性混合效应模型（小区身份为随机截距）揭示，处理、采样时间及其交互作用对AGB和BGB均产生极显著影响（均p < 0.001；表S4）。极高的偏η²值（>0.96）表明这些因素共同解释了植物生物量总变异的96%以上。MI在所有采样时间持续增加AGB和BGB。相对于MD和SD，MD + MI使AGB增加3.7–11.4%，BGB增加3.7–9.4%；SD + MI表现出更大的时间变异，AGB从−8.3%（I-30时不显著，p > 0.05）至20.1%，BGB从−2.8%至37.5%。MI逐步缩小了与未退化（ND）参照小区的生产力差距：MD + MI相对ND的AGB亏缺从45.6%降至39.0%，SD + MI的BGB亏缺在三年接种周期后从62.9%缩小至55.1%。这些植物响应伴随土壤理化性质的显著变化（表S5–S10）。相对于MD和SD，MI施用显著降低土壤pH和EC，同时增加SWC、总养分库（SOC、TN、TP）和所有有效养分组分（DOC、AN、NH4+-N、NO3−-N、AP）（p < 0.001）。土壤肥力与植物生产力的同步改善凸显了MI在修复退化高寒草甸生态系统中的综合作用。
【数据】处理效应：p < 0.001；偏η²：>0.96；MD + MI AGB增幅：3.7–11.4%；MD + MI BGB增幅：3.7–9.4%；SD + MI AGB变化：−8.3%至20.1%；SD + MI BGB变化：−2.8%至37.5%；MD + MI AGB亏缺：45.6%降至39.0%；SD + MI BGB亏缺：62.9%降至55.1%；土壤指标变化：p < 0.001
3.2 MI增强土壤微生物生物量并调控胞外酶活性和微生物碳利用效率
连续三年施用MI显著提升两种退化草甸土壤的微生物生物量碳（MBC）、氮（MBN）和磷（MBP）（表S11–S12）。相对于MD和SD，MD + MI使MBC增加9.6–32.4%，MBN增加4.6–36.1%，MBP增加7.8–22.0%；SD + MI范围较窄，分别为3.0–16.3%、1.5–36.1%和3.1–11.7%，反映严重初始碳耗竭制约了微生物扩增。所有微生物生物量库在晚季采样（I-90、II-90、III-90）达到峰值，与植物物候一致。胞外酶活性呈现处理特异性模式（图2；表S13–S14）。
MI在生长高峰期显著增强碳获取酶（BG、BX）和氮获取酶（NAG、LAP），同时抑制磷获取酶碱性磷酸酶（ALP）活性。通过Sinsabaugh et al.（2016）化学计量模型估算的微生物碳利用效率（CUE）在MI施用后下降：MD土壤从0.31降至0.22（p < 0.001），SD土壤从0.28降至0.19（p < 0.001；图S2；表S14）。
【数据】MD + MI MBC增幅：9.6–32.4%；MBN增幅：4.6–36.1%；MBP增幅：7.8–22.0%；SD + MI MBC增幅：3.0–16.3%；MBN增幅：1.5–36.1%；MBP增幅：3.1–11.7%；CUE：MD从0.31降至0.22（p < 0.001），SD从0.28降至0.19（p < 0.001）；酶活性：BG、BX、NAG、LAP增强，ALP抑制
3.3 MI通过化学计量调整缓解磷限制但加剧碳限制
矢量长度（碳限制代理）和矢量角度（氮 vs. 磷限制代理）受处理、采样时间及其交互作用显著影响（图3；表S14）。
在所有采样时间，ND土壤始终表现出最大矢量长度和最小矢量角度，表明相对于退化土壤，其固有碳限制更强、磷限制更弱。MI施用诱导两个指标的方向性变化：与MD和SD相比，MD + MI和SD + MI显著增加矢量长度同时降低矢量角度（图3a和b）。在III-90生长高峰期，MD + MI矢量长度增加19.5%（从1.18至1.41），矢量角度降低13.3%（从62.4°至54.1°）；SD + MI变化更大，矢量长度增加26.2%（1.07至1.35），矢量角度降低19.6%（65.7°至52.8°）。关键的是，所有处理和采样时间的矢量角度均保持在45°以上，证实即使经过三年MI施用，磷限制仍是这些高寒草甸土壤的背景约束（图4a）。
矢量长度与矢量角度之间观察到强负相关（R² = 0.840，p < 0.001；图4b），表明磷限制的缓解始终伴随碳需求的加剧。矢量长度和矢量角度分别与微生物CUE呈显著负相关和正相关（R² = 0.502，p < 0.001和R² = 0.538，p < 0.001；图4c和d）。此外，矢量长度与AGB（R² = 0.270，p < 0.001）和BGB（R² = 0.666，p < 0.001）呈显著正相关，而矢量角度与AGB（R² = 0.406，p < 0.001）和BGB（R² = 0.681，p < 0.001）呈显著负相关（图S3）。总体而言，这些关系与磷限制抑制有利于植物生产力的观点一致。
【数据】III-90 MD + MI矢量长度：1.18→1.41（+19.5%）；矢量角度：62.4°→54.1°（−13.3%）；III-90 SD + MI矢量长度：1.07→1.35（+26.2%）；矢量角度：65.7°→52.8°（−19.6%）；矢量角度：所有处理>45°；矢量长度与角度：R² = 0.840，p < 0.001；矢量长度与CUE：R² = 0.502，p < 0.001；矢量角度与CUE：R² = 0.538，p < 0.001；矢量长度与AGB：R² = 0.270，p < 0.001；与BGB：R² = 0.666，p < 0.001；矢量角度与AGB：R² = 0.406，p < 0.001；与BGB：R² = 0.681，p < 0.001
3.4 MI驱动土壤细菌群落向富营养型生活史策略转变
α多样性指数（Shannon和Richness）处理效应有限，仅在I-60（两个指数）和II-90（仅Richness）出现显著差异；I-30、I-90或III-90未检测到处理效应（图5a和b）。
相比之下，富营养型-寡营养型比值在I-90、II-90和III-90受显著影响，但在I-30或I-60不受影响（图5c）。总体而言，显著微生物群落变化在60天后出现，并在后期更为明显，表明时间依赖性重构而非即时响应。通过PERMANOVA和ANOSIM的β多样性分析证实所有五个采样时间均存在显著处理诱导的群落分化（图6）。
ANOSIM R值范围为0.391至0.704，表明处理间存在实质性分离。PCoA解释了各阶段总变异的36.85–49.05%。处理解释力（R²）从0.482（I-30）增至0.536（III-90），表明群落分化随重复接种逐步增强，在第三年生长高峰期达到最大分化。在门水平，Proteobacteria、Actinobacteriota、Acidobacteriota、Chloroflexi、Firmicutes、Methylomirabilota和Bacteroidota在所有处理和采样时间的细菌群落中占主导，合计相对丰度>85%（图S4；表S15）。Proteobacteria和Actinobacteriota始终是最丰富的两个门，合计丰度超过60%，并呈现时间和处理依赖性变化。虽然退化或接种下优势门未发生根本改变，表明核心类群高度稳定，但其比例丰度在处理和阶段间显著偏移，表明退化和微生物接种重塑了群落结构而非替换关键成员。在属水平，重复年度MI施用诱导MD和SD土壤中Pseudomonas相对丰度逐渐、时间依赖性增加，从I-60起出现显著富集，在第三年稳态（III-90）达到最大（表S16）。在III-90，Pseudomonas丰度从未接种MD的1.90%增至MD + MI的2.19%，从SD的2.25%增至SD + MI的3.07%。值得注意的是，在五个采样时间检测到的所有优势属中，仅Pseudomonas在MI施用下表现出一致、持续的增加趋势，为引入菌群在三个年度接种周期后成功定殖根际提供了直接田间证据。
【数据】α多样性：仅I-60（Shannon和Richness）和II-90（Richness）显著；富营养型-寡营养型比值：I-90、II-90、III-90显著；ANOSIM R：0.391–0.704；PCoA解释变异：36.85–49.05%；处理R²：0.482（I-30）→0.536（III-90）；优势门相对丰度：>85%；Proteobacteria + Actinobacteriota：>60%；Pseudomonas丰度：MD 1.90%→MD + MI 2.19%；SD 2.25%→SD + MI 3.07%
3.5 退化依赖的调控路径主导植物生物量恢复
采用层次分割法解析土壤性质、养分库、微生物生物量和资源限制对各处理地上（AGB）和地下（BGB）生物量恢复的相对贡献。层次分割……
【数据】原文未详述具体数值
讨论与解读
三年田间试验表明，本地假单胞菌接种剂能够改变退化高寒草甸土壤养分有效性、增加微生物生物量并改变胞外酶活性，最终提高地上和地下植物生物量。然而，修复机制并不统一，强烈依赖于土壤初始退化状态。结果仅部分支持假设H1：MI确实缓解了磷限制，但加剧而非降低了碳限制。假设H2得到支持，微生物生物量增加和酶分配变化是关键路径，但其相对重要性在MD和SD土壤间存在差异。16S rRNA基因扩增子测序数据和化学计量CUE估算的整合，使研究超越了相关性推断，为观察到的化学计量响应背后的微生物群落变化提供了直接证据。
MI施用虽改善土壤肥力和植物生长，却持续增加矢量长度（微生物碳限制指标）。简单化学计量解释会认为MI使微生物群落碳更受限，但这与观察不一致：MI施用后SOC和DOC均增加。这一矛盾反映了微生物生活史策略的转变——富营养型细菌比例上升，其高代谢活性导致单位生物量的碳需求增加，从而在碳供给改善的同时表现出更强的碳限制信号。这一发现提示，退化高寒草甸的微生物修复不能仅以养分供给为指标，还需关注群落代谢策略的定向演替。
编译者解读：该研究将本地假单胞菌接种剂与酶化学计量、CUE和群落生活史策略耦合分析，方法框架清晰，田间三年验证增强了结论的实践可信度。其价值在于揭示了退化梯度上“磷限制缓解—碳限制加剧—富营养型菌群扩张”的连锁机制，为高寒草甸精准修复提供了可操作的生物资源方案。局限在于菌株定殖的长期稳定性、非目标微生物群落的潜在扰动以及成本效益尚未评估，距离规模化牧场应用仍有距离。
参考来源
期刊：Journal of Environmental Management, 2026
DOI: 10.1016/j.jenvman.2026.130880
DOI: 10.1016/j.jenvman.2026.130880