畜牧兽医学报  2019, Vol. 50 Issue (9): 1746-1758. DOI: 10.11843/j.issn.0366-6964.2019.09.003    PDF    
去势对淮南公猪背最长肌转录组的影响
王璟1, 滑留帅1, 陈俊峰1, 张家庆1, 任巧玲1, 白红杰1, 郭红霞1, 徐照学1, 邢宝松1, 白献晓1, 曹海2     
1. 河南省农业科学院 畜牧兽医研究所 河南省畜禽繁育与营养调控重点实验室, 郑州 450002;
2. 河南兴锐农牧科技有限公司, 信阳 465550
摘要:旨在分析去势对猪背最长肌基因表达的影响,探讨去势对肉质性状的分子调控机制。本研究将6头健康公猪分为两组,去势组和非去势组(对照组),在体重达到130 kg时(大约300~315日龄)采集背最长肌样品,利用Illumina HiSeq 2000高通量RNA-seq测序技术对两组猪背最长肌进行转录组测序,使用DESeq软件筛选差异表达基因,并在GO和KEGG数据库中对这些差异表达基因进行功能注释和富集分析。结果显示,去势和非去势淮南公猪分别获得83 613 018和83 746 508条可用读段,与猪参考基因组(Sscrofa10.2)的比对率分别为73.78%和74.09%。与非去势组相比,去势后共有差异表达基因935个,其中上调基因503个,下调基因432个(P < 0.05且|log2(Fold Change)|>0.8)。其中与肌肉发育相关上/下调最显著的基因分别是胶原蛋白XI型α1(COL11A1)和血管紧张素Ⅱ1型受体(AGTR1),而与脂肪代谢相关上/下调最显著的基因分别是硬脂酰辅酶A去饱和酶1(SCD-1)和脂肪特异性磷脂酶A2(PLA2G16)。GO和KEGG通路分析发现,两组差异表达基因显著富集于胰岛素、脂解、脂肪酸延长等脂质代谢相关通路,以及肾上腺素能、心肌收缩、磷酸腺苷活化的蛋白激酶等肌肉发育相关通路。本研究表明,去势后公猪脂肪沉积能力增强,肌肉发育改变,肉质性状随之改变。SCD-1、激素敏感脂酶(HSL)、葡萄糖转运体4(GLUT4)、丝氨酸/苏氨酸激酶(AKT)、促分裂原活化蛋白激酶(MEK)等12个差异表达基因能同时或部分参与激素分泌、脂肪沉积和肌肉发育的调控,这些基因与去势后肉质性状的改变显著相关,值得进行进一步的功能验证。
关键词淮南猪    肉质    转录组测序    肌肉发育    脂肪沉积    
Effect of Castration on Gene Expression in Longissimus dorsi Muscle of Huainan Male Pig by Transcriptome Analysis
WANG Jing1, HUA Liushuai1, CHEN Junfeng1, ZHANG Jiaqing1, REN Qiaoling1, BAI Hongjie1, GUO Hongxia1, XU Zhaoxue1, XING Baosong1, BAI Xianxiao1, CAO Hai2     
1. Henan Key Laboratory of Farm Animal Breeding and Nutritional Regulation, Institute of Animal Husbandry and Veterinary Science, Henan Academy of Agricultural Sciences, Zhengzhou 450002, China;
2. Henan Xing Rui Agricultural and Animal Husbandry Technology Co., LTD, Xinyang 465550, China
Abstract: The aim of this study was to detect the influence of castration on gene expression in porcine longissimus dorsi (LD) muscle and reveal the molecular regulation mechanism of castration on porcine meat quality. Six healthy Huainan male pigs were divided into 2 groups, castrated group and intact group (control group), the LD muscle was collected when their body weight achieved 130 kg (300-315 days of age). The LD muscle transcriptomes of pigs in the 2 groups were analyzed by Illumina HiSeq 2000 high-throughput RNA sequencing system. The differentially expressed genes were screened using DESeq software. All of differentially expressed genes were annotated by using GO and KEGG databases. The results showed that 83 613 018 and 83 746 508 reads were obtained from the castrated and intact groups, respectively. And the ratios of these reads to porcine reference genome (Sscrofa10.2) were 73.78% and 74.09%, respectively. Compared with intact pigs, there were 935 differentially expressed genes in LD muscle of castrated pigs, of which 503 up-regulated genes and 432 down-regulated genes (P < 0.05 and|log2(Fold Change)|>0.8). The top up-and down-regulated differentially expressed genes associated with muscle development were collagen type XI alpha 1 chain (COL11A1) and angiotensin Ⅱ receptor type 1 (AGTR1) in castrated group vs intact group. Similarly, the top up-and down-regulated differentially expressed genes associated with lipid metabolism were stearoyl-coA desaturase 1 (SCD-1) and phospholipase A2 group XVI (PLA2G16). GO and KEGG analysis showed that the differentially expressed genes were mainly associated with lipid metabolism (like insulin, lipolysis, fatty acid elongation) and muscular development (such as adrenergic signaling, muscle contraction, AMPK signaling pathway) related pathways. In conclusion, castration can increase male pigs' lipid deposition, change muscular development and meat quality. There were 12 differentially expressed genes, including SCD-1, hormone-sensitive lipase(HSL), glucose transporter type 4(GLUT4), serine/threonine kinase(AKT), mitogen-activated protein kinase(MEK), and so on. It was reported that these genes could regulate hormone secretion, muscle development and lipid deposition. These genes were significantly associated with the alteration of meat quality in porcine after castration, so they were worth of further functional verification.
Key words: Huainan pig     meat quality     transcriptome sequencing     muscle development     lipid deposition    

目前,去势已经成为规范化生产流程中不可缺少的一步。去势不仅可以减少雄性激素的分泌,减少公猪膻味,实现环境友好型生产。其次,去势后公猪性情更加温和,可为猪场减少不必要的经济损失[1]。睾酮作为性激素,除了与生殖相关,还广泛参与机体各项生命活动,例如蛋白质合成、脂质代谢、肌肉生长、神经系统发育、骨量维持、心血管系统保护等[2]。目前,关于去势对动物肉质性状影响的研究已有很多,有研究表明,去势导致小尾寒羊羔羊日增重降低26.9%(P < 0.05),耗料增重比上升50.5%(P < 0.05),肌肉中亚油酸、花生四烯酸等6种脂肪酸含量显著降低(P < 0.05)[3]。杜环利等[4]发现,公猪去势后背膘厚及肌内脂肪含量显著提高(P < 0.05)。姜锦鹏等[5]研究表明,去势公猪日粮添加脱氢表雄酮后平均背膘厚显著降低(P < 0.05),眼肌面积显著增大(P < 0.05),肌肉pH显著提高(P < 0.05),但对胴体质量、屠宰率、胴体长均无显著影响(P>0.05)。关中黑猪去势试验的结果证明,去势后公猪生长速度快、脂肪沉积能力强、肉质优于母猪,但是瘦肉率不如母猪[6]。同样,去势对公鸡肌肉品质也有影响,试验证明,去势明显增强贵州黄鸡脂肪沉积能力(P < 0.05),显著增加其宰前体重、腿肌重、失水率和肌肉嫩度(P < 0.05),而且其肌肉pH更接近优质肌肉的标准[7]。在牛上也有类似的结果,郭同军等[8]试验证明,去势能提高17~21月龄西门塔尔牛肌红蛋白含量(P < 0.01),下调背最长肌pH(P < 0.01),同时提高大理石纹等级,增加牛肉脂肪、干物质、油酸和不饱和脂肪酸含量(P < 0.05)。综上结果来看,关于去势对肉质性状的研究结果比较一致,不同动物去势后脂质沉积均增加,虽然所测肉质性状的具体指标不太相同,但是脂质沉积增多,肌内脂肪相关的肉质性状指标都有所改善。

目前已有一些关于去势对肉质性状分子调控机制的研究,例如蔡兆伟等[9]通过荧光定量PCR技术,比较了去势和非去势公猪脂肪组织中脂肪代谢关键基因的表达量,结果提示,去势公猪体内脂肪沉积增加可能与过氧化物酶体增殖物激活受体(PPARγ)、脂肪酸合成酶(FAS)、乙酰辅酶A羧化酶(ACC)、硬脂酰辅酶A去饱和酶1(SCD-1)和苹果酸酶1(ME1)基因的表达量上调有关。Bai等[10]以皮特兰×长白二元杂交猪为研究对象,发现去势后皮下脂肪中上调miRNAs的靶基因集中在抑制脂肪形成的MAPK、Wnt通路;其中miR-30a表达量显著上升,其靶基因是雄激素受体(AR),推测去势后的脂肪沉积可能与miR-30a对AR表达量的抑制作用有关;部分下调miRNAs的靶基因是转化生长因子β1受体2(TGFBR2),该基因是脂肪沉积相关TGF-β通路的重要组成部分。在长白猪的试验中,与非去势公猪相比,去势公猪皮下脂肪组织有18个miRNAs表达量变化达到显著水平,功能分析发现它们与脂肪酸代谢相关[11]。Cai等[12]通过比较去势和非去势公猪背最长肌miRNA表达谱,结果发现7个差异表达miRNAs,这7个miRNAs与骨骼肌收缩和脂质代谢相关。

为了验证去势对我国地方猪肉质的影响,笔者之前以河南优良地方品种—淮南黑猪为研究对象,通过屠宰试验分析了去势对淮南黑猪公猪肉质性状的影响。结果发现,与同等体重非去势公猪相比,去势后淮南公猪脂肪细胞平均横截面积显著增大(P < 0.05),背膘厚和肌内脂肪含量显著上调(P < 0.05),瘦肉率和眼肌面积显著减小(P < 0.05),大理石纹等级升高[13-14]。为了进一步研究去势对猪肉质性状影响的分子调控机制,本研究以淮南猪为试验对象,通过高通量测序技术分析去势对淮南公猪背最长肌基因表达的影响,筛选寻找去势调控肉质性状的关键基因,旨在为揭示去势影响公猪肉质性状的分子机制提供相应的理论依据。

1 材料与方法 1.1 供试动物及样品采集

为了比较去势对猪肌肉发育和肌内脂肪沉积的影响,选择体重相近的淮南仔公猪3对,每对来源于同一窝,出生1周以内每对同窝仔猪随机选择1头手术去势(共3头去势),剩余3头在腹部同样位置切口做伪处理,保证2组所受手术应激一致。按照商品猪的饲养标准和流程进行饲养,体重达到130 kg屠宰。屠宰30 min内采集背最长肌样品,剪成黄豆大小,置于冻存管中液氮保存。

1.2 猪背最长肌RNA提取、文库构建及质检

使用TRIzol reagent(Invitrogen,Carlsbad,USA)提取猪背最长肌总RNA。然后使用琼脂糖凝胶电泳、Nano-Drop ND-2000分光光度计(NanoDrop products,Wilmington,USA)和Agilent 2100生物分析仪(Agilent Technologies,Massy,France)检测所提RNA的完整性、纯度和质量。使用Qubit RNA BR试剂盒检测所提RNA浓度。所有样品所提RNA的RIN值(RNA integrity number)大于8的用于高通量测序建库。

1.3 猪背最长肌转录组高通量测序及生物信息学分析

利用RiboZeroTM(Epicentre,Madison,WI,USA)试剂盒去除所提RNA中的rRNA,然后使用片段化缓冲液将RNA片段化成150~200 bp。使用SuperScript Ⅱ reverse transcriptase(Life Techonologies, Saint Aubin, France)合成双链cDNA,末端修复和磷酸化后,然后使用Agencourt AMPure® XP beads (Beckman Coulter, Villepinte, France)纯化。经Agilent 2100 Bioanalyzer定量分析后,将cDNA文库在Illumina HiSeq 2000平台测序。

1.4 猪背最长肌基因结构优化和新转录本的预测

测序所得序列通过Cufflinks软件进行转录本组装,将小于180 bp的过滤掉,然后与已公布的参考基因组(Sscrofa10.2)注释信息进行比对。其中重叠基因2个200 bp的侧翼区域,通过2-way比较法进行blast,重叠度大于等于0.9的就是基因结构优化的结果,会延长已注释基因的5′和/或3′端序列。位于基因间隔区的转录本,且距离现有注释基因200 bp以上的可用于候选新转录本的分析。使用Coding Potential Calculator (CPC,http://cpc.cbi.pku.edu.cn/)对新转录本的编码能力进行预测。

1.5 猪背最长肌基因SNPs以及可变剪接

基于测序数据,使用ASprofile软件对所得样品中的可变剪切事件进行分析和统计,包括单内含子保留(intron retention,IR)、模糊边界单内含子保留(approximate IR,XIR)、多内含子保留(mulit IR-MIR)、模糊边界多内含子保留(approximate mulit-IR,XMIR)、可变末尾外显子(alternative last exons,TTS)、可变第一外显子(alternative first exons,TSS)、单外显子跳跃(skipped exons,SKIP)、模糊边界单外显子跳跃(approximate SKIP,XSKIP)、多外显子跳跃(mulit-exon SKIP,MSKIP)、可变5′或3′端剪切(alternative exon ends,AE)以及模糊边界5′或3′端剪切(approximate AE,XAE)。同时还统计了每个样品的单核苷酸多态性(single nucleotide polymorphisms,SNP)位点和插入缺失突变(insert-deletion,INDEL)。

1.6 猪背最长肌差异表达基因分析

使用bowtie2(http://bowtie-bio.sourceforge.net/bowtie2/manual.shtml)和eXpress(http://www.rna-seqblog.com/express-a-tool-for-quantification-of-rna-seq-data)软件,用序列相似性比对的方法计算各转录本在各样本中的表达丰度。使用FPKM(fragments Per kb per Million reads)来表示基因的表达量。差异表达基因的筛选条件为P < 0.05且|log2(Fold Change)|>0.8。

1.7 猪背最长肌差异表达基因的GO和KEGG分析

为了进一步分析背最长肌差异表达基因的功能,对这些基因进行GO富集分析,对其功能进行描述。同时利用KEGG数据库对差异表达基因进行通路分析,并用超几何分布检验的方法计算每个通路条目中差异基因富集的显著性。同时围绕激素分泌、肌肉发育和脂质代谢3个方面选取富集到的通路,将这些通路中差异表达基因进行分析,寻找能同时参与调控其中2个或3个方面的基因,利用Cytogrape软件进行绘图。

2 结果 2.1 猪背最长肌转录组测序结果评估及数据初步分析

通过去除核糖体RNA法建库,在Illumina HiSeq 2000平台进行测序,去势组和非去势组共得到1.70亿raw reads,去除空读、低质量和接头序列后共获得1.67亿clean reads,去势和非去势淮南公猪分别获得83 613 018和83 746 508条可用读段,序列的有效率高达98.00%以上,2组样品GC含量均为52.50%。

2.2 猪背最长肌测序结果基因组比对分析

与猪的基因组序列进行比对,去势组和非去势组的总比对率分别为73.78%和74.09%,比对率高说明测序所获得的序列可信度高,可应用于进一步的分析。其中有大于13.00%的序列在基因组上有多个比对位置,60.00%左右的序列在参考基因组上具有唯一的比对位置。其中比对到“+”链和“-”链的均为30.00%左右,比对到2个外显子上的序列,也就是junction reads为24.00%~25.43%,而整段比对到外显子的序列为35.00%左右。测得数据已提交到NCBI数据库,登录号为GSE87680。

2.3 猪背最长肌基因结构优化和新转录本的预测

通过比对分析,发现共有486个基因的507个区域得到了基因结构优化,其中5′端191个,3′端316个,其中位于正链239个,位于负链268个。除了Y染色体,所有染色体都有分布,其中1号染色体最多共60个,16号染色体最少仅5个,X染色体有21个。其中5′端扩展序列范围为11~289 272 bp,平均27 963 bp,而3′端扩展序列为64~370 812 bp,平均39 163 bp。

本研究中,2个样品测序共发现11 377个基因的12 550个新转录本,其中6 408个位于“-”链,6 142个位于“+”链,平均外显子数2个,转录本长度范围为180~36 480 bp,平均长度为1 541 bp,经编码潜能软件预测分析,其中3 863个具有编码潜能,剩余8 687个为非编码序列。

2.4 猪背最长肌基因SNPs检测以及可变剪接分析

2个样品共检测到16 049个基因的357 992个SNPs位点,位于所有染色体和线粒体,其中A-G和C-T转换型SNP位点占大多数,约为72.56%,二者具体含量差异不大,颠换型SNP位点(G-C、T-G、A-C、A-T)仅占27.44%,其中A-C数量最多,A-T最少。还检测到8 866个基因的30 172个INDEL,这些INDEL位点位于所有染色体。

图 1可以看出,2个样品中均包含12种可变剪切模式,而且趋势相似,均为TSS(15 561个和14 770个)最多,TTS次之(14 567个和13 924个),接着是SKIP(4 310个和4 908个),剩余几种可变剪切模式均极少,其中XMIR仅为10个和4个。

图 1 去势组和非去势组背最长肌可变剪切类型及数量 Fig. 1 Types and number of alternative splicing events of longissimus dorsi muscle in castrated and intact groups
2.5 去势和非去势猪背最长肌差异表达基因分析

去势组和非去势组共有表达基因16 170个,去势组特有1 468个基因,非去势组特有基因1 071个。使用FPKM值来表示基因的表达量,总体来看去势组和非去势组全部基因表达量变化不明显(图 2)。两个样品中表达量最高的5个基因分别是肌球蛋白重链1(MYH1)、骨骼肌α-肌动蛋白1(ACTA1)、肌球蛋白重链2(MYH2)、线粒体细胞色素c氧化酶第Ⅰ亚基(COX1)和肌酸激酶基因(CKM)。同一基因在2个样品中的表达量符合|log2(Fold Change)|≥0.8,且P<0.05即为差异表达基因,与非去势组相比,去势组背最长肌共有503个上调基因,432个下调基因。其中与肌肉发育相关的上下调最显著的20个基因见表 1,下调的10个基因log2(Fold Change)平均值为-2.27,下调最显著的是血管紧张素Ⅱ-1型受体(AGTR1)。上调的10个基因log2(Fold Change)平均值为2.60,上调最显著的是胶原蛋白XI型α1(COL11A1)。这些基因主要与肌细胞增殖、分化和迁移,骨骼肌的发育、肥大、收缩、肾上腺素能信号等有关。而与脂肪代谢相关的显著变化的前20个基因如表 2所示,下调最显著的是脂肪特异性磷脂酶A2(PLA2G16),而上调最显著的是SCD-1。这些基因参与调控糖酵解/糖异生、脂肪酸的合成、延长、降解、甘油酯代谢等多个方面。

图 2 去势组和非去势组背最长肌基因FPKM值箱线图 Fig. 2 The box plots of FPKM value of genes expressed in longissimus dorsi muscle in castrated and intact groups
表 1 与非去势组相比,去势组与肌肉发育相关上调/下调最显著的20个基因 Table 1 The top 20 up- or down-regulated differentially expressed genes related to muscle development in castrated group compared to intact group
表 2 与非去势组相比、去势组与脂质代谢相关上调/下调最显著的20个基因 Table 2 The top 20 up- or down-regulated differentially expressed genes related to lipid metabolism in castrated group compared to intact group
2.6 去势和非去势猪背最长肌差异表达基因GO富集和KEGG通路分析

对去势组和非去势组差异表达基因进行GO富集分析,图 3显示的是显著富集的前20条通路。其中蛋白质在线粒体的定位、代谢过程、能量动态平衡等属于能量代谢方面的通路。去势组和非去势组在脂质代谢方面有较大差异,例如对于糖脂代谢非常重要的胰岛素受体、脂肪酸合成以及甘油三酯合成等信号通路。去势也引起肌肉发育方面的差异,尤其是心肌收缩通路相关基因。KEGG富集分析结果与GO分析结果一致(图 4),差异表达基因极显著富集于能量代谢、脂肪沉积和肌肉发育相关通路,其中能量代谢相关通路主要涉及蛋白质、氨基酸和淀粉以及糖的代谢,而脂肪沉积相关通路包括脂解的调控、脂肪酸的延长,而心肌细胞肾上腺素能信号、心肌收缩、磷酸腺苷活化的蛋白激酶(AMPK)等均为肌肉发育相关通路。

图 3 差异表达基因GO富集分析 Fig. 3 The GO analysis of the differentially expressed genes
图 4 差异表达基因KEGG富集分析 Fig. 4 The KEGG analysis of the differentially expressed genes
3 讨论

研究表明,去势通过激素分泌、脂质代谢和肌肉发育等3个方面间接调控肉质性状[1]。本研究比较了去势猪和非去势猪背最长肌的基因表达差异,如图 5所示,在这3个方面确实富集到较多差异表达基因。在激素方面,去势最直接影响的是睾酮分泌量,睾酮可抑制脂蛋白酯酶(LPL)活性,促进脂肪分解,抑制脂肪沉积,达到促进能量消耗的作用[15]。本研究结果中,差异表达基因显著富集于雌激素代谢通路而不是雄激素代谢通路。已有研究表明,雄激素去芳香化可转化为雌激素[16],由此推断去势后雄激素通过雌激素代谢通路间接发挥作用。甲状腺激素在脂质代谢方面也很重要,其通过促进生长激素合成、增强脂肪组织对儿茶酚胺和胰高血糖素的敏感性以及拮抗胰岛素的作用,达到促进脂质分解的目的[17]。在本研究中,与非去势组相比,甲状腺激素信号通路大部分基因表达量下调,符合去势后脂质沉积增多的现象。与脂肪代谢密切相关的还有胰岛素,其通过刺激脂肪细胞吸收利用葡萄糖,促进脂肪合成,在本研究中, 胰岛素信号通路中关于葡萄糖转运、脂肪酸合成、糖酵解、糖原合成、脂解等多个方面的基因表达量均有显著变化。其中葡萄糖转运体4(GLUT4)是调控肌肉葡萄糖代谢的关键因子[18],研究表明,在C2C12成肌细胞中过表达广西巴马小型猪Glut4,能增加细胞的葡萄糖摄入量[19]。在莱芜猪群体中,Glut4表达量升高能够促进肌内脂肪含量增加[20]。这些结果提示,去势能通过调控雌激素、甲状腺激素以及胰岛素等的分泌参与机体能量代谢的调控。

灰色.生物过程;红色.信号通路;绿色.基因 Gray. Biological processes; red. Signaling pathways; green. Genes 图 5 与激素、脂质沉积和肌肉发育相关的差异表达基因及通路 Fig. 5 The differentially expressed genes and their pathways associated with hormone, lipid deposition and muscle development

去势后脂肪沉积增多,本研究中, 差异表达基因主要富集在脂质代谢方面,这些基因主要与脂肪酸、甘油酯以及葡萄糖等代谢相关。KEGG富集分析中,与脂质代谢最相关的是淀粉和糖代谢通路,该通路中磷酸二酯酶1(ENPP1)等6个基因在去势组的表达量显著上调,已有研究表明,ENPP1尤其是其K121Q多态位点与二型糖尿病[21]和青少年肥胖[22]等糖脂代谢相关疾病显著相关。其次富集到的是脂肪酸延长通路,该通路中乙酰辅酶A酰基转移酶2(ACAA2)在去势后表达量显著下调,已知该基因主要参与调控脂肪酸的β-氧化和胆固醇的生物合成[23]。已有研究表明,鹅在填饲后,肝中甘油三酯的合成速度急速增加,肝中脂肪酸的β-氧化代偿性增加,ACAA2表达量也显著上调,但随着填饲的持续,增多的脂肪酸氧化后产生的活性自由氧导致线粒体损伤,ACAA2表达量又下降[24]。由此推断,去势后ACAA2表达量降低的原因可能与此类似。值得关注的还有PPAR信号通路(图 6),SCD-1是单不饱和脂肪酸合成的限速酶,通过调控饱和脂肪酸和单不饱和脂肪酸的比例来影响机体的脂质代谢[25]。本研究中,去势促进SCD-1表达量上调,说明去势能够通过上调SCD-1促进脂肪酸的合成进而增加肌内脂肪沉积。围脂滴蛋白(PLIN1)是调控脂解的重要蛋白,为定位于脂滴表面的高磷酸化蛋白,既可通过阻止脂肪酶与脂滴的接触降低脂解,也可促进激素刺激的脂解[26]。去势组PLIN1表达量是非去势组的20倍左右,说明去势组脂滴都被保护起来用于能量储存。研究表明,雌激素能够通过抑制PLIN1减少肝细胞内的脂质沉积[27],去势后雄性激素分泌减少,进而芳香化的雌激素也变少,所以其对PLIN1的抑制作用也降低了,这与去势后脂肪沉积增多相符。脂肪酸转位酶(FAT/CD36)是油酸、棕榈油酸和棕榈酸这3种长链脂肪酸转运的主要受体,而长链脂肪酸的摄取是细胞利用脂肪酸的关键环节[28]。去势后FAT/CD36表达量显著下调,说明睾酮表达量下降以后,机体整体能量合成增多,能量利用减少,所以转运体也相应减少。

红色.上调;绿色.下调。下图同 Red. Up-regulation; green. Down-regulation. The same as below 图 6 PPAR信号通路 Fig. 6 PPAR signaling pathway

去势后肌肉发育也发生了改变,本研究中,AMPK信号通路被显著富集(图 7),AMPK表达量显著下调。已知AMPK既能通过调控骨骼肌脂肪酸的氧化,调控骨骼肌甘油三酯的沉积[29],又能通过促进GLUT4[18]、糖酵解(磷酸化激活6-磷酸果糖激酶-2,PFK-2)[30]和抑制糖原合成(Gs)3个方面来调控骨骼肌的能量代谢,此外AMPK/SIRT1/PGC-1α通路还参与调节骨骼肌的质量和骨骼肌纤维类型转化[31]。去势后AMPK表达量显著降低,肌肉的增殖分化、能量代谢、肌纤维类型、肌内脂肪沉积等都随之发生改变。本研究富集到关于肌肉发育的信号通路还有RAP1[32]、MAPK[33]和PI3K-AKT[34]等。

图 7 AMPK信号通路 Fig. 7 AMPK signaling pathway

去势后的差异表达基因主要富集于激素分泌、肌肉发育和脂质代谢3个方面,为了筛选去势调控肉质性状的关键靶基因,围绕这3个方面将显著富集的相关通路中的差异表达基因进行统计分析,筛选出12个能同时调控其中2个或3个方面的差异表达基因。这12个基因除了前面已经讨论过的SCD-1、GsAktMEKGlut4等,值得关注的还有HSL和促分裂原活化蛋白激酶(MAPK)信号转导通路,HSL与猪的肌内脂肪含量显著相关[35],MAPK通路能够抑制猪骨骼肌成肌细胞的分化[36]。推测这12个基因在去势相关的肉质性状调控中起到了关键作用,值得进一步在细胞水平进行功能验证。

4 结论

综上所述,去势后猪背最长肌共有935个差异表达基因,其中上调基因503个,下调基因432个。GO和KEGG富集分析发现,这些差异表达基因主要与激素分泌、脂肪沉积和肌肉发育相关。综合分析,得到12个能同时参与调控这几个方面的差异表达基因,这些基因可能参与调控了去势导致的肉质性状改变,可作为进一步功能验证的关键靶点。

参考文献
[1] 蔡兆伟.阉割影响猪脂肪沉积的分子机理研究[D].杭州: 浙江大学, 2011.
CAI Z W.Study on molecular mechanisms of the effect of castration on fat deposition in male pigs[D].Hangzhou: Zhejiang University, 2011.(in Chinese) http://www.wanfangdata.com.cn/details/detail.do?_type=degree&id=Y2048417
[2] 张阳海, 李永, 曹迪, 等. 睾酮对动物生殖和生长发育影响的研究进展[J]. 家畜生态学报, 2018, 39(1): 1–7.
ZHANG Y H, LI Y, CAO D, et al. Research progress on effect of testosterone on the animal reproduction and growth[J]. Acta Ecologae Animalis Domastici, 2018, 39(1): 1–7. DOI: 10.3969/j.issn.1673-1182.2018.01.001 (in Chinese)
[3] 李宇, 赵张晗, 段春辉, 等. 去势对小尾寒羊公羔育肥及肉品质的影响[J]. 中国畜牧杂志, 2018, 54(12): 103–108.
LI Y, ZHAO Z H, DUAN C H, et al. Effect of castration on fattening and meat quality of male Small-tail Han lambs[J]. Chinese Journal of Animal Science, 2018, 54(12): 103–108. (in Chinese)
[4] 杜环利, 尚伟伟, 周振金, 等. 去势对公猪生长性能和胴体品质的影响[J]. 畜牧与兽医, 2014, 46(6): 27–32.
DU H L, SHANG W W, ZHOU Z J, et al. Effects of castration on growth performance and carcass qualities in the boar[J]. Animal Husbandry & Veterinary Medicine, 2014, 46(6): 27–32. (in Chinese)
[5] 姜锦鹏, 宁康健, 应如海. 脱氢表雄酮对去势公猪胴体品质与抗氧化能力的影响[J]. 安徽科技学院学报, 2018, 32(4): 1–7.
JIANG J P, NING K J, YING R H. Effects of dietary dehydroepiandrosterone supplementation on carcass quality and anti-oxidative capability in barrows[J]. Journal of Anhui Science and Technology University, 2018, 32(4): 1–7. (in Chinese)
[6] 谢宝财.性别对关中黑猪生长性能、胴体性状及肉品质的影响[D].杨凌: 西北农林科技大学, 2017.
XIE B C.Effects of gender on growth performance, carcasstraits and meat quality in Guanzhong black pigs[D].Yangling: Northwest A & F University, 2017.(in Chinese) http://cdmd.cnki.com.cn/Article/CDMD-10712-1017091394.htm
[7] 林家栋, 张福平, 龚俞, 等. 中草药去势对公鸡肌肉品质的影响[J]. 贵州畜牧兽医, 2018, 42(2): 9–13.
LIN J D, ZHANG F P, GONG Y, et al. Herbal castration on rooster meat quality[J]. Guizhou Journal of Animal Husbandry & Veterinary Medicine, 2018, 42(2): 9–13. DOI: 10.3969/j.issn.1007-1474.2018.02.003 (in Chinese)
[8] 郭同军, 臧长江, 王连群, 等. 去势对西门塔尔牛牛肉品质的影响[J]. 草业科学, 2017, 34(1): 152–160.
GUO T J, ZANG C J, WANG L Q, et al. Effect of castration on meat quality of simmental cattle[J]. Pratacultural Science, 2017, 34(1): 152–160. (in Chinese)
[9] 蔡兆伟, 陈民利, 张立凡, 等. 去势对猪脂肪沉积及其关键基因表达的影响[J]. 农业生物技术学报, 2013, 21(10): 1203–1209.
CAI Z W, CHEN M L, ZHANG L F, et al. Effects of castration on fat deposition and its related gene expression in adipose tissues of male pigs (Sus scrofa)[J]. Journal of Agricultural Biotechnology, 2013, 21(10): 1203–1209. DOI: 10.3969/j.issn.1674-7968.2013.10.009 (in Chinese)
[10] BAI Y, HUANG J M, LIU G, et al. A comprehensive microRNA expression profile of the backfat tissue from castrated and intact full-sib pair male pigs[J]. BMC Genomics, 2014, 15: 47. DOI: 10.1186/1471-2164-15-47
[11] CAI Z W, ZHANG L F, CHEN M L, et al. Castration-induced changes in microRNA expression profiles in subcutaneous adipose tissue of male pigs[J]. J Appl Genet, 2014, 55(2): 259–266. DOI: 10.1007/s13353-014-0194-0
[12] CAI Z W, ZHANG L F, JIANG X L, et al. Differential miRNA expression profiles in the Longissimus dorsi muscle between intact and castrated male pigs[J]. Res Vet Sci, 2015, 99: 99–104. DOI: 10.1016/j.rvsc.2014.12.012
[13] WANG J, HUA L S, CHEN J F, et al. Identification and characterization of long non-coding RNAs in subcutaneous adipose tissue from castrated and intact full-sib pair Huainan male pigs[J]. BMC Genomics, 2017, 18(1): 542. DOI: 10.1186/s12864-017-3907-z
[14] XING B S, BAI X X, GUO H X, et al. Long non-coding RNA analysis of muscular responses to testosterone deficiency in Huainan male pigs[J]. Anim Sci J, 2017, 88(9): 1451–1456. DOI: 10.1111/asj.12777
[15] SHEN M Q, SHI H F. Sex hormones and their receptors regulate liver energy homeostasis[J]. Int J Endocrinol, 2015, 2015: 294278.
[16] 赵冉, 袁中尚, 巩光平, 等. 性激素对脂代谢的影响及其临床应用进展[J]. 山东医药, 2015, 55(37): 98–100.
ZHAO R, YUAN Z S, GONG G P, et al. The effect of sex hormones on lipid metabolism and its clinical application progress[J]. Shandong Medical Journal, 2015, 55(37): 98–100. DOI: 10.3969/j.issn.1002-266X.2015.37.040 (in Chinese)
[17] JUNG K Y, AHN H Y, HAN S K, et al. Association between thyroid function and lipid profiles, apolipoproteins, and high-density lipoprotein function[J]. J Clin Lipidol, 2017, 11(6): 1347–1353. DOI: 10.1016/j.jacl.2017.08.015
[18] SHAMNI O, COHEN G, GRUZMAN A, et al. Regulation of GLUT4 activity in myotubes by 3-O-methyl-D-glucose[J]. Biochim Biophys Acta Biomembr, 2017, 1859(10): 1900–1910. DOI: 10.1016/j.bbamem.2017.06.013
[19] 刘嘉琪, 程晓芳, 周亭亭, 等. 广西巴马小型猪Glut4基因克隆和真核表达[J]. 基因组学与应用生物学, 2016, 35(7): 1690–1695.
LIU J Q, CHENG X F, ZHOU T T, et al. Cloning and eukaryotic expression of Glut4 gene from Guangxi Bama mini-pig[J]. Genomics and Applied Biology, 2016, 35(7): 1690–1695. (in Chinese)
[20] 熊信威.猪12号染色体上增加肌内脂肪含量的因果基因及其因果突变的鉴别与作用机理研究[D].南昌: 江西农业大学, 2017.
XIONG X W.Identification of the causative gene and mutation for increasing IMF content at SSC12 and its mechanism study[D].Nanchang: Jiangxi Agricultural University, 2017.(in Chinese)
[21] HSIAO T J, LIN E. The ENPP1 K121Q polymorphism is associated with type 2 diabetes and related metabolic phenotypes in a Taiwanese population[J]. Mol Cell Endocrinol, 2016, 433: 20–25. DOI: 10.1016/j.mce.2016.05.020
[22] DI J Y, DAI M L, ZHANG Z X. ENPP1 K121Q (rs1044498 C > A) genetic polymorphism confers a high risk of susceptibility to coronary heart disease:a PRISMA-compliant article[J]. Medicine (Baltimore), 2018, 97(27): e11236. DOI: 10.1097/MD.0000000000011236
[23] MILTIADOU D, HAGER-THEODORIDES A L, SYMEOU S, et al. Variants in the 3' untranslated region of the ovine acetyl-coenzyme A acyltransferase 2 gene are associated with dairy traits and exhibit differential allelic expression[J]. J Dairy Sci, 2017, 100(8): 6285–6297. DOI: 10.3168/jds.2016-12326
[24] 王倩倩, 杨彪, 夏丽丽, 等. 鹅乙酰辅酶A酰基转移酶2基因的克隆及其在鹅肥肝形成过程中的表达变化[J]. 畜牧兽医学报, 2016, 47(4): 700–708.
WANG Q Q, YANG B, XIA L L, et al. Cloning of goose acetyl-coenzyme A acyltransferase 2 and its expression pattern in the development of fatty liver[J]. Acta Veterinaria et Zootechnica Sinica, 2016, 47(4): 700–708. (in Chinese)
[25] POLONI S, BLOM H J, SCHWARTZ I V D. Stearoyl-CoA desaturase-1:is It the link between sulfur amino acids and lipid metabolism?[J]. Biology (Basel), 2015, 4(2): 383–396.
[26] SZTALRYD C, BRASAEMLE D L. The perilipin family of lipid droplet proteins:gatekeepers of intracellular lipolysis[J]. Biochim Biophys Acta Mol Cell Biol Lipids, 2017, 1862(10): 1221–1232. DOI: 10.1016/j.bbalip.2017.07.009
[27] 赵志武, 王君实, 马敏, 等. 下调Perilipin 1基因表达对3T3-L1细胞脂解的影响[J]. 中国生物工程杂志, 2016, 36(3): 17–22.
ZHAO Z W, WANG J S, MA M, et al. Effect of down-regulated Perilipin 1 gene expression on lipolysis of 3T3-L1 adipocytes[J]. China Biotechnology, 2016, 36(3): 17–22. (in Chinese)
[28] GLATZ J F C, LUIKEN J J F P. From fat to FAT (CD36/SR-B2):understanding the regulation of cellular fatty acid uptake[J]. Biochimie, 2017, 136: 21–26. DOI: 10.1016/j.biochi.2016.12.007
[29] WU W C, FENG J, JIANG D H, et al. AMPK regulates lipid accumulation in skeletal muscle cells through FTO-dependent demethylation of N6-methyladenosine[J]. Sci Rep, 2017, 7: 41606. DOI: 10.1038/srep41606
[30] SHEN Q W, MEANS W J, UNDERWOOD K R, et al. Early post-mortem AMP-activated protein kinase (AMPK) activation leads to phosphofructokinase-2 and -1 (PFK-2 and PFK-1) phosphorylation and the development of pale, soft, and exudative (PSE) conditions in porcine longissimus muscle[J]. J Agric Food Chem, 2006, 54(15): 5583–5589. DOI: 10.1021/jf060411k
[31] 明鹏飞, 黄莹莹, 董妍丽, 等. LKB1-AMPKα-SIRT1信号通路在奶牛脂肪组织脂代谢中的调控作用[J]. 生物技术通报, 2019, 35(2): 176–181.
MING P F, HUANG Y Y, DONG Y L, et al. Regulation of LKB1-AMPKα-SIRT1 signal pathway in lipid metabolism in the adipose tissue of dairy cows[J]. Biotechnology Bulletin, 2019, 35(2): 176–181. (in Chinese)
[32] LI Q, TENG Y F, WANG J, et al. Rap1 promotes proliferation and migration of vascular smooth muscle cell via the ERK pathway[J]. Pathol Res Pract, 2018, 214(7): 1045–1050. DOI: 10.1016/j.prp.2018.04.007
[33] WILLKOMM L, GEHLERT S, JACKO D, et al. p38 MAPK activation and H3K4 trimethylation is decreased by lactate in vitro and high intensity resistance training in human skeletal muscle[J]. PLoS One, 2017, 12(5): e0176609. DOI: 10.1371/journal.pone.0176609
[34] MATOBA A, MATSUYAMA N, SHIBATA S, et al. The free fatty acid receptor 1 promotes airway smooth muscle cell proliferation through MEK/ERK and PI3K/Akt signaling pathways[J]. Am J Physiol Lung Cell Mol Physiol, 2018, 314(3): L333–L348. DOI: 10.1152/ajplung.00129.2017
[35] ZHAO S M, REN L J, CHEN L, et al. Differential expression of lipid metabolism related genes in porcine muscle tissue leading to different intramuscular fat deposition[J]. Lipids, 2009, 44(11): 1029–1037. DOI: 10.1007/s11745-009-3356-9
[36] WILL K, KUZINSKI J, KALBE C, et al. Effects of leptin and adiponectin on the growth of porcine myoblasts are associated with changes in p44/42 MAPK signaling[J]. Domest Anim Endocrinol, 2013, 45(4): 196–205. DOI: 10.1016/j.domaniend.2013.09.002