The Potential Biological Functions of Circular RNAs During the Initiation of Atresia in Pig Follicles

The specific expression profile and function of circular RNAs (circRNAs) in mammalian ovarian follicles, especially during the atresia process, are unclear. In this study, genome-wide deep circRNA sequencing was applied to screen circRNAs in healthy and early atretic antral follicles in pig ovaries. A total of 40,567 distinct circRNAs were identified in follicles, among which 197 circRNAs (108 upregulated and 89 downregulated) were significantly shifted during the early atresia process. Most differentially expressed circRNAs (DECs) lacked protein-coding potential. Annotation analysis of the DECs revealed 162 known host genes, or noncoding RNAs, and 10 intergenic regions. The key pathways in which these host genes are involved include the focal adhesion-PI3K-Akt-mTOR signaling pathway, vascular endothelial growth factor A (VEGFA)-vascular endothelial growth factor receptor 2 signaling pathway and transforming growth factor-beta signaling pathway. Further comparison analysis between host genes of DECs and the differentially expressed linear messenger RNA transcripts revealed the cotranscription of circRNAs and their linear mRNAs in inhibin beta units (INHBA and INHBB), glutathione S-transferase (GSTA1), and VEGFA. In addition, we predicted 196 pairs of potential circRNA-micro RNA (miRNA) interactions among 77 DECs and 101 porcine miRNAs. We have identified 16 functional miRNAs by comparing the 101 miRNAs to the functional miRNAs reported in mammal ovarian follicle atresia and granulosa cell apoptosis studies. Our study adds new knowledge to circRNA distribution profiles in pig ovarian follicles, offers a valuable reference for transcriptomic profiles in the initiation of follicular atresia, highlights warranted circRNAs for further functional investigation, and provides possible biomarkers for ovarian dysfunctions.