中国食物链中多重耐药沙门氏菌的高风险抗菌药耐药基因组特征：至2050年的预测

《Nature Communications》：High-risk antimicrobial resistance genomic traits in multidrug-resistant Salmonella across China’s food chain: projections to 2050

　　

抗菌药耐药性（antimicrobial resistance, AMR）在全球上升，但塑造其基因组进化的因素仍知之甚少。研究人员调查了社会经济、气候与健康因素如何与中国食物链中多重耐药（multidrug-resistant, MDR）肠炎沙门氏菌（Salmonella enterica）的AMR传播相关联，绘制热点并预测基因组特征至2050年。研究人员分析了来自15个省份、跨度16年的2536个MDR基因组，覆盖患者、食品与非食品工作者、肉类、家畜及环境来源。机器学习框架将基因组监测与风险评估整合，以识别高风险AMR特征（临床相关、MDR、可移动）。回归模型将特征流行率与社会经济、气候、健康和抗菌药物消费数据关联，用于2050年预测。84个高风险特征，包括碳青霉烯、黏菌素和β-内酰胺耐药，预计随抗生素使用、人口增长、医疗支出和城市化上升。预测具特征特异性：blaCTX-M-64与blaTEM-1可增加至20%，而blaCTX-M-65与blaCTX-M-14低于5%。氟喹诺酮耐药出现机制转变：gyrA 87突变耐药下降，质粒介导基因qnrB17、qnrS1上升。食品工作者成为潜在储库，特征流行率高于食品来源，独特携带最后一线抗菌药耐药基因blaTEM-30、blaNDM-5、blaCTX-M-64。将基因组与社会经济监测整合，可使AMR干预由被动转为主动预防。

研究背景方面，抗菌药耐药性（AMR）被视作全球重大公共卫生威胁，预计到2050年每年可造成约1000万死亡。AMR的出现与扩散由MDR基因组特征、可移动遗传元件（mobile genetic elements, MGEs）、跨物种传播及One Health互作共同驱动，并深受抗生素使用、城市化、人口密度、气候、医疗支出与卫生条件等社会经济和健康因素影响。在亚洲尤其中国，农业抗生素广泛使用、快速城市化和季节性气候变异使这一动态更突出。食源性感染中，非伤寒沙门氏菌（nontyphoidal Salmonella, NTS）负担沉重，但AMR如何跨越人类、动物、环境、食品这一One Health食物链传播仍不清楚。既有研究多限于单城市、表型耐药或缺乏风险分层与预测框架，因此研究人员开展该项横跨15省、16年、覆盖全食物链的高分辨率基因组与多维指标融合研究，论文发表于《Nature Communications》。

关键技术方法上，研究人员构建2536个MDR肠炎沙门氏菌全基因组序列队列，来源包括无症状食品与非食品工作者、有症状患者、猪肉与鸡肉零售肉品、家禽与养殖场环境；另设306个非MDR分离株对照。方法分为五阶段：生物信息学识别AMR特征与MGE位置；风险评估框架按传播可能性和临床相关性将特征分为Rank I–V，高危为Rank IV–V；贝叶斯混合效应逻辑回归与PopPUNK谱系分层机器学习，以抗生素药敏表型为输出，随机效应控制谱系、省份、年份和来源；将205个高危特征与96个国家层面社会经济、气候、健康、抗生素消费指标做Pearson相关；用回归、蒙特卡洛模拟、平稳性检验（Augmented Dickey-Fuller、Kwiatkowski-Phillips-Schmidt-Shin）将特征流行率预测至2050。

研究结果部分保留小标题如下。一、A wide array of Salmonella serovars with phylogenetic intermixing are spread throughout the food chain：研究人员经系统发育分析发现88个血清型、80个序列型（sequence type, ST），以S. Typhimurium单相变体4,[5],12:i:-和S. Typhimurium为主；按血清型聚类明显，但来源、MDR状态、地理位置和年份无显著聚类，说明遗传背景在食物链中高度混合。

二、Antimicrobial resistance is disseminated throughout the food chain, via ARGs residing primarily in mobile elements：研究人员识别出421个AMR特征，88%主要位于可移动位置，包括质粒、插入序列（insertion sequences, ISs）相关结构；40个特征存在于全部七类来源，食品工作者特征丰富度最高，并独有134个基因组特征，含黏菌素mcr-1.1_IS1A、碳青霉烯bla_NDM-5等最后一线耐药特征，且93.2%为可移动。

三、Genomic risk scoring reveals food workers and symptomatic humans (diarrhoea patients) carry the highest-risk AMR traits：研究人员经风险评分与重抽样校正发现，95%的高危特征出现于食品工作者，高于食品来源；非食品工作者也高于食品来源但效应较小，有症状患者反而较低。效应梯度为食品工作者7%＞非食品工作者4.3%＞有症状患者-5.9%，支持食品相关职业暴露梯度，但受广西采样地理混杂限制，应视为流行病学关联而非因果源效应。33个Rank V特征中31个见于食品工作者。

四、Machine learning predicts that high-risk AMR traits are strong predictors of observed AMR resistance/susceptibility profiles：研究人员训练18个抗生素特异模型，AUC均＞0.7；省份与年份随机效应方差最大。经OR≥2、95%置信下限≥1、交叉验证稳定度≥80%筛选后得339个唯一特征，其中205个为高危，98%至少在某分离株中可移动。模型在2021—2022内部验证集平均AUC 0.94，在广东独立外部队列中仍具判别力，但类别不平衡和数据集偏移会影响精度与校准。

五、Health, governance, socioeconomic factors, and antibiotic usage are associated with high-risk AMR trait prevalence, highlighting emerging global threats in the food chain：研究人员发现165个高危特征至少与一个指标｜r｜＞0.5，主要涉及社会经济、健康、抗生素使用和气候；101个可预测，84个预计上升。fluoroquinolone耐药中gyrA 87下降，qnrB17、qnrS1上升；tet(A)、tet(M)上升，tet(B)下降；bla_CTX-M-64、bla_TEM-1至2050年可达20%，bla_CTX-M-65低于5%。同一抗生素类内基因轨迹分化，说明需从类水平监测走向基因水平监测。

六、Indicator-indicator correlation and robustness analyses：研究人员对96个指标做相关与PCA，平均｜r｜仅0.237，同类指标聚团，跨类弱相关；主成分无单一主导轴，敏感分析保留95%特征与方向，证明指标冗余下结论稳健。

讨论部分总结：研究人员指出已有框架如Zhang等组学框架与WHO AWARE分类被纳入其临床相关性维度，但本研究框架新增AMR特征—MGE组合、实证跨宿主与跨血清型流动性、以及中国食物链情境校准三项贡献。食品工作者作为潜在储库但并非因果证据；氟喹诺酮由染色体突变转向质粒介导耐药，意味着仅管控促生长用抗生素不足，须针对PMQR基因监测。同一类中基因命运不同，政策不能只看抗生素类别。局限包括来源—省份—年份非平衡、国家层指标难捕捉省份轨迹、相关非因果、短读长质粒重建误差、职业筛查重复采样等。

研究结论翻译：将沙门氏菌AMR基因组监测与全国层面社会经济、气候、健康数据整合，可识别未来25年中国食物链中易扩散的临床相关AMR特征及关键指标。同一抗生素类并非同等风险：bla_CTX-M-64与bla_TEM-1至2050年最高增至20%，bla_CTX-M-65低于5%；tet(A)、tet(M)增而tet(B)减；氟喹诺酮中gyrA 87突变耐药降，qnrB17、qnrS1等质粒介导耐药升。食品工作者携带更多高危特征并独有最后一线耐药基因，是潜在传播热点。AMR特征多位于可移动遗传元件，可水平转移给其他革兰阴性食源性病原体，提示威胁超越肠炎沙门氏菌本身。研究把描述性食品工作者带菌图景转为可预测、有机制依据的靶向干预框架，支持将基因组与社会经济监测融合，使AMR治理由被动反应转向主动预防。