土壤越老越贫瘠，微生物反而会重新组队来维持养分循环

土壤在很长时间里越变越老、越变越酸、越变越缺养分，很多微生物确实会减少，多样性也会下降。但留下来的微生物不是各自孤立地硬撑，而是重新组织关系。细菌的部分养分动员功能变弱，真菌变得更重要，尤其是能和植物形成菌根关系、能分解难降解植物材料、能释放有机磷和有机氮的真菌。病毒数量和多样性整体也下降，但有些噬菌体可能携带帮助宿主应对缺磷的基因。也就是说，老化贫瘠土壤里的微生物世界看起来成员更少，却可能连接更紧、分工更集中，靠真菌和跨域互作维持基本养分循环。

这篇文章研究的是生态系统退化演替对土壤微生物群落的影响。这里的退化演替不是短期污染或人为破坏，而是指一个生态系统在几千年到几百万年里长期缺少新的扰动补给，土壤不断风化、发育，养分慢慢流失，最后变得越来越贫瘠。文章强调，磷会被淋失，或变成植物和微生物难以利用的闭蓄态、难降解形态；氮也会因为磷限制抑制固氮，以及被植物凋落物和腐殖质中的多酚复合物固定而变得不易利用。结果是植物生物量下降，净初级生产力降低，凋落物分解变慢，土壤团聚体稳定性、孔隙度和保水能力也下降。

长期贫瘠化不仅改变植物，也会深刻改变地下微生物食物网。退化演替环境中的植物往往更节约资源，并把更多能量投入结构和化学防御，例如产生更多酚类化合物，或在组织中积累硅。这会改变进入土壤的有机质质量，进一步影响碳、氮、磷循环。以前在澳大利亚 Jurien Bay、新西兰 Haast dune 和 Franz Josef Glacier 等长期年代序列中的研究已经发现，随着磷、pH 和养分下降，细菌和真菌群落会发生变化，有些系统还会出现从细菌主导向真菌主导的转变。本文选择美国加州 Mendocino 的 Ecological Staircase，就是因为这里的海阶由相同母质形成，但土壤年龄、酸度和肥力形成了清楚梯度，适合用空间代替时间来推断长期生态系统变化。

在方法上，作者沿着 Mendocino Ecological Staircase 的海阶年代序列采集土壤。研究地点位于美国加州 Fort Bragg 附近的 Jug Handle State Natural Reserve。这个年代序列包含约 10 万到 50 万年 的海阶，由海蚀和构造抬升形成。作者重点比较年轻海阶，也就是演替建立阶段的 terrace 1，以及处于退化演替阶段的 terrace 3 和 terrace 4，分别称为 young、mid-age 和 old。每个海阶采集 9 个独立土壤样品，总共 27 个样品。采样前去除表层有机层，只采矿质土壤；年轻海阶采 0–20 cm 矿质层，中年和老年海阶采 0–10 cm 较薄的矿质层。

作者测定了土壤 pH、有机质、总磷、总有机碳和总氮等性质，并用较大质量的土壤样品提取 DNA，以尽量覆盖原生生物等低丰度真核微生物。分子方法包括扩增子测序、qPCR、宏基因组测序、基因中心分析、宏基因组组装基因组分析和病毒识别。扩增子测序分别针对原核生物、真菌和原生生物；qPCR 用来估算这些类群的基因拷贝数；宏基因组测序每个海阶做 3 个重复，每个样品约 20 Gb；随后用 MEGAHIT 组装、VEBA 注释、MetaBAT2 和 CONCOCT 分箱，识别原核、微真核和病毒 MAG，并用 VirSorter2、CheckV、geNomad 等工具分析 DNA 病毒。网络分析使用随机矩阵理论构建共现网络，并进一步计算相对模块化、鲁棒性、脆弱性和 cohesion。

主要结果显示，退化演替显著改变了土壤环境和所有主要微生物类群。年轻海阶的 pH、总磷、土壤有机质、总有机碳和总氮最高；中年和老年海阶更酸、更缺养分。微生物 alpha 多样性随海阶年龄显著下降，原核生物、真菌和原生生物的 Shannon 指数都显著降低。原核生物丰度从年轻到中年海阶下降约 74.8%，从中年到老年又下降约 70.9%；原生生物丰度也从年轻到中年下降约 67.1%，从中年到老年下降约 64.6%。真菌的表现不一样，它们在中年海阶达到峰值，比年轻海阶高 117%，而中年和老年之间没有显著差异。

群落组成也被海阶年龄强烈塑造。PERMANOVA 显示，海阶年龄分别解释原核生物、原生生物和真菌群落变化的 95.2%、82.0% 和 80.8%，合并所有微生物类群后，海阶年龄解释总群落组成变化的 86%。把土壤有机质、总氮、总有机碳、总磷、植物 β 多样性和植物 Shannon 指数加入模型后，解释度提高到 91.4%。Mantel 分析显示，总微生物群落组成与 pH、植物群落组成、总氮和总磷关系最强，其中 pH 的相关性最高，r = 0.86，p = 0.001。这说明长期土壤发育不是一个单独因素，而是把酸化、贫养分、植物群落变化和微生物重组连在一起。

随着退化演替加深，微生物网络的 相对模块化 RM 普遍升高。真菌网络 RM 从 0.66 增加到 0.88，原生生物从 0.20 增加到 0.74，原核生物呈 U 型变化，在中年海阶下降后，到老年海阶回升到 0.65；合并网络则从 0.64 增加到 1.33。模块化升高意味着群落被分成更清楚的子群或功能模块，通常被认为有助于稳定性，因为扰动不容易一下子传遍整个网络。作者进一步做随机节点移除和靶向节点移除，发现老年海阶的网络鲁棒性更高，单个节点删除造成的网络脆弱性更低。

真菌在这个重组网络中扮演关键稳定者。作者比较移除真菌节点和移除原生生物节点后网络鲁棒性的变化，发现真菌节点被移除对鲁棒性的影响显著大于原生生物节点。这说明真菌不是老化土壤中普通增加的一类成员，而是维持跨域网络连接的重要结构成分。cohesion 分析也支持这一点：腐生真菌与正 cohesion 显著相关，r = 0.61，p < 0.01；食细菌型和杂食型原生生物也与正 cohesion 有关，但真菌和合并网络在老年海阶的负 cohesion 降低，提示真菌可能推动群落从强竞争向更有协作性的组织方式转变。

群落组装机制也随退化演替变化。作者用 iCAMP 分析同质选择、扩散限制、漂变等过程。结果显示，同质选择、扩散限制和漂变是影响跨域群落组装的主要机制。同质选择显著塑造中年海阶的原核生物群落和老年海阶的真菌群落，也影响年轻海阶的原生生物群落。扩散限制在较老海阶的原核生物和真菌中更明显，在中年海阶的原生生物中也显著。作者认为，老化贫瘠土壤保水能力较低，可能限制微生物移动和定殖，从而加强扩散限制。漂变在年轻海阶的原核生物和真菌组装中更重要，但在退化演替阶段的原生生物中影响更大。

功能基因结果显示，退化演替并不是简单让所有功能都变弱，而是改变了谁来承担养分循环。原核生物方面，固氮酶基因在老年海阶显著更丰富，12 个 MAG 携带固氮酶基因，其中 Beijerinckia.4764、Rhodomicrobium.2215 和 Rhodoblastus.2093 在老年海阶丰度较高。相反，原核生物来源的脲酶、氨单加氧酶、漆酶和磷酸酶基因在年轻海阶更丰富。作者本来预计老年低磷土壤中原核磷酸酶基因会更多，但结果相反，说明在严重退化演替土壤里，原核生物可能不再是磷动员的主要承担者。

真菌功能基因呈现相反趋势。真菌来源的纤维素酶和木聚糖酶基因在老年海阶显著增加，说明真菌在老化酸性土壤中可能更重要地参与植物聚合物分解。真菌来源的植酸酶、酸性和碱性磷酸酶、脲酶基因也在中年和老年海阶显著高于年轻海阶。这意味着，随着土壤变老、变酸、缺氮缺磷，养分获取策略逐渐从原核生物主导转向真菌主导。真菌不仅可能释放有机磷，也可能通过脲酶从有机氮或尿素相关底物中释放氮，还能降解富含木质纤维素和芳香族化合物的植物残体。

这些真菌功能变化与土壤化学高度相关。真菌磷酸酶、植酸酶和脲酶基因与总磷、总氮和 pH 呈负相关，相关系数约为 –0.7 到 –0.9，说明它们在更酸、更贫养分的土壤中更丰富；同类原核基因则与这些土壤指标呈正相关，r = 0.7–0.96。结构方程模型进一步显示，总氮下降主要通过改变微生物群落组成产生间接影响，而总磷和 pH 则是直接限制功能基因潜力的重要环境因素。模型解释了群落组成 95% 的变异和功能基因指数 82%–99% 的变异。作者也提醒，这个模型应被视为探索性的，因为宏基因组重复数有限，而且基因丰度只能代表功能潜力，不等于真实过程速率。

宏基因组组装还帮助作者识别了可能参与养分动员的真菌成员。研究共重建了 15 个真菌 MAG，完整度为 33%–75%，污染度低于 10%。这些真菌 MAG 在中年海阶更丰富，与 qPCR 显示真菌丰度在中年海阶较高一致。其中 Hyaloscypha.39 同时含有碱性和酸性磷酸酶、脲酶、纤维素酶和木聚糖酶基因。其他多个 Hyaloscypha、Endocarpon、Pezizales 和 Helotiales 相关 MAG 也含有磷酸酶、植酸酶或脲酶等基因，并与杜鹃花科植物或针叶树的菌根关联有关。

病毒结果为这篇文章增加了一个常被忽略的维度。作者从宏基因组中识别出 1636 个 DNA 病毒，其中 1621 个属于 Caudoviricetes 噬菌体，另有 15 个属于 Nucleocytoviricota 巨型病毒。DNA 病毒的 Shannon 多样性从年轻到较老海阶显著下降，病毒群落组成也按海阶明显分化。Caudoviricetes 在年轻海阶更丰富，这与年轻海阶更高的原核生物多样性和丰度一致；而 Nucleocytoviricota，包括感染微藻的 Algavirales 和其他巨型病毒类群，在较老海阶更常见，与老年海阶光合原生生物增加相一致。更有意思的是，虽然总体噬菌体多样性下降，但一些与磷代谢有关的病毒基因在贫磷环境中可能变得更重要。作者发现一个 Caudoviricetes 噬菌体携带 phoD，该基因编码碱性磷酸酶，蛋白结构建模支持这个注释。另一个在老年海阶富集的 Caudoviricetes.10593003 编码 PhoH，这是与磷饥饿响应相关的蛋白；iPHoP 预测它的宿主是 Gemmataceae.1082，置信度为 93.9%。这个宿主 MAG 本身含有 PafA、PhoH 和 PhoU 等磷获取和稳态相关基因。噬菌体和宿主的覆盖度在样品内相关显著，r = 0.73，p = 0.026，按海阶均值相关更强，r = 1.00，p = 0.007。作者没有把这说成已经证实的机制，而是谨慎认为这支持一种可能的噬菌体—宿主关联，可能影响低磷土壤中的磷代谢，仍需实验验证。

结论

生态系统退化演替会降低许多微生物类群的多样性和丰度，但不会简单削弱微生物系统。相反，随着土壤变老、变酸、缺氮缺磷，微生物群落会发生跨域重组。原核生物、真菌、原生生物和病毒之间的网络更模块化、更鲁棒，真菌成为维持网络连接和稳定性的关键类群。功能上，贫瘠老化土壤中养分获取潜力从细菌侧更多转向真菌侧，尤其是真菌携带更多与磷动员、氮获取、纤维素和木质纤维素降解、多酚和芳香族化合物处理有关的基因。固氮基因、光合原生生物和携带磷代谢辅助基因的噬菌体，也可能共同支持贫瘠土壤中的养分循环潜力。

这篇文章提醒我们判断生态系统退化不能只看微生物多样性下降，也要看群落是否发生了功能重组和网络重连。某些系统里，多样性降低并不一定立刻意味着功能完全崩溃，因为剩余成员可能通过更强的跨域互作维持功能潜力。但这种稳定性也可能是脆弱的，因为它更依赖少数关键类群，尤其是真菌。未来研究如果能把长期年代序列、宏基因组、宏转录组、稳定同位素示踪、酶活性测定、真菌分离培养、噬菌体—宿主实验和网络扰动实验结合起来，就能更清楚地区分有功能基因与真实发生养分循环之间的关系，也能帮助我们理解和修复贫瘠化、酸化和长期退化的土壤生态系统。