小尾寒羊間情期和發(fā)情期microRNAs差異表達分析

段新崇，魏彥輝，李　陽，李相運，周榮艷，錫建中

(河北農(nóng)業(yè)大學動物科技學院，保定 071001)

?

小尾寒羊間情期和發(fā)情期microRNAs差異表達分析

段新崇，魏彥輝，李陽，李相運，周榮艷*，錫建中

(河北農(nóng)業(yè)大學動物科技學院，保定 071001)

摘要：本研究通過構(gòu)建和分析小尾寒羊間情期和發(fā)情期卵巢組織microRNA(miRNA)表達譜，篩選兩個時期差異表達的microRNAs，為研究microRNA調(diào)控小尾寒羊繁殖過程提供相應的理論基礎。利用活體手術(shù)法在發(fā)情期(第1天)和間情期(第13天)分別采集一側(cè)卵巢。從卵巢中提取總RNA，利用illumina Hiseq2000測序平臺獲取RNA數(shù)據(jù)，對表達譜和差異表達microRNAs進行生物信息學分析。實時熒光定量PCR(qRT-PCR)驗證差異表達的microRNAs在小尾寒羊卵巢中表達水平。結(jié)果，成功地構(gòu)建出小尾寒羊間情期和發(fā)情期microRNA的表達譜，oar-miR-99a和oar-miR-143分別是間情期和發(fā)情期表達量最高的microRNA，并篩選出在兩個時期間3個顯著差異表達的microRNAs，分別為oar-miR-200a、oar-miR-200b和oar-miR-200c。利用qRT-PCR對隨機選擇的2個顯著差異表達的microRNAs進行驗證，其表達水平和RNA-Seq分析結(jié)果一致。microRNA表達譜為后續(xù)的綿羊卵巢microRNA研究提供更詳盡的信息。結(jié)合靶基因預測及通路富集分析，推測差異表達的microRNAs是通過代謝和免疫途徑調(diào)控卵巢周期性活動。

關鍵詞：小尾寒羊；間情期；發(fā)情期；卵巢；microRNA

microRNA是一類長約22個核苷酸的非編碼單鏈RNA分子，在哺乳動物的卵巢、輸卵管、子宮、胎盤等繁殖組織器官中廣泛表達[1-3]。它通過與靶mRNA特異性堿基互補配對，引起靶mRNA降解或者抑制其翻譯，發(fā)揮基因轉(zhuǎn)錄后水平調(diào)控作用，從而參與細胞分化、增殖、凋亡等一系列生物學過程，而這些生物學過程正是卵泡發(fā)生、排卵以及黃體生成所必須的過程[4]。研究發(fā)現(xiàn)microRNA參與動物的多項繁殖活動。首先，microRNA影響激素合成：如對體外培養(yǎng)的人顆粒細胞進行miR-15a、Let-7b、Let-7c、miR-17-3p等36種microRNAs的轉(zhuǎn)染，孕酮釋放量提高了1.3～2.0倍；miR-16、miR-24、miR-25、miR-122等10種microRNAs轉(zhuǎn)染使孕酮釋放量減少了5倍以上；miR-15a、miR-24、miR-25、miR-26a等51種microRNAs對雌激素分泌有抑制作用，最高達4.5倍[5]。在豬有腔卵泡生長過程中，顆粒細胞中的miR-378表達水平上升，抑制芳香酶基因表達和雌激素的分泌[6]。其次，microRNA調(diào)節(jié)顆粒細胞、卵泡、黃體發(fā)育。調(diào)控顆粒細胞凋亡的microRNAs主要有miR-145、miR-23a和miR-26b，其中miR-145通過靶作用于激活素受體IB和細胞周期蛋白D2基因而抑制顆粒細胞增殖[7]。LH 峰后顆粒細胞內(nèi)的miR-21、miR-212表達量增加，miR-21有防止排卵卵泡的顆粒細胞在黃體化過程中發(fā)生凋亡的作用，其調(diào)控的機理尚不清楚[8]；miR-212能夠上調(diào)受體轉(zhuǎn)錄因子CTBP-1在小鼠顆粒細胞核內(nèi)的表達水平[9]。F.Muramatsu等發(fā)現(xiàn)miR-125b、miR-145、miR-31、miR-503和miR-21在反芻動物卵泡與黃體中表達差異顯著[10]。最近R.Di等首次確定了灘羊非發(fā)情季節(jié)卵巢microRNAs表達譜，表達譜中共202個microRNAs，包括63個已知的，136個保守的和3個新預測的microRNAs，miR-n-142是非發(fā)情季節(jié)表達最豐富的microRNA，其靶基因主要富集在氧化磷酸化、甘油脂代謝和磷脂酰肌醇信號通路中[11]?？梢妋icroRNAs在雌性哺乳動物生殖過程中具有極其重要的作用。

小尾寒羊是冀、魯、豫、蘇、皖等地的優(yōu)良品種，在繁殖性能中具有較為突出的表現(xiàn)，如成熟早、四季發(fā)情、產(chǎn)羔率高。其間情期與發(fā)情期卵巢經(jīng)過了一系列變化過程，并最終決定了高繁殖率。建立小尾寒羊間情期與發(fā)情期microRNAs的表達譜并篩選兩個時期差異表達的microRNAs，為研究microRNAs在調(diào)節(jié)卵泡發(fā)育與排卵生物學過程中的作用機制提供理論依據(jù)。

1材料與方法

1.1試驗材料

選擇3只經(jīng)產(chǎn)小尾寒羊母羊，用公羊試情判斷母羊發(fā)情。利用活體外科手術(shù)法在3只母羊的間情期和發(fā)情期(間隔12 d)分別采取一側(cè)卵巢，立即投入液氮速凍，帶回實驗室置于-80 ℃保存。

1.2總RNA的提取、文庫構(gòu)建與測序

采用動物組織總RNA提取試劑盒(CWbio.Co.Ltd，China)，參照試劑盒說明書提取總RNA。并對樣本總RNA進行以下各項質(zhì)量檢測：利用1%瓊脂糖凝膠電泳檢測RNA降解程度以及是否有污染；RNA純度由NanoPhotometer(IMPLEN，CA，USA)進行檢測；Qubit?2.0(Life Technologies，CA，USA)檢測RNA濃度；最后使用Agilent(Agilent Technologies，CA，USA)對RNA的完整性進行檢測。

樣品檢測合格后，小RNA(sRNA)文庫及高通量測序由北京諾禾致源生物信息有限公司完成。每個樣品需3 μg的總RNA，使用Small RNA Sample Pre Kit(NEB，USA)構(gòu)建文庫，利用sRNA的3′及5′端特殊結(jié)構(gòu)，以總RNA為起始樣品在sRNA兩端加接頭并反轉(zhuǎn)錄成cDNA，隨后經(jīng)過PCR擴增，PAGE膠電泳分離目標DNA片段，切膠回收得到的即為cDNA文庫。按照有效濃度及目標數(shù)據(jù)量的要求pooling后進行RNA-Seq測序，得到的序列原始數(shù)據(jù)(Raw Data)用于基本數(shù)據(jù)的分析。

1.3生物信息學分析鑒定sRNA

對原始數(shù)據(jù)進行去除低質(zhì)量、去除接頭、去除污染、去除堿基信息不明確的序列，得到干凈的的序列；對sRNA的長度分別進行統(tǒng)計，長度分布峰可以幫助我們判斷sRNA的種類，篩選一定長度范圍的sRNA進行后續(xù)分析；用bowtie軟件將sRNA定位到參考序列上，將序列比對到RepeatMasker和Rfam本物種的數(shù)據(jù)庫，去除sRNAs中相應的重復序列以及rRNA、tRNA、SnRNA和SnoRNA等非編碼的sRNA；剩余的sRNA對比到miRBase20.0尋找到已知的microRNAs；運用mirdeep2軟件預測新的microRNAs；對各樣本中已知和新的microRNAs進行表達量統(tǒng)計，并進行表達量TPM歸一化處理；將microRNAs表達水平分析中得到的readcount數(shù)據(jù)，采用基于負二項分布的DESeq[12-13]進行分析，從差異倍數(shù)(Fold change)和校正后的顯著水平(q value<0.01)進行評估，對差異microRNAs進行篩選。本研究進行生物信息分析時所參考的數(shù)據(jù)庫和分析程序信息見表1。

表1　生物信息學分析軟件

1.4real-time PCR驗證差異表達的microRNAs

對小尾寒羊兩個時期內(nèi)顯著差異表達的oar-miR-200a和oar-miR-200b進行qRT-PCR檢測。以U6基因為內(nèi)參，microRNA及內(nèi)參基因探針見表2。使用microRNA第一條鏈合成試劑盒(CWbio.Co.Ltd，China) 合成的cDNA作為qRT-PCR的模板。PCR反應體系20 μL：2 μL cDNA、40 μmol·L-1正反引物、10 μL的SYBR Premix Ex Taq(寶生物工程有限公司，大連)。PCR反應條件：95 ℃預變性10 min；95 ℃變性15 s，60 ℃退火30 s，共14個循環(huán)；72 ℃延伸3 min；4 ℃保存。相對表達量結(jié)果采用2-ΔΔCt法進行計算分析。

表2　microRNAs定量檢測的引物序列

2結(jié)果

2.1純化原始序列及sRNA序列長度分布

原始數(shù)據(jù)經(jīng)過去除低質(zhì)量、接頭、污染等過程后得到干凈序列(clean reads)，發(fā)情期(PXW)和間情期(DXW)分別得到clean reads為11 467 290和10 667 271個。sRNA長度分布主要是21～23 bp，測序數(shù)據(jù)中間情期和發(fā)情期長度為22 bp的序列分別占44.55%和43.49%(圖1)。統(tǒng)計兩樣品間公共序列和特有序列的種類(uniq)及數(shù)量(total)，兩個品種間公共序列的種類僅占全部序列的13.05%，而公共序列的數(shù)量卻占全部數(shù)量的93.39%(圖2)。

2.2sRNA與參考序列進行比對

使用Bowtie將長度篩選后的sRNA reads定位到參考序列上，分析sRNA在參考序列上的分布情況，可以看出大部分的sRNA都可以匹配到參考序列上，間情期和發(fā)情期沒有匹配上的sRNA分別僅占12.03%和18.76%(表3)。

2.3microRNAs在間情期和發(fā)情期的表達譜分析

間情期與發(fā)情期卵巢組織sRNA與miRbase數(shù)據(jù)庫進行比對，得到匹配上的已知microRNAs分別有129、135個。兩個時期表達量經(jīng)TPM標準化處理后，排在前10位的microRNAs相同(圖3)，分別為oar-miR-143、oar-miR-21、oar-miR-99a、oar-miR-26a、oar-miR-148a、oar-miR-10b、oar-let-7i、oar-let-7f、oar-let-7g、oar-miR-199a-3p。而間情期和發(fā)情期表達量TPM值最高的microRNA分別是oar-miR-99a(173 266)和oar-miR-143(187 028)。

表3　與參考基因組比對信息統(tǒng)計表

圖1　小RNA長度分布Fig.1　Length distribution of small RNA sequencing

A.兩個時期公共及特有序列種類 (Uniq sRNA)；B.兩個時期公共及特有序列數(shù)量統(tǒng)計(Total sRNA)A.Venn chart for uniq sRNAs in 2 periods；B.Venn chart for total sRNAs in 2 periods圖2　兩個時期公共及特有的序列數(shù)量及種類Fig.2　Venn chart for total and uniq sRNAs in 2 periods

圖3　小尾寒羊間情期和發(fā)情期卵巢組織前20位的microRNAs表達水平柱形圖Fig.3　Bar diagram of top 20 microRNAs expressed in estrus and diestrus ovary in the Small Tailed Han sheep

2.4差異表達microRNAs的篩選

小尾寒羊卵巢組織microRNAs表達量經(jīng)TPM標準化后，獲得3個顯著差異表達的microRNAs，在發(fā)情期較間情期全部表達下調(diào)(log2Fold change<-1)(圖4)，分別是oar-miR-200a由403.6下調(diào)至120.4；oar-miR-200b由119.3下調(diào)至26.7；oar-miR-200c由22.4下調(diào)至3.1，三者均屬于miR-200家族。

方塊表示顯著差異表達的microRNAs (log2Fold change<-1，p value<0.01，q value<0.01)square shape indicate significant difference in expression of microRNAs(log2Fold change<-1，p value<0.01，q value<0.01)圖4　顯著差異表達microRNAs分析火山圖Fig.4　The volcano plot of significant differentially expressed microRNAs

2.5新microRNAs的預測

使用microRNA預測軟件miREvo[12]和mirdeep2[13]進行新microRNA的分析，截取一定長度sRNA比對上的參考序列，通過對其二級結(jié)構(gòu)及Dicer酶切位點信息、能量等特征進行分析，預測樣品中新microRNA，此次測序中小尾寒羊間情期和發(fā)情期共發(fā)現(xiàn)68個新microRNAs：間情期57個，發(fā)情期59個。Novel_67在兩個時期中表達量TPM值最高，分別是201.4和200.9(圖5)。

圖5　小尾寒羊間情期和發(fā)情期卵巢組織新microRNAs表達水平柱形圖Fig.5　Bar diagram of novel microRNAs expressed in estrus and diestrus ovary in the Small Tailed Han sheep

2.6qRT-PCR驗證結(jié)果

在小尾寒羊間情期和發(fā)情期卵巢組織顯著差異表達的microRNAs中選取oar-miR-200a和oar-miR-200b進行熒光定量PCR(qRT-PCR)驗證。結(jié)果表明oar-miR-200a和oar-miR-200b在小尾寒羊卵巢中皆有表達(表4)，并且二者在發(fā)情期較間情期均表達下調(diào)，這與RNA-Seq測序結(jié)果基本一致。

表4　qRT-PCR與RNA-Seq結(jié)果比較

2.7差異microRNAs的靶基因預測和 GO 功能注釋

根據(jù)已有數(shù)據(jù)，oar-miR-200a、oar-miR-200b、oar-miR-200c預測得到的靶基因個數(shù)分別為：540、583和687個。對這些靶基因進行GO富集分析，結(jié)合圖6和在線分析結(jié)果可知，在細胞組分分類中，差異表達microRNAs的靶基因絕大部分參與形成細胞成分(占70%)，其次參與形成細胞內(nèi)組分(占35%)；在分子功能分類中，絕大部分靶基因參與連接分子類(占71%)，其次是催化活性分子類(占55%)；而在生物學過程分類中，參與有機物代謝過程的靶基因較多(占52%)，然后是初級代謝過程(占51%)。

2.8差異microRNAs靶基因KEGG富集分析

在生物體內(nèi)不同基因相互協(xié)調(diào)行使其生物學功能，通過通路顯著性富集能確定基因參與的最主要生化代謝途徑和信號轉(zhuǎn)導途徑。對差異microRNAs靶基因進行KEGG富集分析(P<0.05，F(xiàn)DR<0.4)，結(jié)果表明，差異microRNAs的靶基因在免疫、凝血、脂肪代謝以及糖代謝等19個通路達到顯著富集水平。X.Y.Miao等對不同發(fā)情周期小尾寒羊和道賽特羊的轉(zhuǎn)錄組進行測序，根據(jù)KEGG通路分析提出代謝、免疫等途徑可能調(diào)控家畜繁殖力[14]，并且繁殖力與能量平衡、代謝關系密切[15-16]，我們篩選出4條與卵巢繁殖活動關系較大的通路分別是：脂肪酸代謝通路、I型糖尿病通路、補體和凝血級聯(lián)反應、細胞粘附分子通路(表5)。差異表達microRNAs的靶基因及參與的KEGG通路如圖7所示，所有差異表達microRNAs通過相同靶基因FGN、C1S參與補體凝血級聯(lián)反應，而且通過相同靶基因CD274、ITGA9、ITGA6、NEO1參與細胞粘附分子通路；DYA與OVAR-DQA2作為oar-miR-200c和oar-miR-200b共同的靶基因同時參與細胞粘附分子和I型糖尿病通路，并且oar-miR-200c和oar-miR-200b還通過共同的靶基因ACOX1、CPT1A、OXSM、SCD和ACACB參與脂肪酸代謝通路(圖7)。

1.DNA復制RNA合成；2.DNA復制；3.DNA代謝過程；4.生物合成過程；5.細胞生物合成過程；6.有機物質(zhì)生物合成的過程；7.細胞大分子生物合成過程：8.大分子生物合成過程；9.核酸代謝過程；10.氮的化合物代謝過程；11.初級代謝過程；12.細胞代謝過程；13.蛋白質(zhì)N-鏈糖基化；14.含堿基化合物代謝；15.細胞氮化合物代謝；16.雜環(huán)代謝過程；17.細胞大分子代謝過程；18.有機物質(zhì)代謝過程；19.有機環(huán)狀化合物代謝過程；20.細胞芳香族化合物代謝；21.高爾基體運輸復合體；22.細胞組分；23.胞質(zhì)部分；24.有膜細胞器；25.細胞內(nèi)的有膜細胞器；26.細胞；27.細胞部分；28.細胞外區(qū)域；29.細胞外部分；30.線粒體；31.細胞質(zhì)；32.高爾基體部分；33.細胞器；34.細胞內(nèi)；35.細胞內(nèi)的細胞器；36.線粒體膜部分；37.RNA-直接DNA聚合酶的活動；38.DNA聚合酶活動；39.核酸轉(zhuǎn)移酶活性；40.RNA結(jié)合；41.轉(zhuǎn)移酶活性；42.催化活性；43.核酸結(jié)合；44.β-1，4-甘露糖基-糖蛋白 4 -β-乙酰氨基葡萄糖轉(zhuǎn)移酶活性；45.轉(zhuǎn)移酶活性；46.乙酰葡糖氨基轉(zhuǎn)移酶活性；47.雜環(huán)化合物結(jié)合；48.有機環(huán)狀化合物結(jié)合；49.UDP葡萄糖基轉(zhuǎn)移酶活性；50.戊糖基轉(zhuǎn)移酶活性；51.一碳基團轉(zhuǎn)移酶活性；52.結(jié)合；53.受體結(jié)合1.RNA-dependent DNA replication；2.DNA replication；3.DNA metabolic process；4.Biosynthetic process；5.Cellular biosynthetic process；6.Organic substance biosynthetic process；7.Cellular macromolecule biosynthetic process；8.Macromolecule biosynthetic process；9.Nucleic acid metabolic process；10.Nitrogen compound metabolic process；11.Primary metabolic process；12.Cellular metabolic process；13.Protein N-linked glycosylation；14.Nucleobase-containing compound metabolic；15.Cellular nitrogen compound metabolic；16.Heterocycle metabolic process；17.Cellular macromolecule metabolic process；18.Organic substance metabolic process；19.Organic cyclic compound metabolic process；20.Cellular aromatic compound metabolic；21.Golgi transport complex；22.Cellular component；23.Cytoplasmic part；24.Membrane-bounded organelle；25.Intracellular membrane-bounded organelle；26.Cell；27.Cell part；28.Extracellular region；29.Intracellular part；30.Mitochondrion；31.Cytoplasm；32.Golgi apparatus part；33.Organelle；34.Intracellular；35.Intracellular organelle；36.Mitochondrial membrane part；37.RNA-directed DNA polymerase activity；38.DNA polymerase activity；39.Nucleotidyl transferase activity；40.RNA binding；41.Transferase activity；42.Catalytic activity；43.Nucleic acid binding；44.Beta-1，4-mannosylglycoprotein 4-beta-N-acetylglucosaminyltransferase activity；45.Transferase activity，transferring hexosyl groups；46.Acetylglucosaminyl transferase activity；47.Heterocyclic compound binding；48.Organic cyclic compound binding；49.UDP-glycosyltransferase activity；50.Transferase activity，transferring pentosyl groups；51.Transferase activity，transferring one-carbon groups；52.Binding；53.Receptor binding圖6　差異microRNAs靶基因的 GO 功能分類圖Fig.6　Gene Ontology annotation plot for target genes of microRNAs differentially expressed

表5　4個顯著富集的KEGG 代謝途徑

圖7　差異表達microRNAs靶基因與4個富集KEGG通路關系圖Fig.7　The relationship between target genes of differential expression of microRNAs and the 4 enriched KEGG pathways

3討論

繁殖性狀是綿羊的重要經(jīng)濟性狀，而繁殖性狀是受多基因控制的，遺傳力低，傳統(tǒng)的選擇方法對其遺傳進展的改良非常緩慢，培育新的高繁殖率品種是相當繁瑣甚至是不可能實現(xiàn)的。小尾寒羊是世界上具有非季節(jié)性發(fā)情和多胎特性的高繁殖率綿羊品種之一，研究其microRNA對卵泡發(fā)育和排卵的調(diào)控機制，開展分子標記輔助選擇，為解決這一問題提供了新的方向。

本試驗通過對小尾寒羊卵巢組織不同時期microRNA研究，豐富了綿羊microRNA數(shù)據(jù)庫的信息，發(fā)情期與間情期共有188種microRNAs表達：其中164種microRNAs在兩個時期共同表達，9種在間情期卵巢特異表達，15種在發(fā)情期卵巢中特異表達。間情期與發(fā)情期表達譜中表達量最多的microRNA分別是miR-99a、miR-143，除此之外一些表達豐富的microRNAs如miR-26a、let-7、miR-21，在人、牛[17]、豬[18]、成年和新生小鼠[19]的卵巢內(nèi)表達也是非常豐富的。這些保守的microRNAs在哺乳動物生殖階段高表達，說明它們在卵巢功能中扮演著非常重要的角色。在以往的研究中對這些高表達microRNA的部分功能進行了研究與分析：miR-21確定能夠促進排卵期卵泡細胞的存活[20]，而miR-143是小鼠原始卵泡形成的關鍵[21]以及l(fā)et-7b被證明是黃體正常發(fā)育過程中所必須的[22]。miR-99a的高表達結(jié)果和胎牛與小鼠卵巢microRNA克隆分析研究中miR-99a表達豐富相一致[19，23]。此外，KEGG通路對microRNA的靶基因富集，這對卵巢功能分析也很重要。例如，miR-143在促性腺激素釋放激素信號通路和孕激素介導的卵母細胞成熟途徑上顯著富集，這樣通過內(nèi)分泌功能、卵泡發(fā)育影響卵巢活動[11]。另外，let-7家族的靶基因參與細胞凋亡、Wnt信號通路和C21甾體激素代謝途徑，這對顆粒細胞增殖、卵泡的生長發(fā)育、黃體形成和退化很重要[24]。以上研究均證明了這些保守的高表達microRNAs在動物卵巢功能中的重要作用，同時也是因為這些作用才使得這些microRNAs在此次卵巢測序中高表達，與本試驗結(jié)果相對應。

本研究中，間情期和發(fā)情期卵巢差異表達的microRNAs(oar-miR-200a、oar-miR-200b、oar-miR-200c)均屬于miR-200家族。根據(jù)GO以及KEGG通路分析得到，差異表達microRNAs靶基因參與免疫與代謝過程與小尾寒羊卵巢周期性活動有關：細胞粘附分子通路，參與卵泡與黃體顆粒細胞的凋亡[25-26]；補體和凝血級聯(lián)反應則是因為排卵后卵泡破裂引起卵巢的一系列應答反應[27-28]；而脂肪酸代謝通路中靶基因富集，可能是由于兩個時期卵巢內(nèi)激素水平的變化對應的結(jié)果；糖代謝通路推測與兩個時期卵巢活動如卵泡發(fā)育、破裂所需的能量代謝有關。而以往對miR-200家族的功能研究多見于腫瘤方面，但是探索其對動物繁殖方面的作用也是很有意義的。ZEB1作為miR-200a、miR-200b、miR-200c共同的預測靶基因，在卵巢的表達模型中直接參與調(diào)節(jié)卵泡的作用，切除垂體的倉鼠卵巢內(nèi)ZEB1 mRNA顯著增加，而ZEB1蛋白卻顯著減少，導致mRNA和蛋白差異表達的原因正是由于miR-200轉(zhuǎn)錄后的調(diào)節(jié)作用[29]。H.Hasuwa等發(fā)現(xiàn)，缺乏miR-200b和miR-429的雌性小鼠不排卵，而miR-200b和miR-429通過抑制轉(zhuǎn)錄抑制因子ZEB1的表達調(diào)節(jié)小鼠繁殖過程[30]。在其他的報道中，miR-200b在小鼠的子宮、肌肉、腎中表達量較高，通過雙熒光素酶報告載體系統(tǒng)驗證其靶基因，推測miR-200b通過與靶基因Fog2的相互作用參與小鼠胚胎發(fā)育過程[31]。另外，miR-200a在懷孕6天和4天的小鼠子宮內(nèi)膜中相比較表達量顯著下調(diào)，表明miR-200a在胚胎著床前后對子宮內(nèi)膜起到重要的調(diào)節(jié)作用[32]。

4結(jié)論

通過成功地構(gòu)建小尾寒羊間情期和發(fā)情期microRNA的表達譜，為后續(xù)的綿羊microRNA研究提供更多更詳盡的序列信息。篩選出兩個時期差異表達的microRNAs，結(jié)合靶基因預測及通路富集分析結(jié)果，推測差異表達的microRNAs通過參與代謝與免疫過程調(diào)節(jié)卵泡的發(fā)育進而影響卵巢的周期性活動。

參考文獻(References)：

[1]HU S J，REN G，LIU J L，et al.MicroRNA expression and regulation in mouse uterus during embryo implantation[J].JBiolChem，2008，283(34)：23473-23484.

[2]CARLETTI M Z，CHRISTENSON L K.MicroRNA in the ovary and female reproductive tract[J].JAnimSci，2009，87(14)：29-38.

[3]徐盛玉，王定越，吳德.MicroRNA及其對哺乳動物繁殖的影響[J].畜牧獸醫(yī)學報，2011，42(6)：747-753.

XU S Y，WANG D Y，WU D.Effects of MicroRNA on the Mammalian Reproduction[J].ActaVeterinariaetZootechnicaSinica，2011，42(6)：747-753.(in Chinese)

[4]CHRISTENSON L K.MicroRNA control of ovarian function[J].AnimReprodSci， 2010，7(3)：129-133.

[5]SIROTKIN A V，OVCHARENKO D，GROSSMANN R，et al.Identification of microRNAs controlling human ovarian cell steroidogenesis via a genome-scale screen[J].JCellPhysiol，2009，219(2)：415-420.

[6]XU S，LINHER-MELVILLE K，YANG B B，et al.Micro-RNA378 (miR-378) regulates ovarian estradiol production by targeting aromatase[J].Endocrinology，2011，152(10)：3941-3951.

[7]YAN G J，ZHANG L，F(xiàn)ANG T，et al.MicroRNA-145 suppresses mouse granulosa cell proliferation by targeting activin receptor IB[J].FEBSLett，2012，586(19)：3263-3270.

[8]DAVIS B N，HILYARD A C，LAGNA G，et al.SMAD proteins control DROSHA-mediated microRNA maturation[J].Nature， 2008，454(7200)：56-61.

[9]VO N，KLEIN M E，VARLAMOVA O，et al.A cAMP-response element binding protein-induced microRNA regulates neuronal morphogenesis[J].ProcNatlAcadSciUSA，2005，102(45)：16426-16431.

[10]MURAMATSU F，KIDOYA H，NAITO H，et al.microRNA-125b inhibits tube formation of blood vessels through translational suppression of VE-cadherin[J].Oncogene，2013，32(4)：414-421.

[11]DI R，HE J，SONG S，et al.Characterization and comparative profiling of ovarian microRNAs during ovine anestrus and the breeding season[J].BMCGenomics，2014，15(1)：899.

[12]WEN M，SHEN Y，SHI S，et al.miREvo：An integrative microRNA evolutionary analysis platform for next-generation sequencing experiments[J].BMCBioinformatics，2012，13(1)：140.

[13]FRIEDLANDER M R，MACKOWIAK S D，LI N，et al.miRDeep2 accurately identifies known and hundreds of novel microRNA genes in seven animal clades[J].NucleicAcidsRes，2011，40(1)：37-52.

[14]MIAO X Y，LUO Q M，QIN X Y.Genome-wide transcriptome analysis of mRNAs and microRNAs in Dorset and Small Tail Han sheep to explore the regulation of fecundity[J].MolCellEndocrinol，2015，15(402)：32-42.

[15]BARTEL D P.MicroRNAs：target recognition and regulatory functions[J].Cell，2009，136(2)：215-233.

[16]VALEGGIA C，EILLSON P T.Interactions between metabolic and reproductive functions in the resumption of postpartum fecundity[J].AmJHumBiol，2009，21(4)：559-566.

[17]HOSSAIN M M，GHANEM N，HOELKER M，et al.Identification and characterization of miRNAs expressed in the bovine ovary[J].BMCGenomics，2009，10(1)：443.

[18]LI M Z，LIU Y K，WANG T，et al.Repertoire of porcine microRNAs in adult ovary and testis by deep sequencing[J].IntJBiolSci， 2011，7(7)：1045-1055.

[19]MISHIMA T，TAKIZAWA T，LUO S S，et al.MicroRNA (miRNA) cloning analysis reveals sex differences in miRNA expression profiles[J].Reproduction，2008，136(6)：811-822.

[20]CARLETTI M Z，F(xiàn)IEDLER S D，CHRISTENSON L K.MicroRNA 21 blocks apoptosis in mouse periovulatory granulosa cells[J].BiolReprod，2010，83(2)：286-295.

[21]ZHANG J F，JI X W，ZHOU D D，et al.miR-143 is critical for the formation of primordial follicles in mice[J].FrontBiosci，2012(1)，18：588-597.

[22]OTSUKA M，ZHENG M，HAYASHI M，et al.Impaired microRNA processing causes corpus luteum insufficiency and infertility in mice[J].JClinInvest，2008，118(5)：1944-1954.

[23]TRIPURANI S K，XIAO C，SALEM M，et al.Cloning and analysis of fetal ovary microRNAs in cattle[J].AnimReprodSci，2010，120(1)：16-22.

[24]ROUSH S，SLACK F J.The let-7 family of microRNAs[J].TrendsCellBiol，2008，18(10)：505-516.

[25]MAKRIGIANNAKIS A，COUKOS G，CHRISTOFIDOU-SOLOMIDOU M，et al.N-cadherin-mediated human granulosa cell adhesion prevents apoptosis：a role in follicular atresia and luteolysis?[J].AmJPathol，1999，154(5)：1391-1406.

[26]MORA J M，F(xiàn)ENWICK M A，CASTLE L，et al.Characterization and significance of adhesion and junction-related proteins in mouse ovarian follicles[J].BiolReprod，2012，86(5)：153.

[27]AMBEKER A S，NLIRUJOG R S，SRIKANTH S M，et al.Proteomic analysis of human follicular fluid：a new perspective towards understanding folliculogenesis[J].JProteomics，2013，87(7)：68-77.

[28]EBISCH I M W，THOMAS C M G，WETZELS A M M，et al.Review of the role of the plasminogen activator system and vascular endothelial growth factor in subfertility[J].FertilSteril，2008，90(6)：2340-2350.

[29]WANG C，ROY S K.Gonadotropins and steroids regulation of ZEB1 expression in the hamster ovary[J].BiolReprod，2011，85(1 Supplement)：661.

[30]HASUWA H，UEDA J，IKAWA M，et al.MiR-200b and miR-429 Function in mouse ovulation and are essential for female fertility[J].Science，2013，341(6141)：71-73.

[31]任紅艷，劉楠，陶聰，等.小鼠體型控制相關microRNA-200b(miR-200b)的表達譜及靶標分析[J].農(nóng)業(yè)生物技術(shù)，2013，21(1)：47-54.

REN H Y，LIU N，TAO C，et al.The expression profile and target gene analysis of mouse(Musmusculus) body size REgulation related microRNA-200b(miR-200b)[J].JournalofAgriculturalBiotechnology，2013，21(1)：47-54.(in Chinese)

[32]魯之中，陳雪梅，曾蘭，等.早孕小鼠胚胎著床前后子宮內(nèi)膜microRNAs差異表達[J].細胞生物學，2009，31(6)：848-852.

LU Z Z，CHEN X M，ZENG L，et al.The different expression of microRNAs in endometrium of early pregnant mouse during embryo implantation[J].ChineseJournalofCellBiology，2009，31(6)：848-852.(in Chinese)

(編輯郭云雁)

doi:10.11843/j.issn.0366-6964.2016.07.004

收稿日期：2015-11-13

基金項目：河北省高等學?？茖W技術(shù)研究青年基金項目(QN2015162)；河北省首批青年拔尖人才支持計劃(2013-2015)；河北農(nóng)業(yè)大學青年學術(shù)帶頭人支持項目(2015-2017)

作者簡介：段新崇(1989-)，女，河北邯鄲人，碩士生，主要從事動物遺傳育種與繁殖研究，E-mail: 15733202716@163.com *通信作者：周榮艷，副教授，博士，主要從事動物遺傳育種研究，Tel: 0312-7528451，E-mail: rongyanzhou@126.com

中圖分類號：S826；S813.3

文獻標志碼：A

文章編號:0366-6964(2016)07-1324-09

Analyzing the Differential Expression of microRNAs in Estrus and Diestrus of Small Tailed Han Sheep

DUAN Xin-chong，WEI Yan-hui，LI Yang，LI Xiang-yun，ZHOU Rong-yan*，XI Jian-zhong

(CollegeofAnimalScienceandTechnology，HebeiAgriculturalUniversity，Baoding071001，China)

Abstract:To explore the theoretical basis of the microRNAs in regulating the reproductive process in Small Tail Han sheep，the ovarian microRNAs profiles were constructed and analyzed in estrus and diestrus.One ovary was collected with surgery method at estrus (the 1st day) and diestrus (the 13th day) from the same sheep，respectively.The RNA sequencing data were obtained with illumina Hiseq2000 sequencing platform.The expression profiles and differential expression of microRNAs were analyzed by bioinformatics method.The real-time quantitative PCR (qRT-PCR) was used for validating the differentially expressed microRNAs between the 2 periods in ovary.The microRNAs expression profiles of Small Tailed Han sheep was constructed in the period of estrus and diestrus in ovaries.The highest expression of microRNA was observed for oar-miR-99a and oar-miR-143 in the period of diestrus and estrus，respectively.Three microRNAs，oar-miR-200a，oar-miR-200b and oar-miR-200c，were significantly differentially expressed between estrus and diestrus.The expression level of 2 randomly selected microRNAs by qRT-PCR was the same as the result by RNA-seq data.These microRNAs expression profiles will be used to provide informations for further studying the microRNAs in sheep ovary.The differentially expressed microRNAs may regulate the ovarian periodic activity through the metabolism and immunity pathway by the analysis of predicted target genes and enriched pathways.

Key words:Small Tailed Han sheep；diestrus；estrus；ovary；microRNA