脑膜炎败血伊丽莎白菌(Elizabethkingia miricola)的传播与基因组学特征：一种具内在耐药性和院感暴发潜力的动物源性病原体

摘要：脑膜炎败血伊丽莎白菌(Elizabethkingia miricola，E. miricola)是一种新现的机会性病原体，临床分离率日益升高，其感染带来重大治疗挑战。研究人员对5株临床E. miricola分离株进行全基因组测序，并与全球来源菌株进行比较基因组学分析，以探究其耐药机制、致病性、系统发育关系及传播模式。系统发育分析显示E. miricola的进化与样本来源（临床/环境/动物）、地理起源或地区无严格关联；两栖动物来源分离株聚类于临床谱系内，提示环境储源（可能为水槽积水）是感染来源。地理分布上，致病性E. miricola呈全球播散，在中国和欧洲流行率较高。推定传播网络提示院内播散及跨国循环，存在院感暴发风险。人与蛙来源菌株共享耐药基因（含BlaB和GOB β-内酰胺酶(β-lactamase)）及毒力谱，凸显其内在耐药与致病潜能。米诺环素(minocycline)和多西环素(doxycycline)体外抗菌活性强，是潜在治疗选择。上述发现提示E. miricola可通过院内及潜在国际水平传播构成人畜共患威胁，预防暴发需加强该病原体的基因组监测。

论文解读：

《One Health》发表的该研究针对新现机会性革兰氏阴性杆菌——脑膜炎败血伊丽莎白菌(Elizabethkingia miricola，E. miricola)，针对其临床分离罕见、基因组数据匮乏、进化史与传播路径不明、动物源性潜力未定及广泛耐药导致治疗困难等现状，研究人员对浙江大学医学院附属第一医院2013–2018年临床分离保存的5株E. miricola进行全基因组测序（Whole Genome Sequencing, WGS），并结合NCBI数据库中21株全球来源（临床、环境、两栖动物来源）E. miricola全基因组序列开展比较基因组学分析，阐明其核心基因组与泛基因组特征、系统发育关系与地理分布、潜在传播网络、耐药基因型与表型关联及毒力因子(virulence factor, VF)组成，证实E. miricola具内在多重耐药与固有致病潜能，两栖动物来源株与临床株亲缘关系密切提示人畜共患风险，存在院内及跨国传播可能，米诺环素与多西环素为其敏感治疗药物，为临床防控与新发病原监测提供依据。

主要关键技术方法：研究纳入2013–2018年来自浙江大学医学院附属第一医院菌株样本库经基质辅助激光解吸电离飞行时间质谱（Matrix-Assisted Laser Desorption/Ionization Time-of-Flight Mass Spectrometry, MALDI-TOF MS）初筛并平均核苷酸一致性（Average Nucleotide Identity, ANI＞95%）确认种的5株临床E. miricola（痰3株、血1株、尿1株），联合NCBI下载的21株全球来源公开全基因组序列（临床/环境/两栖动物来源）；采用Illumina NovaSeq平台进行全基因组测序与de novo组装，Glimmer与GeneMarkS预测编码序列(coding sequence, CDS)、tRNAscan-SE与rRNAmmer预测非编码RNA并作KEGG功能注释；FastANI计算ANI评估种间相似性，CMG-Biotools做核心/泛基因组(pan-genome/core-genome)分析，RedDog流程比对与单核苷酸多态性(Single Nucleotide Polymorphism, SNP)过滤后用PhyML构建最大似然(Maximum Likelihood, ML)系统发育树并以随机中点生根法(rooting)；最小生成树(Minimum Spanning Tree, MST)估算遗传距离，≤10个SNP差异归为同一推定传播群；综合抗生素耐药数据库(Comprehensive Antibiotic Resistance Database, CARD)用Resistance Gene Identifier(RGI)预测获得性及固有耐药基因，毒力因子数据库(Virulence Factors of Pathogenic Bacteria Database, VFDB)注释毒力因子；药敏试验采用琼脂稀释与肉汤稀释法按临床和实验室标准协会(Clinical and Laboratory Standards Institute, CLSI)及美国食品药品监督管理局(Food and Drug Administration, FDA)（四环素类新药）折点判读，大肠杆菌ATCC 25922等作质控。

研究结果：

3.1 E. miricola物种基因库组成（Gene repertoire of E. miricola species）：通过对纳入菌株渐进式分析获保守核心基因2607个，泛基因组呈开放性并持续获得附属基因(accessory genes)；菌株SKLX011872含415个特有基因，其余4株临床分离株仅含1～5个株特异性基因，提示多数临床谱系遗传同质性强。

3.2 系统发育推断与地理分布（Phylogenetic inference and geographical distribution）：基于全基因组SNP构建ML系统发育树可将26株分为4个主要谱系(lineage)，各含亚簇，临床与环境株大体分叉但部分两栖动物及环境株嵌入临床分支，进化与样本来源、地理或区域无显著关联；亚洲（尤中国）基因型多样性最高，菌株呈广域分布（亚、欧、北美），提示跨境传播能力。

3.3 遗传聚类与传播浮现（Genetic clusters and transmission emergence）：以≤10个SNP界定假定传播群，发现同一医院不同患者分离株（SKLX069005与SKLX070046）基因组高度一致提示院内同源传播；德产蛙源三株聚为一簇提示动物内传播；美、德、日三株跨国株SNP差异极小提示共同起源与潜在国际播散，表明E. miricola具局部院感暴发与跨国循环能力。

3.4 抗菌药物敏感性及耐药决定因子（Antimicrobial susceptibility and resistance determinants）：5株临床分离株对青霉素类、头孢菌素类、碳青霉烯类、氨基糖苷类、多肽类、氯霉素类、多数四环素类及磷霉素共≥16种抗生素高水平耐药，符合广泛耐药(extensively drug-resistant, XDR)表型；左氧氟沙星(levofloxacin)、米诺环素(minocycline)、多西环素(doxycycline)及替加环素(tigecycline)具良好体外抗菌活性。所有基因组均携带固有碳青霉烯酶基因Bla_B与GOB β-内酰胺酶、头孢菌素耐药基因CME、氨基糖苷类耐药基因RanA/RanB/aadS及氟喹诺酮类耐药相关基因sdiA，证实属水平内在(intrinsic)多重耐药；耐药基因分布与来源及地域无关。

3.5 毒力基因与潜在致病性（Virulence genes and potential pathogenicity）：检出至多38类毒力相关基因，包括荚膜多糖合成、鞭毛/菌毛、脂多糖(lipopolysaccharide, LPS)/脂质A生物合成、去铁胺(desferrioxamine)转运、血红素(heme)合成(hemL普遍存)、抗氧化应激（过氧化氢酶 catalase、过氧化物酶 peroxidase、热休克蛋白 heat shock protein）、双组分调控系统及粘附、分泌系统、溶血素等；22个毒力基因为所有菌株共有，人与两栖动物来源株毒力谱无来源或地域特异差异，支持E. miricola具固有致病潜能。

讨论与结论翻译：E. miricola广泛存在于自然环境与医疗机构，为新现威胁性机会致病菌，可致高危人群侵袭性感染甚至ICU聚集暴发与高病死率，常被误诊或漏检为污染，16S rRNA基因难以区分伊丽莎白菌属(Elizabethkingia)种，建议rpoB与mutT基因qPCR、重组酶聚合酶扩增结合侧流层析(recombinase polymerase amplification combined with lateral flow dipstick, RPA-LFD)、荧光探针RPA或基质辅助激光解吸电离傅里叶变换红外光谱(Matrix-Assisted Laser Desorption/Ionization Fourier-Transform Infrared Spectroscopy, FTIR)及宏基因组二代测序(metagenomic next-generation sequencing, mNGS)提高检出。本菌具内在广泛耐药，碳青霉烯类等通常无效，米诺环素与多西环素体外敏感可作治疗参考，氟喹诺酮类需警惕耐药，中药提取物亦具研究前景。局限性为全球性样本代表性待提升及耐药机制需基因编辑验证。结论：E. miricola是可致人类危及生命感染的重要病原菌，不同来源株共享耐药基因与毒力因子分布证实其内在多重耐药与致病性；研究初步表明E. miricola具人畜共患潜力，能经院内及跨国途径传播，存在全球扩散威胁，结果为制定针对性防控策略提供基础。