Abstract
Circular RNA (circRNA), as an emerging class of noncoding RNA, has been found to play key roles in many biological processes. However, its expression profile in the hypothalamus, a powerful organ initiating the reproductive process, has not yet been explored. Therefore, we used RNA sequencing to explore the expression of circRNAs in the hypothalamus of sheep with the FecB ++ genotype. We totally identified 41,863 circRNAs from sheep hypothalamus, in which 333 (162 were upregulated, while 171 were downregulated) were differentially expressed in polytocous sheep in the follicular phase versus monotocous sheep in the follicular phase (PF vs. MF), moreover, 340 circRNAs (163 were upregulated, while 177 were downregulated) were differentially expressed in polytocous sheep in the luteal phase versus monotocous sheep in the luteal sheep (PL vs. ML). We also identified several key circRNAs including oar_circ_0018794, oar_circ_0008291, oar_circ_0015119, oar_circ_0012801, oar_circ_0010234, and oar_circ_0013788 through functional enrichment analysis and oar_circ_0012110 through a competing endogenous RNA network, most of which may participate in reproduction by influencing gonadotropin-releasing hormone (GnRH) activities or affecting key gene expression, indirectly or directly. Our study explored the overall expression profile of circRNAs in sheep hypothalamus, which potentially provides an alternative insight into the mechanism of sheep prolificacy without the effects of FecB mutation.