《Aquaculture Reports》：Ploidy-associated glycogen differences in
Crassostrea gigas: Insights from the glycogen-metabolism effector layer
编辑推荐：
摘要：倍性变异影响牡蛎的生长、繁殖和营养性状，但其与糖原积累的关系仍知之甚少。本研究以二倍体、三倍体和四倍体太平洋牡蛎（Crassostrea gigas）为材料，采用糖原含量测定、糖原合酶活性检测、mRNA转录组学、miRNA测序及lncRNA profil
摘要：倍性变异影响牡蛎的生长、繁殖和营养性状，但其与糖原积累的关系仍知之甚少。本研究以二倍体、三倍体和四倍体太平洋牡蛎（Crassostrea gigas）为材料，采用糖原含量测定、糖原合酶活性检测、mRNA转录组学、miRNA测序及lncRNA profiling进行了比较研究。在三个倍性组中，三倍体牡蛎表现出最高的糖原含量和糖原合酶活性，代表了最显著的高糖原表型。该表型伴随基因表达及通路富集的显著差异，涉及碳代谢、糖酵解/糖异生、淀粉和蔗糖代谢以及糖原代谢。趋势分析鉴定出与倍性依赖性糖原变异相关的表达模式。关键糖原代谢相关基因，包括Hex-t2、Glys、GBE1和PYGM，在不同倍性组间表现出不同的表达模式。此外，差异表达miRNA显示出与糖原代谢相关基因的潜在调控关系，而lncRNA相关的ceRNA网络提示了在碳水化合物代谢中的额外调控参与。综上，这些发现表明三倍体牡蛎中升高的糖原含量与糖原代谢效应层的协调变化相关，其转录模式提示糖原合成增强且糖原动员可能改变。本研究为牡蛎中倍性相关糖原调控的分子基础提供了见解，并为种质改良和品质育种提供了理论参考。
论文解读：倍性相关糖原差异的分子机制——以太平洋牡蛎为模型
研究背景与科学问题
糖原作为葡萄糖的高度分支聚合物，是动物体内主要的能量储存多糖。在贝类中，糖原不仅是重要的碳水化合物储备，还参与生长、繁殖及环境胁迫适应。牡蛎的糖原含量与组织质地、风味形成、色泽及整体营养品质密切相关，因此是兼具生理意义和食品品质价值的重要性状。中国是牡蛎养殖大国，年产量超过600万吨，其中太平洋牡蛎（Crassostrea gigas）是经济价值最高、养殖最广的物种之一。已有研究表明，牡蛎糖原含量受营养状态、环境条件和生殖发育等多种因素影响，并在个体和群体间表现出稳定而显著的差异。这种可重复的糖原相关表型差异为探索糖原积累的分子机制提供了理想模型。
多倍体育种已成为贝类养殖中改良产量和品质的重要策略。三倍体牡蛎因快速生长、性腺发育延迟和生殖投入减少，通常在营养积累尤其是糖原储存方面具有明显优势。然而，倍性变异影响糖原积累的分子机制远未被完全理解，特别是直接参与糖原合成与分解的关键效应基因是否发生系统性表达变化，以及这些变化是否受到非编码RNA的协同调控，仍不清楚。尽管已有研究关注不同倍性牡蛎在能量代谢上的广泛改变，但对糖原代谢效应层及其上游非编码RNA和信号通路调控的系统分析仍不充分。
研究内容与主要结论
本研究以二倍体（Cg_2N）、三倍体（Cg_3N）和四倍体（Cg_4N）太平洋牡蛎为实验材料，整合糖原含量测定、糖原合酶活性检测、mRNA转录组、miRNA测序和lncRNA分析，系统探究不同倍性水平下糖原积累的机制。结果显示，三倍体牡蛎表现出最高的糖原含量（平均值144.09 mg/g），分别是二倍体（21.08 mg/g）的6.83倍和四倍体（38.33 mg/g）的3.76倍；四倍体糖原含量为二倍体的1.82倍，差异均极显著（p<0.001）。糖原合酶活性也呈现相同趋势：三倍体（708.45 U/g）>四倍体（310.67 U/g）>二倍体（202.31 U/g），差异极显著（p<0.001）。这一结果表明三倍体牡蛎具有最明显的高糖原表型，且该表型与糖原合成能力增强密切相关。
转录组分析发现，三倍体与二倍体及四倍体之间存在大量差异表达基因（DEGs），分别为7893个和9696个，而二倍体与四倍体之间仅682个，表明三倍体经历了更广泛的转录组重编程。GO和KEGG富集分析显示，差异基因主要富集于碳代谢、糖酵解/糖异生、淀粉和蔗糖代谢、氨基酸和核苷酸糖代谢等与糖原前体供给相关的通路。趋势分析进一步鉴定出两个与高糖原表型明确对应的共表达模块：profile 5在三倍体中峰值表达，富集于半乳糖代谢、果糖和甘露糖代谢及磷酸戊糖途径，提示有利于糖原合成的底物供应；profile 2则在二倍体和四倍体中高表达而三倍体中显著降低，主要富集于糖酵解/糖异生、丙酮酸代谢、柠檬酸循环（TCA循环）和氧化磷酸化，提示碳水化合物氧化利用相对减弱。这一“合成增强-分解减弱”的定向模式构成了高糖原表型的核心分子基础。
在糖原代谢效应层，关键基因表达差异显著：葡萄糖利用入口基因Hex-t2上调，直接参与糖原合成的Glys（糖原合酶）和GBE1（糖原分支酶）上调，而糖原动员相关基因PYGM（糖原磷酸化酶）下调。同时，负调控糖原合酶的GSK3β下调，糖原起始蛋白GYG1及糖原结合蛋白Gbs-76A表达呈上升趋势。这些变化共同指向碳流更稳定地导向糖原储存。
关键技术方法
研究人员使用了以下主要技术方法：①流式细胞术鉴定倍性；②糖原含量测定采用碱水解法结合蒽酮显色法；③糖原合酶活性测定使用商品化试剂盒；④mRNA转录组测序基于去rRNA文库构建，利用HISAT2比对至长牡蛎参考基因组（GCA_011032805.1），DESeq2鉴定差异表达基因；⑤miRNA测序通过small RNA文库构建，miRDeep2预测新miRNA，edgeR分析差异表达；⑥lncRNA测序与mRNA共用文库，StringTie重构转录本，Gffcompare区分已知和新型lncRNA；⑦靶基因预测联合RNAhybrid、miRanda和TargetScan，ceRNA网络通过hypergeometric检验构建；⑧qRT-PCR验证。样本来源于山东烟台崆峒岛养殖基地，二倍体亲本为审定品种“鲁益1号”（批准号GS-01-006-2020），三倍体由二倍体雌性×四倍体F1雄性杂交获得，四倍体为F2代。所有牡蛎在自然浅海吊养两年，实验取性腺组织进行分析。
结果分述
1. 倍性鉴定与糖原表型差异：流式细胞术确认三种倍性材料准确；三倍体糖原含量最高，糖原合酶活性也最高，为高糖原表型提供酶学证据。
2. 转录组与糖原相关差异通路：三倍体与另外两组相比DEGs数量多，富集通路涉及碳代谢、糖酵解/糖异生等，关键糖原代谢基因存在方向性差异，表明高糖原表型源于系统性碳流分配改变。
3. 趋势分析：将DEGs分为8种表达模式，其中profile 5（三倍体高表达）和profile 2（三倍体低表达）与糖原表型对应，功能富集显示二者分别偏向糖原合成底物供给和碳水化合物氧化利用。
4. miRNA调控：差异表达miRNA（Cg_2N_vs_Cg_3N:184个，Cg_2N_vs_Cg_4N:96个，Cg_3N_vs_Cg_4N:229个）的靶基因富集于碳代谢和糖原相关通路，预测到与Glys、GBE1、Hex-t2、PYGM的miRNA-mRNA配对关系，提示miRNA可能参与糖原代谢基因的转录后调控。
5. lncRNA与ceRNA网络：差异表达lncRNA数量在包含三倍体的比较中显著较多（Cg_2N_vs_Cg_3N:760个，Cg_3N_vs_Cg_4N:1307个），靶基因功能富集同样涉及碳代谢。以Hex-t2为中心的ceRNA网络显示MSTRG.21754.1、MSTRG.4965.1和MSTRG.5934.1等lncRNA可能在糖原积累中发挥调控作用，网络中的mRNA多与代谢稳态、应激反应和信号转导相关，提示lncRNA可能通过塑造有利于糖原储存的细胞环境来维持高糖原状态。
6. qRT-PCR验证：代表性mRNA、miRNA和lncRNA的表达模式与测序结果一致，验证了数据的可靠性。
讨论总结
研究人员提出，倍性差异导致的糖原含量变化不能简单概括为“能量代谢重编程”，而应更准确地描述为倍性驱动下糖原代谢效应层的定向转变：三倍体中糖原合成相关基因表达增强且糖原分解相关基因表达相对减弱，同时糖原合酶活性升高，共同促成高糖原表型。miRNA可能作为转录后水平的“稳定器”，通过差异调控维持效应层基因的稳态表达；lncRNA相关ceRNA网络则更可能通过调节代谢环境、应激反应和信号转导背景，为高糖原状态提供上游支持。
结论翻译
本研究在表型、糖原合酶活性、mRNA转录组、miRNA及lncRNA水平上系统比较了二倍体、三倍体和四倍体太平洋牡蛎的糖原差异。结果表明，三倍体牡蛎表现出最显著的高糖原表型，这更好地由糖原代谢效应层的定向转变来解释。一方面，葡萄糖输入、糖原合成相关基因表达和糖原合酶活性增强；另一方面，与糖原动员和分解相关的过程相对减弱，从而将整体碳流更强烈地导向储存。同时，miRNA可能在转录后水平维持糖原代谢关键节点的稳定表达，而lncRNA相关ceRNA网络更可能通过调节代谢稳态、应激反应和信号转导背景来为高糖原状态提供上游支持。这些发现为理解牡蛎中倍性相关糖原积累的分子机制提供了新的解释框架，也为未来高品质牡蛎育种和多倍体种质资源的高效利用提供了理论基础。该论文发表于《Aquaculture Reports》。
