IF50.5｜霍乱全球源头再定位：恒河盆地而非三角洲

孟加拉国霍乱演化独特，但全球传播多由印度输出。

这篇按生信解读号的读法拆：先看背景和问题，再看作者怎么做、做出了哪几条硬发现，最后把数据和代码收走。

研究背景

第七次霍乱大流行由霍乱弧菌O1 El Tor型7PET谱系引起，自1960年代起从孟加拉湾周边地区向全球扩散。此前研究认为恒河三角洲及孟加拉国是霍乱的全球源头，但缺乏对恒河盆地（横跨印度和孟加拉国）的长期纵向基因组监测。已知7PET在印度和孟加拉国存在两个主要谱系BD1和BD2，它们在核心基因组上差异较小，但在移动遗传元件（MGEs）上差异显著，包括噬菌体防御系统、抗生素抗性和毒力因子。然而，这些MGEs的演化动态及其对疾病严重性和全球传播的影响尚不清楚。

科学问题

卡在哪里： 此前缺乏对恒河盆地霍乱弧菌的长期纵向基因组数据，无法确定孟加拉国和印度在霍乱全球传播中的各自角色，也不清楚噬菌体防御元件的快速得失如何影响霍乱的临床严重性和跨区域传播。

本文要回答： 通过大规模纵向基因组监测，揭示恒河盆地霍乱弧菌的演化动态，确定全球传播的主要来源国，并阐明噬菌体防御元件与疾病严重性和传播潜力的关系。

技术路线

作者主要用了：全基因组测序、系统发育分析、祖先状态重建、宏基因组筛选、多变量逻辑回归。

作者收集并测序了2014-2023年孟加拉国（n=1,516）和2002-2023年北印度（n=794）的霍乱弧菌分离株，构建了包含全球3,112个已发表基因组的系统发育树。通过SNP距离划分亚谱系，利用TreeTime进行时间尺度和祖先状态重建，追踪谱系在国家和区域间的传播。同时，对关键基因和MGEs进行存在/缺失分析，并结合临床表型（米泔水样便、严重脱水）进行多变量逻辑回归，评估其与疾病严重性的关联。此外，对13,105个肠道宏基因组进行筛选，分析ICP1噬菌体的分布和共演化。

图 1　Fig. 1。图1展示了孟加拉国2014-2018年监测研究中各亚谱系和基因/MGE谱的时间丰度变化。BD2谱系在2014-2017年占主导，2018年被sBD1替换。基因/MGE的丢失事件（如PLE1、ICETET、ddmAB）在BD2中依次发生并固定。系统发育地理学显示，sBD1.3.1和BD2.2.2等亚谱系从达卡扩散，而PLE1- BD2.4.3和ΔVCA0455-VCA0459 BD2.4.1从吉大港扩散。

核心亮点

亮点 1｜孟加拉国霍乱演化受国界限制

系统发育分析显示，尽管孟加拉国和印度地理相邻且共享恒河三角洲，但霍乱弧菌在两个国家内的演化基本独立。BD2谱系在孟加拉国从2013至2017年占主导，但其中10个亚谱系仅在该国境内发现，未在其他国家出现。相反，印度加尔各答的菌株与北印度更相似，而非邻近的孟加拉国库尔纳。这表明霍乱传播主要受人口流动和国界限制，而非水文特征。

图 2　Fig. 2。图2比较了孟加拉国、北印度和加尔各答的亚谱系时间动态，显示两国谱系替换模式不同。南亚内部的传播事件中，孟加拉国和印度之间有40次传播，但多数未能持续。从南亚到其他地区的传播事件中，印度来源的事件远多于孟加拉国，且主要涉及携带ctxB7的sBD1谱系。

亮点 2｜孟加拉国霍乱弧菌快速丢失噬菌体防御元件

在孟加拉国，BD2谱系在2014-2018年间经历了快速的基因和MGE丢失，包括hlyA外毒素基因、ddmABC操纵子（噬菌体防御）、wbeT基因（血清型转换）、PLE1抗噬菌体岛和SXT-ICETET元件。这些丢失事件与疾病严重性增加相关：ddmA、wbeT和ICETET的存在与米泔水样便风险降低相关，而VCA0455-VCA0459和PLE1的存在与严重脱水风险降低相关。宏基因组分析显示，ICETET丢失后，霍乱弧菌在粪便中的相对丰度从16.8%增加到40.5%。

图 3　Fig. 3。图3展示了基因/MGE与临床严重性的关联。ddmA、wbeT和ICETET的存在与米泔水样便风险降低相关，而VCA0455-VCA0459和PLE1与严重脱水风险降低相关。宏基因组分析显示ICETET丢失后霍乱弧菌丰度增加。ICP1噬菌体类型在2018年从CRISPR-Cas型转变为odn型。PLE在传播事件中显著不足。

亮点 3｜2018年谱系替换伴随ICP1噬菌体类型转换

2018年，孟加拉国BD2谱系被来自印度的sBD1谱系完全替换。与此同时，共分离的ICP1裂解性噬菌体从携带CRISPR-Cas抗防御系统的类型转变为携带odn抗防御基因的类型。这种几乎同步的转换表明霍乱弧菌与ICP1之间存在紧密的共演化军备竞赛。sBD1携带的ICEGEN元件含有BREX抗噬菌体系统，可能最初对当地ICP1有效，但随后ICP1通过获得odn基因克服了该防御。

亮点 4｜印度而非孟加拉国是全球传播的主要源头

祖先状态重建显示，2003-2023年间，从印度传出的霍乱弧菌传播事件（n=359）远多于从孟加拉国传出的（n=11），尽管孟加拉国的样本数量是印度的两倍。sBD1谱系在印度获得ctxB7等位基因后，传播至尼日利亚、肯尼亚和海地，引发2010年海地疫情。PLE抗噬菌体岛在孟加拉国以外的传播事件中显著缺失，可能限制了孟加拉国菌株的全球传播潜力。

亮点 5｜PLE抗噬菌体岛限制霍乱弧菌全球传播

PLE1和PLE11抗噬菌体岛在孟加拉国以外的传播事件中显著不足。从印度传出的26个霍乱弧菌中没有一个携带PLE。此外，Inaba血清型（wbeT基因失活）的菌株也很少传播出孟加拉国。这表明PLE等噬菌体防御元件的维持可能以降低传播适合度为代价，从而将霍乱弧菌限制在孟加拉国境内。

生信可带走

这一块是给做分析的人用的，不是科普点缀。

• 方法栈： 全基因组测序、系统发育分析、祖先状态重建、宏基因组筛选、多变量逻辑回归
• 公开数据： 原文未给出公开组学登录号
• 代码： 原文未给出公开 GitHub/GitLab 仓库
• 谁该点开原文： 从事病原体基因组流行病学、噬菌体-宿主共演化或全球健康研究的科研人员。

带走一句

纵向基因组监测揭示霍乱弧菌在局部地区的快速MGE得失与全球传播潜力相关，提示应关注印度等上游地区的监测以预警全球疫情。