Abstract
Ovarian follicles undergo a series of dynamic changes following the ovulatory surge of luteinizing hormone including cumulus expansion, oocyte maturation, ovulation, and luteinization. Post-transcriptional gene regulatory events are critical for mediating LH follicular responses, and among all RNA isoforms, circular RNA (circRNA) is one of the most abundant forms present in cells, yet they remain the least studied. Functionally, circRNA can act as miRNA sponges, protein sponges/decoys, and regulators of transcription and translation. In the context of ovarian follicular development, the identity and roles of circRNA are relatively unknown. In the present study, high throughput RNA sequencing of granulosa cells immediately prior to and 4-h after the LH/hCG surge identified 42,381 circRNA originating from 7712 genes. A total of 54 circRNA were identified as differentially expressed between 0-h and 4-h time points (Fold Change ± 1.5, FDR ≤ 0.1), among them 42 circRNA were upregulated and 12 circRNA were downregulated. All differentially expressed circRNA between the 0-h and 4-h groups were subjected to circinteractome analysis and identified networks of circRNA-protein and circRNA-miRNA were further subjected to "micro-RNA target filter analysis" in Ingenuity Pathway Analyses, which resulted in the identification of miRNA targeted mRNAs. A comparison of these circRNA target mRNAs with LH-induced mRNAs identified Runx2, Egfr, Areg, Sult1el, Cyp19a1, Cyp11a1, and Hsd17b1 as targets of circKif2, circVcan, circMast4, and circMIIt10. These newly identified LH/hCG-induced circRNA, their target miRNA and protein networks provide new insights into the complex interactions associated with periovulatory follicular development.