ESR基因多態性與湖羊產羔性能的關系

郭海燕+金銀+李擁軍+張昊+程國虎

doi：10.15889/j.issn.1002-1302.2017.13.035[HT9.]

摘要：利用PCR-單鏈構象多態性（PCR-SSCP）技術及基因測序技術對湖羊進行了ESR基因的多態性檢測及其與湖羊產羔數的統計分析。測序結果表明，等位基因G與C相比發生了1處堿基突變（G363C）。ESR基因在湖羊中存在多態性（CC、CG和GG）；湖羊中C等位基因的頻率為0.743，G等位基因的頻率為0.257。CC、CG、GG基因型頻率分別為0.564、0.357、0.079。湖羊CG、GG、CC基因型的平均產羔數分別為1.98、1.80、1.96只。數據分析結果表明，ESR基因的多態性對湖羊產羔數的影響差異不顯著（P≥0.05）。

關鍵詞：湖羊；ESR基因；產羔數

中圖分類號： S826.3文獻標志碼： A[HK]

文章編號：1002-1302（2017）13-0126-03[HS）][HT9.SS]

[HJ1.4mm]

收稿日期：2017-01-03

基金項目：江蘇省農業科技自主創新資金[編號：CX（15）1007]。

作者簡介：郭海燕（1992—），女，河南洛陽人，碩士研究生，研究方向為遺傳育種與繁殖。E-mail：1937750042@qq.com。

通信作者：李擁軍，博士，教授，研究方向為養羊生產與羊的繁育。Tel：（0514）87996481；E-mail：liyj@yzu.edu.cn。[HJ]

[ZK）]

湖羊是我國著名的地方多胎品種羊，也是我國一級保護地方畜禽品種，主要分布在江、浙一帶。湖羊1年2胎，每胎多羔，其產羔率能達到200%～250%[1]。綿羊的繁殖力是影響其經濟效益的重要指標之一，而湖羊又是我國著名的多胎品種，因此研究其多胎的分子機制對于湖羊多胎資源的開發利用、新品種的培育以及生產實踐都有著重要意義。

湖羊高產的繁殖性能是得到世界公認的，因此探究并確定影響湖羊多胎的主效基因一直是國內外的研究熱點。綿羊的多胎性狀是由多種基因控制的，且不同品種間控制多胎的主效基因也是不同的。有研究發現，品種之間繁殖性能的差異主要取決于產仔母羊的基因，與后代個體基因無太大關系，且推斷湖羊具有多胎的主基因，并呈顯性遺傳[2-3]。很多學者已證實了[WTBX][STBX]BMPR-IB[WTBZ][STBZ]可能是湖羊多胎性的主效基因，而[WTBX][STBX]BMP15、GDF9[WTBZ][STBZ]等基因可能也參與并影響了湖羊的多胎機制[4-6]。有研究證明湖羊存在FecB基因，且湖羊全部為FecB基因純合攜帶者[7-8]。

雌激素受體（estrogen receptor，簡稱ESR）是核受體超家族的一員，是一種配體依賴性的轉錄因子[9]。ESR能夠轉錄調控蛋白質，參與雌性動物性腺基因的表達與調控。ESR結合配體之后與雌激素應答元件相互作用，可以改變受雌激素調控基因的轉錄，從而影響胚胎發育和系統分化[10]。

許多研究都表明，ESR基因在動物繁殖過程中發揮著重要作用，而且已經確定ESR基因是影響豬產仔數的主效基因[11-17]。研究發現，將小鼠的ESR基因敲除后，小鼠不能排卵，這也證明了ESR基因與繁殖性能的重要聯系[18]。張小雪等對綿羊的ESR基因進行了詳細的生物信息學分析，發現綿羊的ESR基因主要位于細胞核中，并且參與機體的轉錄調控過程[19]。有研究認為ESR可能是控制湖羊和小尾寒羊多胎性能的主效基因或與之存在緊密的連鎖[9，20]。賈立華等對小尾寒羊的ESR基因外顯子4序列進行克隆與序列分析，發現此段序列有很強的保守性[21-22]。馮濤等也得出了相同結論，發現哺乳動物ESR基因外顯子1、4序列的保守性較強，因此此段區域可能不是影響山羊高繁殖力的功能結構域[23]。還有研究發現，ESR基因第5外顯子Pvu Ⅱ多態性與動物的繁殖性能密切相關，但該位點在不同群體中的效應不一致[24-27]。此外，ESR基因對胚胎發育及黃體期的影響及作用也一直是國內外研究熱點[28-30]。本試驗通過研究湖羊ESR基因的多態性及其與產羔數的關系，旨在為篩選湖羊多胎性候選基因提供理論依據。

1材料與方法

1.1試驗動物及血樣采集

選取具有第1胎和第2胎生產記錄的湖羊，共342只。采集湖羊頸靜脈血液約3 mL，并記錄耳標號，羊全部來自江蘇省太倉金倉湖羊場。采血用一次性采血盛血器，并且提前在采血管中加入ACD抗凝劑（主要含枸椽酸、枸椽酸三鈉、葡萄糖）。血樣用冰盒帶回實驗室后于4 ℃保存，以供提取DNA。

1.2DNA的提取及檢測

DNA提取方法參照天根生化科技（北京）有限公司的血液基因組DNA提取試劑盒說明書，提取后于-20 ℃保存。

1.3PCR的擴增

ESR的引物參照畢曉丹的方法[9]，由鉑尚生物技術（上海）有限公司合成。引物序列如下：

上游引物：5′-TGCACCAGATCCAAGCCAACGA-3′；

下游引物：5′-CGGGTACCTGTAGAAGGCGGGAG-3′。

PCR反應體系25 μL，其中DNA 1 μg，上下游引物（10 μmol/L）各1 μL，2×Taq PCR MasterMix[天根生化科技（北京）有限公司]12.5 μL，ddH2O補至25 μL。PCR反應程序：95 ℃預變性4 min；94 ℃變性1 min，57 ℃復性90 s，72 ℃延伸1 min，34個循環；72 ℃終延伸7 min。產物用2%瓊脂糖凝膠檢測。

[HTK]1.4單鏈構象多態性（SSCP）檢測多態性[HT]endprint

取5 μL PCR產物，并與5 μL上樣緩沖液（90%去離子甲酰胺9 mL，0.25%溴酚藍，0.25%二甲苯青FF，1×TBE 1 mL）混合后置于PCR管內，98 ℃變性10 min，迅速放入 -20 ℃ 冰箱10 min，使之保持變性狀態。然后用8%非變性聚丙烯酰胺凝膠電泳檢測（120 V，3 h），電泳結束銀染顯色，拍照分析。

1.5公司測序

將提取的DNA全部送到無錫翼和應用生物技術有限公司進行測序及基因型的鑒定，并整理數據結果。

1.6統計分析

收集整理湖羊的生產記錄，計算湖羊ESR不同基因型及基因的頻率。用SPSS 22.0處理數據，對不同ESR基因型的湖羊產羔數進行LSD多重比較，從而分析比較不同基因型與產羔數之間的關系。

2結果與分析

2.1PCR擴增結果

對PCR產物進行瓊脂糖電泳檢測，結果表明其特異性良好（圖1），目的條帶長度符合預期，為419 bp左右，可直接進行SSCP分析。

2.2SSCP檢測結果

將PCR產物進行SSCP分析，結果表明ESR基因存在多態性，擴增片段存在3種基因型，將其定義為CC、CG和GG（圖2）。其中第10道為CG，第5道為GG，其余幾道為CC。

2.3不同基因型測序結果

為確定ESR基因的突變位點，將CC、CG和GG基因型的PCR產物送無錫翼和應用生物技術有限公司測序。根據測序結果可知，等位基因G與C相比，發生了1處突變（G363C），測序結果見圖3。[FL）]

[FL（2K2]2.4湖羊ESR基因型及基因頻率

檢測結果表明，在342只湖羊樣本中，193只是CC型，122只是CG型，27只是GG型。經計算，C等位基因的頻率為0.743，G等位基因的頻率為0.257。CC基因型頻率為 0.564，CG基因型頻率為0.357，GG基因型頻率為0.079。

2.5湖羊ESR不同基因型的產羔數的均值及標準誤

由表1可知，根據第1胎產羔記錄，GG與CC、CG產羔數差異顯著，CC、CG之間差異不顯著。根據第2胎湖羊產羔記錄，CC、CG、GG基因型的產羔數差異不顯著。從整體水平看，CC、CG、GG基因型的產羔數差異不顯著。[FL）]

3討論與結論

畢曉丹通過研究高繁殖力綿羊品種與低繁殖力綿羊品種的ESR基因的多態性，發現ESR基因外顯子1的第363處發生了C-G突變，在小尾寒羊中AB、BB基因型比AA基因型多產0.51、0.70只羔羊，并認為ESR基因可能是控制小尾寒羊多胎的主效基因或與之存在緊密的連鎖[9]。董文艷等對湖羊ESR基因第一外顯子部分序列進行多態性檢測[20]，結果與畢曉丹的研究結果[9]相似，ESR基因也發生了C-G的突變，并認為ESR基因可能是控制湖羊多羔性能的1個主效基因或與之存在緊密遺傳連鎖的1個標記。李廣錄以中國美利奴羊、湖羊、羅米麗羊和羅米麗×中國美利奴（新疆軍墾型）為研究對象，檢測其ESR基因的5′非翻譯區片段的多態性，結果顯示4種綿羊均出現多態性[10]。狄冉等采用PCR-[JP2]SSCP技術分析ESR基因外顯子4在高繁殖力綿羊品種（小尾寒羊和湖羊）和低繁殖力綿羊品種（特克塞爾、中國美利奴、考力代和杜泊）中的單核苷酸多態性，結果顯示ESR基因在這6個綿羊品種中均不存在多態性，說明所檢測的ESR基因外顯子4序列可能不是影響綿羊高繁殖力的功能結構域[22]。[JP]

本試驗使用的所有湖羊均來自江蘇省太倉金倉湖羊場，此羊場是新建羊場，羊群剛購進不久，且未經過選育，羊場只有湖羊第1、2胎的生產記錄。經統計分析，本次試驗羊第1胎平均產羔1.86只，第2胎試驗羊平均產羔2.04只，此羊場湖羊的整體生產水平未達到通常文獻中報道的每胎 2.30 只羔羊。因此該羊場湖羊的群體構成及其胎次因素可能會對本試驗結果造成一定影響。

本試驗結果顯示，檢測的342只湖羊中ESR基因存在多態性，分別為CG、CC和GG。經檢測發現ESR基因外顯子1的第363處發生了C-G突變。根據湖羊產羔數的統計分析結果可知，CG、CC、GG之間的產羔數并無顯著差異。而這一結果與畢曉丹等的研究結果[9，20]相悖，他們的研究結果表明，等位基因B（B等同于本試驗中的等位基因G）與綿羊的產羔數呈正相關。筆者認為造成這一現象的原因，與上述提及的該羊場的羊群結構以及胎次因素可能有關。湖羊是世界著名的多胎品種，其多胎性能是由多個主效基因控制的，本試驗只是單一地分析了ESR基因與湖羊多胎性能的關系。因此，要確定ESR基因與湖羊繁殖性能的確切關系，還應進行多基因的聯合分析，以尋找各個候選基因之間的聯系。就本研究而言，筆者認為可以在后續試驗中進一步獲取該羊場更多的胎次數據以及擴大樣本數量，從而進一步驗證ESR基因與湖羊產羔數的關系。

[JP3]本試驗表明，ESR基因多態性對湖羊產羔性能無顯著影響。[JP]

[HS2]參考文獻：

[1][ZK（#]史洪才，武堅. 我國對湖羊和小尾寒羊多胎基因的研究[J]. 草食家畜，2004（4）：14-16.

[2]姜運良，吳常信. Booroola Merino綿羊多胎基因FecB的研究進展[J]. 中國畜牧雜志，1999（4）：51-53.

[3]石國慶. 湖羊多胎機制研究[D]. 南京：南京農業大學，2006.

[4]吳麗麗，趙曉楓，張立凡，等. 湖羊[WTBX][STBX]BMPR-IB、BMP15和GDF9[WTBZ][STBZ]基因的RFLP分析[J]. 中國草食動物，2011，31（4）：9-12.

[5]王啟貴，鐘發剛，李輝，等. 綿羊BMPR-IB基因多態性與其產羔數的相關研究[J]. 草食家畜，2003（2）：20-23.endprint

[6]孫紅霞，田秀娥，王永軍. [WTBX][STBX]BMPR-IB，BMP15和GDF9[WTBZ][STBZ]基因作為灘羊繁殖性狀主效候選基因的研究[J]. 西北農業學報，2009，18（5）：17-21.

[7]王根林，毛鑫智，George H D，等. DNA分析發現我國湖羊和小尾寒羊存在Booroola（FecB）多胎基因[J]. 南京農業大學學報，2003，26（1）：104-106.

[8]管峰，艾君濤，劉守仁，等. [WTBX][STBX]BMPR-IB和BMP15[WTBZ][STBZ]基因作為湖羊多胎性候選基因的研究[J]. 家畜生態學報，2005，26（3）：9-12.

[9]畢曉丹. 小尾寒羊高繁殖力候選基因ESR的研究[D]. 延吉：延邊大學，2005.[ZK）]

[10][ZK（#]李廣錄. 綿羊繁殖性狀相關的幾個功能基因的多態性分析[D]. 石河子：石河子大學，2008.

[11]胡建偉，楊曉瓊，余永華. 中外豬種ESR基因多態性與繁殖性狀相關研究[J]. 上海畜牧獸醫通訊，2014（2）：20-23.

[12]武建亮，徐利，周海深，等. 豬繁殖性狀主效基因及其合并基因型對產活仔數的影響[J]. 豬業科學，2016（2）：118-121.[JP]

[13]管玲玉，馬力鵬，劉麗娟，等. FSH-β和ESR基因多態性對大白豬產活仔數的影響[J]. 豬業科學，2015（1）：114-115.

[14][JP2]李俊，陳美娟，許鐘峯，等. 環江香豬FSH-β，ESR及[WTBX][STBX]ZAR1[WTBZ][STBZ]基因多態性分析[J]. 基因組學與應用生物學，2015（3）：495-499.[JP]

[15]Rothschild M，Jacobson C，Vaske D，et al. The estrogen receptor locus is associated with a major gene influencing litter size in pigs[J]. Proceedings of the National Academy of Sciences of the United States of America，1996，93（1）：201-205.

[16]Green S，Walter P，Kumar V，et al. Human oestrogen receptor cDNA：sequence，expression and homology to v-erb-A[J]. Nature，1986，320（658）：134-139.

[17]丁家桐，葛紅山，姜勛平，等. 豬雌激素受體基因與產仔數和乳頭數的關系研究[J]. 南京農業大學學報，2002，25（2）：79-81.

[18]Hewitt S C，Korach K S. Oestrogen receptor knockout mice：roles for oestrogen receptors alpha and beta in reproductive tissues[J]. Reproduction，2003，125（2）：143-149.

[19]張小雪，潘香羽，李發弟，等. 綿羊ESR基因生物信息學分析[J]. 甘肅農業科技，2014（9）：30-33.

[20]董文艷，陳阿琴，王爭光，等. 湖羊高繁殖力候選基因ESR的研究[J]. 浙江農業學報，2009，21（6）：561-564.

[21]賈立華，儲明星，陳宏權，等. 小尾寒羊雌激素受體基因外顯子4的克隆與序列分析[J]. 中國畜牧獸醫，2008，35（2）：60-63.

[22]狄冉，賈立華，儲明星，等. 綿羊雌激素受體基因外顯子4多態性分析[J]. 中國畜牧獸醫，2008，35（12）：89-92.

[23][JP3]馮濤，趙有璋，狄冉，等. 山羊雌激素受體（ESR）基因部分外顯[JP2]子多態性分析[J]. 農業生物技術學報，2009，17（2）：237-242.[JP]

[24]張淑君，熊遠，鄧昌彥，等. 豬產仔數分子標記及其效應分析[J]. 生物技術通報，2003（2）：36-39.

[25]田勇，李學偉，李芳瓊，等. 四川省外種豬ESR基因對繁殖及生長性狀的影響[J]. 黑龍江畜牧獸醫，2004（10）：12-14.

[26][JP3]彭淑紅，徐寧迎，孫萬元，等. 金華豬雌激素受體基因的分布及其對頭胎產仔數的影響[J]. 中國畜牧雜志，2004，40（3）：30-31.[JP]

[27]孟慶利，劉鐵錚. 豬雌激素受體基因PvuⅡ酶切片段多態性與產仔性能的關系[J]. 江蘇農業學報，2005，21（1）：49-52.

[28]Whittle W L，Holloway A C，Lye S，et al. The pattern of glucocorticoid and estrogen receptors may explain differences in steroid dependency of intrauterine prostaglandin production at parturition in sheep[J]. Journal of the Society for Gynecologic Investigation，2006，13（7）：506-511.

[29]Reynolds L P，Haring J S，Johnson M L，et al. Placental development during early pregnancy in sheep：estrogen and progesterone receptor messenger RNA expression in pregnancies derived from in vivo-produced and in vitro-produced embryos[J]. Domestic Animal Endocrinology，2015，53：60-69.

[30]Tasende C，Rodríguez-Pión M，Acua S，et al. Corpus luteum life span and pituitary oestrogen and progesterone receptors in cyclic and gonadotrophin-releasing hormone-treated anoestrous ewes[J]. Reproduction Fertility and Development，2005，17（7）：721-726.endprint