红树林、盐沼和海草床把大量碳锁在泥里，真正干活的是一群肉眼看不见的微生物。环境变化如果换掉一批微生物，碳封存会不会跟着崩？科学家沿着中国海岸取泥，给这些微生物做了一次大规模“户口普查”。

红树林、盐沼和海草床被称为蓝碳生态系统，面积不算大，单位面积却能长期储存大量碳。微生物在这里分解有机质、驱动氮和硫的氧化还原，也参与碳固定。前人通过区域采样已经知道，这些湿地里的微生物种类特殊，耐低氧、耐盐、能处理复杂有机质。七十年来的综述也整理过古菌、细菌和真菌的多样性。但全球大规模微生物基因组调查里，蓝碳生态系统一直代表性不足。过去的方法精细程度有限，常常只能回答“谁比较多”，很难知道“谁在干什么”，更难判断环境压力会把整个群落推向哪里。这项研究想补的，就是大陆尺度上基因组级别的空白：微生物怎么分布，功能怎么分配，遇到生态和人为压力时能不能稳住。

研究者从2017年3月到10月，沿着中国海岸选了121个站点，其中116个有完整数据，包括34个红树林、57个盐沼和25个海草床站点。他们取的是表层0到20厘米沉积物。这一层新鲜有机质输入多，氧化还原变化剧烈，是微生物活动最活跃的地方。每个站点取5个子样混合，总共产生约1.5太字节测序数据，平均每个样品约13.3吉字节。测序就是读取环境DNA片段里的碱基序列，把泥里混杂的遗传信息变成可分析的文本。接下来，研究者把这些DNA碎片拼成一个个微生物基因组草图，学术上叫宏基因组组装基因组。拼出来的初稿有25326个，但很多要么残缺太多，要么是同一物种被不同站点反复拼出，像同一个人被不同相机拍了好几张几乎一样的证件照。研究者把低质量的剔除，再把相似度极高的合并成一个代表，这个过程叫去冗余。最后留下1954个不重复的中高质量草图，这些草图大致是1954个微生物基因组代表。你可以把每个草图想成一本从泥里拼出的微生物“说明书”，记录它可能会做什么。

研究者还整理了13292类功能基因。不同微生物里负责同一件事的基因，即使物种亲缘很远，只要功能相同，就被归到同一组，这叫功能基因直系同源组。13292个组就是13292类功能。

为了判断谁分布广、谁挑剔，研究者计算了每个草图的莱文生态位宽度。做法是先数每个基因组草图在多少站点、多少种环境里露面，露面越广、越均匀，说明越不挑环境。但光看露面次数还不够，因为常见菌天然就容易到处出现。于是研究者把站点标签打乱，让“谁在哪个站点”变成随机组合，重复999次，生成999套假数据。如果某个菌在真实数据里的分布宽度高于绝大多数假数据，说明它确实不挑环境，是泛化种；如果低于绝大多数假数据，说明它比随机情况更挑剔，是特化种。结果只有0.61%是泛化种，88.9%是特化种。这说明海岸环境过滤很强，大多数微生物只适应局部条件，换个地方就可能活不好。

从分类看，假单胞菌属最多，占38.7%，拟杆菌门占21.5%，放线菌门只有4.9%。40.8%的草图在三种生态系统里都能找到，盐沼独有种类最多，占15.9%。地理上，32.9%的草图在渤海、黄海、东海和南海四个海区都出现，29.8%只在单一海区出现。采样到98个站点时，已经覆盖95%的草图。这说明生态系统类型和南北空间共同塑造了微生物名单。

真正反直觉的是功能。研究者检查了19条代谢通路（代谢通路就是一串按顺序进行的生化反应，像工厂流水线，把原料一步步变成产品），其中17条在所有站点都能检测到。氮和硫的氧化还原最丰富，符合潮汐湿地经常缺氧的特点。碳固定方面，卡尔文循环最普遍（卡尔文循环是一种固定二氧化碳的流水线，把无机碳变成细胞能用的有机物，植物和很多微生物都用）。研究者只统计三种有可靠基因标记的碳固定途径，共找到1334个相关基因拷贝，其中卡尔文循环占998个。芳香化合物降解能力也很广。不同站点之间，微生物名单差异很大，典型两站点的名单几乎完全不同；但功能基因组合差异很小，典型两站点约79%相似。这种“不同物种干同样活”的现象，就是功能冗余。

不过，这里有两个统计口径要分清。布雷-柯蒂斯相异度同时看“有哪些物种”和“各占多少”。索伦森相异度只看“有或无”，不管多少。只改变数量，比如某个菌多了一点少了一点，功能可能不变。若一批菌彻底消失、另一批彻底取代，功能才可能大改。论文发现，按布雷-柯蒂斯算功能差异只有0.21，但按索伦森算会升到0.88。所以功能稳定主要怕“谁在”被换掉，不太怕“谁多谁少”。

碳固定是功能冗余的典型例子。10个碳固定潜力最高的草图分别属于脱硫杆菌门和绿弯菌门，亲缘关系很远。它们共享的功能基因只有25%，主要集中在能量保存和蛋白质合成等核心模块，附件功能则分散在不同类群中。负责碳固定的微生物种类在不同站点差别很大，但它们带的碳固定基因组合差别很小。碳固定功能差异小、分类差异大，二者明显脱钩。

环境压力方面，研究者为每个站点测了30种环境条件，包括pH、盐度、总氮、总磷、碳氮比、有机碳、温度、降水、经纬度、人口密度、农业用地比例等。26种与微生物组成显著相关，意思是这些条件不同，微生物名单也明显不同。地球化学条件解释分类周转的41.2%。所谓分类周转，就是沿着环境梯度，一个地点的微生物名单被另一个地点名单替换的程度，不是数量增减，而是“换人”。比如从北到南，温度、盐度逐渐变化，微生物物种也跟着逐渐替换。总磷、碳氮比、pH和盐度都是关键因素。农业用地在集水区增加约25%，分类差异就能从接近0升到0.2，说明轻微农业压力已让群落名单开始变。但功能基因组合被替换的程度一直很低。这意味着环境变、物种名单变，但“能做什么”变化小。总磷、碳氮比、农业用地和人口密度影响都不大，只有温度超过20摄氏度后出现轻微功能变化。这说明人为压力先改变物种名单，核心功能暂时还能稳住。

研究者还构建了三种生态系统的微生物互作网络，节点是基因组草图，边代表共现关系。他们随机移除节点，模拟物种丧失，观察网络效率怎么下降。盐沼在物种齐全时最不脆弱，但超过三分之二物种丧失后，网络最碎片化。到75%物种丧失阈值后，盐沼脆弱性增加30%。全球沿海湿地20世纪约损失63%，指的是湿地面积相比20世纪初减少约63%，这个模拟结果提醒我们，功能冗余不是无限保险。需要说明，这是计算机模拟预测，不是现场控制实验。

这项研究给蓝碳碳汇的稳定性增加了一份信心：即使物种名单变化，碳固定等核心功能仍可能维持。但它也表明，盐沼在严重生物多样性丧失下最危险。研究局限同样清楚：宏基因组草图只能代表可拼出的部分，稀有微生物大多还没被表征；测序深度有限，宏基因组测序像从大池子里随机捞DNA片段，捞得少，低丰度微生物容易漏掉，基因组拼不完整；空间比较和模型预测不能直接说成因果。未来需要更深测序和单细胞基因组等技术。保护红树林、盐沼和海草床，减少农业径流和人口压力，不仅是在保物种，也是在保它们背后的功能。

原文链接：https://www.nature.com/articles/s41467-026-77867-5?sessionid=