2. 广州飞禧特生物科技有限公司, 广州 510640;
3. 北京奥特奇生物制品有限公司, 北京 101499
2. Guangzhou Fishtech Biotechnology Co., Ltd., Guangzhou 510640, China;
3. Beijing Otqi Biological Products Co., Ltd., Beijing 101499, China
鱼类属于变温动物，水温是影响鱼类生理、行为及进化过程中最重要的环境因素[1]。低温应激时，动物可能有不同的营养策略以提高自身的低温适应能力[2]。硒(Se)是动物机体必需的微量元素，具有促生长、抗氧化、提高免疫等功能。硒在体内主要以硒蛋白的形式发挥作用，硒蛋白具有类胰岛素的作用，影响内分泌激素[3]，对促进葡萄糖的摄入、调节糖酵解、糖异生以及脂肪酸的合成具有重要的作用[4]。我们前期研究发现，摄食酵母硒为硒源饲料的黄颡鱼(Pelteobrus fulvidraco)血液中葡萄糖、甘油三酯含量显著高于摄食亚硒酸钠为硒源饲料的黄颡鱼[5]，但是到目前为止这2种硒源发挥作用的分子机制尚不清楚。
转录组是指在特定的发育时期或生理状态下，细胞内所有转录本(transcript)的集合，它能够从整体水平揭示特定生物学过程发生的分子机理，为我们探究动物营养调控与抗环境胁迫机制提供了必要的分子手段[6]。因此，本试验采用转录组测序(RNA-seq)和实时荧光定量PCR(RT-qPCR)技术挖掘低温应激时摄食酵母硒和亚硒酸钠的黄颡鱼肝脏组织的关键信号通路和基因，为深入研究不同硒源的调控机理奠定基础，也为黄颡鱼抗低温功能饲料的研制提供理论依据。
1 材料与方法 1.1 试验饲料以酪蛋白为蛋白质源，玉米淀粉为糖源，鱼油、磷脂油为脂肪源配制基础饲料，其组成及营养水平见表 1。在基础饲料中分别添加亚硒酸钠和酵母硒，硒的实际添加量均为0.23 mg/kg[5]，配制2种试验饲料，添加亚硒酸钠的试验饲料硒含量实测值为0.22 mg/kg，添加酵母硒的试验饲料硒含量实测值为0.21 mg/kg。饲料原料粉碎后过60目筛，各原料采用逐级扩大法混合，然后于NH-10型捏合机中加入鱼油、磷脂油再进行二次混合，使用B20强力搅拌机搅拌均匀至成团，采用SLX-80型双螺杆挤压机制成粒径为2.0 mm的饲料，然后在G-500型造粒机中制粒，制粒完成后放入55~60 ℃烘箱中烘干，随后冷却，筛除粉末，用封口袋分装、标记，置于-20 ℃冰箱中待用。
|
表 1 基础饲料组成及营养水平(干物质基础) Table 1 Composition and nutrient levels of the basal diet (DM basis) |
养殖试验在广东省农业科学院动物科学研究所水产研究室的室内循环水养殖系统中进行。将5.7 g左右的黄颡鱼随机分成2组，每组4个重复，每个重复35尾，以重复为单位放养于水体容积300 L的养殖缸中。2组试验鱼分别投喂添加亚硒酸钠(T1组)和酵母硒(T2组)的试验饲料，每天饱食投喂2次(09:00和18:00)，持续投喂42 d。养殖期间，每天排污1次，每2周换水1次，水温26~30 ℃，pH 7.0~7.5，氨氮含量≤0.2 mg/L，亚硝酸盐含量≤0.02 mg/L，自然光照，24 h不间断曝气，每天记录投喂量和死亡情况。待养殖试验结束，将上述每缸中20尾试验鱼转移至可控温生态水槽系统的养殖缸(每缸水体容积150 L)中，降温前测得水温为26 ℃，降温速度设定为1 ℃/h，连续降温，当水温降到13 ℃时采样，低温应激期间禁食。
1.3 样品采集低温应激试验结束时每缸随机取3尾试验鱼，用40 mg/L MS-222麻醉后尾静脉取血，用于分析生化和抗氧化指标；然后在无菌环境下解剖取肝脏，装入无酶管，液氮速冻，于-80 ℃保存，用于转录组和基因表达分析。
1.4 RNA的提取及转录组文库的建立采用TRIzol(Invitrogen)法提取黄颡鱼肝脏中的总RNA，并使用DNase Ⅰ(TaKaRa)去除基因组DNA。分别采用2100 Bioanalyser(Agilent)、ND-2000(NanoDrop Technologies)检测RNA样品的质量，以保证使用合格的样品[OD260 nm/OD280 nm=1.8~2.2，OD260 nm/OD230 nm≥2.0，RNA完整值(RIN)≥6.5，28S/18S≥1.0，总量>2 μg]进行转录组测序。RNA文库的建立采用TruSeqTM RNA Sample Preparation Kit(Illumina)试剂盒。首先利用带有Oligo(dT)的磁珠从1 μg总RNA中富集有poly-A尾的mRNA。再加入Fragmentation Buffer，将mRNA随机断裂成200 bp左右的小片段。接着采用SuperScript Double-Stranded cDNA Synthesis Kit(Invitrogen)试剂盒，加入六碱基随机引物(Illumina)，以mRNA为模板反转录合成第1链cDNA，随后进行第2链合成，形成稳定的双链结构。双链的cDNA结构为黏性末端，加入End Repair Mix将其补成平末端，随后在3’末端加上1个A碱基，用于连接Y字形的接头，具体步骤参见试剂盒说明书。cDNA经过PCR富集后，使用DNA Clean Beads筛选200~300 bp的条带。经TBS380(Picogreen)定量后，文库使用Ilumina Noveseq 6000测序平台进行高通量测序，测序读长为PE 150。以上过程由上海美吉生物医药科技有限公司完成测序。
1.5 转录组数据分析将Ilumina Noveseq 6000测序平台得到的图像文件经CASAVA碱基识别分析转化为原始数据(raw data)，raw data经过去接头、低质量数据过滤筛选等流程后，得到高质量数据(clean data)。对mapped reads进行拼接，并与原有的基因组注释信息进行比较。本研究采用TPM(transcripts per million reads，即每百万reads中来自于某转录本的读段数)方法计算基因的表达量。再利用DEGseq软件包筛选差异表达基因，筛选阀值为差异倍数(fold change)>2和调整后P值(P-adjust) <0.05，且P-adjust越小基因表达量差异越显著。使用Goatools软件对差异表达基因进行GO功能注释和KEGG通路富集分析，且P < 0.05时表示差异表达基因显著富集。
1.6 RT-qPCR根据差异表达基因对低温应激的响应程度，本研究采用RT-qPCR试验对关键通路中的差异表达基因进行验证，以此评价转录组测序技术鉴定差异表达基因的可靠性。使用Primer 3.0对各差异表达基因进行引物设计(表 2)，并在上海美吉生物医药科技有限公司完成引物的合成。使用反转录试剂盒(BestarTM qPCR RT Kit)将各样品的RNA反转录为cDNA，并以反转录产物为模板。以18S RNA为内参，并使用ChamQ SYBR Color qPCR Master Mix (2×)试剂，采用ABI7300荧光定量PCR仪进行定量检测，每个样品3个平行，以减小误差。采用2-△△Ct的方法计算出各候选基因的相对表达量。
|
|
表 2 RT-qPCR引物 Table 2 Primers for RT-qPCR |
血清总蛋白、白蛋白、球蛋白、胆固醇、葡萄糖、甘油三酯含量采用日立全自动生化分析仪进行测定。
血清抗氧化指标均采用南京建成生物工程研究所生产的试剂盒测定，按照使用说明进行操作。
1.8 生长性能指标计算
|
生长性能、血清生化和抗氧化指标数据采用SPSS 17.0软件进行t检验，结果以平均值±标准差(mean±SD)表示，以P < 0.05为显著差异水平。
2 结果与分析 2.1 生长性能由表 3可知，与摄食以亚硒酸钠为硒源的饲料相比，摄食以酵母硒为硒源的饲料的黄颡鱼的生长性能无显著变化(P>0.05)，但增重率、特定生长率和存活率有升高的趋势。
|
|
表 3 摄食不同硒源饲料黄颡鱼的生长性能 Table 3 Growth performance of Pelteobagrus fulvidraco fed different selenium source diets (n=4) |
由表 3可知，与亚硒酸钠相比，低温应激后酵母硒显著提高了黄颡鱼血清甘油三酯含量以及抑制与产生超氧阴离子自由基能力(P < 0.05)，显著降低了血清丙二醛含量(P < 0.05)。
|
|
表 4 低温应激后摄食不同硒源饲料黄颡鱼的血清生化和抗氧化指标 Table 4 Serum biochemical and antioxidant indexes of Pelteobagrus fulvidraco fed different selenium source diets after low temperature stress (n=4) |
对T1组和T2组黄颡鱼肝脏组织进行转录组测序，共获得42.00 Gb clean data，各样品的clean data均达到6.68 Gb。测序样品的数据质量及比对结果见表 5，各样品的Q30占比均在93%以上，GC含量均在46%以上。将各样品的高质量reads(clean reads)与黄颡鱼基因组进行序列比对，比对率均在66%以上，表明此次测序数据质量较高，可用于后续分析。
|
|
表 5 测序数据质量和序列比对 Table 5 Sequencing data quality and sequence alignment |
基于表达量指标TPM，以fold change的绝对值>2，P-adjust < 0.05为阈值，对同一基因在不同组的表达量进行统计分析。与T1组相比，T2组有76个基因表达量上调，71个基因表达量下调，差异表达基因如图 1所示。
|
图 1 差异表达基因火山图(T1组vs. T2组) Fig. 1 Volcano plot of differentially expressed genes (T1 group vs.T2 group) |
Unigene在GO分类中被注释到生物过程、细胞组分和分子功能等3大类别中。与T1组相比，T2组肝脏基因表达谱主要与代谢过程、催化活性、结合、单组织过程、细胞进程等生物功能有关(图 2)，显著富集到蛋白质水解作用、水解酶活性、肽酶活性、肽链内切酶活性、羧肽酶活性、蛋白质代谢过程(图 3)。
|
X轴为注释的基因/转录本数量，Y轴为GO term，图中显示的是丰度前20个GO term。X-axis was the annotated number of genes/transcripts, and Y-axis was GO term. The figure showed the top 20 GO terms of abundance. Molecular function: 分子功能; Cellular component: 细胞组分; Biological process: 生物过程; Negative Regulation of biological process: 生物过程的负调控；Immune system process: 免疫系统过程; Developmental process: 发展过程; Locomotion: 移动; Response to stimulus: 刺激反应; Multicellular organismal process: 多细胞生物过程; Regulation of biological process: 生物过程调节; Biological regulation: 生物调节; Localization: 定位; Cellular process：细胞过程; Single-organism process: 单组织过程; Metabolic process: 代谢过程; Cell: 细胞; Cell part: 细胞部分; Membrane part: 膜部分; Extracellular region: 细胞外区域; Membrane: 膜; Molecular function: 分子功能调控; Binding: 结合; Catalytic activity: 催化活性。 图 2 差异表达基因的GO功能注释(T1组vs. T2组) Fig. 2 GO functional annotation of differentially expressed genes (T1 group vs. T2 group) |
|
X轴为富集的基因/转录本数量，Y轴为GO term，图中显示的是P-adjust最小的前20个GO term。X-axis was the enriched number of genes/transcripts, and Y-axis was GO term. The figure showed the top 20 GO terms with the smallest P-adjust. Proteolysis: 蛋白质水解过程; Serine-type endopeptidase activity: 丝氨酸肽链内切酶活性; Catalytic activity, acting on a protein: 作用于蛋白质的催化活性; Peptidase activity, acting on L-amino acid peptides: 作用于L-氨基酸肽的肽酶活性; Hydrolase activity: 水解酶活性; Endopeptidase activity: 肽链内切酶活性; Serine hydrolase activity: 丝氨酸水解酶活性; Serine-type peptidase activity: 丝氨酸肽酶活性; Hydrolase activity, acting on acid phosphorus-nitrogen bonds: 作用于酸性磷氮键的水解酶活性; Peptidase activity: 肽酶活性; Metallopeptidase activity: 金属肽酶活性; Metallocarboxypeptidase activity: 金属羧肽酶活性; Metalloexopeptidase activity: 金属外肽酶活性; Carboxypeptidase activity: 羧肽酶活性; Protein metabolic process: 蛋白质代谢过程。 图 3 差异表达基因的GO富集分析(T1组vs. T2组) Fig. 3 GO enrichment analysis of differentially expressed genes (T1 group vs. T2 group) |
Unigene在KEGG分类中被注释到代谢、环境信息处理、细胞过程、生物体系统、人类疾病等五大类别中。与T1相比，T2组差异表达基因主要注释到消化系统、信号分子与相互作用、信号传导、脂代谢、内分泌系统、免疫系统、糖代谢、氨基酸代谢、运输和分解代谢通路(图 4)，显著富集到胰腺分泌、蛋白质消化吸收、脂肪消化吸收、花生四烯酸代谢及醚脂类代谢通路(图 5)。
|
X轴为KEGG通路，Y轴为注释的基因/转录本数量。图中显示的是注释到的KEGG通路。X-axis was KEGG pathways, and Y-axis was the annotated number of genes/transcripts. The figure showed the annotated KEGG pathways. Metabolism: 代谢; Environmental information processing：环境信息处理; Cellular process: 细胞过程; Organismal system: 生物体系统; Human diseases: 人类疾病; Amino acid metabolism: 氨基酸代谢; Carbohydrate metabolism: 碳水化合物代谢; Lipid metabolism: 脂质代谢; Metabolism of cofactors and vitamins: 辅因子和维生素代谢; Nucleotide metabolism: 核苷酸代谢; Signal transduction: 信号传导; Signaling molecules and interaction: 信号分子及相互作用; Cell growth and death: 细胞生长与死亡; Cellular community-eukaryotes: 细胞群落真核生物; Transport and catabolism: 运输和分解代谢; Aging：老化; Circulatory system: 循环系统; Development: 发展; Digestive system: 消化系统; Endocrine system: 内分泌系统; Environmental adaptation: 环境适应; Immune system: 免疫系统; Nervous system: 神经系统; Cancers: Overview: 癌症：概述; Cancers: Specific types: 癌症: 特定类型; Cardiovascular diseases: 心血管疾病; Drug resistance: Antineoplastic: 耐药性: 抗肿瘤药物; Endocrine and metabolic diseases: 内分泌和代谢疾病; Immune diseases: 免疫疾病; Infectious diseases: Bacterial: 传染病: 细菌; Infectious diseases: Parasitic: 感染性疾病: 寄生; Infectious diseases: Viral: 感染性疾病: 病毒; Neurodegenerative diseases: 神经推行性疾病; Substance dependence: 药物依赖。 图 4 差异表达基因的KEGG功能注释(T1组vs. T2组) Fig. 4 KEGG functional annotation of differentially expressed genes (T1 group vs. T2 group) |
|
X轴为富集的基因/转录本数量，Y轴为KEGG通路。图中显示的是所有富集的KEGG通路。X-axis was the annotated number of genes/transcripts, and Y-axis was KEGG pathways. The figure showed all the KEGG pathways. Pancreatic secretion: 胰腺分泌; Protein digestion and absorption: 蛋白质消化吸收; Influenza A: 甲型流感; Fat digestion and absorption: 脂肪消化与吸收; Neuroactive ligand-receptor interaction: 刺激中枢神经的交互作用; Arachidonic acid metabolism: 花生四烯酸代谢; Ether lipid metabolism: 醚脂类代谢。 图 5 差异表达基因的KEGG富集分析(T1组vs. T2组) Fig. 5 KEGG enrichment analysis of differentially expressed genes (T1 group vs. T2 group) |
如表 6所示，胰腺分泌、蛋白质消化与吸收、脂肪消化与吸收、花生四烯酸代谢、醚脂代谢等通路的类胰蛋白酶-1、羧肽酶A1、羧肽酶B、磷脂酶A2等关键基因的表达都显著上调。利用RT-qPCR对蛋白质消化与吸收通路差异表达基因硒蛋白M、磷脂酶A2、羧肽酶A1、类胰蛋白酶-1、胰凝乳蛋白酶样弹性蛋白酶家族成员2A的结果进行验证，结果中表达显著上调的基因其验证结果也为表达显著上调，表明基于转录组测序数据的差异表达基因分析结果可信(表 7)。
|
|
表 6 低温应激后相关差异表达基因的表达 Table 6 Expression of related differential expressed genes after low temperature stress |
|
|
表 7 差异表达基因的qRT-PCR分析 Table 7 RT-qPCR analysis of differentially expressed genes |
生长性能结果显示，酵母硒组黄颡鱼的生长性能优于亚硒酸钠组，这一结果与我们前期开展的8周养殖试验获得的结果[5]一致，同样，低温应激后血清谷胱甘肽过氧化物酶、丙二醛以及抑制与产生超氧阴离子自由基所反映的抗氧化能力提高也与前期研究结果[5]一致。
温度作为对鱼体最重要的影响因素之一，对鱼体自身的新陈代谢起着至关重要的作用。随着水温的持续降低以及应激时间的延长，鱼体内环境出现紊乱。研究显示，部分鱼类承受低温应激能力较强，最终可通过自身调节系统使机体适应外界环境的变化，以维持或重建鱼体内环境的稳态，从而减少或避免应激造成的损伤；但一部分鱼类自身抵抗力差，持续的低温应激对自身机体耗能过多，其自身的调节系统不足以使机体内环境再次达到稳定平衡的状态，最终鱼类自身调节系统衰竭，造成死亡[7]。鱼类的内分泌机能以及蛋白质、脂肪和糖类等营养物质的代谢对低温应激敏感[2, 8-9]，低温应激从多方面影响鱼类的生理活动，显然单独研究某个或某些基因不能揭示信号调控网络的复杂性，转录组测序技术为我们从整体水平揭示特定生物学过程发生的分子机理提供了研究手段。
本研究通过转录组分析发现酵母硒组与亚硒酸钠组之间共产生了147个差异表达基因，GO功能注释表明，差异表达基因的功能主要在蛋白质水解作用、催化活性(作用于蛋白质)、水解酶活性、肽酶活性、水解酶活性(作用于酸性磷氮键)、肽链内切酶活性、蛋白质代谢过程、脂质分解过程富集，KEGG功能注释显示差异表达基因主要注释到消化系统、信号分子与相互作用、信号传导、脂代谢、内分泌系统、免疫系统、碳水化合物代谢、氨基酸代谢、运输和分解代谢通路，并显著富集到胰腺分泌、蛋白质消化吸收、脂肪消化吸收、花生四烯酸代谢、醚脂类代谢。上述结果表明，与亚硒酸钠比较，酵母硒激活了内分泌，蛋白质、氨基酸、脂肪和糖类代谢等多个生物学过程，以调节黄颡鱼适应低温应激。鱼类为了适应环境胁迫，机体表现出神经内分泌活动的改变，低温应激改变了水产动物的内分泌系统[10-11]，当虹鳟(Oncorhynchus mykiss)处于急性低温胁迫时，正常的甲状腺功能降低[12]。硒具有调节机体内分泌激素的作用，硒蛋白P和谷胱甘肽过氧化物酶与胰岛β细胞的功能有关，具有降低血糖和调控胰岛素介导的代谢过程等拟胰岛素作用[13]，而碘甲腺原氨酸脱碘酶能够介导甲状腺激素的形成和降解，具有调节甲状腺激素代谢平衡的作用[14]。本试验中，低温应激后酵母硒组的黄颡鱼肝脏中23个差异表达基因富集在胰液分泌通路且显著上调，这一结果也证明硒参与了黄颡鱼低温应激时的内分泌调控，与前述研究结果一致。有机硒的吸收效率和生物利用率高于无机硒[5, 15-16]，这可能是导致酵母硒组黄颡鱼内分泌通路显著上调的原因，这也恰恰说明当黄颡鱼处于低温应激状态时，饲料中添加酵母硒能够增强黄颡鱼低温应激时的内分泌机能，提高黄颡鱼的抗寒能力。
胰腺分泌的有机物是多种消化酶，可作用于糖类、脂肪和蛋白质的转运及消化过程[17-18]，其中胰腺脂肪酶包括甘油三酯脂酶、磷脂酶A2等，是主要的脂肪消化酶[19-20]。李敏芬[21]研究发现，急性低温胁迫时乌鳢(Channa argus)和乌斑鳢(Channa maculata)的脂肪酶活性显著上升，推测乌鳢及乌斑鳢对低温的适应性可能较好。此外，在暗纹东方鲀(Takifugu fasciatus)[22]、斑马鱼(Danio rerio)[23]、牙鲆(Paralichthys olivaceus)[24]、金头鲷(Sparus aurata[25]、鲤鱼(Cyprinus carpio L.)[26]等鱼类中发现低温下脂类代谢增强，脂类代谢通路在鱼类应对低温时发挥了关键作用。本试验中，酵母硒提高了胰液分泌通路上胰腺三酰甘油、磷脂酶A2、羧基酯脂肪酶的表达，同时脂肪消化与吸收通路显著富集，血清甘油三酯含量显著升高。羧基酯脂肪酶主要参与消化甘油三酯，为生命提供能量的物质，说明酵母硒提高了黄颡鱼的脂肪分解供能能力，为机体积极应对低温应激提供热量。研究显示，酵母硒本身具有显著提高鳡鱼(Elopichthys bambusa)[27]、草鱼(Ctenopharyngodon idellus)[28]、刺参(Apostichopus japonicas)[29]脂肪酶活性的能力，本研究结果与此一致。然而，在我们前期进行的8周养殖试验的低温应激结果中，酵母硒组和亚硒酸钠组黄颡鱼血清甘油三酯含量并未出现显著差异[5]。鱼体的应激调控系统是一个多层次、多功能的结构体系[30]，动物的分子响应与代谢产物产生的时空特性还有待于进一步研究。
蛋白质作为体内重要的营养物质，首先在胃和肠道中被胃蛋白酶、胰蛋白酶以及肽酶等降解生成肽类物质和游离氨基酸。动物在面临低温应激时自动分解的氨基酸，一方面用于蛋白质的合成与周转，另一方面还可能承担着抗应激的功能[2]。研究发现，凡纳滨对虾(Litopenaeus vannamei)在低温应激时肝胰腺启动了蛋白质的分解机制以弥补体内氨基酸的缺乏[31]，鳢在短时间低温胁迫时胰蛋白酶活性显著升高，通过消耗机体蛋白质来应对低温胁迫[21]。本试验结果显示，有18个差异表达基因显著富集在蛋白质消化与吸收通路，该通路中类胰蛋白酶-1、胰凝乳蛋白酶样弹性蛋白酶家族成员2A、羧肽酶A1和羧肽酶B的表达量在酵母硒组显著上调，表明酵母硒促使黄颡鱼在低温应激时蛋白质分解作用增强。本试验中酵母硒组黄颡鱼血清总蛋白含量高于亚硒酸钠组，但差异不显著。目前，硒影响动物蛋白质代谢的报道主要集中在正常养殖环境下硒促进蛋白质合成方面，例如，酵母硒通过提高虹鳟肌肉的生长以及肌肉蛋白质的含量实现生长[32-33]；硒激活胰岛素信号级联和调节蛋白质代谢，有利于骨骼肌发育[34]；母体补硒提高肉鸡雉鸡胸肌含量以及雷帕霉素的基因表达等[35]。但有关硒影响蛋白质水解生成氨基酸的研究鲜有报道。
糖类是低温应激下鱼类的重要能源物质，机体为了抵御寒冷，可能会加快糖的分解利用。研究认为，在应激初期，糖原转化成葡萄糖的速率加快，血糖水平上升；在应激后期，血糖被大量分解生成ATP，以给机体提供能量，此时血糖水平会降低[36]。鲤鱼、在急性低温胁迫时糖酵解/糖质新生过程被显著富集，关键酶葡萄糖-6-磷酸脱氢酶(G6PC)、果糖-1, 6-二磷酸酶(FBP)和磷酸烯醇式丙酮酸羧化酶(PCK)的基因表达量在鲤低温适应过程中均显著升高[37]。然而，本试验条件下未筛选到糖代谢通路。在我们的前期研究中发现当温度下降到13℃时酵母硒组黄颡鱼血清葡萄糖含量高于亚硒酸钠组且达到显著水平[5]，在此次试验中酵母硒组血清葡萄糖含量也高于亚硒酸钠组，但未达到显著水平。因此，本试验中导致糖代谢通路未被注释和富集的原因需要进一步研究。
硒蛋白M是一种内质网硒蛋白，具有维持机体氧化还原水平、保护神经细胞和调节钙稳态平衡等重要作用。研究发现，凡纳滨对虾通过RNA干扰RNAi沉默硒蛋白M会降低鳃部过氧化物酶的活性，硒蛋白M是调节鳃氧化应激的关键酶，可以作为肝胰腺的第二信使调节因子在氧化还原系统中发挥作用[38]。当外源性能源的获取受到抑制时，机体需要动员更多的内源性能源物质或潜在产能底物来提供能量应对低温环境，更多参与转录正调控、信号转导、能量代谢、内分泌等生理活动的基因表达激活、上调都有助于对低温环境的适应[38]。本研究中，酵母硒组黄颡鱼肝脏硒蛋白M的表达显著上调，推测酵母硒可能通过提高硒蛋白M的表达进而上调胰腺分泌、蛋白质消化与吸收、脂肪消化与吸收等通路关键基因的表达，发挥酵母硒优于亚硒酸钠抗低温应激的效果。
4 结论与亚硒酸钠相比，酵母硒可能通过提高黄颡鱼硒蛋白M的表达，上调胰腺分泌、蛋白质消化与吸收、脂肪消化与吸收、花生四烯酸代谢、醚脂代谢等通路中胰蛋白酶、羧肽酶、磷脂酶、羧基酯脂肪酶等关键基因的表达，进而提高黄颡鱼的抗低温应激能力。
| [1] |
陈超, 施兆鸿, 薛宝贵, 等. 低温胁迫对七带石斑鱼幼鱼血清生化指标的影响[J]. 水产学报, 2012, 36(8): 1249-1255. CHEN C, SHI Z H, XUE B G, et al. Influence of low-temperature stress on serum biochemical parameters in juvenile Epinephelus septemfasciatus[J]. Journal of Fisheries of China, 2012, 36(8): 1249-1255 (in Chinese). |
| [2] |
周萌, 王安利, 苗玉涛, 等. 虾蟹低温应激研究进展[J]. 饲料工业, 2011, 32(22): 45-48. ZHOU M, WANG A L, MIAO Y T, et al. Research of shrimp and crab response at low temperature[J]. Feed Industry, 2011, 32(22): 45-48 (in Chinese). DOI:10.3969/j.issn.1001-991X.2011.22.013 |
| [3] |
STEINBRENNER H, SPECKMANN B, PINTO A, et al. High selenium intake and increased diabetes risk: experimental evidence for interplay between selenium and carbohydrate metabolism[J]. Journal of Clinical Biochemistry and Nutrition, 2011, 48(1): 40-45. |
| [4] |
ZENG H W, COMBS GF JR. Selenium as an anticancer nutrient: roles in cell proliferation and tumor cell invasion[J]. The Journal of Nutritional Biochemistry, 2008, 19(1): 1-7. DOI:10.1016/j.jnutbio.2007.02.005 |
| [5] |
胡俊茹, 王国霞, 孙育平, 等. 亚硒酸钠和酵母硒对黄颡鱼幼鱼生长性能、抗氧化能力及抗低温应激的影响[J]. 水产学报, 2019, 43(11): 2394-2404. HU J R, WANG G X, SUN Y P, et al. Effects of dietary sodium selenite and selenoyeast on growth performance, antioxidant responses and low temperature stress resistance of juvenile yellow catfish (Pelteobagrus fulvidraco)[J]. Journal of Fisheries of China, 2019, 43(11): 2394-2404 (in Chinese). |
| [6] |
罗辉, 叶华, 肖世俊, 等. 转录组学技术在水产动物研究中的运用[J]. 水产学报, 2015, 39(4): 598-607. LUO H, YE H, XIAO S J, et al. Application of transcriptomics technology to aquatic animal research[J]. Journal of Fisheries of China, 2015, 39(4): 598-607 (in Chinese). |
| [7] |
贾明亮. 低温胁迫对奥尼罗非鱼的生长、肌肉组成和血液生理生化指标的影响[D]. 硕士学位论文. 湛江: 广东海洋大学, 2010: 1-2. JIA M L. The effect of low temperature stress on growth, muscle composition and hematological and biochemical indices of tilapia (Oreochromis niloticus×O. aures)[D]. Master's Thesis. Zhanjiang: Guangdong Ocean University, 2010: 1-2. (in Chinese) |
| [8] |
HSIEH S L, CHEN Y N, KUO C M. Physiological responses, desaturase activity, and fatty acid composition in milkfish (Chanos chanos) under cold acclimation[J]. Aquaculture, 2003, 220(1/2/3/4): 903-918. |
| [9] |
SILBERGELD E K. Blood glucose: a sensitive indicator of environmental stress in fish[J]. Bulletin of Environmental Contamination and Toxicology, 1974, 11(1): 20-25. DOI:10.1007/BF01685023 |
| [10] |
PICKERING A D. Rainbow trout husbandry: management of the stress response[J]. Aquaculture, 1992, 100(1/2/3): 125-139. |
| [11] |
PAVLIDIS M A, MYLONAS C C. Sparidae: biology and aquaculture of gilthead seabream and other species[M]. [S. l. ]: Blackwell Publishing Ltd., 2011.
|
| [12] |
管标. 急性温度胁迫对虹鳟血液成分和关键代谢的影响[D]. 硕士学位论文. 青岛: 中国海洋大学, 2014: 38-39. GUAN B. The effect of acute temperature stress on blood components and key metabolism mechanisms in rainbow trout (Oncorhynchus mykiss)[D]. Master's Thesis. Qingdao: Ocean University of China, 2014: 38-39. (in Chinese) |
| [13] |
STEINBRENNER H. Interference of selenium and selenoproteins with the insulin-regulated carbohydrate and lipid metabolism[J]. Free Radical Biology & Medicine, 2013, 65: 1538-1547. |
| [14] |
倪银星. 硒蛋白、硒与内分泌激素的关系研究进展[J]. 国外医学(卫生学分册), 2002, 29(1): 38-41, 57. NI Y X. Research progress on the relationship between selenium protein, selenium and endocrine hormone[J]. Foreign Medical Sciences(Section Hygiene), 2002, 29(1): 38-41, 57 (in Chinese). |
| [15] |
WANG C, LOVELL R T. Organic selenium sources, selenomethionine and selenoyeast, have higher bioavailability than an inorganic selenium source, sodium selenite, in diets for channel catfish (Ictalurus punctatus)[J]. Aquaculture, 1997, 152(1/2/3/4): 223-234. |
| [16] |
LIN Y H. Effects of dietary organic and inorganic selenium on the growth, selenium concentration and meat quality of juvenile grouper Epinephelus malabaricus[J]. Aquaculture, 2014, 430: 114-119. DOI:10.1016/j.aquaculture.2014.03.048 |
| [17] |
尤宏争. 鱼类消化系统及消化酶的研究进展[J]. 河北渔业, 2012(3): 53-61. YOU H Z. Research progress on digestive system and digestive enzymes of fish[J]. Hebei Fisheries, 2012(3): 53-61 (in Chinese). |
| [18] |
杨瑞斌. 黄颡鱼仔稚鱼消化系统发育及摄食特性的研究[D]. 博士学位论文. 武汉: 华中农业大学, 2010: 43-44. YANG R B. Study on ontogenetic development of the digestive system and feeding in yellow catfish Pelteobagrus fulvidraco larvae[D]. Ph. D. Thesis. Wuhan: Huazhong Agricultural University, 2010: 43-44. (in Chinese) |
| [19] |
覃川杰. 瓦氏黄颖鱼(Pelteobagrus vachelli)脂肪代谢相关基因cDNA的克隆及表达分析[D]. 博士学位论文. 上海: 华东师范大学, 2010: 8-9. QIN C J. The lipid metabolism-related gene from darkbarbel catfish, Pelteobagrus vachelli juvenile: cDNA cloning characteration and mRNA expression pattern[D]. Ph. D. Thesis. Shanghai: East China Normal University, 2010: 8-9. (in Chinese) |
| [20] |
周姗姗, 李晓南. 胰腺脂肪酶相关研究进展[J]. 生理科学进展, 2015, 46(1): 6-10. ZHOU S S, LI X N. The progress of research about pancreatic lipase[J]. Progress in Physiological Sciences, 2015, 46(1): 6-10 (in Chinese). |
| [21] |
李敏芬. 低温胁迫对乌鳢、斑鳢及其杂交种部分生理指标的影响[D]. 硕士学位论文. 上海: 上海海洋大学, 2016: 23-24. LI M F. Effects of low-temperature stress on some physiological indexes in Channa argus, Channa maculata and their hybrids[D]. Master's Thesis. Shanghai: Shanghai Ocean University, 2016: 23-24. (in Chinese) |
| [22] |
文鑫. 暗纹东方鲀(Takifugu fasciatus)应对低温胁迫的生理响应和分子机制研究[D]. 博士学位论文. 南京: 南京师范大学, 2019. WEN X. The molecular regulatory mechanism and physiological response of Takifugu fasciatus to low-temperature stress[D]. Ph. D. Thesis. Nanjing: Nanjing Normal University, 2019. (in Chinese) |
| [23] |
鲁栋梁. 饥饿提高斑马鱼(Danio rerio)冷应激耐受的机制探究[D]. 博士学位论文. 上海: 华东师范大学, 2019. LU D L. Mechanisms of fasting enhance cold tolerance in zebrafish (Danio rerio)[D]. Ph. D. Thesis. Shanghai: East China Normal University, 2019. (in Chinese) |
| [24] |
HU J W, YOU F, WANG Q, et al. Transcriptional responses of olive flounder (Paralichthys olivaceus) to low temperature[J]. PLoS One, 2014, 9(10): e108582. DOI:10.1371/journal.pone.0108582 |
| [25] |
MATEUS A P, COSTA R, GISBERT E, et al. Thermal imprinting modifies bone homeostasis in cold-challenged sea bream (Sparus aurata)[J]. The Journal of Experimental Biology, 2017, 220(Pt 19): 3442-3454. |
| [26] |
DIETRICH M A, HLIWA P, ADAMEK M, et al. Acclimation to cold and warm temperatures is associated with differential expression of male carp blood proteins involved in acute phase and stress responses, and lipid metabolism[J]. Fish & Shellfish Immunology, 2018, 76: 305-315. |
| [27] |
许友卿, 李太元, 丁兆坤, 等. 添加酵母硒对鳡鱼消化酶活性与饲料转化率的影响[J]. 水产科学, 2013, 32(7): 391-395. XU Y Q, LI T Y, DING Z K, et al. Effects of dietary selenium on digestion and feed conversion in yellowcheek carp Elopichthys bambusa[J]. Fisheries Science, 2013, 32(7): 391-395 (in Chinese). DOI:10.3969/j.issn.1003-1111.2013.07.004 |
| [28] |
苏传福, 罗莉, 文华, 等. 硒对草鱼生长、营养组成和消化酶活性的影响[J]. 上海水产大学学报, 2007, 16(2): 124-129. SU C F, LUO L, WEN H, et al. Effects of dietary selenium on growth performance, quality and digestive enzyme activities of grass carp[J]. Journal of Shanghai Fisheries University, 2007, 16(2): 124-129 (in Chinese). |
| [29] |
刘金凤. 富硒酵母的制备及在刺参养殖中的应用[D]. 硕士学位论文. 大连: 大连工业大学, 2016: 25-26. LIU J F. Preparation of selenium enriched yeast and applicationin in aquaculture of sea cucumber Apostichopus japonicas[D]. Master's Thesis. Dalian: Dalian Polytechnic University, 2016: 25-26. (in Chinese) |
| [30] |
康玉军. 虹鳟肝脏响应高温胁迫的蛋白质组学与代谢组学研究[D]. 博士学位论文. 兰州: 甘肃农业大学, 2020: 16. KANG Y J. Proteomics and metabolomics analysis of rainbow trout(Oncorhynchus mykiss)liver responses to heat stress[D]. Ph. D. Thesis. Lanzhou: Gansu Agricultural University, 2020: 16. (in Chinese) |
| [31] |
ZHOU M, WANG A L, XIAN J A, et al. Variation of free amino acid and carbohydrate concentrations in white shrimp, Litopenaeus vannamei: effects of continuous cold stress[J]. Aquaculture, 2011, 317(1/2/3/4): 182-186. |
| [32] |
WANG L, ZHANG D F, LI S, et al. Dietary selenium promotes somatic growth of rainbow trout (Oncorhynchus mykiss) by accelerating the hypertrophic growth of white muscle[J]. Biological Trace Element Research, 2021, 199(5): 2000-2011. DOI:10.1007/s12011-020-02282-w |
| [33] |
WANG L, WU L, LIU Q, et al. Improvement of flesh quality in rainbow trout (Oncorhynchus mykiss) fed supranutritional dietary selenium yeast is associated with the inhibited muscle protein degradation[J]. Aquaculture Nutrition, 2018, 24(4): 1351-1360. |
| [34] |
ZENG M S, LI X, LIU Y, et al. A high-selenium diet induces insulin resistance in gestating rats and their offspring[J]. Free Radical Biology & Medicine, 2012, 52(8): 1335-1342. |
| [35] |
GAO J, NIE W, WANG F L, et al. Maternal selenium supplementation enhanced skeletal muscle development through increasing protein synthesis and SelW mRNA levels of their offspring[J]. Biological Trace Element Research, 2018, 186(1): 238-248. |
| [36] |
常玉梅, 曹鼎臣, 孙效文, 等. 低温胁迫对鲤血液学和血清生化指标的影响[J]. 水产学杂志, 2006, 19(2): 71-75. CHANG Y M, CAO D C, SUN X W, et al. Changes of serum biochemical indices of common carp affected by cold temperatures[J]. Chinese Journal of Fisheries, 2006, 19(2): 71-75 (in Chinese). |
| [37] |
刘思嘉, 田菲, 张存芳, 等. 鲤在低温胁迫下肝胰腺转录组测序分析[J]. 生物技术通报, 2018, 34(11): 168-178. LIU S J, TIAN F, ZHANG C F, et al. Transcriptome sequencing and analysis of hepatopancreas from carps under cold stress[J]. Biotechnology Bulletin, 2018, 34(11): 168-178 (in Chinese). |
| [38] |
GARCÍA-TRIANA A, PEREGRINO-URIARTE A B, YEPIZ-PLASCENCIA G. Selenoprotein M gene expression, peroxidases activity and hydrogen peroxide concentration are differentially regulated in gill and hepatopancreas of the white shrimp Litopenaeus vannamei during hypoxia and reoxygenation[J]. Comparative Biochemistry and Physiology. Part A, Molecular & Integrative Physiology, 2016, 199: 14-20. |