Co-evolutionary dynamics of defense systems and anti-defense genes in chromosomes and mobile genetic elements of Enterobacteriaceae

技術的要約：腸内細菌科における防御システムと抗防御遺伝子の共進化ダイナミクス

問題提起
細菌の防御システム（DEF）と可動遺伝因子（MGE）は、複雑な進化の軍拡競争を展開している。DEFがしばしばMGE（プラスミド、プロファージ、統合共役要素など）内にコードされていることは確立されているが、これら防御システムの進化的なターンオーバーに対するゲノムコンパートメント（染色体 vs MGE）の具体的な寄与については依然として不明である。さらに、抗防御遺伝子（anti-DEF：宿主の免疫を抑制するもの）の定量的進化ダイナミクスと、その標的となるDEFとの関係、およびそれらの軌跡が連動しているかどうかについても、大規模な系統ゲノムデータセットを用いた体系的な特徴付けがなされていない。本研究は、Enterobacteriaceae（腸内細菌科）におけるDEFおよびanti-DEFの遺伝子獲得、喪失、拡張、および減少率に関する、異なるゲノム背景間での定量的な比較の欠如に対処するものである。

方法論
著者らは、Enterobacteriaceaeの213種から得られた6,907の高精度ゲノムを解析した。ワークフローは以下の通りである：

1. ゲノム処理およびアノテーション： ゲノムは完全性（>99%）およびコンタミネーション（<1%）に基づいてフィルタリングされた。遺伝子の予測およびアノテーションにはpyrodigal-gvおよびgeNomadを用いた。MGE（プロファージ、プラスミド、ICE/IME、およびサテライト）は、特定のツール（geNomad、Platon、SatelliteFinder、ICEfinder-opt）を用いて同定され、染色体由来の遺伝子とMGE関連の遺伝子を区別するための厳格なクロスバリデーション基準が適用された。
2. 防御システムの同定： DEFおよびanti-DEFはDefenseFinderを用いて同定された。オーソログ遺伝子ファミリー（OG）はOrthoFinderを用いて構築された。
3. 進化速度モデリング： Countに実装された系統的出生・死滅モデルを用い、4つの進化プロセス（遺伝子獲得、遺伝子喪失、拡張（コピー数増加）、および減少（コピー数減少））を定量化した。これは、ゲノム上の位置（染色体、プロファージ、プラスミド/ICE/IME）によって分類されたOGに対して適用された。
4. 統計解析：
   • ターンオーバーの比較： ガンマ分布を用いた一般化線形混合効果モデル（GLMM）を適合させ、DEFと非防御（ND）遺伝子、および異なるゲノムコンパートメント間の速度を比較した。
   • 結合解析： 分岐長を制御し、系統的ブラウン運動シミュレーションを通じて評価された偏スペアマン相関を用い、DEFとanti-DEFの獲得と喪失が系統間で進化的に結合しているかどうかを判定した。
   • 層別化： 解析は全ゲノムレベル、および特定のゲノムコンパートメントや防御システムサブタイプ（例：CRISPR-Cas、RM、CBASS）ごとに層別化して行われた。

主な結果

1. ゲノム背景がターンオーバーのバイアスを駆動する： DEFのターンオーバーは、それらを保持するゲノム要素の背景ダイナミクスとは有意に異なる。

   • 染色体： 染色体上のDEF（cDEF）は、染色体上の非防御遺伝子と比較して、最も高いターンオーバー率（獲得、喪失、拡張、および減少）を示した。
   • MGE： 対照的に、プロファージ上に位置するDEFは最も低いターンオーバー率を示し、周囲のプロファージ遺伝子プールよりも高い安定性と、少ないコピー数変動を示した。プラスミドおよびICE/IMEに関連するDEFは、中間的なターン需要を示した。
   • バランス： 絶対的な速度の違いにもかかわらず、獲得/喪失および拡張/減少の相対的なバランスは、染色体とMGEの間で同様であった。

2. 抗防御ダイナミクス： anti-DEFは、一般に標的となるDEFよりも高い喪失、拡張、および減少率を示し、特にMGE上に位置する場合にその傾向が顕著であった。

   • プロファージおよびプラスミド/ICE/IME上において、anti-DEFはDEFよりも有意に高い喪失率とコピー数可塑性を示した。
   • 染色体上では、anti-DEFはDEFよりも獲得率が低かったが、喪失およびコピー数変化の速度は同程度であった。
   • 層別解析により、特定のシステム（例：anti-CBASSおよびanti-RM）において、標的よりも高い喪失およびコピー数変化率を示すという傾向が確認された。

3. 進化的結合： DEFとanti-DEFの獲得と喪失は、細菌の系統間で正の相関がある。

   • この結合は、全ゲノムレベルおよび特定のMGEコンパートメント（プロファージおよびプラスミド）内で観察された。
   • 正の関連性は、「一致する（拮抗的な）」ペア（例：RM Type I–ArdA/ArdU、CRISPR-Cas I-E–AcrIE9）および「一致しない」ペアの両方で見られた。
   • MGE関連のDEFおよびanti-DEFのターンオーバーは、より広範なMGE含有量の変化に伴う機能的遺伝子カテゴリー（COG）のターンオーバーとも相関しており、これは彼らのダイナミクスがしばしばMGEコンテンツのより大規模な変化の一部であることを示唆している。

意義および主張
本論文は、ゲノム背景（染色体 vs MGE）が防御システムの流動性の主要な決定要因であると主張している。著者らは、染色体上のDEFの高いターンオーバーは、細菌が染色体上の防御を失い、MGE由来の防御へと優先的に移行するという「防御のアウトソーシング」戦略を反映しているか、あるいは染色体上の防御アイランドが、ラテラル・トランスダクションや可動要素への近接性により、本質的に可塑性が高いことを反映しているのではないかと提案している。

本研究は、DEFとanti-DEFの進化的な軌跡が結合しているという仮説に対して実証的な支持を与えている。この結合は、直接的な軍拡競争ダイナミクス（一致するペア間の生態学的結合）だけでなく、一致しない防御およびカウンター防御システムを運ぶMGEのより広範な背景的ターンオーバーによっても駆動されている。これらの知見は、プロファージ関連のDEFの保持が選択的な維持またはドメスティケーション（家畜化）の結果である可能性を示す一方で、染色体上のDEFの急速な喪失は、安定したゲノム環境における免疫維持のコスト・ベネフィットのトレードオフを浮き彫りにしている。

著者らは、ゲノム数が少ない種に対する統計的パワーが限定的であることや、解析がEnterobacteriaceaeに限定されていることなど、結論に対して慎重な姿勢を保っている。また、一部のMGE関連DEFの機能的状態（例：それが完全なシステムなのか、あるいは縮小した派生物なのか）については、さらなる直接的な解析が必要であることも認めている。