《Microbiology Spectrum》：A global genomic analysis of Streptococcus equi subsp. zooepidemicus reveals clade-specific pathogenic and antimicrobial resistance profiles with several mobile genetic elements as key adaptive drivers
编辑推荐：
专业术语翻译（摘要部分）：
马链球菌兽疫亚种（Streptococcus equi subsp. zooepidemicus，S. zooepidemicus）是一种影响多种哺乳动物的人畜共患病原体，然而其基因组流行病学仍鲜有了解。为应对这一问题，研究人员分
专业术语翻译（摘要部分）：
马链球菌兽疫亚种（Streptococcus equi subsp. zooepidemicus，S. zooepidemicus）是一种影响多种哺乳动物的人畜共患病原体，然而其基因组流行病学仍鲜有了解。为应对这一问题，研究人员分析了来自6个国家（1993—2024年）的428株S. zooepidemicus菌株。全基因组单核苷酸多态性（single-nucleotide polymorphism，SNP）分析揭示了四个分支（clade），遗传多样性主要由谱系驱动，而非宿主物种或地理来源。研究结果支持S. zooepidemicus为宿主泛嗜性（host generalist）病原体的结论，其中马为主要宿主。所鉴定的毒力相关基因（virulence-associated genes，VAGs）呈现分支特异性分布模式：分支1缺乏黏附因子基因srtC.2/3/4，分支2富集IgG降解酶基因ideZ，分支4富集毒力相关偶联蛋白基因virD4的特定等位基因。抗菌药物耐药基因（antimicrobial resistance genes，ARGs）在18.2%的菌株中检出，未发现多重耐药菌株。ARG流行率与分支4相关，其中四环素耐药最为常见。移动遗传元件（mobile genetic elements，MGEs）在塑造辅助基因组（accessory genome）图谱中起关键作用：前噬菌体（prophage）是外毒素基因的主要携带者，而整合性接合元件（integrative and conjugative elements，ICEs）则是ARGs的主要载体。值得注意的是，侵袭相关的virD4等位基因与ARGs在ICESz1上的共定位主要见于分支4分离株。总之，S. zooepidemicus的致病潜力与耐药谱由MGE驱动的特定进化谱系所塑造。这些元件不仅作为遗传货物，还作为提供接合机制以共传播毒力与耐药性状的驱动力，加速人畜共患谱系的进化。
论文解读文章：
马链球菌兽疫亚种（Streptococcus equi subsp. zooepidemicus，S. zooepidemicus）是一种革兰氏阳性β-溶血性条件致病菌，最常见于马的上呼吸道和下生殖道，也可引起多种哺乳动物包括人的败血症、脑膜炎、心内膜炎和关节炎，具有明确的人畜共患潜力。尽管该菌在兽医学和公共卫生领域具有重要意义，其基因组流行病学长期缺乏系统认识。既往基因组研究多聚焦于特定暴发事件，对该物种的群体结构、进化驱动力和辅助基因组组成缺乏全局性理解。同时，虽然青霉素仍为治疗疑似S. zooepidemicus感染的一线药物，但抗菌药物耐药性（antimicrobial resistance，AMR）的报道日益增多，而相关基因组数据仍然有限。因此，亟需通过整合基因组分析阐明其种群结构、进化机制、毒力相关基因（virulence-associated genes，VAGs）与抗菌药物耐药基因（antimicrobial resistance genes，ARGs）的分布规律，为防控策略提供依据。
研究人员对从美国国家生物技术信息中心（NCBI）序列读长档案（Sequence Read Archive，SRA）下载的910株S. zooepidemicus全基因组测序原始数据进行严格质控和从头组装，最终保留428个高质量基因组用于后续分析。这些菌株分离自1993年至2024年，主要来自美国（n=293）和冰岛（n=106），另有少量来自英国、荷兰、瑞典和爱尔兰；宿主以马为主（n=389），另有犬、人、猫、牛、豚鼠、猪和羊等。核心基因组系统发育分析表明，S. zooepidemicus可分为四个系统发育分支（clade 1–4）。泛基因组分析显示其泛基因组呈开放型，共鉴定出8064个基因簇，其中核心基因1385个（17.2%），软核心基因56个（0.7%），壳基因748个（9.3%），云基因5875个（72.9%）。云基因显著富集于“复制、重组与修复”“防御机制”和“一般功能预测”等COG功能类别，提示频繁的水平基因转移（horizontal gene transfer，HGT）事件。重组分析显示重组与突变多态性比率（r/m）约为7.49，表明重组是该菌谱系多样化的重要进化驱动力。基于全基因组核心SNP的系统发育树显示四个分支，分支内平均SNP差异为15763至24292，分支间差异为24833至25974，且菌株未按宿主或来源国形成单一单系群。辅助基因组主成分降维分析（UMAP）进一步显示，辅助基因组的聚类结构主要与系统发育分支相关，而宿主物种和地理来源的解释力极弱；PERMANOVA分析确认分支是辅助基因组变异的主要驱动因素（R2=0.148，F=24.68），宿主物种（R2=0.031）和地理来源（R2=0.082）解释力可忽略。
研究人员利用自制链球菌毒力因子数据库、ResFinder和CARD-RGI等工具对VAGs和ARGs进行注释，并利用PHASTEST、ICEscreen、ICEfinder2、IslandViewer4、PlasmidFinder和MobSuite等多种生物信息学工具鉴定前噬菌体、整合性接合元件（integrative and conjugative elements，ICEs）、基因组岛（genomic islands，GIs）、整合性可移动元件（integrative mobilizable elements，IMEs）、质粒和接合模块等移动遗传元件（mobile genetic elements，MGEs），并分析其携带的VAGs和ARGs。
结果部分：一、S. zooepidemicus的群体结构。428株高质量基因组平均基因组大小为2.0–2.2 Mb，GC含量为40.92%–41.87%。泛基因组拟合幂律函数得到指数γ=0.233，表明泛基因组开放。核心SNP系统发育树将菌株划分为四个分支，菌株来源宿主和国家彼此混杂，支持S. zooepidemicus为宿主泛嗜性病原体。UMAP降维显示辅助基因组按分支分离，而PERMANOVA和平均轮廓宽度（ASW）进一步确认分支为辅助基因组变异的主要来源。二、VAG分布模式。共鉴定出66个推定的VAGs，其中27个为核心VAGs，39个为辅助VAGs。辅助VAGs呈现分支特异性分布：分支1缺乏黏附因子基因srtC.2/3/4，而sclC在分支1（70.6%）和分支4（51.0%）存在但在分支2完全缺失；分支2富集IgG降解酶基因ideZ（61.7%）；分支4则富集侵袭相关偶联蛋白基因virD4的特定等位基因（87.9%）以及srtC.2/3/4。外毒素基因szeN和szeP主要在分支1富集，szeF主要在分支2（76.5%）。vapE仅存在于分支4的12株菌中。免疫调节基因se18.9和zag在分支3和4中较高而分支1和2中缺失或极低。多数VAGs在不同宿主或地理来源间无显著差异，但ideZ在欧洲菌株中更常见，srtC.3/4在北美菌株中更常见。三、ARG分布模式。ARGs在18.2%（78/428）菌株中检出，主要来自马（70株）。共鉴定出8种ARGs，赋予对四环素类（tet(M)、tet(O)、tet(W)）、林可酰胺类（lsa(C)、lnu(C)）、二氨基嘧啶类（dfrF）和氨基糖苷类（aac(6')-Im、aph(2")-IIa）的耐药性。四环素耐药基因最常见（11.4%）。59.0%的ARG阳性菌株属于分支4，且77/78株来自美国。仅一株分支1菌株同时携带aac(6')-Im和aph(2")-IIa。四、MGEs作为VAG和ARG传播的驱动力。共鉴定出33种不同完整前噬菌体，出现在260株菌中，主要携带外毒素基因，如磷脂酶基因slaA；前噬菌体在ARG传播中作用甚微，仅噬菌体Phi_Geobac_E2携带tet(O)、lsa(C)和dfrF。共检测出18种ICESz1型ICE，其中ICESz1最常见（35%），主要存在于分支4，是virD4等位基因的主要携带者；ICESthJIM8232-1（31%）主要与分支2相关。ICESz1携带大多数量ARGs，包括tet(M)、tet(O)、tet(W)、lnu(C)、lsa(C)和dfrF。ICE携带了65.7%的所有ARGs。GIs在388株菌中检出，主要分布在分支3和4，也主要携带virD4等位基因和四环素耐药基因tet(W)。其他元件如IMEs、质粒和接合模块的作用有限，质粒上未检出VAGs或ARGs。五、ICESz1介导的VAG和ARG共定位。在46株菌中，VAGs和ARGs位于同一MGE上，其中89%（41/46）为ICESz1。最常见的共定位组合为侵袭相关virD4等位基因与四环素耐药基因（tet(W)、tet(O)）或林可酰胺耐药基因lsa(C)。
讨论部分总结：S. zooepidemicus具有开放的泛基因组和大量辅助基因，与克隆性且宿主限制的S. equi形成鲜明对比，这种差异可能与其栖息更多样化的黏膜环境有关，促进了与其他细菌的频繁互作和HGT。云基因中“防御机制”和“复制、重组与修复”类别的显著富集进一步支持频繁的HGT事件。高r/m比率表明重组在谱系多样化中发挥重要作用。VAGs在核心与辅助基因组中的分布差异提示S. zooepidemicus维持稳定的核心毒力机制以保障基本致病性和环境持久性。不同分支间VAG谱的差异表明谱系与辅助基因组协同进化。superantigen基因szeF、szeN、szeP的分支特异性分布可能提示分支相关的疾病表型遗传基础，但受限于缺乏临床元数据。ideZ在分支2中的富集可能赋予宿主适应优势。vapE仅存在于分支4，挑战了其仅与猪宿主毒力相关的既往观点，更可能为分支4的遗传标记。黏附因子基因的相反分布提示不同分支可能进化出不同的宿主定植策略。ARG流行率（18.2%）及无多重耐药的结果与部分表型耐药报道存在差异，可能源于基因型预测与表型表达之间的差异，以及菌株收集和检测方法的不同。ARG阳性菌株高度集中于美国，提示地区性抗菌药物选择性压力差异。青霉素耐药基因未检出，支持青霉素仍可作为一线治疗药物。四环素耐药基因最为常见，与既往报道一致。lsa(C)检出于45株菌。aac(6')-Im和aph(2")-IIa在单一IMEs上的共定位及其与肠杆菌科质粒同源基因的高相似性（>96%）提示氨基糖苷类耐药基因可能从革兰氏阴性菌水平转移至S. zooepidemicus。MGEs分析显示，前噬菌体是外毒素基因的主要载体，且大多前噬菌体GC含量低于宿主基因组，提示为较近期获得；而ICESz1是ARGs的主要载体，其携带的侵袭相关virD4等位基因与ARGs共定位，可能提供致病优势及接合传播机制。ICESz1 GC含量变异最大，提示该元件处于活跃的基因获得与丢失进化过程。GIs在多数菌株中检出，主要分布于分支3和4，并在分支4菌株中高频携带virD4等位基因和tet(W)。质粒在S. zooepidemicus中的作用似乎可忽略。
研究结论：S. zooepidemicus群体由四个系统发育分支构成，每个分支具有独特的辅助基因库。存在分支特异性VAGs，但其是否影响细菌毒力仍有待阐明。人源菌株与其他菌株混杂，凸显其人畜共患风险和采用“一体化健康”（One Health）视角的重要性。MGEs，主要是前噬菌体和ICESz1，对S. zooepidemicus基因组可塑性和分支特异性辅助基因组有重要贡献。前噬菌体参与外毒素VAGs的传播，而ICESz1和GIs主要携带侵袭相关且具有接合转移功能的virD4等位基因和四环素耐药基因，且具有明确的分支关联。ARGs与VAGs在ICESz1上的共定位表明，抗菌药物使用不仅可能选择毒力，还可能促进其水平传播。该研究为S. zooepidemicus提供了新的基因组图谱，强调MGEs是其进化的毒力与耐药谱的关键驱动力。
