IF50.5｜基因组学揭示西非猴痘的人源持续传播与动物源反复溢出

喀麦隆病例源于反复动物传人，尼日利亚疫情则已建立持续人传人链。

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

研究背景

猴痘病毒（MPXV）分为两个主要分支，其中分支IIb历史上局限于西非动物种群。自1970年首次在人类中发现至2017年，猴痘病例多为散发，人际传播有限。然而，2017年尼日利亚和2018年喀麦隆在三十余年后再次报告病例，尼日利亚疫情被认定为持续人传人流行，而喀麦隆的病例驱动因素尚不明确。2022年，分支IIb的B.1谱系引发全球多国猴痘疫情，其基因组表现出APOBEC3胞苷脱氨酶驱动的突变特征，提示病毒在人类中持续传播并积累适应性突变。此前研究估计MPXV于2016年在尼日利亚进入人群并隐秘传播，但缺乏足够基因组数据来解析两国的人畜共患与持续人际传播动态。

科学问题

卡在哪里： 该地区MPXV全基因组序列有限，导致几个关键问题悬而未决：尼日利亚和喀麦隆的人传人程度如何？动物源溢出事件在两国的角色是什么？尼日利亚持续人流行病的驱动因素和地理起源尚不清楚。此外，最接近尼日利亚人流行谱系（hMPXV-1）的动物源外群一直未被找到，限制了对病毒出现时间和地点的推断。

本文要回答： 通过测序2018-2023年尼日利亚和喀麦隆的118个MPXV基因组，解析两国人畜共患与持续人际传播的相对贡献，追溯hMPXV-1的动物源祖先、出现时间和地理起源，并揭示尼日利亚境内病毒传播的时空动态。

技术路线

作者主要用了：全基因组测序、系统发育分析、APOBEC3突变特征分析、贝叶斯分子钟、离散与连续系统地理学、广义线性模型。

作者首先对2018-2023年尼日利亚（n=109）和喀麦隆（n=9）的MPXV临床样本进行全基因组测序，获得118个近全长基因组。随后将这些序列与公共数据库中所有高质量分支IIb基因组合并，构建时间分辨的系统发育树。通过分析APOBEC3特征突变的比例，区分人传人谱系与动物源感染。进一步采用贝叶斯两阶段模型估计hMPXV-1的动物源祖先时间和人类出现时间。利用离散和连续系统地理学重建病毒在尼日利亚各州间的扩散路径，并通过广义线性模型评估影响传播的协变量。

图 1　Fig. 1。图1展示了采样地理分布和系统发育框架。a-b显示喀麦隆样本主要来自西南和西北边境地区，尼日利亚样本集中在南部各州。c显示2017年后南部为疫情中心，2022年起南北均出现病例回升。d-e的系统发育树显示105个尼日利亚序列散布在hMPXV-1的6个亚谱系中，而所有喀麦隆序列形成独立的动物源分支。f显示hMPXV-1内部APOBEC3突变比例高达74%，而动物源分支仅9%，直观区分了两种传播模式。

核心亮点

亮点 1｜喀麦隆病例均为动物源感染

系统发育分析显示，所有9个喀麦隆序列形成一个与hMPXV-1及其已知动物源外群KJ642617并列的基部分支，表明它们来自独立的动物源溢出事件，而非人传人传播。APOBEC3突变特征分析进一步证实：hMPXV-1内部约74%的重建SNP为APOBEC3型突变，而喀麦隆序列和4个尼日利亚动物源序列所在分支的APOBEC3突变比例仅为9%。在总数据集中，13个动物源感染占约11%，其中喀麦隆病例100%为动物源，尼日利亚约3.7%为动物源。这一结果确认了喀麦隆近期病例激增并非由持续人传人驱动，而是反复的动物传人事件。

图 2　Fig. 2。图2呈现时间分辨的系统发育和关键时间节点估计。a显示分支IIb的全球时间树，新发现的动物源外群Zx标注在hMPXV-1旁。b详细展示Zx与hMPXV-1的关系，估计Zx与hMPXV-1的共同祖先在2013年11月，向持续人传人的转变在2014年8月，hMPXV-1的最近共同祖先在2015年8月。这些时间点将病毒出现时间提前了约13.5个月，并明确了动物源祖先的流行窗口。

亮点 2｜跨境森林生态系统共享动物病毒库

两个尼日利亚动物源序列（来自南部阿夸伊博姆州）与喀麦隆序列聚集在同一新发现的动物源谱系中，该谱系与hMPXV-1形成姐妹群。阿夸伊博姆序列与喀麦隆西南部Mbongue序列相差25-29个SNP，但共享共同祖先。贝叶斯重建估计该喀麦隆谱系与KJ642617的共同祖先于1966年2月左右在动物种群中流行，而阿夸伊博姆序列与Mbongue序列的分化时间约为2009年7月。这些结果支持病毒在尼日利亚-喀麦隆边境的动物种群中历史性和近期双向传播，跨境森林生态系统可能承载着共享的动物病毒库。

图 3　Fig. 3。图3展示hMPXV-1在尼日利亚的时空传播。a的最大分支可信树按来源地区着色，显示河流州为起源。b显示河流州贡献了75次病毒输出，主要流向其他南南州、东南和西南地区。c显示早期引入均源自河流州，北部州的引入主要发生在2022年疫情回升后。d的连续系统地理学重建显示病毒在2017年9月前已扩散至500公里外的多个州，包括拉各斯和联邦首都区。

亮点 3｜hMPXV-1动物源祖先追溯至尼日利亚南部

两个来自阿比亚州的尼日利亚动物源序列（命名为Zx）与hMPXV-1形成姐妹谱系，将hMPXV-1与KJ642617之间的长茎枝从27个SNP缩短至8个SNP。利用Zx序列作为外群，两阶段贝叶斯模型估计Zx与hMPXV-1的共同祖先于2013年11月下旬在动物种群中流行，hMPXV-1向持续人传人的转变发生在2014年8月，比先前估计早约13.5个月。hMPXV-1的最近共同祖先时间为2015年8月。Zx与KJ642617均在阿比亚州采样，时间跨度50年，提示hMPXV-1前体谱系在尼日利亚南部动物种群中循环了约50年。

图 4　Fig. 4。图4分析各州传播链的持续性。a显示河流州的传播链从出现起持续存在，而其他州的传播链多由河流州引入。b显示河流州传播链持续时间最长，三角洲和巴耶尔萨的早期传播链持续约4.5和3年。c显示2022年疫情回升期间传播链数量显著增加。d显示南南地区传播链持续性动态变化，但河流州始终以持续传播为主，其他州则依赖反复引入。

亮点 4｜河流州是尼日利亚疫情源头和主要扩散中心

离散和连续系统地理学重建均强烈支持hMPXV-1起源于河流州（后验概率0.97）。河流州也是州际病毒输出的主要来源，估计贡献了75次病毒引入（95% HPD: 70-84）。早期疫情中所有引入均源自河流州，2016年病毒已传播至巴耶尔萨、伊莫和拉各斯等州。到2017年9月疫情宣布时，病毒已扩散至河流州以外500多公里的11个州。河流州的传播链从出现起持续存在，而其他州的局部疫情主要由河流州的反复引入驱动。广义线性模型显示，起源于河流州是病毒扩散的最强预测因子（贝叶斯因子>50），目的地人口密度也有正向关联。

图 5　Fig. 5。图5展示驱动病毒扩散的因素。a的广义线性模型系数显示，起源于河流州是病毒扩散的最强预测因子（贝叶斯因子>50），目的地人口密度也有正向关联（贝叶斯因子>15）。b的系统地理学重建突出河流州作为起源和主要扩散中心，其他州显示为灰色。这一结果排除了采样偏差的影响，确认了河流州在疫情中的核心地位。

亮点 5｜hMPXV-1在检测前已隐秘传播三年

结合时间估计和系统地理学重建，hMPXV-1于2014年8月进入人群，但直到2017年9月才在巴耶尔萨州首次报告病例，期间在尼日利亚南部隐秘传播和多样化超过三年。病毒在河流州持续循环两年以上才被检测到，并在疫情宣布前已建立多个州际传播链。这一发现强调了在动物-人类交界地区加强监测和早期诊断能力的必要性，以缩短从出现到检测的时间差。

生信可带走

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

• 方法栈： 全基因组测序、系统发育分析、APOBEC3突变特征分析、贝叶斯分子钟、离散与连续系统地理学、广义线性模型
• 公开数据： 原文未给出公开组学登录号
• 代码： 原文未给出公开 GitHub/GitLab 仓库
• 谁该点开原文： 从事病毒基因组学、系统发育分析、传染病流行病学或人畜共患病研究的科研人员。

带走一句

用APOBEC3突变特征区分人传人与动物源感染，结合系统地理学可精准定位疫情源头和扩散路径，为监测和干预提供依据。