《JOURNAL OF MASS SPECTROMETRY》：Decoding Biothreats With FT-ICR-MS: Metabotyping of Bacillus cereus Spores Through Untargeted Metabolomics
编辑推荐：
基因组测序是微生物鉴定的金标准，但无法区分遗传特征高度相似的近缘物种。蜡样芽孢杆菌（B. cereus）和细胞毒性芽孢杆菌（B. cytotoxicus）即为典型代表——二者同属蜡样芽孢杆菌群（B. cereus group），却具有截然不同的致病谱。B. c
基因组测序是微生物鉴定的金标准，但无法区分遗传特征高度相似的近缘物种。蜡样芽孢杆菌（B. cereus）和细胞毒性芽孢杆菌（B. cytotoxicus）即为典型代表——二者同属蜡样芽孢杆菌群（B. cereus group），却具有截然不同的致病谱。B. cereus作为已知的机会性致病菌，可引发肠道及肠外感染，产生组织破坏性外酶溶血素BL（HBL）、非溶血性肠毒素（NHE）和细胞毒素K-2（CytK-2），以及质粒编码的催吐毒素蜡样环肽（Ces）；而耐热种B. cytotoxicus则主要产生与严重腹泻暴发相关的细胞毒素K-1（CytK-1）。尽管存在这些差异，单纯基因组数据常不足以明确区分两种细菌。此外，这些细菌在环境样品中主要以高抗性孢子形式存在，进一步增加了传统方法鉴定的难度，阻碍准确分类和风险评估——这一问题在生物战背景下尤为关键。蜡样芽孢杆菌群包含若干携带与炭疽芽孢杆菌（B. anthracis）相似遗传元件的菌株，部分B. cereus菌株甚至可通过获得毒性质粒引发类炭疽感染。若在这些高抗性细菌的孢子水平发生误判或延迟鉴别，在生物袭击或意外事件中可能导致严重后果，此时快速准确检测对于有效及时响应至关重要。本研究表明，基于傅里叶变换离子回旋共振质谱（FT-ICR-MS）的非靶向代谢组学能够实现B. cereus与B. cytotoxicus的鉴别，同时揭示孢子特异性代谢组学特征，为改进微生物法医鉴定和生物监测策略提供关键依据。
生物战剂（bioweapon）的使用历史几乎与人类文明同样悠久，而根据规模和意图可分为生物战（biowarfare）、生物恐怖（bioterrorism）和生物犯罪（biocrime）。并非所有致病微生物都适合作为生物战剂，致病性或毒性本身并不决定其武器化潜力，组织损伤能力、可获得性及向宿主有效递送的效率才是评估关键。美国疾病预防控制中心（CDC）将生物战剂分为三类：A类因易于传播、致死率高、可造成重大公共卫生和社会扰乱而风险最高；B类传播较易、致病率中等、致死率低；C类则包含未来可被工程化大规模传播的新发病原体。炭疽芽孢杆菌（B. anthracis）和蜡样芽孢杆菌生物变种炭疽（B. cereus biovar anthracis）被列为A类生物战剂，纳入美国精选制剂和毒素清单（Select Agents and Toxins List）的Tier 1子集。尽管B. cereus本身未被赋予"充分证实的病原地位"，但其作为食源性机会致病菌可导致催吐和腹泻综合征，也能引起免疫缺陷患者的严重致死性感染，且与B. anthracis同属蜡样芽孢杆菌广义种（B. cereus s.l.）复合群——该复合群包含从重要病原体到农业工业用菌的众多谱系。作为产孢细菌，这些物种能适应极端环境条件，在有利条件恢复时可快速生长繁殖，使细菌孢子成为生物武器应用的金标准。
然而，蜡样芽孢杆菌广义种的分类长期存在挑战和争议。1998年基于表型和经典分子方法将该群分为6个物种；全基因组测序（WGS）和平均核苷酸同一性（ANI）的采用又陆续识别出更多物种，2013年新增B. toyonensis和B. cytotoxicus，2016年新增B. wiedmannii，2017年再增9个物种，使复合群扩展至至少18个公认物种。但ANI分类在95%–96%阈值下的应用对所选截断值和参考基因组敏感，常导致重叠的基因组物种聚类和不一致的分类结果。为此提出了多种命名系统，NCBI分类学虽提供广泛使用的分级框架，但作为动态资源易受错误分类影响。实践中，表型特征分析仍是初步分类的主要步骤，基因组方法用于解决模糊情况，但表型特征往往分辨率不足，基因组分类则存在变异性和可重复性问题，加之高频水平基因转移进一步增加了物种鉴定的复杂性。
在孢子水平开发替代鉴定方法尤为迫切。某些芽孢杆菌属（B. anthracis、B. cereus和B. thuringiensis）的孢子难以通过光镜或电镜区分，同一物种内菌株间几乎相同。核心组分（吡啶二羧酸、小酸溶性蛋白[SASPs]和钙离子）及孢子衣蛋白虽保守，但提取困难耗时。孢子形成过程中，营养细胞经历广泛代谢重编程，形成代谢休眠但非代谢无活性的特殊实体，保留高度特异性的代谢组。这种战略性重组支持长期存活，并确保萌发初期代谢所需的内源性资源。然而孢子代谢组的表征远未完善，却代表着孢子生物分型的强有力工具——细微遗传变异可转化为显著的下游代谢改变，反映生物体的实际生理状态。
本研究旨在利用傅里叶变换离子回旋共振质谱（FT-ICR-MS）对五种B. cereus狭义种（B. cereus s.s.，包括两个野生株和三个参考株）及一种B. cytotoxicus参考株的孢子代谢组进行表征，以实现物种和菌株水平的鉴别。研究样本来源于葡萄牙陆军生物防御实验室（Laboratório de Defesa Biológica）保藏的两株野生B. cereus菌株（分别分离自人体呕吐物[J3]和食物[J271]）以及购自NCTC培养物保藏中心的三株参考菌株[J630/ATCC 10987，源自乳制品；J648/F4810/72，催吐性疾病相关；J649/F528/94，腹泻性疾病相关]和一株B. cytotoxicus参考菌株[J650/NVH0391-98，致死性腹泻相关，分离自食物]。
关键技术方法方面，研究人员采用70%:30%（v/v）甲酸:乙腈混合溶剂同步进行孢子灭活和代谢物提取，验证完全灭活后于?80°C保存待测。分析采用15 T SolariX XR FT-ICR质谱仪（Bruker Daltonics）进行电喷雾电离直接进样，流动速率2.0 μL/h，累积100次瞬态信号，质量范围m/z 90–1500，分辨能力足以分离m/z 340.34850与340.34611的两个峰。数据处理使用Bruker Data Analysis 5.0软件进行峰查找，后续通过基于Python和Metabolinks的自制软件进行m/z与峰强度列表对齐、分子式归属、化合物注释和统计分析，包括缺失值填补（数据集最小强度值的1/5）、总和归一化、广义对数（glog）转换和均值中心化缩放。化合物注释采用人类代谢组数据库（HMDB）和微生物代谢组数据库（MiMeD）。富集分析采用超几何检验评估HMDB化学亚类的过代表性，多重检验校正采用Benjamini–Hochberg错误发现率（FDR）法；无监督分析包括主成分分析（PCA）和基于Bray–Curtis距离、单连接的层次聚类分析（HCA）；有监督分析采用500棵树的随机森林（RF）模型（随机种子42），以Gini重要性值排序，保留前2%最高得分化合物，并通过单因素方差分析（ANOVA）评估统计显著性。
研究结果部分，"孢子代谢谱分析"显示，各菌株谱图差异明显，表明代谢特异性显著。B. cereus野生株J3和J271因分离来源不同，检测到最多m/z峰（分别为502和178个）；B. cytotoxicus参考株J650较B. cereus参考株显示更高代谢多样性。Van Krevelen图分析揭示脂质、氨基酸/肽、碳水化合物、萜类和芳香氧化合物五大化学类群构成孢子代谢组核心架构，其中碳水化合物类群稀疏分布。野生株在脂质、氨基酸/肽和萜类区域代谢多样性显著高于参考株，而J648分子多样性最低且无菌株特异性代谢物，提示其孢子代谢组较为精简。B. cytotoxicus J650与食物来源的J271共享最多m/z峰（212个），暗示分离来源对代谢组组成的影响。
"脂质和萜类的功能意义"部分，过代表性分析在463个共有m/z峰（注释后保留约12%）中识别出脂肪酸及结合物、脂肪酰胺和脂肪酸酯三个显著富集的亚类。除典型膜结构相关脂肪酸（如C15:0、C16:0、C17:0、C18:0及其支链异构体）外，还检测到乳产品相关的不寻常脂肪酸如 vaccenic acid，支持孢子形成环境对孢子代谢组的影响。特别值得注意的是13-docosenamide（erucamide，芥酸酰胺）和多种甲基十四烷酰肉碱异构体的识别：erucamide作为抗菌长链脂肪酸酰胺，在致病性参考株J648和J649中丰度较高，提示其可能作为调节宿主-病原体相互作用的早期信号分子；（3–13）-甲基十四烷酰肉碱则在野生株J3中强度最高，可能与呕吐基质中唾液/胃液分泌物及肉/乳制品消化产生的肉碱丰富来源有关。总体而言，野生株脂肪酸及结合物相对丰度显著高于参考株，腹泻相关菌株（J649和J650）显示相似脂质丰度特征，而食物来源和呕吐来源菌株呈现可区分脂质和酰基肉碱特征。
"支持芽孢杆菌孢子鉴别的代谢特征"部分，PCA和HCA揭示三个主要聚类：代谢分化最高的野生呕吐源J3株；紧密聚类的B. cereus参考株；以及包含B. cytotoxicus和野生食物源J271的第三组。第三主成分（PC3，12.6%方差）可将野生株与参考株分离，而第一、二主成分（55.8%方差）反映的分离来源是主要驱动因素。随机森林模型识别出50个Gini重要性最高的m/z质量，其中15个获得注释，仅3个符合细菌/孢子生物学框架，分别属于脂肪酸酯和二萜类。具体而言，肉碱脂酰酯表现出普遍性作用：十四烷酰肉碱（C14）在各菌株中均检测到，J3中水平特高；而十五烷酰（C15）和十七烷酰（C17）肉碱在致病性B. cereus参考株和B. cytotoxicus中缺失，但在非致死乳制品相关参考株J630和野生食物源J271中高度富集，提示较长链酰基肉碱（C15–C17）可能与非致病性或环境适应菌株相关，而较短链物种（C14）广泛分布，酰基肉碱脂肪酸链长度可能反映菌株的代谢特征并与致病潜力相关。此外，研究还鉴定了两类新型氧合C25 β-apo-12′-胡萝卜素异构体（C25H34O3），在乳制品相关参考株J630和腹泻参考株J649中特异性富集，可能是食物基质中膳食β-胡萝卜素摄入后经类胡萝卜素加氧酶氧化切割的产物，可能增强孢子对紫外线和氧化应激的抗性。
讨论与结论部分，研究强调细菌物种固有相似性带来的生物分型挑战，特别是在预防和控制措施实施中。B. cereus群作为产孢细菌属，因其广泛的致病和非致病物种、传统方法难以遗传区分而尤为突出。本工作证明，基于FT-ICR-MS的非靶向代谢组学显示蜡样芽孢杆菌广义种孢子具有由孢子形成环境塑造的核心代谢组，同时保留足够的代谢特异性以实现物种和菌株鉴别。脂质、氨基酸/肽、碳水化合物、萜类和芳香氧化合物在各物种中一致检出，提示存在支持孢子抗性的共享代谢框架。野生分离株J3和J271表现出最高代谢丰富度和最复杂化学谱，而参考株J648多样性最低且缺乏菌株特异性代谢物；B. cytotoxicus J650显示中间代谢谱，与食物来源株J271聚类更近；J649因脂质谱特征与B. cytotoxicus聚类，J630则以中长链酰基肉碱升高和独特二萜类特征为标志。这些特征共同实现了基于可重复代谢特征的B. cereus菌株和B. cytotoxicus鉴别。
富集分析强调脂质不仅是膜结构组分，还作为环境感知和应激防御介质。外源性脂质掺入孢子结构可能反映原始分离环境的代谢印记，为区分自然来源菌株和实验室产生菌株提供潜在指标。此外，中长链酰基肉碱的检测及其与致病性的关联提示了毒力评估的可能策略。总体而言，这些发现证明了非靶向代谢组学作为细菌孢子生物分型前景策略的潜力，所识别的特定代谢物组合为孢子分类提供了坚实基础，并为向更广泛可用的质谱仪器进行方法学转移铺平了道路。
