Construction and Analysis of miRNA–mRNA Interaction Network in Ovarian Tissue of Wanxi White Geese Across Different Breeding Stages

Abstract

Ovarian development significantly influences the laying performance of geese. In this study, the transcriptome analysis was conducted on the ovarian tissues of Wanxi White Geese during the pre-laying (KL), laying (CL), and ceased-laying period (XL). Short Time-series Expression Miner (STEM) analysis and miRNA–mRNA regulatory network construction were performed to identify the key genes and miRNAs regulating laying traits. Comparative analysis of KL vs. CL, CL vs. XL, and XL vs. KL groups resulted in the identification of 337, 136, and 525 differentially expressed genes (DEGs), and 258, 1131, and 909 differentially expressed miRNAs (DEMs), respectively. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis (p < 0.05) revealed that the main enrichment pathways of DEGs and DEMs at different breeding periods were Neuroactive ligand–receptor interaction, GnRH signaling pathway and Wnt signaling pathway, all associated with ovarian development. According to the three groups of common pathways, four DEGs were screened out, including INHBB, BMP5, PRL, and CGA, along with five DEMs, including let-7-x, miR-124-y, miR-1-y, and miR-10926-z, all of them may affect ovarian development. A miRNA–mRNA regulatory network was constructed through integrated analysis of DEGs and DEMs, revealing nine miRNAs highly associated with ovarian development: miR-101-y, let-7-x, miR-1-x, miR-17-y, miR-103-z, miR-204-x, miR-101-x, miR-301-y, and miR-151-x. The dual-luciferase reporter gene verified the target relationship between WIF1 and miR-204-x, suggesting that these miRNAs may influence ovarian development in Wanxi White Goose by regulating the expression levels of their target genes within ovarian tissue. This study provides a theoretical foundation for analyzing the mechanisms of ovarian development across different breeding periods and accelerating the cultivation of new breeds through post-transcriptional regulation levels.