南京大学任鑫坤/王瑾丰团队ES&T:基于生态学原理设计合成菌群,实现抗生高效降解与稳定定植

推荐理由

图文摘要

文章简介

图文导读

研究背景：

本研究开篇指出现有废水处理工艺的痛点：传统活性污泥法主要针对常规有机物和营养物设计，对磺胺甲恶唑等新兴污染物缺乏专门的分解代谢酶系和选择压力，导致大量此类化合物穿透处理工艺，在环境中持续检出，对生态和健康构成风险。研究进一步指出，当前合成微生物群落设计多呈"功能驱动型"，侧重体外降解能力而严重忽略生态兼容性，导致工程菌群投入复杂生物系统后因扩散、环境过滤及与土著群落竞争而难以建群或维持活性。因此，作者提出需将分解代谢功能与生态适应性深度整合——通过宏基因组识别候选降解菌，评估其群体感应能力（增强环境过滤韧性）和代谢互补网络中的关键物种地位（维持群落稳定性），最终理性设计出兼具高降解效能与生态韧性的合成微生物群落。这一框架为突破生物强化定植瓶颈提供了全新思路。

图1 本研究的总体概念框架与工作流程

1. SMX降解微生物的鉴定及生态学特性

本研究首先构建了针对SMX初始降解酶（SaDA、SaDB、SaDC）的同源注释管道（图2B）。通过系统优化BLASTp比对参数，确定在氨基酸一致性阈值为50%、期望值阈值为1×10⁻¹⁵时，F1分数和马修斯相关系数达到最优平衡（0.78和0.72）（图3C、3D），有效提升了宏基因组中功能基因检出的敏感性与特异性。

利用该管道，从全球活性污泥高质量宏基因组组装基因组中鉴定出403个编码完整SMX初始降解酶系的“初级降解菌”，主要分布于放线菌门（59.5%）和假单胞菌属（36.6%）（图S1、S2）。系统发育分析表明，sadABC基因簇在这些菌门中主要通过垂直遗传传递（图4A），基因水平转移罕见，提示该分解能力在传代中具有高度基因组稳定性。同时，筛选出17个能处理下游中间产物1,2,4-三羟基苯的“次级降解菌”。

进一步评估生态兼容性发现，具有群体感应功能的降解菌在群落中的结构关键性显著高于非群体感应菌（p < 0.01），所有群体感应菌均排在群落成员影响力的前50%（图4B）；且呈现更窄的生态位宽度（p < 0.05）（图4C），符合“窄谱资源利用型专性菌”特征，可减少资源竞争，增强群落功能韧性。

图4 SMX降解菌中 sadABC基因簇的系统发育与群体感应功能特征

2. 关键物种引导的候选菌筛选与代谢交叉喂养网络

基于基因组尺度代谢模型重建了全球活性污泥微生物组的代谢互补网络（包含360个初级和12个次级SMX降解菌，解析为17个非随机模块）（图S4），并识别出网络关键物种（图5A）。将“SMX降解能力”、“群体感应功能”和“网络关键物种地位”三重筛选标准相交，最终从89个MAGs中降重获得5个非冗余菌种（图5B），分别注释为Nocardiopsis corallilcola、Gordonia hydrophobica、Haloglobus pacificus、Mycobacteriumsp020438785和Dokdonella fugitiva。

代谢物交换预测显示，这5个物种间存在高度耦合的交叉喂养网络，共涉及37种可交换代谢物（图5C）。氨基酸及其类似物占主导（半胱氨酸、甘氨酸、谷氨酸和赖氨酸），构成了互惠共生的主要公共物品（表S3）；次级为无机离子（Ca²⁺、Mg²⁺、亚硝酸盐、磷酸盐）（表S4）和碳水化合物（D-葡萄糖、L-阿拉伯糖等）（表S5）。氨基酸介导的代谢互养有效减少了竞争排斥，建立了种间依赖性，结合垂直稳定的分解基因簇、群体感应及关键物种属性，确立了这5个物种为构建SMX降解SynCom的最佳基石。

图5 合成微生物群落（SynCom）候选菌株的筛选

3. 代谢模型引导的SMX降解SynComs最小化与优化

为在性能与简约性间取得平衡，研究利用GSMMs模拟了5个候选菌所有可能的2、3、4、5成员组合（共15种）的SMX摄取通量（图6A）。结果表明，二成员组合平均通量为21.70 ± 2.73 mmol gDW⁻¹ h⁻¹；加入第三个初级降解菌后，平均通量显著提升25%至27.12 ± 3.57 mmol gDW⁻¹ h⁻¹（p < 0.01）；进一步扩增至四或五成员则无统计显著增益，呈现功能饱和现象。

在所有测试组合中，SynCom11（G. hydrophobica、H. pacificus、N. corallilcola、D. fugitiva）预测通量最高（31.71 mmol gDW⁻¹ h⁻¹），其次为SynCom5（G. hydrophobica、H. pacificus、D. fugitiva，通量30.70 mmol gDW⁻¹ h⁻¹）（图6A）。机制在于，次级降解菌D. fugitiva需G. hydrophobica供给半胱氨酸；而G. hydrophobica与H. pacificus构成紧密互养对（双向交换半胱氨酸、酪氨酸、甘氨酸和缬氨酸）；加入N. corallilcola提供了两者所需的异亮氨酸，关闭关键代谢缺口（图6B）。值得关注的是，代谢互作潜力与资源重叠的比值（MIP/MRO）与SMX摄取通量呈显著正线性相关（R² = 0.42，p < 0.05）（图6D），表明高降解性能源于合作与竞争间的有利平衡，而非单纯消除竞争。

图6 代谢模型引导的高性能SMX降解SynComs的设计与互作分析

4. 模型预测SynComs的实验验证

纯培养验证中，G. hydrophobica、H. pacificus和N. corallilcola在72小时内分别降解90.9%、86.5%和88.1%的SMX并积累THB；而D. fugitiva虽不降解SMX，但可去除94.6%的THB，确认了分工代谢路径（图7A）。SynCom5和SynCom11在纯培养中均实现近完全SMX去除，且SynCom11效率略高，与模型预测高度一致（图7B）。

在活性污泥微宇宙中接种SynCom11，7天后SMX去除率达91.3%，显著高于未接种对照的25.8%（图7C）。16S rRNA测序证实，Gordonia（从约0.001%升至>0.2%）和Dokdonella（从<0.6%升至1.1%）成功定植，总定植率达50%（图7D）。α多样性分析显示，SynCom11接种组Shannon指数高于非接种组（图7E），表明快速去除SMX缓解了抗生素对土著群落的胁迫，保存了微生物多样性。功能预测显示接种增强了细菌运动性、群体感应、信号转导和中心碳代谢通路（图7F）。研究同时指出，其余两株菌未检出可能源于竞争排斥、原生动物捕食或生理条件不匹配等因素。

图7  两个最优SynComs（SC5和SC11）的实验验证

主要结论

本研究通过构建"功能-生态"双驱整合框架，成功设计了针对磺胺甲恶唑的高效合成微生物群落。宏基因组引导的降解潜力鉴定结合群体感应与网络关键物种筛选，确保了菌株兼具分解能力与生态适应性。GSMMs模拟与实验验证一致表明，SynCom5（三成员）和SynCom11（四成员）表现出最高降解效能，其性能源于互补的氨基酸交叉喂养所建立的代谢互赖网络，以及合作与竞争间的有利平衡（MIP/MRO比值）。在活性污泥微宇宙中，SynCom11实现91.3%的SMX去除和50%的成员定植率，并通过缓解抗生素胁迫保护了土著微生物多样性。该策略规避了高浓度选择压力带来的抗性风险，为新兴污染物的生物强化提供了从生态理论到应用的普适性解决方案。

往期推荐