Assembly of the Peganum harmala microbiome across plant compartments and growing seasons

技術要約：植物の各コンパートメントおよび生育シーズン間における Peganum harmala マイクロバイオームのアセンブリ

問題提起
北アフリカ、中東、中央アジアの乾燥および半乾燥地域に自生する多年生薬用植物である Peganum harmala L. は、$\beta$-カルボリンアルカロイドを生成することから、重要な民族薬理学的価値を有している。本種は過酷な環境条件（塩分、アルカリ性、栄養不足）の下でも生育することで知られているが、それに随伴するマイクロバイオームの形成と構造を制御するメカニズムについては、依然として十分に解明されていない。これまでの研究は、主に個別の植物コンパートメントや特定の微生物界（細菌または真菌）に焦点を当てた培養依存的な手法を用いてきた。根圏土壌から根、そして葉の内生エンドスフェアへと続く連続体の中で、細菌および真菌コミュニティがいかにして同時に動員、フィルタリング、および構造化されているのか、また、乾燥した生態系においてこれらのコミュニティが季節を通じてどのように変化するのかを理解することには、決定的な空白が存在する。

方法論
本研究では、モロッコのベンゲリール地域における野生の P. harmala 個体群のマイクロバイオームを、2つの連続する生育シーズン（2022年6月および2023年6月）にわたって包括的に特徴付けるため、培養に依存しないアプローチを採用した。

• サンプリング: 研究者は、各シーズンにおいて根圏土壌、根、および表面殺菌した葉を、32個体ごとに採取した。
• 物理化学的分析: 土壌サンプルは、テクスチャー、pH、電気伝導度、有機物、および主要栄養素（P, K, Ca, Mg, Na, N）ならびに微量栄養素（Cu, Fe, Mn, Zn）について分析された。
• シーケンシング: 全DNAを抽出し、ハイスループット・アンプリコンシーケンシングに供した。細菌コミュニティは16S rRNA遺伝子シーケンシング（プライマー799F/1115R）を用いてプロファイリングし、真菌コミュニティはITS2領域を用いてプロファイリングした。
• バイオインフォマティクスおよび統計: 生データはDADA2パイプラインで処理され、アンプリコン配列バリアント（ASV）を生成した。分類学的割り当てには、4つのデータベース（SILVA, Greengenes2, RDP, GTDB）を用いた階層的優先順位を利用した。統計解析には以下を用いた：
  • アルファおよびベータ多様性の評価（Chao1, Shannon, Simpson, UniFrac, Bray-Curtis）。
  • コミュニティ組成の駆動要因を特定するためのPERMANOVA。
  • コンパートメント特異的なバイオマーカーを特定するためのLEfSeおよび指標種解析。
  • 微生物間の相互作用を探求するための共起ネットワーク解析（SparCC）。
  • 微生物のフローを追跡するためのソーストラッキング（サンキーダイアグラム）。
  • PICRUSt2（細菌）およびFUNGuild（真菌）による機能予測。

主な結果

1. 階層的なコンパートメント化:
   マイクロバイオームは明確な階層構造を示した。根圏土壌は最も高い微生物多様性を保持しており、主要なリザーバー（貯蔵庫）として機能していた。植物のコンパートメント（根および葉）は、微生物のより小さなサブセットを段階的に選択した。コンパートメントの同一性は、細菌（PERMANOVA $R^2 = 0.204$）および真菌（$R^2 = 0.187$）の両方において、季節的影響を大きく上回る、コミュニティ組成の主要な駆動要因であった。

2. 分岐する多様性パターン:

   • 細菌: 多様性は、両年ともに根圏から根、そして葉のエンドスフェアへと一貫して減少した。しかし、エンドスフェアの細菌多様性は2年目に大幅に増加し（根で+75%、葉で+43%）、エンドファイトの漸進的な蓄積を示唆した。
   • 真菌: 真菌の多様性パターンは特異であった。初年度は、葉のエンドスフェアは単一の属である Leveillula（約90%）によって支配されていた。2年目には、葉の真菌多様性が著しく増加し、Alternaria、Fusarium、Aspergillus を含む、より多様な集合体へとシフトした。

3. 時間的ダイナミクスとコア・マイクロバイオーム:
   顕著な年間の分類学的ターンオーバーにもかかわらず、安定した細菌コア・マイクロバイオームが季節を超えて持続していた。それは Pseudonocardiaceae、Bacillaceae、および Sphingomonadaceae のメンバーで構成されていた。ソーストラッキングの結果、環境からの動員は初年度の約50%から2年目の約64%へと増加しており、マイクロバイオームが垂直伝播のみによって維持されるのではなく、土壌プールから動的に再構築されていることを示している。

4. バイオマーカーとネットワーク・トポロジー:

   • 細菌: 葉のエンドスフェアは、Tepidiforma、Kallotenue、Arsenophonus といったストレス耐性を持つタクサ（分類群）に富んでいた。また、安定したシアノバクテリアの指標種として、両年ともに葉のエンドスフェアにおける Oscillatoria が見出された。
   • 真菌: 初年度は葉における極端な専門化（Leveillula）を特徴としていたが、2年目はより広範な真菌の多様化を示した。
   • ネットワーク: 細菌ネットワークは緻密でモジュール的であり、「コネクター」ハブ（例：Patescibacteria、Verrucomicrobiota）がモジュール間を連結していた。真菌ネットワークはより単純で結合性が低く、時間的変化に対して高い感受性を示した。

5. 環境駆動要因:
   土壌の物理化学的特性（例：粘土含有量、CEC、pH、CaCO$_3$）は、特定の微生物タクサと相関していた。例えば、Microvirida と Truepera は粘土および湿度と相関し、Fusarium と Penicillium はCaCO$_3$ およびリンと相関していた。しかし、植物のコンパートメントは土壌特性よりも強力なフィルターとして機能していた。

意義および主張
著者らは、本研究が、複数のコンパートメントおよび時点にわたって細菌および真菌コミュニティを統合した、P. harmala マイクロバイオームの最初の包括的な評価を提供すると主張している。主な貢献は以下の通りである：

• メカニズムへの洞察: 本研究の知見は、P. harmala のマイクロバイオームのアセンブリが、メタコミュニティの枠組みの中で機能する、強力なコンパートメント・フィルタリングと動的な環境動員によって駆動されていることを明らかにしている。
• 界（Kingdom）特異的なアセンブリ: 本研究は、対照的なアセンブリ・ルールを浮き彫りにした。すなわち、細菌コミュニティはより決定論的にフィルタリングされ安定しているのに対し、真菌コミュニティはより高い時間的可塑性を示し、環境分散に対してより敏感である。
• 機能的安定性: 年間の間で大幅な分類学的ターンオーバーにもかかわらず、予測された機能プロファイルは安定していた。これは、異なるタクサが異なる生態学的役割（例：ストレス耐性、栄養循環）を果たす機能的冗長性が存在することを示唆している。
• バイオテクノロジーへの可能性: 乾燥適応したホストの中に、ストレス耐性を持つ微生物（窒素固定能を持つ可能性のある Bradyrhizobium や Mesorhizobium を含む）が持続していることは、P. harmala が極限環境に適応した微生物の貴重なリザーバーであることを位置づけている。

限界と謙虚さ
論文は、その結論に関して控えめなトーンを維持している。著者らは以下のことを明示している：

• 機能プロファイルはアンプリコンデータ（PICRUSt2/FUNGuild）に基づく予測であり、測定された代謝活性を表すものではない。
• 共起ネットワークは統計的な関連性を示すものであり、直接的な生物学的相互作用を示すものではない。
• P. harmala の二次代謝物（アルカロイド）がマイクロバイオームの形成に果たす具体的な役割は未解決であり、さらなる実験的検証を必要とする。
• 特定されたタクサの生態学的機能（非マメ科における窒素固定など）は予備的なものであり、それらを確定させるためには分離および接種実験が必要である。

要約すると、本研究は、P. harmala が、乾燥条件に適応するために、安定したコア機能と多様な環境タクサの動員とのバランスを取りながら、構造化されつつも動的なマイクロバイオームを維持する、選択的なホストであることを確立している。