鲍曼不动杆菌（Acinetobacter baumannii）中东地区群体基因组学与机器学习揭示独特区域基因组指纹

《Microorganisms》：Comparative Population Genomics and Machine Learning Reveal a Distinct Regional Genomic Fingerprint of Middle Eastern Acinetobacter baumannii

　　

鲍曼不动杆菌（Acinetobacter baumannii）是重要多重耐药（MDR）医院获得性病原体，其中东地区种群结构与进化动态仍知之甚少。研究人员比较中东与全球分离株以识别区域基因组特征。研究人员分析199个高质量A. baumannii基因组（100个中东、99个全球），采用泛基因组重建、核心基因组系统发育、耐药组与毒力组分析、全基因组关联研究（GWAS）及机器学习。泛基因组含10971个基因，其中1102个为中东西亚分离株特有。与全球株相比，中东基因组携带较少抗菌药耐药（AMR）基因，但移动相关基因丰度更高。一个以序列型2（ST2）为主的不同区域簇占中东分离株66%。GWAS识别区域相关标记，尤其ISAba1变体与tufB基因。随机森林（Random Forest）分类器中东与全球株准确率为79.9%。两组毒力基因组成相似，耐药基因负担随时间稳定。中东A. baumannii分离株呈现独特区域基因组指纹，特征是附属基因组保守、基因组异质性降低、ISAba1相关标记富集及移动相关基因含量增加。这些特征使区域种群区别于全球播散株，尽管两队列中ST2/国际克隆2（IC2）谱系流行率相似。

研究背景方面，鲍曼不动杆菌（Acinetobacter baumannii）是全球重点监控医院病原体，具备强环境存活力、生物被膜形成能力与水平基因转移（HGT）驱动的基因组可塑性。中东地区多重耐药（MDR）与碳青霉烯耐药负担重，IC2/ST2谱系广泛流行；但既有研究多限于单国、单院、短时间窗与小规模样本，缺乏将泛基因组、移动遗传元件（MGE）、纵向流行病学和机器学习整合的比较群体基因组框架，因此区域种群是否具有区别于全球株的基因组特征长期未明。研究人员以199个公开组装基因组为对象，回答中东株是否在克隆扩张背景下形成可识别区域基因组指纹，并评估机器学习用于区域溯源与监测的潜力。论文发表于《Microorganisms》。

主要关键技术方法上，研究人员从NCBI获取并质控人源A. baumannii组装基因组，形成中东队列100株（阿联酋、黎巴嫩、沙特、土耳其、伊朗、伊拉克、科威特）与全球队列99株的对比队列；用Prokka统一注释、QUAST评估组装质量，Panaroo严格模式重建泛基因组并做Jaccard距离与PCA，IQ-TREE 2建核心基因组最大似然树，Abricate比对CARD与VFDB做耐药组与毒力组分析，mlst判定序列型，Scoary做GWAS，Random Forest用泛基因组存在缺失矩阵分类中东与全球来源。

研究结果如下。3.1 泛基因组组成与附属基因组多样性：泛基因组10971个非冗余基因，核心1495个，云基因占61.0%；中东特有1102个、全球特有1970个，区域特有基因主要在附属与云部分，说明中东株附属层更收敛。3.2 耐药基因负担与移动组比较：中东株平均AMR基因30.09±4.53，全球37.63±9.67；中东移动相关基因均值20.42高于全球17.28，呈耐药负担较低而移动相关含量较高的反向关系。3.3 中东特有附属基因功能刻画：1102个特有基因中978个初始为假设蛋白；915个高质量序列经CD-Search后424个有保守结构域，多为代谢/酶功能，其次原噬菌、移动DNA与质粒稳定、毒素-抗毒素、应激与DNA修复、转录调控，显示HGT相关功能是主要可注释类别。3.4 基于附属基因组的种群结构：3354个附属/壳层基因聚类与PCA显示中东株成紧密块与致密簇，全球株碎片化、分散大，反映中东ST2主导下附属基因组保守。3.5 GWAS识别区域标记：火山图显示ISAba1变体、tufB等中东富集，ant1_2~aadA、pitA_1、thyA等全球富集，区域间患病率差达40%–50%。3.6 区域种群机器学习分类：Random Forest五折交叉验证准确率79.89%，Gini重要性最高为ISAba1变体与ampC β-内酰胺酶基因。3.7 移动遗传元件时间动态：2008–2016中东移动组密度稳定且高于全球，全球方差更大并有高移动离群基因组。3.8 毒力因子与序列型：毒力基因中位数约119，中东分布更窄；ST2在两组均约66%，说明ST2共有但中东株基因组更同质。3.9 耐药基因负担时间动态：中东2008–2016耐药基因数25–35为主，年际重叠，无显著积累趋势。3.10 耐药决定簇分布：bla_OXA-23两组均有，bla_OXA-69、bla_OXA-72、bla_ADC-191、bla_PER-7等更见中东代表，bla_OXA-24、bla_OXA-95、bla_ADC-222等更见全球代表，ade家族外排泵广泛存在。

讨论部分总结：研究人员指出中东种群不是因ST2比例更高而特殊，而是ST2/IC2背景下附属基因组保守、ISAba1富集与区域特有基因共同构成区域指纹。ISAba1可作插入序列也可上游启动子增强bla_OXA-23等表达，使区域株依赖既有耐药机制调控优化而非持续扩增AMR基因。较高MGE含量意味着基因组重塑、重排与未来选择压力响应能力，不一定等同更多耐药基因。毒力组总量两组接近，但中东生物膜与分泌系统相关基因均值更高，全球铁获取相关更高；acinetobactin的bas/bau基因在两队列均高流行。全球株反映多区域、多谱系、持续多样化，中东株反映瓶颈式克隆扩张与已适配医疗环境的稳定配置。该区域指纹可支撑分子监测、暴发调查与地理来源归因，但需更大独立队列验证；大量假设蛋白提示中东特有附属层仍有未表征适应潜能。

结论部分翻译：总之，该研究证明中东A. baumannii种群在全球盛行的ST2/IC2谱系中具有独特区域基因组指纹。尽管ST2在中东与全球队列均占优势，中东种群表现出基因组多样性降低、独特附属基因组轮廓，以及ISAba1插入序列显著富集。因此中东种群主要区别特征不是ST2流行率本身，而是其降低的基因组多样性、保守的附属基因组结构、ISAba1相关标记富集及区域特有附属基因。该区域基因组指纹还表现为抗菌药耐药决定簇、移动相关特征与毒力相关基因的差异。中东株维持较保守基因组配置，包括高丰度移动相关特征与保守毒力决定簇集合，同时总体AMR基因少于全球队列。ISAba1富集是区域基因组指纹的重要组分，可能促进相邻耐药决定簇调控，但本研究未评估功能表达。同样，毒力相关基因保守并不证明表型适合度权衡不存在。总体而言，这些发现表明中东A. baumannii种群在ST2/IC2主导种群中具有独特且相对保守的基因组配置；持续基因组监测联合功能与流行病学研究十分重要。