重水（D2O）胁迫下大肠杆菌代谢组学特征解析：菌株特异性耐受的生化基础

《Metabolomics》：Metabolic profiling of Escherichia coli under deuterium oxide-induced stress

　　

摘要 引言：部分细菌能够耐受重水（D2O）并高效生长，而另一些细菌则表现出较慢的生长表型。因此，对某些细菌而言，D2O通过动力学同位素效应（kinetic isotope effects，KIEs）干扰依赖氢的酶促过程，但由此产生的细菌代谢适应性尚未被完全理解。尤其是，菌株特异性D2O耐受性背后的生化因素仍不清楚。 目标：本研究旨在表征D2O敏感型与D2O耐受型大肠杆菌（Escherichia coli）分离株的全局代谢响应，并鉴定与同位素胁迫耐受性相关的代谢通路和代谢物。 方法：研究人员将5株尿路致病性大肠杆菌（uropathogenic E. coli，UPEC）分离株及参考菌株E. coli K-12 MG1655，在添加葡萄糖和酪蛋白氨基酸并含0%、40%和80% D2O的最低培养基中培养。监测生长曲线，并在菌株和条件特异性的、大致对应于所观察到的最大OD600中点的时点提取胞内代谢物。在正离子（ESI+）和负离子（ESI?）模式下进行非靶向液相色谱-质谱联用（LC-MS）分析，并使用成熟的代谢组学工作流程处理数据。统计分析包括主成分分析（PCA）、主成分-判别函数分析（PC-DFA）、双因素方差分析（two-way ANOVA，FDR<0.05）、倍数变化评估及代谢通路富集分析。 结果：生长曲线分析显示，随着D2O浓度的增加，生长抑制呈菌株依赖性且逐步增强，其中敏感型分离株受影响更为显著。代谢组学分析鉴定出跨越9条代谢通路的44种显著改变的代谢物，包括中心碳代谢、氨基酸降解、核苷酸代谢和氧化还原过程。敏感型菌株在谷胱甘肽和多胺代谢、烟酸酯和烟酰胺代谢、核苷酸代谢及中心碳相关过程中表现出显著的代谢物变化。相比之下，耐受型菌株在这些代谢物组中通常表现出较为温和的丰度变化。敏感型分离株的额外响应包括氨基酸、二肽和多胺相关代谢物的变化，而耐受型分离株通常表现出较为温和的响应。 结论：D2O暴露与E. coli中参与中心碳相关、氨基酸、核苷酸、烟酸酯和烟酰胺、谷胱甘肽和多胺代谢的代谢物丰度显著变化相关。敏感型与耐受型分离株的代谢轨迹显著分化，揭示了菌株特异性的适应策略。这些发现为同位素胁迫耐受性提供了生化基础，并增进了对细菌在重水暴露下重塑代谢机制的理解。

论文解读

一、研究背景与科学问题

稳定同位素探针技术（stable isotope probing，SIP）是连接微生物系统发育身份与其生态功能的强大工具，其中重水（D₂O）因水参与几乎所有生化反应而成为通用示踪剂。然而，氢被氘取代后会增强化学键强度并降低零点振动能，由此产生的初级和次级动力学同位素效应（kinetic isotope effects，KIEs）会减缓质子耦合电子转移和氢负离子转移等依赖氢的酶促过程，从而扰乱反应动力学和平衡，导致代谢失衡。尽管已知高浓度D₂O会抑制植物和动物生长，部分微生物（包括细菌、藻类和酵母）却能耐受较高水平的重水，但即使在这些相对耐受的物种中，D₂O暴露仍会造成生长减缓与代谢调整。

此前研究已从单酶体外水平量化了D₂O诱导的KIEs，系统层面的研究也揭示了中心碳通量在70% D₂O条件下向回补途径（anaplerosis）重定向、翻译机器下调和氧化还原酶上调等现象。然而，蛋白质组数据无法捕捉KIEs最先表达的小分子代谢物层面，D₂O胁迫下细菌代谢组的系统变化以及菌株特异性耐受的生化基础仍不清楚。为此，本研究利用非靶向液相色谱-质谱联用（LC-MS）代谢组学，系统表征D₂O敏感型与耐受型大肠杆菌菌株的全局代谢响应，以揭示同位素胁迫适应背后的代谢机制。

二、主要技术方法

研究人员选用5株尿路致病性大肠杆菌（uropathogenic E. coli，UPEC）分离株及参考菌株E. coli K-12 MG1655，在含0%、40%和80% D₂O的M9最低培养基（补充葡萄糖和酪蛋白氨基酸，MMCAA）中培养。采用Bioscreen C自动生长分析系统监测24小时生长曲线，据此将菌株分为D₂O敏感型（分离株99、133、134）和耐受型（分离株25、148及MG1655）。在菌株和条件特异性时点（大致对应最大OD₆₀₀的中点）采集样品，以?48°C预冷甲醇淬灭代谢后提取胞内代谢物，使用UHPLC-Q Exactive Orbitrap质谱仪在正离子（ESI+）和负离子（ESI?）模式下进行非靶向LC-MS分析。数据处理采用Compound Discoverer 3.3进行峰检测、对齐和氘标记分析，经质量过滤后使用MetaboAnalyst进行统计分析和通路富集，并利用MetExplore基于KEGG数据库构建E. coli K-12 MG1655代谢网络进行网络分析。

三、研究结果

3.1 重水培养基中UPEC分离株的生长曲线

生长曲线分析显示，6株大肠杆菌在含0%、20%、40%、60%和80% D₂O的MMCAA培养基中呈现两种响应模式。分离株99、133和134对D₂O敏感性较高，其OD₆₀₀增长随D₂O浓度增加而逐步延迟，部分菌株在20% D₂O时即出现生长变化，高浓度下最终OD₆₀₀值降低。相比之下，分离株25、148和参考菌株MG1655在不同D₂O浓度下生长曲线变化较小，OD₆₀₀增长发生时间更早，整体生长抑制效应较不显著。

3.2 LC-MS代谢谱分析

LC-MS分析在ESI+模式初始检测到7,008个代谢特征（含195个推定注释代谢物），ESI?模式检测到2,800个特征（含52个推定注释代谢物）。经数据预处理和质量过滤后，ESI+保留1,864个特征（含60个注释代谢物），ESI?保留1,337个特征（含37个注释代谢物）。PCA分析显示QC样品紧密聚集，表明数据质量高、重现性好。PC-DFA分析中，DF1轴对应D₂O暴露梯度，对照组位于负侧、40% D₂O居中、80% D₂O位于正端；DF2轴区分敏感型与耐受型菌株，且两表型在40% D₂O条件下的分离最为明显。双因素方差分析（FDR<0.05）鉴定出44种显著改变的MSI Level 2注释代谢物，其中11种在两种电离模式下均被检出。菌株水平敏感性分析显示，仅6.06%的代谢物-条件检验存在显著的菌株效应，支持表型水平合并分析的合理性。

3.2.1 LC-MS代谢物谱数据解读

中心碳代谢：在糖利用和糖磷酸代谢方面，敏感型菌株在0% D₂O时胞内葡萄糖丰度高于耐受型，40% D₂O时葡萄糖和磷酸烯醇式丙酮酸（PEP）显著下降，80% D₂O时葡萄糖部分回升。葡萄糖-1-磷酸（G1P）在40% D₂O时于敏感型中升高，80% D₂O时两组均下降。TCA循环相关代谢物方面，敏感型菌株在40% D₂O时琥珀酸升高、耐受型降低；苹果酸在两组中均随D₂O浓度增加而升高且敏感型更高；黄素单核苷酸（FMN）在耐受型40% D₂O时降幅较大。谷氨酸在耐受型40% D₂O时降低，敏感型80% D₂O时降低。整体而言，耐受型菌株在D₂O暴露下维持了相对稳定的NAD相关辅因子和TCA相关代谢谱，而敏感型菌株在琥珀酸、苹果酸和辅因子相关代谢物上出现更早的变化。

氨基酸和二肽相关代谢：敏感型菌株在40% D₂O时缬氨酰苯丙氨酸（FC=0.10）和丙氨酰酪氨酸（FC=0.08）急剧下降，同时游离缬氨酸升高；脯氨酰亮氨酸也仅在敏感型中于40% D₂O时显著下降（FC=0.20）。耐受型菌株则在40% D₂O时二肽丰度升高。80% D₂O时，耐受型中两种二肽均下降，敏感型较40%部分恢复但仍低于基线。支链氨基酸方面，苏氨酸在敏感型中始终较高且两组均随D₂O浓度增加而下降；异亮氨酸在敏感型40% D₂O时急剧升高（FC=1.51）而在耐受型中下降，可能反映敏感型菌株因辅因子可用性降低和电子传递链活性受损导致下游分解代谢受阻。

3.3 通路富集与网络分析

通路富集分析确定了9条由检测代谢物代表的通路类别，包括谷胱甘肽代谢、精氨酸和脯氨酸代谢、精氨酸生物合成、烟酸酯和烟酰胺代谢、嘌呤和嘧啶代谢、糖酵解、丙酮酸代谢、三羧酸（TCA）循环以及缬氨酸、亮氨酸和异亮氨酸降解。基于KEGG的MetExplore网络分析显示，大多数显著改变的代谢物定位于中心碳代谢（包括糖酵解和TCA循环）以及氨基酸和氧化还原相关通路。

四、讨论与结论

讨论：本研究的代谢组学发现与既往文献一致。Hochuli等人的研究曾表明D₂O暴露主要影响E. coli的中心碳代谢，抑制呼吸和TCA循环，而糖酵解和磷酸戊糖途径几乎不受影响，这与本研究发现能量产生和碳代谢相关通路受影响最显著的结果相符。Opitz等人的蛋白质组学研究显示D₂O暴露改变多种酶系统并下调翻译相关蛋白，为细胞对D₂O的响应提供了更广泛的背景。值得注意的是，通路富集结果需谨慎解读，因为非靶向代谢组学存在固有的代谢物覆盖不完整问题，且缺乏明确的背景代谢物集。通路归属并非互斥，中心代谢中间体可参与多条生化通路。研究团队将通路富集视为系统层面的组织框架，而非通路激活、抑制或代谢通量改变的直接证据。

结论：D₂O暴露与E. coli中参与中心碳相关、氨基酸、核苷酸、烟酸酯和烟酰胺、谷胱甘肽和多胺代谢的代谢物丰度显著变化相关。敏感型与耐受型分离株的代谢轨迹显著分化，揭示了菌株特异性的适应策略。这些发现为同位素胁迫耐受性提供了生化基础，增进了对细菌在重水暴露下重塑代谢机制的理解，对代谢工程、氘代化合物开发以及微生物学中SIP方法的优化具有潜在意义。