造礁蜂窝蠕虫在广泛纬度梯度下的高基因流中的适应性变异

《Molecular Ecology》：Adaptive Variation Amid High Gene Flow Across a Broad Latitudinal Gradient in the Reef-Building Honeycomb Worm, Sabellaria alveolata

　　

气候变化已经在改变，并将继续重塑海洋种群的持久性和空间结构。因此，理解集合种群连通性和局部适应对于识别脆弱种群和具有适应韧性的种群至关重要。研究人员应用海景基因组学方法研究了蜂窝蠕虫的连通性模式和适应反应——这是一种支持高生物多样性并提供基本生态功能的潮间带生态系统工程师。对来自跨越整个纬度范围的17个地点的286个个体进行了基因分型，以评估种群结构和局部适应。异常分析（BayeScan和OutFlank）识别出6个适应性单核苷酸多态性（SNP）和2117个中性SNP。中性SNP分析揭示了从低到中等的遗传分化，表明种群连通性良好。在南部边缘和比斯开湾种群中与可能的祖先基因库的混合表明独特遗传多样性的持续存在。虽然连通性可能是由自然幼虫扩散驱动的，但由支援船只、水产养殖、海洋碎片以及不断扩大的海上和沿海基础设施导致的人类介导的运输也可能促进基因流，从而塑造了观察到的遗传模式。尽管存在这种高连通性，研究人员检测到了对温度、盐度和净初级生产力的适应遗传特征，适应性SNP与参与这些环境因素的分子和生理反应的基因相关联。研究结果表明，Sabellaria alveolata可能对气候变化和其他人为压力具有韧性，使其成为研究海洋种群韧性的遗传和生理基础的宝贵模型。

研究背景与引言
群体连通性和基因流对物种的适合度与长期存活至关重要，它们有助于维持遗传多样性，并促进增强抵御环境变化、过度开发和疾病等压力的韧性。传统观点认为，具有较长浮游幼虫期（PLD）的海洋广播式产卵者通常表现出高度的遗传混合。然而，近期研究对这一假设提出了挑战，发现即使在高扩散潜力的物种中，细尺度遗传结构和适应性分化也比此前认为的更为普遍。因此，研究重心正从单纯评估中性遗传多样性，转向同时识别受选择的适应性位点。在更新世冰川作用的影响下，许多陆地和海洋物种在更新世后发生了种群瓶颈和重新定殖，导致北部边缘遗传多样性降低，南部边缘遗传多样性较高。以往的研究多集中于商业或保护性关注物种，且空间尺度有限，忽视了对物种整个地理分布范围的边缘种群的探索，特别是较难触及的南部边缘种群。这种信息的缺失阻碍了研究人员全面评估自然种群的适应潜力。边缘种群受环境不稳定性塑造，作为增强极端条件韧性的适应性特征进化的前哨，具有更高的遗传变异性，是生物遗传资源的宝贵储库。这对于潮间带生态系统工程师尤为重要，它们缓冲温度、干燥和缺氧等极端环境条件，支撑高生物多样性。如果它们对快速环境变化适应不良，可能导致局部灭绝，破坏生态系统结构和功能，进而对生物多样性和整体生态稳定性产生级联影响。

蜂窝蠕虫是一种从摩洛哥北部到苏格兰西南部广泛分布的造礁潮间带多毛类。作为生态系统工程师，它将粗砂和生物碎屑巩固成硬化管状结构，为无数分类群提供复杂的生物栖息地。此前对该物种的研究多集中在法国比斯开湾的北部和中部种群，揭示了显著的遗传结构，但南部边缘种群仍属研究不足。此外，尽管生理学变异（如膜脂组成和繁殖生理性状的差异）暗示了局部适应的存在，但目前的遗传学研究尚未提供确凿的局部适应证据，且多关注升温，对低温、食物和悬浮颗粒物可用性等选择压力关注较少。针对这些问题，研究人员开展了研究，旨在揭示群体结构、局部适应以及通过基因型-环境关联驱动选择的非生物和生物因素。研究得出结论：尽管存在高连通性，蜂窝蠕虫仍表现出对温度、盐度和净初级生产力的局部适应。该研究发表在《Molecular Ecology》上，为理解该关键生态系统工程师的分布范围内遗传结构、局部适应和生态韧性提供了宝贵见解，对理解S. alveolata如何响应环境变化具有重要意义。

关键技术方法
样本队列来源于2018年3月至2019年9月期间，从摩洛哥至英国北部跨越其已知纬度分布的17个地点采集的286个成年Sabellaria alveolata个体。研究人员主要应用了2b-RAD基因组扫描技术对个体进行基因分型，并利用Illumina NextSeq550进行测序。通过比对参考基因组，使用Stacks等软件识别单核苷酸多态性（SNP）。采用OutFlank和BayeScan进行异常分析以区分中性位点和适应性位点。利用主成分分析（PCA）、判别分析主成分分析（DAPC）和STRUCTURE软件分析群体结构与遗传聚类。通过偏冗余分析（pRDA）进行基因型-环境关联分析，并结合BLASTp和SWISS-MODEL进行基因功能注释和蛋白质结构同源性建模。

研究结果
3.1 遗传多样性与群体结构
研究人员通过计算近交系数（F_IS）和全局固定指数（F_ST）并利用PCA、DAPC和STRUCTURE分析，发现所有种群基本为远交，整体表现出低至中等的遗传分化。大部分种群基因型高度相似，但在南部边缘种群（摩洛哥的Essaouira和Taghazout）和比斯开湾部分种群（La Fontaine aux Bretons和Oléron）中观察到显著的遗传分化和独特的遗传多样性。这表明这些区域受局部海洋过程（如季节性上升流）和种群历史因素影响，基因流较弱。

3.2 与环境参数的关联
研究人员通过偏冗余分析（pRDA）检验适应性等位基因频率与环境参数的关系，发现在最终模型中，平均海水温度是最具影响力的变量，解释了21.4%的总方差。平均潮汐振幅、平均净初级生产力和海水盐度也分别解释了13.7%、13.73%和5.91%的方差。流速和悬浮无机颗粒物浓度无统计学显著性。这表明温度是主要的选择压力，但水动力学和资源可用性等环境因素也塑造了该物种的适应性变异。

3.3 受选择位点的假定基因功能
研究人员通过序列和结构同源性分析，对六个受定向选择的SNPs进行了功能注释。SNP1与转录起始因子TFIID同源，与温度强烈相关，可能参与热胁迫下的转录调控；SNP2与信号黏蛋白HKR1同源，与温度和盐度呈负相关，可能参与渗透压胁迫响应；SNP4和SNP6与高净初级生产力相关，分别编码尿苷5'-一磷酸合酶和金属蛋白酶，可能涉及有机物降解和细胞代谢生物合成；SNP5编码tRNA修饰酶，可能在多种环境胁迫下维持蛋白质合成。此外，约60%的SNPs受到平衡选择，为物种在高度动态的潮间带环境中提供了“基因组缓冲”，有助于应对环境波动。

讨论与结论
在讨论部分，研究表明蜂窝蠕虫在其分布范围内保持高水平连通性，这可能是通过踏脚石扩散、冰期后混合以及人类介导的运输（如船舶压载水、水产养殖设施和海洋废弃物）共同维持的。南部边缘和比斯开湾种群保留的独特遗传多样性反映了局部海洋学过程的影响。研究首次提供了该物种适应性的遗传证据，海水温度是主要的适应性变异驱动因素。受选择的基因功能与其在温度、盐度、净初级生产力等环境胁迫下的生理和分子响应相符。高比例的平衡选择位点作为该物种在潮间带环境中的“基因组缓冲”，为其提供了应对快速环境转变所需的进化基础。

在结论部分，全基因组范围内的SNP分析表明，蜂窝蠕虫的现有种群在欧洲大西洋沿岸保持紧密连通。基因流跨越长距离发生，可能由人类介导的扩散促进。没有种群表现出遗传多样性降低或隔离迹象。尽管存在潜在的破坏性基因流，仍观察到了局部适应。六个SNPs显示出选择特征，主要与温度、盐度和净初级生产力相关。与这些SNPs相关的基因功能证实了对这些环境选择驱动因素的适应。研究结果表明，蜂窝蠕虫能够对环境因素产生局部适应，且可能对环境变化具有韧性，突显了其作为研究潜在气候变化“赢家”模式物种的价值。