基于RNA-Seq識別山羊肉品質(zhì)候選基因

孟憲然，杜 琛，王 靜，付紹印，鄭竹清，張文廣*，李金泉*

(1.內(nèi)蒙古農(nóng)業(yè)大學(xué)動物科學(xué)學(xué)院,動物遺傳育種與繁殖自治區(qū)重點實驗室，呼和浩特 010018； 2.內(nèi)蒙古醫(yī)科大學(xué)附屬醫(yī)院 婦產(chǎn)科生殖醫(yī)學(xué)中心，呼和浩特 010050)

基于RNA-Seq識別山羊肉品質(zhì)候選基因

孟憲然1#，杜 琛2#，王 靜1，付紹印1，鄭竹清1，張文廣1*，李金泉1*

(1.內(nèi)蒙古農(nóng)業(yè)大學(xué)動物科學(xué)學(xué)院,動物遺傳育種與繁殖自治區(qū)重點實驗室，呼和浩特 010018； 2.內(nèi)蒙古醫(yī)科大學(xué)附屬醫(yī)院 婦產(chǎn)科生殖醫(yī)學(xué)中心，呼和浩特 010050)

通過研究不同年齡和不同性別絨山羊背最長肌之間的差異表達(dá)基因，篩選影響絨山羊肉品質(zhì)的候選基因。本研究基于RNA-Seq技術(shù)對4個絨山羊背最長肌的轉(zhuǎn)錄組進(jìn)行高通量測序，在測序評估的基礎(chǔ)上，對差異表達(dá)基因進(jìn)行篩選、功能注釋和富集分析。本研究通過CLC Genomics Workbench 6.0軟件共找到263個候選基因，其中包含123個高表達(dá)有利基因和140個高表達(dá)有害基因。GO功能注釋結(jié)果顯示，高表達(dá)有利基因主要與骨骼肌的生長發(fā)育、細(xì)胞器的形成和蛋白結(jié)合功能有關(guān)；高表達(dá)有害基因主要與脂質(zhì)代謝、細(xì)胞骨架以及結(jié)合功能有關(guān)。利用KEGG數(shù)據(jù)庫作為參考，這些基因主要參與糖酵解/糖異生、絲裂原活化蛋白激酶、凝血-補體級聯(lián)反應(yīng)和色氨酸代謝等通路中。結(jié)合其他家畜基因組學(xué)研究，ADIPOQ、PDK4和CD36可能作為參與肉品質(zhì)重要的候選基因。該結(jié)果為絨山羊肉品質(zhì)的改善以及候選基因的研究提供了理論依據(jù)。

絨山羊；背最長??；肉品質(zhì)；RNA-Seq；候選基因

羊肉是人類飲食中主要的蛋白質(zhì)來源之一，已被證實具有良好的營養(yǎng)價值，其肉質(zhì)特性與現(xiàn)在消費者對營養(yǎng)肉的需求相一致，具有令人滿意的食用品質(zhì)[1]。羊肉品質(zhì)易受去勢、性別、年齡以及品種等多種因素影響。遠(yuǎn)在殷商時代中國已有家畜去勢的記載，目的是消除家畜的性欲，使家畜的肉質(zhì)得以改善、產(chǎn)量提高。去勢能夠增加肌內(nèi)脂肪沉積，降低瘦肉率、屠宰率和剪切力值，提高肉的嫩度，同時能夠增加血液中FSH和LH水平，降低促性腺激素和雄激素睪酮濃度，減少肉的膻味，進(jìn)而改善公畜肉品質(zhì)[2-3]。鄭中朝等[4]對波雜肉羊的屠宰性能以及肉品質(zhì)進(jìn)行了研究，發(fā)現(xiàn)隨著年齡的增長，成年羊的屠宰性能優(yōu)于羔羊，肌內(nèi)脂肪含量增加，但成年羊肉的剪切力值較高，肉的嫩度變差；母羊的生長優(yōu)勢顯著低于公羊，但成年母羊的肌內(nèi)脂肪含量顯著高于成年公羊，剪切力值顯著降低，成年母羊肉的嫩度明顯優(yōu)于成年公羊。王夢霖等[5]研究表明，隨著年齡的增長，周歲羊的產(chǎn)肉性能明顯優(yōu)于羔羊，剪切力值和纖維直徑顯著增加，嫩度降低；母羊的屠宰性能顯著低于公羊，但隨著肌內(nèi)脂肪沉積增加，肥肉率顯著高于公羊，肉的嫩度增加。由此可知，公羊肉因其肉質(zhì)特性較差，膻味重等特點很難被消費者接受。

轉(zhuǎn)錄組是指特定組織或細(xì)胞在某一發(fā)育階段或功能狀態(tài)下轉(zhuǎn)錄出來的所有基因轉(zhuǎn)錄產(chǎn)物的總和，不僅包括mRNA，而且還包括非編碼RNA，是連接基因組遺傳信息與蛋白質(zhì)組生物功能的橋梁[6-7]。轉(zhuǎn)錄組技術(shù)不僅能夠分析轉(zhuǎn)錄本結(jié)構(gòu)和表達(dá)水平，同時也能發(fā)現(xiàn)未知轉(zhuǎn)錄本和稀有轉(zhuǎn)錄本，精確地識別可變剪切位點，提供更為全面的轉(zhuǎn)錄組信息[8-9]。近年來，隨著高通量測序技術(shù)的不斷發(fā)展，特別是RNA高通量測序技術(shù)的應(yīng)用，使得轉(zhuǎn)錄組學(xué)研究得到了長足的進(jìn)步，并應(yīng)用于不同領(lǐng)域的多方位研究[10]。RNA高通量測序技術(shù)(RNA-Seq)，即利用大規(guī)模測序技術(shù)直接對cDNA序列進(jìn)行測序，產(chǎn)生數(shù)以千萬計的reads，從而使得一段特殊的基因組區(qū)域的轉(zhuǎn)錄水平可以直接通過比對到該基因組區(qū)域的reads數(shù)來衡量[11]。目前主要應(yīng)用于差異基因表達(dá)分析、功能基因挖掘、等位基因特異性表達(dá)、RNA剪接和融合基因轉(zhuǎn)錄子等相關(guān)研究中[12-13]。RNA-Seq對于真核生物的基因表達(dá)調(diào)控、癌癥治療和遺傳育種等方面的研究具有不可估量的潛力。隨著轉(zhuǎn)錄組技術(shù)的迅猛發(fā)展，已經(jīng)得到人類和許多動植物的轉(zhuǎn)錄組信息，這為進(jìn)一步研究生物的功能基因組學(xué)奠定了基礎(chǔ)。

目前關(guān)于山羊肉品質(zhì)的研究還不像豬肉和牛肉等那樣廣泛且深入，利用轉(zhuǎn)錄組學(xué)技術(shù)研究山羊肉品質(zhì)的報道更是少之又少。因此，本研究基于RNA-Seq技術(shù)在基因水平上對山羊背最長肌基因表達(dá)進(jìn)行比較研究，以期篩選出與山羊肉品質(zhì)相關(guān)的候選基因，從而為山羊肉品質(zhì)的改善以及功能基因組學(xué)的研究提供理論依據(jù)。

1 材料與方法

1.1 材料

樣品來自內(nèi)蒙古白絨山羊種羊場。選取成年公羊、成年母羊、斷奶公羔和斷奶母羔的背最長肌為試驗樣品，采集后置于液氮中迅速冷凍，之后置于-80 ℃冰箱中保存?zhèn)溆谩?/p>

1.2 背最長肌RNA提取、建庫及測序

按照Trizol法提取絨山羊背最長肌的總RNA，然后用無RNA酶(RNase)的DNA酶(DNase-I)處理去除總RNA中可能存在的痕量基因組DNA，最后再通過瓊脂糖凝膠電泳檢測總RNA的完整性，用核酸測定儀測定RNA濃度。以1 mg總RNA為起始量，參照Illumina建庫試劑盒進(jìn)行cDNA文庫構(gòu)建，具體操作見Truseq mRNA sample prep kit v2說明書。轉(zhuǎn)錄組測序工作由本實驗室完成。

1.3 轉(zhuǎn)錄組數(shù)據(jù)組裝及候選基因分析

利用CLC Genomics Workbench 6.0軟件對4個絨山羊背最長肌轉(zhuǎn)錄組的原始測序數(shù)據(jù)結(jié)果進(jìn)行過濾和拼接組裝，篩選樣品間的差異表達(dá)基因，然后對各個差異表達(dá)基因進(jìn)行注釋和功能分類。以人類基因組為參考，通過DAVID軟件對絨山羊肉品質(zhì)的候選基因進(jìn)行GO和KEGG功能富集分析。

2 結(jié) 果

2.1 絨山羊背最長肌轉(zhuǎn)錄組的組裝

獲得RNA-Seq的原始數(shù)據(jù)后去除其中低質(zhì)量的數(shù)據(jù)，結(jié)果共獲得52 415 590個reads片段，其中成年公羊、成年母羊、斷奶公羔和斷奶母羔背最長肌轉(zhuǎn)錄組的reads數(shù)分別為12 807 070、12 927 794、13 825 310和12 855 416個，平均讀長為90 bp，共包含4 717 403 100個核苷酸序列信息，GC含量為49.4%。由此可以看出此次轉(zhuǎn)錄組測序的結(jié)果較好，可為后續(xù)的數(shù)據(jù)組裝提供較好的原始數(shù)據(jù)。利用DeNovo方法對52 415 590個reads進(jìn)行拼接共獲得65 329個contigs，N50值為1 059 nt，其長度分布如圖1所示。然后將clean reads映射到contigs上，結(jié)果見表1。

圖1 絨山羊背最長肌轉(zhuǎn)錄組的contig長度分布Fig.1 Distribution of contig length for transcriptome of longissimus dorsi in cashmere goats

表1 絨山羊背最長肌轉(zhuǎn)錄組reads的Mapping統(tǒng)計結(jié)果

Table1 Reads mapping statistics of RNA-Seq oflongissimusdorsiin cashmere goats

樣本名稱Sample片段Reads比對上的Reads百分比/%PercentofReadsmapped核苷酸序列堿基數(shù)Nucleotides成年公羊(A)Adultrams(A)1280707091．011152636300成年母羊(B)Adultewes(B)1292779491．401163501460斷奶公羔(C)Malelambsweaned(C)1382531091．361244277900斷奶母羔(D)Femalelambsweaned(D)1285541690．331156987440合計Total524155904717403100

2.2 候選基因的篩選

首先篩選出絨山羊背最長肌4個轉(zhuǎn)錄組的差異表達(dá)基因，將所有差異表達(dá)基因取并集作為可能影響山羊肉品質(zhì)性狀的候選基因?；虮磉_(dá)量的計算使用RPKM(Reads Per Kb per Million Reads)法,其計算公式：

其中，RPKM為特定生物樣品中某一基因的表達(dá)水平；R為匹配到該基因上的原始reads數(shù)目；N為該樣品測序得到的reads總數(shù)；L為該基因的長度(bp)。本研究利用CLCGenomicsWorkbench6.0軟件共找到474個成年公羊與成年母羊，成年公羊與斷奶公羔，成年母羊與斷奶母羔及斷奶公羔與斷奶母羔背最長肌間的差異表達(dá)基因(P<0.05,FDR<0.001 和FoldChange≥2)，其中成年公羊與成年母羊背最長肌間共201個，成年公羊與斷奶公羔背最長肌間共120個，成年母羊與斷奶母羔背最長肌間共100個，斷奶公羔與斷奶母羔背最長肌間共53個。取并集去重后共篩選到263個絨山羊肉品質(zhì)候選基因，其中將成年公羊與另外3種羊基因表達(dá)量的比值大于2的基因定義為高表達(dá)有害基因，共計140個；表達(dá)量比值小于0.5的基因定義為高表達(dá)有利基因，共計123個，結(jié)果見表2。

表2 絨山羊肉品質(zhì)候選基因

Table2 Candidate genes of meat quality in cashmere goats

分類Classification基因名Genesymbol有利基因BeneficialgenePDK4ADIPOQANKRD1ARID5BGPD1MBBuLA?DQBPGM1FLNCIGFBP?5GAPDHCSRNP1MIDNovar?MHCI?F9NFIL3DDIT4SLNMYBPC2RRADHSP70．1LOC101108462ID3ENO3KLF4PDLIM7JUNBDHRS7CLOC101398826Sdc4PGK1IFRD1NDUFB10ALDH2TUBB6LOC101009376JUNDFOSTPM4CD300LGDNAJA1HSPH1LOC101116336CD52ID1MAFFTMEM52KLF2DUSP1LOC100153946Cib2FOSBPER1GADD45GKLF10ABRALOC101388924EMP1CremTb24C26H8orf4THBDZFP36LOC101121285ASB5CST6MYL1BLA?DQBMYOM2GIGYF2LOC100430941HINT1IGF2RASD1PPP1R3CHSPA6RPS21LOC101364826DQA1CLIC4NR4A3Ovar?DRB1DRASPARCLOC100335469TMP1NR4A2ACTG1MACROD1CD74EGR1LOC101385063SRSF3UQCRHRCAN1TRYPTASE?1BTG2JUNLOC101107541TOB1METRNLIGFBP6C10H5orf13PMP22AK1LOC101108849DHDHALDOAAKAP12OMYHC2AXIRP1TIMP3LOC100851389HSPB8ROCK2CHCHD10C21H14orf2CNN1PAI?1LOC101106131USMG5GADL1BoLA?NCTNFRSF12A有害基因DeleteriousgeneCD36LPIN1ADIPOR2C1QCHRCC3EIF4EBP1NRAPCD63CASQ2MYOM3iglC1QBMAP1LC3BDESNFICIMPA2PSAP089IGKCSRP3Iggamma?1PFN2PRM1HSPB2FBXO32IDO1PSMD4METTL21CCOX1USP13RNF115ACADMCUTCSTAT5BVlambda1bACY1Hspb6LMOD3RBM38ASB2SCPEP1TMEM120ASBDSSYNPOSORT1SESN1SELKSEPP1C25H10orf71EIF5DDIT4LSMTNL1TWF2ATGLAMPD3LOC787392CSDAMYBPHFCER1GWIPI1SMPXTTNLOC781267DOK5THBS4UNC45BPPDPFOVARSHISA5LOC101424NOL3EIF2S3PACSIN3LIMS1FTH1EGLN3LOC101112287UGP2REEP5ANKRD2ARRDC2YPEL3IFITM1LOC100157000COX3WARSSERINC2SCHIP1TACC2MYOZ2LOC101798204OPTNTHOP1PSMD13MYL6BEGLN1BNIP3LOC101105397CST3TIAF1SMTNL2HSPB1GLRXCTSSLOC100153853C1QAASB11MUSTN1LMOD2PCBD1SAP1LOC100685697YBX3PALMDTMEM159PM20D2ACTR1BTXLNBLOC100101238APODSYNPO2SYNPO2LSLC29A2CXCL12HADHALOC101121454GMPRSTAT5AQUAKINGHOMER2MLLT11MAPRE2LOC101119652FHL1SYNGR2SLC43A2CDKN1APDLIM3GATSL2LOC101111776

2.3 候選基因的GO功能分類

GO(Gene Ontology)是一個國際標(biāo)準(zhǔn)化的基因功能分類體系，將基因按照其參與的生物學(xué)過程、細(xì)胞組分和分子功能3個角度進(jìn)行分類，進(jìn)而對基因進(jìn)行限定和描述。本研究利用David軟件中的“GO Functional Annotation Chart”模塊分別對高表達(dá)有利基因和高表達(dá)有害基因進(jìn)行GO功能分類，篩選P<0.05的通路并統(tǒng)計每一類基因的數(shù)量，結(jié)果如圖2所示。在生物學(xué)過程中，高表達(dá)有利基因主要與發(fā)育相關(guān)，主要參與骨骼肌的生長發(fā)育過程。已有研究顯示，ADIPOQ(脂聯(lián)素)能夠加快畜禽肌肉組織中葡萄糖的轉(zhuǎn)運、游離脂肪酸的β氧化和改善胰島素抗性等生物學(xué)作用，其主要是通過活化AMP激活的蛋白激酶(AMPK)和絲裂原活化蛋白激酶(MAPK)途徑激活β氧化和氧化磷酸化過程中相關(guān)基因的表達(dá)來實現(xiàn)的[14-16]。另外，脂聯(lián)素受體能夠通過調(diào)控CD36的表達(dá)影響脂肪酸的跨膜轉(zhuǎn)運。高表達(dá)有害基因主要與調(diào)控和脂質(zhì)代謝相關(guān)。在細(xì)胞組成分中，高表達(dá)有利基因主要參與細(xì)胞器的形成；高表達(dá)有害基因主要與細(xì)胞骨架相關(guān)。在分子功能中，高表達(dá)有利基因和高表達(dá)有害基因主要與結(jié)合功能相關(guān)。

a.高表達(dá)有利基因；b.高表達(dá)有害基因a.Highly expressed beneficial genes； b.Highly expressed deleterious genes圖2 高表達(dá)有利和有害基因GO功能分類Fig.2 Gene Ontology of highly expressed beneficial and deleterious genes

2.4 候選基因的KEGG代謝通路富集分析

在生物體內(nèi)，不同基因相互協(xié)調(diào)行使其生物學(xué)功能，根據(jù)KEGG的分析有助于進(jìn)一步分析差異表達(dá)基因中存在的顯著性富集的代謝通路注釋。本研究使用David軟件中的“KEGG Functional Annotation Chart”模塊分別對高表達(dá)有利基因和高表達(dá)有害基因進(jìn)行KEGG通路富集分析，結(jié)果共得到10條通路，其中顯著富集的通路有7條，主要參與到糖酵解/糖異生、絲裂原活化蛋白激酶、凝血-補體級聯(lián)反應(yīng)、朊病毒病、ErbB受體、色氨酸代謝和系統(tǒng)性紅斑狼瘡等通路中(表3)。糖酵解/糖異生過程對于維持動物機體的能量代謝平衡和骨骼肌的發(fā)育有重要的影響，其中PDK4(丙酮酸脫氫酶激酶4)在這一過程中起著重要的作用。有研究表明，PDK4對于脂肪沉積有促進(jìn)作用，其主要是通過抑制自身的活性激活PDC(丙酮酸脫氫酶復(fù)合體)進(jìn)而增加葡萄糖向脂肪酸轉(zhuǎn)化來實現(xiàn)。

3 討 論

本研究通過RNA-Seq方法共篩選到263個對絨山羊肉品質(zhì)有影響的候選基因，其中包含123個高表達(dá)有利基因和140個高表達(dá)有害基因。其中高表達(dá)有利基因PDK4(丙酮酸脫氫酶激酶4)是PDK/BCKDK蛋白激酶家族成員之一，在調(diào)節(jié)骨骼肌糖代謝與脂肪酸轉(zhuǎn)換的過程中發(fā)揮重要的調(diào)控作用，與動物骨骼肌的生長發(fā)育密切相關(guān)。N.De Jager[17]對牛骨骼肌的基因表達(dá)進(jìn)行研究，發(fā)現(xiàn)西澳大利亞牛骨骼肌中的PDK4表達(dá)量遠(yuǎn)高于新南威爾士牛，并且西澳大利亞牛骨骼肌的大理石紋等級優(yōu)于新南威爾士牛，這表明PDK4表達(dá)量越高，肌內(nèi)脂肪沉積越多，肉品質(zhì)越好。本研究發(fā)現(xiàn)，隨著年齡的增長，母羊骨骼肌中的PDK4表達(dá)量隨之升高，而公羊骨骼肌的PDK4表達(dá)量隨之降低，結(jié)果會導(dǎo)致母羊骨骼肌中沉積的脂肪含量高于公羊，這與王夢霖等[5]研究結(jié)果一致。因此，PDK4可作為影響絨山羊肉品質(zhì)的一個重要候選基因，但需要進(jìn)一步深入研究。

表3 高表達(dá)有利基因和有害基因KEGG代謝途徑分類

Table3 KEGG classification of highly expressed beneficial and deleterious genes

分類Classification通路Pathway基因數(shù)GenenumberP值P?value通路IDPathwayID有利基因BeneficialgenesGlycolysis/Gluconeogenesis64．13E?05hsa00010MAPKsignalingpathway77．68E?03hsa04010有害基因DeleteriousgenesComplementandcoagulationcascades41．88E?02hsa04610Priondiseases33．34E?02hsa05020ErbBsignalingpathway43．44E?02hsa04012Tryptophanmetabolism34．26E?02hsa00380Systemiclupuserythematosus44．76E?02hsa05322

ADIPOQ是脂肪細(xì)胞特異性分泌的一種細(xì)胞因子，在骨骼肌和脂肪組織中調(diào)控生物體的糖代謝、脂肪酸代謝和能量平衡，是迄今為止所發(fā)現(xiàn)的惟一與肥胖呈負(fù)相關(guān)的脂肪細(xì)胞特異性蛋白，屬于2型糖尿病保護(hù)因子，具有增加胰島素敏感性、抗糖尿病和抗動脈粥樣硬化等作用。目前人們對小鼠、人、豬、雞、牛和山羊的ADIPOQ基因都有研究，但主要集中在豬、家禽和牛等動物的基因和其受體的克隆以及基因多態(tài)性分析等方面[18]。戴麗荷等[19]對大白豬、長白豬和梅山豬脂聯(lián)素基因內(nèi)含子進(jìn)行突變檢測，發(fā)現(xiàn)內(nèi)含子2中985位的A/G發(fā)生突變，在其F2資源家系群體中進(jìn)行性狀關(guān)聯(lián)分析表明，A等位基因有利于增加眼肌面積和瘦肉率，減少背膘厚，同時降低了肌內(nèi)脂肪含量，肌肉口感變差，最終影響了肉品質(zhì)。值得注意的是，脂聯(lián)素是通過調(diào)控受體AdipoR1和AdipoR2發(fā)揮其重要的生理作用。AdipoR1主要在骨骼肌中表達(dá)并介導(dǎo)激活A(yù)MPK信號通路；而AdipoR2主要在肝和脂肪組織中表達(dá)并介導(dǎo)激活PPARα信號通路，通過誘導(dǎo)FATP(脂肪酸轉(zhuǎn)運蛋白)和CD36(脂肪酸轉(zhuǎn)位酶)的表達(dá)進(jìn)而促進(jìn)脂肪酸的跨膜轉(zhuǎn)運[20-21]。由此可見，脂聯(lián)素及其受體對畜禽肉品質(zhì)具有重要影響。

CD36是參與脂肪酸跨膜轉(zhuǎn)運和脂肪沉積的重要載體蛋白，在脂肪細(xì)胞和骨骼肌細(xì)胞中參與脂肪酸的運輸[22]。束剛等[23]對黃羽肉雞進(jìn)行CD36重組融合蛋白免疫，發(fā)現(xiàn)首次免疫15 d后便產(chǎn)生較高的抗體效價，而且特異性地抑制公雞腹脂的沉積水平，說明體內(nèi)抗CD36抗體能夠有效地中和內(nèi)源性脂肪細(xì)胞膜CD36，通過部分阻斷轉(zhuǎn)運脂肪酸的功能，從而降低肌內(nèi)脂肪的沉積。L.Love-Gregory等[24]研究表明，CD36可以促進(jìn)血漿中游離脂肪酸的吸收，在脂肪組織中參與脂肪酸的貯存。但血漿中高濃度的脂肪酸會產(chǎn)生炎癥性反應(yīng)而影響胰島素信號使機體產(chǎn)生胰島素抗性，進(jìn)而導(dǎo)致機體糖代謝減少，脂肪代謝增加，最終導(dǎo)致脂肪沉積減少。因此，CD36可以作為肉品質(zhì)研究的一個重要候選基因。

在生物體內(nèi)，各生物學(xué)功能的完成都是基于各個基因間系統(tǒng)地相互協(xié)調(diào)而進(jìn)行的，通過對基因進(jìn)行KEGG通路富集分析，可以更深入地了解該基因的生物學(xué)功能以及與其他基因間的相互作用。KEGG通路數(shù)據(jù)庫儲存了基因功能的相關(guān)信息，通過圖形來表示細(xì)胞內(nèi)的生物學(xué)過程，例如代謝、膜運輸、信號轉(zhuǎn)導(dǎo)和細(xì)胞的生長周期等[25]。本研究中高表達(dá)有利基因主要參與到糖酵解/糖異生和絲裂原活化蛋白激酶信號通路中，其中糖酵解/糖異生過程不僅與骨骼肌的生長發(fā)育和動物機體能量代謝平衡的維持相關(guān)，而且與肌內(nèi)脂肪沉積相關(guān)，這一結(jié)果可以為肉品質(zhì)的研究提供一種全新的視角和理論依據(jù)。

4 結(jié) 論

本研究利用Illumina Miseq高通量測序平臺對絨山羊背最長肌的4個轉(zhuǎn)錄組進(jìn)行測序后，結(jié)果共獲得52 415 590個reads片段，包含了4 717 403 100個核苷酸序列信息，利用DeNovo方法對reads片段進(jìn)行拼接，共得到65 329個contigs，N50值為1 059 nt。利用CLC Genomics Workbench 6.0軟件共篩選到263個與山羊肉品質(zhì)相關(guān)的候選基因，其中包含123個高表達(dá)有利基因和140個高表達(dá)有害基因。通過David軟件分別對高表達(dá)有利基因和有害基因進(jìn)行GO功能注釋和KEGG富集分析發(fā)現(xiàn)，糖酵解/糖異生過程對于維持動物機體的能量代謝平衡和骨骼肌的生長發(fā)育有重要的影響，并且與肌內(nèi)脂肪沉積相關(guān)。結(jié)合其他家畜基因組學(xué)研究，ADIPOQ、PDK4和CD36可能作為參與肉品質(zhì)重要的候選基因。高通量RNA-Seq提供了更為全面的山羊背最長肌全轉(zhuǎn)錄組信息，為山羊肉品質(zhì)的改善及功能基因組研究提供了基礎(chǔ)數(shù)據(jù)和參考。

[1] 李洋靜，謝愷舟，李擁軍，等.海門山羊肉品質(zhì)特性的研究[J].上海畜牧獸醫(yī)通訊，2009(6)：10-12. LI Y J，XIE K Z，LI Y J，et al.Research of meat quality characteristics in Haimen goat[J].ShanghaiJournalofAnimalHusbandryandVeterinaryMedicine，2009(6)：10-12.(in Chinese)

[2] 陳 勝，程廣龍，朱德建，等.能量水平及去勢對雜交公山羊日增重和胴體組成的影響[J].當(dāng)代畜牧，2009(10)：21-24. CHEN S，CHENG G L，ZHU D J，et al.Effect of energy levels and castrated on daily gain and carcass composition of male goat hybrid[J].ContemporaryAnimalHusbandry，2009(10)：21-24.(in Chinese)

[3] 袁亞利，李蛟龍，張衛(wèi)輝，等.免疫去勢和手術(shù)去勢方法對豬肉中氨基酸肌苷酸和脂肪酸含量的影響[J].肉類研究，2011(6)：5-8. YUAN Y L，LI J L，ZHANG W H，et al.A comparison of the effects of two different castration methods on the contents of amino acids，inosine monophosphate and fatty acids in pork[J].MeatReserch，2011(6)：5-8.(in Chinese)

[4] 鄭中朝，張耀強，張 力，等.波雜肉羊屠宰性能及肉品質(zhì)研究[J].中國草食動物，2007，27(5)：21-24. ZHENG Z C，ZHANG Y Q，ZHANG L，et al.Studies on the slaughter performance and meat quality of Boer crossbred goats[J].ChinaHerbivores，2007，27(5)：21-24.(in Chinese)

[5] 王夢霖，雒秋江，楊開倫，等.年齡和性別對陶賽特×小尾寒羊F1代羔羊屠宰性能與肉品質(zhì)的影響[J].中國畜牧獸醫(yī)，2009(2)：152-155. WANG M L，LUO Q J，YANG K L，et al.Effects of age and sex on slaughter characteristics and meat quality of Poll Dorset × Small Tail Han sheep F1 lamb[J].ChinaAnimalHusbandry&VeterinaryMedicine，2009(2)：152-155.(in Chinese)

[6] COSTA V，ANGELINI C，DE F I，et al.Uncovering the complexity of transcriptomes with RNA-Seq[J].JBiomedBiotechnol，2010，2010：1-19.

[7] WANG Z，GERSTEIN M，SNYDER M.RNA-Seq：a revolutionary tool for transcriptomics[J].NatRevGenet，2009，10(1)：57-63.

[8] MORTAZAVI A，WILLIAMS B A，MCCUE K，et al.Mapping and quantifying mammalian transcriptomes by RNA-Seq[J].NatMethods，2008，5(7)：621-628.

[9] 祁云霞，劉永斌，榮威恒.轉(zhuǎn)錄組研究新技術(shù)：RNA-Seq及其應(yīng)用[J].遺傳，2011，33(11)：1191-1202. QI Y X，LIU Y B，RONG W H.RNA-Seq and its applications：a new technology for transcriptomics[J].Hereditas，2011，33(11)：1191-1202.(in Chinese)

[10] 蘭道亮，熊顯榮，林寶山，等.基于微量RNA高通量測序技術(shù)的牦牛MII期卵母細(xì)胞轉(zhuǎn)錄組研究[J].畜牧獸醫(yī)學(xué)報，2014，45(5)：722-732. LAN D L，XIONG X Z，LIN B S，et al.Transcriptome analysis of yak metaphase II (MII) Oocytes by a micro-transcriptome sequencing method[J].ActaVeterinariaetZootechnicaSinica，2014，45(5)：722-732.(in Chinese)

[11] ANSORGE W J.Next-generation DNA sequencing techniques[J].NBiotechnol，2009，25(4)：195-203.

[12] LIPSON D，RAZ T，KIEU A，et al.Quantification of the yeast transcriptome by single-molecule sequencing[J].NatBiotechnol，2009，27(7)：652-658.

[13] 周 華，張 新，劉騰云，等.高通量轉(zhuǎn)錄組測序的數(shù)據(jù)分析與基因發(fā)掘[J].江西科學(xué)，2012，30(5)：607-611. ZHOU H，ZHANG X，LIU T Y，et al.Data processing and gene discovery of high-throughput transcriptome sequencing[J].JiangxiScience，2012，30(5)：607-611.(in Chinese)

[14] KADOWAKI T，YAMAUCHI T.Adiponectin receptor signaling：a new layer to the current model[J].CellMetab，2011，13(2)：123-124.

[15] BERMUDEZ V J，ROJAS E，TOLEDO A，et al.Single-nucleotide polymorphisms in adiponectin，AdipoR1，and AdipoR2 genes：insulin resistance and type 2 diabetes mellitus candidate genes[J].AmJTherap，2013，20(4)：414-421.

[16] FERGUSON J F，PHILLIPS C M，TIERNEY A C，et al.Gene-nutrient interactions in the metabolic syndrome：single nucleotide polymorphisms in ADIPOQ and ADIPOR1 interact with plasma saturated fatty acids to modulate insulin resistance[J].AmJClinNutr，2010，91(3)：794-801.

[17] DE JAGER N.gene expression studies in bovine skeletal muscle：investigations into the effects of hormone growth promotant treatment，environment and tenderness genotype associated with tenderness and marbling[D].Queensland：School of Chemistry & Molecular Biosciences，The University of Queensland，2012.

[18] 陳小玲，黃志清.脂聯(lián)素基因功能的研究進(jìn)展及其應(yīng)用[J].湖北農(nóng)業(yè)科學(xué)，2012，50(23)：4777-4779. CHEN X L，HUANG Z Q.Research advance in function of adiponectin gene and its application[J].HubeiAgriculturalSciences，2012，50(23)：4777-4779.(in Chinese)

[19] 戴麗荷，焦青貞，熊遠(yuǎn)著.豬脂聯(lián)素基因內(nèi)含子2的A/G突變檢測及關(guān)聯(lián)分析[C].第十次全國畜禽遺傳標(biāo)記研討會，2006. DAI L H，JIAO Q Z，XIONG Y Z.Detection of A/G polymorphism in porcine adiponectin gene and association analysis with traits[C].Tenth National Symposium on Animal Genetic Markers，2006.(in Chinese)

[20] VAN B A M，GARNIER A，BECKERS P，et al.Functional adiponectin resistance at the level of the skeletal muscle in mild to moderate chronic heart failure[J].Circulation：HeartFailure，2010，3(2)：185-194.

[21] YAMAUCHI T，NIO Y，MAKI T，et al.Targeted disruption of AdipoR1 and AdipoR2 causes abrogation of adiponectin binding and metabolic actions[J].NatMed，2007，13(3)：332-339.

[22] YU K，SHU G.Fatty acid and transcriptome profiling of longissimus dorsi muscles between pig breeds differing in meat quality[J].IntJBiolSci，2013，9(1)：108.

[23] 束 剛，馮嘉穎，余凱凡，等.FAT/CD36融合蛋白的表達(dá)及其對雞腹脂沉積的特異性調(diào)控[J].中國農(nóng)業(yè)科學(xué)，2009，42(2)：650-656. SHU G，F(xiàn)ENG J Y，YU K F，et al.Prokaryotic expression of FAT/CD36 fusion protein and the specific effects on the deposition of visceral fat in cockerel chicks[J].ScientiaAgriculturaSinica，2009，42(2)：650-656.(in Chinese)

[24] LOVE-GREGORY L，ABUMRAD N A.CD36 genetics and the metabolic complications of obesity[J].CurrOpinClinNutrMetabCare，2011，14(6)：527-534.[25] 杜 琛，付紹印，高鴻雁，等.絨山羊肌內(nèi)脂肪細(xì)胞成熟前后比較轉(zhuǎn)錄組分析[J].畜牧獸醫(yī)學(xué)報，2014，45(5)：714-721. DU C，F(xiàn)U S Y，GAO H Y，et al.Transcriptome analysis of intramuscular preadipocytes and matureadipocyte in cashmere goats[J].ActaVeterinariaetZootechnicaSinica，2014，45(5)：714-721.(in Chinese)

(編輯 郭云雁)

RNA-Seq Approach for Identifying Candidate Genes of Meat Quality in Goats

MENG Xian-ran1#，DU Chen2#，WANG Jing1，F(xiàn)U Shao-yin1，ZHENG Zhu-qing1， ZHANG Wen-guang1*，LI Jin-quan1*

(1.KeyLaboratoryofAnimalGenetics，BreedingandReproductionofInnerMongoliaAutonomousRegion，CollegeofAnimalScience，InnerMongoliaAgriculturalUniversity，Hohhot010018，China；2.DepartmentofObstetricsandGynecology，AffiliatedHospitalofInnerMongoliaMedicalUniversity，Hohhot010050，China)

To screen candidate genes which effect on meat quality of cashmere goats，we explored differentially expressed genes oflongissimusdorsiin cashmere goats with different ages and genders.The high-throughput sequencing for transcriptome oflongissimusdorsiin cashmere goats was performed by RNA-Seq.Differentially expressed genes were selected as candidates genes and enriched based on GO and KEGG database after quality control of RNA-Seq.The 263 candidate genes including 123 beneficial and 140 deleterious genes were selected by CLC Genomics Workbench 6.0 software.The results of GO enrichment showed that the highly expressed beneficial genes were primarily related to growth and development of skeletal muscle，organelle formation and protein binding.The highly expressed deleterious genes were primarily related to lipid metabolic process，cytoskeleton and binding.According to the KEGG database，these genes were mainly involved in glycolysis/gluconeogenesis’mitogen-activated protein kinase，complement and coagulation cascades and tryptophan metabolism.As example，ADIPOQ，PDK4 andCD36 might be important candidate genes involved in meat quality，when combining with other livestock genomics study.These results provided theoretical basis for improving meat quality and studying candidate genes in cashmere goats.

cashmere goat；longissimusdorsi；meat quality； RNA-Seq；candidate gene

10.11843/j.issn.0366-6964.2015.08.004

2014-11-10

國家自然科學(xué)基金(30960246)；內(nèi)蒙古自然科學(xué)基金(2010Zd11)；“內(nèi)蒙古自治區(qū)草原英才”與“內(nèi)蒙古自治區(qū)高等學(xué)校創(chuàng)新團(tuán)隊發(fā)展計劃”(NMGIRT1106)；國家絨毛用羊現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術(shù)體系(CARS-40-05)

孟憲然(1988-)，男，河北承德人，碩士生，主要從事動物遺傳育種與繁殖研究，E-mail：mxr1103@126.com；杜 琛(1986-)，女，博士， E-mail：duchen198607@126.com。 孟憲然和杜琛同為第一作者

*通信作者：張文廣，E-mail：actgnmbi@aliyun.com.cn；李金泉，E-mail:lijinquan_nd@126.com

S827.2

A

0366-6964(2015)08-1300-08