楊 酸,郭小江,楊紅文,熊 力,李 晨,譚元成,王春源,張依裕,*
(1.貴州大學 動物科學學院,高原山地動物遺傳育種與繁殖教育部重點實驗室,貴州 貴陽550025; 2.貴州省農(nóng)業(yè)農(nóng)村廳,貴州 貴陽550001)
柯樂豬屬云貴高原烏蒙烏金豬類群,是貴州省優(yōu)良的地方豬種質(zhì)資源,具有肉質(zhì)優(yōu)良、耐粗飼、適宜放養(yǎng)、生長速度緩慢、繁殖性能差等特點,是烏蒙喀斯特山區(qū)真正的“烏金豬”,也是貴州省不可多得的寶貴豬種資源[1-2]。一直以來,當?shù)匮芯克鶎聵坟i進行了種質(zhì)資源的保種、選育和利用等研究工作,但均未開展柯樂豬品系選育,研究缺乏系統(tǒng)性。目前對于該豬種的研究集中在肉質(zhì)性狀[3-4]、毛色[5]、飼養(yǎng)方式[6]和耐粗飼特性[7]等方面,在繁殖性狀方面的研究甚少。近年來,豬的繁殖性狀逐漸作為重點改良目標性狀進而快速提高豬的重要經(jīng)濟性狀,且母豬繁殖性狀的改良已是提高母豬年斷奶仔豬數(shù)的首要考慮因素,于母豬的繁殖性能而言,產(chǎn)仔數(shù)是評定其繁殖力的一個最重要的指標。也有研究表明,增加窩產(chǎn)活仔數(shù)可顯著提高初生窩重和斷奶窩重[8],初生窩重又是斷奶窩重的基礎,斷奶窩重是母豬繁殖性能各項指標的總體現(xiàn)[9]。這些常見的繁殖性能指標遺傳力偏低[10],受不同因素的影響,人們偏向于采用全基因組關(guān)聯(lián)研究來探究豬繁殖性狀的遺傳結(jié)構(gòu)相關(guān)問題[11],其中,單核苷酸多態(tài)性(single nucleotide polymorphism, SNP)被用來評估動物之間的親緣關(guān)系,即使用全基因組SNP面板作為標記和統(tǒng)計方法,同時捕捉大量SNP的影響[12]。目前,大多數(shù)SNP發(fā)現(xiàn)的候選基因或QTL(數(shù)量性狀位點)已被確定,這些數(shù)據(jù)將成為突破性育種技術(shù)的推動力,并提高養(yǎng)豬業(yè)的生產(chǎn)效率。由此,對于產(chǎn)仔數(shù)少、繁殖力低的地方品種,從分子層面探究遺傳性狀的差異,并確定控制表型性狀的基因或與該基因緊密連鎖的遺傳標記,對該品種直接進行基因型選擇或標記輔助選擇,能夠有效提高地方品種的選育效率。
催乳素(prolactin,PRL)是一種蛋白質(zhì)激素,具有多種生物學功能,包括適應性應激反應、滲透調(diào)節(jié)、子宮功能、乳腺的增殖和分化、泌乳等方面[13]。PRL通過其多種形式的膜相關(guān)受體發(fā)揮作用,催乳素受體(prolactin receptor,PRLR)基因較長,其總長度大于100 kb,至少由10個外顯子和9個內(nèi)含子組成,且第一外顯子是可變的[14]。PRLR基因早已被證明與母豬繁殖性狀有顯著關(guān)聯(lián)[15],在大白豬、長白豬、杜洛克豬[16]、藏豬[17]、波蘭豬[18]、定遠豬[19]和松遼黑豬[20]等多個中外豬種上針對其繁殖性狀已有不同程度的研究。PRLR是否能作為柯樂豬在繁殖性能方面的關(guān)鍵候選基因至今未能證明。因此,本文以柯樂豬為研究對象,利用Sanger直接測序法查找PRLR基因的SNP位點,并與柯樂豬的7個繁殖性狀做關(guān)聯(lián)性分析,進一步尋找與柯樂豬繁殖性狀相關(guān)的遺傳標記,為其繁殖性能改良和在養(yǎng)豬生產(chǎn)中的應用提供科學依據(jù)。
試驗豬群來自貴州省赫章縣柯樂豬保種場,隨機選擇柯樂豬母豬204頭,采集約0.25 g耳組織,置于75%乙醇的離心管中,放入冰盒帶回實驗室,轉(zhuǎn)入-20 ℃冰箱保存待用。記錄該豬群母豬頭胎的總產(chǎn)仔數(shù)、產(chǎn)活仔數(shù)、初生窩重、初生個體重、斷奶窩仔數(shù)、斷奶窩重和斷奶個體重等數(shù)據(jù)。
參照陳文娟[21]的苯酚-氯仿-異戊醇抽提法提取母豬耳組織DNA,用1%瓊脂糖凝膠電泳和核酸濃度檢測儀(NanoDrop2000)聯(lián)合檢測基因組DNA的濃度和純度,置于-20 ℃保存。
根據(jù)NCBI中豬PRLR基因的DNA參考序列(登錄號:NC_010458.4),結(jié)合在線軟件Primer Premier 3.0(http://primer3.ut.ee/)設計豬PRLR引物序列(表1),P1引物位置在第2外顯子內(nèi),P2引物位置在第8外顯子內(nèi),P3引物的位置在第3內(nèi)含子和第4外顯子之間,P4引物的位置在第6、第7外顯子,以及第6、第7部分內(nèi)含子之間。引物由生工生物工程(上海)股份有限公司合成。
表1 PRLR基因擴增引物信息
PCR反應體系(20 μL):1 μL DNA,10 μL 2×TaqPCR master mix,上、下游引物各1 μL (10 mmol·L-1),7 μL ddH2O。PCR反應程序:95 ℃預變性5 min;95 ℃變性30 s,最適溫度退火40 s,72 ℃延伸30 s,共35個循環(huán);72 ℃延伸5 min,最后4 ℃保存。PCR產(chǎn)物經(jīng)1.2%的瓊脂糖凝膠電泳檢測,將檢測合格的樣本送生工生物工程(上海)股份有限公司進行雙向測序,測序結(jié)果采用DNAStar軟件中的MegAlign和Editseq程序進行SNP位點篩選,同時進行基因分型。
通過RNAfold web server (http://rna.tbi.univie.ac.at/cgi-bin/RNAWebSuite/RNAfold.cgi)預測mRNA二級結(jié)構(gòu),并對比參考序列和柯樂豬PRLR基因SNPs的mRNA二級結(jié)構(gòu);利用Excel 2010軟件計算PRLR基因各個位點的基因型頻率、等位基因頻率、遺傳雜合度 (heterozygosity, He)、有效等位基因數(shù) (number of effective alleles, Ne)、多態(tài)信息含量 (polymorphism information content, PIC)和Hardy-Weinberg平衡;用SHEsis在線軟件 (http://analysis.bio-x.cn/)對SNP位點連鎖不平衡(linkage disequilibrium, LD)和單倍型進行計算分析;采用SPSS 22.0軟件的一般線性模型(general linearmodel, GLM)分析PRLR基因SNPs與繁殖性狀的關(guān)聯(lián)性,以及PRLR基因SNPs雙倍型與繁殖性狀的差異顯著性,分析時采用以下模型:性狀觀察值=群體均值+基因型效應+隨機殘差效應,再用Duncan’s法進行多重比較,結(jié)果用平均值±標準誤表示,P<0.05表示差異顯著。
瓊脂糖凝膠電泳檢測結(jié)果表明,4對引物的PCR擴增條帶清晰,無冗余雜帶(圖1),與預設片段大小相符,可用于進一步分析。
M,DL 2000 DNA marker;1~3,P1擴增產(chǎn)物;4~6,P2擴增產(chǎn)物;7~9,P3擴增產(chǎn)物;10~12,P4擴增產(chǎn)物。
在柯樂豬PRLR基因第8外顯子檢測到1個SNP突變位點(g.20655220 C>T)(圖2),在第4內(nèi)含子檢測到2個SNP突變位點(g.20648859 C>T和g.20648931 C>T),3個SNP位點均存在3種基因型,即CC、CT、TT。其中g(shù).20655220 C>T突變導致絲氨酸變成脯氨酸,為錯義突變。
柯樂豬PRLR基因的多態(tài)位點遺傳多態(tài)性分析結(jié)果(表2)表明,g.20655220 C>T、g.20648859 C>T和g.20648931 C>T優(yōu)勢基因型均為CC,其頻率分別是0.583、0.657和0.755,優(yōu)勢等位基因均為C,頻率分別是0.762、0.784和0.853;遺傳指數(shù)分析中,g.20655220 C>T和g.20648859 C>T屬于中度多態(tài)(0.25
表2 柯樂豬PRLR基因SNPs群體遺傳信息
利用在線軟件SHEsis對檢測到的3個SNPs位點進行連鎖不平衡分析(表3),根據(jù)D′值和相關(guān)系數(shù)(r)判斷位點之間是否存在連鎖,位點g.20655220 C>T與g.20648859 C>T之間的D′值為0.111,r為0.001;g.20655220 C>T與g.20648931 C>T之間的D′值為0.022,r為0,g.20648859 C>T與g.20648931 C>T之間D′值為1.000,r為0.627。由此可得,g.20648859 C>位點與g.20648931 C>T位點之間有較強的連鎖不平衡效應(D′>0.85,r>0.33)。
表3 PRLR基因SNP位點的連鎖不平衡分析
對柯樂豬PRLR基因3個SNP位點進行單倍型和雙倍型分析,結(jié)果(表4)顯示,3個SNP位點在204頭柯樂豬群體中共檢測到5種單倍型和9種雙倍型,其中,優(yōu)勢單倍型是H1(CCC),頻率為0.610;劣勢單倍型是H5(TTT),頻率為0.049;優(yōu)勢雙倍型是H1H1(CCCCCC),頻率為0.382,劣勢雙倍型是H2H3(TTTTTT),頻率為0.039。
表4 PRLR基因單倍型和雙倍型分析
RNAfold web server預測結(jié)果表明,g.20655220 C>T導致PRLR基因mRNA序列二級結(jié)構(gòu)發(fā)生一定改變,且mRNA二級結(jié)構(gòu)最小自由能有所提高,由-648.40 kcal·mol-1變?yōu)?647.90 kcal·mol-1,可能會使得mRNA二級結(jié)構(gòu)穩(wěn)定性降低。其突變前跟突變后的mRNA二級結(jié)構(gòu)見圖3。
圖3 柯樂豬PRLR基因mRNA二級結(jié)構(gòu)
利用SPSS18.0軟件中的單因素方差分析來分析PRLR基因3個SNPs位點與柯樂豬繁殖性狀的關(guān)聯(lián)性,結(jié)果(表5)顯示,g.20655220 C>T位點CC基因型母豬的初生個體重、斷奶窩重和斷奶個體重均顯著(P<0.05)高于CT基因型,但與TT基因型差異不顯著(P>0.05);g.20648859 C>T位點不同基因型在各繁殖性狀上均無顯著差異(P>0.05);g.20648931 C>T位點CT基因型母豬的初生個體重、初生窩重和斷奶窩重均顯著(P<0.05)高于CC基因型和TT基因型,且該位點在初生個體重上,CC基因型顯著(P<0.05)高于TT基因型。
表5 柯樂豬PRLR基因SNP位點與繁殖性狀的關(guān)聯(lián)性
將PRLR基因3個SNPs組成的雙倍型與可樂豬繁殖性狀進行關(guān)聯(lián)性分析,結(jié)果(表6)顯示,雙倍型H1H5(CTCTCT)個體的總產(chǎn)仔數(shù)顯著(P<0.05)多于H2H4(CTCTCC)雙倍型個體,且產(chǎn)活仔數(shù)均多于其余雙倍型;雙倍型H1H3(CCCTCT)個體的初生窩重和初生個體重均顯著(P<0.05)高于H1H4(CTCCCC)個體和H3H5(CTTTTT)個體,雙倍型H2H3(CCTTCT)的初生個體重、斷奶窩重和斷奶個體重均顯著(P<0.05)高于H3H5(CTTTTT)。
表6 柯樂豬PRLR基因3個SNP位點雙倍型與繁殖性狀的關(guān)聯(lián)性
與大部分國內(nèi)外豬種相比,柯樂豬具有產(chǎn)仔數(shù)低、性成熟早、繁殖率低的生理特性,經(jīng)采集,本試驗研究群體的平均總產(chǎn)仔數(shù)和產(chǎn)活仔數(shù)分別為8.63頭和7.92頭,而近些年雜交豬種的頻繁應用,地方豬種的生存空間進一步縮小,大概是由于地方豬種繁殖力較低下,且生長緩慢、育種周期更長[22],因此,提高地方豬種繁殖力是當前研究的一個重點及熱點。PRLR基因作為母豬繁殖性狀的重要候選基因之一,在大腦中的廣泛表達突出了催乳素在調(diào)節(jié)神經(jīng)元功能中的多效性作用,能夠促進乳腺發(fā)育、維持泌乳功能,并對家畜性腺發(fā)育和母體行為產(chǎn)生影響[23]。本研究在柯樂豬PRLR基因第8外顯子檢測到1個SNP突變位點(g.20655220 C>T),在第4內(nèi)含子檢測到2個SNP突變位點(g.20648859 C>T和g.20648931 C>T),3個位點呈現(xiàn)不同的多態(tài)性:g.20655220 C>T和g.20648859 C>T屬于中度多態(tài)(0.25
PRLR基因定位于豬的16位染色體,與母豬的繁殖性狀密切相關(guān),已有大量學者對PRLR基因多態(tài)性和母豬產(chǎn)仔性能展開研究。Sun等[26]研究表明,PRLR基因在5個山東地方豬種和3個國外豬種的受試群體中存在多態(tài)性位點,且產(chǎn)仔數(shù)性狀的效應明顯,遺傳多態(tài)性豐富,并證實了該基因比RBP4基因的遺傳效應更大;Birgitte等[27]早在2002年就已證實PRLR基因是大白×梅山F2雜交母豬初情年齡、產(chǎn)活仔數(shù)和產(chǎn)仔數(shù)的主效基因或標記。本研究在柯樂豬PRLR基因上新發(fā)現(xiàn)的3個多態(tài)位點中,g.20655220 C>T位點對柯樂豬母豬初生個體重、斷奶窩重和斷奶個體重的影響達顯著性水平,該位點的CC基因型母豬的初生個體重、斷奶窩重和斷奶個體重均顯著高于CT基因型。與此類似,Vashi等[28]研究顯示,PRLRCC基因型與印度本土豬的斷奶窩重呈正相關(guān);范一萍等[29]在大白豬群體上的驗證顯示,其純合基因型BB基因型個體的斷奶窩重顯著高于雜合基因型AB基因型,且其內(nèi)含子上的g.20648931 C>T位點對其初生個體重、初生窩重和斷奶窩重的影響也達到顯著水平,其中,g.20648931 C>T位點CT基因型母豬的初生個體重、初生窩重和斷奶窩重均顯著高于CC基因型和TT基因型,且該位點在初生個體重上,CC基因型又顯著高于TT基因型。出現(xiàn)此類情況的原因可能與樣本量較少,且均為頭胎繁殖數(shù)據(jù)相關(guān),尤其在地方豬種群體中,普遍呈現(xiàn)性成熟較早的情況,而有些后備母豬的身體可能還沒有完全發(fā)育成熟,從而造成頭胎繁殖數(shù)據(jù)的一定局限性。因此,后期有必要進一步展開研究,對經(jīng)產(chǎn)母豬繁殖性狀進行關(guān)聯(lián)分析。
在產(chǎn)仔數(shù)方面,本研究中PRLR基因3個SNPs位點與柯樂豬產(chǎn)仔數(shù)和產(chǎn)活仔數(shù)的關(guān)聯(lián)性均未達到顯著水平,這與胡慧艷等[30]的研究結(jié)論相同,其表示PRLR基因?qū)Ξa(chǎn)仔數(shù)不存在顯著影響。將PRLR基因3個SNPs組成的雙倍型與可樂豬繁殖性狀進行關(guān)聯(lián)性分析發(fā)現(xiàn),H1H5(CTCTCT)雙倍型個體在總產(chǎn)仔數(shù)上顯著高于H2H4(CTCTCC)雙倍型個體,且其產(chǎn)活仔數(shù)均高于其余雙倍型,由此初步認為雙倍型H1H5(CTCTCT)為有利雙倍型,在產(chǎn)仔數(shù)的培育過程中可作為分子標記輔助選擇的參考。另外,雙倍型H1H3(CCCTCT)個體的初生窩重和初生個體重均顯著高于H1H4(CTCCCC)個體和H3H5(CTTTTT)個體;在斷奶窩重和斷奶個體重方面,雙倍型H2H3(CCTTCT)為最優(yōu)組合,均高于其他雙倍型,且在初生個體重、斷奶窩重和斷奶個體重上均顯著高于H3H5(CTTTTT)。綜上所述,g.20655220 C>T位點的CC基因型是改善柯樂豬母豬繁殖性狀,特別是初生個體重、斷奶窩重和斷奶個體重的有利基因型。雙倍型H1H5(CTCTCT)是產(chǎn)仔數(shù)和產(chǎn)活仔數(shù)的有利雙倍型,H2H3(CCTTCT)在斷奶窩重和斷奶個體重上為較優(yōu)雙倍型。