林鵬飛+任麗群+李平+顧麗菊+李謙+燕志宏

實驗室/貴州大學(xué)動物科學(xué)學(xué)院，貴州貴陽 550025； 3.貴州省種畜禽種質(zhì)測定中心，貴州貴陽 550018）摘要：為驗證BMPR-IB基因單核苷酸多態(tài)性（SNP）位點與宗地花豬繁殖性狀的關(guān)聯(lián)性，從而提高宗地花豬繁殖性能，提供有價值的分子輔助標記，推進宗地花豬的選育進程，以貴州宗地花豬經(jīng)產(chǎn)母豬為試驗材料，利用DNA池和DNA直接測序方法檢測BMPR-IB基因的多態(tài)性，對各突變位點不同基因型個體與總產(chǎn)仔數(shù)、產(chǎn)活仔數(shù)、仔豬平均初生質(zhì)量和母豬泌乳力性狀進行關(guān)聯(lián)性分析。在受試群體中發(fā)現(xiàn)，基因BMPR-IB有2個突變位點，G595C、C643T均位于第6外顯子，經(jīng)產(chǎn)母豬總產(chǎn)仔數(shù)、產(chǎn)活仔數(shù)各基因型均表現(xiàn)為AA>AG>GG、BB>BT>TT，產(chǎn)活仔數(shù)為7頭的母豬泌乳力表現(xiàn)為AA基因型個體與AG、GG基因型個體間差異極顯著（P<0.01）；各基因型仔豬平均初生質(zhì)量表現(xiàn)為 TT>BB>BT，但差異不顯著；母豬AA基因型總產(chǎn)仔數(shù)多于AG與GG基因型，AA基因型為優(yōu)勢基因型，A為優(yōu)勢等位基因；BB基因型在總產(chǎn)仔數(shù)、產(chǎn)活仔數(shù)的繁殖性狀中多于BT、TT基因型，BB為優(yōu)勢基因型，B為優(yōu)勢等位基因。推測等位基因A和等位基因B可能為宗地花豬總產(chǎn)仔數(shù)、產(chǎn)活仔數(shù)的有利等位基因，且AA基因型和BB基因型為有利基因型。

關(guān)鍵詞：宗地花豬；BMPR-IB基因；多態(tài)性；SNP位點；繁殖性狀；關(guān)聯(lián)性分析

中圖分類號： S828.3文獻標志碼： A文章編號：1002-1302（2017）17-0035-04

通信作者：燕志宏，教授，研究方向為動物遺傳育種與種質(zhì)資源創(chuàng)新。E-mail：yzh611127@sina.com。骨形態(tài)發(fā)生蛋白受體（bone morphogenetic proteins receptor，簡稱BMPR）是骨形態(tài)發(fā)生蛋白（BMPs）發(fā)揮生理功能的基礎(chǔ)，BMPR屬于TGF-β受體超家族中的一個亞家族成員。BMPR可以分為Ⅰ和Ⅱ2種亞型，其中Ⅰ型受體可以分為BMPR-IA（ALK-3）、BMPR-IB（ALK-6）、ActR-IA（ALK-2），Ⅱ型受體則分為BMPR-Ⅱ、ActR-ⅡA、ActR-ⅡB等3種[1]。骨形態(tài)發(fā)生蛋白受體IB基因?qū)Ψ敝承誀畹挠绊懽钤缡怯蒞ilson等于1980年發(fā)現(xiàn)的，BMPR-IB基因突變后會引起母羊Booroola Merino高產(chǎn)。后來，BMPR-IB基因被改名為FecB基因[2]。

研究表明，BMPR-IB基因是與母豬的生殖系統(tǒng)、卵泡的生長發(fā)育、胚胎發(fā)育等密不可分的基因。試驗證明BMPR-IB基因在豬的顆粒細胞中有表達[3]。宋一萍對8個豬品種的BMPR-IB基因進行研究發(fā)現(xiàn)，在該基因的編碼區(qū)第 369 bp 處C突變成T，山東地方豬種中在該位點的AA基因型總產(chǎn)仔數(shù)、產(chǎn)活仔數(shù)分別比BB基因型多0.45、0.51頭/胎，但差異不顯著；在引進豬種中BB基因型母豬的總產(chǎn)仔數(shù)、產(chǎn)活仔數(shù)分別比AA型多1.05、 0.90頭/胎（P<0.05）[4]。孟洪莉以長白豬、大白豬、杜洛克豬、槐豬為研究對象，發(fā)現(xiàn)BMPR-IB基因在位點A1081C處，長白豬、大白豬和杜洛克豬的經(jīng)產(chǎn)母豬的總產(chǎn)仔數(shù)和產(chǎn)活仔數(shù)均表現(xiàn)為AA>CC>AC，其中大白母豬AA基因型個體的總產(chǎn)仔數(shù)顯著高于AC基因型（P<0.05），AA基因型個體產(chǎn)活仔數(shù)極顯著高于AC基因型（P<0.01），杜洛克母豬AA基因型個體的產(chǎn)活仔數(shù)顯著高于CC基因型（P<0.05）；在該基因T19588A處，經(jīng)產(chǎn)和初產(chǎn)長白母豬的TA基因型個體總產(chǎn)仔數(shù)和產(chǎn)活仔數(shù)顯著高于TT基因型（P<0.05）[5]。

宗地花豬是貴州省優(yōu)良地方豬種之一，皮厚而糯、肉嫩味香，為人們所喜愛。但飼養(yǎng)周期長、產(chǎn)仔數(shù)少、生產(chǎn)效率低、成本高、養(yǎng)殖投入高、收益低等使宗地花豬養(yǎng)殖群體規(guī)模發(fā)展緩慢。針對宗地花豬的保種選育，相關(guān)機構(gòu)和學(xué)者做了很多研究，推進了宗地花豬的開發(fā)利用。馮文豪等對宗地花豬的胴體肉質(zhì)性狀進行測定，結(jié)果表明雖然宗地花豬瘦肉率明顯低于大白豬，但肉質(zhì)性狀優(yōu)于大白豬[6]。燕志宏等研究宗地花豬1世代的生長發(fā)育情況[7]，隨后在2011年研究了宗地花豬及其與長白豬、大白豬的雜種1代豬哺乳期和保育期的生長性能[8]。燕志宏等還對產(chǎn)區(qū)保種選育基礎(chǔ)群進行了繁殖性能的測定，隨后開始對宗地花豬繁殖性能的研究[9]。隨著分子生物學(xué)的蓬勃發(fā)展，燕志宏等分別對影響宗地花豬肉質(zhì)和繁殖性狀的一些相關(guān)基因進行了研究，取得了一定的研究成果[10-14]。

1材料與方法

1.1試驗動物樣本

66頭宗地花豬經(jīng)產(chǎn)母豬來自貴州紫薇畜牧業(yè)開發(fā)有限公司宗地花豬飼養(yǎng)場，各母豬均健康，在相同的環(huán)境和管理條件下進行飼養(yǎng)，用耳缺鉗采取每頭母豬耳根邊緣耳組織，放入裝有70%乙醇的1.5 mL離心管中，放入冰盒中帶回實驗室，于-20 ℃低溫冰箱中保存?zhèn)溆谩?/p>

1.2試劑與儀器

主要試劑和儀器：Ezup柱式動物基因組DNA提取試劑盒、核酸染料、DL 2000 DNA Marker、 2×Taq Master Mix、瓊脂糖等，均購自生工生物工程（上海）股份有限公司；去離子水、無水乙醇、TAE緩沖液等試劑為實驗室自備。PCR擴增儀（PTC200，美國Bio-Rad公司）、凝膠成像系統(tǒng)（BioSens SC 710，美國Bio-Rad公司）、微量移液器（ACURA825，SOCOREX公司）。

1.3DNA提取

用基因組DNA提取試劑盒提取DNA，用含有核酸染料的1%瓊脂糖凝膠電泳檢測純度，并用紫外分光光度計測量每個樣品DNA濃度，共測3次，取其平均值，-20 ℃保存?zhèn)溆谩?/p>

1.4引物設(shè)計及PCR擴增

參照美國國立生物技術(shù)信息中心（National Center for Biotechnology Information，簡稱NCBI）提供的豬BMPR-IB基因序列（登錄號為NC-010450.3），利用Premier 5.0軟件結(jié)合Oligo 6.0軟件設(shè)計引物，應(yīng)用NCBI Blast軟件檢測所設(shè)計引物的特異性（表1）。引物由英濰捷基（上海）貿(mào)易有限公司合成。1.5統(tǒng)計分析

利用DNAStar等軟件分析測序結(jié)果，篩選SNPs位點，分析不同位點在受試群體中的基因型頻率和遺傳特異性，并進行Hardy-Weinbery（哈迪-溫伯格）平衡檢測；根據(jù)最小二乘線性模型，分析不同基因型與樣本母豬繁殖性能[總產(chǎn)仔數(shù)、產(chǎn)活仔數(shù)、仔豬平均初生質(zhì)量、泌乳力（即平均泌乳量）]的關(guān)聯(lián)效應(yīng)。

2結(jié)果與分析

2.1PCR擴增

對66頭宗地花豬經(jīng)產(chǎn)母豬進行PCR擴增，經(jīng)1%瓊脂糖凝膠電泳檢測，均顯示為1條特異性條帶，其片斷長度分別為398、389 bp，與預(yù)期目標片段相符，且條帶清晰、整齊無拖尾（圖1），可用于測序。

2.2BMPR-IB基因DNA池PCR測序結(jié)果

將所得PCR產(chǎn)物送往北京諾賽基因組研究中心有限公司進行純化切膠回收測序，對所得序列運用SeqMan、BioXM等軟件進行分析比對，在BMPR-IB基因上發(fā)現(xiàn)2個SNPs位點，均位于第6外顯子上。測序峰如圖2所示。

2.3宗地花豬個體BMPR-IB基因測序結(jié)果分析

運用SeqMan軟件對所有樣品的測序結(jié)果分析發(fā)現(xiàn)，第6外顯子G595C位點處發(fā)生A→G堿基突變，由圖3可知，有3種基因型，分別命名為AA、AG、GG基因型。第6外顯子C643T位點處發(fā)生B→T堿基突變，發(fā)現(xiàn)3種基因型，分別命名為BB、BT、TT基因型，有2個等位基因B和T。通過序列比對發(fā)現(xiàn)，AA、BB基因型個體序列與GenBank上所提交的豬的基因序列一致，定義為野生型，GG、TT型則定義為突變型，AG、BT基因型定義為雜合型。

2.4基因型頻率及遺傳特異性分析

從群體遺傳學(xué)角度對宗地花豬BMPR-IB基因各SNP位點的基因型和等位基因頻率進行分析。由表2可見，在G595C位點上，AG 基因型頻率最高， 為優(yōu)勢基因型， A等位基因頻率為52.27%，表現(xiàn)為優(yōu)勢等位基因；對于C643T位點，優(yōu)勢基因型為BB基因型，B為優(yōu)勢等位基因。通過對

2.5不同基因型與繁殖性能的關(guān)聯(lián)分析

由表3可知，在宗地花豬母豬群中，總產(chǎn)仔數(shù)、產(chǎn)活仔數(shù)各基因型均表現(xiàn)為AA﹥AG﹥GG，AA基因型為優(yōu)勢基因型，總產(chǎn)仔數(shù)和產(chǎn)活仔數(shù)均最多，A為優(yōu)勢等位基因。方差分析結(jié)果表明，宗地花豬母豬總產(chǎn)仔數(shù)表現(xiàn)為AA基因型個體顯著高于AG、GG基因型個體，分別多產(chǎn)1.13、1.37頭/胎，AG基因型個體與GG基因型個體間總產(chǎn)仔數(shù)、產(chǎn)活仔數(shù)差異均不顯著；AA基因型個體產(chǎn)活仔數(shù)顯著高于AG、GG基因型個體，分別多產(chǎn)1.15、1.56頭/胎；GG基因型的仔豬平均初生質(zhì)量最高，GG基因型為優(yōu)勢基因型，3種基因型個體間仔豬平均初生質(zhì)量差異不顯著；產(chǎn)活仔數(shù)為7頭/胎的母豬泌乳力表現(xiàn)為AA基因型個體與AG、GG基因型個體間差異極顯著，分別比AG、GG基因型高17.297、22.917 kg。由表3還可知，在宗地花豬母豬群中，總產(chǎn)仔數(shù)和產(chǎn)活仔數(shù)各基因型均表現(xiàn)為BB﹥BT﹥TT，BB基因型為優(yōu)勢基因型，總產(chǎn)仔數(shù)和產(chǎn)活仔數(shù)均最多，B為優(yōu)勢等位基因。方差分析結(jié)果表明，宗地花豬母豬總產(chǎn)仔數(shù)表現(xiàn)為BB基因型個體顯著高于BT、TT基因型個體，分別多產(chǎn) 1.06、2.06頭/胎，BT基因型個體顯著高于TT基因型個體，多產(chǎn)1.00頭/胎，BB基因型個體產(chǎn)活仔數(shù)顯著高于BT、TT基因型個體，分別多產(chǎn)1.04、1.56頭/胎；仔豬平均初生質(zhì)量各基因型表現(xiàn)為TT﹥BB﹥BT，但各基因型間差異不顯著；產(chǎn)活仔數(shù)為7頭/胎的母豬泌乳力表現(xiàn)為BB基因型個體高于BT基因型個體，由于樣本數(shù)有限，無法判斷該位點TT基因型個體的平均泌乳力與BB、BT基因型之間的差異。

3討論

本試驗在宗地花豬母豬BMPR-IB基因第6外顯子上發(fā)現(xiàn)2個SNP位點，分別是分別是G595C和C643T，2個位點均為錯義突變，G595C導(dǎo)致編碼的天冬氨酸（aspartic acid，簡稱Asp）變?yōu)榻M氨酸（histidine，簡稱His），C643T位點導(dǎo)致編碼的精氨酸（arginine，簡稱Arg）變?yōu)樯彼幔╰ryptophan，簡稱Trp）。該基因G595C位點在宗地花豬群體中等位基因A、G的基因頻率分別為52.27%、47.73%，多態(tài)信息含量為 0.374 5，屬于中度多態(tài)（0.250.05）。結(jié)果表明，G595C位點在該試驗群體的選擇中均不易受到外界人為因素的干擾。同時這也表明，宗地花豬在該多態(tài)位點的遺傳變異較高，在以后的工作中應(yīng)加強后備母豬基礎(chǔ)群的選育，進行有效的選擇與選育。該基因C643T位點在宗地花豬群體中等位基因B、T的基因頻率分別為83.33%、16.67%，等位基因B的頻率明顯高于等位基因T，該位點多態(tài)信息含量為 0.239 2，屬于低度多態(tài)（PIC<0.25），該位點在試驗群體中基因突變頻率低，該位點可提供的遺傳信息性較差，經(jīng)卡方檢驗，該位點偏離了Hardy-Weinberg平衡（P<0.05），這可能是由該飼養(yǎng)場的宗地花豬群體樣本數(shù)不夠大引起的。C643T位點屬于低度多態(tài)，表明在以后宗地花豬繁殖性能的選種選育時，還須進一步擴大試驗群體進行試驗。

結(jié)果表明，在G595C位點上，宗地花豬母豬總產(chǎn)仔數(shù)、個體產(chǎn)活仔數(shù)均表現(xiàn)為AA基因型個體顯著高于AG、GG基因型個體，產(chǎn)活仔數(shù)為7頭/胎的母豬泌乳力表現(xiàn)為AA基因型個體與AG、GG基因型個體間差異極顯著（P<0.01）。綜合本試驗繁殖性狀指標，AA基因型母豬產(chǎn)仔數(shù)高，且哺乳能力強。在C643T位點上，宗地花豬母豬總產(chǎn)仔數(shù)表現(xiàn)為BB基因型個體顯著高于BT、TT基因型個體，BT基因型個體總產(chǎn)仔數(shù)又顯著高于TT基因型個體（P<0.05）；BB基因型個體產(chǎn)活仔數(shù)顯著高于BT、TT基因型個體。綜合本試驗繁殖性狀指標，BB基因型的母豬為高產(chǎn)母豬。

結(jié)果表明，BMPR-IB基因的G595C、C643T突變，均導(dǎo)致突變型經(jīng)產(chǎn)母豬的產(chǎn)仔數(shù)下降，初步說明BMPR-IB基因的突變基因型可能會對經(jīng)產(chǎn)母豬的產(chǎn)仔數(shù)有一定的抑制作用，這與宋一萍的研究結(jié)果[4]一致。Guan等研究表明，F(xiàn)ecB基因?qū)Ω嵫虻某跎|(zhì)量沒有顯著影響[15]；楊華等對雜交羊的研究結(jié)果[16]與Guan等的結(jié)果一致。目前，還未見BMPR-IB基因?qū)ψ胸i初生質(zhì)量的影響報道，本試驗結(jié)果表明，BMPR-IB基因?qū)ψ胸i初生質(zhì)量沒有顯著影響。BMPR-IB基因?qū)ψ诘鼗ㄘi母豬泌乳力有顯著影響（P<0.05），但目前還未見BMPR-IB基因與母豬泌乳力關(guān)系的相關(guān)性報道。

由于母豬的繁殖性狀受到多方面因素的影響，在生產(chǎn)實踐過程中，除了品種、品系，母豬的斷奶時間、發(fā)情時間、配種時間、母豬的斷奶膘情、營養(yǎng)水平和飼養(yǎng)管理都會影響母豬的繁殖性狀[17-20]。本試驗雖在相同的營養(yǎng)水平和飼養(yǎng)管理條件下進行，但母豬自身的膘情、發(fā)情期、配種期，以及配種公豬的不同等均影響著試驗結(jié)果，這造成母豬繁殖性狀數(shù)據(jù)的收集存在很大差異，使結(jié)果不夠準確。下一步將繼續(xù)擴大試驗群體數(shù)據(jù)的采集，進一步提供宗地花豬母豬繁殖性能的有效標記。

參考文獻：

[1]Rosenzweig B L，Imamura T，Okadome T，et al. Cloning and characterization of a human type Ⅱ receptor for bone morphogenetic proteins[J]. Proceedings of the National Academy of Sciences，1995，92（17）：7632-7636.

[2]楊燕燕. 影響蒙古羊雙羔群體繁殖性狀候選基因的研究[D]. 呼和浩特：內(nèi)蒙古農(nóng)業(yè)大學(xué)，2009.

[3]Quinn R L，Shuttleworth G，Hunter M G. Immunohistochemical localization of the bone morphogenetic protein receptors in the porcine ovary[J]. Journal of Anatomy，2004，205（1）：15-23.

[4]宋一萍. 豬BMP15和BMPR-IB基因的遺傳變異及其與產(chǎn)仔性能關(guān)系的研究[D]. 泰安：山東農(nóng)業(yè)大學(xué)，2010.

[5]孟洪莉. 豬繁殖相關(guān)基因多態(tài)性檢測及其關(guān)聯(lián)分析[D]. 福州：福建農(nóng)林大學(xué)，2012.

[6]馮文豪，燕志宏，毛同輝. 貴州宗地花豬胴體肉質(zhì)性狀測定初報[J]. 貴州畜牧獸醫(yī)，2008，32（2）：9-10.

[7]燕志宏，田松軍，申紅艷，等. 宗地花豬1世代的生長發(fā)育研究[J]. 安徽農(nóng)業(yè)科學(xué)，2010，38（24）：13107-13108.

[8]燕志宏，張蕓，田松軍. 宗地花豬及其雜種1代豬哺乳期和保育期生長性能的研究初報[J]. 養(yǎng)豬，2011（3）：47-48.

[9]燕志宏，毛同輝，馮文豪. 宗地花豬ESR、FSHβ基因的多態(tài)性及其對繁殖性能的影響[J]. 貴州農(nóng)業(yè)科學(xué)，2010，38（7）：134-137.

[10]燕志宏，田松軍，韋玉美. 貴州宗地花豬最佳屠宰體重的研究[J]. 豬業(yè)科學(xué)，2010，27（7）：102-104.

[11]李平，燕志宏，陳偉，等. 宗地花豬H-FABP基因mRNA的表達與肉質(zhì)性狀的相關(guān)研究[J]. 中國畜牧獸醫(yī)，2015，42（3）：644-649.

[12]顧麗菊，任麗群，燕志宏，等. 宗地花豬和從江香豬ADRP基因多態(tài)性及生物信息學(xué)分析[J]. 基因組學(xué)與應(yīng)用生物學(xué)，2015，34（11）：2395-2401.

[13]任麗群，燕志宏，顧麗菊，等. 貴州宗地花豬的GDF9基因多態(tài)性及生物信息學(xué)分析[J]. 基因組學(xué)與應(yīng)用生物學(xué)，2015，34（11）：2402-2407.

[14]閆雷，林家棟，燕志宏，等. 宗地花豬生長激素基因TaqI多態(tài)性及其對屠宰性狀的影響[J]. 貴州農(nóng)業(yè)科學(xué)，2011，39（5）：170-172.

[15]Guan F，Liu S R，Shi G Q，et al. Polymorphism of FecB gene in nine sheep breeds or strains and its effects on litter size，lamb growth and development[J]. Animal Reproduction Science，2007，99（1）：44-52.

[16]楊華，楊永林，何其宏. 綿羊BMPR-IB基因多態(tài)性及其與杜泊羊、陶賽特羊和湖羊雜交F2羔羊初生重的相關(guān)性分析[J]. 中國草食動物科學(xué)，2012（增刊1）：208-210.

[17]付林. 影響母豬繁殖性能的因素分析[J]. 黑龍江畜牧獸醫(yī)，2013（22）：49-52.

[18]張茂，鄒發(fā)蘭，林長春，等. 二元母豬不同時期背膘對繁殖性能的影響[J]. 江蘇農(nóng)業(yè)科學(xué)，2016，44（12）：272-274.

[19]胡新旭，卞巧，劉維綱，等. 影響母豬繁殖性能的因素及其作用機理的研究進展[J]. 中國飼料，2015（4）：12-14，20.

[20]劉笛秋，郭盈魁，賀國安. 影響母豬繁殖性能的因素及解決措施[J]. 中國豬業(yè)，2013（6）：41-43.邵珠田，姜立娜，郭小菲，等. 45份菜用大黃的SRAP標記引物篩選及親緣關(guān)系分析[J]. 江蘇農(nóng)業(yè)科學(xué)，2017，45（17）：39-42.

doi：10.15889/j.issn.1002-1302.2017.17.010