极端氨胁迫下厌氧微生物组的适应性进化揭示代谢可塑性与功能冗余

《Environmental Science & Technology》：Adaptive Evolution Reveals Metabolic Plasticity and Functional Redundancy in an Anaerobic Microbiome under Extreme Ammonia Stress

　　

氨毒性是控制城市固体废弃物有机组分厌氧消化过程中微生物群落结构和性能的主要生化瓶颈。然而，微生物适应长期高氨水平的机制基础仍未得到充分表征。本研究采用在逐步升高的氨浓度（350–1500 mgN L^?1）下进行的长期连续富集策略，并结合以基因组为中心的宏基因组学（genome-centric metagenomics）和宏转录组学（metatranscriptomics），解析了有机废弃物降解微生物组在240天内的响应。升高的氨压力诱导产甲烷作用逐步下降，并伴随挥发性脂肪酸（volatile fatty acids, VFA）的积累，尤其是乙酸盐。尽管存在这些抑制压力，甲烷产量相对于初始基线仅减半，反映出产甲烷群落的韧性。这种稳定性由微生物组的结构重组驱动，其中不同分类群间的功能冗余（functional redundancy）维持了核心代谢功能。关键的适应性响应包括碳固定途径的重构，具体表现为Wood–Ljungdahl途径（Wood–Ljungdahl pathway, WL）变体与甘氨酸裂解系统（glycine cleavage system, GCS）偶联，作为替代的乙酸氧化途径，以及持续的渗透保护剂生物合成。细胞稳态通过多种能量转换复合物介导的H⁺补充和K⁺内流以维持阳离子–质子平衡得以保持。综上，这些发现表明代谢可塑性和核心代谢功能的维持是氨韧性的主要决定因素，在抑制条件下维持甲烷产生。

**论文解读：极端氨胁迫下厌氧微生物组的代谢可塑性与功能冗余**

**一、研究背景与问题提出**

全球城市固体废弃物年产生量约20亿吨，其中有机组分（organic fraction of municipal solid waste, OFMSW）占42%–69%，其高效处理是环境可持续发展的关键挑战。厌氧消化（anaerobic digestion, AD）作为成熟的生物转化技术，可将有机废弃物转化为沼气等可再生能源。AD过程包括水解、产酸、产乙酸和产甲烷四个阶段，其中产甲烷阶段对操作胁迫最为敏感，尤其是高氨氮负荷。在氮含量丰富的底物处理中，蛋白质降解释放活性氮，总氨氮（total ammonia nitrogen, TAN）在水中以铵离子（NH₄⁺）和游离氨（free ammonia nitrogen, FAN）的形式存在平衡，FAN可被动扩散透过细胞膜，破坏胞内pH并抑制酶活性。文献报道的氨抑制阈值差异极大（27–1450 mg FAN L^?1），这主要源于FAN对pH和温度的强依赖性。在核心AD微生物群中，产甲烷古菌最为敏感，尤其是乙酸营养型产甲烷菌。在高氨条件下，乙酸营养途径常被抑制，转而有利于互营乙酸氧化（syntrophic acetate oxidation, SAO）——乙酸被转化为CO₂和H₂，随后供给氢营养型产甲烷菌。尽管对氨抑制的过程层面后果已有较多认识，但微生物群落如何在长期氨胁迫下维持功能性的生态学和代谢机制仍不完全清楚。既往研究已鉴定出耐氨分类群、渗透适应策略和向SAO的转变，但群落更替、代谢重组与转录活性之间的联系仍缺乏解析。为此，本研究旨在通过长期富集策略结合多组学方法，阐明微生物群落组成变化如何转化为功能适应，以及哪些代谢途径在适应氨的群落中被积极维持。

**二、研究方法概述**

研究人员利用来自意大利北部废弃物处理厂的OFMSW作为碳源，采用经BA培养基活化的工业微生物接种物，在120 mL血清瓶中进行三组重复的富集实验。实验以24天为一个世代，在10个世代（G1–G10）共240天内，通过逐步添加NH₄Cl实现FAN浓度从350 mgN L^?1到1500 mgN L^?1的四倍以上递增，并通过定期调节pH至8.0±0.1维持目标FAN水平。研究整合了以基因组为中心的宏基因组学（共组装获得71个宏基因组组装基因组（metagenome-assembled genomes, MAGs），其中51个为高质量）和宏转录组学分析，采用SPAdes共组装、多工具分箱、GTDB-Tk分类注释、eggNOG-mapper功能注释、KEGG模块完整性分析、dbCAN碳水化合物活性酶分析、PhyloPhlAn系统发育推断以及Strobealign转录本比对和TPM定量等方法，在群落和功能两个层面系统解析了微生物组的适应性进化轨迹。

**三、研究结果**

**3.1 渐进性氨负荷对甲烷产量和挥发性脂肪酸波动的影响**

累积甲烷产量在G1–G3保持稳定（约275 mL/24天），但在G4（780 mgN L^?1 FAN）出现超过50%的减产。尽管抑制环境持续加剧，甲烷产量在FAN高达1393 mgN L^?1（G9）时仍维持在约100 mL，表现出显著的微生物抗性。至G10系统达到终末状态，甲烷产量急剧下降并伴随延长的恢复滞后期。CO₂产生轨迹与甲烷相似，而H₂在G6前被完全消耗无积累，从G7起H₂浓度显著增加，表明发酵微生物对高FAN的耐受性优于产甲烷菌。VFA谱显示，从G6起乙酸盐在每世代末期持续积累，反映产甲烷菌群处理中间代谢产物的能力下降，与乙酸营养型产甲烷向SAO-氢营养途径的转变一致。

**3.2 分类学演替驱动乙酸营养型产甲烷菌被耐氨氢营养型替代**

宏基因组重建获得71个MAGs（62个细菌、9个古菌），以梭菌纲（Clostridia）为主（约59%）。G1时群落呈异质组成，具有九种产甲烷古菌，氢营养型占约21%、乙酸营养型约11%、甲基营养型约2.6%。从G2起，Methanosarcina sp. EV05成为优势产甲烷菌（G3时占36%），但随FAN升高，Methanoculleus bourgensis EV08的丰度从G3起逐步上升，伴随Methanosarcina sp. EV05的逐渐衰退。从G6（FAN≥1026 mgN L^?1）起，M. bourgensis EV08成为优势产甲烷菌（RA>15%），而Methanosarcina sp. EV05在G9末完全消失。细菌群落方面，早期优势菌Proteiniphilum massiliensis EV33、Syntrophomonadaceae sp. EV26和Clostridium sartagoforme EV58（RA分别为22.4%、6.5%和4.1%）被耐氨特化菌群替代，至G10时Tepidimicrobiaceae sp. EV16和Xylanivirgaceae sp. EV69占优势（RA>25%）。G10的基因表达谱证实优势分类群的RA与转录活性高度对应，M. bourgensis EV08的表达水平占总古菌活性的93.5%。

**3.3 韧性互营菌的基因组蓝图鉴定**

KEGG模块分析显示，Wood–Ljungdahl（WL）途径在系统中可能双向运行——早期支持乙酸产生，随选择压力增加则转向SAO方向。多个优势分类群编码WL途径的非典型变体：Tepidimicrobiaceae谱系缺乏典型的羰基分支，甘氨酸裂解系统（GCS）被预测作为代谢旁路，将乙酸经由甘氨酸通道化，甘氨酸裂解产生CO₂和还原性C1单元，可进一步生成H₂。Tepidimicrobiaceae sp. EV16、Schnuerera sp. EV29和Anaerosalibacter sp. EV71缺乏典型的乙酸活化酶（pta和acs），但编码丁酰磷酸转移酶（ptb）作为替代策略。丙酸降解途径在所有MAGs中均不完整，唯一具有强丙酸产生潜力的P. massiliensis EV33在G5后的消失与丙酸浓度显著下降在时间上吻合。丁酸代谢途径则广泛分布，在Tepidimicrobiaceae、Schnuerera、Clostridium和Anaerosalibacter等属中完整存在。渗透保护方面，多数耐氨成员依赖谷氨酸积累（hutH、hutU、hutG、putA、putB基因），Virgibacillus promii EV40是唯一编码甘氨酸甜菜碱合成机制的物种（gbsA、gbsB），而M. bourgensis EV08编码β-赖氨酸N6-乙酰转移酶（ablB）并携带ABC型季胺转运系统（ProVWX）和Trk/Ktr钾摄取系统。

**3.4 离子转运体和渗透保护系统在G10的转录激活**

G10转录组分析表明，Tepidanaerobacteraceae谱系主要表达经典WL基因，而多个推定SAOB（Tepidimicrobiaceae成员、Schnuerera sp. EV29、S. ultunensis EV37、Anaerosalibacter sp. EV71）强烈表达GCS相关基因（grdE、grdB、gcvPA、gcvPB、gcvT、gcvH），且ptb基因中等至高表达。这些菌株同时展现完整且活跃转录的丁酸降解途径，提示高度的代谢多功能性。胁迫响应方面，ATP依赖性钾转运系统（Trk/Ktr）和Na⁺/H⁺反向转运体（Nha）高表达，与关键呼吸复合物（rnfA、rfnB、rfnC、nouG）的上调相耦合。Rnf型黄素蛋白复合物在多个群落成员中高表达，将铁氧还蛋白还原与反向电子传递偶联以维持质子动势。渗透保护方面，ProVWX甘氨酸甜菜碱摄取系统在Tepidimicrobiaceae sp. EV16、Anaerosalibacter sp. EV71和M. bourgensis EV08中高表达，提示细胞外甘氨酸甜菜碱被关键群落成员主动同化。M. bourgensis EV08还表达albB基因（N6-乙酰-β-赖氨酸产生）和arcC基因（氨甲酰激酶，催化NH₃和碳酸氢盐可逆转化为氨甲酰磷酸），参与胞内氮解毒。

**3.5 功能冗余维持AD核心代谢功能**

群落核心功能由四类代谢行会维持：有机物降解菌、SAO、同型产乙酸菌和产甲烷菌（分为乙酸营养型和氢营养型）。早期阶段，Methanosarcina sp. EV05和Methanoculleus sp. EV03介导乙酸营养和氢营养产甲烷，P. massiliensis EV33作为关键初级降解菌兼具水解和乙酸代谢双重功能。随FAN升高，早期定殖者被功能等效但更耐受的分类群替代：Xylanivirgaceae sp. EV69取代P. massiliensis EV33成为主要复杂碳水化合物降解者，Tepidimicrobiaceae谱系（EV16和EV50）处理蛋白质和纤维素组分，Anaerosalibacter sp. EV71和Tepidimicrobiaceae sp. EV15贡献互补的淀粉分解和果胶分解活性。从G6起，M. bourgensis EV08成为优势古菌，其多种保护机制确保高效稳态响应。Tepidimicrobiaceae谱系成员和Anaerosalibacter sp. EV71将丁酸降解与WL/GCS旁路偶联，在乙酸营养途径被抑制时建立VFA产生与乙酸消耗之间的代谢桥梁。群落水平上，早期阶段具有高功能冗余，随氨水平升高冗余度逐步精简，选择少数高度特化的耐胁迫微生物，但在G10功能冗余度再次增加，提示适应群落内的二次多样化。

**四、讨论与结论**

本研究表明，极端FAN压力下AD过程的持续不取决于分类学稳定性，而主要由功能冗余决定。随氨浓度升高，胁迫敏感种群被代谢等效但更耐受的微生物替代，驱动代谢网络重组，最终乙酸营养途径被SAO偶联氢营养产甲烷取代，使甲烷产生在严重抑制条件下得以延续。这种向氢营养产甲烷的转变以广泛的代谢重布线为基础，包括Tepidimicrobiaceae谱系利用非经典WL/GCS模块和ptb介导的乙酸活化。尽管这种基因组可塑性维持了SAO并防止系统完全崩溃，但离子稳态和渗透保护所需的高转录投入施加了严格的能量负担，最终导致甲烷产量下降。研究结果揭示了功能冗余和代谢重布线而非分类学稳定性是AD长期耐氨的主要驱动力，为处理富氮底物的全规模消化器运行提供了实践指导：逐步递增的FAN负荷可显著扩展耐受范围，乙酸盐积累和H₂分压上升可作为即将抑制的早期预警信号，氢营养和乙酸氧化分类群的丰度可作为成功适应群落的微生物生物标志物。此外，补充甘氨酸甜菜碱等相容性溶质可能缓解渗透和离子稳态的能量负担。研究也存在一定局限性：使用受控批次系统和单一接种物简化了全规模反应器的固有变异性；NH₄Cl的使用导致离子强度随富集过程增加，观察到的适应轨迹可能反映高FAN与渗透胁迫的联合响应；宏转录组分析仅限于终末阶段；途径重建依赖基因组和转录代理而非直接通量测量。未来的研究应通过¹³C代谢通量分析、靶向蛋白质组学和候选菌株分离验证非经典WL/GCS乙酸氧化途径，并采用连续流反应器测试适应性轨迹在实际运行条件下的适用性。