Integrated analysis of miRNA and mRNA expression profiles associated with wooden breast myopathy in broiler chickens

Wooden breast (WB) myopathy has raised a worldwide concern among broiler industry during the past decade. Despite progress in understanding its etiology from transcriptional regulation, post-transcriptional mechanisms including the regulation of microRNAs (miRNAs) remain largely unknown. In the current study, we described an integrative analysis between mRNA and miRNA expression profiles of pectoralis major muscle from normal and WB myopathic broilers. A total of 1983 differentially expressed mRNAs (DEmRNAs) and 155 DEmiRNAs were identified in WB. We screened crucial biological processes and core DEmRNAs enriched in functional pathways, and established the protein-protein interaction network. DEmiRNAs and negatively correlated DEmRNAs regulatory networks were constructed, including 44 exist DEmiRNAs and 478 DEmRNAs, forming 772 miRNA-mRNA pairs. Upregulated DEmiRNAs including gga-miR-21-3p, gga-miR-460a-5p and gga-miR-6631-5p, as well as downregulated DEmiRNAs including gga-miR-182-5p, gga-miR-183 and gga-miR-96-5p were identified as hub miRNAs. Meanwhile, functional enrichment analysis indicated that upregulated DEmRNAs in the network were enriched in biological processes of response to stimulus, inflammatory response, extracellular matrix organization, whereas downregulated DEmRNAs were enriched in carbohydrate, amino acid and nucleotide metabolic processes. Collectively, our integrative miRNA and mRNA analysis highlighted candidate miRNAs and mRNAs, as well as potential miRNA-mRNA regulatory mechanisms involved in WB myopathy in broiler chicken.