Comprehensive Analysis of lncRNAs and circRNAs Reveals the Metabolic Specialization in Oxidative and Glycolytic Skeletal Muscles

The biochemical and functional differences between oxidative and glycolytic muscles could affect human muscle health and animal meat quality. However, present understanding of the epigenetic regulation with respect to lncRNAs and circRNAs is rudimentary. Here, porcine oxidative and glycolytic skeletal muscles, which were at the growth curve inflection point, were sampled to survey variant global expression of lncRNAs and circRNAs using RNA-seq. A total of 4046 lncRNAs were identified, including 911 differentially expressed lncRNAs (p < 0.05). The cis-regulatory analysis identified target genes that were enriched for specific GO terms and pathways (p < 0.05), including the oxidation-reduction process, glycolytic process, and fatty acid metabolic. All these were closely related to different phenotypes between oxidative and glycolytic muscles. Additionally, 810 circRNAs were identified, of which 137 were differentially expressed (p < 0.05). Interestingly, some circRNA-miRNA-mRNA networks were found, which were closely linked to muscle fiber-type switching and mitochondria biogenesis in muscles. Furthermore, 44.69%, 39.19%, and 54.01% of differentially expressed mRNAs, lncRNAs, and circRNAs respectively were significantly enriched in pig quantitative trait loci (QTL) regions for growth and meat quality traits. This study reveals a mass of candidate lncRNAs and circRNAs involved in muscle physiological functions, which may improve understanding of muscle metabolism and development from an epigenetic perspective.