《Journal of Clinical Microbiology》：Genomic epidemiology of Enterococcus faecium bloodstream infections during a VanB-type VRE peak reveals an oligoclonal scenario: an observational study at a German university hospital (2017–2022)
编辑推荐：
本研究针对德国一所大学医院在2017-2022年间发生的VanB型耐万古霉素屎肠球菌（VRE）感染先激增后锐减的异常流行病学现象，通过整合高分辨率基因组学与患者时空移动数据，深入解析了血流感染分离株的克隆结构、传播动态及van基因型分布，为区分真正暴发与多克隆聚集、优化院感防控策略提供了精准的分子流行病学证据。
进入21世纪，耐药细菌的传播已成为全球公共卫生的重大挑战。其中，耐万古霉素屎肠球菌（Vancomycin-resistant Enterococcus faecium, VRE）作为一种臭名昭著的医院获得性病原体，因其对“最后防线”抗生素的耐药性和在医疗机构中的持续传播能力，备受关注。在德国，VRE的流行病学图谱并非一成不变：2010年代末，一种特定类型——携带vanB耐药基因的VRE（VanB型VRE）——的感染病例经历了戏剧性的“过山车”式变化，先是急剧攀升至顶峰，随后又快速回落。这个神秘的、短暂的高峰背后，究竟隐藏着怎样的克隆传播秘密？是单一优势克隆的集中暴发，还是多种菌株的“群雄并起”？传统的感染控制监测方法，能否准确捕捉到这些微观层面的种群动态变化？为了解开这些谜团，一支研究团队深入挖掘了一家德国大学医院长达五年的临床数据宝库。
这项研究旨在深入剖析2017年至2022年间，一所德国大学医院内与血流感染相关的屎肠球菌种群结构及其动态变化，特别是聚焦于那场来去匆匆的VanB型VRE感染高峰。研究人员采用了一种整合的基因组流行病学方法，对在此期间收集的344株血流感染分离株进行了回溯性分析。他们不仅进行了传统的分子分型，如多位点序列分型（MLST）和核心基因组多位点序列分型（core genome MLST, cgMLST），还应用了分辨率更高的split k-mer分析（SKA），以绘制出精细的菌株遗传图谱。更重要的是，他们将这份基因图谱与通过机器学习辅助分析的、详细的纵向患者移动数据相结合，试图在基因相似性的基础上，重建可能的患者间传播链条，从而区分真正的院内传播事件与仅仅是表型巧合的基因多样化聚集。相关成果已发表在微生物学领域的权威期刊《Journal of Clinical Microbiology》上。
为开展此项研究，作者运用了多项关键技术。首先，他们从德国汉诺威医学院2017年10月至2022年2月间的血流感染病例中，系统收集了屎肠球菌分离株作为研究队列。其次，对所有菌株进行了全基因组测序，并利用生物信息学工具进行质量控制和基因组组装。在分子分型方面，核心工作包括：1）基于7个看家基因的多位点序列分型（MLST）和基于1,423个核心基因位点的cgMLST分型，用于界定菌株的序列型（ST）和复合型（CT）；2）高分辨率的split k-mer分析（SKA），用于在k-mer水平进行更精细的聚类。此外，研究还通过生物信息学方法筛查了vanA、vanB等耐药基因以及一系列毒力基因。最后，通过编程自动查询医院信息系统，提取并分析了患者的住院和病房转移数据，以寻找可能的流行病学联系。
研究结果
The rise and fall of VanB-type VRE at Hannover Medical School
数据分析显示，在汉诺威医学院，VanB型VRE病例数在2017年开始急剧上升，于2019年达到顶峰，随后稳步下降，到2022年已回落至2017年之前的水平。而VanA型VRE的病例数在整个研究期间保持相对稳定。这表明观察到的VRE流行趋势变化主要是由VanB型菌株驱动的。
Isolate collection strategy
研究最终纳入了333株经过质量控制的血流感染屎肠球菌分离株进行分析。其中，238株不携带van基因（van-阴性），76株为vanB阳性，19株为vanA阳性。菌株数量在年份间波动，2021年收集到的菌株最多。
Temporal shifts in clonal lineages and van-genotype prevalence
种群分析揭示了克隆谱系和van基因型随时间的动态变化。vanB阳性分离株的比例在2019年下半年后从30-40%急剧下降至约10%。vanB感染高峰主要由ST117/CT71、ST117/CT36和ST117/CT1917等克隆谱系驱动。与此同时，van-阴性菌株（主要是ST117/CT929和ST117/CT2505）逐渐成为优势种群。vanA阳性分离株则从2019年才开始出现，且数量较少，主要与ST80/CT1470相关。
Prevalence of van genotypes and virulence genes across diverse cgMLST clusters
基于cgMLST的最小生成树分析将菌株分为17个簇和50个单例。四个主要的cgMLST簇（簇1-4）与特定的ST/CT及van基因型相关：簇1（ST117/CT929和CT2505）主要为van-阴性；簇2（ST117/CT71）主要为vanB阳性；簇4（ST80/CT1470）中含有较高比例的vanA阳性菌株。毒力基因分析显示，医院相关变异体在所有主要克隆谱系中均占主导地位，且超过90%的分离株携带了可能与定植优势相关的细菌素基因T8。
SKA provides higher resolution clustering than cgMLST
研究比较了不同分型方法的分辨率。MLST将菌株分为16个ST型，cgMLST进一步细分为70个CT型，并以此定义了20个可能的暴发相关簇。而基于k-mer的SKA方法由于涵盖了附属基因组和基因间区，提供了更高的分辨率，将cgMLST定义的簇进一步分解为40个更小的簇和144个单例，显示出更强的区分能力。
Longitudinal analysis of E. faecium ST117/CT71 transmission
研究深入分析了作为vanB高峰主要驱动力的cgMLST簇2（ST117/CT71）。SKA将该簇进一步分解为8个子簇和31个单例。流行病学数据分析显示，在SKA定义的8个子簇中，有5个存在明确的流行病学联系（如同病房30天内接触）。值得注意的是，即使在SKA定义为单例的菌株中，通过分析患者移动数据，也发现了大量（涉及29名患者）跨越长时间的潜在流行病学联系，这些联系的中位时间间隔长达106天。这表明，仅凭严格的基因分型阈值可能会遗漏由长期、低频传播事件构成的传播链。
研究结论与意义
本研究通过整合高分辨率基因组学和详细的流行病学数据，清晰描绘了德国一所大学医院内屎肠球菌血流感染种群在VanB型VRE高峰期间的复杂动态。核心结论是，这场高峰并非由单一克隆的暴发引起，而是一个涉及多个流行谱系（如ST117/CT71、ST117/CT36、ST117/CT1917、ST80、ST262等）的“寡克隆”场景。不同的谱系具有独特的van基因型，且随时间此消彼长。研究证实，van阳性与van阴性种群之间的基因重叠极小，表明耐药基因的水平转移在此次流行中并非主要驱动力，克隆扩张是主导机制。
这项研究的重要意义在于方法论和实践应用两个层面。在方法学上，它证明了对于解析复杂的、长期的医院内传播动态，单一的分型方法或固定阈值存在局限。cgMLST（使用≤15个等位基因差异的阈值）具有较好的敏感性，能捕获可能的传播簇，但可能高估关联性；而高分辨率的SKA（使用7个SNP的阈值）特异性更强，能区分紧密相关的克隆，但可能因过于严格而割裂实际的传播链。因此，采用cgMLST进行初步筛查，再结合SKA进行高分辨率验证，并整合患者时空移动数据进行分析的“阶梯式”整合策略，能更精准地区分真正的单克隆暴发与寡克隆传播聚集。
在实践意义上，该研究为医院的感染预防与控制团队提供了关键见解。它表明，即使在没有典型“暴发”信号的情况下，院内也可能存在由优势克隆导致的持续性低水平传播。准确识别这些传播模式，有助于将有限的感染控制资源精准投入到经基因和流行病学证据支持的真实传播事件上，从而避免因误判多克隆聚集为暴发而实施的、可能造成不必要干扰的广泛措施。此外，研究发现即使是不携带van基因的敏感菌株（如ST117/CT929、ST117/CT2505），也可能形成医院适应的克隆并传播，提示对这些“敏感”但可能具有高传播性的菌株，也应保持警惕并采取相应的感染控制措施。总之，这项研究强调了持续、精准的分子流行病学监测与多维度数据分析相结合，对于理解和控制VRE等多重耐药菌在医院环境中的传播具有不可替代的价值。
