Integrated metabolomic and genomic characterization of antifungal Streptomyces halstedii S14 from Melipona mondury larval food

技術要約：Streptomyces halstedii S14の統合的なメタボロミクスおよびゲノム特性解析

問題提起
ハリナシバチ属（Melipona spp.）の幼虫食には、天然物探索の可能性を秘めた微生物コミュニティが含まれていることが知られているが、特にMelipona monduryの幼虫食に関連する放線菌の化学的およびゲノム的な潜在能力については、依然として十分に解明されていない。これまでの研究は、本種の腸内細菌叢や蜂蜜のプロファイルに焦点を当てており、このニッチにおける培養可能な放線菌の生合成能に関する理解には空白が存在している。さらに、Streptomycesと真菌（Trichodermaなど）との共培養による相互作用が代謝変化を誘発することは知られているが、これらM. mondury由来の分離株における化学的シグネチャーの空間的分布と、ゲノム潜在能との直接的な相関関係は、体系的に特徴付けられていない。

方法論
本研究では、空間メタボロミクス、分子ネットワーク解析、およびゲノムマイニングを組み合わせた統合的な「メタボロゲノミクス」アプローチを用い、Streptomyces halstedii S14（M. monduryの幼虫食から単離）と糸状菌Trichoderma atrovirideとの間の拮抗作用を調査した。

1. 単離および特性評価： 放線菌を、シクロヘキシミドおよびニスタチンを添加したISP-9、NYDB、およびISP-2培地を用いて濃縮した。菌株S14を単離し、形態学的特性（肉眼的観察、グラム染色、SEM）を評価した後、全ゲノムシーケンシング（Illumina NextSeq 550）および平均塩基配列同一性（ANI）解析により分類学的な同定を行った。
2. 拮抗試験： NYDAプレートを用い、中央にS14を、周辺部にT. atrovirideを接種する共培養アッセイを実施した。阻止円（IZ）の形成を測定した。
3. 空間メタボロミクス： 共培養の5つの異なる領域（放線菌遠位：AD、放線菌近位：AP、阻止円：IZ、T. atroviride近位：TP、T. atroviride遠位：TD）から、アセプティックな手法で寒天片を回収した。コントロールとして、単独培養および培地のみのプレートを用いた。
4. UHPLC-HRMS/MS解析： Nexera X2 UHPLCシステムとLCMS-9030 qTOF質量分析計を組み合わせて抽出物を分析した。データ処理には特徴検出のためにMZmine 4を用い、続いて統計解析において、生物学的変動を考慮するためにプレートブロック線形モデル（limma）を用いたMetaboAnalystを使用した。
5. 分子ネットワーク解析： 特徴ベースの分子ネットワーク（FBMN）をGNPS2プラットフォーム上で実施し、分子ファミリーを整理した後、スペクトルライブラリ照合（GNPS/MassBank）およびイン・シリコ予測（SIRIUS/CANOPUS）による推定注釈を行った。
6. ゲノムマイニング： S14のゲノムをantiSMASH 8.0.4を用いて解析し、生合成遺伝子クラスター（BGC）を特定した。MIBiG 3.1データベースとの比較解析をBiG-SCAPE 2.0.3を用いて行い、BGCの類似性と遺伝子クラスターファミリー（GCF）を決定した。シンテニー解析にはClinkerを用いた。

主な結果

• 分類学的同定： 菌株S14は、参照ゲノムGCA_000720535.1に対して96.5%のANIを示すStreptomyces halstediiと同定された。
• 抗真菌活性： S14は一貫してT. atrovirideを抑制し、平均幅1.28 ± 0.38 cmの阻止円を形成した。
• 空間的メタボロミクスの不均一性： 初期の特徴量1,690個のうち、771個を解析用に保持した。プレートブロック統計モデルにより、515個の特徴量が共培養の勾配に沿って有意な空間的変動を示すことが明らかになった。
  • 阻止円（IZ）候補： 7つの堅牢な特徴量がIZに関連していた。特に、m/z 295.1294（グルタミン酸フェニルアラニンと推定）は、他の領域と比較してIZで3.04倍増加した。
  • 真菌界面（TP）候補： 11個の堅牢な特徴量がTP領域に関連していた。m/z 522.3565（1-オレオイルリゾホスファチジルコリンと推定）は、TPで29.20倍増加した。
• ゲノムとメタボロームの収束：
  • カボキサマイシン（Caboxamycin）： 最も強い収束は、m/z 256.0608（カボキサマイシンと推定）とNODE_598の間で観察された。NODE_598は、参照カボキサマイシンBGC（BGC0001444）に対して0.03のBiG-SCAPE距離を持つPKS様BGCである。遺伝子レベルの解析では、参照配列と96.6–97.6%の高い配列同一性と保存されたシンテニーが示された。しかし、この代謝物は主に放線菌側（AD/AC）に関連しており、IZには関連していなかった。
  • バフィロマイシン（Bafilomycin）： ファミリーレベルの関連性が、m/z 743.399（バフィロマイシンC1と推定）とNODE_118（バフィロマイシンB1参照に対しBiG-SCAPE距離0.34）の間で見出された。ゲノムおよびメタボロミクスによる証拠は異なる類縁体を示唆しており、結論はファミリーレベルの所属に限定された。
  • その他のBGC： NRPS、PKS、テルペン、RiPPsを含む59個のBGCが同定された。NODE_8はポリ環式テトラメートマクロラクタム（PTM）BGCと類似性を示したが、対応するメタボロミクス上のマッチは見出されなかった。

意義および主張
本論文は、ハリナシバチ関連の放線菌の文脈において、領域的なメタボロミクスとゲノムマイニングを統合することで、生合成の仮説を優先順位付けする戦略を実証していることを主張している。

• 空間解像度： 本研究は、S. halstediiとT. atrovireの化学的相互作用が空間的に構造化されていることを強調している。プレート全体を均質化するのではなく、特定の領域（IZおよびTP）をサンプリングすることにより、著者らは拮抗表現型に関連する明確な化学的シグネチャーの特定に成功した。
• 因果関係を伴わない優先順位付け： 著者らは、自らの知見が「空間的に構造化された化学的相互作用」を明らかにしていると明記しており、単一の化合物に抗真菌表現型を帰属させることはしていない。本研究は、生物学的に関連する領域（IZおよびTP）における濃縮度に基づき候補（例：グルタミン酸フェニルアラニンおよび1-オレオイルリゾホスファチジルコリン）を優先順位付けしているが、因果関係を確定するには機能検証（単離および生物活性試験）が必要であることを認めている。
• ゲノム潜在能： カボキサマインに非常に類似したクラスターを含む59個のBGCの同定は、S. halstedii S14の生合成能を裏付けている。ゲノムデータ（NODE_598）とメタボロミクスデータ（m/z 256.0608）の収束は、カボキサマイシン生成の収束的な仮説を提供するが、この代謝物が観察された抗真菌活性において果たす役割は未証明である。
• 生態学的文脈： 本研究は、M. monduryの幼虫食が特殊な二次代謝産物のバイオプロスペクティングのための実行可能なソースであることを確立し、腸内細菌叢や蜂蜜に焦点を当てた従来の研究を超えて、昆虫関連放線菌の既知の化学的多様性を拡大した。

総括すると、本論文は、微生物相互作用における化学的不均一性をマッピングし、それらの空間的パターンをゲノム潜在能に結びつけるための方法論的枠組みを提示しており、現在の抗真菌因果関係の証拠を誇張することなく、構造および機能特性評価のための洗練された候補群を提供している。