Integrated Transcriptome Analysis Reveals the Crucial mRNAs and miRNAs Related to Fecundity in the Hypothalamus of Yunshang Black Goats during the Luteal Phase

SIMPLE SUMMARY: The hypothalamus plays an important role in goat reproduction. In this research, RNA sequencing technology was used to detect and identify microRNAs (miRNAs) related to goat fertility in the goat hypothalamus. The results showed that several key miRNAs may affect the expression of their target genes in different ways with the high-fecundity goats during the luteal phase (LP-HF) vs. low-fecundity goats during the luteal phase (LP-LF), thereby affecting the key reproductive processes. These results reveal the mechanism on the reproductive performance of different litter sizes during the luteal phase of goats, from the perspective of the hypothalamus. ABSTRACT: A normal estrus cycle is essential for the breeding of goats, and the luteal phase accounts for most of the estrus cycle. The corpus luteum (CL) formed during the luteal phase is a transient endocrine gland that is crucial for the reproductive cycle and pregnancy maintenance, and is controlled by many regulatory factors. However, the molecular mechanism of the hypothalamus effect on the reproductive performance of different litter sizes during the luteal phase of goats has not been elucidated. In this study, RNA-sequencing was used to analyze the mRNA and miRNA expression profiles of the hypothalamic tissues with the high-fecundity goats during the luteal phase (LP-HF) and low-fecundity goats during the luteal phase (LP-LF). The RNA-seq results found that there were 1963 differentially expressed genes (DEGs) (890 up-regulated and 1073 down-regulated). The miRNA-seq identified 57 differentially expressed miRNAs (DEMs), including 11 up-regulated and 46 down-regulated, of which 199 DEGs were predicted to be potential target genes of DEMs. Meanwhile, the functional enrichment analysis identified several mRNA-miRNA pairs involved in the regulation of the hypothalamic activity, such as the common target gene MEA1 of novel-miR-972, novel-miR-125 and novel-miR-403, which can play a certain role as a related gene of the reproductive development in the hypothalamic–pituitary–gonadal (HPG) axis and its regulated network, by regulating the androgen secretion. While another target gene ADIPOR2 of the novel-miR-403, is distributed in the hypothalamus and affects the reproductive system through a central role on the HPG axis and a peripheral role in the gonadal tissue. An annotation analysis of the DE miRNA-mRNA pairs identified targets related to biological processes, such as anion binding (GO:0043168) and small molecule binding (GO: 0036094). Subsequently, the KEGG(Kyoto Encyclopedia of Genes and Genomes) pathways were performed to analyze the miRNA-mRNA pairs with negatively correlated miRNAs. We found that the GnRH signaling pathway (ko04912), the estrogen signaling pathway (ko04915), the Fc gamma R-mediated phagocytosis (ko04666), and the IL-17 signaling pathway (ko04657), etc., were directly and indirectly associated with the reproductive process. These targeting interactions may be closely related to the reproductive performance of goats. The results of this study provide a reference for further research on the molecular regulation mechanism for the high fertility in goats.