Circular RNAs (circRNAs) play a significant regulatory role during skeletal muscle development. To identify circRNAs during postnatal skeletal muscle development in chickens, we constructed 12 cDNA libraries from breast muscle tissues of Chinese Gushi chickens at 6, 14, 22, and 30 weeks and performed RNA sequencing. In total, 2112 circRNAs were identified, and among them 79.92% were derived from exons. CircRNAs are distributed on all chromosomes of chickens, especially chromosomes 1-9 and Z. Bioinformatics analysis showed that each circRNA had an average of 38 miRNA binding sites, 61.32% of which have internal ribosomal entry site (IRES) elements. Furthermore, in total 543 differentially expressed circRNAs (DE-circRNAs) were identified. Functional enrichment analysis revealed that DE-circRNAs source genes are engaged in biological processes and muscle development-related pathways; for example, cell differentiation, sarcomere, and myofibril formation, mTOR signaling pathway, and TGF-β signaling pathway, etc. We also established a competitive endogenous RNA (ceRNA) regulatory network associated with skeletal muscle development. The results in this report indicate that circRNAs can mediate the development of chicken skeletal muscle by means of a complex ceRNA network among circRNAs, miRNAs, genes, and pathways. The findings of this study might help increase the number of known circRNAs in skeletal muscle tissue and offer a worthwhile resource to further investigate the function of circRNAs in chicken skeletal muscle development.
Keywords: Gushi chickens; ceRNA; circular RNAs; regulatory network; skeletal muscle.