加勒比海濒危造礁珊瑚麋角珊瑚（Acropora palmata）的种群结构与基因流：基于SNP微阵列的全分布区基因组评估

《Conservation Genetics》：Population structure and gene flow in the endangered Caribbean reef-building coral, Acropora palmata

　　

摘要专业翻译 加勒比海造礁珊瑚麋角珊瑚（Acropora palmata）的种群自21世纪初其遗传结构首次被表征以来已急剧衰退。以往的分析对整个加勒比海及西北大西洋的珊瑚群体进行了全面采样，但受限于所检测位点数量，基因组分辨率有限。这些分析表明该物种存在通过断裂进行的广泛无性繁殖、基因单元（genet）水平上的高异交率，以及东西向的种群分裂。为推进研究并指导这一濒危物种的管理，研究人员利用SNP微阵列（SNP microarray）对来自12个区域的约3,900份样本（代表约1,400个基因单元）在18,898个单核苷酸多态性（single nucleotide polymorphism, SNP）位点上进行了更新后的种群基因组评估。样本由30余个研究和修复团队提供。分析鉴定出12个区域间存在9个空间结构化的遗传聚类（genetic cluster），平均成对FST值较低，介于0.01至0.125之间。值得注意的是，来自佛罗里达的遗留基因单元在两个聚类之间表现出混合祖先（admixture），其中一个聚类主导于西侧的中美洲珊瑚礁带（Mesoamerican Reef Tract），另一个聚类见于从古巴向南的基因单元中。迁移表面分析（migration surface analysis）突显了洋流对基因流的影响。沿大安的列斯群岛（Greater Antilles）存在明显的距离隔离（isolation by distance），但沿佛罗里达珊瑚礁带（Florida Reef Tract）则较弱。野生基因单元在各样点间的亲缘关系（kinship）较低，表明局部亲缘关系有限；然而，修复工作中的辅助有性繁殖导致了独特的亲缘关系模式。这些发现通过解析额外的弱分化遗传聚类，细化了研究人员对Acropora palmata种群结构的理解，并强调了标准化高分辨率基因分型工作流程在生成保护相关数据方面的优势。

论文解读

一、研究背景与意义

加勒比海和西北大西洋的主要造礁珊瑚物种数量已急剧下降，导致群落组成发生显著变化。曾经占主导地位的快速生长、结构复杂的鹿角珊瑚属（Acropora）物种——包括麋角珊瑚（Acropora palmata）和鹿角珊瑚（A. cervicornis）——如今已远不如从前丰富，取而代之的是自然补充仍然成功的物种如Siderastrea siderea和Porites astreoides。这些物种无法维持珊瑚礁框架构建所需的高净钙化速率，导致生态系统从快速生长、结构复杂的Acropora主导系统向更简单、恢复力更弱的珊瑚礁组合转变。最初的衰退由富营养化、沉积、病害暴发、过度捕捞及其他直接人为影响驱动；近年来，不断升高的海表温度加剧了珊瑚白化和死亡，2023年和2024年全球性事件达到灾难性规模。

佛罗里达珊瑚礁带遭受了尤为严重的损失，特别是Acropora palmata。野生种群的基因型多样性自21世纪初以来持续下降。历史估计表明曾有数万至数十万个创始基因单元，但到2020年，野外仅剩73个，另有约120个保存在人工养护环境中。截至2024年，A. palmata和A. cervicornis在佛罗里达珊瑚礁带的野外已功能性灭绝。

两个鹿角珊瑚物种的重大损失伴随着修复活动的增加。常见策略包括生物库保存、将现有群体碎片化用于原位苗圃培育后移植，以及通过配子采集和异交受精进行辅助有性繁殖。配子还被用于在区域间转移等位基因以增强遗传和表型多样性。这些方法需要详细的种群遗传结构知识——包括基因型多样性、亲缘关系和分化程度——以最大化长期恢复力。

此前基于五个微卫星标记的全分布区分析揭示了高遗传多样性、异交以及以波多黎各附近生物物理屏障为界的东西分裂。然而，南加勒比海（如哥伦比亚、委内瑞拉）的有限采样阻碍了对东西区域南部边界的精确划定。基于SNP数据的研究虽然支持区域内位点间低遗传分化的总体模式，但受限于采样位点数量。此外，在种群和珊瑚礁尺度上，任意两个随机采样的野生A. palmata群体之间检测不到亲缘关系，表明近交风险低。然而，该物种的自然有性繁殖现在已失败，需要人工辅助繁殖。加勒比鹿角珊瑚还通过无性繁殖广泛复制：枝条从成体群体上断裂后可在特定条件下重新附着，形成几乎基因相同的大片群体（分株，ramet），这些分株共享同一基因组，可追溯至一次有性繁殖事件，构成一个基因单元。

为支持保护和修复，研究人员开展了这项全分布区种群基因组分析，整合了公开的STAGdb数据库资源，利用标准化高分辨率基因分型工作流程，提供了迄今任何大西洋珊瑚物种中最详细的Acropora palmata种群结构评估。

二、主要技术方法

研究人员利用公开可用的STAGdb数据库中3,925份A. palmata样本的微阵列数据，样本由30余个研究和修复团队贡献，覆盖12个区域约223个用户提交的采样点/珊瑚礁。研究采用基于杂交的珊瑚-藻类基因型微阵列（coral-algal genotype array）和标准化鹿角珊瑚基因分型工具（Standard Tools for Acroporid Genotyping, STAG），在18,898个SNP位点上进行分析。关键方法包括：将样本按基因单元分组（3,925份样本归并为1,432个基因单元）；进行严格的数据过滤（包括基因分型探针质量、缺失率、连锁不平衡剪枝、亲缘关系去除）；采用STRUCTURE和ADMIXTURE进行无监督聚类分析；使用FEEMS进行有效迁移表面分析；利用StAMPP估算成对F_ST；通过Mantel检验评估距离隔离；使用KING方法估算亲缘关系；并整合HYCOM海洋模型的海表流速数据进行生物物理关联分析。

三、研究结果

种群结构分析： 基于地理平衡的子样本数据集（554个基因单元，3,215个SNP位点），主成分分析（PCA）的前两个轴分别解释了11%和4%的变异，第一轴沿东西方向分离采样地点，中美洲珊瑚礁地点（墨西哥、伯利兹、洪都拉斯）大部分重叠。STRUCTURE和ADMIXTURE分析均确定K=9为最优聚类数，鉴定出9个空间结构化的遗传聚类。聚类大致对应地理区域，但存在显著例外：佛罗里达的基因单元表现出混合祖先，反映了一个在佛罗里达和古巴发现的聚类与一个沿中美洲珊瑚礁带发现的聚类之间的贡献。来自古巴的有限基因单元自身也表现出混合，其祖先来源包括巴哈马主导的聚类和与佛罗里达共享的聚类。多米尼加共和国的基因单元在巴哈马/牙买加聚类中具有不同程度的祖先，但分配概率远低于50%。波多黎各和美属维尔京群岛各自形成聚类，库拉索、哥伦比亚大陆和圣安德烈斯群岛的基因单元对其各自聚类具有高分配概率，表明来自其他聚类的基因流较低。

遗传分化与迁移表面： 尽管聚类分析支持全分布区存在基因流屏障，但样点间遗传分化（F_ST）总体较低，范围从伯利兹与洪都拉斯之间的0.001到伯利兹与美属维尔京群岛之间的0.125。有效迁移表面分析显示，从中美洲珊瑚礁带到佛罗里达的基因流几乎没有屏障，而佛罗里达与巴哈马之间的直接扩散受阻，进入佛罗里达的扩散可能经由古巴来自伊斯帕尼奥拉和牙买加。多米尼加共和国与波多黎各之间（莫纳通道）存在基因流屏障，库拉索和哥伦比亚大陆位点各自孤立。这些遗传数据模式与8月（A. palmata产卵期）的主要表层流速和方向高度吻合。

杂合度与近交系数： 所有区域的中位全基因组杂合度估计相似（39.3%–41.6%），哥伦比亚最低（39.3%），牙买加最高（41.6%）。各地点中位近交系数（F_IS）较低，范围从牙买加的?0.008到佛罗里达的0.0446。佛罗里达较高的F_IS值可由两个聚类在佛罗里达珊瑚礁带上的混合（Wahlund效应）解释。

亲缘关系分析： 基于KING亲缘系数分布，所有样点的分布均以0（无关）为中心。除多米尼加共和国两个位点外，大多数推定野生且未受干扰的A. palmata位点不含彼此密切相关的基因单元。然而，已知有修复活动的位点（如库拉索海洋水族馆）表现出更长的尾部和多峰分布，表明人类干预导致高于预期的亲缘关系。

距离隔离分析： 沿大安的列斯群岛（多米尼加共和国、波多黎各、美属维尔京群岛），地理距离与遗传距离（线性化F_ST）之间存在显著关系（Mantel r=0.91，p<0.001，n=91个位点对），但沿佛罗里达珊瑚礁带则不然（Mantel r=0.22，p=0.31，n=45个位点对）。平衡重采样敏感性分析确认了这一差异并非由采样努力不均驱动，而是代表生物学信号。

四、讨论与结论

研究讨论指出，Acropora palmata是加勒比海和西北大西洋仅有的五种具有高净钙化速率的造礁物种之一。基于3,215个全基因组SNP位点，研究将554个基因单元解析为9个遗传聚类，数量多于此前报道，这可归因于本研究使用的位点、样本和采样点数量。聚类大致对应地理区域，但佛罗里达基因单元的混合祖先反映了来自中美洲珊瑚礁带和古巴/佛罗里达聚类的贡献。生物物理建模支持中美洲珊瑚礁带是进入佛罗里达的中美洲/佛罗里达聚类的起源，因为表层流驱动的扩散从中美洲珊瑚礁带到佛罗里达是强单向性的。哥伦比亚-巴拿马涡旋在局部滞留幼虫方面发挥作用，莫纳通道的基因流屏障再次得到验证。沿大安的列斯群岛存在强距离隔离模式，表明邻近位点间的成功繁殖和补充更可能发生，且该区域种群可能处于遗传漂变与基因流之间的平衡状态。佛罗里达珊瑚礁带缺乏距离隔离可能部分由该区域的混合样本和复杂的流型解释。尽管存在生物物理屏障和距离隔离模式，成对F_ST值表明整个加勒比海和西北大西洋的总体分化程度较低，这与该物种幼虫具有较长预感受态期并在缺乏定居线索时持续向上游动的自然史相符。

研究结论指出，A. palmata种群在加勒比海范围内存在9个弱分化的遗传聚类，总体遗传分化低但存在由洋流驱动的空间结构。佛罗里达基因单元表现出中美洲和古巴/佛罗里达聚类的混合祖先，表明历史上存在基因流。野生群体中基因单元间亲缘关系低，而修复活动中的辅助有性繁殖产生了独特的亲缘关系模式。这些发现细化了对该物种种群结构的理解，并强调了标准化高分辨率基因分型工作流程在生成保护相关数据方面的价值，为这一濒危物种的保护和修复策略提供了重要的遗传学依据。