IF50.5｜北美H5N1大流行的生态与传播驱动因素

野生迁徙鸟类是2021-2023年北美H5N1大流行的核心驱动者，反复将病毒引入家禽。

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

研究背景

自2021年底以来，高致病性H5N1分支2.3.4.4b病毒在北美引发了前所未有的动物大流行，导致大量野生鸟类、家禽和哺乳动物死亡。与2014/2015年北美H5N8疫情不同，此次疫情并未因大规模扑杀家禽而终止，且波及许多非典型宿主，如猛禽、猫头鹰和哺乳动物。此前在欧洲和亚洲的研究提示，野生鸟类在分支2.3.4.4b病毒的全球维持和传播中可能扮演越来越重要的角色，但北美地区野生与家养鸟类在传播中的相对贡献尚未得到系统评估。

科学问题

卡在哪里： 尽管北美H5N1疫情已造成巨大经济损失和生态影响，但病毒在北美大陆的引入次数、传播路径以及野生鸟类与家禽之间的传播动态仍不清楚。特别是，非典型宿主（如猛禽、鸣禽和哺乳动物）是否建立了新的传播链，以及后院禽类是否在疫情中充当了早期预警信号，这些问题缺乏基于系统发育动力学的定量证据。

本文要回答： 本文旨在利用贝叶斯系统发育地理学方法，重建2021-2023年北美H5N1病毒的引入、扩散和宿主间传播动态，明确野生鸟类在家禽疫情暴发中的作用。

技术路线

作者主要用了：贝叶斯系统发育地理学、离散性状扩散模型、马尔可夫跳跃、PACT持续性分析、BaTS聚类检验。

作者从GISAID下载了1,824条北美H5N1 HA序列，结合宿主、地理和迁徙行为等元数据，构建了多个数据集。首先通过贝叶斯系统发育重建和离散性状扩散模型，推断病毒从欧洲和亚洲引入北美的次数和时间。随后，利用迁徙飞行路线和宿主分类作为离散性状，量化病毒在不同地理区域和宿主类群间的传播速率和方向。为了评估采样偏差的影响，作者设计了滴定分析，逐步改变野生与家养鸟类序列的比例，并采用随机打乱检验（tip-shuffle）验证结果的稳健性。

图 1　Fig. 1。图1展示了全球1,927条H5N1序列的贝叶斯系统发育树，按大洲着色。主要的大西洋引入事件被标注，来自亚洲的多次引入在放大图中清晰可见，其中一次演化为D1.1谱系。该图直观呈现了北美疫情的多次独立起源。

核心亮点

亮点 1｜至少九次独立引入北美

系统发育分析显示，北美H5N1疫情由至少九次独立引入驱动。其中一次主要引入来自欧洲（2021年9-10月），随后在北美广泛传播；另外七次引入来自亚洲（2022年2-9月），主要进入太平洋飞行路线，但大多短暂存在。值得注意的是，其中一次亚洲引入在2024年底重新出现，形成了D1.1谱系，凸显了太平洋沿岸监测的重要性。

图 2　Fig. 2。图2a显示按迁徙飞行路线着色的系统发育树，插入的PACT分析显示大西洋和太平洋飞行路线中病毒谱系持续最久。图2b的马尔可夫跳跃速率显示密西西比到中部飞行路线的传播最频繁。图2c的地图箭头表明传播以东西向为主。

亮点 2｜迁徙飞行路线驱动快速扩散

病毒在北美内部的传播与候鸟迁徙飞行路线高度相关。序列按飞行路线显著聚类（AI=10.563, P=0.00199），相邻飞行路线间的传播速率约为远距离飞行路线间的10倍。传播方向以从东向西为主，东向西的跳跃频率约为西向东的4.4倍。野生迁徙鸟类被推断为最主要的传播源，其向其他类群的传播速率具有极高的贝叶斯因子支持（BF>3000）。

图 3　Fig. 3。图3a为655条序列的系统发育树，按宿主目着色，显示雁形目在多个深部节点具有高后验支持。图3b的马尔可夫跳跃速率显示雁形目向鸡形目和鸮形目的传播速率最高，而猛禽、鸣禽和哺乳动物主要作为接收端。

亮点 3｜Anseriformes是核心传播宿主

在七个宿主目类群中，雁形目（Anseriformes，包括鸭、鹅、天鹅）被推断为根部状态的概率比随机打乱数据高2-3倍，表明其在持续传播中的关键作用。雁形目向鸡形目（Galliformes）和鸮形目（Strigiformes）的传播速率最高（分别为17.81和13.51次/年，BF>100）。相反，猛禽、鸣禽和哺乳动物等非典型宿主主要作为病毒接收端（sink），传播链短且不持续。

图 4　Fig. 4。图4a为1,080条序列的系统发育树，按野生/家养状态着色，显示家禽序列形成多个小簇，散布在野生鸟类背景中。图4b显示从野生到家禽的传播次数远高于反向传播。图4c的PACT分析显示野生鸟类中病毒谱系持续更久。

亮点 4｜家禽疫情由野生鸟类反复引入

通过滴定分析，作者发现即使野生鸟类序列被大幅下采样，野生鸟类仍被推断为家禽疫情的主要来源。在最终数据集中，推断出106次从野生鸟类到家禽的引入事件，而家禽向野外的回传仅4次。家禽中的病毒谱系平均持续约4.5个月，而野生鸟类中的谱系持续约10个月。这表明2022-2023年的家禽疫情与2014/2015年由少数引入后在家禽间高效传播的模式截然不同。

图 5　Fig. 5。图5a为942条序列的系统发育树，按野生、商业家禽和后院禽类着色，显示后院和商业家禽的序列大多独立成簇。图5b的爆炸视图展示了从野生鸟类到后院和商业家禽的独立引入。图5c-e显示后院禽类的引入次数略多且时间更早。

亮点 5｜后院禽类感染早于商业家禽

对2022年1-5月数据的分析显示，后院禽类（backyard birds）的感染时间平均比商业家禽早约9.6天，且引入次数略多（后院43次，商业39次）。系统发育树中后院禽类序列更常位于基部，但大多数商业家禽的引入并非由后院禽类介导，而是独立来自野生鸟类。这一发现提示后院禽类可能作为野生鸟类病毒活动的早期预警信号。

生信可带走

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

• 方法栈： 贝叶斯系统发育地理学、离散性状扩散模型、马尔可夫跳跃、PACT持续性分析、BaTS聚类检验
• 公开数据： 原文未给出公开组学登录号
• 代码： 原文未给出公开 GitHub/GitLab 仓库
• 谁该点开原文： 从事病毒系统发育动力学、人兽共患病监测和禽流感防控的生信与流行病学研究者。

带走一句

做疫情溯源和传播动力学分析时，不能只看家禽数据；野生鸟类序列的时空分布和宿主分类信息对揭示传播链至关重要，且需用滴定和随机化检验来评估采样偏差。