《Proceedings of the National Academy of Sciences》：Microbial niche partitioning explains N2O distributions and cycling in marine oxygen minimum zones
编辑推荐：
最小含氧带（Oxygen Minimum Zones, OMZs）是海洋中部的低氧区域，低氧重塑了氮循环并促进了强效温室气体及臭氧消耗物质——氧化亚氮（N?O）的排放。然而，塑造N?O分布的生态机制仍不确定，限制了对环境变化下排放响应的预测。为弥补这一空白，研
最小含氧带（Oxygen Minimum Zones, OMZs）是海洋中部的低氧区域，低氧重塑了氮循环并促进了强效温室气体及臭氧消耗物质——氧化亚氮（N?O）的排放。然而，塑造N?O分布的生态机制仍不确定，限制了对环境变化下排放响应的预测。为弥补这一空白，研究人员将基于性状的生态系统模型（该模型将微生物生长与化学能量可利用性联系起来）嵌入到涡分辨率的OMZ模拟中。该模型再现了观测到的N?O剖面，以及诸如亚硝酸盐积累等典型OMZ特征，并揭示了交叉摄食（cross-feeding）、竞争和资源限制如何创造微生物生态位，从而将N?O的生产与消耗区域分离开来。研究结果识别了在低氧水域中构建氮循环并调节温室气体通量的微生物相互作用。
最小含氧带（OMZs）是强效温室气体氧化亚氮（N?O）循环与排放的海洋热点区域。陡峭的氧化还原梯度孕育了多样化的微生物群落，这些群落既产生又消耗N?O，形成了特征性的垂直分布：缺氧核心区浓度较低，而周围氧跃层（oxyclines）中则积累至过饱和状态。尽管这种积累通常被归因于O?对N?O消耗的抑制作用强于对其产生的作用，但这种差异抑制并非始终被观测到，提示氮循环代谢之间的相互作用可能提供了额外的机制。为评估N?O剖面的生态控制因素，研究人员将基于性状的生态系统模型嵌入到东热带南太平洋（Eastern Tropical South Pacific, ETSP）OMZ的涡分辨率模拟中。在该框架下，氧化还原能量学、代谢物交换、竞争和摄食共同构建微生物生长，形成了微生物功能类型的涌现生态位，特别是好氧氨氧化菌和介导N?O循环的多种厌氧反硝化菌。在缺氧核心区，交叉摄食紧密耦合了产生和消耗N?O的反硝化菌，维持了低浓度。随着氧跃层中O?增加，与好氧代谢的竞争以及加剧的摄食逐渐排除了反硝化菌。然而，由于利用了更丰富的底物，产生N?O的反硝化菌比N?O消耗者能在更高的O?水平下持续存在，这反映了依赖中间代谢物种群的脆弱性。这种生态解耦驱动了N?O的积累，并在不施加生理性O?不耐受阈值的情况下捕捉了特征性的双峰垂直结构。更广泛地说，通过揭示生态位分化和生产者-消费者耦合如何塑造OMZ生物地球化学的其他特征，这些结果为预测未来海洋氮损失和N?O排放的变化提供了基础。
一、研究背景与科学问题
海洋最小含氧带（Oxygen Minimum Zones, OMZs）是海洋中部的低氧区域，仅占海洋体积不足0.1%，却因其邻近海岸上升流区而贡献了高达一半的海洋N?O自然排放。N?O既是强效温室气体，又是臭氧消耗物质。在OMZ中，N?O呈现出独特的垂直分布：缺氧核心区浓度低且常低于饱和，而上下氧跃层中则积累至过饱和，形成两个明显的峰值。长期以来，这种分布被归因于O?对N?O消耗步骤的抑制强于对其产生步骤的抑制这一生理学机制，但近期的实验和观测证据表明，这种严格的O?敏感性排序并非一致可观测到，提示微生物生态相互作用可能提供了额外的机制。
反硝化作用是OMZ中N?O的主要来源和汇，其过程为NO??→NO??→NO→N?O→N?。宏基因组调查显示，多数微生物仅编码反硝化途径的部分基因，表明完整的反硝化过程往往由不同种群分工完成。这些代谢模块可由不同种群执行，因此不一定占据相同的生态位。代谢中间产物如NO??和N?O成为连接上游生产者与下游消费者的核心纽带。理解N?O的垂直结构，需要解析N?O生产与消耗作为独立微生物功能在动态OMZ中如何空间组织，以及何种生态和环境条件允许N?O在核心区被消耗至低浓度，同时在氧跃层中积累。
二、研究内容与结论
研究人员将基于性状的生态系统模型嵌入到东热带南太平洋（Eastern Tropical South Pacific, ETSP）OMZ的涡分辨率三维模拟中。模型在原有浮游植物-浮游动物和元素循环框架基础上，增加了十种氧化还原活性的微生物功能类型，包括好氧氨氧化菌、好氧亚硝酸盐氧化菌、厌氧氨氧化菌（anammox）、好氧异养菌，以及六种构成模块化反硝化网络的反硝化菌类型。每种类型的生长受电子供体或受体的Monod动力学约束，生物量产量由底物氧化还原反应的热力学决定。模型不预设任何O?不耐受阈值，N?O的垂直结构完全从微生物生态相互作用中涌现。
研究结果表明，模型成功再现了观测到的N?O双峰垂直剖面、亚硝酸盐积累（次级亚硝酸盐最大值，SNM）等OMZ典型特征。在缺氧核心区，交叉摄食紧密耦合了产生和消耗N?O的反硝化菌，维持了N?O的低浓度。随着氧跃层中O?增加，好氧代谢的竞争和加剧的摄食逐渐排除反硝化菌，但由于产生N?O的反硝化菌利用更丰富的底物（NO??和NO??），它们比依赖中间代谢物N?O的消耗者能在更高的O?水平下持续存在。这种生态解耦驱动了N?O的积累，并自然产生了双峰垂直结构，无需假设生理性O?不耐受阈值。
三、关键技术与方法
研究采用区域海洋模型系统（Regional Ocean Modeling System, ROMS）耦合生物地球化学元素循环模型（Biogeochemical Elemental Cycling model, BEC），在ETSP OMZ的涡分辨率气候态配置中运行。模型嵌入基于性状的生态系统框架，包含十种微生物功能类型，其生长由氧化还原能量学、Monod动力学和Liebig最小定律约束。生物量产量由各代谢途径的Gibbs自由能决定。模型经过25年旋转达到准平衡。观测数据来自PANGAEA数据库汇编的ETSP区域数据集，按位势密度分箱以对齐水团界面。
四、研究结果
微生物功能类型。 模型包含化能自养菌（好氧氨氧化菌、好氧亚硝酸盐氧化菌、anammox）和异养菌（好氧异养菌和六种反硝化菌）。反硝化菌包括三种单步还原类型（NO??→NO??、NO??→N?O、N?O→N?）和三种多步还原类型（NO??→N?O、NO??→N?、NO??→N?），多步类型内部还原中间产物而不与环境交换。
资源竞争层级。 单步反硝化菌中，产量随还原步骤增加而提高，使N?O还原菌成为有机质的最强竞争者。多步反硝化菌的DIN产量随途径长度增加，但以降低有机质效率为代价。通过资源生存浓度（R*）框架分析，好氧异养菌在O?存在时可竞争排除所有反硝化菌。N?O还原菌虽理论上具有最高产量优势，但始终受N?O供应限制，其种群动态与局部N?O供给耦合。
ETSP涌现生物地球化学结构。 模型再现了观测的O?、N?O、NO??和NO??垂直分布。N?O在缺氧核心区低浓度，在上下氧跃层形成双峰。NO??在SNM处积累。这些结构从微生物生态位分化中涌现，而非预设生理阈值。
五、讨论与结论
研究揭示了生态位分化和生产者-消费者耦合如何塑造OMZ氮循环和N?O分布。N?O产生者因利用更丰富的底物而能在较高O?下持续，N?O消耗者则因依赖中间代谢物而更脆弱，这种生态解耦驱动了N?O积累。研究结果提供了预测未来环境变化下海洋氮损失和N?O排放变化的机制基础，表明微生物相互作用而非单纯的生理O?敏感性是控制OMZ中N?O分布的关键因素。
