BackgroundFertility 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.ResultsThe 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.ConclusionsThis 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.