《Aquaculture Reports》：Development and validation of a 40?K Liquid SNP array for ark shell (
Scapharca subcrenata)
编辑推荐：
魁蚶（Scapharca subcrenata）是一种具有重要商业价值的海产双壳贝类，广泛分布于中国南北沿海。近年来，该物种面临种质资源衰退、潜在群体混杂以及育种进展缓慢等问题。SNP芯片因其大规模分析能力、成本效益高和精度优异而成为基因分型领域的基石技术，在
魁蚶（Scapharca subcrenata）是一种具有重要商业价值的海产双壳贝类，广泛分布于中国南北沿海。近年来，该物种面临种质资源衰退、潜在群体混杂以及育种进展缓慢等问题。SNP芯片因其大规模分析能力、成本效益高和精度优异而成为基因分型领域的基石技术，在现代基因组学中不可或缺。为推进该物种的分子育种计划，研究人员首次开发并验证了一种名为“Hanxin No.?Ⅰ”的魁蚶40?K SNP基因分型芯片。研究人员从中国五个野生群体中共采集520个个体，经过对全基因组重测序数据的严格过滤后，获得了5,684,437个高置信度SNP。随后，研究人员筛选出一组在全基因组范围内均匀覆盖的40,558个SNP用于芯片开发，其中包括采用宽松GWAS阈值选出的898个与推定生长相关性状相关的标记，以及1,895个与群体分化相关的标记。该SNP芯片展现出稳健的基因分型能力，在65个魁蚶样本中总体检出率达到97.68%，其中次要等位基因频率（MAF）>?0.1的位点占总数的94.77%。此外，为验证该芯片在群体遗传学中的有效性和准确性，研究人员检测了八个不同群体的遗传多样性和群体结构。研究结果表明，北方与南方群体之间存在明显的遗传分层。综上所述，该SNP芯片为魁蚶的基因分型提供了一个可靠的高分辨率平台，显著推进了种质资源保护的进程，并可能为未来的基因组选择或育种研究提供有用的标记资源。
一、研究背景与科学问题
魁蚶（Scapharca subcrenata）是中国、韩国和日本沿海潮间带广泛分布的重要经济贝类，其肉质鲜美且富含蛋白质和维生素，具有较高的药用价值。然而，近几十年来，由于过度捕捞、围海造陆及海洋环境退化，魁蚶的野生资源显著下降。同时，随着水产养殖规模的迅速扩张及跨省苗种频繁调运，不同地理群体之间可能出现混杂，导致遗传多样性降低的风险升高。尽管此前已有关于魁蚶基因组、转录组及代谢组等多组学研究基础，但其育种工作仍主要依赖野生采捕亲本和群体选择，新品种培育速度受到严重制约。因此，亟需开发高效的SNP标记和高通量基因分型工具，以加速魁蚶的全基因组育种进程。
二、研究方法与技术路线
研究人员从中国五个沿海地区（浙江宁海、浙江三门、辽宁大连、山东即墨、河北唐山）采集520个野生魁蚶个体，测量其壳长、壳宽、壳高和体重四个生长性状，提取肌肉组织DNA后进行全基因组重测序（Whole Genome Sequencing, WGS），平均测序深度为7.09?×。利用标准化生物信息流程进行质量控制和变异检测，经严格过滤获得5,684,437个高质量SNP。在此基础上，研究人员采用GenoBaits液相芯片技术，筛选出40,558个SNP位点开发40?K SNP芯片“Hanxin No.?Ⅰ”。芯片包含四类位点：第1层为采用宽松GWAS阈值（-log10P > 5）选出的898个推定生长相关标记；第2层为1,895个群体分化鉴定标记；第3层为34,270个调控元件相关标记；第4层为3,495个补充覆盖空缺的标记。随后，研究人员从8个不同群体（包括韩国群体）随机选取65个个体进行芯片基因分型验证，并通过重复样本和WGS比对评估芯片的一致性与准确性。
三、研究结果
3.1 全基因组重测序与SNP识别
研究人员对520个野生魁蚡进行全基因组重测序，共获得4,717.87?Gb原始数据，平均测序深度7.09?×，平均Q30达91.54%。经变异检测和严格过滤，最终获得5,684,437个高置信度SNP，均匀分布于所有19条染色体上，染色体SNP密度为0.14–18.58 SNPs/kb。
3.2 基于WGS的主成分分析
主成分分析（Principal Component Analysis, PCA）将520个个体清晰划分为两个遗传类群，PC1和PC2分别解释5.75%和0.32%的遗传方差。连锁不平衡（Linkage Disequilibrium, LD）衰减分析显示即墨群体（JM）具有较快的LD衰减速率，表明其遗传变异性较高。
3.3 生长性状的全基因组关联分析
研究人员对壳长、壳宽、壳高和体重四个生长性状进行全基因组关联分析（Genome-Wide Association Study, GWAS）。壳长检测到2个显著关联SNP和9个暗示性关联SNP；壳宽检测到11个暗示性关联SNP；壳高检测到1个显著关联SNP和9个暗示性关联SNP；体重检测到2个显著关联SNP和22个暗示性关联SNP。位于13号染色体14,669,077?bp处的SNP是四个性状共同的暗示性关联位点。
3.4 “Hanxin No.?Ⅰ”芯片特征分析
该芯片包含40,558个SNP位点，均匀覆盖19条染色体。94.77%的SNP次要等位基因频率（Minor Allele Frequency, MAF）大于0.1。核苷酸替换类型中转换（A/G和C/T）最为常见。功能分类包括2,018个非同义突变、6,776个同义突变、10个终止获得和1个终止丢失突变，以及大量内含子区和基因间区变异。
3.5 芯片基因分型性能评估
65个样本靶向捕获测序平均获得12,879,953.65条reads，平均比对率为96.68%。芯片平均检出率达97.68%。批次间重复样本基因型一致性为96.18%，同批次三次重复一致性为96.71%。芯片与WGS的基因型一致率为93.63%。
3.6 群体遗传结构分析
基于39,617个SNP的群体结构分析显示，最优K值为2，北方群体（河北乐亭、辽宁东港、山东东营、山东乳山、山东潍坊、天津）聚为一支，浙江台州群体独立成一支，韩国群体表现出与北方群体的混合祖先成分。PCA和邻接法系统发育树均支持这一分组模式。
3.7 遗传多样性与分化分析
多态信息含量（Polymorphism Information Content, PIC）范围为0.184（山东潍坊）至0.235（韩国）。观察杂合度（Ho）范围为0.213（浙江台州）至0.352（韩国），期望杂合度（He）范围为0.230（山东潍坊）至0.291（韩国）。近交系数（Fis）范围为-0.165（韩国）至0.057（浙江台州），基因组近交系数（FROH）范围为0.003（韩国）至0.025（浙江台州）。群体间固定指数（Fst）分析显示，除浙江台州外，其余群体间Fst为0.007–0.06，属于低至中度分化；而浙江台州与其他群体间的Fst为0.155–0.223，表明存在显著遗传分化。
四、讨论与结论
讨论部分指出，高通量SNP芯片是水生生物遗传改良的关键工具。“Hanxin No.?Ⅰ”芯片基于520个来自中国五个代表性区域的野生个体的WGS数据设计，具有高多态性和均匀基因组覆盖。与同类芯片相比，其MAF > 0.1的位点比例（94.67%）优于泥蟹、鲈鱼和鲍鱼的相关芯片。芯片检出率（97.68%）与皱纹盘鲍、硬壳蛤、大黄鱼和海参等物种的芯片相当。芯片与WGS的一致性（93.63%）略低于部分报道，研究人员认为这可能与水生动物普遍较高的杂合度有关。群体遗传分析揭示了北方与南方群体之间的显著遗传分化，这与地理分布格局一致，而韩国群体与北方群体共享部分祖先成分，可能源于黄海两侧幼虫的频繁迁移。遗传多样性参数表明韩国群体多样性最高，山东潍坊群体最低，后者可能存在一定程度的近交。
研究结论：研究人员首次开发了魁蚶40?K液相SNP芯片“Hanxin No.?Ⅰ”。该芯片包含40,558个均匀分布于19条染色体的SNP位点，具有高位点多态性、可接受的基因分型性能，以及在种质鉴别和未来育种研究中的潜在应用价值。随着SNP芯片的持续改进，它可能为未来的基因组选择或育种研究提供有用的标记资源。
