Pseudomonas 属呼吸链组成的进化多样化与 Cyo bo3 型末端氧化酶调控分化

《Microbiological Research》：Evolutionary diversification of respiratory chain composition and the Cyo bo?-type terminal oxidase regulation across the genus Pseudomonas

　　

研究人员对 Pseudomonas 属呼吸链组成和 Cyo bo₃ 型末端氧化酶调控的进化多样化进行了研究。Pseudomonas 属好氧电子传递链（ETC）已发生广泛多样化；Cyo 末端氧化酶（terminal oxidase）调控在 Pseudomonas 物种间存在差异；Cyo 氧化酶缺失导致广泛的转录重编程（transcriptional reprogramming）；Cyo 末端氧化酶启动子（promoter）演化为五个调控类别；启动子分化驱动物种特异性 Cyo 表达模式。

Pseudomonas 属呼吸链组成和 Cyo bo₃ 型末端氧化酶调控的进化多样化解读

研究背景与问题。Pseudomonas 是 Gamma-proteobacteria 中系统发育多样、生态分布广泛的类群，包括植物或动物机会致病菌如 Pseudomonas syringae、Pseudomonas aeruginosa，以及植物或土壤有益菌如 Pseudomonas fluorescens、Pseudomonas putida、Pseudomonas rhizosphaerae。其生态成功很大程度上归因于代谢多样性，而需氧电子传递链（ETC）的灵活性是其中关键。许多细菌具有分支 ETC，含多个末端氧化酶，这些酶在氧亲和力、氧化还原特性和质子泵效率上不同，细胞通过调节其相对贡献优化能量保存和氧化还原平衡。P. aeruginosa 和 P. putida 编码五个好氧末端氧化酶：两个泛醌氧化酶 Cyo 和 CIO，以及三个细胞色素 c 氧化酶 Aa₃、Cbb₃-1 和 Cbb₃-2。P. aeruginosa 中 Cyo、CIO、Aa₃ 低氧亲和力，Cbb₃-1 和 Cbb₃-2 高氧亲和力；CIO 在其他氧化酶被氰化物抑制时关键。调控方面，P. aeruginosa 多个系统协调氧化酶表达，包括氧响应调节因子 ANR、氧化还原敏感双组分系统 RoxS/RoxR、稳定期 sigma 因子 RpoS、氰化物响应蛋白 MpaR 和 ZapAL；P. putida 中 ANR 抑制 Cyo 和 CIO 并诱导 Cbb₃-1，但不影响 Aa₃ 和 Cbb₃-2。与 P. putida 不同，P. aeruginosa Cyo 不响应 ANR。Cyo 在富培养基高氧条件重要，氧张力降低或进入静止期时 cyoABCDE 表达下降，并响应 H₂O₂、酸 pH 和铁限制。除氧可利用性外，呼吸链组成与碳代谢存在功能联系，但机制不明。Cyo 丢失增加 CIO 和 Cbb₃-1 表达并损害多种碳氮源生长；Cyo 在 Pseudomonas 种间是否具类似角色尚不确定。cyoABCDE 启动子区在 P. putida 和 P. aeruginosa 间显著不同，尽管酶高度保守，提示调控分化与生态位适应相关。因此，研究人员开展联合系统发育基因组学、启动子水平和功能分析，覆盖 249 个 Pseudomonas 物种，以定义 ETC 组分分布、表征 cyoABCDE 调控区多样化、实验验证启动子分化是否导致不同调控行为，并分析 Cyo 丢失对 P. putida KT2440 转录组的影响。

关键技术方法概括。研究人员使用 248 个 Pseudomonas 模式菌株加 P. putida KT2440，共 249 个基因组；基因组来自 NCBI，经 CheckM 过滤 ≥90% 完整度且 ≤5% 污染，Bakta 重注释。以 493 个核心基因簇推断核心基因组系统发育。末端氧化酶位点通过 curated multi-reference homology-search 识别，使用 DIAMOND 数据库，阈值 ≥30% 氨基酸同一性和 ≥70% query coverage，按邻近关系聚类并分类完整、部分或缺失。提取 208 个含完整 cyoABCDE 菌株的 cyoA 启动子 ?500 至 +50 区，MAFFT 比对、IQ-TREE 建树，定义五组启动子。构建 P_cyoA-gfp 报告质粒，在五株代表菌中测绿色荧光蛋白（GFP）。对 P. putida KT2440 与 KT2440cyoB 突变体做 RNA 测序（RNA-seq）和实时定量 RT-qPCR，并进行竞争实验。相关作者为 Bárbara Pérez-Arroyo、Zulema Udaondo、Rubén Delicado-Mateo、Luis Yuste、Renata Moreno 和 Fernando Rojo。

研究结果。3.1 电子传递链组成在 Pseudomonas 物种间变化。通过 249 个基因组分析，cyoABCDE 存在于 208 株（83.5%），41 株缺失（16.5%）。CIO 为 98.4%，Aa₃ 为 87.1%，Cbb₃-1 为 91.6%，Cbb₃-2 为 60.2%。部分菌株有 cbb₃-1 和 cbb₃-2 不完整基因簇，31 株含 divergent CIO-II，23 株 CIO-II 不完整。多数缺 Cyo 菌株含 Aa₃ 和 CIO 或 CIO-II；所有缺 Aa₃ 菌株有 Cyo 和 CIO，除 P. fluvialis 以 CIO-II 替代。175/249（70.3%）与 P. putida 或 P. aeruginosa 共享同一末端氧化酶组，9 株另含 CIO-II。Mantel 检验（Mantel test）显示氧化酶谱受系统发育显著结构化，ordinal r=0.534，P=0.0001；binary r=0.502，P=0.0001，各家族均有显著系统发育信号，Cyo 最强，Cbb₃-2 最弱。生态元数据不足关联特定氧化酶组与栖息地。

3.2 在 P. putida KT2440 中，Cyo 表达受碳源强烈调节。通过 P_cyoA-gfp 报告系统，在富 LB 中，P_cyoA 在指数早期强诱导，mid-exponential 后急剧下降；高通气下结果不变。最小盐培养基含 casamino acids、acetate 或 malate 也短暂诱导，但荧光低 4.4 至 5.7 倍；fumarate、succinate、protocatechuate、p-coumarate 更低，glucose、citrate、pyruvate、octanoate 进一步低。最大 P_cyoA 活性与倍增时间无显著线性相关（Pearson r=?0.198，P=0.517，n=13）。benzoate 在早期指数期强诱导，接近静止期降半后恒定，独立于生长速率，提示条件特异性调控。

3.3 cyoABCDE 启动子区序列在 Pseudomonas 物种间不同。比较 208 个完整 cyoABCDE 菌株的 cyoA AUG 前 ?500 至 +50 区，包含 P_cyoA 和长 5'-UTR。系统发育和序列比较定义五组启动子 A-E。组内平均同一性：A 为 84.6%，B 为 81.1%，C 为 69.1%，D 为 74.4%，E 为 42.8%；组间 B-C 为 55.1%，B-D 为 58.2%，C-D 为 52.2%；E 内部异质性高。A 组含 P. putida KT2440 等 39 种，通常土壤腐生或根际；B 组含 P. fluorescens/P. protegens 簇，植物相关；C、D 相似更低，D 含 P. syringae 相关植物病原；E 组 25 株含 P. aeruginosa PAO1，与 A 组无可检测相似性。构建各组代表菌 P_cyoA-gfp 在 LB、succinate、citrate、glucose、malate 中比较。与 P. putida 不同，B-E 四株在 LB 中未重现强早期指数诱导，仅 P. rhizosphaerae 有弱爆发，P. fluorescens 或 P. syringae 几乎无活性；P. aeruginosa 指数期低，静止期缓慢增加。P. rhizosphaerae 在 malate 或 succinate 中活性较高；有些菌静止期增加。P_cyoA 活性随生长阶段、碳源和菌株显著变化，无明显生长速率直接关系，支持 cyo 基因调控在谱系间广泛分化。

3.4 Cyo 缺失触发全局转录重编程。通过比较 P. putida KT2440 与 cyoB 灭活突变体 KT2440cyoB 高通气 LB 早期指数期转录组，发现 271 基因上调、271 基因下调 ≥2 倍（Padj<0.05）。ETC 中 CIO 泛醌氧化酶基因强增，cioA 为 13 倍，cioB 为 19 倍，Padj<1.3E-231，RT-qPCR 确认；bc1 复合物上调，petA 为 2.8、petB 为 2.4、petC 为 2.2 倍；Cbb₃-1 基因中度上调，ccoN1 为 1.8、ccoO1 为 1.7、ccoQ1 为 1.8、ccoP1 为 2.0 倍，Padj<7.9E-13，低于 2 倍阈值；Aa₃ 和 Cbb₃-2 不变。电子载体 azurin、cytochrome c4、cytochrome c5 mRNA 增加 3、3、2.2 倍，nuo 基因下降 2 至 2.8 倍。碳摄取和代谢改变：葡萄糖摄取代谢基因下调，三羧酸（TCA）循环关键酶下降，氨基酸和铁摄取改变，NtrB-NtrC 和谷氨酰胺合成酶下调 >2 倍。小 RNA RsmZ 增加，CrcZ 和 CrcY 减少；应激蛋白、热激蛋白和烷基过氧化氢还原酶上调。高通气摇瓶 cyoB 突变倍增时间增加至 37%；竞争实验中 20 代后 cyo/wt 比例从 1:1 降至约 0.2:1。

讨论与结论翻译。讨论指出，Pseudomonas 属好氧呼吸链显著多样化，单株有不同末端氧化酶组合。多数物种编码 Cyo、CIO、Aa₃ 和至少一种 Cbb₃ 型氧化酶，但有些仅保留两种，如 Cyo/CIO 的 P. baltica、P. luteola 或 CIO/Aa₃ 的 P. oryzagri。变异提示呼吸在适应多样栖息地中有重要进化作用；氧化酶谱与核心基因组系统发育显著关联，说明呼吸链组成与谱系历史强相关。生态元数据不足建立特定组合与栖息地稳健关联。CIO 几乎无处不在，可能居中心生理角色；CIO 无铜，耐氰化物和硫化氢，质子转运效率较低，这种耐受性可能是关键适应。Cyo 在多数基因组存在但多种条件下非必需，表明 ETC 功能冗余。主要发现是 Cyo 调控广泛分化：启动子序列变异大，Cyo 表达模式即使氧化酶谱相似也差异大。P. putida、P. aeruginosa 和 P. fluorescens 编码相同好氧末端氧化酶但 Cyo 表达显著不同。启动子组与生态特化相关：A 组多为土壤腐生或根际细菌；B/C 组含非致病植物相关 P. fluorescens 和 P. protegens 簇；D 组植物病原如 P. syringae；E 组含 P. aeruginosa 和 P. luteola 等动物机会致病菌，且 cyoA 启动子序列无关。A 组高序列同一性提示大多数含 ANR 结合位点，如 P. putida KT2440；B、D、E 组缺乏可识别 ANR box。B-E 组控制转录的调控元件不清楚。Cyo 调控变异表明呼吸控制可独立于结构组分进化，可能与生态特化相关。数据还揭示呼吸调控与碳代谢紧密耦合：P. putida KT2440 中 Cyo 表达随碳源和生长阶段强变化，与生长速率无相关，提示细胞内代谢状态如 ATP、NADH 或氧化还原平衡可能决定 ETC 配置；机制未知。RoxSR 和 HskA 等响应氧化还原信号，可能影响某些末端氧化酶表达。苯甲酸强诱导 Cyo，提示芳香化合物胁迫下 Cyo 维持氧化还原稳态和支持能量需求应激反应。cyoB 失活导致广泛转录组重组，细胞通过强诱导 CIO 及较小程度 bc1、Cbb₃-1 和电子载体补偿；氨基酸摄取增加，葡萄糖摄取代谢、TCA 酶和 NADH dehydrogenase I 下降，NtrB-NtrC 和谷氨酰胺合成酶下调提示整体代谢减速。cyoB 失活降低富培养基高氧生长速度，突出 Cyo 重要性。研究结论翻译如下：总之，这些发现表明 Pseudomonas 物种共享一个可变的呼吸骨架，通过谱系特异性调控程序差异部署，可能允许独立调节呼吸性能以响应谱系特异性生态压力，包括氧和营养可利用性变化、氧化还原胁迫和宿主相关生活方式。在此框架中，Cyo 可能代表一个灵活模块，其呼吸贡献取决于环境背景和物种特异性调控架构。cyoABCDE 调控与中心代谢紧密整合将使菌株平衡能量效率、氧化还原稳态和胁迫耐受性，从而促进生态多样化。然而，研究人员注意到 cyoABCDE 表达分析仅在五个代表性菌株中进行，cyoB 突变体转录组分析仅限于 P. putida；因此，将这些观察外推到 Pseudomonas 属时应谨慎。实际意义是，调节呼吸链活性可优化工程菌能量代谢和氧化还原平衡，理解碳源和调控回路如何塑造电子传递可指导增强 Pseudomonas 生物技术平台代谢效率和产量。