Genome-wide analysis of circular RNA-mediated ceRNA regulation in porcine skeletal muscle development
BMC Genomics volume 24, Article number: 196 (2023)
As a diverse and abundant class of endogenous RNAs, circular RNAs (circRNAs) participate in various biological processes including cell proliferation and apoptosis. Nevertheless, few researchers have investigated the role of circRNAs in muscle development in cultivated pigs.
In this study, we used RNA-seq to construct circRNA expression profiles in skeletal muscle of Jinfen White pigs at the age of 1, 90, and 180 days. Among the 16,990 identified circRNAs, 584 circRNAs were differentially expressed. Moreover, the enrichment analysis of DE circRNA host genes showed that they were mainly involved in muscle contraction, muscle organ development and muscle system processes, as well as AMPK and cAMP-related signal pathways. We also constructed a circRNA–miRNA–mRNA co-expression network to find key circRNAs which many involved in the regulation of porcine skeletal muscle development through the competitive endogenous RNA (ceRNA) mechanism. It is noteworthy that circ_0018595/miR-1343/PGM1 axis may play a regulatory role in the development of porcine skeletal muscle.
This study identified the circRNAs and present the circRNA expression profile in the development of pigs, revealed that DE circRNA host genes participate in different cell fates and enriched the porcine ceRNA network. Thus, this work will become a valuable resource for further in-depth study of the regulatory mechanism of circRNA in the development of porcine skeletal muscle.
Skeletal muscle is the main tissue of the animal body, and it has important functions in exercise, heat production, maintenance of body shape, storage of protein, and protection of organs [1, 2]. Muscle fiber is the basic functional unit of skeletal muscle, and its quantity remains unchanged after the birth of mammals. The growth and development of skeletal muscle is achieved through the hypertrophy of existing muscle fiber . In the process of myogenesis, skeletal muscle satellite cells gradually form mature muscle tissue after undergoing myoblasts, multinucleated myotubes, muscle fibers and other stages, and play an indispensable role in the growth of skeletal muscle [4, 5]. Regulation of the activity of myogenic stem cells involves many transcription factors, such as myoblast determination protein 1 (MyoD1), myogenin (MyoG), myogenic factor 5 (Myf5), and myogenic factor 6 (Myf6), as well as other myogenesis-related regulatory transcription factors, including the paired box 3 (PAX3), paired box 7 (PAX7), and myocyte enhancer factor 2 (MEF2) families [6,7,8,9]. Among these, MyoD1 and Myf5 are involved in the first stage of skeletal muscle development, which promoting the proliferation and differentiation of myogenic progenitor cells into myoblasts, while MyoG determines the generation of myotubes, and Myf6 is involved in cell differentiation and cell fate [10,11,12,13]. Furthermore, the expression of MyoG is regulated by MRF4 and MyoD .
In addition to the above roles of transcription factors, a growing number of studies have established that non-coding RNAs (ncRNAs) are also involved in the regulation of myogenesis. These include microRNAs (miRNAs), long non-coding RNAs (lncRNAs), and circular RNAs (circRNAs), whose regulatory roles in the growth of skeletal muscle are gradually being revealed [15,16,17]. With the development of high-throughput sequencing technology, circRNA has received increasing attention in recent years. As the most recently emerging class of ncRNA, circRNA is essentially a closed RNA transcript that is produced by reverse splicing of precursor mRNA (pre-mRNA). The 3′ and 5′ ends are covalently linked, and they differ from their linear counterparts in the absence of a 5′ cap and polyadenylated [poly(A)] tail . Moreover, circRNAs have been identified in most species, and are generally expressed at low levels, but partly more abundantly than their host genes (also referred to as source genes) [19,20,21]. Compared with lncRNAs, circRNAs display highly conserved expression patterns and transcript sequences among species, and they are widely expressed in various tissues and organs [22, 23]. Much current research on circRNAs focuses on aspects of cancer, including cell proliferation, apoptosis, invasion, and chemoresistance [24,25,26,27].
The research on circRNA in skeletal muscle mainly focuses on the proliferation and differentiation of myoblasts. CircRBFOX2 and circFGFR2 can regulate the development of chicken myoblasts, which can sponge and negatively regulate miR-206 and miR-133a-5p respectively [28, 29]. CircFAM188B encoded a novel protein named circFAM188B-103aa that promoted the proliferation of chicken skeletal muscle satellite cells but suppressed differentiation . CircFUT10 sponge miR-133a leaded to inhibition of cell proliferation and ultimately enhanced myoblastic differentiation . At present, the research on circRNA regulation of porcine skeletal muscle growth mainly depends on the in-depth mining of sequencing data. Sun constructed a pig skeletal muscle injury model and screened out the differentially expressed circRNA circSCDE1. Because it is highly conservative in pigs and mice, C2C12 was used to verify the specific mechanism of circCSDE1 regulating the proliferation and differentiation of myoblasts . Ma identified and predicted circRNA in Duroc pigs with different daily weight gain . Yan screened the highly conservative circFGFR2, and explained that circFGFR2 regulates myogenesis and skeletal muscle regeneration through feedback loop on C2C12 . Jinfen White pigs are cultivated from four parental breeds, which have the advantages of good meat quality, strong resistance to stress, and rapid growth [35,36,37]. However, in Jinfen White pigs, the characteristic and mechanism of regulation by circRNAs at different growth stages are still unclear. Thus, further exploration is needed.
Therefore, in this study, RNA-seq data in the longissimus dorsi muscle of Jinfen White pigs were sampled at three growth stages (9 samples: 1 d, 90 d and 180 d). We investigated the feature, expression and potential function of identified circRNAs in Jinfen White pigs and constructed a ceRNA network to find the key circRNAs in the development of pigs. This will advance our knowledge of the feature and predicted function of the circRNAs during the developmental stages of pigs.
Overview of sequencing data
To understand the expression characteristics of circRNA in porcine muscle at different developmental stages, we performed ribosomal RNA-depleted RNA-seq on muscle sampled at three time points. Firstly, we constructed nine non-coding RNA libraries named JFW_1d_1, JFW_1d_2, JFW_1d_3, JFW_90d_1, JFW_90d_2, JFW_90d_3, JFW_180d_1, JFW_180d_2, and JFW_180d_3. Secondly, we purified and sequenced the RNA using the Illumina paired-end RNA-seq method. For these libraries, an average of 100,961,956 raw reads were obtained from the platform, and after filtering, a final dataset of 100,294,750 clean reads was obtained (Table 1). The JFW_180d_3 group contained the greatest number of reads (106,567,824), and the JFW_1d_1 group had the least (91,942,766). The Q30 quality threshold (Phred quality score > 30) was met for 95.2% of the data. Taken together, these results confirmed that the sample and sequencing data generated in this study were of good quality and were reliable for subsequent bioinformatic analysis.
Characterization of circRNAs in porcine skeletal muscle
After mapping the clean data to the pig genome (Sus scrofa11.1) using the find_circ and CIRI software, the intersection selected by the two software is used as the circRNAs for final identification. We identified a total of 16,990 highly credible circRNAs (expressed in at least three or more samples) from the nine libraries, which were transcribed from 4775 source genes (summarized in more detail in Supplementary Table S1). In addition, we mapped the chromosomal distribution of these circRNAs, finding them to be transcribed from all chromosomes but with an uneven distribution (Fig. 1A). In agreement with other earlier work , we observed that the highest number of circRNAs originated from chromosomes 1, 6, and 13 in this study (Fig. 1A). Furthermore, since circRNAs are primarily derived from annotated exons , we also found that > 90% identified circRNAs in this study contained exon sequences (Fig. 1B). The length distribution of the predicted circRNAs showed that most of them were less than 1,000 nt in length (Fig. 1C), which was also consistent with previous reports .
Differential expression of circRNAs
To better understand the differences in circRNA expression patterns in pig muscle at 1, 90, and 180 days, we used DESeq software to identify DE circRNAs between sample pairs (i.e., JFW_90d vs. JFW_1d, JFW_180d vs. JFW_1d, and JFW_180d vs. JFW_90d). We also performed hierarchical clustering analysis of the DE circRNAs. This revealed distinct differential expression patterns for the identified circRNAs between the three time points in muscle, with JFW_1d showing significant differences compared to the other two groups (Fig. 2A). In total, 584 DE circRNAs were detected using the above pairwise comparisons. The number of differentially expressed genes was 477, 63, and 255 for the JFW_180d vs. JFW_1d, JFW_180d vs. JFW_90d, and JFW_90d vs. JFW_1d comparisons, respectively (Fig. 2B). We further constructed a Venn diagram of the identified DE circRNAs in the three pairwise comparisons to find the shared circRNAs (Fig. 2C). The expression level of two kinds of circRNA (circ_0009188 and circ_0018581) increased significantly with age in the three growth stages of Jinfen White pigs. One of these, circ_0009188, is generated by circularization of two exons of tubulin tyrosine ligase-like 7 (TTLL7) and is 311 nt in length. The other, circ_0018581, arises through circularization of two exons of glutamate decarboxylase-like protein 1 (GADL1), and has a length of 608 nt.
Validation of DE circRNAs at different developmental stages of muscle
To validate the reliability of temporal circRNA expression profiles derived from RNA-seq data, divergent primers were designed for eight randomly selected DE circRNAs, and the expression levels of which were quantified by qRT-PCR (Fig. 3A). The results were consistent with the expression patterns of these circRNAs obtained through RNA-seq data. In addition, the expected size PCR products were obtained by amplification using cDNA as the PCR template. The post-splice sites were verified by Sanger sequencing (Fig. 3B). These results suggested that the identified circRNAs in this study are credible.
Functional enrichment analysis of DE circRNA host genes
Previous studies have shown that circRNAs could regulate host gene transcription by competing with linear pre-mRNA splicing [41, 42]. To explore the potential functions of the identified 584 DE circRNAs in the development of skeletal muscle, we performed GO and KEGG pathway enrichment analyses of their host genes . Notably, we found that these host genes mainly have known roles in muscle biology, including muscle contraction, muscle organ development, muscle system processes, and muscle structure development (Fig. 4A). The host genes CHD2, SMAD3, BMPR1A, HOMER1, DMD, and MYL2 are involved in muscle organ development; HOMER1, DMD, and ATP11A are involved in myotube differentiation; and NFATC1, BMPR1A, and PDCD4 are involved in the negative regulation of muscle cell differentiation. The KEGG pathway enrichment analysis showed that the DE circRNA host genes were enriched in top 25 pathways, including skeletal muscle fiber-related signaling pathways, and the cAMP and AMPK signaling pathways and so on (Fig. 4B). PDE4B, NFATC1, PLCE1, ATP2B4, PPP1R12A, and CAMK2B were enriched in the cAMP signaling pathway, and EEF2K, FBP2, PPP2CB, PFKFB1, and PRKAA2 were enriched in AMPK signal pathway.
Construction of a potential circRNA-miRNA-mRNA regulatory network
CircRNAs often act as miRNA molecular sponges to indirectly regulate the expression of miRNA target genes, thereby meaning that they participate in the related biological processes . Therefore, we integrated miRNA and mRNA library data (The sequencing results obtained at the same time, the accession number is consistent with the circRNA) with the results from the present study to constructed a ceRNA network and identified key circRNAs associated with pork quality and regulation of cell proliferation and differentiation. This ceRNA network contained 229 circRNAs, 55 miRNAs and 143 mRNAs (Supplementary Table S2). Because the ceRNA network is too large, we notice the top 20 DE circRNAs with the most obvious up- or down-regulation profiles to present two ceRNA network maps (Fig. 5). In ceRNA network, circ_0018595 is predicted to function as a sponge for ssc-miR-1343, which binds to the mRNAs transcribed from the PGM1, TNNT3, and IFNAR2 genes (Fig. 5A).
To explore the functional significance of mRNAs in the circRNA–miRNA–mRNA regulatory network, we performed functional enrichment analysis on 143 targeted mRNAs. As a result the GO enrichment analysis showed that these target mRNAs were significantly enriched in pyruvate metabolism, hexose metabolism, striated muscle contraction, glycolysis process, muscle system process, and muscle contraction. More specifically, PGAM2, ALDOA, ACTN3, PKM, and GPI are involved in the process of glycolysis, and ACTN3, PI16, WDR1 are involved in striated muscle development. It is worth noting that ACTN3 is also involved in the positive regulation of fast-twitch skeletal muscle fiber contraction, skeletal muscle fiber development, skeletal muscle tissue growth, and other important processes (Supplementary Figure S1A). Furthermore, KEGG enrichment analysis showed that the same target mRNAs were mainly enriched in carbon metabolism, glycolysis/gluconeogenesis, and regulation of the actin cytoskeleton. Specifically, PGAM2, ALDOA, GPI, and EFCAB7 are involved in glycolysis and gluconeogenesis, HOMER3 and PDPK1 are involved in the FOXO signaling pathway, and PDPK1, ITGA5, and ITGB4 are associated with the PI3K-AKT signaling pathway (Supplementary Figure S1B). These results suggested that these circRNAs in the ceRNA network may play important role in the muscle development.
Experimental validation of circ_0018595
From the above network, we selected circ_0018595 for follow-up study, speculating that it might be involved in regulating growth and development of porcine skeletal muscle via the predicted circ_0018595/miR-1343/PGM1 axis (Fig. 5A). Circ_0018595 is an exonic circRNA generated from exons 3 and 4 of STT3B. It was amplified from longissimus dorsi cDNA using specific divergent primers, and the product was subjected to Sanger sequencing to validate the ligation site (Fig. 6A). We also processed the total RNA extract with RNase R and performed qRT-PCR, revealing that circ_0018595 was more resistant to RNase R than STT3B mRNA and 18S rRNA as expected (Fig. 6B). The circ_0018595 was expressed in various tissues, especially in lung, liver, skeletal muscle, and kidney (Fig. 6C). We also examined the expression of circ_0018595 in muscle of Jinfen White pigs at 1, 90, and 180 days, and found it was significantly associated with developmental stage, with strong upregulation between 1 and 90 days, and maintained at 180 days (Fig. 6D). Furthermore, the results of a nucleocytoplasmic separation assay showed that circ_0018595 is primarily localized in the cytoplasm, which was consistent with its proposed role as a sponge for miRNAs (Fig. 6E). Finally, the RNAhybrid website found that circ_0018595 and miR-1343 had potential binding sites, which was also found between miR-1343 and PGM1 mRNA transcript (Fig. 6F). The function and mechanism experiment of circ_0018595 in the skeletal muscle development will be carried out in the future work.
CircRNAs are abundant in transcriptomes and participate in various biological activities by regulating transcription, splicing, mRNA stability, and translation. At present, many studies have investigated circRNAs as potential therapeutic targets for various diseases [45, 46]. With the ongoing development of related research strategies and methods, many biological functions of circRNAs in animals have also been emerging, and circRNAs have been found to regulate the proliferation and differentiation of animal skeletal muscle cells by acting as ceRNAs and encoding small peptides . For example, Shen et al. , found that circTMTC1 can bind miR-128-3p and regulate expression of myostatin (MSTN), thereby inhibiting the proliferation and differentiation of chicken skeletal muscle satellite cells. Pigs are important farm animals that provide meat for humans and are often used for disease models in medical research. In 2017, Liang and coworkers constructed the first publicly available Sus scrofa circRNA database through analysis of nine organs of Guizhou pigs at three developmental stages . Current research on circRNAs in porcine skeletal muscle mainly relies on the mining and discussion of sequencing data, and the specific mechanisms of circRNA involvement in skeletal muscle growth and fat deposition are rarely studied [49,50,51]. Here, we used RNA-seq to construct circRNA expression profiles for Jinfen White pig skeletal muscle at three developmental stages (1, 90, and 180 d). Compared with previous studies, our sequencing data are robust, with > 97% of bases meeting the Q20 threshold.
RNA-seq is widely used in the study of circRNA in model animals and non-model animals. We use the library construction method of linear RNA for RNA-seq identification. The circRNA detected by different software and algorithms is different. We use the intersection of circRNA predicted by find_circ and CIRC software can improve the accuracy of circRNA. This method has been used in mice , cattle , sheep , chicken  and other species . We identified a total of 16,990 expressed circRNAs, which is much higher than the number reported in one recent study. We also analyzed the chromosomal distribution, composition, and length of these circRNAs. qRT-PCR confirmed the expression patterns derived from RNA-seq data for eight examples, further establishing the accuracy of our sequencing results, along with the developmental stage specificity of circRNA expression profiles. These findings are also consistent with previous studies [49,50,51]. Our results therefore enrich the known porcine circRNA library and confirm that circRNAs are consistent across pig species. In addition, we identified a total of 584 circRNAs that were differentially expressed between the three developmental stages. This may be helpful in identifying circRNAs that regulate skeletal muscle growth and development, and opening up new avenues of future exploration to elucidate the specific mechanisms involved.
Usually, circRNAs are generated by back-splicing of pre-mRNAs. They have notable structural stability and spatiotemporal specificity, and can positively or negatively regulate transcription of their parent genes . We performed enrichment analysis on the host genes of all DE circRNAs identified in this study, and found that these host genes were mainly involved in the AMPK and cAMP signaling pathways. The AMPK signaling pathway is known as a regulator of oxidative metabolism in skeletal muscle, and it promotes the conversion of glycolytic muscle fibers into oxidative fibers [58, 59]. For example, in C2C12 myotubes, AMPK signaling drives expression of skeletal slow-twitch muscle fiber genes . On the other hand, the cAMP signaling pathway affects protein synthesis rate, promotes skeletal muscle hypertrophy, and is expected to become a new target for the treatment of muscle mass loss during atrophic diseases [61,62,63]. Interestingly, more than one circRNA can be transcribed from a single gene due to alternative splicing mechanisms. In this study, we found that Lim-domain only protein 7 (Lmo7) produced 15 circRNAs, but among these, only circ_0016243 and circ_0016253 had high expression levels, suggesting that these two novel circRNAs may play important roles in cell differentiation in Jinfen White pigs. LMO7 is essential for skeletal muscle development and encodes a transcription factor that binds to the PAX3, MyoD, and Myf5 promoters to activate expression of key myogenic differentiation genes in C2C12 cells . Thus, downregulation of LMO7 was shown to inhibit myogenesis in chicken myoblasts . In addition, circLMO7 (formed by circularization of exons 3–5 of bovine LMO7), can target miR-378a-3p to inhibit the differentiation and apoptosis of myoblasts, and the sequence of bovine circLMO7 is completely different from that of porcine circ_0016243 and circ_0016253 in our study .
As mentioned above, circRNAs can act as miRNA sponges to regulate the expression of target genes, through the ceRNA mechanism. We constructed a circRNA-miRNA-mRNA co-expression network using high-throughput sequencing of mRNA and miRNA expression network data. For some of the miRNAs in this ceRNA network, a role in regulating proliferation and differentiation of skeletal muscle satellite cells has already been reported in other animals. These include miR-370 , miR-214 , miR-133 , and miR-885 . We predicted that circ_0000576 could act as a sponge for ssc-miR-15a, ssc-miR-98, and ssc-miR-885-3p, modulating the target genes of these miRNAs to regulate muscle growth and development. In addition, circ_0004213 is an exon-derived circRNA from the polypyrimidine tract binding protein 1 (PTBP1) gene, which itself encodes a splicing factor involved in regulatory processes related to cell differentiation, cell cycle and apoptosis, cell motility, cell metabolism, and the immune response . We also predicted that circ_0004213 could target miR-370, which in turn targets α-actinin-3 (ACTN3). Involvement of miR-370 has been reported in many cellular functions, including cell proliferation, migration, and differentiation , while ACTN3 encodes the sarcomeric α-actinin-3 protein that is a component of the Z line in mammalian skeletal muscle fibers . We also identified circ_0018595, which is generated by cyclization of exons 3–4 of STT3B and is predicted to bind miR-1343. This miRNA is involved in cell proliferation and apoptotic processes, as well as affecting the pluripotency of stem cells. In the ovaries of highly fertile sows, miR-1343 suppresses expression of transforming growth factor-β receptor type 1 (TGFBR1) to inhibit pig granule cell proliferation and promote apoptosis, thus affecting sow fecundity . Furthermore, miR-1343 specifically binds to the 3´UTR of Orthodenticle homeobox 2 (OTX2) and inhibits endogenous OTX2 expression in pig pluripotent stem cells (piPSCs), enhancing the expression of pluripotency genes and thus maintaining piPSCs pluripotency in pigs . Finally, we chose to focus specifically on PGM1, which has the highest expression level in porcine longissimus dorsi and biceps femoris. It is involved in the glycolytic pathway, which may be related to pork quality and productivity [74, 75], since reducing the proportion of glycolytic fibers in porcine longissimus dorsi can slow down the rate and extent of muscle pH decline after slaughter . Thus, PGM1 may be involved in regulating skeletal muscle growth and lipid deposition. Given these considerations, the circ_0018595/miR-1343/PGM1 axis may play a key role in the growth and development of skeletal muscle in Jinfen White pigs.
We also have some limitations of this study. All circRNAs and their target miRNAs or mRNAs were predicted computationally, and the specific mechanisms of these circRNAs in regulating skeletal muscle still require experimental verification. A series of complex potential mechanisms, such as circRNAs interacting with proteins, regulating parental gene transcription, and participating in protein coding, which remain to be explored.
This study identified and revealed the expression profile and potential role of circRNA in the longissimus dorsi muscle of Jinfen White pigs by using RNA-seq data. Based on enrichment analysis, the host genes of DE circRNA are mainly AMPK and cAMP pathway. In the constructed ceRNA network, a new axis in the network circ_0018595/miR-1343/PGM1 may participate in the regulation of skeletal muscle growth. In conclusion, we have identified a series of candidate circRNAs, which have potential regulatory effects on the growth and development of porcine skeletal muscle and provide important references for future research.
Experimental animals and samples
All experimental animals and procedures in this experiment were approved by the Institutional Animal Care and Use Committee of Shanxi Agricultural University (Shanxi, China). The 9 Jinfen White pigs used in the experiment came from Datong Pig Breeding Farm (Shanxi, China) and were raised under standard conditions without restriction on feeding and drinking. Three heads per group (denoted as samples JFW_1d_1, JFW_1d_2, JFW_1d_3; JFW_90d_1, JFW_90d_2, JFW_90d_3; JFW_180d_1, JFW_180d_2, JFW_180d_3) were sacrificed by corona and neck exsanguination on days 1, 90, and 180, respectively. Tissues such as heart, liver, spleen, lung, kidney, pancreas, longissimus dorsi and subcutaneous fat were quickly collected and immediately snap-frozen in liquid nitrogen, and stored at -80 °C for later use.
Library preparation and RNA sequencing
Total RNA was isolated and purified using the Trizol reagent (Takara, Japan), following the manufacturer’s instructions. The concentration and purity of isolated RNA samples were tested by electrophoresis and further examined with a nucleic acid analysis system (Experion System, Bio-Rad, USA). The A260/A280 ratio was between 1.9 and 2.1, indicating that the extracted RNA samples were of good quality. Sequencing was done by Beijing Novogene, using the Illumina HiSeq 2500 sequencing platform for two-end (pair-end), high-throughput sequencing (2 × 150 bp), and the amount of sequencing data (clean reads) per library was about 10 Gb.
Identification of differentially expressed circRNAs
To obtain high quality clean reads for subsequent bioinformatic analysis, raw data were filtered to remove low quality reads, adapter contamination, high unknown base content, and all rRNA-mapping reads. The Q20 and Q30 scores of the clean reads were calculated simultaneously, and circRNAs were then identified by aligning clean reads to the pig reference genome (Sscrofa 11.1) using find_circ and CIRI. The true expression levels of junction reads and circRNAs extracted from the sequencing data were affected by other factors, such as sequencing depth. Therefore, in order to accurately estimate true circRNA expression levels, the readcount data for predicted circRNAs were corrected using the TPM method. The filtered sequences were compared with the pig reference genome, the genomic features of each identified circRNA were analyzed, then the wider genomic and sequence characteristics of circRNAs in Jinfen White pigs were analyzed; mainly including the distribution among chromosomes, genomic location, and exon number of circRNAs.
Gene Ontology (GO) and pathway enrichment analysis (KEGG)
Using GO analysis (http://www.geneontology.org/) and KEGG pathway enrichment analysis (https://www.kegg.jp/), the host genes of DE circRNAs at the three growth and developmental stages were categorized in order to understand their biological functions. In the GO analysis, using the pig reference genome as background, the host genes were mapped to different functional entries. Thereafter, the enrichment significance of each entry was expressed as the corrected q-value, where entries with P < 0.05 were considered to be significantly enriched. In the KEGG pathway enrichment analysis, the hypergeometric test was applied—again using the pig reference genome as the analysis background—to screen the significantly enriched pathways (P < 0.05) for the host genes, with the results plotted as a bubble chart.
Target miRNA and gene prediction, and network analysis
To reveal functional interactions between ncRNAs and mRNAs, we used circRNA–miRNA interaction network analysis to construct a competing endogenous RNA network. Here, the public project TargetScan 7.0 (http://www.microrna.org/microrna/home.do) and miRanda software were used to predict potential regulatory relationships between circRNA, mRNA, and miRNA. Cytoscape (http://cytoscape.org/) was used for network visualization.
Quantitative Reverse Transcription Polymerase Chain Reaction (qRT-PCR)
Total RNA was extracted from the nine muscle samples used in the RNA-seq experiments, then reverse transcribed using PrimeScript ® RT Master Mix (Takara, Japan) with divergent primers designed using Premier 3.0 (including qRT-PCR primers and PCR primer-specific alignment by NCBI Primer-BLAST). Details of primers used in this study are listed in Supplementary Table S3. The expression of each circRNA, as well as other parameters, was detected using a CFX96 quantitative PCR instrument (Bio-Rad). The relative expression of target circRNAs was calculated using the 2−ΔΔCt method.
Circularization site verification of circRNA
Eight selected circRNAs were amplified by PCR using the designed qRT-PCR primers and the longissimus dorsi cDNA as template. Then, 2% agarose gel electrophoresis and Sanger sequencing of the PCR products were used to confirm reverse splicing of each circRNA by identifying the junction sequence.
RNase R enzyme digestion test
To test the circular nature of selected circRNAs, total RNA from porcine longissimus dorsi tissue was treated with RNase R enzyme (Epicentre) followed by qRT-PCR analysis.
Data processing and production of graphics were carried out using SPSS and GraphPad Prism 6 software. To test for statistically significant differences between two averages, data were analyzed by Student’s t test, and P < 0.05 was considered to indicate statistical significance .
Availability of data and materials
The circRNA, miRNA and mRNA raw sequence data reported in this paper have been uploaded to NCBI and are publicly accessible at PRJNA867525.
Kyoto Encyclopedia of Genes and Genomes
Competitive endogenous RNA
Long non-coding RNAs
Pig pluripotent stem cells
Myoblast determination protein 1
Myogenic factor 5
Myogenic factor 6
Paired box 3
Paired box 7
Myocyte enhancer factor 2
Tubulin tyrosine ligase like 7
Glutamate decarboxylase-like protein 1
Lim-domain only protein 7
Polypyrimidine tract binding protein 1
Transforming growth factor-β receptor type 1
Orthodenticle homeobox 2
Horak M, Novak J, Bienertova-Vasku J. Muscle-specific microRNAs in skeletal muscle development. Dev Biol. 2016;410(1):1–13.
Yin H, Price F, Rudnicki M. Satellite cells and the muscle stem cell niche. Physiol Rev. 2013;93(1):23–67.
Kuang S, Kuroda K, Le Grand F, et al. Asymmetric self-renewal and commitment of satellite stem cells in muscle. Cell. 2007;129(5):999–1010.
Verma N, Rettenmeier A, Schmitz-Spanke S. Recent advances in the use of Sus scrofa (pig) as a model system for proteomic studies. Proteomics. 2011;11(4):776–93.
Sanger H, Klotz G, Riesner D, et al. Viroids are single-stranded covalently closed circular RNA molecules existing as highly base-paired rod-like structures. Proc Natl Acad Sci USA. 1976;73(11):3852–6.
Brand-Saberi B. Genetic and epigenetic control of skeletal muscle development. Ann Anat. 2005;187(3):199–207.
Buckingham M, Relaix F. The role of Pax genes in the development of tissues and organs: Pax3 and Pax7 regulate muscle progenitor cell functions. Annu Rev Cell Dev Biol. 2007;23:645–73.
Sabourin L, Rudnicki M. The molecular regulation of myogenesis. Clin Genet. 2000;57(1):16–25.
Bi P, Ramirez-Martinez A, Li H, et al. Control of muscle formation by the fusogenic micropeptide myomixer. Science. 2017;356(6335):323–7.
Ito Y, Kayama T, Asahara H. A systems approach and skeletal myogenesis. Comp Funct Genomics. 2012;2012: 759407.
Rudnicki M, Schnegelsberg P, Stead R, et al. MyoD or Myf-5 is required for the formation of skeletal muscle. Cell. 1993;75(7):1351–9.
Kassar-Duchossoy L, Gayraud-Morel B, Gomès D, et al. Mrf4 determines skeletal muscle identity in Myf5: Myod double-mutant mice. Nature. 2004;431(7007):466–71.
Hasty P, Bradley A, Morris J, et al. Muscle deficiency and neonatal death in mice with a targeted mutation in the myogenin gene. Nature. 1993;364(6437):501–6.
Myer A, Olson E, Klein W. MyoD cannot compensate for the absence of myogenin during skeletal muscle differentiation in murine embryonic stem cells. Dev Biol. 2001;229(2):340–50.
Li J, Zhao W, Li Q, et al. Long Non-Coding RNA H19 promotes porcine satellite cell differentiation by interacting with TDP43. Genes. 2020;11(3):259–71.
Geng T, Liu Y, Xu Y, et al. H19 lncRNA promotes skeletal muscle insulin sensitivity in part by targeting AMPK. Diabetes. 2018;67(11):2183–98.
Raza S, Kaster N, Khan R, et al. The role of microRNAs in muscle tissue development in beef cattle. Genes. 2020;11(3):295–308.
Vicens Q, Westhof E. Biogenesis of circular RNAs. Cell. 2014;159(1):13–4.
Yuan X, Diao J, Du A, et al. Circular RNA expression profiles and features in NAFLD mice: a study using RNA-seq data. J Transl Med. 2020;18(1):476–88.
Hao Z, Zhou H, Hickford J, et al. Identification and characterization of circular RNA in lactating mammary glands from two breeds of sheep with different milk production profiles using RNA-Seq. Genomics. 2020;112(3):2186–93.
Li M, Xie X, Zhou J, et al. Quantifying circular RNA expression from RNA-seq data using model-based framework. Bioinformatics. 2017;33(14):2131–9.
Li X, Yang L, Chen L. The Biogenesis, functions, and challenges of circular RNAs. Mol Cell. 2018;71(3):428–42.
Memczak S, Jens M, Elefsinioti A, et al. Circular RNAs are a large class of animal RNAs with regulatory potency. Nature. 2013;495(7441):333–8.
Tao M, Zheng M, Xu Y, et al. CircRNAs and their regulatory roles in cancers. Mol Med. 2021;27(1):94–116.
Shi X, Yang J, Liu M, et al. Circular RNA ANAPC7 inhibits tumor growth and muscle wasting via PHLPP2-AKT-TGF-β signaling axis in pancreatic cancer. Gastroenterology. 2022;162:2004–17.
Wang Y, Yan Q, Mo Y, et al. Splicing factor derived circular RNA circCAMSAP1 accelerates nasopharyngeal carcinoma tumorigenesis via a SERPINH1/c-Myc positive feedback loop. Mol Cancer. 2022;21(1):62–82.
Fan H, Chen Z, Chen X, et al. METTL14-mediated mA modification of circORC5 suppresses gastric cancer progression by regulating miR-30c-2-3p/AKT1S1 axis. Mol Cancer. 2022;21(1):51–65.
Ouyang H, Chen X, Wang Z, et al. Circular RNAs are abundant and dynamically expressed during embryonic muscle development in chickens. DNA Res. 2018;25(1):71–86.
Chen X, Ouyang H, Wang Z, et al. A Novel Circular RNA generated by FGFR2 gene promotes myoblast proliferation and differentiation by sponging miR-133a-5p and miR-29b-1-5p. Cells. 2018;7(11):7110199.
Yin H, Shen X, Zhao J, et al. Circular RNA CircFAM188B encodes a protein that regulates proliferation and differentiation of chicken skeletal muscle satellite cells. Front Cell Dev Biol. 2020;8.
Li H, Yang J, Wei X, et al. CircFUT10 reduces proliferation and facilitates differentiation of myoblasts by sponging miR-133a. J Cell Physiol. 2018;233(6):4643–51.
Sun D, An J, Cui Z, et al. CircCSDE1 regulates proliferation and differentiation of c2c12 myoblasts by sponging miR-21–3p. Int J Mol Sci, 2022, 23(19):12038.
Ma L, Chen W, Li S, et al. Identification and functional prediction of circular RNAs related to growth traits and skeletal muscle development in duroc pigs. Front Genet. 2022;13: 858763.
Yan J, Yang Y, Fan X, et al. CircRNAome profiling reveals circFgfr2 regulates myogenesis and muscle regeneration via a feedback loop. J Cachexia Sarcopenia Muscle. 2022;13(1):696–712.
Cai B, Ma M, Zhou Z, et al. CircPTPN4 regulates myogenesis via the miR-499-3p/NAMPT axis. J Anim Sci Biotechnol. 2022;13(1):2–17.
Shen X, Wei Y, Liu W, et al. A novel circular RNA circITSN2 targets the miR-218-5p/LMO7 axis to promote chicken embryonic myoblast proliferation and differentiation. Front Cell Dev Biol. 2021;9.
Legnini I, Di Timoteo G, Rossi F, et al. Circ-ZNF609 is a circular RNA that can be translated and functions in myogenesis. Mol Cell. 2017;66(1):22–37.
Liang G, Yang Y, Niu G, et al. Genome-wide profiling of Sus scrofa circular RNAs across nine organs and three developmental stages. DNA Res. 2017;24(5):523–35.
Kristensen L, Andersen M, Stagsted L, et al. The biogenesis, biology and characterization of circular RNAs. Nat Rev Genet. 2019;20(11):675–91.
Das A, Shyamal S, Sinha T, et al. Identification of potential circRNA-microRNA-mRNA regulatory network in skeletal muscle. Front Mol Biosci. 2021;8.
Zhang Y, Zhang X, Chen T, et al. Circular intronic long noncoding RNAs. Mol Cell. 2013;51(6):792–806.
Li Z, Huang C, Bao C, et al. Exon-intron circular RNAs regulate transcription in the nucleus. Nat Struct Mol Biol. 2015;22(3):256–64.
Kanehisa M, Furumichi M, Sato Y, et al. KEGG for taxonomy-based analysis of pathways and genomes. Nucleic Acids Res. 2023;51:D587–92.
Chen B, Yu J, Guo L, et al. Circular RNA circHIPK3 promotes the proliferation and differentiation of chicken myoblast cells by sponging miR-30a-3p. Cells. 2019;8(2):177–86.
Cheng X, Ai K, Yi L, et al. The mmu_circRNA_37492/hsa_circ_0012138 function as potential ceRNA to attenuate obstructive renal fibrosis. Cell Death Dis. 2022;13(3):207–18.
Deng W, Zhou X, Zhu K, et al. Novel circular RNA circ_0086722 drives tumor progression by regulating the miR-339-5p/STAT5A axis in prostate cancer. Cancer Lett. 2022;533: 215606.
Wu X, Zhou J, Zhao L, et al. CircCYP24A1 hampered malignant phenotype of renal cancer carcinoma through modulating CMTM-4 expression via sponging miR-421. Cell Death Dis. 2022;13(2):190–291.
Shen X, Liu Z, Cao X, et al. Circular RNA profiling identified an abundant circular RNA circTMTC1 that inhibits chicken skeletal muscle satellite cell differentiation by sponging miR-128-3p. Int J Biol Sci. 2019;15(10):2265–81.
Li B, Yang J, He J, et al. Spatiotemporal regulation and functional analysis of circular RNAs in skeletal muscle and subcutaneous fat during pig growth. Biology. 2021;10(9):841–62.
Li M, Zhang N, Zhang W, et al. Comprehensive analysis of differentially expressed circRNAs and ceRNA regulatory network in porcine skeletal muscle. BMC Genomics. 2021;22(1):320–32.
Qi K, Liu Y, Li C, et al. Construction of circRNA-related ceRNA networks in longissimus dorsi muscle of Queshan black and large white pigs. Mol Genet Genomics. 2022;297(1):101–12.
Jakobi T, Czaja-Hasse L, Reinhardt R, et al. Profiling and validation of the circular RNA repertoire in adult murine hearts. Genomics Proteomics Bioinformatics. 2016;14:216–23.
Yang Xinran, Ma X, Mei C, et al. A genome-wide landscape of mRNAs, lncRNAs, circRNAs and miRNAs during intramuscular adipogenesis in cattle. BMC Genomics, 2022, 23: 691.
Wang J, Chen H, Zhang Y, et al. Comprehensive analysis of differentially expressed circRNAs in the ovaries of low- and high-fertility sheep. Animals, 2023, 13:236.
Zhang G, Zhang J, Wu P, et al. Transcriptome sequencing analysis of circRNA in skeletal muscle between fast- and slow-growing chickens at embryonic stages. Animals. 2022;12(22):3166.
Chen X, Wang Z, Chen Y, et al. Transcriptome analysis of differentially expressed circRNAs miRNAs and mRNAs during the challenge of coccidiosis. Front Immunol. 2022;13: 910860.
Barrett SP, Salzman J. Circular RNAs: analysis, expression and potential functions. Development. 2016;143(11):1838–47.
Witczak C, Sharoff C, Goodyear L. AMP-activated protein kinase in skeletal muscle: from structure and localization to its role as a master regulator of cellular metabolism. Cell Mol Life Sci. 2008;65(23):3737–55.
Lee Y, Yun M, Kim H, et al. Exogenous administration of DLK1 ameliorates hepatic steatosis and regulates gluconeogenesis via activation of AMPK. Int J Obesity. 2016;40(2):356–65.
Xu M, Chen X, Huang Z, et al. Procyanidin B2 promotes skeletal slow-twitch myofiber gene expression through the ampk signaling pathway in C2C12 myotubes. J Agric Food Chem. 2020;68(5):1306–14.
Lauterbach N, Gonçalves D, Silveira W, et al. Urocortin 2 promotes hypertrophy and enhances skeletal muscle function through cAMP and insulin/IGF-1 signaling pathways. Molecular Metabolism. 2022;60: 101492.
Machado J, Manfredi L, Silveira W, et al. Calcitonin gene-related peptide inhibits autophagic-lysosomal proteolysis through cAMP/PKA signaling in rat skeletal muscles. Int J Biochem Cell Biol. 2016;72:40–50.
Tong T, Kim M, Park T. α-Ionone attenuates high-fat diet-induced skeletal muscle wasting in mice via activation of cAMP signaling. Food Function. 2019;10(2):1167–78.
Dedeic Z, Cetera M, Cohen T, et al. Emerin inhibits Lmo7 binding to the Pax3 and MyoD promoters and expression of myoblast proliferation genes. J Cell Sci. 2011;124:1691–702.
Wei X, Li H, Yang J, et al. Circular RNA profiling reveals an abundant circLMO7 that regulates myoblasts differentiation and survival by sponging miR-378a-3p. Cell Death Dis. 2017;8(10): e3153.
Zhang P, Du J, Guo X, et al. LncMyoD promotes skeletal myogenesis and regulates skeletal muscle fiber-type composition by sponging miR-370-3p. Genes. 2021;12(4):589–95.
Du Y, Wang Y, Li Y, et al. MiR-214-5p regulating differentiation of intramuscular preadipocytes in goats via targeting KLF12. Front Genet. 2021;12:748629.
Feng Y, Niu L, Wei W, et al. A feedback circuit between miR-133 and the ERK1/2 pathway involving an exquisite mechanism for regulating myoblast proliferation and differentiation. Cell Death Dis. 2013;4: e934.
Elsaeid Elnour I, Dong D, Wang X, et al. Bta-miR-885 promotes proliferation and inhibits differentiation of myoblasts by targeting MyoD1. J Cell Physiol. 2020;235(10):6625–36.
Zhu W, Zhou B, Rong L, et al. Roles of PTBP1 in alternative splicing, glycolysis, and oncogenesis. J Zhejiang Univ Sci B. 2020;21(2):122–36.
Harada N, Gotoda Y, Hatakeyama A, et al. Differential regulation of Actn2 and Actn3 expression during unfolded protein response in C2C12 myotubes. J Muscle Res Cell Motility. 2020;41:199–209.
Hu H, Fu Y, Zhou B, et al. Long non-coding RNA TCONS_00814106 regulates porcine granulosa cell proliferation and apoptosis by sponging miR-1343. Mol Cell Endocrinol. 2021;520:111064.
Xie Y, Cao H, Zhang Z, et al. Molecular network of miR-1343 regulates the pluripotency of porcine pluripotent stem cells via repressing OTX2 expression. RNA Biol. 2019;16(1):82–92.
Yu J, Shao S, Xiong Y. Molecular characterization, expression patterns, and promoter activity analysis of PGM1 in pigs. Genet Mol Res. 2015;14(1):2816–24.
Xing K, Wang K, Ao H, et al. Comparative adipose transcriptome analysis digs out genes related to fat deposition in two pig breeds. Sci Rep. 2019;9(1):12925.
Kim G, Jeong J, Yang H, et al. Differential abundance of proteome associated with intramuscular variation of meat quality in porcine longissimus thoracis et lumborum muscle. Meat Sci. 2019;149:85–95.
We thank Novogene for providing helpful technical support in bioinformatics analysis. We also thank all the staff of the Datong Pig Breeding Farm.
This research was jointly supported by the National Natural Science Foundation of China (31872336), Special Funds for Scholars Support Program of Shanxi Province (2016; 2017), the Basic Research Project of Shanxi Province (20210302124639), Science and Technology Innovation Project of Shanxi Agricultural University (2020BQ56), Shanxi Province Excellent Doctoral Work Award-Scientific Research Project (SXBYKY2021026), University Science and Technology Innovation Project of Shanxi Province (2021L158).
Ethics approval and consent to participate
The nine pigs used in this study were all from the Datong Pig Breeding Farm (Shanxi, China). We obtained the written consent of the director of Datong Pig Breeding Farm and provided a proof. All methods are reported in accordance with ARRIVE guidelines (https://arriveguidelines.org) for the reporting of animal experiments, and all experimental procedures were carried out in accordance with the guidelines of the Regional Animal Ethics Committee and approved by the Animal Ethics Committee of Shanxi Agricultural University (Approval No. SXAU-EAW-2021Sus.AB.0926001).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Yun, J., Huang, X., Liu, C. et al. Genome-wide analysis of circular RNA-mediated ceRNA regulation in porcine skeletal muscle development. BMC Genomics 24, 196 (2023). https://doi.org/10.1186/s12864-023-09284-7
- Skeletal muscle