Persistence of the coral Porites harrisoni in the Persian/Arabian Gulf is underpinned by strong selection on few genomic loci

技術要約：ペルシャ湾／アラビア湾における Porites harrisoni の生存

問題提起
造礁サンゴは熱的限界に近づいており、加速する海洋温暖化の下で壊滅的な損失が生じることが予測されている。ペルシャ湾／アラビア湾（PAG）は、夏季の海水温が34°Cを超え、他の熱帯のリーフの熱限界を2°C以上上回る極限環境である。このような過酷な条件下にもかかわらず、造礁サンゴの一種である Porites harrisoni はPAGにおいて生存し続けている。しかし、この熱耐性のゲノム基盤、およびPAGの個体群と近隣のより穏やかな環境（オマーン湾［GO］など）を結ぶ集団遺伝学的構造、そしてこの文脈における宿主と藻類の共生の特異性は、依然として十分に理解されていない。これらのメカニズムを理解することは、他地域での回復力を高めるための「支援的遺伝子流動（assisted gene flow）」のような能動的な介入戦略を評価する上で極めて重要である。

方法論
本研究では、熱的に極端な南PAG（Al Saada ReefおよびSaadiyat Reef）の2地点と、環境的に穏やかなオマーン湾（Al Aqah）の1地点における P. harrisoni を調査するために、集団ゲノミクス・アプローチを用いた。

• サンプリングおよびシーケンシング： 2022年に計119個体のサンゴコロをサンプリングした。Illumina NovaSeq6000を用いた全ゲノムショットガン（WGS）シーケンシングを実施し、平均カバレッジは45Xであった。
• バリアント・コーリングおよびフィルタリング： リードの処理はGATK Best Practicesに従った。2つの異なるSNPフィルタリング・パイプラインを適用した。一つは集団構造解析用（低頻度のバリアントおよびクローンを除去）、もう一つは選択的スイープ検出用の、より寛容なパイプライン（固定シグナルを捉えるために稀なバリアントを保持）である。
• 集団構造： 解析には、主成分分析（PCA）、祖先推定のためのADMIXTURE、ペアワイズ $F_{ST}$ 計算、および人口統計学的履歴再構築のためのStairway Plot v2を用いた。
• 宿主－共生体関連性： 宿主の遺伝子型と藻類共生体の集合（ITS2型プロファイル）との相関を、ジェノタイプ－ジェノタイプ間の関連性を評価するためのマンテル検定を用いて分析した。
• 選択のシグネチャー： 正の選択を受けている候補領域は、Weir & Cockerhamの $F_{ST}$（遺伝的分化）と、集団間拡張ハプロタイプホモ接合性（XP-EHH）統計量を統合することで特定された。アウトライヤーSNPはSnpEffを用いてアノテーションを行い、塩基多様性（$\pi$）および連鎖不平衡（LD）の減衰を評価して、選択的スイープの性質を特徴づけた。

主な結果

• 集団構造： 全体的な遺伝的分化は低く（$F_{ST}$ は0.021から0.041の範囲）、ホルムズ海峡を越えた遺伝子流動を示唆している。しかし、ADMIXTURE解析（$K=2$）により、PAGの個体群はGOと比較して明確に異なる祖先系統を共有していることが明らかになり、最も南のPAG地点（SA）では混合が最小限であった。人口統計学的再構築は、PAGの形成（約6,000年前の完新世）と一致する有効個体数（$N_e$）の著しい減少を示しており、これは創始者効果またはボトルネック現象を示唆している。
• 宿主－共生体特異性： 宿主の遺伝的類似性が共生体の同一性を厳密に決定するわけではなかったが、宿主と Cladocopium の遺伝的距離との間に、中程度ではあるが有意な関連が見られた。PAGのコロニーは主に熱耐性のある Cladocopium thermophilum と関連していたのに対し、GOのコロニーは Durusdinium trenchii および Symbiodinium microadriaticum を宿主としていた。
• ゲノムにおける選択のシグネチャー： 本研究では、GOのサンゴには見られない、PAGのサンゴにおいて強い正の選択を受けている2つの特定のゲノム領域を特定した。
  1. Contig 29： 最も高い $F_{ST}$ 値（最大0.93）と拡張されたハプロタイプを示した。この領域には、細胞外マトリックス（ECM）のリモデリング（例：COL4A1, COL2A1）、TGF-$\beta$ シグナル伝達（例：ACVR2B, BMPR2）、および細胞周期調節（例：CDK2, ORC4）に関連する遺伝子が含まれている。
  2. Contig 62： 非常に長いハプロタイプ（正規化されたXP-EHHは最大10.47）と高い分化（$F_{ST}$ は最大0.69）を示した。この領域は、免疫応答、ストレス応答（例：KDELC1）、アポトーシス（例：PDCD6）、およびRNAプロセシングに関与する遺伝子が豊富である。
• 選択の性質： アレル頻度プロットは、GOのサンゴには見られない、PAGのサンゴにおける特定のアレルのほぼ完全な固定を示した。拡張されたハプロタイプを伴いながらも顕著な塩基多様性（$\pi$）の減少が見られないことは、「ソフト」な選択的スイープ、すなわち単一の新規変異ではなく、既存の遺伝的変異に対して選択が働いたことを示唆している。

意義および主張
本論文は、PAGにおける P. harrisoni の生存は、少数のゲノム遺伝子座に作用する強い方向性選択によって支えられていると主張している。「非交渉的SNP（non-negotiable SNPs）」――すなわちPAG個体群においてほぼ固定へと駆り立てられたアレル――の特定は、極端な熱および塩分ストレスを生き抜くための、これら特定の遺伝的変異の不可欠な役割を浮き彫りにしている。

本研究は、選択的育種を通じた支援的遺伝子流動のような、将来の能動的な介入戦略のためのゲノム標的を特定するための第一歩を提供すると主張している。しかし、著者らは重要な制約についても強調している。すなわち、PAGのサンゴの適応は、湾内の特定の環境条件（高塩分を含む）と密接に結びついている。したがって、これらの遺伝子型を、より低温で低塩分の環境（GOのような）に移転することは、長期的には実現可能ではない可能性がある。これらの知見は、地域の遺伝的多様性が貴重な資源である一方で、極限環境における狭い選択圧が将来の適応能力を制限する可能性があることを強調しており、局所的な介入だけでなく、地球規模の気候緩和策がより緊急であることを示している。