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）之间的种群遗传结构联系，以及在此背景下宿主-藻类共生的特异性，目前仍知之甚少。理解这些机制对于评估主动干预策略（如辅助基因流）以增强其他地区珊瑚韧性的潜力至关重要。

研究方法
本研究采用群体基因组学方法，对三个位点的 P. harrisoni 进行了调查：两个位于热极端环境的南侧 PAG 位点（Al Saada Reef 和 Saadiwad Reef），以及一个位于环境较温和的阿曼湾（Al Aqah）。

• 采样与测序： 2022 年共采集了 119 株珊瑚群体。使用 Illumina NovaSeq6000 进行全基因组鸟枪法（WGS）测序，平均覆盖度为 45X。
• 变异检测与过滤： 按照 GATK 最佳实践处理读取数据。应用了两种不同的 SNP 过滤流程：一种用于群体结构分析（去除低频变异和克隆），另一种是更宽松的流程用于检测选择性扫荡（保留稀有变异以捕捉固定信号）。
• 群体结构： 分析包括主成分分析（PCA）、ADMIXTURE 祖先成分估计、成对 $F_{ST}$ 计算以及用于重建人口统计历史的 Stairway Plot v2。
• 宿主-共生体关联： 使用 Mantel 检验评估基因型-基因型关联，将宿主基因型与藻类共生体组合（ITS2 类型谱）进行相关性分析。
• 选择信号： 通过整合 Weir & Cockerham 的 $F_{ST}$（遗传分化）和跨群体扩展单倍型纯合性（XP-EHH）统计量，识别候选正向选择区域。使用 SnpEff 对离群 SNP 进行注释，并评估核苷酸多样性（$\pi$）和连锁不平衡（LD）衰减，以表征选择性扫荡的性质。

主要结果

• 群体结构： 三个位点之间的整体遗传分化较低（$F_{ST}$ 范围为 0.021 至 0.041），表明通过霍尔木兹海峡存在基因流。然而，ADMIXTURE 分析（$K=2$）显示，PAG 种群与 GO 相比具有独特的祖先谱系，其中最南端的 PAG 位点（SA）表现出极少的混合。人口统计重建表明，有效种群大小（$N_e$）在全新世期间 PAG 形成（约 6,000 年前）时出现了显著下降，这表明存在奠基者效应或瓶颈效应。
• 宿主-共生体特异性： 虽然宿主遗传相似性并未严格决定共生体身份，但研究发现宿主与 Cladocopium 遗传距离之间存在中度但显著的关联。PAG 珊瑚主要与耐热的 Cladocopium thermophilum 相关联，而 GO 珊瑚则托管着 Durusdinium trenchii 和 Symbiodinium microadriaticum。
• 选择的基因组特征： 研究确定了 PAG 珊瑚中存在两个在 GO 珊瑚中未观察到的强正向选择特定基因组区域：
  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 加工的基因。
• 选择的性质： 等位基因频率图显示，特定等位基因在 PAG 珊瑚中趋于近乎固定，但在 GO 珊瑚中并非如此。核苷酸多样性（$\pi$）并未出现明显的降低，却伴随着扩展的单倍型，这表明是一种“软”选择性扫荡，即选择作用于现有的遗传变异而非单一的新发突变。

意义与主张
本文认为，P. harrisoni 在极端 PAG 环境中的持续存在是由作用于一小组基因组位点的强定向选择所支撑的。研究确定的“非协商性 SNP”（被驱动至 PAG 种群近乎固定的等位基因）强调了这些特定遗传变异在生存极端热应力和盐度压力中的核心作用。

研究声称，这是识别未来主动干预策略（如通过选择性育种进行的辅助基因流）以增强珊瑚韧性的第一步。然而，作者强调了一个关键限制：PAG 珊瑚的适应性与海湾特定的环境条件（包括高盐度）紧密相连。因此，将这些基因型转移到较冷、低盐度的环境（如 GO）在长期内可能无法实现。研究结果强调，虽然区域遗传多样性是一项宝贵的资源，但极端环境下的狭窄选择机制可能会限制未来的适应潜力，这使得全球气候缓解策略比单纯的局部干预更加紧迫。