牦牛TMEM-18基因多態(tài)性及其與生產(chǎn)性能關(guān)聯(lián)分析

張歡張全偉王琪馬友記張勇趙興緒，

（1. 甘肅農(nóng)業(yè)大學(xué)動(dòng)物醫(yī)學(xué)院，蘭州 730070；2. 甘肅農(nóng)業(yè)大學(xué)生命科學(xué)技術(shù)學(xué)院，蘭州 730070；3. 甘肅農(nóng)業(yè)大學(xué)動(dòng)物科技學(xué)院，蘭州 730070）

牦牛TMEM-18基因多態(tài)性及其與生產(chǎn)性能關(guān)聯(lián)分析

張歡1張全偉2王琪1馬友記3張勇1趙興緒1，2

（1. 甘肅農(nóng)業(yè)大學(xué)動(dòng)物醫(yī)學(xué)院，蘭州 730070；2. 甘肅農(nóng)業(yè)大學(xué)生命科學(xué)技術(shù)學(xué)院，蘭州 730070；3. 甘肅農(nóng)業(yè)大學(xué)動(dòng)物科技學(xué)院，蘭州 730070）

跨膜蛋白18基因（TMEM-18）是一個(gè)末端為低聚嘧啶的基因，與動(dòng)物生長(zhǎng)發(fā)育密切相關(guān)。為了研究牦牛TMEM-18 基因多態(tài)性與生產(chǎn)性能的關(guān)系。以192頭天祝雌性牦牛血樣NDA構(gòu)建混合池，擴(kuò)增TMEM-18基因內(nèi)含子1和外顯子5的序列。運(yùn)用BLAST和Chromas軟件分析突變位點(diǎn)，運(yùn)用高分辨率熔解曲線分析技術(shù)（HRM）統(tǒng)計(jì)基因型分型，采用SHEsis 軟件進(jìn)行基因型頻率和等位基因頻率計(jì)算，同時(shí)對(duì)多態(tài)位點(diǎn)進(jìn)行配對(duì)連鎖不平衡和單倍型分析，運(yùn)用SPSS21.0分析基因多態(tài)位點(diǎn)與生產(chǎn)性能關(guān)聯(lián)。結(jié)果表明，牦牛TMEM-18基因內(nèi)含子1處存在2個(gè)多態(tài)位點(diǎn)，分別是861（A/G）和1267（C/T），外顯子5處存在1個(gè)多態(tài)位點(diǎn)為4 447（C/T）。關(guān)聯(lián)分析表明，牦牛TMEM-18基因861（A/G）位點(diǎn)的不同基因型與牦牛的體高、體重、胴體重和屠宰率差異性顯著（P＜0.05）；1267（C/T）位點(diǎn)的不同基因型與牦牛的體高、體重和屠宰率差異性顯著（P＜0.05），而牦牛TMEM-18基因4447（C/T）位點(diǎn)的不同基因型與牦牛的體斜長(zhǎng)、管圍、胸圍、體重、胴體重、凈肉重、凈肉率和屠宰率均存在顯著性差異（P＜0.05）。通過單倍型分析發(fā)現(xiàn)群體中存在3種單倍型組合，即ACC單倍型、GTC單倍型和GTT單倍型，其中ACC單倍型組合個(gè)體數(shù)目明顯多于其他單倍型組合，為優(yōu)勢(shì)單倍型。本實(shí)驗(yàn)揭示TMEM-18基因可作為牦牛生產(chǎn)性能開發(fā)的候選分子標(biāo)記，為牦牛遺傳資源的利用、開發(fā)與新品種的選育提供依據(jù)。

牦牛；TMEM-18基因；SNP；HRM；生產(chǎn)性狀

牦牛（Bos grunniens）是長(zhǎng)期放牧于海拔高達(dá)3 000 m以上的高寒地區(qū)特殊牛種，屬于地方類半野生型的種群，是經(jīng)過長(zhǎng)期的自然選擇和人工選育而形成的高原特有畜種資源。牦牛既是我國特有的畜種資源，也是我國稀有而珍貴的遺傳資源［1］。牦牛生長(zhǎng)發(fā)育緩慢、群體規(guī)模小、生產(chǎn)性能低下，從而影響牦牛的經(jīng)濟(jì)價(jià)值和使用價(jià)值［2］。

近年來，人們已對(duì)與牦牛生長(zhǎng)發(fā)育等相關(guān)基因進(jìn)行了分子結(jié)構(gòu)特征、多態(tài)性和表達(dá)分析，且有部分侯選基因的研究相對(duì)比較深入。如肌肉生長(zhǎng)抑制素（MSTN）基因和胰島素樣生長(zhǎng)因子家族（IGF-IR）基因的單核苷酸多態(tài)性（single nucleotide polymorphisms，SNPs）對(duì)不同年齡階段的牦牛生產(chǎn)性能有顯著影響［3］；胰島素樣生長(zhǎng)因子2（IGF-2）基因的SNPs對(duì)牦牛的生長(zhǎng)發(fā)育具有顯著差異［4］；生長(zhǎng)分化因子10（GDF-10）、印度刺猬因子（IHH）、黑素皮質(zhì)素受體-4（MC4R）等基因的SNPs對(duì)不同牦牛種群的生產(chǎn)性狀有顯著影響［5］等。上述基因均可作為牦牛選育的候選基因。

研 究 發(fā) 現(xiàn)，跨膜 蛋 白18（Transmembrane Proteins-18，TMEM-18）是一個(gè)末端為低聚嘧啶的基因，該基因能夠在體內(nèi)和體外參與神經(jīng)干細(xì)胞遷移的調(diào)節(jié)［6］。Hindorff等［7］通過全基因組關(guān)聯(lián)分 析（Genome-wide association study，GWAS） 發(fā)現(xiàn)TMEM-18基因與肥胖密切相關(guān)。此后研究表明TMEM-18基因主要通過大腦中的神經(jīng)系統(tǒng)來調(diào)節(jié)脂肪細(xì)胞分化從而調(diào)控肥胖［8］；而TMEM-18基因除神經(jīng)系統(tǒng)中表達(dá)外，在不同組織中均見表達(dá)［9］。近年來，人類全基因組的研究發(fā)現(xiàn)，TMEM-18基因存在與人肥胖連鎖的單核苷酸多態(tài)性（single nucleotide polymorphism，SNP），該基因的多個(gè)SNPs位點(diǎn)（rs6548238 、rs2867125、rs4854344和rs756131）均可增加人體質(zhì)指數(shù)（Body Mass Index，BMI）和體重［10-15］。早期研究發(fā)現(xiàn)TMEM-18基因可調(diào)節(jié)小鼠的生長(zhǎng)發(fā)育，不僅在小鼠大腦中強(qiáng)烈表達(dá)，而且在小鼠其他組織中也有表達(dá)［16］。在TMEM-18基因的家畜動(dòng)物研究中，大白豬、南陽牛、郟縣紅牛、草原紅牛、魯西牛和秦川牛中均檢測(cè)發(fā)現(xiàn)TMEM-18基因的多個(gè)SNPs位點(diǎn)均對(duì)家畜動(dòng)物的生長(zhǎng)性能有顯著影響［17，18］。因此，我們推測(cè)TMEM-18可能與牦牛生長(zhǎng)發(fā)育及其生產(chǎn)性能有密切關(guān)系。然而，牦牛TMEM-18基因的結(jié)構(gòu)、功能及SNPs與生產(chǎn)性能關(guān)聯(lián)分析的研究迄今尚未見報(bào)道。本實(shí)驗(yàn)擬對(duì)牦牛TMEM-18基因的SNPs進(jìn)行檢測(cè)，初步篩選與牦牛生產(chǎn)性狀相關(guān)的多態(tài)位點(diǎn)，分析其遺傳變異與生長(zhǎng)性能的關(guān)系，從而為今后牦牛遺傳改良和選育以及TMEM-18基因的功能研究奠定基礎(chǔ)。

1 材料與方法

1.1 材料

實(shí)驗(yàn)樣本采自于甘肅省天祝縣某屠宰場(chǎng)，選取2-3歲的牦牛群體采集192頭天祝雌性牦牛血樣。頸靜脈采血3 mL，輕微振蕩混勻后于-80℃保存?zhèn)溆谩?shí)地測(cè)量并記錄其生長(zhǎng)性狀（體斜長(zhǎng)、體高、胸圍、管圍、體重）及屠宰性能（胴體、凈肉、屠宰率、凈肉率）。

1.2 方法

1.2.1 基因組DNA提取 利用血液基因組DNA提取試劑盒（Omega Bio-Tek，D6919-01B）提取基因組DNA，經(jīng)1%的瓊脂糖凝膠電泳檢測(cè)DNA的完整性后使用紫外分光光度計(jì)檢測(cè)DNA的濃度。從上述檢測(cè)合格的DNA樣本中隨機(jī)抽取20 μL稀釋至18-20 ng/μL，用作PCR擴(kuò)增的模板，其余的母液于-80℃保存?zhèn)溆谩?/p>

1.2.2 基因組DNA池構(gòu)建 選取稀釋后的200份DNA樣本，每個(gè)樣本各取1 μL，混合均勻，構(gòu)建牦?；蚪MDNA混合池。

1.2.3 引物設(shè)計(jì)和PCR擴(kuò)增 參考NCBI數(shù)據(jù)庫中牦牛的TMEM-18基因（登錄號(hào)：NW_005393077.1），采用交叉重疊原理［15］運(yùn)用Primer 5.0軟件設(shè)計(jì)3對(duì)引物對(duì)部分內(nèi)含子1和外顯子5區(qū)域進(jìn)行擴(kuò)增（表1），引物由上海生工生物科技有限公司合成。

PCR擴(kuò)增體系（25 μL）：20 ng/μL的基因組DNA 1 μL，10 pmol 的上下游引物各1 μL，2×Taq PCR Master Mix（ 內(nèi) 含Taq DNA 聚 合 酶、Mg2+、dNTPs 等）（北京，百泰克）12.5 μL，滅菌超純水9.5 μL。

PCR反應(yīng)程序：95℃預(yù)變性3 min；94℃變性30 s，復(fù)性30 s（退火溫度見表1），72℃延伸30 s，35個(gè)循環(huán)；72℃延伸10 min；4℃保存。

表1 鑒定牦牛TMEM-18基因突變位點(diǎn)所用引物信息表

1.2.4 DNA測(cè)序分析 PCR擴(kuò)增產(chǎn)物經(jīng)瓊脂糖凝膠回收試劑盒（北京，全式金）回收純化后送至蘇州金唯智生物科技有限公司進(jìn)行測(cè)序，測(cè)序結(jié)果通過DNAMAN軟件和Chromas軟件比對(duì)分析。

1.2.5 HRM小片段法基因分型 針對(duì)所發(fā)現(xiàn)的多態(tài)位點(diǎn)，利用Light Scanner Primer Design Software（Idaho公司，美國）設(shè)計(jì)出用于高分辨率熔解曲線（High Resolution melting，HRM）分析多態(tài)位點(diǎn)的引物（表2），對(duì)試驗(yàn)群體樣本進(jìn)行基因型分析［19，20］。合成2對(duì)高低溫內(nèi)標(biāo)［21］（表 3）對(duì)溶解曲線進(jìn)行校正，提高分型的精確度。

表2 HRM 分析多態(tài)位點(diǎn)的引物信息表

表3 高低溫內(nèi)標(biāo)序列

HRM-PCR擴(kuò)增體系（11 μL）：20 ng/μL的基因組DNA 1 μL，10 pmol 的上下游引物各0.2 μL，2×Taq PCR Master Mix（ 內(nèi) 含Taq DNA 聚 合 酶、Mg2+、dNTPs 等）（北京，百泰克）5 μL，滅菌超純水3.6 μL，10×LC Green 飽和染料（美國，Idaho）1 μL。

HRM-PCR反應(yīng)程序：95℃預(yù)變性5 min；94℃變性20 s，Tm（℃）（表2）復(fù)性20 s，72℃延伸20 s，35個(gè)循環(huán)；72℃延伸10 min。

高低溫內(nèi)標(biāo)稀釋方法：1 OD內(nèi)標(biāo)單鏈加 400 μL水，退火體系：飽和氯化鈉1 μL，內(nèi)標(biāo)互補(bǔ)雙鏈各1 μL，補(bǔ)水至10 μL，95℃水浴3 min，緩慢降至室溫。

HRM熒光信號(hào)采集：將稀釋好的高低溫內(nèi)標(biāo)各1 μL加入HRM-PCR產(chǎn)物單孔中，95℃水浴30 s；25℃水浴30 s；放入LightScanner 96（美國，Idaho）中收集熒光信號(hào)。

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

根據(jù)群體遺傳學(xué)理論［22］計(jì)算基因型和等位基因頻率，統(tǒng)計(jì)牦牛群體TMEM-18基因多態(tài)位點(diǎn)的純和度Ho、雜合度He、有效等位基因數(shù)Ne、多態(tài)信息含量PIC［23-25］，并進(jìn)行Hardy-Weinberg平衡檢驗(yàn)。采用SHEsis軟件進(jìn)行配對(duì)連鎖不平衡分析。采用SPSS 21.0軟件對(duì)TMEM-18基因多態(tài)性位點(diǎn)在牦牛群體中的基因型分布進(jìn)行卡方檢測(cè)，并對(duì)不同基因型與生產(chǎn)性狀間的關(guān)聯(lián)性進(jìn)行分析。考慮到性別、環(huán)境和年齡因素的影響，采用固定模型分析基因型效應(yīng)對(duì)體尺性狀的影響。

式中，Yijk：個(gè)體表型記錄；u：群體均值；Gi：標(biāo)記基因型效應(yīng)；Sj：性別因素；Ak：年齡效應(yīng)；Eijk：隨機(jī)誤差。

2 結(jié)果

2.1 PCR擴(kuò)增產(chǎn)物

以牦牛血液基因組DNA為模板，根據(jù)設(shè)計(jì)的引物進(jìn)行PCR特異性擴(kuò)增，結(jié)果見圖1。由圖1可知PCR擴(kuò)增產(chǎn)物均為單一特異性條帶。產(chǎn)物經(jīng)過測(cè)序后長(zhǎng)度分別為1 001、953和722 bp，與引物設(shè)計(jì)時(shí)預(yù)期長(zhǎng)度一致。

圖1 TMEM-18基因PCR擴(kuò)增電泳圖

2.2 牦牛TMEM-18基因HRM方法的分型

對(duì)PCR擴(kuò)增產(chǎn)物進(jìn)行測(cè)序，結(jié)果（圖 2）顯示，牦牛TMEM-18基因內(nèi)含子1存在2個(gè)多態(tài)位點(diǎn)，分別是861（A/G）和1267（C/T），外顯子5存在1個(gè)多態(tài)位點(diǎn) 4447（C/T）。通過DNA混合池對(duì)牦牛TMEM-18基因測(cè)序后，應(yīng)用高分辨率熔解曲線法對(duì)3個(gè)多態(tài)位點(diǎn)進(jìn)行基因型分析（圖 3）。HRM 結(jié)果（圖 3-A）顯示，牦牛TMEM-18基因861（A/G）位點(diǎn)處存在3種基因型，分別是AA、AG和GG。統(tǒng)計(jì)分析發(fā)現(xiàn)，AA基因型屬于優(yōu)勢(shì)基因型（n=130），個(gè)體數(shù)明顯多于AG（n=56）和GG（n=6）型；牦牛TMEM-18基因1267（C/T）和 4447（C/T）位點(diǎn)處均存在3種基因型，分別為CC、CT和TT 型（圖3-B和圖 3-C）。統(tǒng)計(jì)分析發(fā)現(xiàn)，CC基因型均屬于優(yōu)勢(shì)基因型。

圖2 TMEM-18基因3個(gè)多態(tài)位點(diǎn)的測(cè)序峰圖

2.3 牦牛TMEM-18基因多態(tài)位點(diǎn)遺傳多態(tài)性的分析

對(duì)牦牛TMEM-18基因多態(tài)位點(diǎn)進(jìn)行基因型頻率和基因頻率以及相關(guān)遺傳特性指標(biāo)計(jì)算，結(jié)果（表4）顯示，牦牛TMEM-18基因 861（A/G）的優(yōu)勢(shì)等位基因?yàn)锳，其基因頻率為0.823；牦牛TMEM-18基因 1267（C/T）和 4447（C/T）的優(yōu)勢(shì)等位基因均為C，其基因頻率均大于0.8 以上。3個(gè)位點(diǎn)均表現(xiàn)為低度多態(tài)（PIC＜0.25），經(jīng)Hardy-Weinberg平衡檢驗(yàn)發(fā)現(xiàn)，3個(gè)位點(diǎn)均未達(dá)到Hardy-Weinberg平衡（P＜0.05）（表 5）。配對(duì)連鎖不平衡分析結(jié)果（表 6）表明，牦牛TMEM-18基因中上述位點(diǎn)間存在強(qiáng)連鎖不平衡（D'＞0.75，r2＞0.33）。多態(tài)位點(diǎn)的單倍型組合的分布及頻率結(jié)果（表 7）顯示，192頭牦牛群體中存在3種單倍型，其中ACC在研究群體中個(gè)體數(shù)最多（＞80%），而GTC和GTT單倍型在群體中只存在少量的個(gè)體。

圖3 TMEM-18基因3個(gè)多態(tài)性位點(diǎn)HRM分型曲線

表4 TMEM-18基因突變位點(diǎn)的基因型頻率及基因頻率

表5 TMEM-18基因不同突變位點(diǎn)的遺傳多態(tài)性

表6 TMEM-18基因多態(tài)位點(diǎn)配對(duì)連鎖不平衡分析

表7 TMEM-18基因多態(tài)位點(diǎn)的單倍型組合的分布及頻率

2.4 TMEM-18基因基因型及單倍型與生產(chǎn)性狀的關(guān)聯(lián)性分析

牦牛TMEM-18基因861（A/G）和1267（C/T）位點(diǎn)的不同基因型與不同個(gè)體的體高均存在顯著差異（P＜0.05），而牦牛體斜長(zhǎng)、管圍和胸圍差異均不顯著（P＞0.05），其中861（A/G）位點(diǎn)的GG 基因型的個(gè)體的生產(chǎn)性能均高于AA和AG基因型的個(gè)體；和1267（C/T）位點(diǎn)的TT 基因型的個(gè)體的生產(chǎn)性能均高于CC、CT基因型個(gè)體；牦牛TMEM-18基因4447（C/T）位點(diǎn)的不同基因型個(gè)體在體斜長(zhǎng)、管圍和胸圍均存在差異顯著（P＜0.05），而體高差異不顯著（P＞0.05），其中CT基因型的個(gè)體的生產(chǎn)性能均高于CC、TT基因型個(gè)體（表8）。

牦牛TMEM-18基因861（A/G）位點(diǎn)的不同基因型個(gè)體的體重、胴體重和屠宰率均存在顯著差異（P＜0.05），而凈肉重和凈肉率差異均不顯著（P＞0.05），其中GG基因型的個(gè)體的屠宰性能均高于AA和AG基因型的個(gè)體；1 267（C/T）位點(diǎn)的不同基因型個(gè)體的體重、屠宰率存在顯著差異（P＜0.05），而胴體重和凈肉重和凈肉率差異均不顯著（P＞0.05），其中TT基因型的個(gè)體的屠宰性能均高于CC和CT基因型的個(gè)體；4 447（C/T）位點(diǎn)的不同基因型個(gè)體的體重、胴體重、凈肉重、凈肉率和屠宰率均存在顯著差異（P＜0.05），其中CT基因型的個(gè)體的屠宰性能均高于CC和TT基因型的個(gè)體（表9）。

表8 牦牛TMEM-18基因多態(tài)位點(diǎn)不同基因型與生產(chǎn)性狀關(guān)聯(lián)分析

表9 牦牛TMEM-18基因多態(tài)位點(diǎn)不同基因型與屠宰性能關(guān)聯(lián)分析

3 討論

本實(shí)驗(yàn)通過DNA混合池測(cè)序得到牦牛TMEM-18基 因3個(gè)SNPs（861（A/G）、1 267（C/T） 和4 447（C/T）），其中1 267（C/T）和4 447（C/T）兩個(gè)位點(diǎn)為TMEM-18基因上新發(fā)現(xiàn)的SNPs，其他物種均未見報(bào)道。牦牛TMEM-18基因861（A/G）位點(diǎn)存在3種基因型，分別是AA、GG和AG 型，且等位基因A是優(yōu)勢(shì)等位基因；1 267（C/T）和4 447（C/T）位點(diǎn)均存在3種基因型，分別是CC、TT和 CT 型，且等位基因C是優(yōu)勢(shì)等位基因。這3個(gè)多態(tài)位點(diǎn)的優(yōu)勢(shì)等位基因頻率均達(dá)到了0.8以上，這可能與樣本的種群和自然環(huán)境有關(guān)。此外，這3個(gè)多態(tài)位點(diǎn)均為低度多態(tài)，遺傳變異低，選擇余地較小，這可能因牦牛長(zhǎng)期生長(zhǎng)在高海拔環(huán)境下，長(zhǎng)年的高寒、缺氧等惡劣的自然環(huán)境下使有效等位基因和多態(tài)信息含量降低。而連鎖不平衡是單倍型基因非隨機(jī)分布的現(xiàn)象。若兩個(gè)連鎖基因座的兩個(gè)等位基因頻率分別為p1和p2，這兩個(gè)等位基因組成的單體型的頻率h如果為預(yù)期值（plxp2），則稱這兩個(gè)基因座是連鎖平衡，如果h高于plxp2就存在連鎖不平衡［26］。本實(shí)驗(yàn)中，Hardy-Weinberg平衡檢驗(yàn)發(fā)現(xiàn)，3個(gè)位點(diǎn)均未達(dá)到Hardy-Weinberg平衡（P＜0.05），這可能與本次實(shí)驗(yàn)的數(shù)據(jù)結(jié)構(gòu)有關(guān)，還需要進(jìn)一步擴(kuò)大樣品進(jìn)行研究，同時(shí)也表明牦牛品種的選育、遷移和遺傳漂變等因素對(duì)TMEM-18基因的突變?cè)斐捎绊?，可能?dǎo)致未達(dá)到平衡。近年來在中國牛上也開始了對(duì)TMEM-18 基因的序列測(cè)序和多態(tài)性研究，馬偉［17］在不同品種的牛TMEM-18基因中共發(fā)現(xiàn)4個(gè)不同SNPs（1085（T/G）、2303（A/G）、3835（G/A）和3865（G/A）），均未達(dá)到Hardy-Weinberg平衡（P＜0.05），與本結(jié)果相似。

牦牛TMEM-18基因內(nèi)含子1處存在2個(gè)多態(tài)位點(diǎn)，分別是861（A/G）和1267（C/T）。內(nèi)含子雖為基因中無編碼功能的序列，但其對(duì)基因表達(dá)調(diào)控有著決定作用［27］。一些基因中的內(nèi)含子對(duì)基因表達(dá)具有正調(diào)控作用，Chang等［28］發(fā)現(xiàn)豬MyHC基因內(nèi)含子中存在重要調(diào)控元件，可以調(diào)控轉(zhuǎn)錄起始來增強(qiáng)基因表達(dá)。而張冉［29］檢測(cè)魯西黃牛TMEM-18基因在起始密碼子上游541 bp和861 bp處存在SNPs位點(diǎn)，且證實(shí)861 bp處的SNPs位點(diǎn)影響啟動(dòng)子活性，表現(xiàn)為負(fù)調(diào)控作用。由此推測(cè)，牦牛TMEM-18基因內(nèi)含子861（A/G）突變也可能參與介導(dǎo)了該基因的轉(zhuǎn)錄調(diào)控。因此，本研究對(duì)了解與分析該基因遺傳特征及其對(duì)重要經(jīng)濟(jì)性狀的影響仍然有重要的意義。

關(guān)聯(lián)分析表明，牦牛TMEM-18基因861（A/G）位點(diǎn)的不同基因型與牦牛的體高、體重、胴體重和屠宰率差異性顯著（P＜0.05）；1267（C/T）位點(diǎn)的不同基因型與牦牛的體高、體重和屠宰率差異性顯著（P＜0.05），而牦牛TMEM-18基因4447（C/T）位點(diǎn)的不同基因型與牦牛的體斜長(zhǎng)、管圍、胸圍、體重、胴體重、凈肉重、凈肉率和屠宰率均存在顯著性差異（P＜0.05）。此外，本實(shí)驗(yàn)在牦牛后代群體中找到了3種單倍型組合，其中ACC單倍型組合個(gè)體數(shù)目高達(dá)80%以上，表明ACC單倍型為主要單倍型。這與目前相關(guān)文獻(xiàn)對(duì)動(dòng)物TMEM-18基因多態(tài)性及其與生產(chǎn)性能關(guān)聯(lián)分析研究報(bào)道的相關(guān)結(jié)論一致。例如，林平［16］在大白豬TMEM-18基因中共發(fā)現(xiàn)了1個(gè)SNPs（91（T/G）），且與豬生長(zhǎng)性狀和腰薦結(jié)合處背膘厚性狀顯著關(guān)聯(lián)（P＜0. 01），說明TMEM-18基因具有在選種中運(yùn)用的價(jià)值。馬偉［17］在不同品種的牛TMEM-18基因中共發(fā)現(xiàn)4個(gè)不同SNPs，對(duì)中國牛的體尺性狀有顯著影響，并說明在其4個(gè)新SNPs中，有3個(gè)可作為中國牛品種體尺和體重等生長(zhǎng)性狀選擇分子育種的候選標(biāo)記。因此本研究推測(cè)TMEM-18基因上的3個(gè)多態(tài)位點(diǎn)可影響牦牛的生長(zhǎng)性狀。綜上所述，牦牛TMEM-18基因可作為牦牛生產(chǎn)性狀有關(guān)的候選分子標(biāo)記，對(duì)提高分子標(biāo)記輔助選擇的準(zhǔn)確性，為今后牦牛分子育種提供主要的理論依據(jù)，為牦牛遺傳育種提供參考。

4 結(jié)論

牦牛TMEM-18基因3個(gè)多態(tài)位點(diǎn)與牦牛的體尺指標(biāo)和屠宰性能顯著相關(guān)，可以嘗試作為牦牛生長(zhǎng)性狀的候選分子標(biāo)記，應(yīng)用于牦牛分子標(biāo)記輔助育種中。

［1］何梅蘭, 馬秀山. 天祝白牦牛產(chǎn)業(yè)發(fā)展現(xiàn)狀及發(fā)展戰(zhàn)略［J］.中國草食動(dòng)物科學(xué), 2015, 35（5）：71-74.

［2］馬蒼. 天祝白牦牛產(chǎn)業(yè)可持續(xù)發(fā)展思考［J］. 家畜生態(tài)學(xué)報(bào), 2008, 29（6）：133-135.

［3］梁春年. 牦牛MSTN和IGF-IR基因的克隆及SNPs與生長(zhǎng)性狀相關(guān)性［D］. 蘭州：甘肅農(nóng)業(yè)大學(xué), 2011.

［4］李亞蘭. 天祝白牦牛IGF-2基因的全長(zhǎng)cDNA克隆, 組織表達(dá)及生產(chǎn)性狀相關(guān)性分析［D］. 蘭州：甘肅農(nóng)業(yè)大學(xué), 2015.

［5］李天科. 牦牛GDF-10, IHH, MC4R基因的多態(tài)與生長(zhǎng)性狀的相關(guān)性研究［D］. 蘭州：甘肅農(nóng)業(yè)大學(xué), 2014.

［6］Jurvansuu J, Zhao Y, Leung DS, et al. Transmembrane protein 18 enhances the tropism of neural stem cells for glioma cells［J］. Cancer Research, 2008, 68（12）：4614-4622.

［7］Hindorff LA, Sethupathy P, Junkins HA, et al. Potential etiologic and functional implications of genome-wide association loci for human diseases and traits［J］. Proc Natl Acad Sci USA, 2009, 106（23）：9362-9367.

［8］Gutierrez-Aguilar R, Kim DH, Woods SC, et al. Expression of new loci associated with obesity in diet-induced obese rats：from genetics to physiology［J］. Obesity（Silver Spring）, 2011, 20, 306-312.

［9］Bernhard F, Landgraf K, Kl?ting N, et al. Functional relevance of genes implicated by obesity genome wide association study signals for human adipocyte biology［J］. Diabetologia, 2013, 56（2）：311-322.

［10］Rask-Andersen M, Jacobsson JA, Moschonis G, et al. Association of TMEM18 variants with BMI and waist circumference in children and correlation of mRNA expression in the PFC with body weight in rats［J］. European Journal of Human Genetics, 2012, 20（2）：192-197.

［11］Hotta K, Nakamura M, Nakamura T, et al. Association between obesity and polymorphisms in SEC16B, TMEM18, GNPDA2, BDNF, FAIM2 and MC4R in a Japanese population［J］. Journal of Human Genetics, 2009, 54（12）：727-731.

［12］Holzapfel C, Grallert H, Huth C, et al. Genes and lifestyle factors in obesity：Results from 12462 subjects from MONICA/KORA［J］. International Journal of Obesity, 2010, 34（10）：1538-1545.

［13］Hong KW, Oh B. Recapitulation of genome-wide association studies on body mass index in the Korean population［J］. International Journal of Obesity, 2011, 36（8）：1127-1130.

［14］Zhao J, Bradfield JP, Li M, et al. The role of obesity-associated loci identified in genome-wide association studies in the determination of pediatric BMI［J］. Obesity（Silver Spring）, 2009, 17（12）：2254-2257.

［15］Almén MS, Jacobsson JA, Shaik JH, et al. The obesity gene, TMEM18, is of ancient origin, found in majority of neuronal cells in all major brain regions and associated with obesity in severely obese children［J］. BMC Medical Genetics, 2010, 11（1）：58.

［16］林平. 大白豬和通城豬肌纖維組織學(xué)特性差異和SH2B1, TEM18, IDH3B基因遺傳效應(yīng)分析［D］. 武漢：華中農(nóng)業(yè)大學(xué), 2013.

［17］馬偉. 牛TMEM18基因克隆, SNP檢測(cè)及其與部分經(jīng)濟(jì)性狀的關(guān)聯(lián)分析［D］. 西安：西北農(nóng)林科技大學(xué), 2012.

［18］郭憲, 李天科, 裴杰, 閻萍等. 南德溫雜交肉牛生長(zhǎng)分化因子10基因多態(tài)性與體尺性狀的關(guān)聯(lián)分析［J］. 中國畜牧獸醫(yī), 2014, 41（12）：197-203.

［19］Michael L, Robert P, Robert P, et al. Genotyping of singlenucleotide polymorphisms by high- resolution melting of small amplicons［J］. Molecular Diagnostics and Genetics, 2004, 50（7）：1156-1164.

［20］Ye MH, Chen JL, Zhao GP, et al. Associations of A-FABP and H-FABP markers with the content of intramuscular-lar fat in Beijing-you chicken［J］. Animal Biotechnology, 2010, 21：14-24.

［21］曾凱. 群體遺傳學(xué)理論在檢測(cè)正向選擇與研究物種形成機(jī)制中的應(yīng)用［D］. 廣州：中山大學(xué), 2007.

［22］Cameron N, Steven F. Base-pair neutral homozygotes can be discriminated by calibrated high-resolution melting of small amplicons［J］. Nucleic Acids Research, 2008, 1093（10）：1-8.

［23］Nei M, Roychoudhury AK. Sampling variances of heterqozygosity and genetic distance［J］. Genetics, 1974, 76（2）：379-390.

［24］Nei M, Li WH. Mathematical model for studying genetic variation in terms of restriction endonucleases［J］. Proc Natl Acad Sci USA, 1979, 76（10）：5269-5273.

［25］Yamashita R, Suzuki Y, Takeuchi N, et al. Comprehensive detection of human terminal oligo-pyrimidine（TOP）genes and analysis of their characteristics［J］. Nucleic Acids Research, 2008, 36（11）：3707-3715.

［26］牛曉亮. 牦牛CAPN4和CAST基因多態(tài)性及與胴體和肉質(zhì)性狀關(guān)聯(lián)研究［D］. 蘭州：甘肅農(nóng)業(yè)大學(xué), 2015.

［27］Caws E, White SN, Riley DG, et al. Assessment of single nucleotide polymorphisms in genes residing on chromosomes 14 and 29 for association with carcass composition traits in Bos indicus cattle［J］. Journal of Animal Science, 2005, 83（1）：13-19.

［28］Chang KC. Critical regulatory domains in intron 2 of porcine sarcomeric myosin heavy chain gene［J］. Journal of Muscle Research and Cell Motility, 2000, 21：451-461.

［29］張冉. 肉質(zhì)性狀相關(guān)基因的表達(dá)及消化道菌群多樣性與動(dòng)物生長(zhǎng)發(fā)育的關(guān)系［D］. 濟(jì)南：山東大學(xué), 2014.

（責(zé)任編輯 李楠）

Polymorphism of TMEM-18 Gene and Its Correlation with Production Traits in Yak

ZHANG Huan1ZHANG Quan-wei2WANG Qi1MA You-ji3ZHANG Yong1ZHAO Xing-xu1，2
（1. College of Veterinary Medicine，Gansu Agriculture University，Lanzhou 730070；2. College of Life Science and Technology，Gansu Agriculture University，Lanzhou 730070；3. College of Animal Science＆Technology，Gansu Agriculture University，Lanzhou 730070）

The TMEM-18（transmembrane protein 18）gene is the one with oligomer in its terminal，and closely related to the growth and development of animal. This study aims to study the correlation between polymorphism of TMEM-18 gene and production traits in yak. Firstly，constructing DNA mixing pool with blood samples from 192 female Tianzhu yaks，and the intron1 and exon5 of TMEM-18 gene were amplified. Then，BLAST and Chromas were used to analyze the mutation sites in TMEM-18 gene，high-resolution melting curve（HRM）to determine the genotype，software SHEsis to calculate and the frequencies of genotypes and alleles and to analyze the linkage disequilibrium and haplotypes in the SNPs of TMEM-18 gene，and software SPSS21.0 to understand the correlation between gene polymorphic loci and product traits. The experimental results showed that the intron1 of yak TMEM-18 gene contained 2 polymorphic loci of 861（A/G）and 1 267（C/T），and the exon5 contained 1 polymorphic locus of 4 447（C/T）. Correlation analysis indicated that the different genotypes at locus 861（A/G）resulted in the significant differences on body height，body weight，carcass weight，and slaughter rate（P＜0.05）；while the varied genotypes at locus 1 267（C/T）caused significant differences on body height，body weight and slaughter rate（P＜0.05）；and the various genotypes at locus 4 447（C/T）let to significant differences on body length，tube girth，chest girth，body weight，carcass weight，meat weight，netmeat rate，and slaughter rate（P＜0.05）. There were 3 haplotype combinations in the yak populations，namely ACC，GTC and GTT. Among them，the number of ACC was significantly more than that of others，i.e.，it was dominant haplotype. The results suggested that the TMEM-18 gene could be used as molecular-marker candidate for developing the production traits of yak，providing a scientific basis for the utilization and development of yak genetic resources and breeding of new varieties.

yak；TMEM-18 gene；SNP；HRM；production traits

10.13560/j.cnki.biotech.bull.1985.2017.02.013

2016-06-04

“十二五”農(nóng)村領(lǐng)域國家科技計(jì)劃課題（2013AA102503-3）

張歡，女，碩士研究生，研究方向：動(dòng)物醫(yī)學(xué)工程；E-mail：1373707490@qq.com

趙興緒，博士，教授，博士生導(dǎo)師；E-mail：zhaoxx@gsau.edu.cn