李 梅，安小婭，裴浩宇，白福霞，王 剛，王延濤，于 杰，黨瑞華*

(1.西北農(nóng)林科技大學(xué)動物科技學(xué)院，陜西 楊凌 712100；2.山東東阿阿膠公司，山東 東阿 252299)

生長性狀是畜牧業(yè)生產(chǎn)中關(guān)注的一個重要的經(jīng)濟(jì)性狀，是衡量育種價值及經(jīng)濟(jì)效益的重要指標(biāo)，尋找相關(guān)的分子標(biāo)記可以縮短對生長性狀的育種年限。SNP是指在基因組上單個核苷酸的變異，包括轉(zhuǎn)換、顛換、缺失和插入形成的遺傳標(biāo)記，其數(shù)量很多，多態(tài)性豐富。對于SNP的研究方法有毛細(xì)管電泳、變形高效液相色譜(DHPLC)、質(zhì)譜分型等，質(zhì)譜分型的發(fā)展基于20世紀(jì)80年代產(chǎn)生的一些新的軟電離方式，包括電噴霧電離質(zhì)譜(electrospray ionization-mass spectrometry，ESI-MS)和基質(zhì)輔助激光解吸電離質(zhì)譜(matri-assisted laser desorption ionization-mass spectrometry，MALDI-MS)等，目前國外已有諸多文獻(xiàn)報道生物質(zhì)譜在SNP分型方面的應(yīng)用，主要采用MALDI-TOFMS方法[1]。

Sirtuins家族是廣泛存在于地球所有生命形態(tài)中的具有共同催化結(jié)構(gòu)的古老蛋白質(zhì)，在包括細(xì)菌、針具、酵母菌、瘧原蟲、后生動物、哺乳動物甚至病毒等生命體中都廣泛保守[2-4]。編碼該酶家族的類sirtuin基因包括SIRT1、SIRT2和SIRT3，是組蛋白去乙酰化酶家族中煙酰胺腺嘌呤二核苷酸(NAD)依賴性家族成員，在衰老、代謝和凋亡中發(fā)揮重要作用[5]。SIRT3[沉默配對型信息調(diào)節(jié)3(silent mating-type information regulation 3)]是SIRT3編碼sirtuin家族III類組蛋白去乙酰化酶，定位于線粒體基質(zhì)中[6]，可通過去乙?；€粒體中與ATP生產(chǎn)相關(guān)的酶從而調(diào)劑ATP的生成[7]。2010年Zhao等[8]通過高通量蛋白質(zhì)組學(xué)證實，參與三羧酸循環(huán)、尿素循環(huán)和脂肪酸的β氧化等代謝的蛋白質(zhì)可被乙?；揎?，而在線粒體中大約1/5的蛋白質(zhì)可被乙?；?，結(jié)合SIRT3的活性，SIRT3可通過調(diào)控線粒體蛋白質(zhì)乙?；絹韰⑴c代謝活動。因此，SIRT3可通過調(diào)節(jié)呼吸鏈、三磷酸腺苷(ATP)的產(chǎn)生以及檸檬酸和尿素循環(huán)在調(diào)節(jié)能量代謝和衰老過程中發(fā)揮重要作用，于代謝旺盛的組織如肌肉、褐色脂肪、腎臟、肝臟、腦和心臟等富含線粒體的組織器官中表達(dá)較高[9]。敲除SIRT3基因會導(dǎo)致機(jī)體的線粒體過度乙酰化，造成細(xì)胞能量代謝及有絲分裂過程發(fā)生異常[10]。

王海濤等人[11]研究發(fā)現(xiàn)，SIRT3在骨骼肌中也高表達(dá)，參與能量代謝并影響機(jī)體的衰老、細(xì)胞死亡和腫瘤發(fā)生，且經(jīng)過驗證，適度的運動會增加骨骼肌中的SIRT3含量并進(jìn)一步影響脂肪代謝。有研究[12]表明，SIRT3啟動子轉(zhuǎn)錄因子結(jié)合位點的多態(tài)性可能影響SIRT3的轉(zhuǎn)錄活性，從而改變?nèi)馀＜?nèi)脂肪含量，且SIRT1和SIRT2基因的基因型組合對秦川牛眼肌面積和肌間脂肪有顯著影響。在心臟、肝臟和腎臟等代謝旺盛的組織中，ATP水平較高，且SIRT3也具有較高的表達(dá)水平，由此推測SIRT3的表達(dá)有利于ATP的生成，從而用于機(jī)體的運動代謝等，使得機(jī)體相應(yīng)部位相對發(fā)達(dá)，所以和動物的體尺性狀有一定的關(guān)聯(lián)[13]。

目前關(guān)于SIRT3基因的多態(tài)性研究較少，在動物體尺性狀方面僅在秦川牛上做了相關(guān)研究，而未見關(guān)于驢的相關(guān)報道，鑒于SIRT3基因在能量代謝和脂質(zhì)代謝過程中的重要作用，我們推測其在驢的生長性狀方面可能有直接或間接地作用。因此，本研究擬以SIRT3基因為目標(biāo)基因，采用質(zhì)譜分型方法，分析不同標(biāo)記基因型與德州驢生長性狀之間的關(guān)系，以期為德州驢分子育種提供可參考的基礎(chǔ)數(shù)據(jù)。

1 材料與方法

1.1 試驗材料

血樣采自山東省德州市東阿阿膠公司驢場，選擇同等飼養(yǎng)條件下的30月齡以上健康的214頭公驢、190頭母驢。對每頭驢頸靜脈采血10 mL，ACD抗凝(V(ACD)∶V(血液)=1∶6)，置于-80 ℃保存。同時采集試驗驢的9個生長性狀指標(biāo)(體高、體長、胸圍、胸寬、胸深、管圍、尻長、尻寬、尻高)。

1.2 基因組DNA的提取及檢測

采用常規(guī)的苯酚—氯仿提取法[14]提取德州驢血樣基因組DNA，提取的DNA用TE緩沖液處理后，采用0.5%瓊脂糖凝膠電泳檢測DNA質(zhì)量，用分光光度計檢測D260/D280值，放入-20 ℃保存?zhèn)溆?，在使用時稀釋至50 ng/μL，保存在4 ℃。

1.3 SIRT3基因多態(tài)性鑒定、基因分型

根據(jù)GenBank(http://www.ncbi.nlm.nih.gov/)公布的驢SIRT3序列(GenBank NC_19763和NC_35148)，針對編碼區(qū)利用Primer 5.0軟件設(shè)計引物，引物均由上海生物技術(shù)有限公司合成。

PCR反應(yīng)體系25 μL：含有核酸染料的dNTPs、TaqDNA聚合酶、10×Buffer的Mix 12.5μL，ddH2O 9.5 μL，上游、下游引物各(10 pmoL/μL)1.0 μL，模板DNA(50 ng/μL)1.0 μL。PCR反應(yīng)程序：95 ℃預(yù)變性5 min；95 ℃ 30 s，退火30 s(引物序列及退火溫度見表1)，72 ℃ 30 s，35個循環(huán)；72 ℃延伸10 min；4 ℃保存。PCR擴(kuò)增產(chǎn)物用10 g/L瓊脂糖凝膠電泳檢測。將所有個體DNA稀釋到適宜濃度并送往康普森生物技術(shù)有限公司進(jìn)行質(zhì)譜分型。每個位點隨機(jī)選取20個樣本進(jìn)行PCR擴(kuò)增并送至上海生物工程股份有限公司進(jìn)行測序，測序結(jié)果使用SeqMan軟件進(jìn)行對比分析，進(jìn)行驗證。

表1 SIRT3-1和SIRT3-2位點的引物信息

1.4 統(tǒng)計學(xué)分析

用Excel對基因遺傳變異質(zhì)譜分型結(jié)果進(jìn)行統(tǒng)計分析，利用哈代—溫伯格平衡檢驗軟件對2個突變位點的基因型頻率、等位基因頻率和哈代—溫伯格平衡(HWE)進(jìn)行了分析。同時，運用在線軟件(http:///www.msrcall.com/)估算純合子(Ho)、雜合子(He)、有效等位基因數(shù)(Ne)和多態(tài)信息含量(PIC)。使用SHEsis在線軟件(http://analysis.bio-x.cn/myAnalases)檢測位點之間的連鎖不平衡強(qiáng)度，采用SPSS 17.0軟件中的單因素方差分析和獨立樣本T檢驗分析分析德州驢基因型與生長性狀的相關(guān)性，結(jié)果用“平均數(shù)±標(biāo)準(zhǔn)誤”表示(P<0.05為差異顯著水平)。

2 結(jié)果與分析

2.1 德州驢SIRT3基因的多態(tài)性檢測

隨機(jī)選取德州驢的DNA樣本作為SIRT3-1和SIRT3-2位點的PCR擴(kuò)增模板，產(chǎn)物加入濃度為1%的瓊脂糖凝膠中進(jìn)行電泳鑒定。同時，對擴(kuò)增產(chǎn)物進(jìn)行測序和驗證。在SIRT3-1(g.19763T>C)發(fā)現(xiàn)了TT、TC、CC 3個基因型，在SIRT3-2(g.35148A>G)發(fā)現(xiàn)了AG、GG 2個基因型，測序結(jié)果如圖1。將404頭德州驢的DNA樣品送到康普生生物技術(shù)有限公司進(jìn)行質(zhì)譜分析，SIRT3-1和SIRT3-2位點分型結(jié)果如圖2所示。

2.2 群體遺傳學(xué)分析、連鎖不平衡分析

從群體遺傳學(xué)角度分別分析這2個SNPs的基因型頻率和等位基因頻率，并進(jìn)行卡方檢驗檢測Hardy-Weinberg平衡狀態(tài)，結(jié)果如表2、表3所示。SIRT3-1位點上，TC基因型頻率最高，PIC范圍約為0.30，說明母驢和公驢具有中度多態(tài)性(0.250.05)。

表2 德州驢SIRT3-1位點基因型、等位基因頻率及群體遺傳參數(shù)

表3 德州驢SIRT3-2位點基因型、等位基因頻率及群體遺傳參數(shù)

在連鎖不平衡分析中，連鎖強(qiáng)度的大小可以用D’和r2表示。r2<0.33表示SIRT3-1和SIRT3-2沒有連鎖關(guān)系，如圖3。

2.3 SIRT3基因型與生長性狀的關(guān)聯(lián)分析

對德州驢214頭公驢和190頭母驢的9個生長性狀分別與SIRT3基因2個SNPs位點不同基因型進(jìn)行關(guān)聯(lián)分析，結(jié)果見表4、表5。在SIRT3-1位點，公驢和母驢的TT型個體在體高和尻長方面均顯著性高于TC型個體均值(P<0.05)，即公母驢存在一致的顯著性差異。此外，公驢中TT型個體均值在尻高方面也顯著高于TC型個體均值(P<0.05)，而CC型個體均值在管圍方面顯著高于TT型個體均值(P<0.05)(如表4)。在SIRT3-2位點，公驢和母驢GG基因型個體均值在胸寬方面顯著高于AG型個體均值(P<0.05)，公驢GG型個體均值在尻寬方面極顯著高于AG型個體均值(P<0.01)(如表5)。

表4 德州驢SIRT3-1位點不同基因型與生長性狀的關(guān)聯(lián)分析

表5 德州驢SIRT3-2位點多態(tài)與生長性狀的關(guān)聯(lián)分析

3 討 論

Shi等人[15]研究發(fā)現(xiàn)，SIRT3可被寒冷刺激誘導(dǎo)并在棕色脂肪組織中表達(dá)冷暴露會增加小鼠棕色脂肪SIRT3 mRNA水平，靜息狀態(tài)下骨骼肌70%的能量來源于脂肪的分解代謝，由此可推測骨骼肌中SIRT3的表達(dá)對脂肪代謝具有重要影響。

研究[16-24]表明，SIRT3基因的突變位點與大圍子豬和湖村黑豬的肉質(zhì)性狀存在顯著性相關(guān)，表現(xiàn)為在SIRT3基因第5外顯子上存在的SNP突變位點顯著影響大圍子豬的24 h滴水損失、失水率和肌內(nèi)脂肪含量(P<0.05)。Gui等人[5]對SIRT1、SIRT2和SIRT3 3個基因的4個SNP位點與牛的體尺性狀進(jìn)行關(guān)聯(lián)分析，發(fā)現(xiàn)僅有SIRT1的基因變異及其相應(yīng)的基因型可以作為牛育種經(jīng)濟(jì)性狀的分子標(biāo)記。鄧冠群等[10]人對秦川牛SIRT3基因的研究表明，不同位點的不同基因型與體斜長、體高、腰角寬、尻長、胸深、胸圍、背膘厚、眼肌面積和肌內(nèi)脂肪含量顯著(P<0.05)或極顯著相關(guān)(P<0.01)。

本試驗通過質(zhì)譜分析對190頭母驢、214頭公驢共404頭德州驢的SIRT3基因進(jìn)行了多態(tài)性分析，發(fā)現(xiàn)了2個多態(tài)位點。通過基因分型在SIRT3-1發(fā)現(xiàn)3種基因型，而SIRT3-2只有兩種基因型，這可能與試驗中選擇的群體過小或者在進(jìn)化過程中自然和人為定向選擇有關(guān)。相對于外顯子區(qū)域的突變，越來越多的研究表明，內(nèi)含子區(qū)域的突變也能對性狀產(chǎn)生顯著的影響。雖然內(nèi)含子的突變不會對編碼的氨基酸序列發(fā)生影響，但其可能對順式作用元件的功能產(chǎn)生影響或者改變細(xì)胞內(nèi)mRNA的可變剪切，進(jìn)而對基因的功能或者表達(dá)產(chǎn)生影響。

4 結(jié) 論

本研究發(fā)現(xiàn)SIRT3基因在德州驢中存在多態(tài)且與生長性狀顯著關(guān)聯(lián)，德州驢群體中SIRT3基因的2個SNPs位點與體高、尻高、尻長、胸寬這4個生長性狀有顯著的相關(guān)性，說明SIRT3基因多態(tài)有潛力作為影響生長性狀的候選分子標(biāo)記用于輔助選擇，這為德州驢的選育工作提供了科學(xué)依據(jù)。