Data describing the effects of Induced knockout of Heterogeneous nuclear ribonucleoprotein K in mouse skeletal muscle satellite cells.

HnRNPK, a prominent member of the heterogeneous nuclear ribonucleoprotein (hnRNP) family, is widely expressed in mammalian tissues and plays a crucial role in animal development. Despite its well-established functions, limited information is available regarding its role in skeletal muscle development and regeneration. To elucidate the functional role of hnRNPK in skeletal muscle, we utilized Pax7CreER; HnrnpkLoxP/LoxP (Hnrnpk pKO) mice as a model, isolated primary mouse skeletal muscle satellite cells (MuSCs), and induced hnRNPK knockout using 4-OTH. Transcriptome sequencing was performed on four distinct groups: Hnrnpk pKO MuSCs undergoing proliferation for 24 h (ethanol 24 h) and 48 h (ethanol 48 h) after treatment with ethanol as the control, as well as Hnrnpk pKO MuSCs undergoing proliferation for 24 h (4-OHT 24 h) and 48 h (4-OHT 48 h) after treatment with 4-OHT as the hnRNPK-induced knockout group. The RNA sequencing data was generated using the Illumina HiSeq 2000/2500 sequencing platform. The raw data files have been archived in the Sequence Read Archive at the China National Center for Bioinformation (CNCB) under the accession number CRA015864. This data article is related to the research paper "Deletion of heterogeneous nuclear ribonucleoprotein K in satellite cells leads to inhibited skeletal muscle regeneration in mice, Genes & Diseases 11: 101,062, DOI: 10.1016/j.gendis.2023.06.031".

RNA-Binding Proteins in the Regulation of Adipogenesis and Adipose Function.

The obesity epidemic represents a critical public health issue worldwide, as it is a vital risk factor for many diseases, including type 2 diabetes (T2D) and cardiovascular disease. Obesity is a complex disease involving excessive fat accumulation. Proper adipose tissue accumulation and function are highly transcriptional and regulated by many genes. Recent studies have discovered that post-transcriptional regulation, mainly mediated by RNA-binding proteins (RBPs), also plays a crucial role. In the lifetime of RNA, it is bound by various RBPs that determine every step of RNA metabolism, from RNA processing to alternative splicing, nucleus export, rate of translation, and finally decay. In humans, it is predicted that RBPs account for more than 10% of proteins based on the presence of RNA-binding domains. However, only very few RBPs have been studied in adipose tissue. The primary aim of this paper is to provide an overview of RBPs in adipogenesis and adipose function. Specifically, the following best-characterized RBPs will be discussed, including HuR, PSPC1, Sam68, RBM4, Ybx1, Ybx2, IGF2BP2, and KSRP. Characterization of these proteins will increase our understanding of the regulatory mechanisms of RBPs in adipogenesis and provide clues for the etiology and pathology of adipose-tissue-related diseases.

ITGß6 Facilitates Skeletal Muscle Development by Maintaining the Properties and Cytoskeleton Stability of Satellite Cells.

Integrin proteins are important receptors connecting the intracellular skeleton of satellite cells and the extracellular matrix (ECM), playing an important role in the process of skeletal muscle development. In this research, the function of ITGß6 in regulating the differentiation of satellite cells was studied. Transcriptome and proteome analysis indicated that Itgß6 is a key node connecting ECM-related proteins to the cytoskeleton, and it is necessary for the integrity of the membrane structure and stability of the cytoskeletal system, which are essential for satellite cell adhesion. Functional analysis revealed that the ITGß6 protein could affect the myogenic differentiation potential of satellite cells by regulating the expression of PAX7 protein, thus regulating the formation of myotubes. Moreover, ITGß6 is involved in muscle development by regulating cell-adhesion-related proteins, such as ß-laminin, and cytoskeletal proteins such as PXN, DMD, and VCL. In conclusion, the effect of ITGß6 on satellite cell differentiation mainly occurs before the initiation of differentiation, and it regulates terminal differentiation by affecting satellite cell characteristics, cell adhesion, and the stability of the cytoskeleton system.

Deletion of Hnrnpk Gene Causes Infertility in Male Mice by Disrupting Spermatogenesis.

HnRNPK is a heterogeneous nuclear ribonucleoprotein (hnRNP) that has been firmly implicated in transcriptional and post-transcriptional regulation. However, the molecular mechanisms by which hnRNPK orchestrates transcriptional or post-transcriptional regulation are not well understood due to early embryonic lethality in homozygous knockout mice, especially in a tissue-specific context. Strikingly, in this study, we demonstrated that hnRNPK is strongly expressed in the mouse testis and mainly localizes to the nucleus in spermatogonia, spermatocytes, and round spermatids, suggesting an important role for hnRNPK in spermatogenesis. Using a male germ cell-specific hnRNPK-depleted mouse model, we found that it is critical for testicular development and male fertility. The initiation of meiosis of following spermatogenesis was not affected in Hnrnpk cKO mice, while most germ cells were arrested at the pachytene stage of the meiosis and no mature sperm were detected in epididymides. The further RNA-seq analysis of Hnrnpk cKO mice testis revealed that the deletion of hnRNPK disturbed the expression of genes involved in male reproductive development, among which the meiosis genes were significantly affected, and Hnrnpk cKO spermatocytes failed to complete the meiotic prophase. Together, these results identify hnRNPK as an essential regulator of spermatogenesis and male fertility.

RNA-seq profiling of white and brown adipocyte differentiation treated with epigallocatechin gallate.

Due to serious adverse effects, many of the approved anti-obesity medicines have been withdrawn, and the selection of safer natural ingredients is of great interest. Epigallocatechin gallate (EGCG) is one of the major green tea catechins, and has been demonstrated to possess an anti-obesity function by regulating both white and brown adipose tissue activity. However, there are currently no publicly available studies describing the effects of EGCG on the two distinct adipose tissue transcriptomes. The stromal vascular fraction (SVF) cell derived from adipose tissue is a classic cell model for studying adipogenesis and fat accumulation. In the current study, primary WAT and BAT SVF cells were isolated and induced to adipogenic differentiation in the presence or absence of EGCG. RNA-seq was used to determine genes regulated by EGCG and identify the key differences between the two functionally distinct adipose tissues. Taken together, we provide detailed stage- and tissue-specific gene expression profiles affected by EGCG. These data will be valuable for obesity-related clinical/basic research.

Assessment of CircRNA Expression Profiles and Potential Functions in Brown Adipogenesis.

Brown adipose tissue (BAT) is specialized for energy expenditure, thus a better understanding of the regulators influencing BAT development could provide novel strategies to defense obesity. Many protein-coding genes, miRNAs, and lncRNAs have been investigated in BAT development, however, the expression patterns and functions of circRNA in brown adipogenesis have not been reported yet. This study determined the circRNA expression profiles across brown adipogenesis (proliferation, early differentiated, and fully differentiated stages) by RNA-seq. We identified 3,869 circRNAs and 36.9% of them were novel. We found the biogenesis of circRNA was significantly related to linear mRNA transcription, meanwhile, almost 70% of circRNAs were generated by alternative back-splicing. Next, we examined the cell-specific and differentiation stage-specific expression of circRNAs. Compared to white adipocytes, nearly 30% of them were specifically expressed in brown adipocytes. Further, time-series expression analysis showed circRNAs were dynamically expressed, and 117 differential expression circRNAs (DECs) in brown adipogenesis were identified, with 77 upregulated and 40 downregulated. Experimental validation showed the identified circRNAs could be successfully amplified and the expression levels detected by RNA-seq were reliable. For the potential functions of the circRNAs, GO analysis suggested that the decreased circRNAs were enriched in cell proliferation terms, while the increased circRNAs were enriched in development and thermogenic terms. Bioinformatics predictions showed that DECs contained numerous binding sites of functional miRNAs. More interestingly, most of the circRNAs contained multiple binding sites for the same miRNA, indicating that they may facilitate functions by acting as microRNA sponges. Collectively, we characterized the circRNA expression profiles during brown adipogenesis and provide numerous novel circRNAs candidates for future brown adipogenesis regulating studies.

Identification of Circular RNA Expression Profiles in White Adipocytes and Their Roles in Adipogenesis.

Obesity and its related metabolic diseases have become great public health threats worldwide. Although accumulated evidence suggests that circRNA is a new type of non-coding RNAs regulating various physiological and pathological processes, little attention has been paid to the expression profiles and functions of circRNAs in white adipose tissue. In this study, 3,771 circRNAs were detected in three stages of white adipogenesis (preadipocyte, differentiating preadipocyte, and mature adipocyte) by RNA-seq. Experimental validation suggested that the RNA-seq results are highly reliable. We found that nearly 10% of genes which expressed linear RNAs in adipocytes could also generate circRNAs. In addition, 40% of them produced multiple circRNA isoforms. We performed correlation analysis and found that a great deal of circRNAs (nearly 50%) and their parental genes were highly correlated in expression levels. A total of 41 differential expression circRNAs (DECs) were detected during adipogenesis and an extremely high ratio of them (80%) were correlated with their parental genes, indicating these circRNAs may potentially play roles in regulating the expression of their parental genes. KEGG enrichment and GO annotation of the parental genes suggesting that the DECs may participate in several adipogenesis-related pathways. Following rigorous selection, we found that many up-regulated circRNAs contain multiple miRNAs binding sites, such as miR17, miR-30c, and miR-130, indicating they may potentially facilitate their regulatory functions by acting as miRNA sponges. These results suggest that plenty of circRNAs are expressed in white adipogenesis and the DECs may serve as new candidates for future adipogenesis regulation.

Role of juvenile hormone receptor Methoprene-tolerant 1 in silkworm larval brain development and domestication.

The insect brain is the central part of the neurosecretory system, which controls morphology, physiology, and behavior during the insect's lifecycle. Lepidoptera are holometabolous insects, and their brains develop during the larval period and metamorphosis into the adult form. As the only fully domesticated insect, the Lepidoptera silkworm Bombyx mori experienced changes in larval brain morphology and certain behaviors during the domestication process. Hormonal regulation in insects is a key factor in multiple processes. However, how juvenile hormone (JH) signals regulate brain development in Lepidoptera species, especially in the larval stage, remains elusive. We recently identified the JH receptor Methoprene tolerant 1 ( Met1) as a putative domestication gene. How artificial selection on Met1 impacts brain and behavioral domestication is another important issue addressing Darwin's theory on domestication. Here, CRISPR/Cas9-mediated knockout of Bombyx Met1 caused developmental retardation in the brain, unlike precocious pupation of the cuticle. At the whole transcriptome level, the ecdysteroid (20-hydroxyecdysone, 20E) signaling and downstream pathways were overactivated in the mutant cuticle but not in the brain. Pathways related to cell proliferation and specialization processes, such as extracellular matrix (ECM)-receptor interaction and tyrosine metabolism pathways, were suppressed in the brain. Molecular evolutionary analysis and in vitro assay identified an amino acid replacement located in a novel motif under positive selection in B. mori, which decreased transcriptional binding activity. The B. mori MET1 protein showed a changed structure and dynamic features, as well as a weakened co-expression gene network, compared with B. mandarina. Based on comparative transcriptomic analyses, we proposed a pathway downstream of JH signaling (i.e., tyrosine metabolism pathway) that likely contributed to silkworm larval brain development and domestication and highlighted the importance of the biogenic amine system in larval evolution during silkworm domestication.

Comparative Proteomics and Phosphoproteomics Analysis Reveal the Possible Breed Difference in Yorkshire and Duroc Boar Spermatozoa.

Sperm cells are of unique elongated structure and function, the development of which is tightly regulated by the existing proteins and the posttranslational modifications (PTM) of these proteins. Based on the phylogenetic relationships of various swine breeds, Yorkshire boar is believed to be distinctly different from Duroc boar. The comprehensive differential proteomics and phosphoproteomics profilings were performed on spermatozoa from both Yorkshire and Duroc boars. By both peptide and PTM peptide quantification followed by statistical analyses, 167 differentially expressed proteins were identified from 1,745 proteins, and 283 differentially expressed phosphopeptides corresponding to 102 unique differentially phosphorylated proteins were measured from 1,140 identified phosphopeptides derived from 363 phosphorylated proteins. The representative results were validated by Western blots. Pathway enrichment analyses revealed that majority of differential expression proteins and differential phosphorylation proteins were primarily concerned with spermatogenesis, male gamete generation, sperm motility, energy metabolism, cilium morphogenesis, axonemal dynein complex assembly, sperm-egg recognition, and capacitation. Remarkably, axonemal dynein complex assembly related proteins, such as SMCP, SUN5, ODF1, AKAP3, and AKAP4 that play a key regulatory role in the sperm physiological functions, were significantly higher in Duroc spermatozoa than that of Yorkshire. Furthermore, phosphorylation of sperm-specific proteins, such as CABYR, ROPN1, CALM1, PRKAR2A, and PRKAR1A, participates in regulation of the boar sperm motility mainly through the cAMP/PKA signal pathway in different breeds, demonstrating that protein phosphorylation may be an important mechanism underlying the sperm diversity. Protein-protein interaction analysis revealed that the 14 overlapped proteins between differential expression proteins and differential phosphorylation proteins potentially played a key role in sperm development and motility of the flagellum, including the proteins ODF1, SMCP, AKAP4, FSIP2, and SUN5. Taken together, these physiologically and functionally differentially expressed proteins (DEPs) and differentially expressed phosphorylated proteins (DPPs) may constitute the proteomic backgrounds between the two different boar breeds. The validation will be performed to delineate the roles of these PTM proteins as modulators of Yorkshire and Duroc boar spermatozoa.

Genome-wide analysis of DNA G-quadruplex motifs across 37 species provides insights into G4 evolution.

G-quadruplex (G4) structures have been predicted in the genomes of many organisms and proven to play regulatory roles in diverse cellular activities. However, there is little information on the evolutionary history and distribution characteristics of G4s. Here, whole-genome characteristics of potential G4s were studied in 37 evolutionarily representative species. During evolution, the number, length, and density of G4s generally increased. Immunofluorescence in seven species confirmed G4s' presence and evolutionary pattern. G4s tended to cluster in chromosomes and were enriched in genetic regions. Short-loop G4s were conserved in most species, while loop-length diversity also existed, especially in mammals. The proportion of G4-bearing genes and orthologue genes, which appeared to be increasingly enriched in transcription factors, gradually increased. The antagonistic relationship between G4s and DNA methylation sites was detected. These findings imply that organisms may have evolutionarily developed G4 into a novel reversible and elaborate transcriptional regulatory mechanism benefiting multiple physiological activities of higher organisms.

Circular RNA Regulation of Myogenesis.

Circular RNA (circRNA) is a novel class of non-coding RNA generated by pre-mRNA back splicing, which is characterized by a closed-loop structure. Although circRNAs were firstly reported decades ago, their regulatory roles have not been discovered until recently. In this review, we discussed the putative biogenesis pathways and regulatory functions of circRNAs. Recent studies showed that circRNAs are abundant in skeletal muscle tissue, and their expression levels are regulated during muscle development and aging. We, thus, characterized the expression profile of circRNAs in skeletal muscle and discussed regulatory functions and mechanism-of-action of specific circRNAs in myogenesis. The future investigation into the roles of circRNAs in both physiological and pathological conditions may provide novel insights in skeletal muscle development and provide new therapeutic strategies for muscular diseases.

Post-translational modification control of RNA-binding protein hnRNPK function.

Heterogeneous nuclear ribonucleoprotein K (hnRNPK), a ubiquitously occurring RNA-binding protein (RBP), can interact with numerous nucleic acids and various proteins and is involved in a number of cellular functions including transcription, translation, splicing, chromatin remodelling, etc. Through its abundant biological functions, hnRNPK has been implicated in cellular events including proliferation, differentiation, apoptosis, DNA damage repair and the stress and immune responses. Thus, it is critical to understand the mechanism of hnRNPK regulation and its downstream effects on cancer and other diseases. A number of recent studies have highlighted that several post-translational modifications (PTMs) possibly play an important role in modulating hnRNPK function. Phosphorylation is the most widely occurring PTM in hnRNPK. For example, in vivo analyses of sites such as S116 and S284 illustrate the purpose of PTM of hnRNPK in altering its subcellular localization and its ability to bind target nucleic acids or proteins. Other PTMs such as methylation, ubiquitination, sumoylation, glycosylation and proteolytic cleavage are increasingly implicated in the regulation of DNA repair, cellular stresses and tumour growth. In this review, we describe the PTMs that impact upon hnRNPK function on gene expression programmes and different disease states. This knowledge is key in allowing us to better understand the mechanism of hnRNPK regulation.

New Insights into the Interplay between Non-Coding RNAs and RNA-Binding Protein HnRNPK in Regulating Cellular Functions.

The emerging data indicates that non-coding RNAs (ncRNAs) epresent more than the "junk sequences" of the genome. Both miRNAs and long non-coding RNAs (lncRNAs) are involved in fundamental biological processes, and their deregulation may lead to oncogenesis and other diseases. As an important RNA-binding protein (RBP), heterogeneous nuclear ribonucleoprotein K (hnRNPK) is known to regulate gene expression through the RNA-binding domain involved in various pathways, such as transcription, splicing, and translation. HnRNPK is a highly conserved gene that is abundantly expressed in mammalian cells. The interaction of hnRNPK and ncRNAs defines the novel way through which ncRNAs affect the expression of protein-coding genes and form autoregulatory feedback loops. This review summarizes the interactions of hnRNPK and ncRNAs in regulating gene expression at transcriptional and post-transcriptional levels or by changing the genomic structure, highlighting their involvement in carcinogenesis, glucose metabolism, stem cell differentiation, virus infection and other cellular functions. Drawing connections between such discoveries might provide novel targets to control the biological outputs of cells in response to different stimuli.

Artificial selection on storage protein 1 possibly contributes to increase of hatchability during silkworm domestication.

Like other domesticates, the efficient utilization of nitrogen resources is also important for the only fully domesticated insect, the silkworm. Deciphering the way in which artificial selection acts on the silkworm genome to improve the utilization of nitrogen resources and to advance human-favored domestication traits, will provide clues from a unique insect model for understanding the general rules of Darwin's evolutionary theory on domestication. Storage proteins (SPs), which belong to a hemocyanin superfamily, basically serve as a source of amino acids and nitrogen during metamorphosis and reproduction in insects. In this study, through blast searching on the silkworm genome and further screening of the artificial selection signature on silkworm SPs, we discovered a candidate domestication gene, i.e., the methionine-rich storage protein 1 (SP1), which is clearly divergent from other storage proteins and exhibits increased expression in the ova of domestic silkworms. Knockout of SP1 via the CRISPR/Cas9 technique resulted in a dramatic decrease in egg hatchability, without obvious impact on egg production, which was similar to the effect in the wild silkworm compared with the domestic type. Larval development and metamorphosis were not affected by SP1 knockout. Comprehensive ova comparative transcriptomes indicated significant higher expression of genes encoding vitellogenin, chorions, and structural components in the extracellular matrix (ECM)-interaction pathway, enzymes in folate biosynthesis, and notably hormone synthesis in the domestic silkworm, compared to both the SP1 mutant and the wild silkworm. Moreover, compared with the wild silkworms, the domestic one also showed generally up-regulated expression of genes enriched in the structural constituent of ribosome and amide, as well as peptide biosynthesis. This study exemplified a novel case in which artificial selection could act directly on nitrogen resource proteins, further affecting egg nutrients and eggshell formation possibly through a hormone signaling mediated regulatory network and the activation of ribosomes, resulting in improved biosynthesis and increased hatchability during domestication. These findings shed new light on both the understanding of artificial selection and silkworm breeding from the perspective of nitrogen and amino acid resources.

Genome-wide epigenetic landscape of pig lincRNAs and their evolution during porcine domestication.

AIM: We aimed to identify previously unreported long intergenic noncoding RNAs (lincRNAs) in the porcine liver, an important metabolic tissue, and further illustrate the epigenomic landscapes and the evolution of lincRNAs. MATERIALS & METHODS: We used porcine omics data and comprehensively analyzed and identified lincRNAs and their methylation, expression and evolutionary patterns during pig domestication. RESULTS: LincRNAs exhibit highly methylated promoter and downstream regions, as well as lower expression levels and higher tissue specificity than protein-coding genes. We identified a batch of lincRNAs with selection signals that are associated with pig domestication, which are more highly expressed in the liver than in other tissues (19:10/8/6/3/2/1/1). Interestingly, the lincRNA linc-sscg1779 and its target gene C6, which is crucial in liver metabolism, are differentially expressed during pig domestication. CONCLUSION: Although they may originate from noisy transcripts, lincRNAs may be subjected to artificial selection. This phenomenon implies the functional importance of lincRNAs in pig domestication.

A survey of transcriptome complexity in Sus scrofa using single-molecule long-read sequencing.

Alternative splicing (AS) and fusion transcripts produce a vast expansion of transcriptomes and proteomes diversity. However, the reliability of these events and the extend of epigenetic mechanisms have not been adequately addressed due to its limitation of uncertainties about the complete structure of mRNA. Here we combined single-molecule real-time sequencing, Illumina RNA-seq and DNA methylation data to characterize the landscapes of DNA methylation on AS, fusion isoforms formation and lncRNA feature and further to unveil the transcriptome complexity of pig. Our analysis identified an unprecedented scale of high-quality full-length isoforms with over 28,127 novel isoforms from 26,881 novel genes. More than 92,000 novel AS events were detected and intron retention predominated in AS model, followed by exon skipping. Interestingly, we found that DNA methylation played an important role in generating various AS isoforms by regulating splicing sites, promoter regions and first exons. Furthermore, we identified a large of fusion transcripts and novel lncRNAs, and found that DNA methylation of the promoter and gene body could regulate lncRNA expression. Our results significantly improved existed gene models of pig and unveiled that pig AS and epigenetic modify were more complex than previously thought.

Identification and Functional Analysis of Long Intergenic Non-coding RNAs Underlying Intramuscular Fat Content in Pigs.

Intramuscular fat (IMF) content is an important trait that can affect pork quality. Previous studies have identified many genes that can regulate IMF. Long intergenic non-coding RNAs (lincRNAs) are emerging as key regulators in various biological processes. However, lincRNAs related to IMF in pig are largely unknown, and the mechanisms by which they regulate IMF are yet to be elucidated. Here we reconstructed 105,687 transcripts and identified 1,032 lincRNAs in pig longissimus dorsi muscle (LDM) of four stages with different IMF contents based on published RNA-seq. These lincRNAs show typical characteristics such as shorter length and lower expression compared with protein-coding genes. Combined with methylation data, we found that both the promoter and genebody methylation of lincRNAs can negatively regulate lincRNA expression. We found that lincRNAs exhibit high correlation with their protein-coding neighbors in expression. Co-expression network analysis resulted in eight stage-specific modules, gene ontology and pathway analysis of them suggested that some lincRNAs were involved in IMF-related processes, such as fatty acid metabolism and peroxisome proliferator-activated receptor signaling pathway. Furthermore, we identified hub lincRNAs and found six of them may play important roles in IMF development. This work detailed some lincRNAs which may affect of IMF development in pig, and facilitated future research on these lincRNAs and molecular assisted breeding for pig.

Integrated analysis of methylome, transcriptome and miRNAome of three pig breeds.

AIM: Integrated analysis of methylome and transcriptome may help understand the molecular basis of the different breeds with different traits of commercial interest. MATERIALS & METHODS: We obtained the first genome-wide methylome with single-base resolution, miRNAome and transcriptome from three swine breeds. RESULTS: We displayed the landscape of the three omics in the whole-genome level. Integrated outcomes of methylome with genetic selection, miRNAome and transcriptome are also provided. Finally, we identified 11 candidate differentially methylated genes associated with phenotype variance in pigs. CONCLUSION: DNA methylation not only suppresses transcriptome but also miRNAome. The different -omics data have complicated interaction in directly or indirectly and exhibited close relations with the distinct phenotypic traits of growth, disease resistance and energy metabolism.

Assessment of myoblast circular RNA dynamics and its correlation with miRNA during myogenic differentiation.

Myoblast differentiation is a highly complex process that is regulated by proteins as well as by non-coding RNAs. Circular RNAs have been identified as an emerging new class of non-coding RNA in the modulation of skeletal muscle development, whereas their expression profiles and functional regulation in myoblast differentiation remain unknown. In the present study, we performed deep RNA-sequencing of C2C12 myoblasts during cell differentiation and uncovered 37,751 unique circular RNAs derived from 6943 hosting genes. The ensuing qRT-PCR and RNA fluorescence in situ hybridization verification were carried out to confirm the RNA-sequencing results. An unbiased analysis demonstrated dynamic circular RNA expression changes in the process of myoblast differentiation, and the circular RNA abundances were independent from their cognate linear RNAs. Gene ontology analysis showed that many down-regulated circular RNAs were exclusive to cell division and the cell cycle, whereas up-regulated circular RNAs were related to the cell development process. Furthermore, interaction networks of circular RNA-microRNA were constructed. Several microRNAs well-known for myoblast regulation, such as miR-133, miR-24 and miR-23a, were in this network. In summary, this study showed that circular RNA expression dynamics changed during myoblast differentiation. Circular RNAs play a role in regulating the myoblast cell cycle and development by acting as microRNA binding sites to facilitate their regulation of gene expression during myoblast differentiation. These findings open a new avenue for future investigation of this emerging RNA class in skeletal muscle growth and development.