Background: Litter size is an important index of mammalian prolificacy and is determined by the ovulation rate. The ovary is a crucial organ for mammalian reproduction and is associated with follicular development, maturation and ovulation. However, prolificacy is influenced by multiple factors, and its molecular regulation in the follicular phase remains unclear. Methods: Ten female goats with no significant differences in age and weight were randomly selected and divided into either the high-yielding group (n = 5, HF) or the low-yielding group (n = 5, LF). Ovarian tissues were collected from goats in the follicular phase and used to construct mRNA and miRNA sequencing libraries to analyze transcriptomic variation between high- and low-yield Yunshang black goats. Furthermore, integrated analysis of the differentially expressed (DE) miRNA-mRNA pairs was performed based on their correlation. The STRING database was used to construct a PPI network of the DEGs. RT–qPCR was used to validate the results of the predicted miRNA-mRNA pairs. Luciferase analysis and CCK-8 assay were used to detect the function of the miRNA-mRNA pairs and the proliferation of goat granulosa cells (GCs). Results: A total of 43,779 known transcripts, 23,067 novel transcripts, 424 known miRNAs and 656 novel miRNAs were identified by RNA-seq in the ovaries from both groups. Through correlation analysis of the miRNA and mRNA expression profiles, 263 negatively correlated miRNA-mRNA pairs were identified in the LF vs. HF comparison. Annotation analysis of the DE miRNA-mRNA pairs identified targets related to biological processes such as “estrogen receptor binding (GO:0030331)”, “oogenesis (GO:0048477)”, “ovulation cycle process (GO:0022602)” and “ovarian follicle development (GO:0001541)”. Subsequently, five KEGG pathways (oocyte meiosis, progesterone-mediated oocyte maturation, GnRH signaling pathway, Notch signaling pathway and TGF-β signaling pathway) were identified in the interaction network related to follicular development, and a PPI network was also constructed. In the network, we found that CDK12, FAM91A1, PGS1, SERTM1, SPAG5, SYNE1, TMEM14A, WNT4, and CAMK2G were the key nodes, all of which were targets of the DE miRNAs. The PPI analysis showed that there was a clear interaction among the CAMK2G, SERTM1, TMEM14A, CDK12, SYNE1 and WNT4 genes. In addition, dual luciferase reporter and CCK-8 assays confirmed that miR-1271-3p suppressed the proliferation of GCs by inhibiting the expression of TXLNA. Conclusions: These results increase the understanding of the molecular mechanisms underlying goat prolificacy. These results also provide a basis for studying interactions between genes and miRNAs, as well as the functions of the pathways in ovarian tissues involved in goat prolificacy in the follicular phase.

Granulosa cell (GC) proliferation provides essential conditions for ovulation in animals. A previous study showed that DENND1A plays a significant role in polycystic ovary syndrome. However, the modulation of DENND1A in GCs remains unclear. Our previous integrated analysis of miRNA–mRNA revealed that the 3'-untranslated region of DENND1A could be a target of chi-miR-324-3p. In this study, we used quantitative reverse transcription polymerase chain reaction (RT-qPCR) to investigate DENND1A expression in ovarian tissues of high- and low-yielding goats. Furthermore, dual-fluorescent reporter vector experiments, Cell Counting Kit-8 (CCK-8) assay, and RT-qPCR were used to elucidate the regulatory pathway of chi-miR-324-3p-DENND1A in GCs. The results revealed an opposite tendency between the expressions of chi-miR-324-3p and DENND1A in the ovaries of high- and low-yielding goats. The CCK-8 assay indicated that chi-miR-324-3p overexpression significantly suppressed GC proliferation, whereas chi-miR-324-3p inhibition promoted GC proliferation. In addition, the expressions of GC proliferation markers LHR, Cylin D2, and CDK4 showed the same tendency. The dual-fluorescent reporter assay revealed that chi-miR-324-3p directly targeted DENND1A, and the RT-qPCR results revealed that DENND1A expression was inhibited by chi-miR-324-3p. In summary, chi-miR-324-3p inhibited the proliferation of GCs by targeting DENND1A.

CircRNAs acting as miRNA sponges play important roles in the growth process of animal individuals. The prolificacy trait of goats is involved in many pathways, however, the variation of circRNA expression profiles in the different phases of the estrus cycle at high and low fecundity groups is still unknown. Here, we analyzed the circRNA profiles of ovarian tissues among high and low fecundity groups in the follicular phase (HF vs LF), high and low fecundity groups in the luteal phase (HL vs LL), and high and low fecundity in the whole estrus cycle (HF vs HL and LF vs LL) using RNA-seq. A total of 283 (114 upregulated and 169 downregulated), 559 (299 upregulated and 260 downregulated), 449 (254 upregulated and 195 downregulated), and 314 (210 upregulated and 104 downregulated) differentially expressed (DE) circRNAs were screened in HF vs LF, HF vs HL, HL vs LL, and LF vs LL groups, respectively. Enrichment analysis suggested that the targeting of DE circRNAs was mainly enriched in oocyte meiosis, the GnRH signaling pathway, and estrogen signaling pathway. After integrating our previous study of miRNA-seq, there were 56 miRNAs that could target to 192 DE circRNAs, including the miR-133 family (including miR-133a-3p and miR-133b), miR-129-3p, and miR-21, which also had important influence on the prolificacy trait of goats. Then, 18 circRNAs with coding potential were obtained by four software predictions, and 9 circRNAs were validated by RT-qPCR. Together, circRNAs play a key role in the prolificacy trait and the transformation of the follicular phase to the luteal phase in the estrus cycle of goats.

Retinol-binding protein 4 (RBP4), a member of the lipocalin family, is a specific carrier of retinol (vitamin A) in the blood. Numerous studies have shown that RBP4 plays an important role in mammalian embryonic development and that mutations in RBP4 can be used for the marker-assisted selection of animal reproductive traits. However, there are few studies on the regulation of reproduction and high-prolificacy traits by RBP4 in goats. In this study, the 5′ flanking sequence of RBP4 was amplified, and a G>C polymorphism in the promoter region -211 bp (g.36491960) was detected. An association analysis revealed that the respective first, second and third kidding number and mean kidding number of nanny goats with CC and GC genotypes (2.167 ± 0.085, 2.341 ± 0.104, 2.529 ± 0.107 and 2.189 ± 0.070 for CC and 2.052 ± 0.047, 2.206 ± 0.057, 2.341 ± 0.056 and 2.160 ± 0.039 for GC) were significantly higher (p < 0.05) than those with the GG genotype (1.893 ± 0.051, 2.027 ± 0.064, 2.107 ± 0.061 and 1.74 ± 0.05). The luciferase assay showed that luciferase activity was increased in C allele individuals compared with that in G allele individuals. A competitive electrophoretic mobility shift assay (EMSA) showed that individuals with the CC genotype had a stronger promoter region binding capacity than those with the GG genotype. In addition, transcription factor prediction software showed that the RBP4 g.36491960G>C mutation added a novel binding site for transcription factor DP-1 (TFDP1). RT–qPCR results showed that the expression of TFDP1 was significantly higher in the high-prolificacy group than in the low-prolificacy group, and the expression of RBP4 was higher in both the CC and GC genotypes than that in the GG genotype. TFDP1 overexpression significantly increased the expression of RBP4 mRNA (p < 0.05) and the expression of the cell proliferation factors cyclin-D1, cyclin-D2 and CDK4 (p < 0.05). The opposite trend was observed after interference with TFDP1. Both the EdU and CCK-8 results showed that TFDP1 expression could regulate the proliferation of goat ovarian granulosa cells. In summary, our results showed that RBP4 g.36491960G>C was significantly associated with fecundity traits in goats. The g.36491960G>C mutation enhanced the transcriptional activity of RBP4 and increased the expression of RBP4, thus improving the fertility of Yunshang black goats.

The oviduct is associated with embryo development and transportation and regulates the pregnancy success of mammals. Previous studies have indicated a molecular mechanism of lncRNAs in gene regulation and reproduction. However, little is known about the function of lncRNAs in the oviduct in modulating goat kidding numbers. Therefore, we combined RNA sequencing (RNA-seq) to map the expression profiles of the oviduct at the luteal phase from high- and low-fecundity goats. The results showed that 2023 differentially expressed mRNAs (DEGs) and 377 differentially expressed lncRNAs (DELs) transcripts were screened, and 2109 regulated lncRNA-mRNA pairs were identified. Subsequently, the genes related to reproduction (IGF1, FGFRL1, and CREB1) and those associated with embryonic development and maturation (DHX34, LHX6) were identified. KEGG analysis of the DEGs revealed that the GnRH- and prolactin-signaling pathways, progesterone-mediated oocyte maturation, and oocyte meiosis were related to reproduction. GSEA and KEGG analyses of the target genes of DELs demonstrated that several biological processes and pathways might interact with oviduct functions and the prolificacy of goats. Furthermore, the co-expression network analysis showed that XLOC_029185, XLOC_040647, and XLOC_090025 were the cis-regulatory elements of the DEGs MUC1, PPP1R9A, and ALDOB, respectively; these factors might be associated with the success of pregnancy and glucolipid metabolism. In addition, the GATA4, LAMA2, SLC39A5, and S100G were trans-regulated by lncRNAs, predominantly mediating oviductal transport to the embryo and energy metabolism. Our findings could pave the way for a better understanding of the roles of mRNAs and lncRNAs in fecundity-related oviduct function in goats.