賈美婷,李顯耀,康 麗,姜運良,樊新忠,唐 輝
(山東農(nóng)業(yè)大學動物科技學院,泰安 271018)
?
壽光雞GPR147基因5′調(diào)控區(qū)多態(tài)性與產(chǎn)蛋性狀的關(guān)聯(lián)分析
賈美婷,李顯耀,康麗,姜運良,樊新忠,唐輝*
(山東農(nóng)業(yè)大學動物科技學院,泰安 271018)
摘要:為了了解GPR147基因多態(tài)性及其與產(chǎn)蛋性狀的關(guān)系,本試驗以壽光雞為試驗材料,用直接測序的方法檢測了壽光雞該基因5′調(diào)控區(qū)的SNPs,并與產(chǎn)蛋性狀進行了關(guān)聯(lián)分析。結(jié)果發(fā)現(xiàn),GPR147基因的5′ 調(diào)控區(qū)在雞群中共有11個SNPs位點,其中-593和-579位點的不同基因型在后期(45~56周)產(chǎn)蛋量和56周總產(chǎn)蛋數(shù)有顯著差異(P<0.05);-593位點TT基因型對應較高產(chǎn)蛋數(shù),-579位點AA基因型對應較高產(chǎn)蛋數(shù)。上述兩個位點的單倍型分析結(jié)果表明,單倍型TA與CG相比,中期(26~44周)產(chǎn)蛋量多了1.5枚(P=0.041 8),后期產(chǎn)蛋量多2.5枚(P=0.032 6),300天產(chǎn)蛋量多1.6枚(P=0.1),56周總產(chǎn)蛋量多4.4枚(P=0.014 4)。綜上,GPR147 5′ 調(diào)控區(qū)存在豐富的SNPs位點,-593和-579位點SNPs及其單倍型對壽光雞的產(chǎn)蛋性狀有顯著影響,是具有潛在應用價值的分子標記。
關(guān)鍵詞:雞;產(chǎn)蛋性能;GPR147基因;多態(tài)性
促性腺激素抑制激素(GnIH)是2000年K.Tsutsui等在日本鵪鶉中發(fā)現(xiàn)的一種下丘腦十二肽,能夠抑制促性腺激素的合成和釋放[1]。研究顯示,GnIH廣泛存在于鳥類、魚類及哺乳動物的下丘腦——垂體——性腺軸(HPG軸)上[2-5],通過多種方式參與生殖調(diào)控[6-7]。T.Ikemoto等和J.A.Bonini等報道雞的GnIH有兩種受體(RFRPR和NPFFR),其中RFRPR又稱為GPR147基因,主要在垂體表達[8-9]。GnIH與其受體GPR147結(jié)合,抑制腺苷酸環(huán)化酶活性,從而降低細胞內(nèi)cAMP水平,干擾G蛋白耦聯(lián)受體(包括GnRH、LH、FSH受體)信號傳導[10]。T.Osugi等研究發(fā)現(xiàn),給麻雀注射GnIH能快速減少血漿LH水平,同時注射GnIH和雞促性腺激素釋放激素I(cGnRH-I),能降低GnRH誘導的血漿LH的升高水平[2]。T.Ubuka等研究表明,白冠雀和鵪鶉腹內(nèi)注射GnIH后,LH濃度降低,鵪鶉LHα、LHβ和FSHα亞單位的mRNA表達量顯著減少[11]。因此推測,GnIH可以通過直接或間接作用抑制促卵泡激素(FSH)和促黃體激素(LH)的合成與釋放,繼而影響雌激素的合成。這些激素共同作用影響生殖器官的發(fā)育,進而對蛋雞的產(chǎn)蛋水平產(chǎn)生影響。
近年來,有關(guān)GnIH及其受體GPR147的研究主要停留在神經(jīng)內(nèi)分泌調(diào)控上,對雞GPR147基因5′調(diào)控區(qū)多態(tài)性及其與產(chǎn)蛋性狀關(guān)系的研究還未見報道。本研究以壽光雞為試驗材料,研究其GPR147基因5′調(diào)控區(qū)的SNPs,并與產(chǎn)蛋性狀進行關(guān)聯(lián)分析,為地方雞產(chǎn)蛋性狀分子標記輔助選擇提供依據(jù)。
1材料與方法
1.1材料
1.1.1試驗動物本試驗雞群來自山東壽光慈倫大雞養(yǎng)殖場,選取有完整的開產(chǎn)日齡及每日產(chǎn)蛋情況記錄的樣本432只。所有參試雞群單籠飼養(yǎng),飼養(yǎng)管理條件相同。
1.1.2主要試劑PJET1.2 Vector、DNA凝膠回收試劑盒購自天根生化科技有限公司,感受態(tài)細胞E.coliDH5α、PrimerSTAR HS DNA Polymerase購自TaKaRa公司。
1.2DNA提取
50周齡時經(jīng)翅靜脈采血后,根據(jù)血液基因組DNA提取試劑盒說明書提取DNA。用紫外分光光度計和1%瓊脂糖凝膠電泳檢測DNA的濃度與純度,于-20 ℃保存。
1.3引物設計
根據(jù)GenBank上發(fā)表的雞GPR147基因序列(登錄號:NC_006093.3),在5′側(cè)翼序列設計引物。引物由生工生物工程(上海)股份有限公司合成,序列:正向引物(F):5′-GCTGCTACTGTCACCTGA-GG-3′;反向引物(R):5′-GGAGAAACCGCAATGAAAGGAG-3′。
1.4PCR擴增和SNPs位點檢測
根據(jù)產(chǎn)蛋記錄,選取30個產(chǎn)蛋性狀差異較大個體的DNA樣品,混池后進行PCR擴增。PCR反應條件:94 ℃預變性5 min;98 ℃變性10 s,63.5 ℃退火15 s,72 ℃延伸1 min,35個循環(huán);72 ℃延伸10 min;4 ℃保存。
PCR產(chǎn)物經(jīng)瓊脂糖凝膠電泳檢測后,用DNA凝膠回收試劑盒回收,并與PJET1.2 Vector連接。將連接產(chǎn)物轉(zhuǎn)化至感受態(tài)細胞E.coliDH5α,通過菌液PCR鑒定陽性克隆,將陽性克隆菌送至北京華大基因有限公司進行測序,檢測多態(tài)性位點。用TFSEARCH在線軟件預測可能的轉(zhuǎn)錄因子結(jié)合位點(可能性大于85%)。
1.5多態(tài)性位點的群體檢測與基因分型
每個個體PCR擴增GPR147基因5′調(diào)控區(qū)序列,PCR產(chǎn)物直接送北京華大基因有限公司測序,檢測其多態(tài)位點基因型。
1.6數(shù)據(jù)處理
根據(jù)產(chǎn)蛋率是否高于60%,將產(chǎn)蛋期分為前期(20~25周齡)、中期(26~44周齡)和后期(45~56周齡)3個階段,統(tǒng)計各階段每只雞產(chǎn)蛋數(shù),300 天產(chǎn)蛋數(shù)及產(chǎn)蛋全期總產(chǎn)蛋數(shù)。
用R軟件的Genetics程序包來計算分析各個多態(tài)位點的群體遺傳學特征。用SAS8.2軟件包的GLM程序?qū)⒉煌鄳B(tài)位點及其單倍型與產(chǎn)蛋性狀進行關(guān)聯(lián)分析,統(tǒng)計模型:Yij=μ+Gi+eij,其中,Yij為性狀表型值,μ為該性狀的總體均值,Gi為基因型效應值,eij為隨機殘差效應[12]。計算基因的加性效應、顯性效應,并進行t檢驗分析。用PHASE軟件(V2.0)估計每個個體多態(tài)位點的單倍型。
2結(jié)果
2.1GPR147基因的多態(tài)性
混池模板PCR擴增后,用1%瓊脂糖凝膠電泳檢測,結(jié)果如圖1所示,片段大小與911 bp相符,特異性好。單克隆測序結(jié)果發(fā)現(xiàn),壽光雞GPR147基因5′調(diào)控區(qū)有11個SNP位點:-660(G→A)、-602(C→T)、-593(T→C)、-579(A→G)、-562(C→A)、-523(T→G)、-516(A→G)、-489(A→G)、-435(C→T)、-427(T→G/C)、-357(G→C/A)。根據(jù)測序結(jié)果可以判斷個體基因型,單峰為純合子,雙峰為雜合子[13](圖2)。
各位點用TFSEARCH在線軟件預測,發(fā)現(xiàn)突變后-660位點失去CDP CR轉(zhuǎn)錄因子結(jié)合位點(得分為88.1),-602位點失去USF轉(zhuǎn)錄因子結(jié)合位點(得分為86.6),-579位點失去GATA-1(得分為91.0)和GATA-2(得分為90.9)轉(zhuǎn)錄因子結(jié)合位點,-562、-489位點失去MZF1(得分分別為89.6、85.2)轉(zhuǎn)錄因子結(jié)合位點,-427位點失去RORalp轉(zhuǎn)錄因子結(jié)合位點(得分為87.5)。其余位點堿基突變后,轉(zhuǎn)錄因子結(jié)合位點沒有變化,其中-593位點的轉(zhuǎn)錄因子結(jié)合位點為GATA-1(得分為98.4)和GATA-2(得分為87.0)。
2.2GPR147基因5′調(diào)控區(qū)SNPs位點在壽光雞群體中的分布
GPR147基因5′調(diào)控區(qū)SNPs位點群體遺傳學分析結(jié)果見表1。χ2檢驗表明,-602、-562、-523三個多態(tài)位點處于Hardy-Weinberg平衡狀態(tài)(P>0.05),其余位點未達到Hardy-Weinberg平衡狀態(tài)(P<0.05)。-660、-593、-579位點的多態(tài)信息含量分別為0.25、0.27和0.26,呈中度多態(tài)(0.25≤PIC<0.5),其余位點均呈低度多態(tài)(PIC<0.25)[14]。
表1GPR147基因5′調(diào)控區(qū)SNPs位點的群體遺傳學分析結(jié)果
Table1The results of population genetic analysis on SNPs in 5′UTR ofGPR147
用SHEsis軟件分析各位點間的連鎖關(guān)系,發(fā)現(xiàn)-602、-562、-523三個位點完全連鎖(D′=1,r2=1)[14],-516、-489兩個位點完全連鎖(D′=1,r2=1),-593和-579位點不完全連鎖(D′=1,r2=0.928)。
2.3單位點多態(tài)性與壽光雞產(chǎn)蛋性狀的關(guān)聯(lián)分析
對壽光雞群體GPR147基因-593和-579兩個多態(tài)位點的單位點與產(chǎn)蛋性狀進行關(guān)聯(lián)分析,結(jié)果見表2。由表2可知,-593和-579位點的不同基因型在后期產(chǎn)蛋量和56周總產(chǎn)蛋數(shù)有顯著差異(P<0.05)。在-593位點,TT型個體后期產(chǎn)蛋數(shù)比CC型多4.3枚,比CT型多4.7枚;TT型個體56周產(chǎn)蛋數(shù)比CC型多6.1枚,比CT型多6.5枚。等位基因T與C相比能使雞后期產(chǎn)蛋數(shù)增加2.2枚,56周產(chǎn)蛋總數(shù)增加3.0枚。-579位點AA基因型對應較高產(chǎn)蛋數(shù),AA型個體后期產(chǎn)蛋數(shù)比AG型多4.1枚,比GG型多5.1枚;56周產(chǎn)蛋數(shù)AA型個體比AG型多4.4枚,比GG型多8.8枚;等位基因A與等位基因G相比能使雞的中期產(chǎn)蛋量增加1.5枚,后期產(chǎn)蛋數(shù)增加2.5枚,56周產(chǎn)蛋總數(shù)增加4.4枚。
其余低度多態(tài)位點,由于基因型頻率差別太大,分析結(jié)果沒有代表性,因此未作關(guān)聯(lián)分析。
表2不同位點基因型對產(chǎn)蛋性狀的最小二乘均數(shù)及等位基因的遺傳效應
Table 2Least squares means and allele genetic effects for laying traits at each locus
P<0.05表示該位點在該性狀不同基因型之間差異顯著;同一位點同一列標有不同字母的兩均數(shù)間差異顯著(P<0.05)。*.P<0.05;**.P<0.01。AE為加性效應的縮寫;DE為顯性效應的縮寫。下表同
P<0.05 indicates that different genotypes at this locus differ significantly in this trait.Value within the same column at the same locus with different superscripts differ significantly (P<0.05).*.P<0.05;**.P<0.01.AE.Additive effect;DE.Dominant effect.The same as below
2.4多位點聯(lián)合基因型及單倍型對壽光雞產(chǎn)蛋性能的效應分析
在該群體中,-593和-579位點有4種聯(lián)合基因型:TA/TA、TA/CG、CA/CG、CG/CG,其中CA/CG型只有10個個體,頻率低于5%,因此不作分析。對其他3種聯(lián)合基因型與產(chǎn)蛋性狀進行關(guān)聯(lián)分析,結(jié)果見表3。由表3可知,-593、-579位點不同聯(lián)合基因型對后期產(chǎn)蛋數(shù)和56周總產(chǎn)蛋量有顯著影響(P<0.05),表現(xiàn)為后期產(chǎn)蛋數(shù)TA/TA型比TA/CG型多4.7枚,比CG/CG多5.1枚,總產(chǎn)蛋數(shù)TA/TA型比TA/CG型多6.5枚,比CG/CG多8.8枚。加性效應及顯性效應分析結(jié)果發(fā)現(xiàn),單倍型TA與單倍型CG相比,中期產(chǎn)蛋量多1.5枚(P=0.041 8),后期產(chǎn)蛋量多2.5枚(P=0.032 6),300天產(chǎn)蛋量多1.6枚(P=0.1),56周總產(chǎn)蛋量多4.4枚(P=0.014 4)。
表3兩個位點聯(lián)合基因型的最小二乘均數(shù)及單倍型的遺傳效應
Table 3Least squares means by combined genotypes and haplotype genetic effects
3討論
本研究在壽光雞GPR147基因5′調(diào)控區(qū)發(fā)現(xiàn)了11個SNPs位點,其中3個為中度多態(tài),8個為低度多態(tài),并且哈代溫伯格平衡檢驗發(fā)現(xiàn)只有3個位點處于Hardy-Weinberg平衡狀態(tài)(P>0.05),原因可能與該群體的育種措施造成的偏離有關(guān)。在群體中發(fā)現(xiàn),-602、-562、-523位點只有兩種基因型,可能是該群體的人工選擇造成的。
本研究發(fā)現(xiàn),-579位點多態(tài)性對56周產(chǎn)蛋量和后期產(chǎn)蛋量有顯著影響(P<0.05)。TFSEARCH在線軟件預測發(fā)現(xiàn),-579位點由堿基A突變成堿基G后,失去GATA-1(得分為91.0)和GATA-2(得分為90.9)轉(zhuǎn)錄因子結(jié)合位點。哺乳動物研究表明,GATA家族對類固醇發(fā)生、性別決定和分化等過程有著不可或缺的作用[15]。小鼠GATA-2基因通過與PIT1相互作用,調(diào)節(jié)生長激素細胞、催乳素細胞、促性腺激素細胞和促甲狀腺激素細胞的生成,并能降低促性腺激素細胞中LHβ的表達[16-19]。因此-579位點很可能是因為失去GATA-1和GATA-2轉(zhuǎn)錄因子結(jié)合位點影響了基因表達。-593位點由堿基T突變成C,轉(zhuǎn)錄因子結(jié)合位點沒有變化,為GATA-1(得分為98.4)和GATA-2(得分為87.0),其多態(tài)性對產(chǎn)蛋性狀的影響可能是通過轉(zhuǎn)錄因子結(jié)合方式的改變或與其他位點(比如-579位點)連鎖造成的。
單位點多態(tài)性分析,-593位點等位基因T與C相比能使雞后期產(chǎn)蛋數(shù)增加2.2枚,56周產(chǎn)蛋總數(shù)增加3.0枚。-579位點等位基因A與G相比能使雞的中期產(chǎn)蛋量增加1.5枚,后期產(chǎn)蛋數(shù)增加2.5枚,56周產(chǎn)蛋總數(shù)增加4.4枚。單倍型分析發(fā)現(xiàn),單倍型TA與單倍型CG相比,中期產(chǎn)蛋量多1.5枚,后期產(chǎn)蛋量多2.5枚,56周總產(chǎn)蛋量多4.4枚。由此可以看出單倍型加性效應與-579位點的加性效應相同,進一步可以推測單倍型效應可能主要是-579位點的作用,-593位點的效應是連鎖作用的結(jié)果。
4結(jié)論
壽光雞GPR147基因5′調(diào)控區(qū)-593和-579位點的多態(tài)性對產(chǎn)蛋數(shù)有影響,單倍型TA個體擁有相對較高的中期產(chǎn)蛋數(shù)、后期產(chǎn)蛋數(shù)、300天產(chǎn)蛋量以及56周總產(chǎn)蛋量。因此,單倍型TA是具有潛在應用價值的分子遺傳標記,可為壽光雞產(chǎn)蛋性狀的標記輔助選擇提供依據(jù)。
參考文獻(References):
[1]TSUTSUI K,SAIGOH E,UKENA K,et al.A novel avian hypothalamic peptide inhibiting gonadotropin release[J].BiochemBiophysResCommun,2000,275(2):661-667.
[2]OSUGI T,UKENA K,BENTLEY G E,et al.Gonadotropin-inhibitory hormone in Gambel’s white-crowned sparrow (Zonotrichia leucophrys gambelii):cDNA identification,transcript localization and functional effects in laboratory and field experiments[J].JEndocrinol,2004,182(1):33-42.
[3]UKENA K,UBUKA T,TSUTSUI K.Distribution of a novel avian gonadotropin-inhibitory hormone in the quail brain[J].CellTissueRes,2003,312(1):73-79.
[4]UBUKA T,KIM S,HUANG Y C,et al.Gonadotropin-inhibitory hormone neurons interact directly with gonadotropin-releasing hormone-I and-Ⅱ neurons in European starling brain[J].Endocrinology,2008,149(1):268-278.
[5]KRIEGSFELD L J,MEI D F,BENTLEY G E,et al.Identification and characterization of a gonadotropin-inhibitory system in the brains of mammals[J].ProcNatlAcadSciUSA,2006,103(7):2410-2415.
[6]BENTLEY G E,TSUTSUI K,KRIEGSFELD L J.Recent studies of gonadotropin-inhibitory hormone (GnIH) in the mammalian hypothalamus,pituitary and gonads[J].BrainRes,2010,1364:62-71.
[7]CALISI R M,RIZZO N O,BENTLEY G E.Seasonal differences in hypothalamic EGR-1 and GnIH expression following capture-handling stress in house sparrows (Passer domesticus)[J].GenCompEndocrinol,2008,157(3):283-287.
[8]IKEMOTO T,PARK M K.Chicken RFamide-related peptide (GnIH) and two distinct receptor subtypes:identification,molecular characterization,and evolutionary considerations[J].JReprodDev,2005,51(3):359-377.
[9]BONINI J A,JONES K A,ADHAM N,et al.Identification and characterization of two G protein-coupled receptors for neuropeptide FF[J].JBiolChem,2000,275(50):39324-39331.
[10]HINUMA S,SHINTANI Y,F(xiàn)UKUSUMI S,et al.New neuropeptides containing carboxy-terminal RFamide and their receptor in mammals[J].NatCellBiol,2000,2(10):703-708.
[11]UBUKA T,UKENA K,SHARP P J,et al.Gonadotropin-inhibitory hormone inhibits gonadal development and maintenance by decreasing gonadotropin synthesis and release in male quail[J].Endocrinology,2006,147(3):1187-1194.
[12]趙雪丹,王慧,康麗,等.褪黑激素受體MTNR1B 5′調(diào)控區(qū)多態(tài)及其與雞產(chǎn)蛋性狀的關(guān)聯(lián)分析[J].中國獸醫(yī)學報,2014(11):1855-1860.
ZHAO X D,WANG H,KANG L,et al.Polymorphisms of the 5′ regulatory region of melatonin receptor gene MTNR1B and their correlations with laying traits in chickens[J].ChineseJournalofVeterinaryScience,2014(11):1855-1860.(in Chinese)
[13]宋倩,王慧,康麗,等.壽光雞褪黑素受體基因(MTNR1A)5′控區(qū)多態(tài)性與產(chǎn)蛋性狀的關(guān)聯(lián)分析[J].農(nóng)業(yè)生物技術(shù)學報,2013,21(5):579-587.
SONG Q,WANG H,KANG L,et al.Polymorphisms of the 5′ regulatory region of melatonin receptor gene (MTNR1A) and their associations with laying traits in Shouguang chickens(Gallusgallus)[J].JournalofAgriculturalBiotechnology,2013,21(5):579-587.(in Chinese)
[14]賈祥捷,王長法,楊桂文,等.中國荷斯坦牛POU1F1基因與PRL基因的多態(tài)性及其聚合效應對產(chǎn)奶性狀的影響[J].遺傳,2011,33(12):1359-1365.
JIA X J,WANG C F,YANG G W,et al.Polymorphism of POU1F1 gene and PRL gene and their combined effects on milk performance traits in Chinese Holstein cattle[J].Hereditas(Beijing),2011,33(12):1359-1365.(in Chinese)
[15]葉凱.尼羅羅非魚GATA因子表達及其功能的初步研究[D].重慶:西南大學,2012.YE K.Preliminary studies on expression patterns and functional analysis of GATA factors in Nile tilapia[D].Chongqing:Southwest University,2012.(in Chinese)[16]DASEN J S,O’CONNELL S M,F(xiàn)LYNN S E,et al.Reciprocal interactions of Pit1 and GATA2 mediate signaling gradient-induced determination of pituitary cell types[J].Cell,1999,97(5):587-598.
[17]GORDON D F,WOODMANSEE W W,BLACK J N,et al.Domains of Pit-1 required for transcriptional synergy with GATA-2 on the TSH beta gene[J].MolCellEndocrinol,2002,196(1-2):53-66.
[18]ROSENFELD M G,BRIATA P,DASEN J,et al.Multistep signaling and transcriptional requirements for pituitary organogenesisinvivo[J].RecentProgHormRes,2000,55:1-14.
[19]張未麗,劉智皓,吳風瑞,等.轉(zhuǎn)錄因子GATA-2在脊椎動物發(fā)育過程中的作用[J].四川動物,2009,28(6):945-948.
ZHANG W L,LIU Z H,WU F R,et al.Role of the GATA-2 transcription factor in vertebrates during development[J].SichuanJournalofZoology,2009,28(6):945-948.(in Chinese)
(編輯郭云雁)
Association of Polymorphisms in the 5′ Regulatory Region ofGPR147 with Laying Performance in Shouguang Chickens
JIA Mei-ting,LI Xian-yao,KANG Li,JIANG Yun-liang,F(xiàn)AN Xin-zhong,TANG Hui*
(CollegeofAnimalScienceandVeterinaryMedicine,ShandongAgriculturalUniversity,Tai’an271018,China)
Key words:chicken;laying traits;GPR147 gene;polymorphism
Abstract:To investigate the single nucleotide polymorphisms (SNPs) ofGPR147 gene and the correlations between SNPs and laying traits,the SNPs in the 5′ regulatory region ofGPR147 of Shouguang chicken were detected by direct sequencing method and their associations with laying traits of Shouguang chickens were analyzed.The results showed that there were 11 SNPs loci in the 5′ regulatory region ofGPR147 gene detected in the population.The -593 and -579 loci had significant effects on egg production from 45- to 56-week-old and total egg number until 56-week-old (P<0.05).Genotype TT at the -593 locus and AA at -579 locus were associated with higher egg production.Haplotype analysis of the 2 loci showed that hens with haplotype TA laid 1.5,2.5,1.6 and 4.4 eggs more than that with haplotype CG during 26-44 weeks of age (P=0.041 8),45-56 weeks of age (P=0.032 6),the period until 300 days of age (P=0.1) and the total period until 56 weeks of age (P=0.014 4),respectively.In summary,the 5′ regulatory region of theGPR147 has highly polymorphic,SNPs and haplotypes at -593 and -579 loci have significant impacts on the egg number of Shouguang chickens,and they can be used as potential and valuable candidate molecular markers.
doi:10.11843/j.issn.0366-6964.2016.01.005
收稿日期:2015-01-21
基金項目:山東省現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術(shù)體系項目(SDAIT-13-011-02)
作者簡介:賈美婷(1989-),女,山東濟寧人,碩士生,主要從事動物遺傳與育種研究,E-mail:jiameiting2008@163.com *通信作者:唐輝,副教授,碩士生導師,E-mail:tanghui@sdau.edu.cn
中圖分類號:S813.2
文獻標志碼:A
文章編號:0366-6964(2016)01-0034-07