Genotype-specific ecological and environmental drivers of HPAI H5N1 spread in wild birds in France, 2021-2023

技術的要約：フランスにおける野生鳥類へのHPAI H5N1拡散の遺伝子型特異的な生態学的および環境的要因、2021-2023年

問題提起
高病原性鳥インフルエンザ（HPAI）H5N1ウイルス（系統2.3.4.4b）は、2020年以降、世界的なパンズーチ（家畜伝染病）を引き起こしており、野生鳥類、家禽、および哺乳類に深刻な影響を与えている。野生鳥類は、地球規模および地域規模の両方においてウイルスの拡散を駆動する中心的存在であると認識されているが、異なるウイルス遺伝子型の拡散を形作る具体的な生態学的および環境的要因については、依然として理解が進んでいない。先行研究では、土地被覆（例：湿地）や宿主グループ（例：カモ科、カモツル科）との一般的な関連性が特定されているが、異なるウイルス系統間での拡散パターンの不均一性や、ドライバーが空間的・時間的にどの程度変化するかについては、十分に特徴付けられていない。フランスは、2021年から2023年にかけて重大なエピゾーチック（動物流行）の影響を受け、かつ高いシーケンシング・カバレッジを維持していたことから、これらのダイナミクスを調査するための重要なケーススタディとなる。

手法
本研究では、連続的な系統地理学的枠組みを用いて、フランス国内で循環している3つの主要なH5N1遺伝子型（EA-2020-C、EA-2021-AB、EA-2022-BB）の時空間的な広がりを再構築した。

• データキュレーション： 2021年10月から2023年7月までの期間の、フランス国立参照研究所（ANSES）およびGISAIDから取得した全ゲノム配列を用いた。データセットには野生鳥類および家禽の配列が含まれており、特に野生鳥類から家禽への一次導入に焦点を当て、その後の農場間伝播は除外した。
• 系統再構築： 遺伝子型はEURLの分類を用いて分類した。最大尤度樹はIQ-TREEを用いて推定された。主要なフランス国内の単系統グループ（単一の導入イベントに起因するクラスター）を特定するために、離散形質解析（BEAST v1.10.4）を用いた。
• 連続的な系統地理学： 特定された各単系統グループに対し、Cauchy relaxed random walk (RRW) モデルの下でBEASTを用いた連続的な系統地理学的推論を行い、ウイルスの拡散を再構築した。これにより、km²/year単位の加重拡散係数（WDC）の推定が可能となった。
• 生態学的および環境的ドライバー： 3km解像度の変数群を分析した。これには以下が含まれる：
  • 宿主変数： 症例の80%以上を占める種の出現確率（eBirdより）および、グループレベルの種豊かさ（カモ科、カモツル科、水鳥、海鳥、猛禽類、陸鳥）。
  • 土地被覆変数： 湿地、農業用地、森林、人工地、重要鳥類・生物多様性領域（IBA）、および家禽農場の位置。
  • リスクゾーン： 既存のフランスのリスクゾーン（野生鳥類による導入リスクを示すZRP、および農場間伝播リスクを示すZRD）。
• 関連性テスト： seraphimパッケージを用い、ウイルス系統の所在（系統樹内のノード）と環境変数との間の正および負の関連性をテストした。証拠の強さを評価するためにベイズ因子（BF）を算出した（BF > 20：強い、BF > 150：非常に強い）。
• リスクマッピング： 遺伝子型間で最も一貫した強い関連性を示した変数を組み合わせ、複合的な生態学的リスクマップを作成した。その後、このマップを2021-2023年および2024-2025年の観察された症例分布に対して検証した。

主な結果
本研究は、H5N1の拡散における空間的および宿主のダイナミクスが、遺伝子型間で顕著な不均一性を示すことを明らかにした。

1. 遺伝子型 EA-2020-C (遺伝子型 C):

   • ダイナミクス： 強い沿岸部のシグネチャーを示し、主に北ガンの繁殖期（2022年6月〜9月）に循環した。
   • 宿主： 第2波は北ガンによって支配されていたが、この系統はカモ科、カモツル科、水鳥、および海鳥の種豊かさと非常に強い正の関連を示した。
   • 環境： 人工地およびIBAと強く関連し、農業用地とは負の関連を示した。
   • 驚き： ウイルスの所在は、この特定の系統において感染報告されていない種（例：マガモ、アジサシ）とも強く関連しており、これは生息地の共有（生態学的インターフェース）または未検出の感染を示唆している。

2. 遺伝子型 EA-2021-AB (遺伝子型 AB):

   • ダイナミクス： フランス全土で最も広範なマルチホスト拡散を示し、3つの異なる単系統グループが異なる拡散率（最大約3,800 km²/週）を示した。
   • 宿主： 幅広い宿主（カモ科、カモツル科、猛禽類、家禽）に影響を及ぼし、特定の宿主への依存は見られなかった。
   • 環境： 湿地およびIBAと強く関連した。最初の単系統グループは、最も広い宿主範囲と最も強い生態学的関連を示した。

3. 遺伝子型 EA-2022-BB (遺伝子型 BB):

   • ダイナミクス： 主にカモツル科（アジサシ類）に関連しており、2023年2月から7月の間に循環した。
   • 宿主： カモツル科（アジサシ類およびカモメ類）に対して高い特異性を示したが、他の水鳥との関連も見られた。
   • 環境： IBA、湿地、およびZRPと強く関連し、農業用地および森林地帯とは負の関連を示した。

4. 生態学的リスクマッピング:

   • 8つの変数が複合リスクマップのために保持された：5種の鳥類（マガモ、オオハクチョウ、アジサシ、ヨーロッパアジサシ、コエジロ）および3つの鳥類グループ（カモ科、カモツル科、水鳥）。
   • 得られたリスクマップは、湿地および主要な河川ネットワークに集中した高リスクエリアを特定した。
   • 検証： このマップは、低リスク閾値（>0.125）において、2021-2023年期の症例の78%、および2024-2025年のエピゾーチックの76%を捉え、現在のフランスのZRP（それぞれ34%と56%を捕捉）を上回る性能を示した。
   • ZRP vs. ZRD： ウイルス系統は、遺伝子型BBにおいてのみZRP（野生鳥類リスク）と関連しており、ZRD（農場間伝播リスク）とは関連していなかった。これは、野生鳥類のダイナミクスが家禽へのスピルオーバーではなく、野生鳥類内の伝播によって駆動されているという仮説を支持している。

意義および主張
本論文は、H5N1の野生鳥類における拡散は高度に不均一であり、遺伝子型ごとに異なる宿主、環境、およびウイルス学的要因の複雑な相互作用によって形作られていると主張している。主な貢献は以下の通りである：

• 遺伝子型特異的なドライバー： 本研究は、異なる遺伝子型が異なる空間的ダイナミクスと宿主との関連性を示すことを実証しており、HPAIの拡散に単一の均一なドライバーが存在するという概念に異を唱えている。
• 生態学的インターフェース： 感染報告されていない種との関連が頻繁に検出されることは、共有された生息地（生態学的インターフェース）がウイルスの循環において重要な役割を果たしており、それが伝播経路を隠蔽しているか、あるいは未検出の感染を示唆していることを示唆している。
• サーベイランスへの示唆： 今回の知見は、主にカモ類の分布に基づいた現在のリスクゾーニング（ZRP）では不十分である可能性を示唆している。著者らは、拡大する宿主範囲や特定の生態学的ドライバー（例：カモツル科、特定の水鳥、沿岸湿地）を考慮した、更新されたリスクゾーンとサーベイランス戦略が必要であると主張している。
• 手法の進展： 連続的な系統地理学と高解像度の生態学的データを統合することで、疫学的データだけでは見逃される可能性のある、祖先ウイルスの移動とそれを駆動する要因を再構築する枠組みを提供している。

著者らは、ウイルスの拡散を予測するには、この不均一性を認識する必要があると結論付けており、サーベイランス戦略は、進化する疫学的パターンや影響を受ける宿主の範囲の拡大を反映して、頻繁に更新されなければならないとしている。