北太平洋九种岩鱼（Sebastes spp. 和 Sebastolobus spp.）的隐存多样性及种群遗传结构模式——基于低覆盖度全基因组重测序

《Molecular Ecology》：Cryptic Diversity and Patterns of Structure in Nine Species of North Pacific Rockfishes (Sebastes spp. and Sebastolobus spp.)

了解栖息地和生活史特征如何影响种群遗传结构，对于评估保护风险和确定野生种群的最优管理至关重要。来源于适应性辐射（adaptive radiation）的物种可作为研究此类模式的有用模型，因为存在一些一致的特征，使得能够更细致地研究差异属性。在此，研究人员研究了岩鱼（Sebastes 和 Sebastolobus 属）的种群结构；这是一个物种丰富、近期分化、并表现出从隐存物种到假定遗传泛交（genetic panmixia）连续统的海洋鱼类类群。本研究分析了来自9种岩鱼的低覆盖度全基因组重测序数据，识别出四种主要类型的遗传分化：隐存多样性（cryptic diversity）、距离隔离（isolation-by-distance, IBD）、无可检测结构（假定泛交），以及无可检测结构但单个基因组区域推测为基础局域适应（local adaptation）。两种岩鱼——暗色岩鱼（Dark Rockfish, Sebastes ciliatus）和北方岩鱼（Northern Rockfish, S. polyspinis）——各自表现出由两个同域分布且在整个基因组所有染色体上分化的遗传类群组成的隐存多样性。阿拉斯加东南部的黑岩鱼（Black Rockfish）也存在两个遗传类群，尽管存在从北美西海岸到阿拉斯加湾（Gulf of Alaska, GOA）的整体IBD信号。五种岩鱼在整个采样范围内显示出遗传泛交的证据，而黄眼岩鱼（Yelloweye Rockfish）在约4000 km范围内（从南加州到GOA）表现出高基因流，但在萨利什海（Salish Sea）与外海岸之间存在强烈分化。在未检测到结构的物种中，仅黑斑岩鱼（Blackspotted Rockfish）具有一个表现出潜在局域适应特征的离群区域。研究结果突显了多种机制可能促进岩鱼的遗传结构，并提示深度分布和离岸栖息地可能与岩鱼更广泛的分化模式存在联系。此外，本研究为东北太平洋商业渔业中最丰富的部分岩鱼提供了遗传数据，增强了可用于指导渔业管理的科学信息。

研究背景与科学问题

海洋鱼类种群结构常由生活史特征与历史过程共同塑造，并因扩散屏障少、有效种群大或近期分化而常呈现低分化甚至接近泛交。岩鱼（rockfishes，Sebastes spp. 和 Sebastolobus spp.）是东北太平洋常见且经济重要的鱼类，作为近期辐射分化、物种丰富的类群，其遗传结构表现出从隐存物种到假定遗传泛交（genetic panmixia）的连续谱。以往研究多基于微卫星等有限标记，并按采集地点估计等位基因频率，可能掩盖个体水平上的隐存分化。为厘清种群结构类型及其与生活史特征的关系，并为渔业管理提供依据，研究人员对9种岩鱼开展了低覆盖度全基因组重测序（low-coverage whole genome resequencing, lcWGS）的个体水平分析。

研究内容与总体结论

研究人员共分析995份岩鱼样本，涵盖Dark Rockfish（Sebastes ciliatus）、Dusky Rockfish（S. variabilis）、Northern Rockfish（S. polyspinis）、Shortraker Rockfish（S. borealis）、Blackspotted Rockfish（S. melanostictus）、Rougheye Rockfish（S. aleutianus）、Yelloweye Rockfish（S. ruberrimus）、Black Rockfish（S. melanops）和Shortspine Thornyhead（Sebastolobus alascanus）。样本来自阿拉斯加渔业科学中心、阿拉斯加渔猎部、NOAA西南渔业科学中心和加拿大渔业与海洋部等调查/研究收藏。总体识别出四类遗传分化：隐存多样性、距离隔离（isolation-by-distance, IBD）、无可检测结构（假定泛交），以及无可检测结构但存在一个推测局域适应（local adaptation）的基因组区域。研究指出，深度分布和离岸栖息地可能与分化模式存在潜在联系，并为东北太平洋商业渔业中多种高丰度岩鱼提供了新的遗传信息。

主要技术方法

每物种独立处理原始序列，经接头修剪后比对到韩国岩鱼（S. schlegelii）染色体级参考基因组（GCA_014673565.1）；用ANGSD计算基因型似然（genotype likelihoods, GLs），并基于GLs进行PCAngsd全基因组主成分分析（principal component analysis, PCA）和NGSadmix混合分析；按PCA聚类或地理区域分组，用realSFS估算成对F_ST（固定指数），并用核平滑基因组扫描识别分化峰；对IBD使用海洋距离与F_ST/(1?F_ST)相关分析；对显著峰进行连锁不平衡（linkage disequilibrium, LD）、局部PCA和等位基因剂量分析。测序深度平均为2.1×–3.3×。

研究结果

3.1 Rougheye Rockfish

基因组PCA显示从俄勒冈至白令海无地理聚集，admixture分析未支持多个祖先群，基因组扫描无显著分化峰，平均F_ST为0.011–0.018，结论为未检测到种群结构。

3.2 Shortraker Rockfish

PCA沿PC2出现两个非地理聚集的遗传类群，经基因组扫描判别为可能位于chr19的染色体倒位；去除该染色体后PCA无聚集，但仍显示测序批次效应；admixture主要受批次影响；基于地理区域的基因组扫描无显著峰，成对F_ST为0.007–0.010。

3.3 Dusky Rockfish

PCA从阿拉斯加东南部至西部GOA无聚集，admixture多为混合祖先，基因组扫描无显著分化区域，平均F_ST为0.008–0.010。

3.4 Shortspine Thornyhead

PCA从南加州至白令海无聚集，admixture结果不支持稳定分组；基因组扫描无显著峰，平均F_ST为0–0.005，为本研究中最低。

3.5 Yelloweye Rockfish

第一批次PCA显示Salish海个体与外海岸个体明显分离；第二批次仅GOA样本单独分析时无地理聚集。排除Salish海后，从南加州至GOA的遗传分化极低，成对F_ST为0.007–0.028，而Salish海与外海岸比较为0.024–0.028。

3.6 Black Rockfish

PCA中聚类间距离与地理距离一致，海洋距离与遗传距离相关分析R²=0.88支持IBD；此外阿拉斯加东南部样本在PCA中分为两个类群，染色体水平和基因组F_ST扫描显示两群间分化遍布整个基因组；最大平均F_ST为西海岸与西部GOA之间的0.065。

3.7 Dark Rockfish

全基因组PCA揭示两个先前未识别的遗传群，并在西部GOA同域分布；两个个体表现居中，admixture显示约50/50祖先比例。染色体PCA显示除一条染色体外，所有染色体均在两群间分离；基因组平均F_ST=0.012。

3.8 Northern Rockfish

初始PCA受三个性别相关染色体（chr20、22、24）影响；去除后PC1分离出两个与采样地点无关的遗传群，PC2存在批次效应；admixture一致支持两个群；染色体PCA显示所有染色体均分离；平均F_ST=0.019。

3.9 Blackspotted Rockfish

PCA从俄勒冈至白令海无聚集，admixture支持K=1；平均F_ST为0.010–0.017；这是在无可检测结构背景下唯一出现显著F_ST峰的物种，峰位于chr14约30 kb，最高F_ST=0.34；LD升高、局部PCA和等位基因剂量显示该峰与地理相关，白令海和阿留申群岛个体中主要等位基因剂量较低，而西海岸、阿拉斯加东南部和GOA个体大多为纯合主要型。

讨论总结

讨论指出，浮游幼虫持续时间（pelagic larval duration, PLD）与基因流的关系并非简单对应，例如Shortspine Thornyhead具长扩散期且无结构，Black Rockfish PLD较短且存在IBD。物理屏障可限制扩散，如Yelloweye Rockfish在Salish海的分化，以及Black Rockfish在阿拉斯加东南部的两个类群可能涉及内水与外太平洋的有效隔离。隐存多样性在岩鱼中可能比以往认识更常见；Dark和Northern rockfishes中隐存群同域分布，但分布比例沿地理空间变化。个体水平分析可揭示按采集地点估计频率所掩盖的结构，例如Northern Rockfish先前的IBD信号实际由两组比例变化造成。lcWGS存在批次效应和参考基因组偏差等限制。管理方面，泛交不代表不存在生态学上重要的种群；隐存多样性可能需作为独立管理单元；Blackspotted Rockfish的适应性离群峰可作为渔区分界的参考。

研究结论

商业化捕捞海洋物种基因组数据质量的持续提高，为表征物种多样性、细化种群信息并用于管理边界提供了机会。就太平洋岩鱼而言，全基因组数据和个体水平分析揭示了隐存多样性以及近似泛交的结构，并提示栖息地、深度、海洋环流和物种总体地理范围可能与岩鱼更广泛的分化模式存在联系。虽然本研究的9个物种和约1000份样本的测序数据已是探索多种遗传结构形成机制的重要数据集，但未来项目宜扩展分析以涵盖更多生活史特征（如潮间带和海藻林物种）、更广地理范围（如西太平洋和北大西洋物种），并可能显式检验物理和环境变量。