Identification of miRNA–mRNA Networks Associated with Pigeon Skeletal Muscle Development and Growth

SIMPLE SUMMARY: Elucidating the molecular mechanisms that regulate the growth and development of skeletal muscle is an essential prerequisite for using molecular breeding technology to improve the meat production of pigeons. In the present study, we characterized the expression of miRNAs and mRNAs in pigeon skeletal muscle at two embryonic stages (E8 and E13) and two growth stages (D1 and D10). A total of 839 DEmiRNAs and 11,311 DEGs were identified among the four groups. An miRNA–mRNA network was constructed by correlation analysis and target prediction between DEmiRNAs and DEGs. Based on the connectivity degree, twenty potential hub miRNAs responsible for pigeon skeletal muscle development and growth were identified, including cli-miR-20b-5p, miR-130-y, cli-miR-106-5p, cli-miR-181b-5p, miR-1-z, etc. GO and KEGG enrichment analysis identified candidate biological processes and pathways related to pigeon skeletal muscle development and growth. Our results provide a database for further investigating the miRNA–mRNA regulatory mechanism underlying pigeon skeletal muscle development and growth. ABSTRACT: The growth and development of skeletal muscle determine the productivity of pigeon meat production, and miRNA plays an important role in the growth and development of this type of muscle. However, there are few reports regarding miRNA regulating the growth and development of skeletal muscle in pigeons. To explore the function of miRNA in regulating the growth and development of pigeon skeletal muscle, we used RNA sequencing technology to study the transcriptome of pigeons at two embryonic stages (E8 and E13) and two growth stages (D1 and D10). A total of 32,527 mRNAs were identified in pigeon skeletal muscles, including 14,378 novel mRNAs and 18,149 known mRNAs. A total of 2362 miRNAs were identified, including 1758 known miRNAs and 624 novel miRNAs. In total, 839 differentially expressed miRNAs (DEmiRNAs) and 11,311 differentially expressed mRNAs (DEGs) were identified. STEM clustering analysis assigned DEmiRNAs to 20 profiles, of which 7 were significantly enriched (p-value < 0.05). These seven significantly enriched profiles can be classified into two categories. The first category represents DEmiRNAs continuously downregulated from the developmental stage to the growth stage of pigeon skeletal muscle, and the second category represents DEmiRNAs with low expression at the development and early growth stage, and significant upregulation at the high growth stage. We then constructed an miRNA–mRNA network based on target relationships between DEmiRNAs and DEGs belonging to the seven significantly enriched profiles. Based on the connectivity degree, 20 hub miRNAs responsible for pigeon skeletal muscle development and growth were identified, including cli-miR-20b-5p, miR-130-y, cli-miR-106-5p, cli-miR-181b-5p, miR-1-z, cli-miR-1a-3p, miR-23-y, cli-miR-30d-5p, miR-1-y, etc. The hub miRNAs involved in the miRNA–mRNA regulatory networks and their expression patterns during the development and growth of pigeon skeletal muscle were visualized. GO and KEGG enrichment analysis found potential biological processes and pathways related to muscle growth and development. Our findings expand the knowledge of miRNA expression in pigeons and provide a database for further investigation of the miRNA–mRNA regulatory mechanism underlying pigeon skeletal muscle development and growth.