技術要約：健康時およびIBDにおけるPhocaeicola vulgatusの温和性および糸状ファージにおける明確な進化的選択パターン

問題提起

腸内ウイルス叢は炎症性腸疾患（IBD）において著しく変化しており、溶菌性ファージから温和性ファージへのシフトを特徴とする。しかし、これらの変化を駆動する具体的な進化メカニズムは依然として十分に理解されていない。機能的な理解における大きな障壁は、腸内ファージの大部分について既知の細菌宿主が欠如していることであり、これが炎症の文脈におけるファージと細菌の相互作用および進化の研究を制限している。標準的なメタゲノム解析手法では、低存在量のファージを捉えたり、特定の宿主集団内における一塩基多型（SNP）や構造的再編成といった微細な集団動態を解明したりすることが困難である場合が多い。

方法論

本研究では、特定の宿主であるPhocaeicola vulgatus（旧名Bacteroides vulgatus。重要なパソビオン）と物理的に相互作用するファージを単離・特性評価するために、**ウイルス・タギング（Viral Tagging: VT）**を用いた。

• サンプル収集: クローン病（CD）患者および、それらと一致する健康な同居家族コントロール（HHC）の糞便サンプルから、ウイルス様粒子（VLP）をプールした。
• 実験デザイン:
  • ウイルス・タギング: VLPプールを、P. vulgatus株（$\Delta$BV01、天然のBV01プロファージを欠損させたもの）とインキュベートし、フローサイトメトリーを用いて蛍光（細菌用SYTO 9、VLP用SYTO 62）に基づきソートした。バルクソートされた細胞に対して、複数置換増幅法（MDA）およびシーケンシングを行った。
  • in vitro濃縮: 並行して、単一のHHCサンプル由来のVLPをP. vulgatus培養物に通すことで、活発に複製しているファージを選択するアプローチを用いた。
  • メタバイローム比較: 結果を、入力されたメタバイロームおよび外部データセット（169の公開されたCD、潰瘍性大腸炎（UC）、およびHHCのバイローム）と比較した。
• バイオインフォマティクス:
  • リードのフィルタリング、アセンブリ（MEGAHIT, metaSPAdes）、および宿主除去のためにカスタムパイプライン（PvVTPhageFinder, VLPVirifer）を使用した。
  • ウイルスコンティグを、95%の平均塩基同一性（ANI）でウイルス操作分類単位（vOTU）へとクラスター化した。
  • 進化解析として、選択圧を評価するためにSNP呼び出し（LoFreq）および$P_n/P_s$（非同義/同義）比の計算を行った。また、リードマッピングとモチーフ特定を通じて、構造的再編成（特に「シャッフロン」組換えイベント）を解析した。
  • ライフサイクル動態として、不一致リードペアおよびスプリットリードを用いて、プロファージの統合およびゲノム環状化率を定量化した。

主な貢献

1. 宿主定義によるファージ回収: 本研究は、標準的なメタバイロームのアセンブリではほとんど捕捉できなかった、P. vulgatusに特異的な43のvOTUを回収することに成功した。これは、VTが超低存在量のウイルス多様性を捉える能力を持つことを示している。
2. 新規ファージの発見:
   • リファレンスファージBV01に近縁な温和性尾部ファージ（vOTU_0）を特定した。
   • PvIno1、すなわちBacteroidaceae科の宿主として初めて記述される糸状ファージ（Inoviridae科）を発見した。
3. 高解像度の進化的動態: 本研究は、コンセンサス配列を超えて集団内の微小進化を解析し、健康状態と疾患状態の間で異なる選択パターンを明らかにした。

結果

1. 方法論的比較と多様性

VTおよびin vitro濃縮は、標準的なメタバイロームには存在しないユニークなファージ集団を捉えた。メタバイロームはCaudoviricetes、Crassvirales、およびMicroviridaeによって支配されていたが、VTは糸状ファージ（Inoviridae）を強く濃縮しており、これらはメタバイロームではほとんど観察されなかった。クロスマッピングの結果、条件間（CD vs HHC）の明白な排他性は、生物学的な不在ではなく、多くの場合、低存在量によって引き起こされていることが判明したが、VTは他の手法では検出不可能な疾患関連の多様性を識別することに成功した。

2. BV01様ファージ（vOTU_0）の進化

• 選択パターン: BV01様ファージは、CDサンプルと比較してHHCサンプルにおいて有意に高い塩基多様性を示した。
  • 健康状態（HHC）: 尾部構造遺伝子（テイルタンパク質1、2、3、およびテープメジャータンパク質）において**多様化選択（diversifying selection）**を示しており、これは宿主との活発な共進化の軍拡競争を示唆している。特定の置換（テイルタンパク質2におけるPro14Ser）がHHC集団において固定されていた。
  • 疾患状態（CD）: ゲノム全体にわたって**浄化選択（purifying selection/保存的選択）**を示し、健康時で見られた多様性が消失していた。
• シャッフロン組換え: リコンビナーゼに挟まれたマルチカセット・シャッフロンを含む1.7 kbのゲ genomic領域が特定された。このシャッフロンは、カセットの反転を介して6つの異なるテイルバリアントを生成する。
  • 6つすべての構成がCDとHHCの両方で検出されたが、その相対的な頻度は条件間で大きく異なっていた。
  • in vitro濃縮では、構成1（天然型）と構成6が優占しており、特定のテイルタイプがin vitroでの維持に有利であることを示唆している。
• ライフサイクル活性: BV01はこれまで培養不可能であったが、VTサンプルでは既知のBV01統合部位における高い統合率が示された。in vitro濃致では、活発な感染、高い環状化率、およびライソジェニシスが確認された。

3. 糸状ファージPvIno1（vOTU_10）の進化

• 疾患との関連: PvIno1は100%のCD VTサンプルで検出され、これらのみから由来しており、HHC背景に対して>58,000倍の濃縮を示した。これは、潜在的な疾患関連バイオマーカーであることを示唆している。
• ゲノムの保存性: このファージは顕著な配列保存性（ANI > 99.85%）を示し、13〜14個の遺伝子をコードするコンパクトなゲノム（~5,994 bp）を有している。
• 複製動態: 正規化された深度プロファイルは、複製起点（$ori+$および$ori-$）に対応する2つの主要な枯渇部位を示した。線状コンカテマー中間体の存在と高い環状化率（42.9%）は、宿主細胞内での活発なローリングサークル複製（RCR）を示している。
• 変異パターン: BV01様ファージとは異なり、PvIno1は非常に低い塩基多様性を示した。特定された数少ないSNPは、構造遺伝子ではなく、予測される$ori-$周辺の調節領域に集中していた。これは、糸状ファージにおける進化圧が、テイルファージで見られるような構造遺伝子の多様化ではなく、調節遺伝子を標的にしていることを示唆している。

重要性と主張

本論文は、炎症が同じ宿主に感染する異なるタイプのファージに対して、ユニークかつ明確な進化的圧力を課していると主張している。

• 温末性尾部ファージ（BV01様）は、宿主の表面変異（CPS）を回避するために、健康時には構造的なテイル遺伝子において多様化選択を受け、炎症時にはその多様性が浄化選択によって崩壊する。
• 糸状ファージ（PvIno1）は根本的に異なる軌跡を辿り、多様化は調節領域に集中しており、疾患状態（CD）との強い関連性を示す。

ウイルス進化をIBDファージオームのエコロジカルモデルに統合することで、本研究はファージベースの治療薬設計のための基礎を確立している。また、ファージの活動（SNP、シャッフロン、調節変異などの微細な適応）を予測するためには、これらを理解することが不可ло欠であることを強調している。本研究は、低存在量の宿主特異的なウイルス動態を研究する上で、標準的なメタゲノミクスの限界を克服するために、VTのような宿主標的型の手法が必要であることを裏付けている。