楊樹猛，馬登錄，郭淑珍，李保明，馬忠濤，才讓鬧日，牛小瑩，格桂花，趙光平，楊秀蘭，張玉珍，尕旦吉，張海濱

（甘肅省甘南州畜牧科學研究所，合作 747000）

遺傳育種

藏羊脂聯(lián)素基因多態(tài)性及其與產肉性能的相關性分析

楊樹猛，馬登錄，郭淑珍，李保明，馬忠濤，才讓鬧日，牛小瑩，格桂花，趙光平，楊秀蘭，張玉珍，尕旦吉，張海濱

（甘肅省甘南州畜牧科學研究所，合作 747000）

利用高分辨率熔解曲線技術對176頭2周歲的甘肅藏羊（歐拉型、甘加型、喬科型）脂聯(lián)素基因SNPs位點進行檢測,運用GLM模型將檢測到的SNPs位點與部分胴體及肉質性狀的相關性進行了分析。結果表明：在脂聯(lián)素第2外顯子發(fā)現(xiàn)+67G>C突變使編碼氨基酸由谷氨酸突變?yōu)楣劝滨０?；GG、GC基因型個體的宰前活重、胴體重均顯著高于CC型（P<0.05）。說明脂聯(lián)素基因該位點可能是影響藏羊胴體及肉質性狀的主效QTL或與之緊密連鎖,可作為藏羊高檔羊肉生產的候選分子標記。

脂聯(lián)素基因；SNPs；產肉性能；藏羊

脂聯(lián)素（adiponectin）是脂肪細胞高效表達并通過特殊運輸途徑分泌進入血液循環(huán)的一種活性物質［1］，又被稱為ACRP30［2］、ADIPOQ［3］、APM1［4］或GBP28［5］，是迄今為止發(fā)現(xiàn)的唯一與肥胖呈負相關的脂肪細胞因子，具有調控葡萄糖、脂肪酸代謝及生物體能量穩(wěn)態(tài)等重要作用。脂聯(lián)素在生物體內的生理作用通過受體ADIPOR1和ADIPOR2介導［6］。人和哺乳動物脂聯(lián)素基因均包括3個外顯子和2個內含子,而禽類脂聯(lián)素基因只包括2個外顯子和1個內含子［3-4,7］。在人類，脂聯(lián)素與胰島素抵抗、動脈粥樣硬化之間的相關性已被許多研究者所關注，但對動物脂聯(lián)素基因的研究相對較少，且主要集中在豬、牛等家畜［8-10］。研究發(fā)現(xiàn)該基因可能是控制畜禽脂肪性狀的主效基因或與控制該性狀的主效基因相連鎖，可作為畜禽脂肪性狀的候選基因。目前關于羊脂聯(lián)素基因的研究很少，僅見羊脂聯(lián)素基因與繁殖性能相關性的報道［11］。

隨著全球消費者對綠色食品重要性的認識和需求的日益增加，生活在高海拔無污染地區(qū)的甘肅藏羊以生產“綠色肉食品”的獨特優(yōu)勢引起國內外的廣泛關注。藏羊肉中氨基酸豐富，種類齊全，接近FAO模式，蛋白質品質優(yōu)良，肌間脂肪含量適中，礦物質含量較高，肉品質量上乘，屬于典型的高蛋白、低脂肪、營養(yǎng)豐富的動物性食品，符合現(xiàn)今消費者對肉食品的選擇性需求。為合理開發(fā)利用和擴大藏羊產品市場，應在提高其生產性能的基礎上，加大對其肉品質量的研究。因此，本研究擬采用高分

辨率熔解曲線分析檢測脂聯(lián)素基因第2、3外顯子在甘肅藏羊中的多態(tài)性，并與肉質性狀進行相關性分析，以期尋找控制藏羊肉品質的分子標記位點，為甘肅藏羊的分子標記輔助選擇提供理論依據(jù)。

1 材料與方法

1.1 材料

1.1.1 材料 供試羊均為2±0.2周歲母羊，包括歐拉型藏羊53只，來自甘肅省甘南州瑪曲縣河曲馬場；甘加型藏羊57只，來自甘南州夏河縣甘加鄉(xiāng)核心育種戶；喬科型藏羊66只，來自甘南州碌曲縣尕海鄉(xiāng)核心育種戶。血樣采自綿羊頸靜脈（10 mL/只），用肝素鈉抗凝，迅速運回實驗室，-70℃凍存?zhèn)溆谩Ｍ涝浊昂髮γ恐谎蜻M行屠體性狀分析，屠宰后取背最長肌處肉樣，將其分割為100 g左右的小塊裝入潔凈保鮮袋中，迅速冷置于-18℃冰箱冷凍保存，用于肉質品質、風味等分析。測定的主要性狀和指標包括宰前活重、胴體重、GR值、眼肌面積、失水力、剪切力，所有檢測均參照《綿羊和山羊生產性能測定技術規(guī)范》（NYT1236—2006）進行。

1.1.2 主要試劑與儀器 2×Taq PCR Master Mix、血液基因組DNA提取試劑盒均購自天根生化科技（北京）有限公司；LCGreen PLUS飽和熒光染料購自美國Idaho公司；石蠟油購自北京鼎國生物技術有限責任公司；其他常規(guī)藥品均為國產分析純。VeritiTM96-Well Thermal Cycler型PCR儀為美國ABI公司產品；Allegra TM21R離心機為美國Beckman公司產品；Lightscanner TMHR-I為美國Idaho公司產品；NanoDrop ND-2000核酸定量儀為美國Nanodrop公司產品。引物合成及測序由寶生物工程（大連）有限公司完成。

1.1.3 PCR引物 根據(jù)羊ADIPOQ序列（GenBank登錄號：100294598）［11］，用Oligo 6.0設計2對引物，擴增藏羊脂聯(lián)素基因的第2、3外顯子的部分片段。引物序列和擴增片段特征見表1。其中HRM內參核苷酸片段3’端氨基-C6封閉，PAGE級別。

表1 引物序列、產物大小及退火溫度

1.2 方法

1.2.1 基因組DNA的提取及定量 按照血液基因組試劑盒說明提取DNA，1%瓊脂糖凝膠電泳檢測，用NanoDrop ND-2000超微量核酸蛋白檢測儀測定其濃度和純度，PCR前調整樣本DNA終濃度至30 ng/μL。

1.2.2PCR擴增 PCR反應體系：10 mmol/L dNTP 0.2 μL，Taq DNA聚合酶0.2 μL，10×PCR反應緩沖液1 μL，0.01 mmol/L上游引物0.1 μL，0.01 mmol/L下游引物0.1 μL，1×LCGreen PLUS飽和熒光染料1 μL，寡核苷酸內參0.2 μL（高溫寡核苷酸內參與低溫寡核苷酸內參各0.1 μL），模板1 μL（待測模板0.5 μL和野生型模板0.5 μL），加去離子水至10 μL。PCR時于每一體系中加入20 μL的石蠟油，防止液體揮發(fā)。PCR反應程序：95℃預變性5 min，然后94℃30 s，58℃退火30 s，72℃30 s，40個循環(huán)；最后72℃延伸10 min，4℃保存。

1.2.3 HRM檢測 將反應后的8 μL PCR反應產物移至96孔伯樂 PCR平板中，然后加入 10 μL礦物油，1 000 r/min離心。LightScanner儀器從45℃開始，以0.3℃/s的斜率采集熔解曲線，到 98℃結束，用LightScanner Call IT軟件對采集后的曲線進行分析。

1.2.4 測序 選擇基因型不同的樣品，送大連寶生物工程有限公司進行PCR產物測序。

1.2.5 數(shù)據(jù)統(tǒng)計

1）等位基因頻率：按下列公式計算。

式中，Pi為第i個等位基因的頻率，i為純合復等位基因，j1、j2……jn為與i共有的第1到第n個等位基因。

2）基因型頻率：基因型頻率＝基因型個體數(shù)/測定群體總數(shù)。

3）多態(tài)性分析：多態(tài)性信息含量(polymorphism information content,PIC)、觀察雜合度（observed heterozygosiy, Ho）、遺傳雜合度（heterozygosity,He）和有效等位基因數(shù)（efective number ofalleles,Ne）等遺傳參數(shù)用POPGEN32軟件統(tǒng)計分析（Nei Met al，1979；Nei M，et al，1974）。

4）最小二乘法統(tǒng)計分析模型：利用SPSS 13.0軟件分析場地和基因型等因素與肉質性狀的相關性。根據(jù)影響藏羊生長性狀的不同因素，采用以下固定模型，同時根據(jù)實際情況進行取舍。固定模型如下：

y=μ+Farm+Genetype+X+e

其中：y為個體表型記錄；μ為群體平均值；Farm為場際效應；Genetype為標記基因型效應；X為各種二級和二級以上互作效應，如Farm×Genetype等；e為隨機誤差。用SPSS13.0軟件對數(shù)據(jù)進行分析，并用最小二乘法擬合線性模型，對各基因型間生長性狀進行差異顯著性分析，分析不同基因型與宰前活重、胴體重、GR值、眼肌面積、系水力、嫩度等性狀的相關性。

2 結果與分析

2.1 HRM突變掃描檢測

經(jīng)優(yōu)化PCR條件后，將用于HRM突變掃描檢測的PCR產物1%凝膠電泳和HRM檢測。由圖1可見，由引物E2F/E2R、E3F/E3R擴增的PCR產物分別為183 bp、258 bp的特異性條單帶，無雜帶，可進行HRM檢測。經(jīng)HRM檢測，E2F/E2R擴增的PCR產物明顯地分為3種類型，測序鑒定發(fā)現(xiàn)這3種類型分別為CC野生型、GG純合型和GC雜合型3種基因型（圖2）。測序結果顯示，在脂聯(lián)素第2外顯子片段出現(xiàn)一處突變，對應于基因組序列（GenBank登錄號：NC_007299）+67G>C，此處的突變位于編碼區(qū)，并引起了氨基酸的改變，使第22位氨基酸的第一位密碼子G變?yōu)镃，氨基酸由原來的谷氨酸改變?yōu)楣劝滨０贰Ｔ谥?lián)素第3外顯子中未發(fā)現(xiàn)多態(tài)性。

圖1 脂聯(lián)素基因第2、3外顯子引物擴增瓊脂糖凝膠電泳檢測結果

圖2 脂聯(lián)素基因第2外顯子+67G>C位點基因分型結果

2.2 群體遺傳分析

2.2.1 藏羊脂聯(lián)素基因第2外顯子+67G>C位點遺傳分析 對甘肅歐拉型、甘加型、喬科型藏羊脂聯(lián)素基因第2外顯子+67G>C位點遺傳分析，G等位基因為優(yōu)勢等位基因，3個群體基因頻率分別為0.821、0.860、0.864。GG基因型為優(yōu)勢基因型，3個群體基因型頻率分別為0.736、0.737和0.773（表2）。進行Hardy-Weinberg平衡檢驗，僅歐拉型藏羊在 +67G>C位點上處于Hardy-Weinberg不平衡狀態(tài)（P<0.05）。進行基因雜合度、有效等位基因數(shù)、多態(tài)信息含量統(tǒng)計分析，結果如表2所示，可見+67G>C位點在3個群體中均為中度多態(tài)，多態(tài)信息含量歐拉型藏羊最高，純合度甘加型藏羊最高，雜合度歐拉型藏羊最高。

表2 脂聯(lián)素基因第2外顯子+67G>C位點在甘肅藏羊中的遺傳分析

2.2.2 +67G>C位點不同基因型與甘肅藏羊產肉和肉質性狀的關聯(lián)分析 應用建立的最小二乘法線性模型，對甘肅歐拉型、甘加型、喬科型藏羊脂聯(lián)素基因第2外顯子+67G>C位點的不同基因型與宰前活重、胴體重、GR值、眼肌面積、系水力、嫩度等性狀進行方差分析（表3）。由表3可以看出，GG、GC基因型個體的宰前活重、胴體重均值均顯著高于CC型個體（P<0.05），GR值、眼肌面積、失水力、剪切力等性狀在不同基因型個體之間差異均不顯著（P>0.05）。

表3 脂聯(lián)素基因第2外顯子+67G>C位點不同基因型與藏羊產肉和肉質性狀的關聯(lián)分析

3 討論

脂聯(lián)素在脂肪代謝中起著重要作用。本研究在脂聯(lián)素第2、3外顯子檢測的區(qū)域中，僅在第2外顯子中檢測到一處突變，為錯義突變，這與楊彥杰等［10］在秦川牛中的結果一致。Morsci等［12］在編碼區(qū)也僅檢測到一處突變，為同義突變（第59位氨基酸的第3位密碼子T>C）。其余突變均在非編碼區(qū)，而且大部分都在啟動子區(qū)。Morsci等［12］在安格斯牛的啟動子區(qū)也檢測到了8個突變，張良志等［13］在秦川牛的啟動子區(qū)也檢測到了11個突變。因此，在今后的研究中將重點開展脂聯(lián)素啟動子區(qū)多態(tài)性研究。

Botstein等提出了衡量基因變異程度高低的多態(tài)信息含量指標（PIC）。當PIC>0.5時,該位點為高度多態(tài)位點；當0.5>PIC>0.25時,該位點為中度多態(tài)位點；當PIC<0.25時,該位點為低度多態(tài)位點。本試驗對脂聯(lián)素基因的第2外顯子+67G>C位點研究發(fā)現(xiàn),歐拉型藏羊雜合度和多態(tài)信息含量最高,說明歐拉型藏羊在本研究中發(fā)現(xiàn)的多態(tài)處選擇的遺傳潛力最大。本研究結果表明，脂聯(lián)素基因在不同的品種中具有豐富的多態(tài)性,這為進一步分析其與脂肪等性狀的關系,進而確定能否用于甘肅藏羊的脂肪沉積的選擇標記奠定了基礎。

Dai等［14］在豬的外顯子2中發(fā)現(xiàn)一處錯義突變A/G（第59位的纈氨酸被替換為異亮氨酸），在內含子2中也發(fā)現(xiàn)了一處突變（第298位A/G），并通過關聯(lián)分析發(fā)現(xiàn)這兩處突變與豬的部分脂肪沉積性狀和胴體性狀（背膘厚，眼肌面積，瘦肉率，骨的比例，平均脂肪量，肩部背膘厚等）有關聯(lián)。楊彥杰等［10］在脂聯(lián)素基因第2外顯子64 bp處和第3外顯子50 bp處發(fā)現(xiàn)的兩處突變都發(fā)生在編碼區(qū),并且導致氨基酸發(fā)生改變。在對這兩個位點突變與胴體及肉質關聯(lián)性分析的結果表明：脂聯(lián)素基因與秦川牛胸圍、眼肌面積、背膘厚、胴體腿臀圍方面密切相關。本研究進一步將脂聯(lián)素基因第2外顯子+67G>C位點不同基因型與甘肅藏羊產肉和肉質性狀的關聯(lián)分析表明，GG、GC基因型個體的宰前活重、胴體重均顯著高于CC型個體（P<0.05），這與楊彥杰等［10］在秦川牛中的研究結果一致。因此，初步推測脂聯(lián)素基因第2外顯子+67G>C位點可能是控制甘肅藏羊宰前活重和胴體重的一個主效基因或是與之存在緊密遺傳連鎖的分子標記，但由于本研究的樣本數(shù)量較少，因此有必要進行更大樣本量或對其他品種進行分析驗證。

［1］Xie L，O'ReillyCP，Chapes SK，et al.Adiponectin and leptin are secreted through distinct trafficking pathways in adipocytes［J］.Biochem Biophys Acta，2008，1782（2）:99-108.

［2］Scherer P E，Williams S，F(xiàn)oglianoM，et al.Anovel serumprotein similar to C1q，produced exclusivelyin adipocytes［J］.J Biol Chem，1995，270（45）：26746-26749.

［3］Hu E，Liang P，Spiegelman B M.AdipoQ is a novel adipose-specific gene dysregulated in obesity［J］.J Biol Chem，1996，271（18）：10697-10703.

［4］Maeda K，Okubo K，Shimomura I，et al.cDNA cloning and expression of a novel adipose specific collagen-like factor，apM1（Adipose most abundant Gene transcript 1）［J］.Biochem Biophys Res Commun，1996，221（2）：286-289.

［5］Nakano Y，Tobe T，Choi-Miura N H，et al.Isolation and characterization of GBP28，a novel gelatin-binding protein purified from human plasma［J］.Journal ofBiochemistry，1996，120（4）：803-812.

［6］Kadowaki T，Yamauchi T.Adiponectin and adiponectin receptors［J］. Endocrine Reviews，2005，26（3）：439-451.

［7］徐國慶，龔道清，儲冬生，等.鵝脂聯(lián)素基因的克隆、序列分析及組織表達［J］.農業(yè)生物技術學報，2008，16（6）：941-946.

［8］張依裕，徐琪，段修軍，等.白羽番鴨脂聯(lián)素基因SNP與肌內脂肪和血清總膽固醇的關聯(lián)分析［J］.中國畜牧雜志，2011（3）：6-10.

［9］俞亞波，王金玉，顧玉萍，等.脂聯(lián)素基因c.249 C>T突變對京海黃雞部分主要經(jīng)濟性狀的影響［J］.西北農業(yè)學報，2011（4）：7-11.

［10］楊彥杰，昝林森，王洪寶，等.秦川牛脂聯(lián)素基因第2外顯子多態(tài)性及其與部分產肉性狀的相關性［J］.西北農林科技大學學報：自然科學版，2009（9）：53-58.

［11］劉重旭，王憑青，張寶云，等.貴州白山羊和古藺馬羊脂聯(lián)素基因多態(tài)性及其與繁殖性能的關聯(lián)研究［J］.中國農業(yè)科學，2011，44（9）：1916-1922.

［12］MorsciNS，SchnabelRD，Taylor JF.Association analysisofadiponectin

and somatostatin polymorphisms on BTA1 with growth and carcass traits in Angus cattle［J］.AnimGenet，2006，37（6）：554-562.

［13］張良志.黃牛脂聯(lián)素基因遺傳特性及在大腸桿菌中的表達研究［D］.陜西楊凌：西北農林科技大學，2008.

［14］Dai L H，Xiong Y Z，Deng C Y，et al.Association of the A-G polymor phism in porcine adiponectin gene with fat deposition and carcass traits［J］.Aian-AustralaianJournalofAnimalScience，2006，19（6）：779-783.

SNPs of Adiponectin Gene and Its Relationship with Meat Traits in Tibetan Sheep

YangShu-meng，Ma Deng-lu，GuoShu-zhen，et al
（Gannan Institute ofAnimal Sciences，Hezuo747000，Gansu，China）

Gansu Tibetan sheep（Euler type，Ganga type，Choco type）（n=176）were used to detect SNPs of adiponectin gene by high-resolution melting（HRM）analysis and toanalyze the correlation ofSNPs with carcass and meat traits usingthe general linear model（GLM）.There was onlyone single nucleotide polymorphism（SNP）identified as SNP+67G>Cin exon2 ofadiponectin，results in a missense mutation from Glu to Gln.Statistical analysis revealed that Tibetan sheep with GG，GC genotype was higher than CC genotype in slaughter weight，carcass weight（P<0.05）.Adiponectin gene was proved to be closely related to carcass and meat traits，which could be used as a candidate molecular marker for production ofhigh-grade meat in Tibetan sheep.

adiponectin gene；SNPs；meat traits；Tibetan sheep

S826.2

A

2095-3887（2014）02-0005-05

10.3969/j.issn.2095-3887.2014.02.001

2013-11-27

楊樹猛（1976-），男，高級畜牧師。

馬登錄（1970-），男，高級畜牧師。