Transcriptome Analysis of Differentially Expressed mRNA Related to Pigeon Muscle Development

Abstract

The mechanisms behind the gene expression and regulation that modulate the development and growth of pigeon skeletal muscle remain largely unknown. In this study, we performed gene expression analysis on skeletal muscle samples at different developmental and growth stages using RNA sequencing (RNA−Seq). The differentially expressed genes (DEGs) were identified using edgeR software. Weighted gene co−expression network analysis (WGCNA) was used to identify the gene modules related to the growth and development of pigeon skeletal muscle based on DEGs. A total of 11,311 DEGs were identified. WGCNA aggregated 11,311 DEGs into 12 modules. Black and brown modules were significantly correlated with the 1st and 10th day of skeletal muscle growth, while turquoise and cyan modules were significantly correlated with the 8th and 13th days of skeletal muscle embryonic development. Four mRNA−mRNA regulatory networks corresponding to the four significant modules were constructed and visualised using Cytoscape software. Twenty candidate mRNAs were identified based on their connectivity degrees in the networks, including Abca8b, TCONS−00004461, VWF, OGDH, TGIF1, DKK3, Gfpt1 and RFC5, etc. A KEGG pathway enrichment analysis showed that many pathways were related to the growth and development of pigeon skeletal muscle, including PI3K/AKT/mTOR, AMPK, FAK, and thyroid hormone pathways. Five differentially expressed genes (LAST2, MYPN, DKK3, B4GALT6 and OGDH) in the network were selected, and their expression patterns were quantified by qRT−PCR. The results were consistent with our sequencing results. These findings could enhance our understanding of the gene expression and regulation in the development and growth of pigeon muscle.