Background: Fertility is an important economic trait in the production of meat goat, and follicular development plays an important role in fertility. Although many mRNAs and microRNAs (miRNAs) have been found to play critical roles in ovarian biological processes, the interaction between mRNAs and miRNAs in follicular development is not yet completely understood. In addition, less attention has been given to the study of single follicle (dominant or atretic follicle) in goats. This study aimed to identify mRNAs, miRNAs, and signaling pathways as well as their interaction networks in the ovarian follicles (large follicles and small follicles) of uniparous and multiple Chuanzhong black goats at estrus phase using RNA-sequencing (RNA-seq) technique. Results: The results showed that there was a significant difference in the number of large follicles between uniparous and multiple goats (P < 0.05), but no difference in the number of small follicles was observed (P > 0.05). For the small follicles of uniparous and multiple goats at estrus phase, 289 differentially expressed mRNAs (DEmRNAs) and 16 DEmiRNAs were identified; and for the large follicles, 195 DEmRNAs and 7 DEmiRNAs were identified. The functional enrichment analysis showed that DE genes in small follicles were significantly enriched in ovarian steroidogenesis and steroid hormone biosynthesis, while in large follicles were significantly enriched in ABC transporters and steroid hormone biosynthesis. The results of quantitative real-time polymerase chain reaction were consistent with those of RNA-seq. Analysis of the mRNA-miRNA interaction network suggested that CD36 (miR-122, miR-200a, miR-141), TNFAIP6 (miR-141, miR-200a, miR-182), CYP11A1 (miR-122), SERPINA5 (miR-1, miR-206, miR-133a-3p, miR-133b), and PTGFR (miR-182, miR-122) might be related to fertility, but requires further research on follicular somatic cells. Conclusions: This study was used for the first time to reveal the DEmRNAs and DEmiRNAs as well as their interaction in the follicles of uniparous and multiple goats at estrus phase using RNA-seq technology. Our findings provide new clues to uncover the molecular mechanisms and signaling networks of goat reproduction that could be potentially used to increase ovulation rate and kidding rate in goat.

MicroRNAs (miRNAs) play a significant role in ovarian follicular maturity, but miRNA expression patterns in ovarian stroma (OS), large follicles (LF), and small follicles (SF) have been rarely explored. We herein aimed to identify miRNAs, their target genes and signaling pathways, as well as their interaction networks in OS, LF, and SF of Chuanzhong black goats at the estrus phase using small RNA-sequencing. We found that the miRNA expression profiles of LF and SF were more similar than those of OS—32, 16, and 29 differentially expressed miRNAs were identified in OS vs. LF, OS vs. SF, and LF vs. SF, respectively. Analyses of functional enrichment and the miRNA-targeted gene interaction network suggested that miR-182 (SMC3), miR-122 (SGO1), and miR-206 (AURKA) were involved in ovarian organogenesis and hormone secretion by oocyte meiosis. Furthermore, miR-202-5p (EREG) and miR-485-3p (FLT3) were involved in follicular maturation through the MAPK signaling pathway, and miR-2404 (BMP7 and CDKN1C) played a key role in follicular development through the TGF-β signaling pathway and cell cycle; nevertheless, further research is warranted. To our knowledge, this is the first study to investigate miRNA expression patterns in OS, LF, and SF of Chuanzhong black goats during estrus. Our findings provide a theoretical basis to elucidate the role of miRNAs in follicular maturation. These key miRNAs might provide candidate biomarkers for the diagnosis of follicular maturation and will assist in developing new therapeutic targets for female goat infertility.

Follicular development and maturation has a significant impact on goat reproductive performance, and it is therefore important to understand the molecular basis of this process. The importance of long non-coding RNAs (lncRNAs) in mammalian reproduction has been established, but little is known about the roles of lncRNAs in different follicular stages, especially in goats. In this study, RNA sequencing (RNA-seq) of large follicles (>10 mm) and small follicles (<3 mm) of Chuanzhong black goats was performed to investigate the regulatory mechanisms of lncRNAs and mRNAs in follicular development and maturation. A total of 8 differentially expressed lncRNAs (DElncRNAs) and 1,799 DEmRNAs were identified, and the majority of these were upregulated in small follicles. MRO, TC2N, CDO1, and NTRK1 were potentially associated with follicular maturation. KEGG pathway analysis showed that the DEmRNAs involved in ovarian steroidogenesis (BMP6, CYP11A1, CYP19A1, 3BHSD, STAR, LHCGR, and CYP51A1) and cAMP signaling play roles in regulating follicular maturation and developmental inhibition respectively. Five target pairs of DElncRNA-DEmRNA, namely, ENSCHIT00000001255-OTX2, ENSCHIT00000006005-PEG3, ENSCHIT00000009455-PIWIL3, ENSCHIT00000007977-POMP, and ENSCHIT00000000834-ACTR3 in co-expression analysis provide a clue in follicular development and maturation of lncRNA-mRNA interaction. Our findings provide a valuable resource for lncRNA studies, and could potentially provide a deeper understanding of the genetic basis and molecular mechanisms of goat follicular development and maturation.