《Aquaculture Reports》：Genome-wide association study and genomic selection reveal the genetic architecture of glycogen content in autotetraploid
Crassostrea gigas
编辑推荐：
专业摘要翻译：
同源四倍体长牡蛎（autotetraploid Crassostrea gigas）是一种具有重要经济价值的海产软体动物，因其能够稳定产生三倍体长牡蛎而备受关注。糖原含量是影响长牡蛎风味和品质的关键营养性状，改良该性状已成为重要的育种和商业目标
专业摘要翻译：
同源四倍体长牡蛎（autotetraploid Crassostrea gigas）是一种具有重要经济价值的海产软体动物，因其能够稳定产生三倍体长牡蛎而备受关注。糖原含量是影响长牡蛎风味和品质的关键营养性状，改良该性状已成为重要的育种和商业目标。本研究首次在同源四倍体长牡蛎中开展了糖原含量的全基因组关联分析（Genome-wide association study, GWAS）和基因组选择（Genomic selection, GS）分析。共鉴定出4个显著SNP和18个提示性SNP，候选基因主要涉及与糖原代谢相关的生物学过程，包括信号转导、葡萄糖转运和蛋白质磷酸化。研究人员系统评估了不同SNP密度和标记选择策略对预测准确性的影响，采用六种预测模型。结果表明，在本研究群体及评估的群体间，5K SNP标记组合足以实现对糖原含量相对较高的预测准确性，且纳入显著SNP可进一步提高预测准确性。值得注意的是，包含显著SNP且在染色体上均匀分布的1K SNP组合表现出相对较高的预测性能。此外，RKHS（Reproducing Kernel Hilbert Space）模型在预测准确性和计算效率之间表现出良好的平衡，可能是同源四倍体长牡蛎糖原含量预测的最优模型。然而，GWAS-GS策略的预测稳定性可能受群体遗传结构影响，需要在额外的独立群体中进行验证。本研究为深入理解同源四倍体长牡蛎糖原含量的遗传和代谢机制奠定了基础，并为开发低密度SNP芯片以及实施高效、经济的高糖原含量同源四倍体长牡蛎育种计划提供了理论支持和实证依据。
论文解读文章：
一、研究背景与科学问题
长牡蛎（Crassostrea gigas）是全球水产养殖中具有重要经济和营养价值的海洋软体动物，中国是全球最大的牡蛎生产国，2022年养殖产量达620万吨。糖原作为一种葡萄糖多糖，在动物和真菌中广泛分布，在能量储存、氧化还原调节及维持内环境稳态等生理过程中发挥重要作用。在长牡蛎中，糖原参与生殖、抗逆性和性腺发育等生理过程，同时也是决定其营养价值和风味品质的关键性状。长牡蛎糖原含量呈现季节性变化规律：夏季产卵期因配子释放而显著下降，秋冬季节则大量积累，导致夏季收获的长牡蛎营养价值和风味品质下降，难以满足市场需求。
三倍体长牡蛎因育性降低，能将更多能量分配给生长发育，在维持较高肥满度的同时持续积累糖原，因此市场接受度不断提高。三倍体长牡蛎的生产主要通过两种途径：化学诱导法存在无法稳定获得三倍体及化学试剂残留的食品安全隐患；而二倍体与同源四倍体杂交产生的三倍体则100%为三倍体且无化学残留。鉴于糖原含量在长牡蛎中具有可遗传性（h2=0.29），培育高糖原含量的同源四倍体长牡蛎群体可显著改善商品三倍体长牡蛎的营养品质，从而推动长牡蛎养殖业发展。
尽管GWAS已应用于长牡蛎糖原含量相关SNP的鉴定，如Meng等（2019）对427个个体鉴定出168个显著相关SNP，Liu等（2019a）通过GWAS鉴定出PPP1R3B为与糖原含量显著相关的基因，但GS和GWAS尚未应用于同源四倍体长牡蛎糖原含量的遗传分析，这限制了糖原代谢相关基因和机制的鉴定，以及高糖原含量同源四倍体长牡蛎的高效育种。
二、研究内容与主要结论
本研究首次在同源四倍体长牡蛎中开展了糖原含量的GWAS和GS分析。研究共鉴定出4个显著SNP和18个提示性SNP，候选基因主要涉及信号转导、葡萄糖转运和蛋白质磷酸化等与糖原代谢相关的生物学过程。研究系统评估了不同SNP密度和标记选择策略对预测准确性的影响，采用六种预测模型进行比较分析。结果表明，5K SNP标记组合足以实现相对较高的预测准确性，纳入显著SNP可进一步提高预测准确性；包含显著SNP且均匀分布的1K SNP组合表现出相对较高的预测性能；RKHS模型在预测准确性和计算效率之间表现出良好平衡，可能是最优预测模型。该研究为深入理解同源四倍体长牡蛎糖原含量的遗传和代谢机制奠定了基础，为开发低密度SNP芯片和实施高效经济的高糖原含量育种计划提供了理论支持和实证依据。
三、关键技术方法
研究人员从山东烟台崆峒岛实业有限公司获取同源四倍体长牡蛎样本，共360个性腺成熟的个体，来自三个群体：Group 1（206个个体，2龄）为崆峒岛长牡蛎群体诱导的同源四倍体群内繁殖获得的F1代；Group 2（125个个体，1龄）为Group 1群内繁殖获得的F2代；Group 3（29个个体，1.5龄）为山东烟台长岛长牡蛎群体诱导的同源四倍体群内繁殖群体。Group 1是Group 2的亲本，Group 3为独立群体。采用基于长牡蛎基因组的"鲁东力新-4N-I"SNP芯片进行基因分型，包含41,270个高质量SNP，筛选MAF>0.05、检出率>90%的双等位SNP用于后续分析，基因型按等位基因剂量编码为0、1、2、3、4。糖原含量采用蒽酮比色法糖原测定试剂盒测量，每个样本测量三次取平均值。主要技术方法包括：基于sommer R包的双变量基因组混合模型分析遗传相关性；基于adegenet R包的PCA分析和Vcfpop群体结构分析；基于ldsep R包的连锁不平衡分析；基于GWASpoly R包的Q+K混合模型GWAS分析；基于SnpEff的SNP功能注释；基于DAVID的GO和KEGG富集分析；六种GS预测模型的比较评估。
四、研究结果
4.1 群体结构与连锁不平衡分析
PCA分析揭示了样本间的遗传分化模式，Vcfpop群体结构分析通过ΔK方法确定最优K值，kinship矩阵进一步阐明了个体间的潜在关系。LD分析计算了SNP对间的r2并可视化LD衰减。
4.2 GWAS鉴定显著SNP
采用GWASpoly R包中的Q+K混合模型进行GWAS分析，显著性阈值和提示性显著性阈值分别设为0.05/SNP数和1/SNP数。共鉴定出4个显著SNP和18个提示性SNP。SnpEff用于SNP注释，选取显著SNP上下游50kb范围内的基因作为候选基因，基于基因潜在功能鉴定出最可能与糖原代谢相关的候选基因。
4.3 候选基因功能富集分析
通过DAVID进行GO和KEGG富集分析，q值<0.05的GO条目和KEGG信号通路被认为显著富集。候选基因主要涉及信号转导、葡萄糖转运和蛋白质磷酸化等与糖原代谢相关的生物学过程。
4.4 基因组选择预测准确性评估
研究系统评估了不同SNP密度和标记选择策略对预测准确性的影响，采用六种预测模型进行比较。结果表明5K SNP组合足以实现相对较高的预测准确性，纳入显著SNP可进一步提高预测准确性。包含显著SNP且均匀分布的1K SNP组合表现出相对较高的预测性能。RKHS模型在预测准确性和计算效率之间表现出良好平衡，可能是同源四倍体长牡蛎糖原含量预测的最优模型。
五、讨论与结论
讨论部分指出，GWAS-GS策略的预测稳定性可能受群体遗传结构影响，需要在额外的独立群体中进行验证。研究结论表明，本研究为深入理解同源四倍体长牡蛎糖原含量的遗传和代谢机制奠定了基础，为开发低密度SNP芯片和实施高效、经济的高糖原含量同源四倍体长牡蛎育种计划提供了理论支持和实证依据，有助于推动长牡蛎养殖业的发展。
