大堡礁磁岛（Yunbenun）广布产卵珊瑚Acropora kenti的遗传隔离与多样性丧失

《Conservation Genetics》：Genetic isolation and loss of diversity in the broadcast spawning coral, Acropora kenti at Yunbenun (Magnetic Island) on the Great Barrier Reef

　　

中文摘要翻译： 大堡礁（Great Barrier Reef, GBR）的近岸珊瑚礁种群正面临着来自环境胁迫因子的日益增长的压力。遗传多样性和连通性是决定珊瑚适应能力的关键因素。Acropora kenti在大堡礁沿岸数百公里范围内形成了至少四个遗传上截然不同的集群。其中一个因地理分布受限而可能较为脆弱的集群位于Yunbenun（磁岛，Magnetic Island）。在此，研究人员对来自Yunbenun周边七个珊瑚礁以及中央大堡礁八个邻近近岸和中陆架珊瑚礁的A. kenti种群进行了遗传评估，以确定其种群结构、遗传多样性和连通性。利用单核苷酸多态性（single nucleotide polymorphisms, SNPs）标记，研究人员证明一个遗传上截然不同（Fst >0.33）的谱系在Yunbenun周围的所有海湾占据主导地位，但在邻近珊瑚礁上却很少发现。对Yunbenun谱系内部进行的精细尺度分析揭示了由有限基因流驱动的距离隔离模式。由水动力模型支持的混合分析表明，Yunbenun与周围珊瑚礁之间存在幼虫扩散，但基因流仍然有限。与邻近珊瑚礁（Ne = 9494）相比，Yunbenun种群的等位基因丰富度（allelic richness）（p < 0.001）更低，有效种群大小（effective population size, Ne）也小得多（北部海湾Ne = 1348；南部海湾Ne = 91）。这些发现表明，尽管A. kenti在整个大堡礁的遗传多样性仍然很高，但诸如Yunbenun这类孤立且可能已局部适应的种群可能面临遗传衰退的风险，这将削弱其适应潜力并增加其保护关注度。

论文解读文章：

一、研究背景与问题

大堡礁（Great Barrier Reef, GBR）作为联合国教科文组织世界遗产地，是全球最大且最多样化的珊瑚礁生态系统之一。然而，气候变化引发的更频繁、更强烈的热浪已导致大规模珊瑚白化事件，对大堡礁构成日益严重的威胁。此外，靠近陆地的近岸珊瑚礁还受到水质恶化、沉积物暴露增加、藻类占优势、气旋和漂白事件等局部压力的影响。珊瑚种群的持续存在依赖于其对环境干扰的耐受性和恢复力，而这在一定程度上通过适应过程以及遗传连通种群间的混合来实现。

在珊瑚遗传学研究中，Acropora kenti（曾与A. tenuis合并为同义词）作为大堡礁上丰富且重要的造礁物种，其种群遗传结构受到广泛关注。已有研究表明，A. kenti在大堡礁沿岸数百公里范围内形成了至少四个遗传上截然不同的集群。其中，位于Yunbenun（磁岛，Magnetic Island）的集群因地理分布受限而可能较为脆弱。Yunbenun位于中央大堡礁的近岸区域，其周围珊瑚礁受到水质差、沉积物暴露、藻类占优势、气旋和漂白事件等多重压力的影响。先前基于全基因组序列的群体遗传学模型显示，Yunbenun的A. kenti种群可能在大约0.5至150万年前与周围珊瑚礁的种群发生分化，经历了历史上的遗传隔离时期。

由于YunbenunA. kenti种群的地理分布范围受限且遗传独特性显著，加上近岸珊瑚礁面临的累积急性与慢性压力，亟需深入了解该种群的遗传结构、遗传多样性和连通性，以评估其对当前和未来胁迫因子的脆弱性。小型、地理上孤立的种群特别容易因种群规模小而丧失多样性，这可能导致近交衰退并降低适应能力。基因流是影响YunbenunA. kenti遗传脆弱性的关键因素，包括Yunbenun各海湾之间的局部基因流以及与邻近广泛分布谱系的基因流。因此，该研究旨在探究Yunbenun周围A. kenti珊瑚的遗传结构（包括先前未检测的采样点）、Yunbenun与邻近珊瑚礁之间的基因流和遗传连通性，以及Yunbenun的A. kenti种群是否在遗传上不同于其他邻近珊瑚礁位置。同时，研究人员还利用水动力模型研究当地和区域洋流如何影响幼虫扩散、迁移模式和潜在的连通性障碍，以帮助细化该种群的地理和遗传多样性、连通性及隔离状况。

二、研究内容与主要结论

研究人员在中央大堡礁的15个珊瑚礁种群中采集了总共512个A. kenti样本，采样时间跨越6年（2017年、2020-2022年）。采样地点分为四类：（1）Yunbenun北部，（2）Yunbenun南部，（3）邻近Yunbenun的中陆架珊瑚礁，（4）邻近Yunbenun的近岸珊瑚礁。Yunbenun北部和南部类别包括来自岛屿周围七个地点的171个样本。通过单核苷酸多态性（SNPs）标记分析，研究人员发现一个遗传上截然不同的谱系在Yunbenun周围所有海湾占主导地位，但在邻近珊瑚礁上很少发现。精细尺度分析揭示了Yunbenun谱系内部的距离隔离模式，主要由有限基因流驱动。混合分析结合水动力模型表明，Yunbenun与周围珊瑚礁之间存在幼虫扩散，但基因流仍然有限。Yunbenun种群的等位基因丰富度显著低于邻近珊瑚礁，有效种群大小也远小于邻近珊瑚礁。这些发现表明，尽管A. kenti在整个大堡礁的遗传多样性仍然很高，但Yunbenun这类孤立且可能已局部适应的种群可能面临遗传衰退的风险。

三、主要技术方法

研究人员采用了多种关键技术方法。首先，利用DArTseq技术进行简化基因组测序（reduced representation sequencing），获得全基因组范围的SNP标记。其次，使用多种群体遗传学分析软件（如Structure、PCA、NewHybrids等）进行种群结构、混合分析和杂交个体鉴定。第三，利用NeEstimator软件基于连锁不平衡（linkage disequilibrium, LD）方法估算有效种群大小（N_e）。第四，利用bcftools识别长纯合片段（long runs of homozygosity, LROH）以评估近交水平。第五，构建二维水动力数值模型（Delft3D FM Suite），模拟珊瑚幼虫的扩散轨迹和连通性。样本来源包括Yunbenun周边七个地点（北部海湾：Huntingfield Bay、Wilson Bay、Maud Bay、Horseshoe Bay；南部海湾：Geoffrey Bay、Picnic Bay、Middle Reef）以及邻近的中陆架珊瑚礁（Bramble、John Brewer、Keeper）和近岸珊瑚礁（Palm Island群岛的Orpheus Island、Pelorus Island、Havannah Reef）。

四、研究结果

（一）种群结构

主成分分析（PCA）和基于Structure的混合分析（K=2）清晰显示了YunbenunA. kenti种群与邻近近岸和中陆架珊瑚礁种群之间的遗传分化。Yunbenun集群包括岛屿周围采样的所有六个海湾以及Middle Reef。邻近近岸和中陆架珊瑚礁的珊瑚形成了一个与Yunbenun明显不同的单一遗传集群。这种强烈的遗传分化在成对F_st值中也得到体现，Yunbenun与邻近珊瑚礁之间的比较值（0.334至0.348）远大于Yunbenun海湾之间（0-0.01）或邻近珊瑚礁群内部（0-0.0028）的值。

混合分析显示，大多数个体几乎完全归属于Yunbenun或邻近珊瑚礁集群。在邻近珊瑚礁（Havannah Reef）中发现了三个具有Yunbenun血统的推定迁移个体，在Yunbenun（Horseshoe、Huntingfield和Wilson Bay）中发现了13个具有邻近珊瑚礁血统的个体。这些迁移个体集中的地点代表Yunbenun北部和地理上最近的近岸邻近珊瑚礁（Havannah Reef）。此外，还鉴定出七个高度混合个体（次要混合比例>0.05），NewHybrids分析显示其中一个来自Horseshoe Bay的个体可能是早期世代杂交个体。

精细尺度分析显示，Yunbenun各海湾之间的遗传结构比邻近珊瑚礁群内部位置之间的结构更强。尽管邻近珊瑚礁群中珊瑚礁之间距离长达75公里，但未发现距离隔离模式，PCA中也没有明显的结构。相比之下，Yunbenun各海湾之间（相距5至15公里）显示出显著但嘈杂的距离隔离模式，主要由北部和南部海湾之间的高F_st值驱动。PCA揭示了三个部分重叠的集群，北部海湾（Huntingfield、Wilson、Maud和Horseshoe Bay）局限于一个中央集群，而南部海湾（Middle Reef、Geoffrey和Picnic Bay）的样本分布在所有三个集群中。

（二）遗传多样性与有效种群大小

稀薄化等位基因丰富度（A_r）在Yunbenun（A_r 1.54）显著低于邻近珊瑚礁（A_r 1.65）（p < 0.001），但Yunbenun北部和南部海湾之间没有差异。期望杂合度（H_e）在Yunbenun和邻近珊瑚礁之间没有显著差异，但在与杂合子过量相关的两个指标（H_o和F_IS）上观察到显著差异。南部Yunbenun的H_o值最高、F_IS值最低，北部Yunbenun居中，邻近珊瑚礁的H_o最低（F_IS最高）。这些模式具有统计学显著性，但数值较小。

基于连锁不平衡的当代有效种群大小（N_e）估计显示，Yunbenun的N_e（北部1352；南部91）比邻近珊瑚礁（9504）低数个数量级。Jackknife 95%置信区间反映了邻近珊瑚礁估计值的置信度较低（最小2217至最大40656），而Yunbenun种群的置信区间较窄（北部445-4078；南部44-190）。这些N_e模式对样本量、样本合并方案、稀有等位基因阈值和排除推定亲属均具有稳健性。

在检测到的21个高Phred评分的长纯合片段（LROH）中，三个最高评分的LROH（反映更长片段和更高纯合性）在Yunbenun发现，但这一模式在统计上不显著，且各群体间总LROH长度没有显著差异。

（三）生物物理连通性

水动力模型显示了在Yunbenun周围十一个海湾和六个邻近珊瑚礁释放的珊瑚产卵在60天内的扩散情况。模型表明，Yunbenun的幼虫扩散通常保持在珊瑚礁地点附近，不与邻近珊瑚礁重叠。大多数北部亚种群的扩散范围不与南部亚种群重叠，反之亦然。少数例外包括Nelly和Geoffrey Bay，它们都向北扩散到Huntingfield Bay。中陆架珊瑚礁的珊瑚产卵向东南方向扩散至Yunbenun，但在60天模拟期内未完全到达该地点。Palm Island近岸珊瑚礁（Havannah Reef、Orpheus和Pelorus Island）释放的幼虫主要向西北方向移动，并与中陆架珊瑚礁混合。

五、讨论与结论

研究确认了Yunbenun的A. kenti种群与邻近珊瑚礁种群在遗传上截然不同，其地理分布范围很可能高度受限，仅限于岛屿周围的裙礁和附近的Middle Reef。尽管存在一些基因流的迹象，但推定迁移个体相对稀少（427个个体中的16个），且集中在Yunbenun北部海湾和Havannah Reef。高度混合个体的数量相对于推定迁移个体数量较少，这可能反映了迁移率的较大波动，或迁移个体因远交衰退或缺乏局部适应性性状而导致繁殖成功率降低。

Yunbenun北部和南部亚种群之间的遗传结构差异可能由多种因素导致，包括水动力模型显示的北部和南部亚种群幼虫扩散范围仅在海岛东北和西北点重叠，以及南部海湾极低的N_e加速了遗传漂变，可能导致随机补充效应驱动的虚假遗传聚类模式。

有效种群大小的分析确认了邻近珊瑚礁种群规模庞大，而Yunbenun的N_e要小得多，特别是南部海湾。敏感性分析支持了南部Yunbenun低N_e的发现。等位基因丰富度显著低于邻近珊瑚礁，而期望杂合度没有显著变化，这表明遗传漂变已导致等位基因丢失，但尚未显著影响杂合度。观察到的杂合度小幅增加可能反映了极小N_e下二项式抽样效应导致的预期增加。

该研究的发现揭示了Yunbenun的A. kenti种群与邻近珊瑚礁在遗传上截然不同，基因流有限，遗传多样性中等。Yunbenun北部和南部亚种群之间的遗传差异强调了在珊瑚礁管理中考虑精细尺度遗传结构的重要性。受限的连通性、低有效种群大小和中等近交水平突出了该种群对环境胁迫因子的脆弱性。Havannah Reef成为北部Yunbenun与中央大堡礁其他近岸和中陆架珊瑚礁之间潜在的遗传踏脚石。这些结果强调了维持高有效种群大小和遗传多样性以确保YunbenunA. kenti种群及其周围种群长期生存力和对环境变化恢复力的重要性。该研究为制定考虑当地珊瑚种群独特遗传特征的有针对性的保护策略奠定了基础，最终有助于这些重要生态系统的恢复力和保护。