《Geoderma》：Litter-type-specific succession of microeukaryotic communities and their associations with litter decomposition
编辑推荐：
为解决微核生物（真菌和原生动物）在凋落物分解过程中的演替规律及其对分解速率的调控机制尚不明确的问题，研究人员利用长读长纳米孔测序技术，对六种不同C/N比、分属四个植物科的凋落物进行为期一年的野外分解实验，系统解析了微核生物群落组成与功能的演替动态及其对凋落物质量损失的预测作用。研究发现微核生物群落的演替具有强烈的凋落物类型特异性，且原生动物相关指标随时间推移对分解速率的预测重要性增加，强调了整合真菌和原生动物动态对理解陆地生态系统碳循环的关键意义。
在森林和农田的土壤表层，落叶、枯枝和其他植物残体构成了凋落物层。这些看似不起眼的“垃圾”实则至关重要，它们的分解过程是全球碳循环和养分周转的核心引擎。传统观点认为，真菌是驱动这一过程的绝对主力军。然而，土壤中还存在着一类数量庞大、种类繁多的微小真核生物——原生动物，它们以捕食细菌和真菌为生。虽然已知真菌与原生动物之间存在互动，但两者如何协同演替，又如何共同影响不同类型凋落物的分解命运，仍然是一个巨大的“黑箱”。理解这些微观生命世界的动态，对于预测生态系统如何响应环境变化、调控土壤碳储量具有深远意义。
为了揭开这一谜团，一个由Yuxin Wang、Geert Smant、Stefan J.S. van de Ruitenbeek和Stefan Geisen组成的研究团队在《Geoderma》上发表了一项研究。他们不再局限于单一凋落物类型，而是选择了六种覆盖作物（大麦、萝卜、三叶草、羊茅、燕麦、万寿菊）的凋落物，这些材料横跨了宽广的碳氮比（C/N）梯度（6–47），并来自禾本科、十字花科、豆科和菊科等多个植物家族，旨在获得更具普适性的发现。研究人员采用尼龙网袋法将这些凋落物埋入田间土壤，在分解开始后的30天（第一阶段）、210天（第二阶段）和360天（第三阶段）分别取样。通过对凋落物残留物进行DNA提取，并利用能覆盖18S、ITS、5.8S和28S rRNA区域（约4.5 kb）的通用引物进行扩增，他们创新性地采用了牛津纳米孔长读长测序技术，以获取微核生物（包括真菌和原生动物）从分类到功能的完整信息。通过计算α-多样性（如Chao1丰富度指数、Shannon多样性指数）和β-多样性（基于Bray-Curtis距离的主坐标分析），并结合随机森林模型分析，系统探究了微核生物群落的演替规律及其与凋落物分解速率的关系。
主要技术方法概述：
本研究在荷兰瓦赫宁根大学的实验农场进行，使用了六种覆盖作物的凋落物。核心方法包括：1) 野外凋落物分解实验设计，将不同凋落物装入尼龙网袋埋于土中，分三个时间点取样；2) 利用长读长纳米孔测序（SSU–ITS–LSU，~4.5 kb）对凋落物样品中的微核生物核糖体RNA操纵子进行测序；3) 使用EukFunk数据库对获得的序列进行功能注释，将其划分为共生营养型、腐生营养型和捕食者等功能群；4) 运用多元统计（如PERMANOVA）和随机森林机器学习模型，分析群落动态并识别预测凋落物质量损失的关键因子。
3.1. 微核生物群落α-多样性的演替
通过分析Chao1指数和Shannon指数，研究人员发现微核生物群落的α-多样性呈现出特定的时间动态。在所有凋落物类型中，分类单元丰富度（Chao1）遵循单峰轨迹，在分解中期达到峰值，而Shannon多样性则随时间持续增加。更重要的是，这些动态具有显著的凋落物类型特异性。例如，在分解初期（第一阶段），具有较高C/N比的燕麦比C/N比较低的大麦表现出更高的α-多样性。而含有抗菌次生代谢物（如萝卜中的硫代葡萄糖苷）的凋落物，其α-多样性在后期阶段（第二、三阶段）增长更为陡峭。
3.2. 微核生物群落组成的演替
群落组成（β-多样性）分析表明，凋落物类型是塑造微核生物群落结构的主要因素，而时间进程则降低了同一凋落物类型内的群落差异。研究人员识别出了阶段特异性指示类群和与特定凋落物持久关联的类群。例如，真菌中的毛霉目和肉座菌目在早期阶段（第一阶段）在所有凋落物中相对丰度均较高；而到了中期（第二阶段），能够降解木质素等顽固基质的粪壳菌目则在燕麦和万寿菊中更为丰富。对于原生动物，萝卜和大麦在早期阶段显示了较低的滑行虫目丰度，而羊茅在中期则显示出较高的尾丝虫目、诺兰虫目和真变形虫目丰度。
3.3. 微核生物群落功能群的演替
功能分析揭示了微核生物功能群随时间的演替规律。总体上，捕食者的相对丰度随时间先增加后下降，而共生营养型生物（特别是植物病原菌）的相对丰度则持续下降。然而，腐生营养型生物的相对丰度未表现出统一的跨凋落物类型的时间变化趋势，其动态更多受到凋落物类型的影响。例如，万寿菊在整个分解过程中都保持了较高的腐生营养型生物丰度，同时具有较低的植物病原菌和捕食者丰度，这可能与其富含抗菌次生代谢物和高C/N比有关。
3.4. 凋落物分解速率的预测因子
随机森林模型分析揭示了预测凋落物质量损失的关键因子。整合所有分解阶段的数据，模型可以解释80%的凋落物质量损失变异。最重要的四个预测因子是原生动物属Paracercomonas、原生动物与真菌的比率、凋落物初始C/N比以及原生动物属Cercomonas。分阶段分析显示，真菌相关指标在分解早期（第一、二阶段）对解释质量损失变异更为重要（分别占52%和41%），而原生动物相关指标的重要性在后期（第二、三阶段）显著增加（分别占37%和47%），在分解末期成为最强的预测因子，解释了47%的变异。
研究结论与意义
这项为期一年的研究清晰地揭示，在凋落物分解过程中，微核生物（真菌和原生动物）群落经历着强烈的、凋落物类型特异性的演替。这种演替不仅体现在分类学组成的动态变化上，也体现在功能群的更替中。研究证实了初始假设的部分内容，如α-多样性中的丰富度呈单峰变化，以及群落组成和功能（向捕食者增加、共生营养型减少）随时间的转变。尤为关键的发现是，原生动物在分解过程中的作用随时间推移而凸显。随机森林模型表明，原生动物相关指标在分解后期对预测凋落物质量损失变得至关重要。
这项研究的意义在于，它将真菌和原生动物这两个常常被分开研究的微核生物类群置于同一框架下，并置于多种凋落物类型的真实生态梯度中进行考察。结果强调了凋落物内在性状（如C/N比、植物家族特异性次生代谢物）通过“自下而上”的调控方式，深刻影响着分解者群落的组建与演替轨迹。这些群落的动态反过来又与分解速率紧密关联。因此，未来要更准确地预测碳循环和生态系统功能，必须采用更全面的土壤食物网视角，整合更多生物类群，并对凋落物性状进行多维度的动态表征。同时，需要通过直接的操控性实验（如原生动物的排除或添加）来验证其在凋落物分解中的因果关系。这项工作为理解地下生态过程及其对全球变化的反馈提供了新的、更整合的见解。
