Sequence-equivalent bacterial ASVs occupy dynamic network roles across amphibian development

技術要約：配列が等価な細菌ASVは両生類の発生過程において動的なネットワーク上の役割を占める

問題提起
宿主に関連する微生物叢は、環境的な微生物ソースと宿主特異的な生態学的フィルタリングの相互作用によって形成される。これらのコミュニティの組成については十分に研究されているが、その構造的組織化を制御するメカニズムについては、未だ理解が進んでいない。特に、特定の発生段階（個体発生）を超えて存続する特定の細菌系統が、どのように微生物ネットワークのトポロジーの再編成に関連しているのかは不明である。先行研究により、両生類の皮膚微生物叢は種、段階、および環境によって異なることが確立されているが、特定の細菌の同一性の存続を、コミュニティネットワーク構造内での変化する役割と結びつけた研究はほとんどない。

方法論
本研究では、分類学的アプローチとネットワーク的アプローチを統合した多段階の解析フレームワークを用い、アンデス産の両生類 Rhinella spinulosa の皮膚微生物叢を特徴付けた。

• データの取得と統合： 著者らは、変態前（Gosnerステージ ≤25）のオタマジャクシおよびその周囲の水柱環境から、新たな16S rRNA（V3–V4領域）シーケンシングデータを生成した。これらのデータは、既報の幼体および成体の R. spinulosa のデータセットと統合された。
• 標準化された処理： 異なる地域からサンプリングされた各段階間の方法論的一貫性を確保するため、すべてのデータセット（オタマジャクシ、幼体、成体）を統一されたQIIME 2ワークフロー（v2024.10）を用いて再処理した。これには、プライマー除去、Deblurによるデノイジング、および199 bpへのトリミングが含まれ、アンプリコン配列バリアント（ASV）を生成した。
• 生態学的分類： オタマジャクシの皮膚ASVを、水柱に対する皮膚での出現に基づき、water-only（水のみ）、shared (non-core)（共有・非コア）、skin-only (non-core)（皮膚のみ・非コア）、core (shared)（コア・共有）、および core (skin-restricted)（コア・皮膚限定）のカテゴリーに分類した。コア微生物叢は、≥80%の個体で存在するASVと定義した。
• ネットワーク推論： 微生物関連ネットワークを、SPIEC-EASI（Sparse Inverse Covariance Estimation for Ecological Association Inference）のMeinshausen-Bühlmann近傍選択法を用い、オタマジャクシ、幼体、および成体の各段階について個別に推論した。
• トポロジー解析： グローバルなネットワーク特性（連結性、クラスタリング、不均一性、中心性）およびノードレベルの指標（次数、媒介中心性）を算出した。ハブは、次数の分布または媒介中心性の分布の上位5%にあるノードとして定義された。
• 存続解析： 配列レベルのマッチングにより、3つの段階すべてで検出されたASV（TJA）、オタマジャクシと幼体のみのASV（TJ）、およびオタマジャクシのみのASV（T）を特定した。本研究では、これらの存続する系統の存在量とトポロジー上の役割（ハブの状態）を、発生過程を通じて追跡した。
• 堅牢性： トポロジーパターンの安定性を評価するために、ブートストラップ・リサンプリングおよび幼体サンプルのランダム・サブサンプリングによる感度分析を実施した。

主な結果

1. 宿主によるフィルタリングと分画： オタマクシの皮膚微生物叢は、水中の微生物プールから構造化されたサブセットを形成していた。多くのASVが水と共有されていたが、かなりの部分が皮膚限定的であった。コアASV（≥80%のオタマジャクシに存在する）は、平均相対存在量および出現率が高いことを示した。決定的なことに、一貫した出現（出現頻度）は、水と皮膚の間での存在量の差に厳密には依存しておらず、これらが補完的なアセンブリの次元であることを示唆している。
2. ネットワーク構造のダイナミクス： ネットワークトポロジーは、発生段階によって有意に変化した：
   • 幼体： 最も高い連結性（ノードあたり22.20の隣接ノード）、クラスタリング（0.470）、不均一性、および中心性を示した。
   • 成体： ネットワークは縮小し（幼体の1,196ノードに対し308ノード）、7つの成分に断片化した。しかし、高いクラスタリングと中心性は維持していた。
   • オタマジャクシ： 中間的な連結性（7.27の隣接ノード）と、単一の連結成分を示した。
   • 存在量とトポロジーの関係： オタマジャクシにおいて、存在量は次数とは無関係であったが、媒介中心性と正の相関があった。幼体および成体においては、存在量は次数と負の相関を示すようになったが、媒介中心性とは引き続き正の相関を示した。
3. コア微生物叢の役割： オタマジャクシにおいて、コアASVは次数の定義によるハブに不釣り合いに多く含まれており、非コアASVよりも高い媒介中心性を示した。これは、一貫して関連する分類群が、構造的に顕著な地位を占めていることを示している。
4. 存続とネットワーク上の役割： 配列レベルのマッチングにより、3つの段階すべてに存在する44のASV（TJA）を特定した。
   • 存在量： TJA ASVは、TJまたはオタマジャクシのみのASVよりも有意に高頻度であった。
   • トポロジーの軌跡： 主要な解析において、TJA ASVはTJ ASVと比較して、次数は低いが媒介中心性は高いという特性を示した。しかし、感度分析により、媒介中心性のパターン（TJAの方が高い）は確認されたものの、次数のパターンを一貫して再現することには失敗した。
   • 役割の流動性： ほとんどの存続するASV（31/44）は、発生過程を通じて非ハブのままであった。他のASVは、特定の段階でハブの状態（次数または媒介中心性）を獲得、あるいは喪失した。したがって、細菌系統の存続は、そのネットワーク上の役割の存続を意味しない。

意義と主張
本論文は、宿主によるフィルタリング、段階に関連したネットワークの再編成、および系統の存続が、両生類の皮膚微生物叢のアセンブリにおける相補的な次元であることを主張している。本研究は、分類学的な連続性（特定の細菌系統の存続）が、コミュニティ構造およびネットワークトポロジーの劇的な変化と共存し得ることを実証している。

著者らは、初期の宿主によるフィルタリングが皮膚関連微生物のサブセットを確立する一方で、これらの系統は宿主の発達に伴い、そのネットワーク上の役割において動的な再編成が行われると結論付けている。具体的には、成体まで存続する系統は、高い媒介中心性（ネットワーク領域を橋渡しする役割）を持ちつつも、低い局所的な連結性（次数）を持つ傾向があり、このパターンはサンプリングの変動に対して堅牢であると考えられる。本研究は、一貫した出現が静的な生態学的役割を意味するわけではないこと、すなわち、微生物系統は存続しながらも、段階特異的なネットワーク内でのトポロジー上の位置を変化させ得ることを強調している。

著者が認める限界事項
著者らは、発生段階が地理的起源（各段階で異なる地域）と混同されているため、個体発生の効果と集団レベルまたは環境レベルの効果を分離できないことを指摘している。また、サンプルサイズが小さいこと、およびネットワーク推論の性質（直接的な生態学的相互作用ではなく統計的な関連性であること）が、グローバルな指標の解釈を制約している。存続は199-bpフラグメントの正確な配列一致によって定義されており、これは連続的な定着と、環境からの配列等価な系統の繰り返しの補充とを区別できない。