《Fishes》：Detecting Natural Reproduction Following Reintroduction Using Complementary Genetic Approaches in the Endangered Zingel asper
编辑推荐：
拟穴青蟹（Scylla paramamosain）第五步足（游泳足）的浆状趾节是一种特化的附肢，但其性状的遗传基础和基因组预测潜力尚不清楚。研究人员利用323个个体的全基因组重测序数据，对双侧趾节长度、宽度、面积及第一主成分（PC1）进行了全基因组关联分析（G
拟穴青蟹（Scylla paramamosain）第五步足（游泳足）的浆状趾节是一种特化的附肢，但其性状的遗传基础和基因组预测潜力尚不清楚。研究人员利用323个个体的全基因组重测序数据，对双侧趾节长度、宽度、面积及第一主成分（PC1）进行了全基因组关联分析（GWAS）和基因组预测。遗传力估计值范围为0.184至0.538。GWAS鉴定出与左侧和右侧趾节长度相关的显著SNP分别为116个和70个，与对应宽度性状相关的分别为73个和68个，与对应面积性状相关的分别为102个和76个，以及与PC1相关的91个。主要信号出现在长度性状的19号和28号染色体，宽度性状的19号和23号染色体，以及复合性状的19号、21号和30号染色体。两个SNP（Chr16:23074328和Chr30:357279）在所有七个表型中共享，并对双侧长度和宽度表现出一致的加性效应。Chr30:357279被注释为SRSF6-like基因中的错义变异。基于GWAS优先SNP子集的GBLUP模型比全基因组SNP提高了预测准确性，其中FDR=0.05子集表现最佳（0.39–0.47）。这些结果为了解拟穴青蟹功能形态性状的遗传解析和基因组预测提供了新见解。
**论文解读：拟穴青蟹游泳足趾节性状的遗传基础与基因组预测研究**
**1. 研究背景与问题**
拟穴青蟹（Scylla paramamosain）是东南亚和中国重要的海水养殖经济物种，其第五步足（游泳足）的趾节特化为浆状结构，在游泳和摄食活动中发挥关键功能。然而，目前对拟穴青蟹的遗传育种研究主要集中于体重、甲壳宽等常规生长性状，而对游泳足趾节这类功能特化附肢的遗传基础知之甚少。尽管已有基于家系的研究利用全基因组重测序构建高密度连锁图谱并定位了相关数量性状位点（QTL），但群体水平的关联模式和基因组预测潜力尚未系统评估。因此，亟需通过群体水平的全基因组关联分析（GWAS）和基因组选择（GS）方法，揭示游泳足趾节形态性状的遗传架构，并评估其分子育种价值。
**2. 研究内容与结论**
本研究利用323个拟穴青蟹个体的全基因组重测序数据（平均深度9×），对游泳足趾节的长度、宽度、面积及第一主成分（PC1）进行了表型测量和遗传分析。通过GWAS鉴定出大量显著关联的单核苷酸多态性（SNP）和候选基因，并基于不同显著性阈值的SNP子集构建了基因组最佳线性无偏预测（GBLUP）模型，评估了基因组预测准确性。主要结论包括：游泳足趾节性状为低至中等遗传力的复杂数量性状（遗传力0.184–0.538），不同表型维度间存在共享和特异的遗传信号；鉴定出两个跨所有表型共享的SNP（Chr16:23074328和Chr30:357279），其中Chr30:357279位于SRSF6-like基因外显子中，为错义变异；基于GWAS优先SNP子集的GBLUP模型（尤其是FDR=0.05子集）在预测准确性上显著优于全基因组SNP模型（0.39–0.47 vs 0.02–0.10）。该研究为拟穴青蟹功能形态性状的遗传解析和基因组选择提供了新证据，论文发表在《Fishes》。
**3. 主要关键技术方法**
- **全基因组重测序与变异检测**：对323个拟穴青蟹个体（来自中国福建省福州市养殖场，178雄145雌，养殖周期180±7天）的肌肉DNA进行Illumina HiSeq X Plus平台测序，平均深度9×，经质控、比对至参考基因组ASM3559412v1、SNP calling和过滤后，获得4,042,299个高质量SNP。
- **全基因组关联分析（GWAS）**：采用线性混合模型（LMM），以PLINK提取的前六个主成分、性别、采样日期和体重作为固定效应，遗传相关矩阵作为随机效应，在GEMMA软件中实现，以p=1.0×10
-5为建议显著性阈值。
- **候选基因识别**：利用R包GALLO提取显著SNP上下游50 kb范围内的蛋白编码基因，并用EggNOG-mapper进行功能注释。
- **基因组最佳线性无偏预测（GBLUP）**：基于GCTA软件构建基因组关系矩阵（GRM），采用限制性最大似然（REML）估计遗传力，并构建GBLUP模型进行10折交叉验证（10次重复），预测准确性为预测值与观测值之间的Pearson相关系数。
**4. 研究结果**
**4.1 表型统计**
游泳足趾节长度、宽度和面积性状均近似正态分布，适合数量遗传分析。长度和宽度性状的变异系数（CV）为14.75%–15.27%，面积性状的CV为29.26%–29.77%。Pearson相关性分析显示面积与长度、宽度呈显著正相关（r=0.933–0.991，p<0.001）。SNP遗传力（h
2）估计值：左侧第5步足趾节长度（5PLL）为0.457，右侧（5PLR）为0.538；左侧宽度（5PWL）为0.213，右侧（5PWR）为0.184；左侧面积（SL）为0.358，右侧（SR）为0.382。
**4.2 长度和宽度性状的GWAS**
对5PLL鉴定出116个显著SNP，主要富集于19、28和1号染色体；5PLR鉴定出70个显著SNP，富集于28、27和24号染色体；25个SNP在双侧长度性状间共享，以28号染色体密度最高。对5PWL鉴定出73个显著SNP，主要分布于19、15和16号染色体；5PWR鉴定出68个显著SNP，富集于21、9、23和30号染色体；13个SNP在双侧宽度性状间共享，以23号染色体富集密度最高。
**4.3 复合性状的GWAS**
左侧趾节面积（SL）和右侧面积（SR）分别鉴定出102和76个显著SNP，其中34个共享，主要富集于30、21和19号染色体。PC1鉴定出91个显著SNP，25个位点在三个复合性状间共享。
**4.4 共享候选基因分析**
从所有GWAS显著SNP上下游50 kb区域共提取出多个候选基因。双侧长度性状共享49个候选基因，包括GRIA1-like、GRIN2B-like、SUZ12A-like、MMP17-like、HMCN2-like等，涉及信号转导和细胞外基质组织。双侧宽度性状共享29个候选基因，包括CLIP1-like、PTGS2-like、HSDL1-like、CFAP251-like、POFUT2-like等，与细胞骨架组织和脂代谢相关。面积与复合性状共享73个候选基因，包括IRS1B-like、ACSBG2-like、MTOR-like、AM1159-like、MMP17-like、HMCN2-like等，涉及生长信号和脂肪酸代谢。六个候选基因在所有七个表型中共享：LOC135116013、BP10-like、SRSF6-like、LOC135116139、LRP2-like、TL5A-like，均位于Chr30:357279周围50 kb区域。Chr30:357279被注释为SRSF6-like基因的错义变异（丙氨酸→苏氨酸替换）。
**4.5 关键SNP的基因型效应**
两个共享SNP（Chr16:23074328和Chr30:357279）在所有四个基础性状中均表现出一致的加性效应：随着替代等位基因T的累积，趾节长度和宽度增加。
**4.6 基于GBLUP的基因组预测准确性**
基于全基因组SNP（4,042,299个）的GBLUP模型预测准确性较低（0.02–0.10）。基于GWAS显著SNP子集（p<1.0×10
-5，n=68–116）的模型准确性为0.25–0.39；FDR=0.01子集（n=571–670）为0.34–0.44；FDR=0.05子集（n=11,035–12,187）为0.39–0.47，表现最佳。其中5PLR准确性最高（0.47），SL最低（0.39）。
**5. 讨论与结论**
**讨论部分总结**：游泳足趾节性状的遗传力差异表明，尽管表型高度相关，但不同维度的遗传基础可能不完全共享。GWAS结果揭示了多个染色体区域涉及共享和特异遗传成分，一致支持其复杂遗传架构。共享候选基因的功能注释提示了细胞外基质重塑、细胞骨架动态、代谢调控、外骨骼结构及转录后调控等生物学过程的潜在参与。两个共享SNP（尤其是Chr30:357279的SRSF6-like错义变异）具有一致的加性效应，是优先待验证的候选位点。GWAS优先SNP子集显著提高了基因组预测准确性，且FDR=0.05子集实现了最佳平衡，表明在样本量有限时，优先信息标记可有效提升预测效率。
**结论部分翻译**：本研究利用全基因组重测序数据，研究了拟穴青蟹（S. paramamosain）第五步足浆状趾节长度、宽度和面积的遗传基础，并评估了基于GWAS优先SNP的GBLUP模型的基因组预测性能。结果显示，这些性状为低至中等遗传力的复杂数量性状，不同表型维度间存在共享和特异遗传成分。GWAS鉴定出多个显著关联位点和候选基因，其中两个共享SNP（Chr16:23074328和Chr30:357279）在四个基础性状中表现出一致的加性效应。值得注意的是，Chr30:357279被注释为丝氨酸/精氨酸富集剪接因子6-like基因编码区的错义变异，使该位点成为优先验证的候选。与全基因组SNP模型相比，基于GWAS优先SNP的GBLUP模型显著提高了预测准确性，以FDR=0.05阈值选择的SNP子集表现最佳。总体而言，本研究为拟穴青蟹游泳足趾节性状的遗传解析和基因组预测提供了新证据，并为这类功能形态性状的分子育种提供了参考。
