The spatio-temporal expression patterns of Circular RNA (circRNA) across organs and developmental stages are critical for its function and evolution analysis. However, they remain largely unclear in mammals. Here, we comprehensively analysed circRNAs in nine organs and three skeletal muscles of Guizhou miniature pig (S. scrofa), a widely used biomedical model animal. We identified 5,934 circRNAs and analysed their molecular properties, sequence conservation, spatio-temporal expression pattern, potential function, and interaction with miRNAs. S. scrofa circRNAs show modest sequence conservation with human and mouse circRNAs, are flanked by long introns, exhibit low abundance, and are expressed dynamically in a spatio-temporally specific manner. S. scrofa circRNAs show the greatest abundance and complexity in the testis. Notably, 31% of circRNAs harbour well-conserved canonical miRNA seed matches, suggesting that some circRNAs act as miRNAs sponges. We identified 149 circRNAs potentially associated with muscle growth and found that their host genes were significantly involved in muscle development, contraction, chromatin modification, cation homeostasis, and ATP hydrolysis-coupled proton transport; moreover, this set of genes was markedly enriched in genes involved in tight junctions and the calcium signalling pathway. Finally, we constructed the first public S. scrofa circRNA database, allowing researchers to query comprehensive annotation, expression, and regulatory networks of circRNAs.

DNA methylation plays a pivotal role in biological processes by affecting gene expression. However, how DNA methylation mediates phenotype difference of skeletal muscle between lean-, obese-, and mini-type pigs remains unclear. We systematically carried out comparative analysis of skeletal muscle by integrating analysis of genome-wide DNA methylation, mRNA, lncRNA and miRNA profiles in three different pig breeds (obese-type Tongcheng, lean-type Landrace, and mini-type Wuzhishan pigs). We found that the differentially methylated genes (DMGs) were significantly associated with lipid metabolism, oxidative stress and muscle development. Among the identified DMGs, 253 genes were related to body-size and obesity. A set of lncRNAs and mRNAs including UCP3, FHL1, ANK1, HDAC4, and HDAC5 exhibited inversely changed DNA methylation and expression level; these genes were associated with oxidation reduction, fatty acid metabolism and cell proliferation. Gene regulatory networks involved in phenotypic variation of skeletal muscle were related to lipid metabolism, cellular movement, skeletal muscle development, and the p38 MAPK signaling pathway. DNA methylation potentially influences the propensity for obesity and body size by affecting gene expression in skeletal muscle. Our findings provide an abundant information of epigenome and transcriptome that will be useful for animal breeding and biomedical research.

Background Circular RNAs (circRNAs) represent a novel class of non‐coding RNAs formed by a covalently closed loop and play crucial roles in many biological processes. Several circRNAs associated with myogenesis have been reported. However, the dynamic expression, function, and mechanism of circRNAs during myogenesis and skeletal muscle development are largely unknown. Methods Strand‐specific RNA‐sequencing (RNA‐seq) and microarray datasets were used to profile the dynamic circRNAome landscape during skeletal muscle development and myogenic differentiation. Bioinformatics analyses were used to characterize the circRNAome and identify candidate circRNAs associated with myogenesis. Bulk and single‐cell RNA‐seq were performed to identify the downstream genes and pathways of circFgfr2. The primary myoblast cells, C2C12 cells, and animal model were used to assess the function and mechanism of circFgfr2 in myogenesis and muscle regeneration in vitro or in vivo by RT‐qPCR, western blotting, dual‐luciferase activity assay, RNA immunoprecipitation, RNA fluorescence in situ hybridization, and chromatin immunoprecipitation. Results We profiled the dynamic circRNAome in pig skeletal muscle across 27 developmental stages and detected 52 918 high‐confidence circRNAs. A total of 2916 of these circRNAs are conserved across human, mouse, and pig, including four circRNAs (circFgfr2, circQrich1, circMettl9, and circCamta1) that were differentially expressed (|log2 fold change| > 1 and adjusted P value < 0.05) in various myogenesis systems. We further focused on a conserved circRNA produced from the fibroblast growth factor receptor 2 (Fgfr2) gene, termed circFgfr2, which was found to inhibit myoblast proliferation and promote differentiation and skeletal muscle regeneration. Mechanistically, circFgfr2 acted as a sponge for miR‐133 to regulate the mitogen‐activated protein kinase kinase kinase 20 (Map3k20) gene and JNK/MAPK pathway. Importantly, transcription factor Kruppel like factor 4 (Klf4), the downstream target of the JNK/MAPK pathway, directly bound to the promoter of circFgfr2 and affected its expression via an miR‐133/Map3k20/JNK/Klf4 auto‐regulatory feedback loop. RNA binding protein G3BP stress granule assembly factor 1 (G3bp1) inhibited the biogenesis of circFgfr2. Conclusions The present study provides a comprehensive circRNA resource for skeletal muscle study. The functional and mechanistic analysis of circFgfr2 uncovered a circRNA‐mediated auto‐regulatory feedback loop regulating myogenesis and muscle regeneration, which provides new insight to further understand the regulatory mechanism of circRNAs.