Analysis of potential regulatory LncRNAs and CircRNAs in the oxidative myofiber and glycolytic myofiber of chickens

Abstract

Abstract SART and PMM are mainly composed of oxidative myofibers and glycolytic myofibers, respectively, and myofiber types profoundly influence postnatal muscle growth and meat quality. SART and PMM are composed of lncRNAs and circRNAs that participate in myofiber type regulation. To elucidate the regulatory mechanism of myofiber type, lncRNA and circRNA sequencing was used to systematically compare the transcriptomes of the SART and PMM of Chinese female Qingyuan partridge chickens at their marketing age. The luminance value (L*), redness value (a*), average diameter, cross-sectional area, and density difference between the PMM and SART were significant (p < 0.05). ATPase staining results showed that PMMs were all darkly stained and belonged to the glycolytic type, and the proportion of oxidative myofibers in SART was 81.7%. A total of 5 420 lncRNAs were identified, of which 365 were differentially expressed in the SART compared with the PMM (p < 0.05). The cis-regulatory analysis identified target genes that were enriched for specific GO terms and KEGG pathways (p < 0.05), including striated muscle cell differentiation, regulation of cell proliferation, regulation of muscle cell differentiation, myoblast differentiation, regulation of myoblast differentiation, and MAPK signaling pathway. Pathways and coexpression network analyses suggested that XR_003077811.1, XR_003072304.1, XR_001465942.2, XR_001465741.2, XR_001470487.1, XR_003077673.1 and XR_003074785.1 played important roles in regulating oxidative myofibers by TBX3, QKI, MYBPC1, CALM2, and PPARGC1A expression. A total of 10 487 circRNAs were identified, of which 305 circRNAs were differentially expressed in the SART compared with the PMM (p < 0.05). Functional enrichment analysis showed that differentially expressed circRNAs were involved in host gene expression and were enriched in the AMPK, calcium signaling pathway, FoxO signaling pathway, p53 signaling pathway, and cellular senescence. Novel_circ_004282 and novel_circ_002121 played important roles in regulating oxidative myofibers by PPP3CA and NFATC1 expression. Using lncRNA-miRNA/circRNA-miRNA integrated analysis, we identified many candidate interaction networks that might affect muscle fiber performance. Important lncRNA-miRNA-mRNA networks, such as lncRNA-XR_003074785.1/miR-193-3p/PPARGC1A, regulate oxidative myofibers. This study reveals that lncXR_003077811.1, lncXR_003072304.1, lncXR_001465942.2, lncXR_001465741.2, lncXR_001470487.1, lncXR_003077673.1, XR_003074785.1, novel_circ_004282 and novel_circ_002121 might regulate oxidative myofibers. The lncRNA-XR_003074785.1/miR-193-3p/PPARGC1A pathway might regulate oxidative myofibers. All these findings provide rich resources for further in-depth research on the regulatory mechanism of lncRNAs and circRNAs in myofibers.