基于基因组学方法区分受威胁石首鱼类加州湾石首鱼（Totoaba macdonaldi）合法圈养繁育与非法野生捕捞

《Aquaculture Reports》：Genomic approach to distinguish between legal captive breeding and illegal wild fishing of threatened totoaba (Totoaba macdonaldi)

　　

遗传多样性的理解对于保护策略与可持续水产养殖实践，以及产品合法来源的判定至关重要，尤其是对受威胁物种而言，例如加州湾石首鱼（Totoaba macdonaldi）这一加利福尼亚湾特有的大型鱼类。研究共对170尾个体进行了双酶切限制性位点相关DNA测序（double-digest restriction-site associated DNA, ddRAD），其中92尾为野生个体，78尾为来自两个野生动物保护与管理单元（Units for Conservation and Management of Wildlife, UMAs：CREMES和UABC）的圈养亲本，最终获得4992个中性单核苷酸多态性位点（single nucleotide polymorphisms, SNPs）。野生采样地点之间未检测到显著遗传分化，表明在本研究所分析的空间尺度上未发现可检测的种群结构。相较之下，野生种群与UMA圈养亲本之间，以及两个UMA之间，均存在显著差异。野生种群内部未检测到亲缘关系，而UMA的F₁与F₂亲本中检测到显著水平的亲缘性。种群遗传学方法能够将全部个体划分至各自来源。以野生种群与合并后的UMAs为来源进行判定时，来源归属准确率总体均值较高，>88%（WILD：不同位点子集为94–100%；UMAs：84–88%）；以野生种群与各UMA亲本分别比较时，准确率>91%（CREMES约74–80%，UABC约90%，WILD为100%）。在控制亲缘关系后，两类分析的总体平均准确率分别略降至>86%和>90%。本研究凸显了SNP在加州湾石首鱼遗传多样性监测中的有效性，并补充了个体来源精准判定所需的信息，这对于合法可追溯性与执法监管具有重要意义。

该文发表于《Aquaculture Reports》，围绕受威胁鱼类加州湾石首鱼（Totoaba macdonaldi）合法养殖来源与非法野生捕捞来源的分子鉴别展开。研究背景在于，保护遗传学（conservation genetics）已成为濒危物种保护、资源管理与可持续利用的重要基础工具。对于兼具高经济价值与严格法律监管属性的物种，野生个体与圈养个体在法律地位上存在显著差异，因此需要可靠的来源追溯体系，以支撑食品安全、执法监管以及对非法、未报告和无管制捕捞（illegal, unreported and unregulated fisheries, IUU）的打击。加州湾石首鱼是加利福尼亚湾特有的大型石首鱼科（Sciaenidae）鱼类，曾是当地极为重要的渔业资源，但因资源急剧下降，自1975年起被永久禁捕，并先后受到《濒危野生动植物种国际贸易公约》（CITES）附录Ⅰ以及美国《濒危物种法》等法规保护。尽管其保护等级已有调整，但鱼鳔高额国际需求持续推动偷捕与非法贸易，使“合法养殖产品”与“非法野生来源产品”的区分成为现实而迫切的治理难题。

在此背景下，研究人员以基因组学（genomics）方法建立野生种群与圈养种群的遗传基线，旨在回答两个核心问题：其一，加州湾石首鱼野生群体在研究空间尺度上是否呈随机交配（panmixia）；其二，能否借助单核苷酸多态性（SNP）标记实现个体来源的高准确率判定，从而服务于可追溯性与执法。文章指出，虽然现有管理体系已采用微卫星（microsatellite）标记开展亲本鉴定和群体研究，但SNP具有更高的实验室间可重复性，也更适合在渔业与水产养殖管理中作为补充或升级工具。因此，本研究的意义不仅在于阐明种群遗传结构，也在于为监管部门和养殖管理者提供可操作的遗传识别依据。

在方法上，研究人员共分析170份样本，其中野生样本92份，来自加利福尼亚湾5个采样点；圈养样本78份，来自两个UMA机构，即CREMES与UABC，并覆盖F₀、F₁和F₂世代。研究采用双酶切限制性位点相关DNA测序（ddRAD）获取全基因组分布的SNP位点，通过STACKS进行de novo组装与分型，在缺乏参考基因组的条件下构建位点目录；随后经离群位点筛查、连锁不平衡过滤、哈代-温伯格平衡（Hardy-Weinberg equilibrium, HWE）检验及缺失数据控制，保留4992个中性SNP。之后，研究人员估计遗传多样性参数、有效群体大小（effective population size, N_e）、种群分化与亲缘关系，并利用判别主成分分析（DAPC）、STRUCTURE聚类及基于支持向量机（support vector machine, SVM）的assignPOP框架开展来源归属分析。样本来源明确，包括2010—2013年的野生群体与2021—2022年的UMA圈养亲本。

研究结果部分可按原文小标题归纳如下。

3.1. Individuals and SNP filtering
研究人员通过两条测序lane共获得1,733,678,472条原始reads，经质量过滤后保留1,492,680,881条。de novo组装初步检测到226,807个RAD位点和215,565个SNP，经群体过滤后形成5691个SNP的全局数据集。进一步剔除低测序深度、高缺失率样本、存在连锁不平衡位点、显著偏离HWE位点及疑似交叉污染相关个体后，最终建立包含145个个体、4992个中性SNP的高质量数据集。这一结果说明研究建立了可用于后续种群遗传分析和来源判定的稳健分子标记基础。

3.2. Genetic diversity
遗传多样性分析显示，野生种群的等位基因丰富度（allelic richness, Ar）最高，为1.96，高于两个圈养群体，提示野生群体保留了更丰富的遗传变异。观察杂合度（H_O）和期望杂合度（H_E）总体处于中等水平，近交系数（F_IS）在野生与圈养群体中均为负值，表明未观察到明显的近交信号。UABC与野生群体和CREMES在H_O与H_E方面存在显著差异。多数SNP位点具有较好的多态信息含量（polymorphic information content, PIC），其中36.2%的位点PIC > 0.30，说明该SNP数据集在群体区分和来源归属方面具有较强信息性。总体而言，野生加州湾石首鱼尽管面临长期捕捞压力，仍保有相对较高的遗传多样性，而圈养群体的多样性略低。

3.3. Estimates of effective population size (N_e)
有效群体大小估计显示，野生群体的N_e远高于圈养群体。采用连锁不平衡法（LD-N_e）时，野生群体估值为9272.2，置信区间较宽；采用杂合子过剩法（He-N_e）时，野生群体结果为无限值。相比之下，UABC和CREMES的圈养群体N_e明显偏低。该结果表明，野生种群在当前样本框架下仍表现出较大的有效群体规模，而圈养亲本群体由于繁殖个体数量有限和家系结构明显，其有效群体规模较小，存在遗传漂变（genetic drift）和近交风险增加的管理隐患。

3.4. Population differentiation
3.4.1. Wild totoaba locations
针对5个野生采样点的群体结构分析中，成对F_ST值非常低且均不显著，DAPC散点图中各地点高度重叠，STRUCTURE分析亦未显示清晰亚群划分。研究人员据此认为，在本研究评估的空间尺度上，野生加州湾石首鱼可视为单一遗传单元，即支持随机交配假说。这一结论意味着加利福尼亚湾范围内存在持续基因流，与该物种迁移、产卵与索饵生活史模式相一致。

3.4.2. Wild totoaba population vs UMAs captive broodstock
与野生地点之间缺乏分化不同，野生群体与两个UMA圈养亲本之间存在显著遗传分化。成对F_ST比较显示，UABC与CREMES之间分化最大，UABC与野生群体之间次之，CREMES与野生群体之间分化相对较低。DAPC识别出与野生、CREMES和UABC相对应的3个遗传簇，STRUCTURE在K = 3时也呈现三群分化格局。两家UMA中的F₀个体因源自野外，与野生群体存在遗传混合信号。亲缘关系分析进一步表明，野生群体内部无显著亲缘关系，而CREMES和UABC中大量F₁与F₂个体表现出一级和二级亲缘关系。这说明圈养群体的独立管理历史和家系结构共同塑造了可识别的基因组差异。

3.5. Population assignment
3.5.1. Wild population vs global UMAs captive broodstock
当将两个UMA合并为一个圈养来源组，与野生组进行比较时，来源归属准确率整体较高，在不同高F_ST位点子集和训练集比例下总体均超过89%。野生个体的归属准确率高于圈养个体，最高可达100%；圈养个体准确率为84%至88%。在控制亲缘关系、即每个家系仅保留一个代表用于模型训练后，总体准确率仍维持在86%以上，但圈养群体准确率有所下降且波动增大。这说明家系结构会对圈养组内部判定产生一定影响，但并未削弱野生与圈养总体区分能力。

3.5.2. Wild population vs individual UMAs captive broodstock
当分别将CREMES或UABC与野生群体进行比较时，整体归属准确率进一步提高，均超过91%。其中，CREMES的准确率约为74%至80%，UABC约为90%，野生群体在所有设置下均达到100%。控制亲缘关系后，总体准确率仍高于90%，但两个UMA内部的准确率均有所下降，尤其在样本量进一步减少时更为明显。这表明若将不同UMA分别建模，来源判定性能更优，同时也反映出样本量和家系非独立性对圈养组预测稳定性具有实质影响。

3.5.3. Assignment probabilities
对于独立设置的未知来源个体，模型给出的归属概率整体较高。合并UMA模型中，UMA个体平均归属概率为99%，野生个体为88%；UMA分组模型中，CREMES为95%，UABC为96%，野生个体为85%。在控制亲缘关系后，这些数值变化很小，说明所构建的SNP数据集和分类模型在真实来源判定场景中具有较强稳健性。研究据此认为，这套标记系统具备发展为精简SNP面板的潜力。

讨论部分集中强调了四方面内容。首先，野生加州湾石首鱼在研究所用历史样本中表现出相对较高的遗传多样性和较大的有效群体规模，这与其濒危地位形成一定反差，但与既往研究结果一致，说明该物种尚未出现明显的遗传衰退信号。其次，野生群体在所检测空间尺度上呈随机交配状态，支持其为单一空间遗传单元的观点。再次，两个UMA圈养群体因家系结构与独立管理而形成清晰的基因组分化，这种分化恰恰为来源追踪提供了判别基础。最后，基于SNP的归属方法在区分野生来源与圈养来源方面表现出高准确率，且在控制亲缘关系后总体性能仍较稳健，因此对合法贸易监管、打击非法捕捞和未来种群遗传监测具有直接应用价值。文章同时指出，若未来扩大圈养群体样本规模，并进一步开发精简SNP面板，来源追踪工具的实用性还可继续提升。

研究结论部分可译为：本研究利用RAD测序产生的SNP标记，揭示了野生与圈养加州湾石首鱼亲本群体的遗传多样性和种群结构。与通常对濒危物种的预期不同，野生Totoaba macdonaldi种群表现出相对较高的遗传多样性，这与既有研究一致，表明尽管经历过度捕捞和栖息地改变，该物种仍维持了显著遗传变异。基于现有野生样本的分析未发现采样地点之间存在显著空间遗传分化，支持在本研究评估尺度上的随机交配假说。相对地，SNP标记能够准确区分野生群体与UMA圈养亲本之间的基因组结构，并对来源种群实现较高精度归属；据此可进一步开发精简SNP面板，以有效识别来自野生、UABC和CREMES的个体。这对于加州湾石首鱼个体合法或非法来源的判定、可追溯性建设及执法监管具有高度重要性。未来工作应聚焦于系统性的基因组监测与管理策略优化，以确保野生群体和圈养亲本群体的长期可持续性。