克罗诺杆菌分离株的核心基因组与全基因组多位点序列分型

《Microbiology Spectrum》：Core genome and whole genome multi-locus sequence typing of Cronobacter isolates

　　

克罗诺杆菌属（Cronobacter）物种，尤其是阪崎克罗诺杆菌（C. sakazakii）和丙二酸盐克罗诺杆菌（C. malonaticus），是机会性病原体，与婴儿严重感染及高病死率相关。研究人员调查了全基因组测序（WGS）分析方法，具体包括7基因多位点序列分型（7-gene MLST）、核心基因组MLST（cgMLST）和全基因组MLST（wgMLST）对克罗诺杆菌分离株进行分型。研究人员分析了源自临床、食品和环境来源的743株克罗诺杆菌综合分离株集。研究人员还评估了高质量单核苷酸多态性（hqSNP）、cgMLST和wgMLST在聚集流行病学相关株及区分散发性阪崎克罗诺杆菌分离株中的表现。结果表明，cgMLST与wgMLST均能准确识别密切相关分离株，且与流行病学发现一致。基于等位基因的分析也与当前金标准hqSNP分析相当。研究人员的工作流同时输出7-gene MLST等位基因调用、克罗诺杆菌序列型及克隆复合体，可用于暴发调查中的历史比对。随着克罗诺杆菌感染近期在美国被列为全国法定报告疾病，研究发现证明了在PulseNet框架下基于WGS的方法可有效改善克罗诺杆菌暴发检测与响应策略。

研究背景方面，克罗诺杆菌属为革兰氏阴性兼性厌氧菌，2007年前归为肠杆菌科阪崎肠杆菌，含至少10种，临床以阪崎克罗诺杆菌与丙二酸盐克罗诺杆菌为主，可致婴儿（尤其低出生体重者）脑膜炎、败血症，病死率可达44%，存活者常有神经后遗症。既往分型依赖质粒谱、生物型、脉冲场凝胶电泳（PFGE），分辨率有限；近年全基因组测序（WGS）成为主流。2022年美国FDA大规模召回奶粉后，2024年起婴儿克罗诺杆菌感染被列为全国法定报告疾病，PulseNet需建立标准化WGS分型方案。本研究即在PulseNet框架内验证cgMLST与wgMLST相较金标准hqSNP的聚类效能，并明确质控阈值。

关键技术方法上，研究人员整合美国CDC肠道病实验室分支历史株、四家PulseNet公共卫生实验室送检株及FDA-CFSAN生物项目序列，共743株（639株阪崎克罗诺杆菌、43株C. turicensis、33株C. dublinensis、19株C. malonaticus、7株C. muytjensii、2株C. universalis），其中73株阪崎克罗诺杆菌来自2014至2023年18起食品与环境调查。Illumina测序后于PulseNet 2.0平台质控过滤与从头组装；采用7-gene MLST（PubMLST）、cgMLST（1865位点）、wgMLST（15739位点，含13867辅助基因组位点）调用等位基因；以Lyve-SET进行hqSNP最大似然树分析；以GrapeTree构建UPGMA树、iTOL注释；R语言ggplot2做质量指标与成对距离回归。

研究结果部分，序列质量评估显示，639株阪崎克罗诺杆菌中34株覆盖度低于40×，其中44%核心位点调用率低于95%但仍过85%门槛；覆盖度≥40×组中仅3%低于95%，证明设覆盖度≥40×可优化cgMLST调用。所有基因组均过Q值≥30、contig数与N50门槛。非阪崎种（如C. muytjensii、C. dublinensis）核心调用率略低（86%至95%）。

流行病学一致性分析表明，18起调查中每起内部遗传差均极小（hqSNP≤8、cgMLST≤5、wgMLST≤7），与流行病学归属完全吻合。

WGS系统发育比较中，cgMLST与wgMLST树均能在6个种水平分开聚类；18起调查各自成独立簇，调查间cgMLST差11至>200、wgMLST差26至>200（除调查10与18交叉外均>40）。组内成对距离显著低于组间（P<0.001）。

遗传成对距离与线性回归显示，cgMLST对hqSNP的R²=0.78、斜率0.62，wgMLST对hqSNP的R²=0.83、斜率0.87；cgMLST对wgMLST的R²=0.79、斜率1.20。Pearson系数cgMLST=0.88、wgMLST=0.91、cg/wg间=0.89，证明等位基因法与国际金标准高度线性相关。

MLST部分，7-gene MLST在阪崎克罗诺杆菌检出66种ST，ST4最常见（197株，占15/18起调查），次为ST1、ST64、ST40、ST8、ST17等；部分株因新等位基因未获ST。ST4克隆复合体为临床优势克隆。

讨论与结论翻译：侵入性克罗诺杆菌感染自2024年在美国列为全国法定报告疾病。为支持PulseNet下全国监测，研究人员以743株基因组将cgMLST、wgMLST对照hqSNP验证，明确Q≥30、覆盖度≥40×、长度4.0至4.9 Mbp、contig≤200、N50≥75000的阈值。三种方法在流行病学相关株上高度一致，支持cgMLST作为PulseNet常规监测标准、wgMLST与hqSNP供暴发高分辨使用。调查内cgMLST差≤10等位基因判为密切相关；若放宽至20或40会合并流行病学无关调查，故暂定10为界。ST4为临床与环境共有优势克隆，但环境分离株与临床株直接溯源链仍缺。研究局限包括仅含已测序株、环境采样不全、无直接环境-临床株确认试验。该工作为PulseNet纳入克罗诺杆菌、实现婴儿感染早检与风险消减提供依据。论文发表于《Microbiology Spectrum》。