陳來(lái)運(yùn) 袁超 孫曉萍

摘要：旨在探究4個(gè)綿羊品種的多胎性狀分子遺傳機(jī)制，以高山美利奴羊、青海細(xì)毛羊、中國(guó)美利奴羊和敖漢細(xì)毛羊?yàn)檠芯繉?duì)象，通過(guò)分析全基因組測(cè)序結(jié)果對(duì)4個(gè)綿羊品種FSHR基因外顯子單核苷酸多態(tài)性進(jìn)行篩選，利用專門軟件預(yù)測(cè)單核苷酸多態(tài)性對(duì)FSHR基因mRNA二級(jí)結(jié)構(gòu)、蛋白質(zhì)二級(jí)結(jié)構(gòu)和三級(jí)結(jié)構(gòu)的影響。研究發(fā)現(xiàn)，高山美利奴羊、中國(guó)美利奴羊、敖漢細(xì)毛羊FSHR基因在外顯子上分別存在3個(gè)SNPs（T904G、C975T、G1572A）、7個(gè)SNPs（G149A、G813A、T817C、T904G、C975T、T1234C、G1572A）、6個(gè)SNPs（G813A、T817C、T904G、C975T、T1234C、G1572A），青海細(xì)毛羊FSHR基因在外顯子上不存在突變位點(diǎn)。生物信息學(xué)分析發(fā)現(xiàn)，試驗(yàn)所檢測(cè)到的7個(gè)SNPs中T817C、T904G、T1234C導(dǎo)致mRNA的二級(jí)結(jié)構(gòu)和最小自由能發(fā)生改變，G813A引起最小自由能改變，而未引起mRNA二級(jí)結(jié)構(gòu)發(fā)生改變，G149A、C975T與G1572A未引起最小自由能與mRNA的二級(jí)結(jié)構(gòu)發(fā)生改變。5個(gè)錯(cuò)義突變位點(diǎn)（G149A、G813A、T904G、C975T、G1572A）均導(dǎo)致編碼蛋白質(zhì)二級(jí)結(jié)構(gòu)與三級(jí)結(jié)構(gòu)發(fā)生改變。

關(guān)鍵詞：綿羊;FSHR基因;多態(tài)性;SNPs;生物信息學(xué)分析

中圖分類號(hào)： S826.2 ?文獻(xiàn)標(biāo)志碼： A ?文章編號(hào)：1002-1302（2019）22-0047-05

卵泡刺激素受體（follicle-stimulating hormone receptor，F(xiàn)SHR）屬于G蛋白偶聯(lián)受體超家族中的糖蛋白亞家族成員，成熟的FSHR是由二硫鍵構(gòu)成的穩(wěn)定低聚體[1]。FSHR與FSH（follicle stimdating hormone，F(xiàn)SH）功能的發(fā)揮有著緊密的聯(lián)系[2]。FSH是由垂體合成并分泌的一種糖蛋白類促性腺激素，在卵泡生長(zhǎng)過(guò)程中發(fā)揮關(guān)鍵作用，但其是一種不能穿過(guò)細(xì)胞膜的生物大分子，必須通過(guò)FSHR特異性介導(dǎo)[3]。因此，F(xiàn)SHR對(duì)動(dòng)物卵泡的發(fā)育與成熟具有重要作用，從而對(duì)動(dòng)物繁殖性狀造成影響。研究證實(shí)，F(xiàn)SHR基因顯著影響動(dòng)物繁殖性狀，吳井生等對(duì)106頭小梅山豬研究發(fā)現(xiàn)，F(xiàn)SHR基因第10個(gè)外顯子存在2個(gè)多態(tài)位點(diǎn)對(duì)小梅山豬產(chǎn)仔數(shù)有顯著影響，在P1位點(diǎn)上，2胎以上小梅山母豬中AA型個(gè)體的總產(chǎn)仔數(shù)、產(chǎn)活仔數(shù)比BB型分別顯著高出1.95、1.66頭（P<001）;在P3位點(diǎn)上，2胎以上及所有胎次的母豬中，DD型個(gè)體的總產(chǎn)仔數(shù)和產(chǎn)活仔數(shù)均極顯著高于CC型（P<0.01）[4]。邢耀潭對(duì)12頭單胎、92頭雙胎天祝白牦牛研究發(fā)現(xiàn)，F(xiàn)SHR基因 P3片段（-1 195 bp處）在天祝白牦牛單胎群體和雙胎群體中，雙胎母牛的突變率顯著高于單胎母牛的突變率，說(shuō)明FSHR基因很可能是控制天祝白牦牛雙胎性狀的主效基因[5]。Chu等統(tǒng)計(jì)252只小尾寒羊不同F(xiàn)SHR基因型產(chǎn)羔數(shù)的最小二乘平均值和標(biāo)準(zhǔn)差，發(fā)現(xiàn)FSHR基因不同基因型產(chǎn)羔數(shù)存在極顯著差異[6]。

Jiang等于2014年公布綿羊新一代從頭測(cè)序的基因組信息，綿羊FSHR基因（GenBank登錄號(hào)為NC_019460）位于第3號(hào)染色體上，基因全長(zhǎng)196 149 bp，包含9個(gè)內(nèi)含子和10個(gè)外顯子，外顯子總長(zhǎng)為2 431 bp，編碼395個(gè)氨基酸[7]。目前，對(duì)動(dòng)物FSHR基因SNP的研究主要集中于豬[4，8]、牛[5，9]、羊[6，10]等的第10個(gè)外顯子及5′ 側(cè)翼區(qū)，但對(duì)FSHR基因其他外顯子分析較少。因此，本研究選取高山美利奴羊、青海細(xì)毛羊、中國(guó)美利奴羊和敖漢細(xì)毛羊?yàn)檠芯繉?duì)象，根據(jù)全基因組測(cè)序的結(jié)果，分析FSHR基因外顯子在4個(gè)綿羊品種上的單核苷酸多態(tài)性（single nucleotide polymorphism，SNP），并對(duì)篩選出的SNPs進(jìn)行生物信息學(xué)分析，以期篩選出與4個(gè)綿羊品種產(chǎn)羔數(shù)相關(guān)的遺傳標(biāo)記，為4個(gè)綿羊品種多胎性狀的研究奠定分子基礎(chǔ)。

1 材料與方法

1.1 材料

采集30只高山美利奴羊（甘肅省綿羊繁殖技術(shù)推廣站）、30只中國(guó)美利奴羊（新疆鞏乃斯種羊場(chǎng)）、30只敖漢細(xì)毛羊（內(nèi)蒙古敖漢種羊場(chǎng)）、30只青海細(xì)毛羊（青海省三角城種羊場(chǎng)）共120只羊血樣，每只羊頸部靜脈采血5 mL于江蘇宇力醫(yī)療器械有限公司EDTA-K2抗凝真空采血管中，每只綿羊分別采集2管血液樣品，于-20 ℃冰箱冷凍保存，用于DNA提取。

1.2 方法

1.2.1 全基因組測(cè)序 4個(gè)綿羊品種共120只，分別對(duì)每只羊取1管血液并送至北京諾禾致源科技股份有限公司，提取DNA樣品（試驗(yàn)過(guò)程中須無(wú)降解、無(wú)污染、無(wú)斷裂等，以保證所提4個(gè)綿羊品種DNA純度及完整性），經(jīng)檢測(cè)合格的4個(gè)綿羊品種的DNA樣品分別每個(gè)樣品取2 μg用于10×genomics平臺(tái)建庫(kù)，基于Illumina Hiseq平臺(tái)進(jìn)行雙末端 150 bp 測(cè)序。根據(jù)2014年公布的綿羊參考基因組Oar_v 40（GCF_000298735.2），利用不同綿羊品種的測(cè)序結(jié)果序列reads，應(yīng)用SOAPdenovo（諾禾自主研發(fā)軟件）進(jìn)行組裝得到長(zhǎng)的序列片段，然后利用lastz軟件獲得4個(gè)綿羊品種的每個(gè)個(gè)體全基因組信息，該試驗(yàn)由北京諾禾致源科技股份有限公司完成。

1.2.2 FSHR基因多態(tài)位點(diǎn)篩選 將120個(gè)未分群未過(guò)濾的SNP vcf文件按照4個(gè)品種（高山美利奴羊、青海細(xì)毛羊、中國(guó)美利奴羊、敖漢細(xì)毛羊）進(jìn)行分群，所用軟件為vcftools_v 0114。將分好群的SNP文件進(jìn)行過(guò)濾，過(guò)濾條件為：?jiǎn)蝹€(gè)樣本測(cè)序覆蓋深度>2、maf（最小等位基因頻率）>0.05、mis（缺失）<0.1。根據(jù)FSHR基因位置信息：Chromosome 2，NC_000002.12 （48953161-49154527，complement），利用 vcftools_v 0.1.14軟件進(jìn)行4個(gè)綿羊品種FSHR基因外顯子SNP位點(diǎn)的提取，然后統(tǒng)計(jì)各個(gè)品種FSHR基因外顯子SNP位點(diǎn)數(shù)，利用A-NNOVAR軟件將各個(gè)品種FSHR基因外顯子SNP位點(diǎn)進(jìn)行注釋。

1.2.3 FSHR基因多態(tài)位點(diǎn)遺傳多樣性分析 通過(guò)生物學(xué)分析軟件POPGEN 1.32軟件統(tǒng)計(jì)4個(gè)綿羊品種FSHR基因各SNPs位點(diǎn)的野生型、雜合型及突變純合型的基因型頻率、野生型和突變型等位基因頻率、Hardy-Weinberg平衡檢驗(yàn)，利用群體多態(tài)信息含量（PIC）計(jì)算程序計(jì)算4個(gè)綿羊品種FSHR基因各SNPs位點(diǎn)遺傳多樣性參數(shù)多態(tài)信息含量（PIC）。

1.2.4 生物信息學(xué)分析 4個(gè)綿羊品種FSHR基因由于SNPs產(chǎn)生了不同的基因型。利用RNA二級(jí)結(jié)構(gòu)分析軟件RNAfold進(jìn)行4個(gè)綿羊品種不同基因型的FSHR基因RNA二級(jí)結(jié)構(gòu)預(yù)測(cè);利用NPSA軟件中的MLRC程序進(jìn)行4個(gè)綿羊品種不同基因型的FSHR基因蛋白質(zhì)二級(jí)結(jié)構(gòu)預(yù)測(cè);利用蛋白質(zhì)三級(jí)結(jié)構(gòu)預(yù)測(cè)軟件SWISS-MODEL進(jìn)行4個(gè)綿羊品種不同基因型的FSHR基因蛋白質(zhì)三級(jí)結(jié)構(gòu)預(yù)測(cè)。

2 結(jié)果與分析

2.1 4個(gè)綿羊品種FSHR基因外顯子SNPs篩選

根據(jù)4個(gè)綿羊品種共120只羊的全基因組測(cè)序結(jié)果與綿羊參考基因組比對(duì)分析4個(gè)綿羊品種FSHR基因外顯子上存在的多態(tài)位點(diǎn)，高山美利奴羊FSHR基因外顯子上共發(fā)現(xiàn)3個(gè)SNPs，包括2個(gè)轉(zhuǎn)換位點(diǎn)（C→T、G→A）和1個(gè)顛換位點(diǎn)（T→G）;中國(guó)美利奴羊FSHR基因外顯子上共發(fā)現(xiàn)7個(gè)SNPs，包括6個(gè)轉(zhuǎn)換位點(diǎn)（3個(gè)G→A、2個(gè)T→C和1個(gè)C→T）和1個(gè)顛換位點(diǎn)（T→G）;敖漢細(xì)毛羊INHA外顯子上FSHR基因外顯子上共發(fā)現(xiàn)6個(gè)SNPs，包括5個(gè)轉(zhuǎn)換位點(diǎn)（2個(gè)G→A、2個(gè)T→C和1個(gè)C→T）和1個(gè)顛換位點(diǎn)（T→G）;青海細(xì)毛羊FSHR基因外顯子上沒(méi)有發(fā)現(xiàn)突變。通過(guò)對(duì)4個(gè)綿羊品種FSHR基因外顯子氨基酸序列進(jìn)行分析發(fā)現(xiàn)，高山美利奴羊發(fā)生了3個(gè)錯(cuò)義突變，中國(guó)美利奴羊存在2個(gè)同義突變T817C、T1234C和5個(gè)錯(cuò)義突變G149A、G813A、T904G、C975T、G1572A，敖漢細(xì)毛羊存在2個(gè)同義突變T817C、T1234C和5個(gè)錯(cuò)義突變G149A、G813A、T904G、C975T、G1572A，4個(gè)綿羊品種FSHR基因外顯子均不存在無(wú)義突變（表1）。

2.2 4個(gè)綿羊品種多態(tài)位點(diǎn)統(tǒng)計(jì)學(xué)分析

根據(jù)4個(gè)品種共120只羊的全基因組測(cè)序結(jié)果分析高山美利奴羊FSHR基因外顯子SNPs位點(diǎn)基因型頻率、基因頻率、Hardy-Weinberg平衡狀態(tài)檢測(cè)（χ2、P值）、群體多態(tài)信息含量（PIC）計(jì)算結(jié)果。高山美利奴羊FSHR基因外顯子上檢測(cè)到的T904G、C975T、G1572A均表現(xiàn)出野生型和突變雜合型2種基因型，野生型基因型頻率均高于突變雜合型基因型頻率，野生型等位基因T、C、G分別為優(yōu)勢(shì)等位基因，均表現(xiàn)為低度多態(tài)。中國(guó)美利奴羊FSHR基因外顯子上檢測(cè)到的G149A、T904G、C975T、T1234C、G1572A表現(xiàn)出野生型和突變雜合型2種基因型，野生型基因型頻率均高于突變雜合型基因型頻率，野生型等位基因G、T、C、T、G分別為優(yōu)勢(shì)等位基因，均表現(xiàn)為低度多態(tài);中國(guó)美利奴羊FSHR基因外顯子上檢測(cè)到的G813A、T817C表現(xiàn)出野生型、突變雜合型和突變純合型3種基因型，G813A野生型基因型頻率高于突變雜合型與突變純合型基因型頻率和等位基因G為優(yōu)勢(shì)等位基因，T817C突變雜合型基因型頻率高于野生型與突變雜合型基因型頻率和等位基因C為優(yōu)勢(shì)等位基因，G813A、T817C都表現(xiàn)為中度多態(tài)。敖漢細(xì)毛羊FSHR基因外顯子上檢測(cè)到的T904G、C975T、G1572A表現(xiàn)出野生型和突變雜合型2種基因型，野生型基因型頻率均高于突變雜合型基因型頻率，野生型等位基因T、C、G分別為優(yōu)勢(shì)等位基因，均表現(xiàn)為低度多態(tài)。敖漢細(xì)毛羊FSHR基因外顯子上檢測(cè)到的G813A、T817C、T1234C表現(xiàn)出野生型、突變雜合型和突變純合型3種基因型，G813A和T1234C野生型基因型頻率高于突變雜合型與突變純合型基因型頻率和等位基因G、T分別為優(yōu)勢(shì)等位基因，T817C突變雜合型基因型頻率高于野生型與突變雜合型基因型頻率和等位基因C為優(yōu)勢(shì)等位基因，G813A、T817C和T1234C都表現(xiàn)為中度多態(tài);高山美利奴羊、中國(guó)美利奴羊和敖漢細(xì)毛羊FSHR基因外顯子上檢測(cè)到的所有突變位點(diǎn)χ2值均未達(dá)到顯著水平（P>0.05），說(shuō)明3個(gè)綿羊品種FSHR基因外顯子上檢測(cè)到的所有突變位點(diǎn)均達(dá)到Hardy-weinberg平衡狀態(tài)（表2）。

2.3 FSHR基因的RNA二級(jí)結(jié)構(gòu)分析

利用RNA二級(jí)結(jié)構(gòu)分析軟件RNAfold預(yù)測(cè)3個(gè)綿羊品種不同基因型的FSHR基因mRNA二級(jí)結(jié)構(gòu)，結(jié)果（圖1），表明G149A、C975T、G1572A位點(diǎn)的突變不會(huì)導(dǎo)致最小自由能的改變，G149A、C975T、G1572A位點(diǎn)的變異也不會(huì)導(dǎo)致mRNA二級(jí)結(jié)構(gòu)的改變;G813A位點(diǎn)的突變最小自由能增加，結(jié)構(gòu)穩(wěn)定性降低，G813A未導(dǎo)致mRNA二級(jí)結(jié)構(gòu)的改變;T817C、T904G和T1234C位點(diǎn)的突變均導(dǎo)致FSHR基因最小自由能改變，其中T817C導(dǎo)致FSHR基因RNA二級(jí)結(jié)構(gòu)最小自由能降低，二級(jí)結(jié)構(gòu)穩(wěn)定性增加，T904G和T1234C位點(diǎn)的突變導(dǎo)致FSHR基因RNA二級(jí)結(jié)構(gòu)最小自由能增加，二級(jí)結(jié)構(gòu)穩(wěn)定性降低，T817C、T904G和T1234C位點(diǎn)的突變導(dǎo)致FSHR基因mRNA二級(jí)結(jié)構(gòu)的改變。

2.4 FSHR基因多態(tài)位點(diǎn)蛋白質(zhì)二級(jí)結(jié)構(gòu)預(yù)測(cè)

利用NPSA軟件中的MLRC程序進(jìn)行4個(gè)綿羊品種不同基因型的FSHR基因蛋白質(zhì)二級(jí)結(jié)構(gòu)的預(yù)測(cè)，結(jié)果表明，F(xiàn)SHR基因編碼蛋白質(zhì)為全α蛋白，不存在β-轉(zhuǎn)角。錯(cuò)義突變G149A位點(diǎn)使編碼蛋白質(zhì)無(wú)規(guī)則卷曲所占比例減少，擴(kuò)展鏈所占比例增加，α-螺旋所占比例不變;錯(cuò)義突變G813A位點(diǎn)使編碼蛋白質(zhì)α-螺旋所占比例減少，無(wú)規(guī)則卷曲與擴(kuò)展鏈所占比例增加;錯(cuò)義突變位點(diǎn)T904G使編碼蛋白質(zhì)無(wú)規(guī)則卷曲與α-螺旋所占比例減少，擴(kuò)展鏈所占比例增加;錯(cuò)義突變位點(diǎn)C975T使編碼蛋白質(zhì)無(wú)規(guī)則卷曲所占比例減少，α-螺旋與擴(kuò)展鏈所占比例增加;錯(cuò)義突變位點(diǎn)G1572A使編碼蛋白質(zhì)無(wú)規(guī)則卷曲與擴(kuò)展鏈所占比例減少，α-螺旋所占比例增加。

2.5 FSHR基因多態(tài)位點(diǎn)蛋白質(zhì)三級(jí)結(jié)構(gòu)預(yù)測(cè)

利用蛋白質(zhì)三級(jí)結(jié)構(gòu)預(yù)測(cè)軟件SWISS-MODEL進(jìn)行3個(gè)綿羊品種不同基因型的FSHR基因蛋白質(zhì)三級(jí)結(jié)構(gòu)的預(yù)測(cè)，結(jié)果表明，3個(gè)綿羊品種各錯(cuò)義突變位點(diǎn)導(dǎo)致FSHR基因蛋白質(zhì)三級(jí)結(jié)構(gòu)發(fā)生明顯改變，與蛋白質(zhì)二級(jí)結(jié)構(gòu)預(yù)測(cè)結(jié)果相一致。

3 討論

羊的繁殖性狀是影響羊養(yǎng)殖業(yè)的重要經(jīng)濟(jì)性狀，直接影響綿羊養(yǎng)殖業(yè)的生產(chǎn)成本與生產(chǎn)效益，限制綿羊養(yǎng)殖業(yè)的發(fā)展，因此提高綿羊的繁殖性能對(duì)畜牧業(yè)的發(fā)展具有十分重要的意義[11]。大量研究表明，F(xiàn)SHR與哺乳動(dòng)物的繁殖性狀有密切的關(guān)系，例如雌性動(dòng)物卵巢功能[12-13]、卵泡生長(zhǎng)發(fā)育[14]和雌性動(dòng)物睪丸發(fā)育[15]、精子發(fā)生[16]、精子活力與精液質(zhì)量[17-19]等。近年，關(guān)于哺乳動(dòng)物FSHR基因多態(tài)性與多胎性狀之間的關(guān)系研究較多。許瑤等研究豬FSHR基因外顯子10的多態(tài)性與產(chǎn)仔數(shù)間的關(guān)聯(lián)分析發(fā)現(xiàn)，F(xiàn)SHR基因可作為提高豬產(chǎn)仔數(shù)的分子標(biāo)記[8]。陳祥等對(duì)220只產(chǎn)羔記錄完整的黔北麻羊FSHR基因SNPs位點(diǎn)的不同基因型與繁殖性狀進(jìn)行分析，在C1246A位點(diǎn)，黔北麻羊群體中BB型個(gè)體的產(chǎn)羔數(shù)與AA型和AB型差異顯著（P<0.05）[10]。Pan等對(duì)FSHR基因的5′側(cè)翼區(qū)與產(chǎn)羔數(shù)進(jìn)行關(guān)聯(lián)分析，發(fā)現(xiàn)綿羊FSHR基因與綿羊產(chǎn)羔數(shù)極顯著相關(guān)（P<0.01），CC的基因型比TC和TT多0.42（P<0.01）和0.53（P<0.01），表明綿羊FSHR基因可以作為提高綿羊產(chǎn)羔數(shù)候選基因[20]。

由上述研究可推測(cè)，F(xiàn)SHR基因多態(tài)性與高山美利奴羊、青海細(xì)毛羊、中國(guó)美利奴羊、敖漢細(xì)毛羊的多胎性狀也存在相關(guān)性。本研究通過(guò)分析4個(gè)品種全基因組測(cè)序結(jié)果，首次對(duì)4個(gè)綿羊品種FSHR的全外顯子的單核苷酸多態(tài)性進(jìn)行研究。研究發(fā)現(xiàn)，高山美利奴羊FSHR基因在外顯子上存在3個(gè)多態(tài)位點(diǎn)（T904G、C975T、G1572A），中國(guó)美利奴羊FSHR基因在外顯子上存在7個(gè)多態(tài)位點(diǎn)（G149A、G813A、T817C、T904G、C975T、T1234C、G1572A）、敖漢細(xì)毛羊FSHR基因在外顯子上存在6個(gè)多態(tài)位點(diǎn)（G813A、T817C、T904G、C975T、T1234C、G1572A）。4個(gè)品種之間檢查結(jié)果與龍威海等對(duì)黔北麻羊FSHR基因的檢測(cè)結(jié)果[21]及王惠娥等對(duì)多浪羊和卡拉庫(kù)爾羊FSHR基因的檢測(cè)結(jié)果[22]存在差異，可能是由地域不同造成的自然選擇不同或者是品種的不同造成的?；蛲怙@子編碼氨基酸序列，外顯子上的單核苷酸突變可能會(huì)造成RNA二級(jí)結(jié)構(gòu)、蛋白質(zhì)二級(jí)結(jié)構(gòu)、蛋白質(zhì)三級(jí)結(jié)構(gòu)改變，進(jìn)而對(duì)其生物學(xué)功能造成影響[23]。本研究所檢測(cè)7個(gè)SNPs中T817C、T904G、T1234C既使mRNA的二級(jí)結(jié)構(gòu)發(fā)生改變，又改變了最小自由能，G813A只改變最小自由能，mRNA二級(jí)結(jié)構(gòu)未發(fā)生改變，G149A與C975T、G1572A既沒(méi)有使最小自由能發(fā)生改變，又沒(méi)有改變mRNA的二級(jí)結(jié)構(gòu)。堿基突變會(huì)導(dǎo)致蛋白質(zhì)二級(jí)結(jié)構(gòu)發(fā)生改變，本研究對(duì)4個(gè)品種所檢測(cè)的7個(gè)多態(tài)位點(diǎn)中5個(gè)錯(cuò)義突變位點(diǎn)（G149A、G813A、T904G、C975T、G1572A）進(jìn)行蛋白質(zhì)二級(jí)結(jié)構(gòu)分析發(fā)現(xiàn)，5個(gè)突變位點(diǎn)會(huì)導(dǎo)致 α-螺旋、無(wú)規(guī)則卷曲、擴(kuò)展鏈的比例發(fā)生改變，進(jìn)一步分析蛋白質(zhì)三級(jí)結(jié)構(gòu)發(fā)現(xiàn)錯(cuò)義導(dǎo)致蛋白質(zhì)三級(jí)結(jié)構(gòu)的改變，對(duì)FSHR基因功能是否造成影響還需要進(jìn)一步的試驗(yàn)研究。

本試驗(yàn)為高山美利奴羊、青海細(xì)毛羊、中國(guó)美利奴羊、敖漢細(xì)毛羊的FSHR基因多態(tài)性與多胎性相關(guān)性的研究提供了參考，后續(xù)試驗(yàn)研究將進(jìn)一步分析試驗(yàn)結(jié)果中的各品種FSHR基因SNPs位點(diǎn)與綿羊多胎性狀相關(guān)性，研究FSHR基因多態(tài)性對(duì)綿羊產(chǎn)羔數(shù)的影響。

參考文獻(xiàn)：

[1]Dattatreyamurty B，Smith R A，Zhang S B，et al. The size of the mature membrane receptor for follicle-stimulating hormone is larger than that predicted from its cDNA[J]. Journal of Molecular Endocrinology，1992，9（2）：115-115.

[2]Casarini L，Moriondo V，Marino M，et al. FSHR polymorphism p.N680S mediates different responses to FSH in vitro[J]. Molecular and Cellular Endocrinology，2014，393（1/2）：83-91.

[3]Jiang X L，Liu H L，Chen X Y，et al. Structure of follicle-stimulating hormone in complex with the entire ectodomain of its receptor[J]. Proceedings of the National Academy of Sciences of the United States of America，2012，109（31）：12491-12496.

[4]吳井生，王金玉. FSHR基因第10外顯子多態(tài)性及其與小梅山豬產(chǎn)仔數(shù)的相關(guān)性[J]. 中國(guó)農(nóng)業(yè)科學(xué)，2012，45（13）：2728-2736.

[5]邢耀潭. 天祝白牦牛FSH、FSHR基因多態(tài)性分析與雙胎性狀的關(guān)聯(lián)研究[D]. 蘭州：甘肅農(nóng)業(yè)大學(xué)，2011：22-26.

[6]Chu M X，Guo X H，F(xiàn)eng C J，et al. Polymorphism of 5′regulatory region of ovine FSHR gene and its association with litter size in Small Tail Han sheep[J]. Molecular Biology Reports，2012，39（4）：3721-3725.

[7]Jiang Y，Xie M，Chen W B，et al. The sheep genome illuminates biology of the rumen and lipid metabolism[J]. Science，2014，344（6188）：1168-1173.

[8]許 瑤，劉 暢，牛 熙，等. 香豬FSHR基因外顯子10的多態(tài)性及其與產(chǎn)仔數(shù)間的關(guān)聯(lián)分析[J]. 中國(guó)畜牧獸醫(yī)，2017，44（3）：799-806.

[9]雷雪芹，魏伍川，陳 宏. 6個(gè)牛品種在FSHR基因位點(diǎn)的遺傳關(guān)系及其多態(tài)對(duì)雙胎性狀的標(biāo)記[J]. 西北農(nóng)林科技大學(xué)學(xué)報(bào)（自然科學(xué)版），2004，32（7）：1-6.

[10]陳 祥，龍威海，孫振梅，等. 貴州地方山羊FSHR基因與繁殖性狀的相關(guān)性研究[J]. 農(nóng)業(yè)生物技術(shù)學(xué)報(bào)，2017，25（1）：94-101.

[11]王建英，曹頂國(guó)，劉昭華，等. 綿羊多羔性狀候選基因的研究進(jìn)展[J]. 草食家畜，2017（3）：1-8.

[12]Laan M，Grigorova M，Huhtaniemi I T. Pharmacogenetics of follicle-stimulating hormone action[J]. Current Opinion in Endocrinology Diabetes and Obesity，2012，19（3）：220-227.

[13]Boudjenah R，Molina-Gomes D，Torre A A，et al. Genetic polymorphisms influence the ovarian response to rFSH stimulation in patients undergoing in vitro fertilization programs with ICSI[J]. PLoS One，2012，7（6）：e38700.

[14]李俊杰，賈 青，張正珊. 雌性哺乳動(dòng)物促卵泡素受體與促黃體素受體的結(jié)構(gòu)與表達(dá)[J]. 生命的化學(xué)，2003，23（4）：289-291.

[15]Lend A K，Belousova A，Haller-Kikkatalo K A，et al. Follicle-stimulating hormone receptor gene haplotypes and male infertility in estonian population and meta-analysis[J]. systems Biology in Reproductive Medicine，2010，56（1）：84-90.

[16]Wu Q Y，Shai Y C，Xia X Y，et al. FSH and FSHR gene polymorphisms and male infertility：an update[J]. National Journal of Andrology，2015，21（11）：1031-1031.

[17]Sang L，Du Q Z，Yang W C，et al. Polymorphisms in follicle stimulation hormone receptor，inhibin alpha，inhibin bata A，and prolactin genes，and their association with sperm quality in Chinese Holstein bulls[J]. Animal Reproduction Science，2011，126（3/4）：151-156.

[18]Lazaros L，Xita N，Takenaka A，et al. Synergistic effect of follicle-stimulating hormone receptor and androgen receptor gene variants on semen quality[J]. Andrologia，2013，45（5）：339-344.

[19]Zhylkova I S，Sotnik N N，Yegunkova O V，et al. Analysis of single nucleotide polymorphisms G919A and a2039G of gene FSHR in infertile men[J]. Cytology and Genetics，2018，52（2）：132-138.

[20]Pan X，Liu S，Li F，et al. Erratum to：molecular characterization，expression profiles of the ovine FSHR gene and its association with litter size[J]. Molecular Biology Reports，2014，41（12）：7749-7754.

[21]龍威海，丁 玫，馮文武，等. 卵泡刺激激素受體基因（FSHR）在貴州地方山羊的多態(tài)性及表達(dá)研究[J]. 農(nóng)業(yè)生物技術(shù)學(xué)報(bào)，2015，23（6）：755-761.

[22]王惠娥，王 永，汪 瀾，等. 南疆綿羊FSHR基因第10外顯子多態(tài)性檢測(cè)及序列分析[J]. 黑龍江畜牧獸醫(yī)，2017，17）（17）：114-117，294.

[23]楊 紅，孫 鵑，裘貴華，等. 3個(gè)貴州地方雞種NKX2-5基因多態(tài)性及生物信息學(xué)分析[J]. 中國(guó)畜牧獸醫(yī)，2014，41（12）：209-214.