發(fā)情相關(guān)基因在小尾寒羊性腺軸組織中的表達譜研究

陳玉林 劉玉芳 儲明星

摘要 [目的]探究綿羊發(fā)情相關(guān)基因KITLG、MAPK1、NOTCH4和NR5A2在小尾寒羊性腺軸組織（大腦、小腦、下丘腦、垂體、子宮、卵巢、輸卵管）中的表達差異，為闡明綿羊高繁殖力的分子機理提供理論依據(jù)。[方法]以FecB BB型小尾寒羊為研究對象，利用實時熒光定量PCR檢測上述4個基因在小尾寒羊發(fā)情周期（黃體期和卵泡期）性腺軸相關(guān)的7種組織中的表達差異。[結(jié)果] KITLG基因在小尾寒羊卵泡期輸卵管組織中的表達量最高，其在黃體期小腦與垂體中的表達量顯著高于卵泡期（P<0.05）;MAPK1基因在小尾寒羊卵泡期大腦組織中的表達量最高且極顯著高于黃體期（P<0.01），其在黃體期下丘腦和小腦組織中的表達量顯著高于卵泡期（P<0.05），而其在卵泡期輸卵管組織中的表達量顯著高于黃體期（P<0.05）;NOTCH4基因在小尾寒羊黃體期下丘腦中的表達量最高，而其在黃體期子宮中的表達量顯著高于卵泡期（P<0.05）;NR5A2基因在小尾寒羊黃體期卵巢組織中的表達量最高，其在黃體期子宮組織中的表達量顯著高于卵泡期（P<0.05）。[結(jié)論]KITLG、MAPK1、NOTCH4和NR5A2基因可作為潛在的候選基因來調(diào)控綿羊繁殖力。該研究結(jié)果為進一步闡明綿羊高繁殖力的分子機理提供了新思路。

關(guān)鍵詞 綿羊;發(fā)情相關(guān)基因;性腺軸;組織表達

中圖分類號 S 826? 文獻標識碼 A? 文章編號 0517-6611（2022）01-0092-04

doi：10.3969/j.issn.0517-6611.2022.01.023

Study on the Expression Profile of Estrus-related Genes in the Gonadal Axis Tissues of Small Tail Han Sheep

CHEN Yu-lin1，2， LIU Yu-fang1，2， CHU Ming-xing1

（1.Key Laboratory of Animal Genetics， Breeding and Reproduction of Ministry of Agriculture and Rural Affairs， Institute of Animal Sciences， Chinese Academy of Agricultural Sciences， Beijing 100193;2. College of Life Science and Food Engineering， Hebei University of Engineering， Handan， Hebei 056001）

Abstract [Objective] To explore the expression difference of estrus-related genes KITLG， MAPK1， NOTCH4 and NR5A2 in the gonadal axis tissues （brain， cerebellum， hypothalamus， pituitary， uterus， ovary， fallopian tube） of Small Tail Han Sheep， and provide the theoretical basis for clarifying the molecular mechanism of high-fertility sheep. [Method]Using FecB BB type Small Tail Han Sheep as the research objects， real-time quantitative PCR was used to detect the expression difference of the four genes in seven kinds of tissues related to the gonadal axis of Small Tail Han Sheep in the estrus period（luteal phase and follicular phase）. [Result]The expression level of KITLG gene in the oviduct tissues of Small Tail Han Sheep in the follicular phase was the highest， its expression levels in the cerebellum and pituitary tissues of Small Tail Han Sheep in the luteal phase were significantly higher than that in the follicular phase （P<0.05）.The expression level of MAPK1 gene in the brain tissues of Small Tail Han Sheep in the follicular phase was the highest， which was extremely significantly higher than that in the luteal phase（P>0.01）. Its expression levels in hypothalamus and cerebellum tissues of Small Tail Han Sheep in the luteal phase were significantly higher than that in the follicular phase （P<0.05）， and its expression level in fallopian tube tissues in the? follicular phase was significantly higher than that in the luteal phase（P<0.05）. The expression level of NOTCH4 gene in the hypothalamus of Small Tail Han Sheep in the luteal phase was the highest， and the expression level in uterus tissues in the? luteal phase? was significantly higher than that in the follicular phase （P<0.05）. [Conclusion]KITLG， MAPK1， NOTCH4 and NR5A2 genes might be used as potential candidate genes to regulate the fecundity of sheep. The research provided a new idea for further clarifying the molecular mechanism of sheep’s high fertility.

Key words Sheep;Estrus-related genes;Gonadal axis;Tissue expression

基金項目 國家自然科學基金項目（31772580）;國家肉羊產(chǎn)業(yè)技術(shù)體系專項（CARS-38）;中國農(nóng)業(yè)科學院科技創(chuàng)新工程項目（CAAS-ZDRW202106，ASTIP-IAS13）;河北省自然科學基金青年項目（C2019402261）。

作者簡介 陳玉林（1997—），男，河南駐馬店人，碩士研究生，研究方向：動物遺傳育種。

通信作者：劉玉芳，講師，博士，從事動物遺傳育種研究;儲明星，研究員，博士，從事羊優(yōu)異繁殖性狀分子機理研究。

收稿日期 2021-05-19

如何提高綿羊的繁殖力一直是畜牧業(yè)關(guān)注的問題[1]。發(fā)情期是動物繁殖的重要時期，分為卵泡期和黃體期，卵母細胞必須經(jīng)過卵泡期的發(fā)育才能成熟進而排卵，黃體期黃體的生長發(fā)育則釋放黃體酮等激素，為母體妊娠提供適宜的環(huán)境，為進一步的繁殖做準備[2]。綿羊的繁殖力受到其發(fā)情方式的直接影響，常年發(fā)情綿羊的繁殖力通常比季節(jié)性發(fā)情綿羊有更大的優(yōu)勢[3]。因此，研究綿羊的發(fā)情機制及其如何影響綿羊的繁殖力一直是一個備受關(guān)注的問題。小尾寒羊作為一種優(yōu)質(zhì)的常年發(fā)情綿羊，其高繁殖力的特點使其成為研究綿羊發(fā)情機制及繁殖力機制的完美模型。筆者篩選出幾個與動物發(fā)情期相關(guān)的基因，旨在為闡明綿羊發(fā)情機制提供一定的理論基礎(chǔ)，為提高綿羊繁殖力提供新的研究思路。

KIT配體（KITLG）是c-KIT跨膜酪氨酸激酶受體（KIT）的配體，也被稱為干細胞因子、鋼鐵因子或肥大細胞生長因子。動物模型顯示，顆粒細胞衍生的KITLG與卵母細胞、卵泡膜細胞衍生的KITLG之間的相互作用涉及卵母細胞和卵泡發(fā)育的各個方面，包括卵巢中原始生殖細胞的形成、原始卵泡的激活、卵母細胞的存活和生長、顆粒細胞的增殖、膜細胞的募集和減數(shù)分裂停滯的維持[4]。MAPK1是MAPK（絲裂原活化蛋白激酶）的成員之一，在MAPK通路中發(fā)揮其作用，是卵泡發(fā)育的一條重要途徑。MAPK通路的激活可以促進瘦素誘導的卵母細胞成熟[5]。NOTCH4是一個保守的跨膜受體家族成員，可通過與一些特定的配體相互作用調(diào)節(jié)細胞，在許多發(fā)育過程中起關(guān)鍵作用，影響細胞分化、增殖和凋亡[6-8]。Hernandez等[9]研究表明NOTCH4通過Notch信號通路發(fā)揮促進黃體細胞活力和類固醇合成的作用。NR5A2也被稱為胎兒蛋白轉(zhuǎn)錄因子（FTF）、肝臟受體同源1（LRH-1），是核激素受體（NR）轉(zhuǎn)錄因子亞家族的成員，NR5A2與哺乳動物的類固醇合成、細胞增殖、卵泡發(fā)育、排卵和生育等多種過程有關(guān)[10]。

筆者利用實時熒光定量PCR（qPCR）技術(shù)對上述發(fā)情期相關(guān)候選基因在綿羊不同發(fā)情期（卵泡期和黃體期）性腺軸組織的表達情況進行檢測，為進一步揭示KITLG、MAPK1、NOTCH4和NR5A2基因的功能和解析綿羊繁殖機制提供理論基礎(chǔ)。

1 材料與方法

1.1 試驗樣品的收集

隨機選取3周歲健康、經(jīng)產(chǎn)的空懷小尾寒羊卵泡期和黃體期母羊各3只，于2017年10月飼養(yǎng)于天津市畜牧獸醫(yī)研究所試驗羊場，所有試驗羊的飼養(yǎng)環(huán)境和飼料均相同。2017年11月對試驗羊進行屠宰，取其大腦、小腦、下丘腦、垂體、子宮、卵巢和輸卵管，取樣后迅速裝入1.8 mL RNase-Free凍存管中（最大樣品量為凍存管體積的2/3）。所有樣品采集要在屠宰后30 min內(nèi)完成，樣品采集完后迅速放入液氮中冷凍保存，用干冰帶回實驗室，放入-80 ℃冰箱中冷凍保存，備用。

1.2 RNA的提取及質(zhì)量檢測

將采集的小尾寒羊卵泡期和黃體期性腺軸7種組織研磨后，用Trizol試劑（Invitrogen，美國）進行裂解，此后按照動物組織RNA提取試劑盒（天根，北京）的說明書進行總RNA的提取，用1.0%瓊脂糖凝膠電泳和Nanodrop 2000檢測提取RNA的質(zhì)量和濃度。經(jīng)檢驗合格的組織總RNA于-80 ℃下保存?zhèn)溆谩?/p>

1.3 引物設(shè)計

根據(jù)GenBank提供的綿羊KITLG、MAPK1、NOTCH4、NR5A2基因序列（登錄號分別為NM_001267888.1、XM_027956867.1、XM_027958770.1、NM_001246280.1）信息，并使用Primer Premier 5.0軟件進行引物設(shè)計，以RPL19為內(nèi)參基因。引物由北京天一輝遠生物科技有限公司合成。各引物濃度均為10 μmol/L，引物序列詳細信息見表1。

1.4 cDNA的合成

利用TaKaRa反轉(zhuǎn)錄試劑盒（PK0445）合成cDNA，反轉(zhuǎn)錄體系如下：1.0 μL反轉(zhuǎn)錄混合引物，1.0 μL Oligo dT引物，1.0 μL隨機引物，4.0 μL熒光定量緩沖液，1.0 μL RNA，12.0 μL雙蒸水，全程在冰上操作。反應(yīng)條件：37 ℃ 15 min，85 ℃ 5 s，將反轉(zhuǎn)錄完成后的cDNA產(chǎn)物稀釋后，用持家基因RPL19進行PCR檢測，質(zhì)量合格后-20 ℃下保存，以備檢測基因mRNA表達。

1.5 實時熒光定量PCR

1.5.1 實時熒光定量PCR體系及反應(yīng)程序。

實時熒光定量PCR反應(yīng)體系（總體積20 μL）：10.0 μL熒光定量Ex Taq引物Ⅱ，0.8 μL上游引物，0.8 μL下游引物，2.0 μL cDNA 模板，6.4 μL雙蒸水。PCR反應(yīng)程序：95 ℃預變性5 min;95 ℃變性5 s，60 ℃ 30 s，40個循環(huán)。反應(yīng)結(jié)束后進行熔解曲線分析。

1.5.2 標準曲線的建立。

將cDNA樣本5倍稀釋后，進行2倍梯度稀釋后獲得5個濃度梯度（1、1/2、1/4、1/8、1/16）的cDNA樣品。用這些cDNA為模板，對目的基因和持家基因進行熒光定量PCR，以濃度梯度的對數(shù)值（以10為底數(shù)）為橫坐標，以檢測所得Ct值為縱坐標，繪制目的基因和持家基因的標準曲線。

1.5.3 實時熒光定量檢測與數(shù)據(jù)分析。

使用Roche Light Cycler480 Ⅱ 型熒光定量PCR儀進行熒光定量檢測，使用Excel軟件進行數(shù)據(jù)處理，采用2-△△Ct的方法計算各基因的相對表達量。獲得的表達量數(shù)據(jù)使用SPSS 19.0統(tǒng)計軟件進行獨立樣本t檢驗。

2 結(jié)果與分析

2.1 RNA的提取

提取后的總RNA用1%瓊脂糖凝膠電泳進行檢測，結(jié)果發(fā)現(xiàn)28S條帶明顯比18S條帶更亮，且條帶完整性較好，表明該RNA質(zhì)量合格，能用于后續(xù)試驗。

2.2 4個基因在各組織中表達分析

2.2.1 KITLG基因組織表達檢測。

qPCR結(jié)果顯示，KITLG基因在7種組織中均有表達，其在小尾寒羊卵泡期輸卵管中的相對表達量最高，其在黃體期小腦與垂體組織中的相對表達量均大于卵泡期且差異顯著（P<0.05），在其他組織中的相對表達量差異均不顯著（P>0.05）（圖1）。

2.2.2 MAPK1基因組織表達檢測。

MAPK1基因在小尾寒羊卵泡期大腦組織中的相對表達量最高且極顯著高于黃體期（P<0.01），在黃體期下丘腦組織中的相對表達量極顯著高于卵泡期（P<0.01）;其在黃體期小腦組織中的相對表達量顯著高于卵泡期（P<0.05），其在卵泡期輸卵管組織中的相對表達量顯著高于黃體期（P<0.05）;其在其他組織中的相對表達量差異均不顯著（P>0.05）（圖2）。

2.2.3 NOTCH4基因組織表達檢測。

如圖3所示，NOTCH4基因在小尾寒羊黃體期下丘腦組織中的相對表達量最高，其在黃體期子宮組織中的相對表達量顯著高于卵泡期（P<0.05），而其在其他組織中的相對表達量差異均不顯著（P>0.05）。

2.2.4 NR5A2基因組織表達檢測。

如圖4所示，NR5A2基因在小尾寒羊黃體期卵巢組織中的相對表達量最高，其在黃體期子宮組織中的相對表達量顯著高于卵泡期（P<0.05），而其在其他組織中的相對表達量均無顯著差異（P>0.05）。

3 討論

動物發(fā)情是研究動物繁殖育種不可忽視的一環(huán)，其中發(fā)情方式、發(fā)情時間均受到多種遺傳因素的影響，研究綿羊發(fā)情期相關(guān)基因調(diào)控機制對于提高綿羊繁殖力有著重大意義[11]。卵泡在發(fā)情期的生長發(fā)育狀況直接決定了綿羊的排卵數(shù)。排卵數(shù)是決定綿羊產(chǎn)羔數(shù)的直接因素，卵泡的生長發(fā)育與排卵數(shù)有著直接的聯(lián)系[12-14]。性腺軸是以下丘腦-垂體-卵巢為主的主要參與調(diào)控卵巢內(nèi)卵泡及黃體生長發(fā)育的相關(guān)組織，研究性腺軸組織中發(fā)情相關(guān)基因的表達為進一步探究提高綿羊產(chǎn)羔數(shù)的方法提供了新思路[15]。

KITLG作為KIT的配體，是細胞表面錨定因子家族成員之一[16]。KIT/KITLG信號通路是生殖細胞生長發(fā)育過程中不可或缺的，其與生殖細胞的發(fā)育密切相關(guān)[17]。研究表明，KITLG基因在多種動物卵巢組織中均有表達，且對排卵前卵泡的生長發(fā)育起重要作用，已有試驗證實KITLG基因多態(tài)性與綿羊產(chǎn)羔數(shù)具有顯著關(guān)聯(lián)[18-19]。該試驗中KITLG基因在性腺軸7種組織中均有表達，其在小腦和垂體中的相對表達量在2個時期間存在顯著差異（P<0.05），其在垂體中的相對表達量較少，其在2個時期卵巢組織中均有表達，但差異未達到顯著水平（P>0.05），推測該基因不參與卵泡期與黃體期的轉(zhuǎn)換調(diào)控過程，其在不同時期的調(diào)控機制有待進一步探究。

絲裂原活化蛋白激酶1（MAPK1）作為MAPK通路中的重要成員，參與生殖過程中的信號轉(zhuǎn)導，可激活卵母細胞的成熟進程，使其經(jīng)歷生發(fā)泡破裂等一系列過程[20-21]。該試驗中MAPK1基因在卵泡期大腦組織中的相對表達量最高且極顯著高于黃體期（P<0.01），其在黃體期下丘腦組織中的相對表達量極顯著高于卵泡期（P<0.01），黃體期小腦組織中該基因的相對表達量顯著高于卵泡期（P<0.05），在卵泡期輸卵管組織中該基因的相對表達量顯著高于黃體期（P<0.05），但除了卵泡期大腦組織外該基因在2個時期性腺軸各組織中的相對表達量均不高，推測其在小尾寒羊性腺軸對發(fā)情期的調(diào)控過程中沒有起到重要的調(diào)控作用。

NOTCH4在NOTCH通路中作為受體發(fā)揮作用，NOTCH信號通路在許多器官與細胞的發(fā)育中非常重要。雷俊川[22]研究表明其在哺乳動物卵泡發(fā)育中發(fā)揮重要作用。該試驗中NOTCH4基因在小尾寒羊不同時期性腺軸7種組織中均有表達。NOTCH4基因在小尾寒羊黃體期子宮組織中的相對表達量顯著高于卵泡期（P<0.05），且其在卵泡期子宮中的相對表達量較低。根據(jù)前人研究結(jié)果，推測出其在子宮妊娠環(huán)境的建立過程中參與調(diào)控[23]。

NR5A2由4個結(jié)構(gòu)域組成：DNA結(jié)合（DBD）和配體結(jié)合（LBD）結(jié)構(gòu)域、連接DBD和LBD的長Hinger（HR）以及修飾結(jié)構(gòu)域（MD）[24]。在成年哺乳動物中，NR5A2基因主要在前脂肪細胞以及肝、腸、胰腺、卵巢和睪丸中均有表達，尤其在卵巢中高表達，該基因是卵泡發(fā)育和排卵的重要和多效性調(diào)節(jié)因子[7]。該試驗中NR5A2基因在小尾寒羊黃體期和卵泡期大腦、小腦和垂體組織中低表達，其在黃體期子宮組織中的相對表達量顯著高于卵泡期（P<0.05），其在卵巢組織中的相對表達量最高，但黃體期與卵泡期之間無顯著差異。結(jié)合前人研究結(jié)果，推測出該基因在小尾寒羊卵泡和黃體生長發(fā)育過程中起到重要作用，其具體作用機制有待進一步研究。

4 結(jié)論

該試驗利用實時熒光定量PCR研究了KITLG、MAPK1、NOTCH4和NR5A2基因在小尾寒羊卵泡期和黃體期性腺軸7種組織中的表達水平，結(jié)果發(fā)現(xiàn)KITLG、MAPK1、NOTCH4和NR5A2基因在卵泡期、黃體期性腺軸相關(guān)組織中的表達量有較大差異，這4個基因可能對綿羊發(fā)情周期的轉(zhuǎn)換具有一定的調(diào)控作用，可為進一步闡明綿羊發(fā)情周期轉(zhuǎn)換的分子機理提供參考。

參考文獻

[1] 樊宇嵐.綿羊多羔性狀候選基因的研究及發(fā)展[J].農(nóng)業(yè)技術(shù)與裝備，2019（11）：119-120.

[2] 田占偉，趙宗勝，林杉，等.新疆哈薩克綿羊在乏情期和發(fā)情期生殖激素的變化規(guī)律[J].江蘇農(nóng)業(yè)科學，2016，44（2）：256-259.

[3] 夏青，劉秋月，王翔宇，等.綿羊季節(jié)性繁殖分子機制及休情季節(jié)誘導綿羊發(fā)情配種技術(shù)[J].遺傳，2018，40（5）：369-377.

[4] AN X P，HOU J X，LEI Y N，et al.Two mutations in the 5′-flanking region of the KITLG gene are associated with litter size of dairy goats[J].Anim Genet，2015，46（3）：308-311.

[5] DI R，HE J N，SONG S H，et al.Characterization and comparative profiling of ovarian microRNAs during ovine anestrus and the breeding season[J].BMC Genom，2014，15（1）：1-15.

[6] CHENG R，CAI X R，KE K，et al.Notch4 inhibition suppresses invasion and vasculogenic mimicry formation of hepatocellular carcinoma cells[J].Curr Med Sci，2017，37（5）：719-725.

[7] ZHANG Z X，BU X，YANG J Y，et al.NOTCH4 regulates colorectal cancer proliferation，invasiveness，and determines clinical outcome of patients[J].J Cell Physiol，2018，233（10）：6975-6985.

[8] PURI R V，YERRATHOTA S，HOME T，et al.Notch4 activation aggravates NF-κB-mediated inflammation in HIV-1 associated Nephropathy[J].Dis Models Mech，2019，12（12）：1-31.

[9] HERNANDEZ F，PELUFFO M C，STOUFFER R L，et al.Role of the DLL4-NOTCH system in PGF2alpha-induced luteolysis in the pregnant rat[J].Biol Reprod，2011，84（5）：859-865.

[10] LI Y X，ZHANG J，QIAN Y，et al.Molecular characterization，expression，polymorphism of NR5A2 and its relationship with litter size in Hu sheep[J].Genet Mol Res，2015，14（4）：12765-12775.

[11] LI X Y，LI C Y，XU Y R，et al.Analysis of pituitary transcriptomics indicates that lncRNAs are involved in the regulation of sheep estrus[J].Funct Integr Genom，2020，20（4）：563-573.

[12] 郭晶，李新宇，李隱俠，等.湖羊TGF-β1基因特征、表達及其與排卵數(shù)的相關(guān)性分析[J].中國農(nóng)業(yè)科學，2013，46（21）：4586-4593.

[13] 王翔宇，狄冉，劉秋月，等.小鼠和綿羊BMP15和GDF9突變在卵泡發(fā)育中作用及其對排卵數(shù)的影響[J].農(nóng)業(yè)生物技術(shù)學報，2020，28（10）：1870-1880.

[14] 徐業(yè)芬，李齊發(fā)，李二林，等.湖羊BMP2、BMP4、BMP6和BMP7基因mRNA表達水平與排卵數(shù)關(guān)系的研究[J].中國農(nóng)業(yè)科學，2009，42（10）：3655-3661.

[15] PLANT T M.60 YEARS OF NEUROENDOCRINOLOGY：The hypothalamo-pituitary-gonadal axis[J].J Endocrinol，2015，226（2）：T41-T54.

[16] 侯佳妮.藍狐KIT、KITLG及EDNRB基因的克隆及與顯性白毛色突變的相關(guān)性分析[D].長春：吉林大學，2014.

[17] 張敏，丁顯平，朱一劍，等.KIT基因單核苷酸多態(tài)性與四川地區(qū)男性不育的相關(guān)性研究[J].四川大學學報（自然科學版），2011，48（4）：961-966.

[18] 張軍霞，樂祥鵬，王維民，等.綿羊KITLG基因的表達及其多態(tài)性與產(chǎn)羔性狀關(guān)聯(lián)分析[J].農(nóng)業(yè)生物技術(shù)學報，2017，25（6）：893-900.

[19] 安小鵬.山羊產(chǎn)羔性狀候選基因的篩選及其多基因聚合效應(yīng)的研究[D].楊凌：西北農(nóng)林科技大學，2013.

[20] PONZA P，YOCAWIBUN P，SITTIKANKAEW K，et al.Molecular cloning and expression analysis of the Mitogen-activating protein kinase 1（MAPK1）gene and protein during ovarian development of the giant tiger shrimp Penaeus monodon[J].Mol Reprod Dev，2011，78（5）：347-360.

[21] GONZALEZ-GARCIA J R，BRADLEY J，NOMIKOS M，et al.The dynamics of MAPK inactivation at fertilization in mouse eggs[J].J Cell Sci，2014，127（Pt 12）：2749-2760.

[22] 雷俊川.Notch信號途徑在哺乳動物卵泡發(fā)育中作用的初步研究[D].西安：第四軍醫(yī)大學，2008.

[23] 王會平，石瑛，張展.Notch信號通路在著床和胎盤形成中的表達及作用[J].廣東醫(yī)學，2016，37（7）：1099-1102.

[24] YUMOTO F，NGUYEN P，SABLIN E P，et al.Structural basis of coactivation of liver receptor homolog-1 by β-catenin[J].PNAS，2012，109（1）：143-148.