Molecular characteristics, fitness, and virulence of high-risk and non-high-risk clones of carbapenemase-producing Klebsiella pneumoniae
产碳青霉烯酶高风险与非高风险肺炎克雷伯菌克隆的分子特征、适应性及毒力
来源：Microbiology Spectrum, February 2024, Volume 12, Issue 2, DOI:10.1128/spectrum.04036-22
《微生物学谱》，2024年2月，第12卷第2期，DOI:10.1128/spectrum.04036-22
摘要
广泛耐药（XDR）肺炎克雷伯菌给全球医疗体系带来沉重负担，产碳青霉烯酶菌株治疗手段极为有限，部分序列型克隆在全球大范围流行，但此类克隆流行成功的分子机制尚不明确。本研究共纳入20株高风险克隆（HiRC）、10株非高风险克隆（non-HiRC）XDR肺炎克雷伯菌，探究驱动部分克隆流行病学优势的相关因素。通过核心基因组多位点分型（cgMLST）解析菌株克隆亲缘关系；利用生物信息学预测荚膜分型、O抗原、毒力基因、耐药基因、前噬菌体、CRISPR-Cas基因座；体外生长实验测定绝对与相对适应性；大蜡螟幼虫体内感染模型评估菌株毒力。结果显示，体外适应性、已知毒力因子、前噬菌体、CRISPR-Cas单独均无法解释特定克隆全球流行优势；但体内实验证实高风险克隆毒力显著高于非高风险克隆。该研究表明，肺炎克雷伯菌克隆的毒力与流行病学优势无法由本研究检测的单一性状解释，仍需后续更多实验深入探究。
关键词
抗菌耐药；适应性；竞争；流式细胞术；大蜡螟；肺炎克雷伯菌；毒力
研究目的
1. 对比产碳青霉烯酶XDR肺炎克雷伯菌高风险、非高风险克隆在基因组层面的差异，包括荚膜、O抗原、毒力基因、耐药基因、前噬菌体、CRISPR-Cas分布特征。
2. 体外生长与竞争实验，量化两类克隆绝对、相对适应性差异。
3. 利用大蜡螟动物模型，在体内水平对比两类克隆毒力强弱。
4. 明确基因组特征、体外适应性、体内毒力哪一因素可以解释高风险克隆全球大规模流行的流行病学优势。
研究思路
1. 菌株筛选与分型：收集30株多国来源产KPC/NDM/VIM型碳青霉烯酶菌株，依据暴发次数、报道国家划分为20株高风险、10株非高风险克隆；全基因组测序，cgMLST划分克隆群。
2. 基因组生物信息分析：使用Kleborate、CRISPR-Cas++、PHASTER工具，注释荚膜分型、脂多糖O抗原、毒力基因、耐药基因、CRISPR-Cas系统、前噬菌体数量与组成。
3. 体外适应性检测：采用芬兰Bioscreen C仪器测定菌株生长曲线，计算代时代表绝对适应性；构建YFP荧光标记高风险菌株，与未标记非高风险菌株共培养，流式细胞术测定选择系数表征相对适应性。
4. 体内毒力评价：大蜡螟幼虫腹腔注射菌液，48小时持续记录幼虫存活时间，采用Kaplan–Meier生存分析对比两类克隆毒力。
5. 关联与统计分析：对比克隆群、高低风险组各项指标，分析代时与毒力相关性，采用多种参数、非参数统计方法完成差异检验。
研究亮点
1. 同时结合全基因组生物信息、体外适应性、体内昆虫模型多维度系统对比高低风险产碳青霉烯酶肺炎克雷伯菌，突破仅单一表型或基因组分析的局限。
2. 证实常规基因组毒力元件、前噬菌体、CRISPR-Cas、纯体外生长适应性均不能解释高风险克隆全球流行，推翻单一性状驱动流行的传统猜想。
3. 首次在统一菌株库中证实高风险克隆体内毒力显著更强，明确体内毒力是区分两类克隆的关键表型差异。
4. 发现CG258流行克隆完全缺失完整CRISPR-Cas系统，该特征利于其摄取外源耐药移动元件，为该克隆广泛传播提供基因组线索。
5. 发现两类克隆生长-毒力关联模式存在本质区别：非高风险菌株毒力与代时负相关，高风险克隆毒力不依赖生长速率，拥有独立致毒通路。
可延伸的方向
1. 添加抗生素选择压力开展体外竞争实验，模拟临床抗生素环境下高低风险克隆适应性差异。
2. 开展小鼠等高等动物体内感染实验，验证大蜡螟模型毒力结论是否适用于哺乳动物。
3. 深度挖掘CG258缺失CRISPR-Cas促进耐药质粒捕获的分子机制，构建基因敲除回补株验证功能。
4. 转录组、代谢组对比高低风险克隆，挖掘现有注释未覆盖的新型毒力、适应性相关基因。
5. 扩大菌株样本，纳入更多地区、不同碳青霉烯酶亚型菌株，构建更大流行病学数据集。
6. 筛选靶向高风险克隆特有毒力通路的抑制剂，为院感防控提供新靶点。
测量的数据及研究意义
1. cgMLST克隆分群、荚膜分型、O抗原、毒力评分、耐药评分、前噬菌体数量、CRISPR-Cas完整度数据，来自表1、表2、图1
研究意义：1. 高低风险克隆荚膜、O抗原、已知毒力、耐药基因总量无统计学差异，证明上述经典分子特征无法区分两类克隆流行潜力；2. 所有菌株均携带多个前噬菌体，但前噬菌体数量、来源基因占比在两组间无差异；3. CRISPR-Cas完整度在不同克隆群存在显著区别，全球流行CG258克隆无完整CRISPR，利于外源耐药基因获取，是其广泛传播潜在基因组基础。
2. Bioscreen仪器测定各菌株代时（绝对适应性），来自图2
研究意义：高低风险克隆体外代时无显著差异，单纯无抗生素体外环境下生长快慢不能解释高克隆的流行病学优势，说明医院抗生素选择压力、体内微环境才是筛选关键。
3. 共培养竞争实验选择系数（相对适应性），来自图3
研究意义：CG258高风险克隆在无药体外培养中对比多数非高风险菌株无竞争优势，甚至被部分非高风险菌株抑制，进一步证明纯体外环境无法复刻临床流行优势，必须结合体内宿主因素。
4. 大蜡螟幼虫生存曲线、标准化毒力值，来自图4、图5
研究意义：1. 高风险克隆感染后幼虫中位存活时间更短、48h存活率仅2.2%，显著低于非高风险克隆，体内毒力存在明确组间差异；2. CG14克隆毒力显著高于CG258等克隆群，说明毒力差异存在克隆层级特异性；3. 仅体内模型可区分两类克隆，体外基因组、生长指标均无法实现区分。
5. 代时与毒力相关性统计数据
研究意义：非高风险菌株增殖越快毒力越强，高风险克隆毒力与生长速度无关，提示高风险克隆进化出不依赖快速增殖的独特致病机制，是其院内暴发重要原因。
结论
1. 荚膜、脂多糖、已知毒力基因、前噬菌体、CRISPR-Cas系统、无抗生素体外生长适应性单一指标，均不能解释产碳青霉烯酶高风险肺炎克雷伯菌全球大规模流行的现象。
2. 大蜡螟体内感染模型证实高风险克隆毒力显著高于非高风险克隆，体内致病性是两类克隆最明确的区分表型。
3. CG258等主流高风险克隆完全缺失完整CRISPR-Cas系统，利于捕获耐药移动遗传元件，是其广泛传播的重要基因组特征，但不足以单独解释其流行优势。
4. 非高风险克隆毒力与生长速率负相关，高风险克隆毒力不受增殖速度调控，二者生存进化策略完全不同。
5. 高风险克隆的流行病学成功是基因组、宿主体内微环境、抗生素压力等多重因素协同作用的复杂结果，单一性状无法主导，后续需结合抗生素暴露、高等动物模型开展深入机制研究。
使用芬兰 Bioscreen仪器测量数据的研究意义
1. Bioscreen C可高通量、24小时不间断自动检测96孔板菌株OD600，每4分钟采集一次数值，完整捕获菌株延滞期、对数期、稳定期全生长动态，精准计算指数生长期代时，客观量化菌株绝对适应性，避免人工单点检测带来误差。
2. 标准化生长条件统一所有菌株培养环境，平行设置生物学重复，实现高风险、非高风险、不同克隆群生长参数客观横向对比，得出两组体外生长无显著差异的可靠结论。
3. 代时数据作为核心体外适应性指标，与后续体内毒力、共培养竞争结果形成对照，明确体外生长优势不等于体内流行优势，区分体外、宿主内两套筛选体系的差异。
4. 代时数值用于相关性分析，挖掘出两类克隆“生长-毒力”关联模式完全不同，揭示高风险克隆拥有独立致病通路，为解释其院内高致病性提供表型证据。
5. 生长曲线数据为后续竞争实验、体内攻毒接种浓度标准化提供依据，统一各菌株接种量，消除生长快慢对竞争、毒力实验的干扰，保障多组实验数据可比性。