《SCIENCE ADVANCES》：Defining epidemiologically relevant units of cholera transmission in sub-Saharan Africa
编辑推荐：
关于非洲霍乱传播动态仍存在大量未知。研究人员联合运用传统系统地理学技术（phylogeographic techniques）与整合流行病学及基因组数据的创新方法，识别出3至4个独立的"霍乱传播单元"（cholera transmission units），单
关于非洲霍乱传播动态仍存在大量未知。研究人员联合运用传统系统地理学技术（phylogeographic techniques）与整合流行病学及基因组数据的创新方法，识别出3至4个独立的"霍乱传播单元"（cholera transmission units），单元内部的跨边境霍乱传播更为频繁。这些单元可归并为两大传播集群（transmission clusters），其一位于西非及毗邻国家，其二横跨中部、东部和南部区域。模拟研究（simulation studies）表明，上述集群可用于评估引入的霍乱菌株（cholera strains）国际传播风险，并基于暴发地点识别国际传播速度的异质性（heterogeneities）。研究方法进一步允许在缺乏测序数据（sequencing data）的情况下推断特定谱系（specific lineages）所致的霍乱发病率（cholera incidence）。该工作对霍乱防控（cholera control）具有明确意义。
**研究背景与问题**
霍乱是由产毒霍乱弧菌（Vibrio cholerae）引起的急性胃肠道感染，以水样腹泻和快速脱水为特征，未经治疗可迅速致死。尽管改善水和卫生基础设施、疫苗以及口服补液盐（oral rehydration salts）有效治疗等防控工具已被证实有效，霍乱仍造成大量发病和死亡，全球每年估计有130万至400万病例和2.1万至14.3万例死亡。该病在孟加拉湾周边国家广泛流行，但全球大部分霍乱死亡集中于撒哈拉以南非洲。自20世纪60年代大流行O1 El Tor型霍乱弧菌从东南亚传入以来，霍乱在撒哈拉以南非洲呈地方性流行（endemic）。该区域发病率呈季节性模式，但大规模的、常发生于非流行季的暴发并不罕见，可造成毁灭性影响。自然灾害和社会动荡常触发重大霍乱疫情，但无法解释霍乱在全非洲的持续存在与传播。为揭示霍乱传播的潜在驱动因素，研究人员转向霍乱弧菌全基因组测序（whole-genome sequencing），发现自1970年以来至少有17次O1 El Tor型霍乱弧菌从东南亚传入非洲（谱系AFR1至AFR17）。近期研究对霍乱传入非洲大陆后的传播方式提供了一定见解，发现邻近国家间频繁跨边境传播、国家内菌株共循环（strain cocirculation）以及流行间歇期谱系的长期持续存在等证据。然而，这些研究在地理和时间范围上有限，部分原因是可用序列数据的覆盖有限且不均衡。因此，关于非洲大陆内霍乱传播的地理结构、谱系持续存在以及新引入谱系可能如何传播等问题仍未解决。研究人员因此开展本研究，旨在界定流行病学相关的霍乱传播单元并深入理解非洲霍乱传播动态。该论文发表在《SCIENCE ADVANCES》。
**关键技术方法**
研究整合了1970至2023年来自GenBank和VibrioWatch的霍乱弧菌序列数据（含国家、年份和谱系标注元数据），以及世界卫生组织（WHO）成员国报告的年度疑似和确诊病例数据。主要方法包括：（1）基于1762条第三波大流行（AFR9及以后）霍乱序列（1989至2023年，来自36个非洲国家）的最大似然系统发育分析（maximum likelihood phylogeny），使用bacpage流程、Gubbins v2.3.4和IQ-TREE2构建，并以BEAST X v10.5.0进行时间校准和离散状态祖先重建（discrete state ancestral reconstruction）；（2）创新性地开发隐马尔可夫模型（hidden Markov model, HMM），整合流行病学数据与序列衍生谱系信息，推断缺失谱系及跨边境连通性（cross-border connectivity）；（3）利用Louvain社区检测算法（community detection）和共识聚类算法（consensus clustering）定义霍乱传播单元；（4）基于HMM后验参数进行霍乱传播模拟（simulation），评估引入谱系的1年暴发风险和传播速度；（5）通过交叉验证和敏感性分析评估模型稳健性。
**研究结果**
**系统地理学分析（Phylogeographic analysis）**
研究人员利用系统地理学技术基于1762条第三波大流行霍乱序列推断非洲内部国际霍乱传播速率。观察到地理转换密度和幅度最高的为空间邻近国家之间，远距离地点间转换支持较少，与既往区域研究一致。国家质心间每增加1000公里距离，有支持转换的概率下降88%（OR 0.12；95% CI 0.07, 0.17）。研究人员识别出4个独立但相互联系的霍乱传播单元，定义为单元内与单元外霍乱传播概率之比。单元内转换率平均为单元间转换率的13.5倍（95% HPD区间：9.7, 17.6）。这些传播单元可进一步归并为西部集群（phyW1和phyW2）与横跨中东部南部集群（phyCE和phyS）。
**整合流行病学与谱系数据的隐马尔可夫建模方法（Combining epidemiologic and lineage data in a hidden Markov modeling approach）**
研究人员开发了创新建模方法，利用流行病学数据和序列衍生霍乱谱系推断缺失谱系并深入理解非洲传播动态。使用HMM综合1970至2023年45个非洲国家的发病率数据（含当地同期序列谱系信息的国家-年份占26%），建模霍乱潜在传播与观察过程，估计谱系存在的隐藏状态概率及谱系在国家间移动的概率。
**推断感染谱系的存在、流行率与持续存在（Inferring presence, prevalence, and persistence of infecting lineages）**
利用该方法推断撒哈拉以南非洲所有报告霍乱病例的感染谱系。估计自1970年以来非洲观察到的霍乱病例中46%（240万例）可归因于两个谱系：AFR10（33% [95% CI: 31, 35]）和AFR9（14% [95% CI: 12, 15]），无其他谱系占报告病例超过10%。分析还揭示局部动态，如提示AFR5在20世纪90年代末和21世纪初被AFR10在东非的置换是跨国现象，并显示AFR10在病例稀少且无序列数据期间在这些国家的持续存在。
**定义霍乱传播单元——HMM模型（Defining cholera transmission units—HMM model）**
利用HMM模型推断的连通性，识别出3个独立霍乱传播单元，与系统地理学分析基本一致：一个涵盖西非及毗邻国家（hmmW），一个横跨中非和东非（hmmCE），一个主要位于南部非洲（hmmS）。这些单元进一步归并为西部集群（hmmW）和横跨中东部南部集群（hmmCE和hmmS），与系统地理学分析发现的主要集群完全相同。国家间谱系移动在传播单元内可能性为单元间的7.8倍（95% CI 6.5, 9.4）。方法间一致性证实了西部霍乱流行与非洲其他地区的独特性，并显示传播模式随时间的稳定性。
**模拟霍乱后续传播（Simulating onward spread of cholera）**
利用HMM模型对新引入谱系传播的模拟进一步支持霍乱传播单元的重要性。对45个国家各2000次引入的模拟发现，模拟引入后第一年内，同一传播单元内国家经历后续霍乱暴发的可能性为单元外国家的10.3倍（95% CI 10.1, 10.5）。新引入菌株到达同一传播单元内国家的平均时间比单元外国家低3.3年（95% CI ?3.2, ?3.2）。hmmCE与hmmS之间的传播也比这些单元与hmmW之间更可能。即使在传播单元内，各国预测传播范围和速度也存在显著差异，如加纳引入预计1年内导致hmmW内4.69个其他国家暴发，而邻国喀麦隆仅预计导致0.58个。
**讨论与结论**
研究结果丰富了对非洲霍乱传播的理解，对霍乱风险评估和国家层面防控具有明确意义。霍乱传播单元的存在为预判非洲大陆其他地区霍乱暴发中哪些国家风险最高以及协调霍乱防控的自然国家分组提供了框架。传播速度模拟提示某些国家（如莫桑比克、刚果民主共和国）的暴发更可能播散至非洲大陆其他地区，而其他国家（如塞内加尔、科摩罗）的暴发可能保持自限性，有助于完善全球疫苗储备分发风险评估。研究局限性包括：系统地理学分析受数据稀疏性限制；HMM推断亦受数据可用性限制且损失了同一谱系暴发间遗传关系的信息；推断基于报告疑似病例，为真实发病率的非完美度量；使用国家级数据和质心距离可能掩盖霍乱动态的关键方面（如刚果民主共和国东部和西部疫情相对独立）；未明确考虑免疫或菌株间竞争；HMM分析覆盖较长时间段，传播动态、监测强度和病例报告实践可能已变化。尽管存在局限，验证实验显示建模方法的稳健性：交叉验证中HMM模型推断检测谱系的召回率为73%（SD = 7），后期更密集采样谱系恢复更一致（AFR10和AFR13召回率分别为94%和87%）；模拟到达时间与观察到达时间呈中度一致（Pearson r = 0.49至0.55）。该工作展示了整合分子与流行病学数据获取疾病生态学重要见解的力量，观察到的明确霍乱传播单元及谱系持续存在和后续传播推断既可指导公共卫生行动也可推动进一步科学探究，并希望为其他病原体研究中分子与流行病学数据的更紧密整合提供模板。研究结论：研究人员联合运用传统系统地理学技术与整合流行病学及基因组数据的创新方法，识别出非洲撒哈拉以南地区3至4个独立霍乱传播单元，可归并为西部与中东部南部两大传播集群；模拟研究表明这些集群可用于评估引入霍乱菌株的国际传播风险并识别传播速度异质性；研究方法允许在缺乏测序数据情况下推断特定谱系所致霍乱发病率；该工作对霍乱防控具有明确意义。
