FlavoTyper: a genome-based in-silico serotyping tool for the fish pathogen Flavobacterium psychrophilum

技術要約：FlavoTyper – Flavobacterium psychrophilum のゲノムに基づくイン・シリコ血清型決定ツール

問題提起
Flavobacterium psychrophilum は、淡水魚に甚大な被害をもたらす病原体であり、世界中の養殖業において大きな経済的損失を引き起こし、野生個体群にも影響を及ぼしている。効果的な疫学調査やワクチン開発は、リポ多糖（LPS）のO-多糖（O-PS）部位の抗原多様性に依存しているが、従来の血清学的手法はコストが高く、労力を要し、動物由来の抗血清を必要とするほか、解釈の主観性や研究室間での標準化の欠如といった課題がある。wzy 遺伝子を標的としたマルチプレックスPCR（mPCR）スキームによって標準化は改善されたものの、これは限定的なゲノム変異（5つの型）しか捉えることができず、R基（R-group）のバリエーションを含むO-抗原（O-Ag）遺伝子座の全容を解明するには不十分である。全ゲノムシーケンシング（WGS）データの蓄積に伴い、O-Ag 構造の多様性の全スペクトルを解明できる、堅牢で再現性のある包括的なイン・シリコ血清型決定ツールの必要性が高まっている。

方法論
本研究では、ゲノムアセンブリから直接血清型を予測するために設計された、オープンソースの Python ベースのバイオインフォマティクス・パイプラインである FlavoTyper を開発した。その手法は、主に以下の3つのフェーズで構成される：

1. バイオマーカーの発見とデータベースの構築：

   • 公開ゲノムアセンブリ（NCBI および PubMLST）の O-Ag 生合成遺伝子座に対して、広範な比較ゲノミクス解析を実施した。
   • O-Ag 遺伝子座を解析し、血清型を定義する3つの可変領域、すなわち O-型（wzy または wzy 様遺伝子によって決定）、S-型（特定の4遺伝子クラスターによって決定）、および R-型（R基の生合成/転送に関与する遺伝子によって決定）を特定した。
   • 配列保存性が低い wzy 変異体を分類するため、著者らは AlphaFold 3 を用いた 3D 構造予測、および DALI/Foldseek による構造類似性検索を利用した。これにより、配列同一性は低いものの、構造的な保存性によって特定の wzy タンパク質（例：Wzy0, Wzy4, Wzy6, Wzy7）が異なるファミリーにグループ化されることが明らかになり、これらが独立した血清型決定因子として分類されることが検証された。
   • O-型（O:0 ～ O:7）、S-型（S0, S1）、および R-型（R0 ～ R4、R1 サブバリアント V1, V2, V3 を含む）の参照配列を含むバイオマーカーデータベースを構築した。

2. ソフトウェアの実装：

   • FlavoTyper は、FASTA 形式のゲノムアセンブリを受け取るモジュール式パイプラインである。以下の処理を行う：
     • 品質管理（QC）： 型株 F. psychrophilum NCIMB 1947T に対する FastANI を用いた種バリデーション（閾値 ≥95% ANI）。
     • タイピング： バイオマーカーデータベースに対するクエリゲノムの BLASTn アライメント。
     • 遺伝子座解析（オプション）： 検出されたマーカーを可視化するための O-Ag 遺伝子座の抽出およびマッピング。
   • 真のマーカーと部分的なヒットを区別するために、経験的に導出された配列同一性 ≥97% かつカバレッジ ≥94% というタイピング規則を確立した。

3. 検証と解析：

   • 本ツールを NCBI の 482 個のアセンブリおよび PubMLST の厳選された 357 個のゲノムに適用した。
   • 334 個の解明された株に対してコアゲノム系統樹を再構築し、血清型とクローナルコンプレックスの関係を分析した。
   • ホスト魚種と血清型の統計的な関連性を、フィッシャーの正確確率検定を用いて分析した。

主な結果

• 多様性の解明： FlavoTyper は、公開データセット全体において 23 の異なる血清型（O:x–S:y–R:z の組み合わせ）を特定することに成功し、従来の mPCR スキームの5つの型と比較して、解像度を大幅に向上させた。最も支配的な2つの血清型は O:2–S0–R1V1 および O:1–S0–R1V1 であった。
• 性能： 本ツールは、NCBI データセットの 77.8%、および厳選された PubMLST データセットの 94% に対して血清型を解明した。また、非 F. psychrophilum のアセンブリ（NCBI データの 15.6%）を効果的に排除し、断片化したゲノムや混合培養の問題をフラグ立てした。
• ホストとの関連性： イン・シリコ・タイピングにより、既知のホスト・血清型関連性が確認および精緻化された：
  • アユ（Plecoglossus altivelis）は O:3 と強く関連していた。
  • コホーサーモン（Oncorhynchus kisutch）は O:0 および O:7 が豊富であった。
  • ニジマス（Oncorhynchus mykiss）は O:1 および O:2 が豊富であった。
• 組換えと進化： 系統解析の結果、コアゲノム系統と O-Ag 遺伝子座との間に顕著な不一致が認められた。O-型は様々なクローナルコンプレックスに分散しており（例：CC-ST10 には 5 つの異なる O-型が含まれる）、O-Ag 遺伝子座が組換えホットスポットとして機能していることを示している。これは、血清型の多様性が垂直継承ではなく、頻繁な水平遺伝子交換によって駆動されていることを示唆しており、O-抗原がコアゲノムのバックボーンから切り離されていることを意味する。

意義と主張
著者らは、FlavoTyper を F. psychrophilum に特化した初のゲノムベースのイン・シリコ血清型決定ツールとして位置づけている。その意義は以下の通りである：

• 標準化： 従来の血清学における主観性とリソース集約的なプロセスを排除する、再現可能で自動化された枠組みを提供する。
• 包括的な解像度： 従来の mPCR では解明できなかった R 基および S 基のバリエーションを含む、O-抗原の全ゲノム変異を捉える。
• 疫学的有用性： 確立されたホスト特異的な傾向を再現するツールの能力を検証しており、監視、アウトブレイクの検出、および非サケ科魚類（ウナギやコイなど）における病原体の出現追跡への利用を支持している。
• アクセシビリティ： 本ツールは PyPI、Bioconda、Galaxy Europe を通じて配布されており、また PubMLST BIGSdb のプラグインとしても提供されているため、日常的な監視やワクチン株の選択を行う研究者や獣医専門家にとって容易に利用可能である。

本論文は、O-Ag 遺伝子座は頻繁にシャッフルされるものの、本ツールは全ゲノムシーケンシング時代における、この種のホスト・病原体相互作用および進化ダイナミクスを理解するための必要な枠組みを提供するものであると結論付けている。