英文题目:Ecological memory of recurrent drought modifies soil processes via changes in soil microbial community
中文题目:反复干旱的生态记忆通过改变土壤微生物群落来调节土壤过程
期刊名称:Nature Communications
影响因子:15.7
作者单位:奥地利维也纳大学微生物学与环境系统科学中心
DOI:https://doi.org/10.1038/s41467-021-25675-4
通过在奥地利高山草甸开展长达10年的野外控制实验,系统探究了反复干旱事件诱导土壤微生物群落形成生态记忆及其功能缓冲效应的机制。研究表明,与单次干旱相比,连续10年的反复干旱显著改变了土壤细菌、古菌、真菌和丛枝菌根真菌的群落组成(Bray-Curtis相异指数显著升高),并重塑了微生物量C:N:P化学计量特征,使其在后续干旱中趋于稳定;长期干旱促使微生物群落优化养分获取策略(显著提高C:P和N:P酶活比),并富集与耐旱性相关的功能基因(如渗透保护物质合成、胞外多糖代谢、活性氧清除等),从而在重新湿润后表现出更快的生长和呼吸恢复能力。最关键的是,土壤多功能性指数在反复干旱处理下显著高于单次干旱组,表明历史干旱经历通过创造“生物遗产”,使微生物群落“记住”胁迫事件并预先适应,有效缓冲了后续干旱对土壤碳氮磷硫循环和生态系统功能的冲击。该研究首次在原位条件下证实,生态记忆的形成是土壤微生物群落应对气候变化、增强生态系统恢复力的关键机制。
生态记忆形成:反复干旱对土壤微生物量化学计量特征的影响
图1展示了干旱对土壤微生物量及其化学计量学特征的影响,揭示了反复干旱处理下的生态记忆效应。研究发现,尽管微生物量碳(MBC)未受显著影响,但微生物量氮(MBN)和微生物量磷(MBP)在不同处理间存在显著差异,且呈现一致模式:单次干旱处理(1年)显著偏离对照组,而连续10年反复干旱处理与对照组数值相似。微生物量的C:N和C:P比同样在单次干旱组显著高于对照组和10年处理组,而N:P比未受显著影响。这一结果证实,长期反复干旱通过改变微生物群落组成,使其形成“生态记忆”,从而在后续干旱中维持更稳定的生物量化学计量特征,表现出对干旱胁迫的缓冲和适应能力。
图1.干旱对微生物量C、N和P的影响。颜色表示土壤处理(对照=灰色,1年=浅蓝色,10年=蓝色)。上图显示微生物量碳(MBC)、微生物量氮(MBN)和微生物量磷(MBP)。MBC未受干旱显著影响(F值:2.02,p=0.18),而MBN和MBP受到显著影响(F值:5.7,p=0.04,n=4;F值:4.44,p=0.06,n=4)。这三个数值均呈现相同模式,即对照和10年处理具有相似值,而1年处理则不同。下图显示微生物量的化学计量比。C:N和C:P比值显示出相似结果,其中1年处理与对照和10年处理均存在显著差异(分别为F值:5.72,p=0.0249,n=4;F值:8.06,p=0.0098,n=4)。N:P比值未受干旱显著影响,10年处理值介于1年处理和对照之间。字母表示双侧Tukey HSD事后检验结果(字母旁的圆点表示p值在0.05阈值上)。箱线图中心线代表中位数,箱体边界代表上下四分位数,须状线代表1.5倍四分位距,而单独的点代表异常值。样本量n代表生物学独立样本。
通过扩增子测序揭示了长期反复干旱对土壤微生物群落结构的深远影响(图2)。非度量多维尺度分析(NMDS)显示,10年反复干旱处理在细菌/古菌、真菌和丛枝菌根真菌(AMF)群落中均与对照组及单次干旱处理形成显著分离(PERMANOVA检验p<0.05),而单次干旱处理与对照组无明显差异。差异丰度分析进一步表明,长期干旱导致大量微生物类群发生显著变化:酸杆菌门、拟杆菌门、担子菌门和多孢目的相对丰度降低,而放线菌门、子囊菌门和球囊菌目的丰度增加;相比之下,单次干旱仅在真菌担子菌门中检测到显著差异。值得注意的是,各处理间的α-多样性指数(丰富度、Shannon指数等)无显著变化,说明群落变异主要由类群替换而非多样性丧失驱动。这一结果从分子层面证实,长期反复干旱通过选择性富集或淘汰特定微生物类群,在群落水平上形成了持久的“生物遗产”,这种结构重组是土壤生态记忆形成的关键机制,使微生物群落能够“记住”历史胁迫并在后续干旱中表现出更强的功能稳定性。
图2.通过扩增子测序评估干旱对微生物群落组成的影响。颜色表示土壤处理(对照=灰色,1年=浅蓝色,10年=蓝色)。上方为NMDS图,代表古菌/细菌、真菌和AMF(丛枝菌根真菌)。所有NMDS图均显示相似模式,干旱有显著影响(分别为p=0.014,p=0.01,p=0.015;PERMANOVA结果,n=4),且10年处理与1年处理和对照明显分离。下图显示对照vs 1年(上排)和对照vs 10年(下排)的差异丰度分析结果。红色圆圈代表与对照相比具有显著差异丰度的分类单元,灰色代表无显著差异。与NMDS图类似,对照vs 1年结果显示差异很小(仅在真菌界的担子菌门),而对照vs 10年在所有水平上均显示出强烈差异(受影响最大的类群包括:酸杆菌门、放线菌门、变形菌门、子囊菌门、担子菌门、多孢目和球囊菌目)。
图3揭示了长期反复干旱通过重塑酶活性策略形成生态记忆的功能机制。研究发现,C循环相关酶在10年处理组中维持对照水平,而单次干旱组则显著降低(p=0.016),表明历史干旱使微生物群落稳定了碳分解功能;N循环酶呈现分化响应,单次干旱组所有酶活性均下降,而10年组在酪氨酸氨肽酶(TAP)恢复至对照水平(p<0.001),显示特定功能的选择性恢复;P/S循环酶(磷酸酶和硫酸酯酶)则对两类干旱均敏感,活性显著降低(p<0.01)。更关键的是,长期反复干旱显著提高了微生物群落的C:P和N:P获取比(p=0.001和p=0.008),而单次干旱无此效应,说明群落通过优先通过碳氮获取策略来适应胁迫。这一结果从功能层面证实,生态记忆不仅体现在群落结构改变,更通过系统性优化养分获取策略和酶促反应模式,在后续干旱中缓冲土壤多功能性损失,增强生态系统气候抗性。

图3.干旱对土壤潜在酶活性的影响。颜色表示土壤处理(对照=灰色,1年=浅蓝色,10年=蓝色)。上图显示所有测量酶的潜在酶活性,按C、N以及P和S循环相关酶分组。C循环酶均显示相似趋势,其中10年处理与对照值相似,而1年处理始终不同(仅在木糖苷酶中达到统计显著:F=6.68,p=0.016,n=4;而在β-葡萄糖苷酶中为部分显著:F=3.2,p=0.089,n=4)。N循环酶显示10年处理的不同响应,其中两种酶在两种干旱处理中均显示下降(内切几丁质酶:F=4.9,p=0.036,n=4;亮氨酸氨肽酶:F=14.58,p=0.0015,n=4),而另外两种仅1年处理显示下降值,10年处理与对照值相似(仅在酪氨酸氨肽酶中显著:F=18.79,p<0.001,n=4)。P和S循环酶显示相同趋势,即两种干旱处理均与对照有显著差异(磷酸酶:F=6.53,p=0.018,n=4;硫酸酯酶:F=14.75,p=0.0014,n=4)。星号代表Tukey HSD事后检验的显著差异,圆点代表部分显著值(p<0.1)。左下图显示2017年潜在酶活性的NMDS分析,所有处理间存在明显分离(R²=0.58,p=0.0003,n=4,stress=0.13)。右下图代表获取策略(潜在酶活性的化学计量比),所有比值均显示相似趋势,尽管仅在C:P和N:P中达到统计显著(分别为F=14.8,p=0.001,n=4;F=8.6,p=0.008,n=4)。字母表示双侧Tukey HSD事后检验结果。箱线图中心线代表中位数,箱体边界代表上下四分位数,须状线代表1.5倍四分位距,而单独的点代表异常值。样本量n代表生物学独立样本。图4揭示了反复干旱对土壤微生物生理响应的塑造作用,体现了生态记忆在群落功能层面的表达。研究显示,长期反复干旱(10年)显著提高了微生物的质量特异性呼吸速率(p=0.006)和生长速率(p=0.0015),而单次干旱处理仅引起微弱响应;尽管碳利用效率(CUE)未受显著影响,但10年处理组呈现出较低的CUE趋势,表明微生物群落将更多资源分配给快速代谢和生长而非碳储存。这种生理特征的适应性转变,说明历史干旱经历筛选出了能够在重新湿润后迅速激活的耐受性菌群,形成了可遗传的“记忆”功能,使其在再次遭遇干旱时能更快恢复生理活性,从而增强生态系统的整体恢复力。
图4.土壤再湿润后的生理参数。颜色表示土壤处理(对照=灰色,1年=浅蓝色,10年=蓝色)。干旱对单位质量特异性呼吸速率和生长速率有显著影响(分别为F=9.67,p=0.006,n=4和F=14.55,p=0.0015,n=4),但对CUE无显著影响。10年干旱在单位质量特异性呼吸速率上与1年干旱存在显著差异。字母表示双侧Tukey HSD事后检验结果。箱线图中心线代表中位数,箱体边界代表上下四分位数,须状线代表1.5倍四分位距,而单独的点代表异常值。样本量n代表生物学独立样本。
通过PICRUSt2功能预测分析揭示了反复干旱对土壤细菌群落功能潜力的深层塑造作用(图5)。与单次干旱相比,10年处理组在氨基酸代谢、膜转运、辅因子代谢等主要功能类群的丰度上发生显著转变,且这些功能变化涉及多个基因同时富集和缺失的复杂模式。更关键的是,长期反复干旱显著富集了一系列与耐旱性直接相关的功能基因,包括渗透保护物质合成(如甜菜碱、脯氨酸)、胞外多糖(EPS)代谢、肽聚糖生物合成、磷脂重塑及活性氧(ROS)清除等相关基因,表明群落通过系统性强化细胞保护机制来适应水分胁迫。同时,不同菌门成员在同一功能类群中同时呈现正负贡献,说明功能适应发生在更精细的分类水平上,而非简单的类群替代。这一结果从功能基因层面证实,生态记忆的形成伴随着微生物群落的重构——将资源从基础代谢转向环境胁迫应对,这种功能层面的预适应是土壤生态系统抵御未来干旱的分子基础。
图5.土壤微生物群落的潜在功能能力和获取策略。三个图显示PICRUSt2对1年干旱与10年干旱的比较结果。揭示了与干旱历史响应相关的微生物群落许多功能性状的潜在富集或消耗(分别为正和负的Log2倍数变化)。左图显示1年和10年处理中细菌群落内KEGG直系同源基因的平均比例。颜色表示土壤处理(1年=浅蓝色,10年=蓝色)。红色星号表示整个功能组的显著差异丰度。中间图代表使用DESeq2富集(黑色圆圈)或消耗(紫色圆圈)的功能数量(分别为正和负的Log2倍数变化)。右图代表对富集或消耗功能有贡献的微生物类群(气泡大小表示每个门内ASV在各功能组中的频率高低)(分别为正或负)。各门的颜色在图例中描述。
生态记忆的量化:反复干旱增强群落差异并缓冲土壤多功能性
通过构建Bray-Curtis相异指数和土壤多功能性指数,从量化和整合的角度验证了生态记忆的形成效应(图6)。左侧图片显示,反复干旱处理组(10年)与对照组之间的微生物群落差异(包括古菌/细菌、真菌和丛枝菌根真菌)显著高于单次干旱组(p<0.01),表明长期干旱历史在群落水平上形成了可累积的“生物遗产”,导致群落结构持续偏离原始状态。右侧图片进一步揭示,单次干旱事件导致土壤多功能性指数显著下降,而经过10年反复干旱的群落,其多功能性指数反而更接近对照组水平(p=0.027),说明历史干旱经历通过重塑微生物群落功能,有效缓冲了后续干旱事件对土壤综合功能的冲击。这两个指数共同证明,生态记忆不仅表现为群落结构的定向转变,更关键的是通过优化功能配置,在生态系统层面增强了气候胁迫下的功能稳定性与恢复力。
图6.干旱对微生物群落和土壤功能的影响。颜色表示土壤处理(对照=灰色,1年=浅蓝色,10年=蓝色)。左侧图片显示箱线图,代表干旱处理与对照之间的Bray-Curtis相异性指数,分别针对古菌/细菌(F=44.969,p=0.0005,n=4)、真菌(F=32.86,p=0.012,n=4)和AMF(F=23.002,p=0.003,n=4)。右侧图片显示箱线图,代表土壤多功能性指数(通过广泛的土壤参数计算,包括酶活性、土壤C和养分库、微生物量C、N和P以及根系生物量生产力)按处理分组(F=5.514,p=0.027,n=4)。字母表示双侧Tukey HSD事后检验结果。箱线图中心线代表中位数,箱体边界代表上下四分位数,须状线代表1.5倍四分位距,而单独的点代表异常值。
该研究通过在奥地利高山草甸开展的长期野外实验,发现10年反复干旱事件使土壤微生物群落(细菌、古菌、真菌和丛枝菌根真菌)的组成发生显著变化,与单次干旱形成鲜明对比,并富集了与耐旱相关的功能基因(如渗透保护剂合成、胞外多糖代谢等),从而形成了生态记忆。这种生物遗留效应使土壤多功能性(涵盖酶活性、养分循环和微生物生理过程)在后续干旱中表现出更强的抵抗力,其功能变化更接近未受干旱的对照水平。反复干旱后的微生物群落具有更快的生长和呼吸速率,但碳利用效率降低,表明其能更迅速地响应水分恢复。研究证实,长期反复干旱通过重塑土壤微生物群落结构和功能策略,缓冲了干旱对土壤过程的负面影响,这对预测气候变化下日益频繁的极端干旱事件及其生态后果具有重要意义。参考文献:Canarini A , Schmidt H , Fuchslueger L ,et al.Ecological memory of recurrent drought modifies soil processes via changes in soil microbial community[J].Nature Communications, 2021, 12.DOI:10.1038/s41467-021-25675-4.