Identification and Functional Analysis of circRNAs during Goat Follicular Development

Abstract

Litter size is a crucial quantitative trait in animals, closely linked to follicular development. Circular RNA (circRNA), a type of single-stranded closed-loop endogenous RNA with stable expression, plays pivotal roles in various biological processes, yet its function in goat follicular development remains unclear. In this study, we collected large (follicle diameter > 3 mm) and small (1 mm < follicle diameter < 3 mm) follicles from black goats in the Chuanzhong region for circRNA sequencing, with the aim of elucidating the functional circRNAs that influence follicle development in goats. Differential analysis revealed that 17 circRNAs were upregulated in large follicles, and 28 circRNAs were upregulated in small follicles. Functional enrichment analysis revealed significant enrichment of pathways related to reproduction, including cellular response to follicle-stimulating hormone stimulus, the PI3K-Akt signaling pathway, the MAPK signaling pathway, and the Notch signaling pathway. Based on the ceRNA mechanism, 45 differentially expressed circRNAs were found to target and bind a total of 418 miRNAs, and an intercalation network including miR-324-3p (circRNA2497, circRNA5650), miR-202-5p (circRNA3333, circRNA5501), and miR-493-3p (circRNA4995, circRNA5508) was constructed. In addition, conservation analysis revealed that 2,239 circRNAs were conserved between goats and humans. Prediction of translation potential revealed that 154 circRNAs may potentially utilize both N6-methyladenosine (m6A) and internal ribosome entry site (IRES) translation mechanisms. Furthermore, the differential expression and circularization cleavage sites of five circRNAs were validated through RT-qPCR and DNA sequencing. Our study constructed a circRNA map in goat follicle development, offering a theoretical foundation for enhancing goat reproductive performance.