miRNA-34基因家族進化分析及其在鵝生殖周期內的表達變化

劉　然，張惠子，云皓琪，陳　陽，2，黃正洋，2，張　揚，2，徐　琪，2，陳國宏，2*

(1.揚州大學 動物科學與技術學院，揚州 225009； 2.江蘇省動物遺傳繁育與分子設計重點實驗室，揚州 225009)

?

miRNA-34基因家族進化分析及其在鵝生殖周期內的表達變化

劉然1，2#，張惠子1#，云皓琪1，陳陽1，2，黃正洋1，2，張揚1，2，徐琪1，2，陳國宏1，2*

(1.揚州大學 動物科學與技術學院，揚州 225009； 2.江蘇省動物遺傳繁育與分子設計重點實驗室，揚州 225009)

摘要：旨在探明miRNA-34基因家族(miRNA-34s)系統進化歷程及其在鵝生殖周期內的表達變化規律。通過文獻和miRBase數據庫檢索脊椎動物的miRNA-34s序列，利用Ensembl數據庫信息，確定miRNA-34s在基因組中的位置，采用MEGA 6.0構建miRNA-34s的系統進化樹；并利用RT-qPCR檢測鵝卵巢組織miRNA-34s在生殖周期內的表達變化。結果在脊椎動物中共搜索70余條miRNA-34s序列，miRNA-34s 3個成員(miRNA-34a、miRNA-34b和miRNA-34c)在脊椎動物中分布廣泛，高度保守，均位于基因間隔區(Intergenic region，IGR )，且miRNA-34b和miRNA-34c兩個基因位于同一染色體上，在進化過程中緊密相連，可能是由于串聯重復所致。RT-qPCR檢測結果顯示，miRNA-34s 3個成員在鵝的開產前期、排卵期、排卵后期和就巢期卵巢組織中表現出相同趨勢的表達，即在排卵期、排卵后期表達水平高于或顯著高于開產前期和就巢期。miRNA-34家族在進化中高度保守，且在調控鵝生殖周期中發揮一定的作用。

關鍵詞：鵝；miRNA-34；進化；定量表達

microRNA(miRNA)作為廣泛影響基因表達的潛在調控因子，在真核細胞中具有高度的保守性。miRNA在細胞生長、維持內環境穩定及病理過程中起重要作用。miRNA主要通過結合目的mRNA的3′非編碼區域的互補序列來發揮作用，并通過加速轉錄的進程來抑制基因表達，或通過與多蛋白RISC的相互作用來抑制翻譯進程[1]。不同物種間的同源miRNA存在高度進化保守性，如miRNA-l、miRNA-34、miRNA-87在無脊椎動物和脊椎動物中高度保守，序列比對分析發現這些保守序列的堿基差異僅為1～2 nt[2]。miRNA在物種間的保守性使其成為探索非編碼RNA進化規律的理想分子。

miRNA-34家族包括3種miRNA(a/b/c)，它們包含相同的種子區域，具有高度保守性。這種保守性與其功能有著密切關系。研究顯示，miRNA-34家族成員在配子形成和早期胚胎發育過程中發揮重要的作用[3-5]。A.Tscherner等[6]研究發現miRNA-34s主要作用于牛配子形成時期，且在精子中miRNA-34s有獨特的表達特征，可以作為雄性的生產力指標。有研究發現miRNA-34a和miRNA-34b在多種組織中廣泛表達[7-8]，但miRNA-34c僅在生殖腺等組織中檢測到[9]。王勃[10]通過RT-qPCR技術比較miRNA-34c在牛胎兒成纖維細胞、精子、MII期卵母細胞中表達，確定miRNA-34c在精子中表達量較高，而在體細胞、卵母細胞中的相對表達量非常低。精子中所攜載的miRNA-34c可以通過誘導表達的形式賦予供體細胞，能提高胚胎發育能力、核重編程及胚胎質量。精子miRNA-34c表達水平可能會影響卵裂期的胚胎質量，對卵裂期的胚胎發育起到一定的積極作用[11]。本課題組前期以浙東白鵝為研究對象，采用Solexa深度測序產蛋和就巢組鵝卵巢miRNA的轉錄組，用生物信息學來確定miRNA的差異表達。鑒定1 328條已知的保守miRNA和22個新的潛在候選miRNA，發現miRNA-34家族在產蛋鵝和就巢鵝之間表達差異達到顯著水平[12]。綜上表明，miRNA-34家族與生殖調控關系密切，且與配子生成有關。

鑒于miRNA-34s與生殖調控及配子生成密切相關，在前期研究的基礎上，為進一步探明miRNA-34s進化分析及其在鵝生殖周期內的表達變化，本研究利用Ensembl數據庫信息，比對幾種常見禽類miRNA-34s前體序列，確定miRNA-34s在基因組中的位置，構建miRNA-34s系統進化樹；并利用RT-qPCR檢測鵝卵巢組織miRNA-34s的表達變化，旨在為進一步揭示鵝miRNA-34s在生殖調控中作用奠定基礎。

1材料與方法

1.1miRNA-34序列的獲得

從miRBase(http：//www.mirbase.org/cgi-bin/browse.pl) 中下載獲得25種脊椎動物的miRNA-34基因成熟序列及前體序列，并從已報道文獻中獲得鴨和鵝的miRNA-34基因成熟序列[12-13]，其前體序列從鴨、鵝genome assembly scaffolds BLAST獲得。

1.2幾種常見禽類的miRNA-34s基因定位

采用Ensembl數據庫BLAST(http：//asia.ensembl.org/index.html)程序，分別用家雞、家鴨、家鵝、斑胸草雀等4個禽類(以人、小鼠)的miRNA-34前體序列對其基因組數據庫進行搜索，以確定miRNA-34s在基因組中的位置。

1.3miRNA-34s系統發生分析

選擇21種脊椎動物和無脊椎動物中每個種屬中代表性物種，采用Clustal W軟件對其miRNA-34基因的前體序列進行多序列比對分析(Multiple sequences alignment，MSA)，利用MEGA 6.0軟件采用基于距離參數的鄰接法(Neighbor-joining，NJ)，并自展分析(Bootstrap) 1 000次，構建miRNA-34s的系統進化樹。

1.4miRNA-34s在鵝生殖周期內的定量表達檢測

鵝的生殖周期分段明顯，包含產蛋前期、產蛋期(就巢期)和休產期等幾個階段。鵝150日齡開產，產蛋一段時間后，體溫升高，被毛蓬松，抱蛋而窩，停止產蛋即進入就巢期。本試驗采用miRcute miRNA Isolation Kit (TIANGEN，China，DP501)分別提取開產前期(120日齡)、產蛋期(380日齡，排卵前)、產蛋期(380日齡，排卵后)和就巢期(380日齡)浙東白鵝母鵝卵巢基質的sRNA，取3 μL sRNA，采用Poly(A)加尾法進行RT-PCR，具體操作按miRcute miRNA First-Strand cDNA Synthesis kit (Tiangen)試劑盒說明書進行。采用SYBR Green qPCR 法檢測miRNA-34a、miRNA-34b和miRNA-34c的表達量，擴增引物見表1，PCR擴增采用兩步法，PCR反應條件參照miRcute miRNA qPCR detection kit試劑盒說明書進行，PCR反應在Applied Biosystems 7500 熒光定量 PCR系統中進行，U6作為內參，每個樣本重復3次，采用2-ΔΔCt計算miRNA相對表達量(以開產前期卵巢組織為對照)。各時期表達量差異采用獨立樣本t檢驗(Independent-samplesttest)進行統計分析。

表1miRNA-34s RT-qPCR擴增引物

Table 1Sequence of the oligonucleotide primers of RT-qPCR

2結果

2.1miRNA-34s及其在基因組中的定位

根據文獻[12-14]和miRBase檢索共搜索27種脊椎動物70余條miRNA-34s序列，其中人類(Homosapiens)、小鼠(Musmusculus)和雞(Gallusgallus)等18種脊椎動物都具有miRNA-34s 3個成員(miRNA-34a、miRNA-34b和miRNA-34c)；絨毛猴(Lagothrixlagotricha)、蜘蛛猴(Atelesgeoffroyi)和平頂猴(Macacanemestrina)僅檢索到miRNA-34a；樹鼩(Tupaiachinensis)和熱帶爪蟾(Xenopustropicalis)擁有miRNA-34a和miRNA-34b，且發現熱帶爪蟾存在xtr-miRNA-34b、xtr-miRNA-34b-1、xtr-miRNA-34b-2、xtr-miRNA-34b-3和xtr-miRNA-34b-4 同源序列；鴨嘴獸(Ornithorhynchusanatinus)和豬(Susscrofa)擁有miRNA-34a和miRNA-34c。

以家雞、家鴨、家鵝、斑胸草雀4個禽類(以人、小鼠)的miRNA-34前體作為查詢序列(Query sequence)，分別對其相關基因組進行BLAST 搜索，確定各序列在基因組中的位置(表2)，基因定位顯示，與其它哺乳動物miRNA主要定位在內含子區域不同的是，禽類miRNA-34s家族3個成員均位于基因間隔區(Intergenic region，IGR )，且miRNA-34a和miRNA-34b兩個基因位于同一染色體上。

2.2miRNA-34s序列分析

對14種典型物種的miRNA-34s的成熟序列進行多序列比對分析(圖1)，結果顯示，miRNA-34s的成熟序列同源性很高，序列長度(22～23 nt)相近，保守堿基17個，分別為第3～10位、14～20位以及22～24位堿基，說明miRNA-34s在不同物種間的高度保守性。比對結果顯示，miRNA-34s在其進化過程中可能發生了缺失，有的物種在第1位、25位和26位缺失了堿基U/C，大多數物種在第13位缺失了堿基A/G；在第12和21位堿基分別出現A>C和A>G。

2.3miRNA-34s系統進化樹的構建

來源于miRBase中miRNA-34基因主要分布在脊椎動物哺乳綱、鳥綱、兩棲綱以及無脊椎動物的昆蟲綱和線蟲綱等，采用MEGA 6.0軟件，以這些綱目中的代表性物種miRNA-34s序列構建基于距離參數的鄰接法系統發育進化樹(圖2)。由圖2可以看出，進化樹大致可分為兩支：miRNA-34b和 miRNA-34c位于第一分支中，這一分支中只有一個例外是spu-miRNA-34，另一分支為 miRNA-34/miRNA-34a。從進化樹的分支和距離上看，miRNA-34a基因的分支出現在進化早期，并且一直延續到現在，因而它很可能是miRNA-34 基因家族中最早出現的基因形式，即祖先基因。在無脊椎動物中，miRNA-34 基因只有一個拷貝(miRNA-34)，而在脊椎動物中 miRNA-34 基因一般都有3種形式 (miRNA-34a、miRNA-34b和miRNA-34c)，且miRNA-34基因家族的3個成員穩定地存在于從魚類到哺乳動物中。

2.4miRNA-34s在鵝生殖周期內的定量表達檢測

采用RT-qPCR對開產前期、產蛋期(排卵)、產蛋期(排卵后)和就巢期鵝卵巢基質miRNA-34家族成員進行檢測，結果見圖3。由圖3可以看出，miRNA-34a、miRNA-34b和miRNA-34c在鵝卵巢組織中均有表達，且開產前期、排卵期、排卵后期和就巢期呈現一定規律表達，即在排卵期、排卵后期表達水平高于或顯著高于開產前期和就巢期。

表2miR-34s在基因組中的定位

Table 2Genomic location of miR-34 gene family

3討論

microRNA作為重要的內源性調控基因，分布范圍廣，參與多種生物學調控過程[15]。miRNA-34s在脊椎動物中分布廣泛，且高度保守，但在無脊椎動物和脊椎動物中miRNA-34s成員數并不相同，在漫長的進化過程中，3個成員組成的miRNA-34 基因家族保守地存在于幾乎所有的脊椎動物中，這是否意味著與無脊椎動物相比，脊椎動物需要更多的miRNA-34 基因來形成一個基因調控網絡，從而成功的行使其調控功能。對miRNA-34s的進化分析顯示，miRNA-34a基因主要出現在進化早期，為祖先基因，且miRNA-34b和mir-34c基因親緣關系很近，且位于同一染色體上聚集成簇，暗示它們在進化過程中緊密相連，這可能在miRNA-34基因進化過程中出現了串聯重復所致。J.Hertel等[16]在后生動物基因組miRNA基因家族之間的系統發育分析，揭示了串聯重復和片段重復是miRNA基因簇進化的主要機制，也因此產生miRNA基因簇和miRNA基因家族的多樣性分布。作為轉錄后水平調控元件，miRNA成熟序列5′端 的第2～8位堿基對于其功能的正常發揮至關重要，被稱為種子序列(Seed sequence)[2]。Blast結果顯示，不同物種的miRNA-34s種子序列保守性為100%，暗示其功能的保守性。本研究中也發現miRNA-34s成員在鵝的開產前期、排卵期、排卵后期和就巢期表現出相同趨勢的表達。miRNA-34家族成員在不同生殖時期的表達變化較大。這提示miRNA-34家族在調控鵝生殖周期中發揮重要作用。有報道miRNA-34家族的成員 (miRNA-34a/b/c)是p53作用的直接靶點[15]。miRNA-34s在小鼠睪丸等組織的高表達在很大程度上依賴p53在這些組織中的表達，并在p53通路中的存在其他作用。在鵝的生殖周期調控中，受生殖激素水平的制約，如孕酮(P4)，開產前P4慢慢上升，并在產蛋高峰期維持高水平，就巢期最低，且排卵周期內P4在產蛋前13 h內分別出現兩次高峰[17]。這一結果與本研究miRNA-34s在在不同生殖時期的表達變化規律一致。孕酮能上調卵巢細胞的p53表達[18]；筆者推測miRNA-34家族成員在不同生殖時期的表達變化仍依賴p53表達。

4結論

4.1miRNA-34家族3個成員(miRNA-34a、miRNA-34b和miRNA-34c)在脊椎動物中分布廣泛，且高度保守；均位于基因間隔區(Intergenic region，IGR )，且miRNA-34b和miRNA-34c兩個基因位于同一染色體上，在進化過程中緊密相連，可能是由于串聯重復所致。

4.2RT-qPCR檢測結果顯示，miRNA-34s成員(miRNA-34a、miRNA-34b和miRNA-34c)在鵝的開產前期、排卵期、排卵后期和就巢期卵巢組織中表現出相同趨勢的表達，即在排卵期、排卵后期表達水平高于或顯著高于開產前期和就巢期。

參考文獻(References)：

[1]GREGORY R I，CHENDRIMADA T P，COOCH N，et al.Human RISC couples microRNA biogenesis and post transcriptional gene silencing[J].Cell，2005，123(4)：631-640.

[2]羅丹丹，羅玉萍，彭娟，等.miRNA-9基因家族進化分析及其靶基因預測[J].中國細胞生物學學報，2011，33(5)：510-519.

LUO D D，LUO Y P，PENG J，et al.Molecular evolution of miRNA-9 gene family and prediction of their target genes[J].ChineseJournalofCellBiology，2011，33(5)：510-519.(in Chinese)

[3]BOUHALLIER F，ALLIOLI N，LAVIAL F，et al.Role of miR-34c microRNA in the late steps of spermatogenesis[J].RNA，2010，16(4)：720-731.

[4]LIN W M，PANG R T，CHIU P C，et al.Sperm-borne microRNA-34c is required for the first cleavage division in mouse[J].ProcNatlAcadSciUSA， 2012，109(2)：490-494.

[5]SONI K，CHOUDHARY A，PATOWARY A，et al.miR-34 is maternally inherited inDrosophilamelanogasterandDaniorerio[J].NucleicAcidsRes， 2013，41(8)：4470-4480.

[6]TSCHERNER A，GILCHRIST G，SMITH N，et al.MicroRNA-34 family expression in bovine gametes and preimplantation embryos[J].ReprodBiolEndocrinol， 2014，12：85.

[7]楊蕊珂，楊麗紅，王淑輝，等.微小核糖核酸-34a對大鼠肥厚心肌組織肌球蛋白重鏈α基因表達的調控[J].中國循環雜志，2014，29(2)：138-141.

YANG R K，YANG L H，WANG S H，et al.Regulation of microRNA-34a on myosin heavy chain-α expression in rats with left ventricular hypertrophic tissue[J].ChineseCirculationJournal，2014，29(2)：138-141.(in Chinese)

[8]鐘金發，湯愛萍.DNA甲基化與miR-34b表達及生物學功能的關系[J].腫瘤防治研究，2013，40(11)：1105-1108.

ZHONG J F，TANG A P.Biological relations between DNA methylation and miR-34b expression[J].CancerResearchJournal，2013，40(11)：1105-1108.(in Chinese)

[9]DUTTA K K，ZHONG Y，LIU Y T，et al.Association of microRNA-34a overexpression with proliferation is cell type-dependent[J].CancerSci，2007，98(12)：1845-1852.

[10]王勃.miR-34c和miR-449b對牛體細胞克隆胚早期發育的影響[D].楊凌：西北農林科技大學，2014.

WANG B.Affects of miR-34c and miR-449b on the development of early bovine somatic cell nuclear transfer embryos[D]Yangling：Northwest A&F University，2014.(in Chinese)

[11]江雯.IVF患者精子miR-21，miR-34c，miR-140，miR-375與卵裂期胚胎質量的相關性研究[D].長沙：中南大學，2013.

JIANG W.The correlation research of the IVF patients’ miR-21，miR-34c，miR-140 and miR-375 expression in the sperm and cleavage stage embryo quality[D].Changsha：Central South University，2013.(in Chinese)

[12]XU Q，ZHANG Y，CHEN Y，et al.Identification and differential expression of microRNAs in ovaries of laying and Broody geese (Ansercygnoides) by Solexa sequencing[J].PLoSOne，2014，9(2)：e87920.

[13]GU L，XU T，HUANG W，et al.Identification and profiling of microRNAs in the embryonic breast muscle of pekin duck[J].PLoSOne，2014，9(1)：e86150.

[14]ZHANG L，NIE Q，SU Y，et al.MicroRNA profile analysis on duck feather follicle and skin with high-throughput sequencing technology[J].Gene，2013，519(1)：77-81.

[15]高佳莉，羅玉萍，李思光.miR-34基因家族的分子進化[J].動物學研究，2007，28(3)：271-278.

GAO J L，LUO Y P，LI S G.Molecular Evolution of miR-34 gene family[J].ZoologicalResearch，2007，28(3)：271-278.(in Chinese)

[16]HERTEL J，LINDEMEYER M，MISSAL K，et al.The expansion of the metazoan microRNA repertoire[J].BMCGenomics，2006，7：25.

[17]李定健.鵝發育期和繁殖周期生殖激素變化的研究進展[J].湖南農業科學，2012(22)：26-28.

LI D J.Progress on reproductive hormones during the development period and the reproductive cycle of geese[J].HunanAgriculturalSciences，2012(22)：26-28.(in Chinese)

[18]BU S Z，YIN D L，REN X H，et al.Progesterone induces apoptosis and up-regulation of p53 expression in human ovarian carcinoma cell lines[J].Cancer，1997，79(10)：1944-1950.

(編輯程金華)

Molecular Evolution Analyzing of miRNA-34 Gene Family and Its Expression in Goose Reprodutive Cycle

LIU Ran1，2#，ZHANG Hui-zi1#，YUN Hao-qi1，CHEN Yang1，2，HUANG Zheng-yang1，2，ZHANG Yang1，2，XU Qi1，2，CHEN Guo-hong1，2*

(1.CollegeofAnimalScienceandTechnology，YangzhouUniversity，Yangzhou225009，China；2.JiangsuKeyLaboratoryofAnimalGenetics&BreedingandMolecularDesign，Yangzhou225009，China)

Key words:goose；miRNA-34；evolution；quantitative expression

Abstract:This experiment was conducted to investigate the miRNA-34 gene family evolutionary process and its expression regularity in goose reproductive cycle.The miRNA-34 genes were searched in literatures and miRBase.It’s genes were located by Ensemble database，and its phylogenetic tree was constructed by MEGA 6.0.In addition，the expression level of miRNA-34s was detected in ovarian tissue of goose reproductive cycle using RT-qPCR.A total 70 sequences in 27 species which indicated the extensive existence of miRNA-34s(miRNA-34a，miRNA-34b and miRNA-34c) in different species.The analysis of gene localization showed all miRNA-34 genes were located on the intergenic region(IGR).miRNA-34b and miRNA-34c shared on the same chromosome and they were closely linked in evolutionary process，which maybe caused by tandem repeats.Besides，the similar expression of miRNA-34s were detected in goose ovary tissues of pre-laying stage，ovulation，oviposition，and the broody phase.More miRNA-34s expression was detected at ovulation，oviposition than that in pre-laying and broody phase period.The miRNA-34s are highly conservative in evolution process，and play a role in regulation reproductive activity in goose.

doi:10.11843/j.issn.0366-6964.2016.02.007

收稿日期：2015-06-10

基金項目：江蘇省高等學校大學生創新創業訓練計劃項目(201411117050Y)；現代農業產業技術體系建設專項CARS-43-3)；江蘇高校優勢學科建設工程資助項目(蘇政辦發[2011]137號)

作者簡介：劉然(1990-)，女，山東濟寧人，碩士生，主要從事動物遺傳育種與繁殖方面研究，E-mail:betterlr@163.com；張惠子(1993-)，女，江蘇連云港人，本科生，主要從事動物科學研究，E-mail：15396752381@163.com。劉然和張惠子為并列第1作者 *通信作者：陳國宏，E-mail：ghchen@yzu.edu.cn

中圖分類號：S835.2

文獻標志碼：A

文章編號:0366-6964(2016)02-0260-08