二倍体、三倍体与四倍体太平洋牡蛎（Crassostrea gigas）硒诱导抗氧化响应的转录组学分析

《Foods》：Transcriptomic Analysis of Selenium-Induced Antioxidant Responses in Diploid, Triploid, and Tetraploid Pacific Oysters (Crassostrea gigas)

　　

富硒酵母是提高水产品营养价值与生理功能的重要有机硒源。然而，倍性是否影响太平洋牡蛎（Crassostrea gigas）在硒补充条件下的分子响应仍不清楚。在本研究中，研究人员对二倍体、三倍体和四倍体太平洋牡蛎投喂富硒饲料10天，采集肝胰腺组织进行抗氧化酶活性测定及转录组学分析。结果显示，三倍体牡蛎表现出更高的超氧化物歧化酶（superoxide dismutase, SOD）、谷胱甘肽过氧化物酶（glutathione peroxidase, GPx）和过氧化氢酶（catalase, CAT）活性，同时丙二醛（malondialdehyde, MDA）水平相对较低，表明其在硒补充条件下具有更强的抗氧化防御能力。转录组学分析显示，所有比较组间均存在大量差异表达基因（differentially expressed genes, DEGs），提示存在显著的倍性依赖性转录响应。功能富集分析表明，这些基因主要与代谢过程、氧化应激响应、膜组分、内分泌相关通路、脂质代谢和信号转导相关。多个氧化还原相关基因，包括谷胱甘肽S-转移酶（glutathione S-transferase, GST）家族基因、CHAC1、SOD家族基因和NRX1，呈现倍性依赖性表达模式。特别是，三倍体牡蛎似乎更依赖于酶促抗氧化防御和谷胱甘肽周转，而四倍体牡蛎则表现出更强的氧化还原信号整合和稳态调节倾向。这些发现为太平洋牡蛎倍性依赖性硒响应提供了新的见解，并为富硒牡蛎产品的开发利用提供了理论基础。

硒（Se）是维持机体正常生理功能所必需的微量元素，其生物学功能主要通过整合入硒蛋白（如谷胱甘肽过氧化物酶和硫氧还蛋白还原酶）并参与其功能而实现，在抗氧化防御、免疫调节、甲状腺激素代谢及维持细胞氧化还原稳态中发挥关键作用。硒的生理效应具有显著的剂量依赖性特征，缺乏和过量摄入均可导致代谢紊乱。因此，筛选安全、稳定且高生物利用度的天然硒源成为功能性水产品开发与营养强化研究的重要方向。

在水生生物中，硒利用效率与物种特异性生理及代谢特征密切相关，这为硒在海洋动物体内的积累与转化提供了生物学基础。牡蛎作为典型滤食性双壳贝类，具有强大的矿质元素积累与转化能力，是开发海洋功能性食品的理想候选物种。研究表明，牡蛎中的硒主要以硒代蛋氨酸（selenomethionine, C₅H₁₁NO₂Se）等有机形式存在，具有较高的生物利用度。在太平洋牡蛎中，硒代蛋氨酸可作为高效膳食硒源促进硒在组织中的积累与留存，并通过参与含硒蛋白合成提高抗氧化能力和应激抵抗力。这一特性赋予牡蛎较高的营养价值，使其成为开发富硒水产品的重要资源。然而，尽管目前已有研究探索了不同倍性牡蛎在生长性能、繁殖特性和环境适应性方面的差异，但关于营养代谢，特别是硒补充条件下硒的吸收、转化、利用及抗氧化调控机制方面的潜在差异研究仍较有限。

从分子角度来看，倍性变异可能影响基因剂量效应和转录调控网络，从而可能改变与硒相关的代谢和应激响应通路。已有研究表明，硒干预可显著改变牡蛎中代谢和抗氧化相关基因的表达，差异表达基因主要富集于抗氧化防御、物质代谢和免疫调节等通路。同时，针对不同倍性牡蛎的转录组学研究显示，倍性改变可重塑与生长、代谢和应激响应相关的多个分子通路。但截至目前，尚缺乏关于不同倍性牡蛎在统一硒营养干预下是否表现出不同转录响应模式的系统比较，尤其不清楚硒代谢相关和抗氧化调控相关通路是否呈现倍性依赖性调控特征。

基于上述背景，该研究以二倍体、三倍体和四倍体太平洋牡蛎为实验材料，在硒补充条件下，结合转录组测序技术系统比较了不同倍性个体间的转录响应，重点筛选与硒代谢、抗氧化防御及相关代谢过程有关的差异表达基因和功能通路，旨在阐明倍性差异对牡蛎硒响应策略及其分子调控基础的影响，为理解太平洋牡蛎倍性依赖性硒响应及其分子机制提供新见解。该成果发表于《Foods》。

研究样本来源于中国山东省烟台市崆峒岛养殖基地的二倍体、三倍体和四倍体太平洋牡蛎。主要关键技术方法包括：（1）流式细胞术进行倍性验证；（2）抗氧化酶活性测定（SOD、GPx、CAT及MDA）；（3）TRIzol法提取总RNA并构建cDNA文库；（4）Illumina NovaSeq 6000平台高通量转录组测序；（5）HISAT2进行参考基因组比对；（6）DESeq2软件进行差异表达基因分析；（7）GO和KEGG功能富集分析。

研究结果部分首先介绍了富硒酵母对太平洋牡蛎酶活性的影响。通过测定不同倍性牡蛎的抗氧化酶活性，发现三倍体牡蛎的SOD、GPx和CAT活性最高，四倍体居中，二倍体最低；而作为脂质过氧化产物的MDA含量则呈现相反趋势，二倍体最高，三倍体最低。这表明三倍体牡蛎在硒补充条件下具有更强的抗氧化防御能力。

在测序数据组装方面，转录组测序数据质量整体较高，二倍体、三倍体和四倍体样本的clean reads分别为43,182,512、44,689,791和44,264,606条，Q20和Q30值分别超过97.52%和92.98%，GC含量在43.40%-43.58%之间，测序数据质量良好，适用于后续分析。

差异表达基因分析显示，在相同富硒酵母干预下，不同倍性牡蛎表现出显著的基因表达差异。Se_3N vs. Se_2N、Se_4N vs. Se_2N和Se_4N vs. Se_3N三个比较组分别鉴定出1,564、1,446和1,407个差异表达基因，其中上调基因分别为806、854和727个，下调基因分别为758、592和680个。三个比较组共有36个核心差异表达基因。

功能注释与富集分析结果表明：GO富集分析显示差异表达基因主要涉及代谢过程、结合功能和膜结构等基本生物学活动，在生物过程方面显著关联生物调控、刺激响应和信号转导；在分子功能方面涉及信号转导和转录调控相关功能；在细胞组分方面主要集中于膜和膜部分、细胞器及细胞外区域等结构。KEGG通路富集分析显示，差异表达基因主要富集于六大功能类别，其中代谢相关通路富集基因数最多，其次是遗传信息处理和膜运输等；进一步分析显示，内分泌系统（如胰岛素信号通路）、脂质代谢（如脂肪酸降解）和信号转导（如MAPK信号通路）相关通路富集显著性较高（-log₁₀(p-value) > 5）。

讨论部分围绕不同倍性牡蛎在富硒酵母处理下的氧化还原调控特征展开深入分析。关于谷胱甘肽代谢相关基因的倍性依赖性调控机制，研究人员发现GST家族基因（GSTO1、GST-8、GSTT1和GSTN等）在不同倍性比较中显著下调，推测倍性差异可能影响谷胱甘肽代谢通路内的亚型特异性功能分工，而非简单诱导全局性抗氧化增强。特别是GSTO1在Se_3N vs. Se_2N和Se_4N vs. Se_2N中均下调，提示高倍性牡蛎在硒营养干预下可能不主要依赖GSTO1介导的经典抗氧化通路，而是通过其他GST亚型、SOD家族成员及相关氧化还原调控分子的协同调控维持细胞稳态，这更可能反映抗氧化通路的功能重组和补偿性调节。

CHAC1基因呈现显著的倍性依赖性差异表达模式：在Se_3N vs. Se_2N中上调，但在Se_4N vs. Se_3N中下调。该基因编码γ-谷氨酰环转移酶家族成员，可特异性降解谷胱甘肽（glutathione, GSH），在调控细胞氧化还原稳态中起关键作用。CHAC1在三倍体中的上调可能反映其在硒营养干预下更高的GSH周转需求，结合三倍体中更高的SOD、GPx和CAT活性，推测三倍体可能通过加速GSH更新与抗氧化酶协同作用，更有效地缓冲氧化压力。而四倍体中CHAC1的下调则表明其更倾向于抑制过度GSH分解以维持稳定的细胞内还原环境。

关于SOD基因家族在不同倍性抗氧化酶系统中的表达特征，研究发现多个胞质Cu/Zn-SOD相关基因（SODC、SOD1和SOD3）在三倍体和四倍体中较二倍体整体呈上调趋势，与三倍体较高的SOD酶活性一致，表明高倍性个体在硒处理下可能具有更强的超氧阴离子清除能力。特别是胞外Cu/Zn-SOD成员SOD3在Se_4N vs. Se_2N中的上调尤为值得关注，提示四倍体牡蛎的抗氧化调控可能不仅限于细胞内ROS清除，还可能延伸至细胞外微环境水平，具有更广泛的抗氧化调控范围。

硫氧还蛋白相关基因的表达变化揭示了不同倍性的调控层级差异。NRX1（nucleoredoxin 1）在Se_4N vs. Se_3N中显著上调，表明四倍体牡蛎在硒营养干预下可能更多依赖硫氧还蛋白相关调控机制。NRX1属于非经典硫氧还蛋白家族，除参与氧化还原反应外，还能调控Wnt/β-catenin和MAPK等多种氧化还原敏感信号通路，在维持细胞稳态和环境适应中发挥上游调控作用。与三倍体相比，四倍体中NRX1的更明显上调，暗示氧化还原信号和细胞稳态调控可能参与其对硒补充的响应。研究人员推断，随着倍性增加，牡蛎在硒处理下的氧化还原调控模式可能从"代谢抗氧化"模式向"基于信号的稳态调控"模式转变。

基于上述分析，研究人员提出了不同倍性太平洋牡蛎在富硒酵母补充下的氧化还原调控整合模型。外源富硒酵母处理可能影响活性氧（reactive oxygen species, ROS）生成与清除之间的平衡，从而触发涉及谷胱甘肽代谢通路、SOD酶系统和硫氧还蛋白相关调控网络的协同响应。在该模型中，GST家族成员（GSTO1-like、GSTT-like）、CHAC1、SOD1/SOD3和NRX1等关键差异表达基因被明确映射到相应的功能模块，突出其在不同倍性背景下调节氧化还原稳态中的潜在作用。结合生理指标与转录组证据，三倍体牡蛎以更高的SOD、GPx、CAT活性和较低的MDA水平为特征，可能通过增强的酶促抗氧化活性和谷胱甘肽周转具有更强的缓解氧化应激的能力；而四倍体牡蛎中NRX1的更明显上调则表明其更依赖氧化还原信号和稳态调节维持细胞平衡。

研究结论指出，在相同富硒酵母营养干预下，不同倍性太平洋牡蛎表现出不同的抗氧化生理响应和转录表达谱。三倍体牡蛎表现出更高的SOD、GPx和CAT活性以及较低的MDA水平，表明其在硒补充下具有更强的酶促抗氧化能力。转录组学分析显示，差异表达基因主要富集于抗氧化防御、物质代谢、免疫调节和信号转导等通路。总体而言，不同倍性牡蛎对富硒酵母采取不同的氧化还原调控策略：三倍体更倾向于依赖直接的酶促抗氧化防御，而四倍体则更强调氧化还原信号协调和细胞稳态维持。这些发现凸显了硒响应性抗氧化调控的倍性依赖性差异，为多倍体太平洋牡蛎的功能适应提供了新见解。该研究目前仅限于转录组水平分析和生理指标检测，未对关键基因或蛋白进行直接功能验证。未来研究应结合蛋白质组学和功能实验进一步验证候选基因（如GST家族成员、CHAC1和NRX1）的作用，并更详细地阐明其潜在调控机制。