畜牧兽医学报 ›› 2017, Vol. 48 ›› Issue (4): 627-636.doi: 10.11843/j.issn.0366-6964.2017.04.005

• 生物技术与繁殖 • 上一篇    下一篇

营养诱导舍饲哈萨克羊非繁殖季节卵巢组织中miRNAs差异表达分析

雷小萍1, 林杉2, 于要升1, 梁慧慧1, 赵宗胜1*   

  1. 1. 石河子大学动物科技学院, 石河子 832003;
    2. 石河子大学生命科学学院, 石河子 832003
  • 收稿日期:2016-10-25 出版日期:2017-04-23 发布日期:2017-04-21
  • 通讯作者: 赵宗胜,教授,主要从事动物分子遗传研究,E-mail:zhaozongsh@shzu.edu.cn
  • 作者简介:雷小萍(1991-),女,甘肃通渭人,硕士生,主要从事动物遗传育种研究,E-mail:leixiaoping0329@163.com
  • 基金资助:

    国家自然科学基金项目(31260537)

Analyzing the Differential Expression of miRNAs in Ovary of Stall-feed Kazakh Sheep Induced by Nutrition in Non-breeding Season

LEI Xiao-ping1, LIN Shan2, YU Yao-sheng1, LIANG Hui-hui1, ZHAO Zong-sheng1*   

  1. 1. College of Animal Science and Technology, Shihezi University, Shihezi 832003, China;
    2. College of Life Sciences, Shihezi University, Shihezi 832003, China
  • Received:2016-10-25 Online:2017-04-23 Published:2017-04-21

摘要:

旨在利用Solexa测序技术分析、挖掘营养诱导非繁殖季节绵羊发情期和乏情期卵巢组织micro RNAs(miRNA),为进一步研究特定miRNA参与非繁殖季节绵羊发情的分子调控机制提供理论基础。以不同营养水平饲喂的哈萨克绵羊作为研究对象,在非繁殖季节前后分别屠宰已鉴定的3只发情(EN)和乏情(AN)母羊,分别采集卵巢(O)组织,构建两个小片段RNA文库(OEN和OAN),进行高通量测序和生物信息学分析,实时荧光定量PCR(qRT-PCR)验证差异表达的miRNAs在哈萨克羊卵巢中表达水平。结果显示,两文库分别获得8 665 978(发情期)和9 071 102条(乏情期)clean reads,且绝大多数小RNA序列长度为21~23 nt。两文库相比较,发现有9个已知和104个未知miRNAs表达差异显著(P<0.01)。利用qRT-PCR随机选择的6个显著差异表达的miRNAs进行验证,其表达水平和测序分析结果一致。KEGG通路分析表明,这些差异miRNAs的靶基因主要参与了甘氨酸、丝氨酸和苏氨酸代谢,氨基糖和核苷酸糖代谢,精氨酸和鸟氨酸代谢及甘油磷脂代谢等营养调控通路。结合靶基因预测及通路富集分析,推测差异表达的miRNAs是通过营养代谢通路对绵羊的非季节性发情起重要作用。

关键词: 哈萨克羊, 发情期, 乏情期, 卵巢, Solexa测序, microRNA

Abstract:

To further explore the theoretical basis of the microRNAs in regulating the non-breeding season estrous in sheep, characterization of ovarian microRNAs were constructed and analyzed in estrous and anestrous of non-breeding season induced by nutrition by Solexa sequencing. In this study, the Kazakh sheep with different nutrition levels were used as the research object. In non-breeding season,3 estrous (EN) and 3 anestrous (AN) sheep were slaughtered with their ovary tissue(O) collected, respectively. Two miRNA libraries were constructed and analyzed by the bioinformatics methods. The real-time quantitative PCR (qRT-PCR) was used for validating the differentially expressed miRNAs betweent the 2 periods in ovary. The results showed that the clean reads of 8 665 978 (OEN) and 9 071 102 (OAN) were obtained in the two libraries, and most of the small RNA sequences were between 21 and 23 nt. Compared the two ovary (OEN vs OAN) libaries, the 9 known and 104 novel miRNAs had significant expression difference (P<0.01). The expression level of 6 randomly selected miRNAs by qRT-PCR was the same as the result by Solexa sequencing. KEGG pathway analysis showed that these target genes of difference miRNAs were primarily involved in nutrition regulatory related pathways, such as glycine, serine and threonine metabolism, amino sugar and nucleotide sugar metabolism, D-arginine and D-ornithine metabolism, glycerol phospholipid metabolism, and so on. The differentially expressed miRNAs may play an important role in non-breeding season estrous of sheep through nutrition metabolism pathway by the analysis of predicted target and enriched pathway analysis.

Key words: Kazakh sheep, estrus, anestrus, ovary, solexa sequencing, microRNA

中图分类号: