miRNA-mRNA network regulation in the skeletal muscle fiber phenotype of chickens revealed by integrated analysis of miRNAome and transcriptome

Skeletal muscle fibers are primarily categorized into oxidative and glycolytic fibers, and the ratios of different myofiber types are important factors in determining livestock meat quality. However, the molecular mechanism for determining muscle fiber types in chickens was hardly understood. In this study, we used RNA sequencing to systematically compare mRNA and microRNA transcriptomes of the oxidative muscle sartorius (SART) and glycolytic muscle pectoralis major (PMM) of Chinese Qingyuan partridge chickens. Among the 44,705 identified mRNAs in the two types of muscles, 3,457 exhibited significantly different expression patterns, including 2,364 up-regulated and 1,093 down-regulated mRNAs in the SART. A total of 698 chicken miRNAs were identified, including 189 novel miRNAs, among which 67 differentially expressed miRNAs containing 42 up-regulated and 25 down-regulated miRNAs in the SART were identified. Furthermore, function enrichment showed that the differentially expressed mRNAs and miRNAs were involved in energy metabolism, muscle contraction, and calcium, peroxisome proliferator-activated receptor (PPAR), insulin and adipocytokine signaling. Using miRNA-mRNA integrated analysis, we identified several candidate miRNA-gene pairs that might affect muscle fiber performance, viz, gga-miR-499-5p/SOX6 and gga-miR-196-5p/CALM1, which were supported by target validation using the dual-luciferase reporter system. This study revealed a mass of candidate genes and miRNAs involved in muscle fiber type determination, which might help understand the molecular mechanism underlying meat quality traits in chickens.