《Current Microbiology》：Genomic Insights into a New Burkholderia cenocepacia Sequence Type Linked to Cepacia Syndrome in Cystic Fibrosis
编辑推荐：
洋葱伯克霍尔德菌（Burkholderia cenocepacia）是囊性纤维化（cystic fibrosis，CF）患者中的重要病原体，可导致慢性肺部感染；在严重情况下，可引发洋葱伯克霍尔德菌综合征（cepacia syndrome，CS），其特征为坏死性
洋葱伯克霍尔德菌（Burkholderia cenocepacia）是囊性纤维化（cystic fibrosis，CF）患者中的重要病原体，可导致慢性肺部感染；在严重情况下，可引发洋葱伯克霍尔德菌综合征（cepacia syndrome，CS），其特征为坏死性肺炎、菌血症和高死亡率。流行性谱系如 B. cenocepacia J2315（ET12）和 ST32 具有高度传播性，并且常与 CS 相关。本研究对1例发生 CS 的儿童 CF 患者痰液和血液样本中分离获得的5株 B. cenocepacia 分离株进行测序，并采用多位点序列分型（multilocus sequence typing，MLST）、系统基因组学（phylogenomics）、基因组相似性指标和比较基因组学开展计算机模拟分析（in silico analysis），以鉴定毒力决定子、耐药基因及移动遗传元件。研究鉴定出一种新的序列型（sequence type，ST2424），该序列型不隶属于任何已知克隆复合群。系统基因组学分析显示，各分离株之间具有高度遗传相关性，提示该感染符合单克隆感染并伴随宿主体内微进化（within-host microevolution）的特征。比较分析揭示出不同的遗传谱：aaiB（T6SS）存在于全部分离株中，而 narG 和 manCcore 仅见于血流分离株。研究还检测到获得性耐药决定子，包括 blaOXA?1043、sul1、aac(6′)-IIc 和 aadA1，而这些基因在参考菌株 B. cenocepacia J2315 和 ST32 中均未发现。此外，还鉴定到移动遗传元件，包括属于 IS3、IS5 和 IS200/IS605 家族的插入序列，以及一个推定可接合质粒，其编码Ⅳ型分泌系统（type IV secretion system，T4SS）的 Tra/Trb 亚基。这些发现描述了一种与致死性 CS 相关的新型 B. cenocepacia 序列型，并强调该物种具有高度基因组可塑性，进一步凸显了开展基因组监测以识别新兴变异株并改善 CF 临床管理的重要性。
该论文发表于《Current Microbiology》，围绕囊性纤维化（cystic fibrosis，CF）患者中洋葱伯克霍尔德菌（Burkholderia cenocepacia）的致病与进化问题展开。Burkholderia cepacia complex（Bcc，洋葱伯克霍尔德菌复合群）是一类非发酵革兰阴性杆菌，在CF患者中可引起持续性呼吸道感染，并与肺功能快速恶化及死亡风险升高密切相关。其临床危害不仅源于对多种抗菌药物、消毒剂的固有耐受，还与生物被膜形成、生态适应范围广及基因组高度可塑性有关。B. cenocepacia 尤其值得关注，因为该物种与更差的预后及洋葱伯克霍尔德菌综合征（cepacia syndrome，CS）密切相关。CS 是一种急性坏死性肺炎综合征，常伴高热、菌血症及迅速进展的呼吸衰竭，病死率较高。既往研究提示，J2315（ET12）和 ST32 等流行株具有较强传播能力和毒力，但对于新发序列型在宿主体内持续、播散及致死性结局中的基因组特征，仍缺乏深入认识。因此，开展本研究的必要性在于：通过对CS病例分离株进行全基因组比较，明确其系统发育位置、毒力与耐药构成以及移动遗传元件特征，从而为理解严重CF感染的分子基础及强化基因组监测提供依据。
研究人员纳入1例来自巴西 IFF/Fiocruz 囊性纤维化参考中心的11岁女性CF患儿，其基因型为 F508del 纯合，并在出现 CS 后进展为呼吸衰竭和死亡。研究对象为住院期间获得的5株 B. cenocepacia 临床分离株，包括1株痰液分离株和4株血液分离株。主要技术方法包括：采用基质辅助激光解吸电离飞行时间质谱（matrix-assisted laser desorption/ionization time-of-flight mass spectrometry，MALDI-TOF MS）进行菌种鉴定；应用 Illumina MiSeq 平台完成全基因组测序；结合多位点序列分型（MLST）、数字DNA-DNA杂交（digital DNA–DNA hybridization，dDDH）、平均核苷酸一致性（average nucleotide identity，ANI）、平均氨基酸一致性（average amino acid identity，AAI）和基于单核苷酸多态性（single-nucleotide polymorphism，SNP）的系统发育分析确定分类与亲缘关系；进一步利用毒力因子数据库、耐药数据库及多种生物信息学平台解析毒力组（virulome）、耐药组（resistome）和移动组（mobilome）。
在“Identification by Mass Spectrometry”部分，研究显示5株分离株经 MALDI-TOF MS 鉴定得分为2.33–2.42，均达到种水平鉴定可信范围，全部被判定为 B. cenocepacia。该结果说明后续全基因组分析具有明确的物种学基础。
在“Whole-Genome Sequencing Analysis”部分，研究人员基于 recA 基因序列分析确认全部分离株均属于 B. cenocepacia genomovar IIIA。通过 PubMLST 进行MLST分析后，在 gltB 位点发现新等位基因 gltB_1015，进而定义出新的序列型 ST2424，且该型别未归属于任何已知克隆复合群。进一步的 dDDH、ANI 和 AAI 分析表明，5株临床分离株与参考菌株 B. cenocepacia J2315 及 ST32 均保持很高相似性，但系统基因组树和TYGS分析显示，这5株菌自身聚为单一分支，并与 J2315、ST32 处于相近但可区分的位置。该结果证明本研究对象属于 IIIA 谱系中的一个新序列型，同时在物种层面上与已知流行株保持紧密亲缘关系。
在同一部分的SNP系统发育分析中，研究人员发现5株临床分离株形成单系群。分离株间两两SNP差异为16–253个，其中3株血流分离株之间差异最小，仅16–32个SNP；痰液分离株与血流分离株之间差异为82–112个SNP；另1株血流分离株3415与其余血流分离株差异较大，为183–253个SNP。与参考株相比，临床株的SNP距离则明显增大。该结果支持单克隆感染在宿主体内持续存在并发生微进化的判断，也提示从肺部到血流的播散过程中伴随有限但可检测的遗传分化。
在“Distribution and Comparison of Virulence Genes”部分，研究共鉴定出100个与毒力相关的基因，按功能归为11类，其中90个构成所有基因组共享的保守核心。主要类别包括黏附相关因子、抗吞噬相关的 wcb 操纵子、侵袭/趋化相关基因群、群体感应系统以及 T3SS/T6SS-1 分泌系统成分。这说明本研究分离株整体上保留了 B. cenocepacia 典型而完整的毒力框架。进一步比较发现，aaiB 为所有临床分离株共有，而在参考流行株中未检出；narG 和 manCcore 则仅存在于4株血流分离株中。由此可见，不同解剖部位来源的分离株在部分毒力与适应性相关基因上存在差异，这为理解血流播散阶段的代谢适应和宿主环境适应提供了直接基因组证据。
在“Resistome Analysis”部分，ResFinder 分析显示，5株临床分离株具有一致的获得性耐药谱，均携带 aac(6′)-IIc、aadA1 和 sul1，而这些基因在 J2315 和 ST32 中均不存在。CARD/RGI 进一步检出 aadA、blaOXA?1043、penB-4 和 qacEΔ1 等决定子，其中 blaOXA?1142 仅见于 J2315。PATRIC AMR 模块共预测43个耐药基因，分布于10类机制之中，包括抗生素靶位、外排泵、调控、通透性、细胞壁修饰、保护、激活和靶替代等，多数在临床株与参考株之间是保守的。最显著的差异集中在抗生素灭活机制相关基因，尤其是 aac(6’)-Ib、aac(6’)-II 和 aadA 在所有临床分离株中均存在。该部分结果表明，ST2424 除具有 B. cenocepacia 核心基因组固有的多重耐药背景外，还额外获得了多个与氨基糖苷类和磺胺类耐药相关的基因，使其耐药谱更具临床挑战性。
在“Mobilome Analysis”部分，PHASTEST 在所有临床分离株中均预测到完整前噬菌体区域，但这些区域主要包含典型病毒模块，未见明确毒力或耐药基因。IslandViewer 预测每个基因组含132–152个基因组岛（genomic island，GI），总长度约2.10 Mb，分布于3条染色体。人工校正证实所有分离株的GI中稳定存在 dihydropteroate synthase type 2 和 ANT(3″)-Ia（aadA 家族）等耐药决定子，提示部分临床相关耐药因子整合于水平获得区域。插入序列方面，研究检出 IS3、IS5 和 IS200/IS605 三个家族成员；在多数分离株中，携带 ISPa18 的 contig 同时包含 aadA1、qacE 和 sul1，3415中还与 aac(6′)-IIc 位于同一 contig。该结果提示插入序列可能参与耐药元件的聚集与转移。质粒分析则在全部分离株中发现一个推定可接合质粒，平均长度约80 kb，编码 MOBF 型松弛酶、MPF_F 型交配形成系统，以及 traI、traD 和多个 Tra/Trb 相关成分。该发现表明，这些菌株具有潜在的接合转移能力，可能促进适应性性状在菌株间乃至种间传播。
讨论部分围绕新序列型 ST2424 的临床与生物学意义展开。研究人员指出，该型别来源于1例发生致死性 CS 的CF患儿，属于 B. cenocepacia IIIA，但不同于既有克隆复合群。基于相似性指标与系统树拓扑，ST2424 与 ET12 谱系代表株 J2315 关系更近，这与 IIIA 谱系具有较高流行潜能的既往认识相一致。SNP分析支持单一感染谱系在宿主体内持续存在并发生微进化，提示肺部环境压力和长期抗菌治疗可能促进遗传分化。毒力层面，绝大多数已知决定子高度保守，但血流分离株特有的 narG 和 manCcore 以及临床分离株共有的 aaiB，显示出与厌氧代谢、脂多糖（lipopolysaccharide，LPS）生物合成和 T6SS 相关的特征差异。耐药层面，研究强调其核心耐药背景之外还叠加了获得性耐药决定子，尤其是 sul1、blaOXA?1043 及多种氨基糖苷修饰酶基因，提示新序列型可能具有更强的治疗逃逸能力。移动组分析进一步表明，该菌基因组中存在大量可移动和不稳定区域，结合插入序列与可接合质粒的存在，共同体现出显著的基因组可塑性。作者同时指出，本研究全部结论均基于 in silico 分析，且样本仅来自单一患者，因此相关决定子的功能及其普遍性仍需更大样本和纵向研究加以验证。
研究结论部分可概括为：本研究为来源于1例发生 CS 并进展为致死性菌血症的CF患者的 B. cenocepacia genomovar IIIA 新序列型 ST2424 提供了重要基因组学认识。基因组分析显示，该序列型与流行性 ET12 谱系高度相似，同时具有兼具保守性与特异性的毒力组、耐药组和移动组特征，提示基因组可塑性和水平基因转移可能共同促进其持续定植与系统性播散。需要强调的是，本研究所有分析均为计算机模拟分析，仍需实验验证以确认所鉴定决定子的功能作用；此外，由于样本仅来自单一患者，研究结果的可推广性受到限制。未来应在更大CF队列及纵向采样框架下评估 ST2424 相关模式的一致性及其临床意义。总体而言，该研究突出了持续开展基因组监测的重要性，有助于发现新兴变异株，并为CF患者感染控制与临床管理策略优化提供支持。
