沙门氏属蒙特维迪亚血清型（Salmonella enterica serovar Montevideo，S. Montevideo）与美国牛及牛肉产品相关的毒力减弱现象——基于基因组与流行病学证据

《Frontiers in Microbiology》：Genomic and epidemiologic evidence for attenuated virulence in Salmonella Montevideo associated with cattle and beef products in the United States

　　

摘要：肠炎沙门氏菌（Salmonella enterica）蒙特维迪亚血清型（S. Montevideo）常从牛、牛源产品及相关环境中分离，但通过该类来源引起人类疾病较为罕见，提示其不同谱系间存在对人毒力及生态位适应的差异。基于研究人员此前鉴定出的包含12个具不同生态关联的进化支（clade）、占主导的S. Montevideo系统发育群，研究人员开展了比较基因组与进化分析，以表征在人源关联及来源适应性上存在差异的进化支。研究人员鉴定出一个非人源关联（non-human-associated, NHA）进化支（clade 10，牛适应型）和三个与人源关联（human-associated, HA）进化支（clade 7，非动物环境适应型；clades 3和6，无明确来源特异性），四者均表现为开放泛基因组（open pan-genome），表明存在持续基因交流。相对于NHA进化支，HA进化支中显著富集（overrepresented）的基因35%–66%位于可移动遗传元件（mobile genetic elements, MGEs）上，而HA进化支中富集的所有5个毒力因子均由染色体编码。与HA进化支相比，NHA进化支在编码毒力因子的基因（涉及上皮细胞入侵及全身性感染）中提前终止密码子（premature stop codons, PMSCs）/无义突变富集达7–19倍，提示功能缺失与宿主适应及毒力降低相关。此外，进化重建显示clade 10起源最近且进化速率最快，反映快速宿主适应；clade 7为古老进化支，具基因组稳定性，提示长期环境持久定殖。Clades 3和6呈中等进化速率，其中clade 3包含两个近期出现的亚支，具克隆扩增特征。本研究强调需针对农业动物及其源性食品中的沙门氏菌采取基于风险的防控策略，以提高公共卫生效益。

论文解读：《Genomic and epidemiologic evidence for attenuated virulence in Salmonella Montevideo associated with cattle and beef products in the United States》发表于《Frontiers in Microbiology》

一、研究背景与意义

肠炎沙门氏菌（Salmonella enterica）蒙特维迪亚血清型（S. Montevideo）是美国牛、牛肉及牛相关环境中高频检出的血清型，但由牛肉和乳制品引发的S. Montevideo人类感染暴发极为少见。这种现象提示S. Montevideo内部存在谱系（cline/clade）水平的异质性——部分谱系可能已发生宿主适应性演化并伴随对人毒力的减弱（attenuated virulence）。此前研究人员（Chen et al., 2024）发现S. Montevideo优势系统发育谱系（S. Montevideo A）含12个进化支（clade），各clade来源分布及与人源分离株的关联程度（odds ratio, OR）差异显著，其中clade 10以牛源为主且人源极少，clade 7则富集于人源及非动物自然环境。本研究在此基础上，通过比较基因组学与分子进化分析，阐明与人关联（HA, human-associated）及非人关联（NHA, non-human-associated）clade之间毒力差异的遗传基础、宿主适应机制及进化轨迹，为沙门氏菌按亚型风险分级管控提供科学依据。

二、主要关键技术方法概述

研究人员选取NCBI病原体检测（Pathogen Detection, PD）数据库中S. Montevideo A的全基因组测序（WGS）数据及元数据，限定美英（US/UK）分离株计算各clade人源分离株比值比（odds ratio, OR）并进行Fisher精确检验与Benjamini–Hochberg校正，按OR>2且p<0.05定义为HA、OR<0.5且p<0.05定义为NHA、其余为中性（Neutral），据此确定HA clades（3、6、7）与NHA clade（clade 10）。从各目标clade按SNP聚类分层随机抽样40株进行泛基因组（pan-genome）分析（Prokka注释、Panaroo推断、Heaps定律评估开放性），通过pan-GWAS（Scoary）筛选HA/NHA富集的附属基因（accessory genes）并与毒力因子数据库（VFDB、PATRIC_VF、Victors）比对；用kSNP4进行无参核心SNP（core SNP）检出并识别提前终止密码子（premature stop codon, PMSC），GWAS分析PMSC在各clade的富集情况；通过Platon和PHASTER判定基因位于染色体/质粒/前噬菌体（prophage）；按采样年份分层随机选100株各clade代表性分离株进行参考基因组依赖SNP检出、重组剔除（Gubbins）后，用BEAST v2.5.2宽松分子钟与贝叶斯天际线（Bayesian Skyline Plot, BSP）模型估算分歧时间（TMRCA）及进化速率，并通过Date-Randomization Test验证时间信号。

三、研究结果

The major phylogenetic lineage of S. Montevideo comprises clades with significant over- or underrepresentation of human isolates（S. Montevideo主系统发育谱系包含人源分离株显著过代表征或低代表征的进化支）

对S. Montevideo A 12个clade进行流行病学分型，确认clades 3、6、7为HA（OR>2, p<0.05），clade 10为NHA（OR<0.5, p<0.05）。系统发育树显示clades 3与6亲缘关系近共享最近共同祖先，clade 7与clade 10分别适应非动物自然环境与牛源，支持S. Montevideo A各clade已发生导致来源特异性适应及人源致病潜能差异的分化。

Source-specific adaptation does not constrain pan-genome expansion among HA and NHA clades within S. Montevideo A（来源特异性适应不限制S. Montevideo A内HA与NHA进化支的泛基因组扩展）

四个目标clade泛基因组大小介于5206（clade 7）至5542（clade 3）基因，Heaps定律拟合α值均<1（clade 7 α=0.944最高，clade 3 α=0.923最低），表明所有clade均为开放泛基因组（open pan-genome），提示来源适应未限制基因获得/丢失动态，各clade在其生境中仍持续进行水平基因转移（horizontal gene transfer, HGT）。

Comparative genomics reveals clade-specific gene content variation and diverse contributions of mobile genetic elements (MGEs)（比较基因组学揭示clade特异性基因含量变异及可移动遗传元件的不同贡献）

相对于NHA clade 10，HA clades 3、6、7分别显著富集174、198、283个HA基因（HA genes），clade 10相对各HA clade分别富集160、123、209个NHA基因。五个定位于染色体且跨HA clade一致或选择性富集的已知毒力因子为STM0520（编码假定硫代乙酸转运蛋白SauU）、hyi（STM0518，羟丙酮酸异构酶）、allD（STM0528，脲基乙醛酸脱氢酶）——三者于全部三个HA clade富集；STM0030（LysR家族转录调节因子LeuO）——于clades 3、6富集；steB（STM1629，III型分泌系统效应蛋白Type III secretion system effector, T3SS effector）——仅于clade 7富集。HA基因中35%–66%位于MGE（质粒或前噬菌体），但上述毒力因子均不在MGE上；clade 3部分HA基因来自IncI1质粒（含IncI1.1复制子），clades 6和7 HA基因以前噬菌体/染色体来源为主。

NHA-clade-specific accumulation of premature stop codons indicates functional genome disruption linked to cattle adaptation and attenuated human virulence（NHA clade特异性提前终止密码子积累提示功能基因组破坏与牛适应及人毒力减弱相关联）

clade 10较各HA clade核心基因组中出现显著更多PMSC（clade 10: 170–204个 vs. HA: 17–29个），其中影响毒力因子的PMSC富集7–19倍。clade 10中13个被破坏毒力相关基因显著富集于该clade，包括：参与牛致死感染的sspA（严谨饥饿蛋白A Stringent starvation protein A）；参与牛肠道定植的tdcE、mglA、acnA、nrdD；以及编码SPI-1（Salmonella pathogenicity island 1）III型分泌系统（Type III Secretion System, T3SS-1）组分的invE（入侵蛋白InvE gatekeeper）、sicP（SptP分子伴侣chaperone）、orgB（侵袭相关蛋白OrgB）。这些PMSC导致关键入侵与全身感染相关基因截短失活，提示clade 10向牛共生/定植方向演化并伴随对人侵染能力减退。

Divergent evolutionary trajectories of HA and NHA clades reflect the accelerated evolution and cattle adaptation of clade 10（HA与NHA进化支的分歧进化轨迹反映clade 10加速演化及牛适应）

分子钟分析显示：clade 10约1905年最近共同祖先（TMRCA≈116年），进化速率最高（1.9×10^?7substitutions/site/year），符合牛肠道内短世代时间及放松净化选择（relaxed purifying selection）下适应假基因化（pseudogenization）；clade 7最古老（TMRCA≈3765年，约1744 BCE），进化速率最低（8.1×10^?9substitutions/site/year），体现非动物自然环境中长期基因组稳定与环境持久性；clades 3（≈1837年，8.9×10^?8）与6（≈1852年，1.1×10^?7）呈中等速率且无明确来源偏好，clade 3含两近期亚支（约2007、2013年）具克隆扩增及有效群体变小迹象，clade 6具地理/序列型（ST195欧洲中东 vs. ST4美国，部分具抗微生物药物抗性AMR）分化。

四、讨论与结论总结（翻译浓缩）

讨论指出S. Montevideo A内clade 10为牛适应、毒力减弱的共生谱系体，其SPI-1 T3SS组分（invE、sicP、orgB）及牛致病/定植相关基因（sspA、tdcE、mglA、acnA、nrdD）中PMSC积累致功能缺失，支持宿主适应伴随毒力衰退假说；HA clades 3、6、7与clade 10的差异毒力相关基因均为染色体编码，MGE主要驱动clade间附属基因组多样化但不贡献已知毒力因子差异；clade 7低进化速率反映自然环境长期定殖，具水源/农产品污染潜在公共卫生风险；clade 10高进化速率源于牛宿主内放松净化选择与适应性基因丢失；clade 3近期亚克隆扩增虽未关联暴发但提示潜在公共卫生威胁上升可能。研究强调S. Montevideo种内谱系毒力不均一，应实施基于谱系风险分级（risk-based approach）的沙门氏菌监测与干预策略而非等同视之，对优化农业食品安全与公共卫生资源配置具重要意义。