利用伪靶向SWATH/IDA-MRM进行耐药性大肠杆菌代谢组学分析：兼顾高覆盖率和精确性

1. 领域背景与文献

文献英文标题：Drug-resistantEscherichia colimetabolomics via pseudo-targeted SWATH/IDA-MRM: Bridging high coverage and precision；

发表期刊：未明确标注；影响因子：未公开；研究领域：微生物代谢组学，细菌抗生素耐药性研究。

领域共识：随着抗生素在临床和农业生产中的广泛应用，细菌耐药性问题不断加剧，耐碳青霉烯类、黏菌素的多重耐药大肠杆菌已对临床抗感染治疗安全构成严重威胁，而新型抗菌药物研发速度始终落后于细菌耐药性演化速度，因此明确细菌耐药过程中的分子机制，尤其是代谢层面的适应性改变，是开发新型抗菌策略的核心方向。已有研究证实，代谢重编程是细菌应对抗生素压力的核心适应性机制，抗生素会干扰细菌代谢稳态，而细菌的固有代谢状态反过来会调控抗生素敏感性，靶向调控代谢通路可重新敏化耐药菌对抗生素的反应，因此系统解析耐药菌的代谢组特征具有重要的应用价值。

当前代谢组学研究主要分为非靶向和靶向两类技术体系，两类技术都存在难以克服的缺陷，未解决的核心问题是现有技术无法同时实现代谢组的高覆盖检测和精准定量，难以满足耐药菌代谢组研究对低丰度差异代谢物检测的需求。非靶向代谢组学虽然可同时检测已知和未知代谢物，适合发现新的代谢特征，但存在动态范围有限、定量准确性差、低丰度代谢物易被掩盖、电离偏差大等问题；靶向代谢组学虽然定量精准、灵敏度高，但只能检测预先定义的代谢物，无法发现新的潜在生物标志物，且依赖标准品，不适合发现阶段研究。

本研究针对现有技术的核心缺陷，开发整合顺序窗口获取所有理论碎片离子/信息依赖采集-多反应监测（SWATH/IDA-MRM）的伪靶向代谢组学方法，同时结合三种不同的色谱分离模式实现不同理化性质代谢物的全覆盖检测，兼顾了非靶向的高覆盖性和靶向的定量精准性，填补了细菌代谢组研究的技术空白，并将该方法应用于携带不同耐药基因的耐药大肠杆菌代谢组分析，明确了耐药相关的代谢重编程特征，为后续代谢靶向抗菌策略开发提供了技术框架和理论基础。

2. 文献综述解析

作者对现有代谢组学研究技术按照技术类型分类展开评述，首先梳理了非靶向代谢组学的核心优势与局限性，非靶向代谢组学作为假说驱动的研究平台，优势在于可同时检测已表征和新型代谢物，能够覆盖更大范围的代谢物种类，适合发现阶段的研究，但局限性也十分突出，其动态范围受限，定量准确性不足，尤其对低丰度代谢物检测效果差，高浓度分析物易引发检测器饱和，掩盖低丰度的关键代谢物，同时复杂基质的干扰会降低离子检测可靠性，代谢物理化性质差异带来的电离偏差也会进一步影响痕量代谢物的检测稳定性。

随后作者评述了靶向代谢组学的技术特点，基于液相色谱-串联质谱的靶向代谢组学是代谢物定量的金标准，优势在于灵敏度高、通过多反应监测实现高特异性、动态范围宽、重现性好，能够满足精准定量的需求，但局限性在于其靶向属性限制了代谢物覆盖范围，只能检测预先设定的代谢物，不可避免会遗漏潜在的重要生物标志物，并且定量依赖标准品，从根本上限制了结构新型或未表征代谢物的检测，不适用于需要全面代谢谱分析的发现阶段研究。

通过对现有技术的评述，作者指出现有两类主流技术都无法同时兼顾高覆盖度和定量精准性，这一技术瓶颈严重限制了细菌耐药代谢组学的发展，本研究在此基础上提出整合非靶向和靶向技术优势的伪靶向策略，创新点在于首次结合SWATH和IDA两种采集模式的优势，整合构建高置信度的多反应监测离子对文库，同时结合反相、亲水相互作用、金属敏感分析物三种液相色谱分离模式，实现了不同理化性质代谢物的全覆盖，在保证高覆盖度的同时兼顾了定量精准性，解决了现有技术的核心缺陷，为细菌代谢组学研究提供了新的技术范式。

3. 研究思路总结与详细解析

本研究的整体研究目标为开发同时具备高代谢物覆盖度和定量精准性的伪靶向代谢组学方法，核心科学问题是如何整合非靶向和靶向代谢组学的技术优势，解析携带不同耐药基因大肠杆菌的代谢重编程特征，技术路线遵循“技术优化-方法构建-性能验证-样本应用”的闭环逻辑，首先优化色谱分离条件，比较不同质谱采集模式的性能，随后整合两种采集模式构建SWATH/IDA-MRM方法，验证方法的分析性能，最后应用该方法解析耐药大肠杆菌的代谢组特征，明确耐药相关代谢重编程。

3.1 实验菌株与样本制备

实验目的为制备携带不同耐药基因的大肠杆菌样本和质控样本，为后续方法开发和代谢组分析提供合格的实验材料。方法细节：实验菌株包括大肠杆菌DH5α和ATCC25922来源的不同基因型菌株，分别为携带空载体的对照株、携带mcr-1、blaNDM-1、blaNDM-5的单耐药株，以及同时携带mcr-1/blaNDM-1的双耐药株；所有菌株在37℃下于LB培养基中培养，收获对数生长期（OD600=0.3±0.05，n=6）的菌体以保证代谢活性一致，所有培养和提取过程平行进行以降低批次效应；采用预冷（-60℃）甲醇/乙二醇（45:55，体积比）快速淬灭菌体，沉淀后用0.85%氯化钠溶液洗涤两次，采用乙醇/水（75:25，体积比）裂解细胞，离心后取上清浓缩，用50%乙腈复溶，合并所有菌株的样本制备混合质控样本，每检测6个样本插入一个质控样本以保证系统稳定性。文献未提及具体实验产品，领域常规使用微生物培养、代谢提取类试剂与仪器。

3.2 色谱条件优化与不同采集模式性能比较

实验目的为优化金属敏感代谢物的色谱分离条件，比较全扫描、IDA、SWATH三种质谱采集模式的性能，明确不同模式的互补性。方法细节：构建整合反相液相色谱、亲水相互作用液相色谱、金属敏感分析物液相色谱三种分离模式，结合正负电离模式共得到6种分析条件，分别用于分离不同极性的代谢物，其中反相分离非极性和中等极性化合物，亲水相互作用色谱保留高极性代谢物，金属敏感分析物液相色谱提升含磷酸和多羧酸代谢物的检测灵敏度；考察流动相中不同浓度（0、1、2、5、10μM）的medronic acid对金属敏感代谢物峰形和信号强度的影响，确定最优添加浓度；随后分别采用三种采集模式采集质控样本数据，比较不同模式的MS1特征数量、质荷比匹配度、MS2覆盖度和谱图质量。结果解读：色谱优化结果显示，当medronic酸浓度为5μM时，乙酰辅酶A、GDP、NADH、GTP等金属敏感代谢物的峰形和信号强度最优，浓度过高（10μM）会引发离子抑制，过低无法改善峰形（文献未明确提供具体数值，基于图表趋势推测）；不同采集模式比较结果显示，全扫描模式不进行前体离子选择，检测到的MS1特征数量最多，质荷比匹配度最高，但无法获得MS2谱图，MS2覆盖度为0；IDA模式需要分配采集时间给高丰度离子的MS2采集，MS1特征数量少于全扫描，仅能对高丰度离子采集MS2，MS2覆盖度有限；SWATH模式采用分段采集方式，MS1特征数量最少，但实现了100%的MS2覆盖，且MS2谱图质量更高，RDP得分大于0.8的代谢物比例显著高于IDA模式；韦恩图分析显示三种采集模式的代谢特征存在部分重叠，每种模式都能检测到独有的特征，具备互补性，整合多种模式可有效提升代谢组覆盖度。

实验所用关键产品：Waters的Acquity I-Class Plus超高效液相色谱系统，SCIEX的X500B四极杆飞行时间质谱仪、Triple Quad 7500三重四极杆质谱仪。

3.3 SWATH/IDA-MRM方法构建与性能验证

实验目的为整合SWATH和IDA的特征离子对，构建伪靶向SWATH/IDA-MRM方法，并系统验证方法的分析性能。方法细节：分别从SWATH和IDA的原始数据中提取MS1和MS2信息，筛选特征离子对，随后按照质荷比偏差≤5ppm、保留时间偏差≤0.05min的标准匹配合并两种方法得到的离子对，然后采用三重四极杆质谱对合并后的离子对进行实验验证，去除信噪比<3或无目标色谱峰的离子对，得到最终的高置信度离子对文库；采用梯度稀释（1倍到1024倍稀释）的质控样本，比较全扫描、IDA、SWATH和SWATH/IDA-MRM四种方法的线性、灵敏度、重现性，通过主成分分析比较不同方法的样本分离效果。结果解读：SWATH模式得到5694个特征多反应监测离子对，IDA模式得到1997个特征离子对，两种模式存在部分重叠，最终验证后得到3529个高置信度离子对，共注释得到2270个代谢物，其中430个代谢物通过标准品验证，符合代谢组学标准倡议一级注释标准；性能验证结果显示，SWATH/IDA-MRM方法中超过60%的代谢物线性决定系数R2>0.95，超过80%的代谢物R2>0.8，其中金属敏感分析物负电离模式中91.83%的代谢物R2>0.8，70.26%的R2>0.95，线性表现显著优于其他三种方法；灵敏度方面，SWATH/IDA-MRM在所有浓度范围都表现出更优的检测能力，对低丰度代谢物的检测优势尤为明显；重现性方面，SWATH/IDA-MRM检测到的代谢物中97.1%±2.0%的相对标准偏差<30%，显著优于常规非靶向方法；主成分分析显示，SWATH/IDA-MRM可实现清晰的组间样本分离，质控样本聚类紧密，方法稳定性和重现性优异。

3.4 耐药大肠杆菌代谢组特征分析

实验目的为应用构建完成的SWATH/IDA-MRM方法，解析携带不同耐药基因大肠杆菌的代谢重编程特征。方法细节：每组耐药株和对应的对照株均设置6个生物学重复，采用构建好的方法进行代谢组检测，原始数据经过归一化后进行主成分分析和正交偏最小二乘判别分析，差异代谢物筛选标准为倍数变化>2或<0.5，学生t检验P<0.05，正交偏最小二乘判别分析的变量重要性投影得分>1，对差异代谢物进行KEGG通路富集分析，明确耐药相关的核心代谢通路改变。结果解读：主成分分析结果显示，耐药株和敏感株在多元空间中完全分离，不同耐药株形成独立的聚类，说明耐药株在获得耐药基因后发生了系统性的代谢组重编程；差异分析结果显示，SWATH/IDA-MRM可以检测到非靶向方法遗漏的多种低丰度差异代谢物，包括N-氨甲酰-L-天冬氨酸、琥珀酸、CDP、尿苷二磷酸、三磷酸胞苷、氨甲酰磷酸等，这些代谢物均参与细菌核心代谢通路，对细菌耐药适应具有重要作用；通路富集分析显示，所有耐药株的差异代谢物主要富集在氨基酸代谢、碳水化合物代谢、能量代谢、聚糖生物合成、核苷酸代谢五大核心代谢模块；不同耐药基因型菌株存在共同和特异性的代谢改变，携带mcr-1的菌株和双耐药菌株的O抗原生物合成、氨基糖代谢发生显著改变，三羧酸循环中间产物水平升高，ppGpp、褪黑素水平显著升高，肌苷一磷酸水平显著降低；携带blaNDM-1和blaNDM-5的菌株糖酵解、糖异生、磷酸戊糖途径、嘌呤代谢发生显著改变，还原型谷胱甘肽水平升高，氧化型谷胱甘肽水平降低；所有耐药株都存在核苷酸代谢和能量代谢的上调，为细菌在抗生素压力下的增殖和存活提供能量和核苷酸前体。

4. Biomarker研究及发现成果

本研究基于开发的伪靶向代谢组学方法，筛选得到一系列与大肠杆菌抗生素耐药表型相关的代谢生物标志物，筛选逻辑为：构建整合高覆盖和精准定量的技术方法→对携带不同耐药基因的耐药株和敏感对照株进行代谢组检测→差异分析结合统计检验获得耐药相关差异代谢物→结合已有功能研究明确候选生物标志物的应用价值，整个逻辑链条完整，从方法开发到临床相关菌株验证，保障了生物标志物的可靠性。

本研究的生物标志物来源于大肠杆菌全细胞代谢物，采用伪靶向代谢组学进行定量检测，通过六次生物学重复的统计检验筛选差异代谢物，最终得到多个不同耐药基因型特异性的候选生物标志物。其中，ppGpp在所有耐药株中均显著升高（P<0.05，n=6），是所有研究耐药株共同的差异代谢物；mcr-1阳性耐药株中，褪黑素水平显著升高（倍数变化>2，P<0.05，n=6），肌苷一磷酸水平显著降低（倍数变化<0.5，P<0.05，n=6）；blaNDM阳性耐药株中，氧化型谷胱甘肽水平显著降低（P<0.05，n=6）。

核心成果方面，本研究首次系统利用整合伪靶向代谢组学方法解析了携带不同临床关键耐药基因的大肠杆菌的代谢特征，明确了多个具有临床应用潜力的代谢生物标志物：ppGpp作为细菌严谨应答的核心信号分子，可调控细菌转录、翻译、DNA复制等过程，重塑细菌代谢促进耐药适应，本研究证实ppGpp在耐药株中显著积累，支持其作为抗菌治疗的新型靶点，统计学结果符合研究要求（P<0.05，n=6）；对于mcr-1介导的黏菌素耐药，褪黑素水平升高和肌苷一磷酸水平降低是特征性代谢改变，已有研究证实外源性补充褪黑素可逆转黏菌素耐药，外源性补充肌苷一磷酸可直接杀死mcr-1阳性耐药菌，这两种代谢物可作为mcr-1阳性耐药菌的候选诊断标志物和治疗靶点；对于blaNDM介导的碳青霉烯耐药，氧化型谷胱甘肽水平降低是特征性改变，已有研究证实外源性补充氧化型谷胱甘肽可恢复耐药菌对美罗培南的敏感性，该代谢物可作为blaNDM阳性耐药的候选生物标志物和治疗靶点；本研究获得的所有差异代谢物都经过严格的统计检验，为后续耐药诊断和代谢靶向抗菌策略开发提供了可靠的候选分子，也为细菌耐药机制研究提供了新的方向。