Integrated analysis of miRNA–mRNA interaction in ovaries of Turpan Black Sheep during follicular and luteal phases

To investigate the regulatory mechanism of the follicular-luteal phase transition in Turpan black sheep (Ovis aries), the genome-wide expression patterns of microRNAs (miRNAs) and genes were investigated in ovaries of six sheep (3 years and single lamb with 3 consecutive births) during follicular and luteal phases of the oestrous cycle. Bioinformatic analysis was used to screen potential miRNAs and genes related to Turpan black sheep ovarian function. RT-qPCR was used to validate the sequencing results. In total, we identified 139 known and 71 novel miRNAs in the two phases with miRNA-seq, and a total of 19 miRNAs were significantly differentially expressed, of which 7 were up-regulated and 12 were down-regulated in the follicular phase compared with luteal phase. A total of 150 genes were significantly differentially expressed, including 63 up-regulated and 87 down-regulated in the follicular phase compared with the luteal phase by RNA-seq data analysis. Those DEGs were significantly enriched in 103 GO terms and several KEGG pathways, including metabolic pathway, ovarian steroidogenesis, steroid hormone biosynthesis and oestrogen signalling pathway. In addition, we created a miRNA-mRNA regulatory network to further elucidate the mechanism of follicular-luteal transition. Finally, we identified key miRNAs and genes including miR-143, miR-99a, miR-150, miR-27a, miR-125b, STAR, STAT1, which might play crucial roles in reproductive hormone biosynthesis and follicular development. The miRNA-mRNA interactive network clearly illustrates molecular basis involving in follicular-luteal transition.