Putative MicroRNA-mRNA Networks Upon Mdfi Overexpression in C2C12 Cell Differentiation and Muscle Fiber Type Transformation

Abstract

Mdfi, an inhibitor of myogenic regulatory factors, is involved in myoblast myogenic development and muscle fiber type transformation. However, the regulatory network of Mdfi regulating myoblasts has not been revealed. In this study, we performed microRNAs (miRNAs)-seq on Mdfi overexpression (Mdfi-OE) and wild-type (WT) C2C12 cells to establish the regulatory networks. Comparative analyses of Mdfi-OE vs. WT identified 66 differentially expressed miRNAs (DEMs). Enrichment analysis of the target genes suggested that DEMs may be involved in myoblast differentiation and muscle fiber type transformation through MAPK, Wnt, PI3K-Akt, mTOR, and calcium signaling pathways. miRNA-mRNA interaction networks were suggested along with ten hub miRNAs and five hub genes. We also identified eight hub miRNAs and eleven hub genes in the networks of muscle fiber type transformation. Hub miRNAs mainly play key regulatory roles in muscle fiber type transformation through the PI3K-Akt, MAPK, cAMP, and calcium signaling pathways. Particularly, the three hub miRNAs (miR-335-3p, miR-494-3p, and miR-709) may be involved in both myogenic differentiation and muscle fiber type transformation. These hub miRNAs and genes might serve as candidate biomarkers for the treatment of muscle- and metabolic-related diseases.

Introduction

Myogenic differentiation is an essential process in myogenesis, regulated by specific transcription factors such as myoblast determination protein (Myod), myogenin (Myog), myogenic factor 5 (Myf5), and myogenic regulatory factor 4 (Mrf4) (Luo et al., 2013). Loss of muscle mass and function due to aging, injury, and muscle diseases decreases quality of life (Argiles et al., 2016). Thus, it is of considerable significance to understand the muscle myogenic differentiation mechanism for treating human muscle-related diseases.

Muscle fiber is the basic unit of skeletal muscle. On the basis of myosin heavy chains (MyHC), muscle fiber can be classified into type I (MyHC I), type IIa (MyHC IIa), type IIx (MyHC IIx), and type IIb (MyHC IIb) (Schiaffino and Reggiani, 2011). Slow-twitch muscles mainly contain type I and type IIa muscle fibers, which have high hemoglobin and mitochondria content, and a strong resistance to fatigue resistance. Fast-twitch muscles mostly contain type IIb muscle fibers, displaying low hemoglobin and mitochondria content, and have a weak resistance to fatigue resistance (Chen X. et al., 2018). Human metabolic diseases, such as insulin resistance, type 2 diabetes, and hypertension, are usually accompanied with the transformation of muscle fiber types (Oberbach et al., 2006; Lu et al., 2017; Aguirre et al., 2020). Therefore, it is essential to elucidate the mechanism of muscle fiber type transformation in order to develop remedies for metabolic-related diseases.

microRNAs (miRNAs) are a class of endogenous, non-coding, short (21–25 nucleotides) RNAs that post-transcriptionally regulate protein-coding genes through interaction with the mRNAs 3′ -untranslated region (Jiao et al., 2018). Many miRNAs play vital roles in skeletal muscle development, including cell proliferation, differentiation, and apoptosis processes (Vienberg et al., 2017). Muscle-specific miRNAs (MyomiRs) are highly expressed in skeletal muscle, including miR-1, miR-133a, miR-133b, and miR-206 (Zhang et al., 2018). They primarily target myogenic transcription factors, such as MyoD (Wu et al., 2019) and serum response factor (SRF) (Kim et al., 2006), to regulate the proliferation and differentiation of skeletal muscle cells. Besides MyomiRs, there are many other miRNAs involved in muscle development regulation. For instance, miR-27a/b act on the myostatin (MSTN) gene to promote myogenic cell proliferation by inhibiting MSTN (McFarlane et al., 2014). miR-29a/b/c and miR-24 target the transcriptional activator SMAD family gene and up-regulate the myogenic differentiation process (Sun et al., 2008; Wei et al., 2013). Myogenic regulatory factor inhibitor (Mdfi) interacts with the MyoD family members. The regulation of its interaction on downstream miRNAs expression is, however, unknown. Therefore, it is of significance to know the relationship between Mdfi and miRNAs for the next step of research on the muscle development regulation mechanism.

Mdfi is a MyoD family inhibitor, also known as I-mfa, which was reported to play a negative role in regulating the MyoD family’s transcription in the myogenic differentiation of NIH-3T3 cells (Chen et al., 1996). Injection of lentivirus vectors containing Mdfi interfering fragments into leg muscles of a 3-week-old mouse showed that interference with Mdfi significantly promoted the hypertrophy of mouse leg muscle fibers (Hou et al., 2018). We previously constructed a Mdfi-overexpressing C2C12 cell line by the CRISPR/Cas9 system and found that overexpression of Mdfi promoted the C2C12 cell myogenic differentiation and muscle fiber transformation from type IIb to type IIa and type I (Huang et al., 2021).

In the present study, we proposed the miRNA-mRNA networks to further explore the function of Mdfi on C2C12 cell myogenic differentiation and muscle fiber type transformation. We selected differentiated C2C12 cells with Mdfi overexpression (Mdfi-OE) and wild-type (WT) for miRNA-seq analysis. We proposed pathways in which differentially expressed miRNAs (DEMs) perform essential functions in myogenic differentiation and muscle fiber type transformation. We also proposed a regulatory network of overexpression of Mdfi on C2C12 cell differentiation and muscle fiber type transformation and identified putative hub miRNAs and genes in them. These findings help to explain how Mdfi regulates myogenic development and muscle fiber type transformation. These hub miRNAs and genes might provide candidate biomarkers for treating muscle- and metabolic-related diseases.

Materials and Methods

Preparation of Cells and Cell Culture

In this study, we selected the mouse skeletal muscle cell line (C2C12), which was preserved in our laboratory (South China Agricultural University, Guangzhou, China) as the experimental cells. The C2C12 cells were cultivated in a basal medium containing Dulbecco’s modified Eagle’s medium (DMEM, Gibco, Grand Island, New York, NY, United States ) with 10% fetal bovine serum. C2C12 cells were seeded in 6-well plates (2 × 10⁵ cells per well). When reached 80–90% confluence, the cells were cultured by myogenic differentiation induction medium (with 2% horse serum) (Hou et al., 2019). We previously constructed an Mdfi stably overexpressing cell line from C2C12 cells by CRISPR/Cas9 (Huang et al., 2021). Our previous studies showed that overexpression of Mdfi promoted C2C12 cell differentiation and positively modulated muscle fiber type transformation (Huang et al., 2021). In the present study, we collected experimental samples (Mdfi-OE) from differentiated C2C12 cells on day 5 and used them for further miRNA-seq analysis.

RNA Extraction and Sequencing

Total RNA was extracted from two groups of six C2C12 cell samples in triplicate with Trizol reagent (Invitrogen, Carlsbad, CA, United States ). The purity of RNA was checked with a NanoDrop ND-2000 spectrophotometer (Thermo Scientific, Waltham, CA, United States ). The integrity of total RNA was assessed with the RNA Nano 6000 assay kit of the Agilent Bioanalyzer 2100 system (Agilent Technologies, Santa Clara, CA, United States ) (Hou et al., 2017). After library preparation, entire libraries from six samples were sequenced according to both 50 bp single-end miRNA methods based on the Illumina Hiseq platform (Caporaso et al., 2012).

Analysis of miRNA Sequencing Data

We used Trimmomatic software to remove the adapter and low-quality sequences from raw data and to clean reads. The sequences from 18 to 32 nt were considered to be high-quality clean reads. After quality control, clean reads were mapped to the mouse reference genome (Mus musculus 10.fa) using BWA software. By using blast, the mapped reads were aligned to the mirbase version 22 (www.mirbase.org) (Kozomara et al., 2019) to identify known miRNAs. The remaining reads were used to detect other small RNA in the Rfam databases (rfam.janelia.org) (Kalvari et al., 2018) by blast software. Differential expression analysis was performed on DESeq2 software, and read counts were the main data. miRNA-seq data was normalized by DESeq2 normalization (Love et al., 2014). The volcano plot and heatmap were generated by (https://software.broadinstitute.org/morpheus) to display an overview of the DEMs. The upset plots were generated by webtool (https://intervene.shinyapps.io/intervene/) (Khan and Mathelier, 2017). The q < 0.05 and |Log2FoldChange| > 1 was set as the cutoff criteria for differentially expressed analysis.

Quantitative Real Time-PCR of the Differentially Expressed miRNAs

To confirm differential expression analysis results, we randomly selected DEMs for quantitative PCR verification. The internal reference gene was U6 in the DEM. The primer sequences designed for the DEMs are listed in Supplementary Table S1. The qRT-PCR reactions were as follows: 40 cycles of 95°C for 20 s, 60°C for 20 s, and 72°C for 20 s. All data were expressed as the mean ± standard error of the mean (S.E.M.). The assumptions of normality of data and homogeneity of variances between the groups were analyzed with SPSS. Significant differences between treatment groups were determined via one-way ANOVA (SPSS 18.0, Chicago, IL, United States ). A value of p < 0.05 was considered to be statistically significant. * stands for p < 0.05, while ** stands for p < 0.01.

Gene Set Enrichment Analysis

Gene Ontology (GO) enrichment analyses of genes were conducted using the Database for Annotation, Visualization, and Integrated Discovery (DAVID, david.ncifcrf.gov) (Huang et al., 2009). Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses of genes were conducted using KOBAS 3.0 (kobas.cbi.pku.edu.cn) (Xie et al., 2011). We used GO enrichment analysis to infer the biological function of DEGs. GO terms including three categories: biological process (BP), cellular component (CC), and molecular function (MF). The KEGG database analysis was used to analyze the enriched pathways for the DEGs. The Benjamini-Hochberg procedure was applied to control the false discovery rate (p < 0.05).

Construction of miRNA–mRNA Interaction Network

We used two strategies to establish miRNA-mRNA regulatory networks: predictive target genes and negative regulatory relationships. Prediction of multiple DEM target genes was performed with the DIANA database (Karagkouni et al., 2018). According to the principle of miRNA’s negative regulatory gene, we overlapped the predicted DEM target genes and DEGs to establish miRNA-mRNA interaction networks. We then utilized Cytoscape 3.6 to analyze and visualize the miRNA-mRNA network.

Results

Overview of miRNA Sequencing Data

The miRNA profiles were generated by sequencing six C2C12 cell samples (three each from Mdfi-OE and WT). The sample’s length, clean reads, Q30 quality scores, GC content, and unique mapped rate were summarized in Supplementary Table S2. After filtration, a total of 63.49 million clean reads were retrieved from the miRNA profiles, of which more than 89% were mapped to the mouse reference genome. The length of the clean reads ranged from 22 to 25 nt, Q30 quality scores of the six samples were above 97%, and the samples’ GC content ranged from 46 to 49%. These data met the standard analysis requirements and could be further analyzed.

Differentially Expressed miRNA Analysis in Mdfi-OE Versus WT

We performed miRNA-Seq to identify microRNAs (miRNAs) differentially expressed between Mdfi-OE and WT groups in triplicate per group. Based on the criteria of q < 0.05 and |Log2FoldChange| > 1, we identified 66 differentially expressed miRNAs (DEMs), among which 27 miRNAs were upregulated, and 39 miRNAs were downregulated by overexpression of Mdfi in C2C12 cells (Supplementary Table S3). The Heatmap and Volcano plot were generated to display an overview of the DEMs (Figures 1A,B). The most abundant upregulated DEM was miR-6356, whose function has not been reported before. The most significant difference downregulated DEM was miR-10a-5p, which belonged to the miR-10 family. The function of miR-10a-5p in skeletal muscle differentiation has not been reported before. However, miRNA-10b, another member of the miRNA-10 family, targeted the nuclear factor of activated T-cells (NFAT) to inhibit myoblast differentiation (Ge et al., 2019). Therefore, the function of miRNA-10 in myogenic differentiation warrants further investigation. The top 10 upregulated and downregulated DEMs with a significant difference (p < 0.05) were listed in Table 1. We also identified four miRNAs (miR-34b-3p, miR-494-3p, miR-182-5p, and miR-351-5p) that have been reported to involve in muscle development and muscle fiber type transformation.

FIGURE 1

Intersection Analysis of DEM, Myogenesis and Myopathies

miRNA is not only involved in the regulation of muscle development (Mok et al., 2017), but also related to the occurrence of muscle diseases (Sayed and Abdellatif, 2011). Therefore, we analyzed the intersection of DEMs, miRNAs involved in myogenesis (Chen et al., 2006), and miRNAs involved in myopathies (Eisenberg et al., 2007). As shown in the Figure 1C, miR-351-5p is muscle-specific miRNA identified in DEMs; In the DEMs, seven miRNAs were associated with muscle disease, such as miR-494-3p, miR-130b, and miR-335-5p. In particular, miR-335 was found upregulated in seven muscle disorders and was the same direction of expression when Mdfi was overexpressed. miR-335 is also the hub miRNA in the miRNA-mRNA interaction network. These hub miRNAs may be the candidate biomarkers for treating muscle-related diseases.

Quantitative Real-Time PCR Validation

To validate the DEMs identified by RNA-seq data, we randomly selected eight DEMs, including two hub miRNAs (miR-335-5p and miR-10a-5p), for analysis by qPCR in Mdfi-OE vs. WT. The qPCR results for DEMs were consistent with the results of the miRNA-seq data, except for miR-190-5p and miR-1258-3p (Figure 2). Thus, the qPCR results had a high consistency rate (75%) with miRNA-seq results, which indicates good reliability of our sequencing data.

FIGURE 2

Prediction and Functional Analysis of Target Genes of DEMs

We herein used DIANA Tools (version 2.0) to predict the target genes of each DEMs (Supplementary Table S4). There were 2074 target genes for 27 upregulated DEMs, consisting of 2454 miRNA-mRNA interaction pairs. We also predicted 2,846 target genes for 39 downregulated DEMs, consisting of 4,083 miRNA-mRNA interaction pairs.

We then performed GO and KEGG enrichment analysis on target genes. GO functional enrichment analyses revealed that 755 and 1189 GO terms were significantly enriched by upregulated and downregulated target genes (p < 0.05), respectively (Supplementary Tables S5, S6). These GO terms were mainly involved in the regulation of transcription, transcription regulator, and nuclear. The top 10 GO terms in each category with a significant difference were shown in Figure 3. There were also 145 KEGG pathways enriched by the target genes of upregulated DEMs (p < 0.05), including the Pathways in cancer, Axon guidance, MAPK signaling pathway, and Regulation of actin cytoskeleton (Supplementary Table S7). There were a total of 184 KEGG pathways enriched by the target genes of downregulated DEMs (p < 0.05), including the pathways in cancer, MAPK signaling pathway, regulation of actin cytoskeleton, and PI3K-Akt signaling pathway (Supplementary Table S7). The top 20 significantly enriched pathways were shown in Figure 3.

FIGURE 3

Integrated Analysis of DEMs and Differentially Expressed Genes in C2C12 Cell With Overexpressed Mdfi

Our previous mRNA-seq analysis identified 434 upregulated DEGs and 1088 downregulated DEGs in Mdfi-OE vs. WT. Mdfi was upregulated 64-fold in the mRNA-seq results (Huang et al., 2021). In the present study, we intersected the prediction target genes of DEMs with these DEGs and formed miRNA-mRNA interaction pairs based on miRNA’s negative regulation genes. A total of 48 genes were overlapped between the 434 upregulated DEGs and 2,846 predicted target genes of downregulated DEMs, and a total of 114 genes were overlapped between the 1,088 downregulated DEGs and 2,074 predicted target genes of upregulated DEMs (Figure 4).

FIGURE 4

We further paired 23 upregulated miRNAs with the above 114 downregulated genes and 25 downregulated miRNAs with the above 48 upregulated genes to form 122 and 71 miRNA-mRNA interactions, respectively. After integrated analysis, we used Cytoscape 3.6 software to display miRNA-mRNA interaction networks. The top 5 significant upregulated miRNAs in the upregulated miRNA-downregulated mRNA interaction network were miR-335-3p (degree = 120), miR-494-3p (degree = 45), miR-6938-3p (degree = 41), miR-709 (degree = 41), and miR-34b-5p (degree = 31) (Supplementary Figure S1). Many downregulated genes were also regulated by multiple miRNAs. Runx1t1 (degree = 8), Onecut2 (degree = 8), and Stc1 (degree = 6) were targeted by three hub miRNAs (Table 2). We further selected the first neighbors of these five hub miRNAs to construct the hub upregulated miRNA-downregulated mRNA interaction network (Figure 5).

FIGURE 5

In the downregulated miRNA-upregulated mRNA interaction network, the top 5 significant downregulated miRNAs were miR-330-5p (degree = 24), miR-9-5p (degree = 23), miR-322-5p (degree = 22), miR-96-5p (degree = 19), and miR-182-5p (degree = 18) (Supplementary Figure S2). The upregulated genes, Ttn (degree = 13) and Slc5a3 (degree = 9), were targeted by all five hub miRNAs (Table 2). Similarly, we constructed the hub downregulated miRNA-upregulated mRNA interaction network (Figure 6).

FIGURE 6

Integrated Analysis of DEMs and DEGs of Overexpressed Mdfi on Muscle Fiber Type Transformation

Overexpression of Mdfi positively modulates fast-to-slow-twitch muscle fibers transformation (Huang et al., 2021). We herein used miRNA-seq analysis to identify the miRNAs involved in muscle fiber type transformation. We selected the 101 DEGs related to the transformation of muscle fiber types in our previous study and matched these DEGs with the DEMs in the miRNA-seq.

Furthermore, we proposed miRNA-mRNA interaction networks related to muscle fiber type transformation (Figure 7) and identified 11 miRNA-mRNA interaction pairs. These interaction pairs included four upregulated miRNAs, four downregulated miRNAs, four upregulated genes, and seven downregulated genes. In the upregulated miRNAs, miR-335-3p, miR-494-3p, and miR-709 were also hub miRNAs in the miRNA-mRNA interaction network of Mdfi overexpressed differentiated C2C12 cells. This result indicates that these miRNAs may be involved not only in the process of myogenic differentiation but also in muscle fiber type transformation.

FIGURE 7

Discussion

We previously constructed an Mdfi-overexpressing C2C12 cell line, performed mRNA-seq on Mdfi-OE and WT C2C12 cells. Overexpression of Mdfi promoted C2C12 cell differentiation and positively modulated muscle fiber type transformation (Huang et al., 2021). Mdfi overexpression may not be the only possibility to regulate myoblasts differentiation. Many factors may regulate myoblasts differentiation, but the regulatory pathways are different. For example, MYOC regulates TGFb by inhibiting CAV1 to promote myoblasts differentiation. Lmod3 promotes differentiation of C2C12 cells by activating the AKT and ERK pathways. However, the regulatory network in Mdfi-overexpressing C2C12 cells has not been revealed. In this study, we performed miRNA-seq on Mdfi-OE and WT C2C12 cells. We drafted miRNA-mRNA interaction networks involved in myogenic differentiation and muscle fiber type transformation regulated by Mdfi, and identified hub miRNAs and genes in networks. These results may help in understanding the mechanisms of Mdfi in myogenic development and muscle fiber type transformation.

Differentially Expressed miRNA Analysis in Mdfi-OE Versus WT

In the comparative analysis of Mdfi-OE vs. WT, we identified 66 DEMs. Among these DEMs, four miRNAs (miR-34b-3p, miR-494-3p, miR-182-5p, and miR-351-5p) have been reported to involve in muscle development and muscle fiber type transformation. miR-34b-3p promoted the myogenic differentiation of C2C12 cells and produced more multicore myotubes (Tang et al., 2017). miR-351-5p targeted lactamase beta (Lactb) to inhibit C2C12 cell differentiation and led to muscle fiber type transformation (Du et al., 2019). The knockout of miR-182-5p promoted muscle fibers’ transformation from type IIa to type I (Zhang et al., 2016). The overexpression of miR-494-3p inhibited myotube formation and the expression of MYH2 during human skeletal myogenesis (Yamamoyo et al., 2012). The expression trend of three miRNAs (miR-34b-3p, miR-351-5p, and miR-182-5p) was consistent with our previous results that overexpression of Mdfi promotes the differentiation of C2C12 cells and muscle fiber type transformation. However, based on the function of miRNA-494 previously reported, the expression trend of miR-494-3p in miRNA-seq was inconsistent with our previous results. We hypothesized that a novel interaction might exist between Mdfi and miR-494-3p.

Functional Analysis of Target Genes of DEMs in Mdfi-OE Versus WT

The enriched GO terms of upregulated and downregulated DEMs target genes were mainly involved in the regulation of transcription, nuclear, and metal ion binding and transport. Mdfi plays essential roles in inhibiting protein transport into the nucleus and ultimately inhibiting its transcriptional activity (Chen et al., 1996). These results suggested that DEM may regulate gene entry into the nucleus and its transcriptional activity in C2C12 cells with overexpressed Mdfi.

The KEGG pathways enriched by upregulated and downregulated DEMs target genes, especially MAPK signaling pathway (Troy et al., 2012), Wnt signaling pathway (Girardi and Le Grand, 2018), PI3K-Akt signaling pathway (Rommel et al., 2001), mTOR signaling pathway (Go et al., 2020), and calcium signaling pathway (Swoap et al., 2000), and are all known to be involved in myogenic differentiation and muscle fiber type transformation. In particular, PI3K/AKT signaling pathway activation has been reported to induce myomiR maturation in C2C12 cells (Briata et al., 2012). mTOR is a crucial pathway signaling pathway of skeletal muscle hypertrophy associated with increased protein synthesis (Go et al., 2020). mTOR regulates mitochondrial function and participate in glycolytic flux and mitochondrial respiration (Tan et al., 2015). The increase of intracellular cAMP level increases the expression of CaMKIV and PGC-1 in myoblast cells, and then regulate the formation of muscle fibers (Handschin et al., 2003). There is feedback regulation between cAMP and intracellular calcium homeostasis, which plays an important role in the regulation of muscle fiber transformation(Walke et al., 1994). miRNA-21 has been reported to inhibit PI3K/Akt/mTOR signaling by targeting TGFBI during skeletal muscle development (Briata et al., 2012). Muscle development is associated with calcium homeostasis, and many miRNAs in calcium signaling appeared to be suppressed in aged skeletal muscles (Yan et al., 2018). These results suggest that DEM may regulate myoblast differentiation and transform muscle fiber types through these signaling pathways.

Identification of Hub miRNAs and Genes in Regulatory Networks Associated With Mdfi Overexpression on C2C12 Cells Differentiation

To investigate how miRNA regulates the gene expression of overexpressed Mdfi, we drafted miRNA-mRNA interaction networks associated with Mdfi overexpression on C2C12 cell differentiation. Among the networks, high degreed-nodes in the regulatory network are also known as hub miRNAs or genes.

In the hub upregulated miRNA-downregulated mRNA interaction network, we identified five hub miRNAs (miR-335-3p, miR-494-3p, miR-6938-3p, miR-709, and miR-34b-5p) and three hub genes [RUNX1 partner transcriptional co-repressor 1 (Runx1t1), one cut homeobox 2 (Onecut2), and stanniocalcin 1 (Stc1)]. These hub miRNAs paired with hub genes to form nine interaction pairs (miR-335-3p-Runx1t1/Onecut2/Stc1, miR-709-RUNX1T1unx1t1/Onecut2/STC1, miR-494-3p/ONECUT2necut2, miR-6938-3p/RUNX1T1unx1t1, and miR-34b-5p/Stc1). MiR-335-3p was induced during the differentiation of C2C12 cells and highly expressed during muscle regeneration (Greco et al., 2009). MiR-709 has been documented to be regulated adipocyte differentiation by targeting glycogen synthase kinase3b (Gsk3b) (Chen et al., 2014), but its function on muscle development has not been reported. MiR-494-3p has been reported to inhibit the differentiation of myoblast and lead to muscle fiber type transformation (Iwasaki et al., 2015). Still, its expression trend was inconsistent with the results of Mdfi promotes C2C12 cell differentiation. Therefore, we hypothesized that a novel interaction might exist between Mdfi and miR-494-3p, thereby regulating myoblast differentiation. To our knowledge, there are no published reports on miR-6938-3p. MiR-34b-5p was previously reported to promote the differentiation of C2C12 cells, and its upregulation is required for myotube formation (Tang et al., 2017). In our study, miR-34b-5p was upregulated, indicating that Mdfi may promote the differentiation of C2C12 cells by upregulating the expression of miR-34b-5p. Runx1t1 inhibited brown adipogenesis and inhibited white adipogenesis in 3T3 cells (Sun et al., 2011). Stc1 was involved in regulating calcium homeostasis in mammals and expressed from myoblast differentiation to myotube fusion. (Jiang et al., 2000). Onecut2 has been discovered as a transcriptional factor and regulated the hepatocyte differentiation (Francius and Clotman, 2010), but its function on muscle development has not been reported. In addition, in our enrichment analysis, these three hub genes were not enriched in any pathways involved in myogenic differentiation. All the above hub miRNAs and genes were significantly changed in our RNA-Seq results. They play essential roles in the process of Mdfi regulating C2C12 cell differentiation.

By contrast, we identified five hub miRNAs (miR-330-5p, miR-9-5p, miR-322-5p, miR-96-5p, and miR-182-5p) and two hub genes [titin (Ttn) and solute carrier family 5 member 3 (Slc5a3)] in the hub downregulated miRNA-upregulated mRNA interaction network. They formed ten interactions pairs (miR-330-5p/miR-9-5p/miR-322-5p/miR-96-5p/miR-182-5p-Ttn/Slc5a3). In the present study, miR-330-5p was proposed to bind STAT3, inhibiting muscle wasting (Mubaid et al., 2019). The overexpression of miR-9-5p significantly inhibited myoblasts differentiation, whereas the inhibition of miR-9-5p had the opposite effect (Yin et al., 2020). In our study, miR-9-5p was downregulated, indicating that Mdfi may promote the differentiation of C2C12 cells by downregulating miR-9-5p. The inhibition of miR-322-5p attenuated Dexamethasone (Dex)-induced muscle atrophy in C2C12 myotubes (Geng et al., 2020). MiR-96-5p has been reported to regulate proliferation, migration, and apoptosis of vascular smooth muscle cells (Tian et al., 2020), but its function on muscle cell differentiation remains unclear. MiR-182-5p was highly expressed in the mouse’s fast-twitch muscles, and knockout of miR-182-5p promoted MyHC I expression and reduced MyHC IIa expression (Zhang et al., 2016). TTN, an abundant filamentous protein, is activated by calcium. TTN promoted muscle development by increasing the expression of MYH6, myosin light chain 1 (MYL1) (Kelly et al., 2020). SLC5A3, also known as SMIT1, is expressed in cardiomyocytes, and its primary function is to detect increased extracellular glucose concentrations (Van Steenbergen et al., 2017).

All the above hub miRNA genes were significantly changed in the miRNA-Seq results. Therefore, we hypothesized that overexpression of Mdfi would promote C2C12 cell differentiation through these 19 interaction pairs, and the role of ten hub miRNAs and five genes contained in interaction pairs needs to be further explored. In particular, the hub miRNAs mentioned above, we identified two miRNAs (miR-182-5p and miR-494-3p), suggesting that overexpression of Mdfi may induce muscle fiber type transformation. Their functions need to be further explored.

Integrated Analysis of Muscle Fiber Type Transformation in Response to Mdfi Overexpression

In the miRNA-mRNA interaction networks, we identified four hub miRNA (miR-335-3p, miR-335-3p, miR-494-3p, and miR-709) and seven hub genes [collagen type I alpha 2 chain (Col1a2), 5-hydroxytryptamine receptor 7 (Htr7), lysophosphatidic acid receptor 4 (Lpar4), prostaglandin E receptor 2 (Ptger2), fms related receptor tyrosine kinase 1 (Flt1), phosphodiesterase 1B (Pde1b), and ring finger protein 152 (Rnf152)]. These hub miRNAs paired with hub genes to form seven interaction pairs (miR-335-3p-Col1a2/Htr7/Lpar4/Ptger2, miR-335-5p-Flt1, miR-709-Pde1b, miR-494-3p-Rnf152). MiR-335-3p was previously reported to involve in the differentiation of myoblasts (Greco et al., 2009), but its function on muscle fiber type transformation has remained unclear. MiR-335-5p was shown to be involved in cardiac differentiation by regulating WNT and TGFb signaling pathways (Kay et al., 2019). To our knowledge, there are no published reports on the function of miR-709 in muscle fiber type transformation. As we described above, miR-494-3p has been reported to inhibit the expression of Myh2, the slow-twitch muscle fiber marker gene during human skeletal myogenesis (Iwasaki et al., 2015). However, Myh2 was upregulated in our previous results (Huang et al., 2021) and miR-494-3p was upregulated in this study, which was inconsistent with the literature reports. The expression trend of miR-494-3p was also inconsistent with our results of Mdfi positively modulates fast-to-slow-twitch muscle fibers transformation. Therefore, we hypothesized that it might be attributed to different species, or that a novel interaction may exist between Mdfi and miR-494-3p, thereby regulating the transformation of muscle fiber types. Col1a2 is a subtype of collagen 1. Collagen type I was proposed to markedly inhibit C2C12 cell differentiation and muscle regeneration, and increase collagen production (Alexakis et al., 2007). Htr7 promoted axon growth and synapse formation of the nerve (Olusakin et al., 2020), but its functions on muscle have not been reported before. Lpar4, also known as GPR23, is a G-protein coupled receptor and has been reported to increase intracellular concentrations of cAMP and calcium (Lee et al., 2007). Flt1 was downregulated in the myogenic differentiation and expressed in regenerating muscle fibers (Germani et al., 2003). Similarly, Flt1 was downregulated and identified as a hub gene in the miRNA-mRNA network of this study, suggesting possible regulation via miR-335-5p and involvement in muscle fiber type transformation. Pde1b is distributed in the human myocardium and activated by calcium/calmodulin (Wennogle et al., 2017). Rnf152 is an essential negative regulator of the mTOR pathway by targeting RagA (Deng et al., 2015). Rnf152 was downregulated and targeted by miR-494-3p in this network. Thus, Mdfi may downregulate Rnf152 via upregulating miR-494-3p, thereby regulating muscle fiber type transformation through the mTOR signaling pathway.

Similarly, we identified four hub miRNAs (miR-130-3p, miR-6540-5p, miR-10a-5p, and miR-330-3p) and four hub genes [phosphatidylinositol-4,5-bisphosphate 3-kinase catalytic subunit beta (Pik3cb), ryanodine receptor 3 (Ryr3), calmodulin-dependent protein kinase II beta (Camk2b), and Prospero homeobox 1 (Prox1)] in the downregulated miRNA-upregulated mRNA interaction network. They formed four interactions pairs (miR-130-3p-Pik3cb, miR-6540-5p-Ryr3, miR-10a-5p-Camk2b, miR-330-3p-Prox1). MiR-130-3p is closely related to the occurrence and development of cancer (Wang et al., 2020). There are no published reports on miR-6540-5p. The function of miR-10a-5p in muscle fiber type transformation has not been reported, but miRNA-10 family member, miRNA-10b, targeted NFAT to inhibit myoblast differentiation (Ge et al., 2019). Therefore, miRNA- 10a-5p may also be related to muscle development. Overexpression of miR-330-3p induced cardiac hypertrophy in mice (Chen Y. et al., 2018). In C2C12 cells, the knockout of Pik3cb delayed differentiation, whereas overexpression of Pik3cb promoted differentiation (Matheny et al., 2015). In skeletal muscle, Ryr3 is a sarcoplasmic reticulum calcium release channel (Protasi et al., 2000), which regulated the release and balance of calcium. Ryr3 was increased in our mRNA-seq results and enriched in the calcium signaling pathway, suggesting that Mdfi may regulate muscle fiber type transformation by regulating calcium concentration. Camk2b, a calcium-dependent kinase, promoted mitochondrial biogenesis, and participated in the formation of slow-twitch muscle fibers (Al-Shanti and Stewart, 2009). The knockout of Prox1 in the skeletal muscle caused slow to fast skeletal muscle fiber transformation (Petchey et al., 2014). In this study, Prox1 was upregulated; therefore, we thus speculated that Prox1 upregulation was involved in Mdfi regulated fast-to-slow- twitch muscle fibers transformation. These hub genes also were enriched in the PI3K-Akt signaling pathway, the cAMP signaling pathway, the calcium signaling pathway, and the mTOR signaling pathway. This result suggests that these hub miRNAs mainly play key regulatory roles in muscle fiber type transformation through these signaling pathways.

All the eight hub miRNAs and eleven genes in the interaction pairs mentioned above play essential roles in the process of Mdfi regulating the transformation of muscle fiber types. Their role needs to be further explored. In addition, three hub miRNAs (miR-335-3p, miR-494-3p, and miR-709) identified in the network of C2C12 cell differentiation were also identified as the hubs genes in the network of muscle fiber type transformation. Their functions need to be further explored.

Conclusion

In conclusion, based on miRNAs’ expression profiles, we performed differential analysis, enrichment analysis, and proposed miRNA-mRNA interaction networks of Mdfi overexpression on C2C12 cell differentiation and muscle fiber type transformation (Figure 8). We have putatively identified the pathways in which DEM performs its functions in myogenic differentiation and muscle fiber type transformation. Furthermore, we putatively identified hub miRNAs and genes in the miRNA-mRNA networks. In particular, three hub miRNAs (miR-335-3p, miR-494-3p, and miR-709) were involved in regulating both myogenic differentiation and muscle fiber type transformation. Our findings may provide insights into the regulatory role of Mdfi in myogenic development and muscle fiber type transformation. The hub miRNAs and genes might provide candidate biomarkers for treating muscle- and metabolic-related diseases. The results warrant our interest in verifying the relationship between Mdfi and hub miRNAs/genes to understand how Mdfi regulates C2C12 cell differentiation and muscle fiber type transformation.

FIGURE 8

References

• 1

  AguirreF.AbrigoJ.GonzalezF.GonzalezA.SimonF.Cabello-VerrugioC. (2020). Protective Effect of Angiotensin 1-7 on Sarcopenia Induced by Chronic Liver Disease in Mice. Int. J. Mol. Sci.21, 3891. 10.3390/ijms21113891

  • CrossRef
  • Google Scholar

• 2

  Al-ShantiN.StewartC. E. (2009). Ca2+/calmodulin-dependent Transcriptional Pathways: Potential Mediators of Skeletal Muscle Growth and Development. Biol. Rev. Camb Philos. Soc.84, 637–652. 10.1111/j.1469-185X.2009.00090.x

  • CrossRef
  • Google Scholar

• 3

  AlexakisC.PartridgeT.Bou-GhariosG. (2007). Implication of the Satellite Cell in Dystrophic Muscle Fibrosis: a Self-Perpetuating Mechanism of Collagen Overproduction. Am. J. Physiol. Cell Physiol.293, C661–C669. 10.1152/ajpcell.00061.2007

  • CrossRef
  • Google Scholar

• 4

  ArgilésJ. M.CamposN.Lopez-PedrosaJ. M.RuedaR.Rodriguez-MañasL. (2016). Skeletal Muscle Regulates Metabolism via Interorgan Crosstalk: Roles in Health and Disease. J. Am. Med. Directors Assoc.17, 789–796. 10.1016/j.jamda.2016.04.019

  • CrossRef
  • Google Scholar

• 5

  BriataP.LinW.-J.GiovarelliM.PaseroM.ChouC.-F.TrabucchiM.et al (2012). PI3K/AKT Signaling Determines a Dynamic Switch between Distinct KSRP Functions Favoring Skeletal Myogenesis. Cell Death Differ.19, 478–487. 10.1038/cdd.2011.117

  • CrossRef
  • Google Scholar

• 6

  CaporasoJ. G.LauberC. L.WaltersW. A.Berg-LyonsD.HuntleyJ.FiererN.et al (2012). Ultra-high-throughput Microbial Community Analysis on the Illumina HiSeq and MiSeq Platforms. ISME J.6, 1621–1624. 10.1038/ismej.2012.8

  • CrossRef
  • Google Scholar

• 7

  ChenC.-M. A.KrautN.GroudineM.WeintraubH. (1996). I-mf, a Novel Myogenic Repressor, Interacts with Members of the MyoD Family. Cell86, 731–741. 10.1016/s0092-8674(00)80148-8

  • CrossRef
  • Google Scholar

• 8

  ChenH.MoD.LiM.ZhangY.ChenL.ZhangX.et al (2014). miR-709 Inhibits 3T3-L1 Cell Differentiation by Targeting GSK3β of Wnt/β-Catenin Signaling. Cell Signal.26, 2583–2589. 10.1016/j.cellsig.2014.07.017

  • CrossRef
  • Google Scholar

• 9

  ChenJ.-F.MandelE. M.ThomsonJ. M.WuQ.CallisT. E.HammondS. M.et al (2006). The Role of microRNA-1 and microRNA-133 in Skeletal Muscle Proliferation and Differentiation. Nat. Genet.38, 228–233. 10.1038/ng1725

  • CrossRef
  • Google Scholar

• 10

  ChenX.GuoY.JiaG.LiuG.ZhaoH.HuangZ. (2018a). Arginine Promotes Skeletal Muscle Fiber Type Transformation from Fast-Twitch to Slow-Twitch via Sirt1/AMPK Pathway. J. Nutr. Biochem.61, 155–162. 10.1016/j.jnutbio.2018.08.007

  • CrossRef
  • Google Scholar

• 11

  ChenY.LiuX.ChenL.ChenW.ZhangY.ChenJ.et al (2018b). The Long Noncoding RNA XIST Protects Cardiomyocyte Hypertrophy by Targeting miR-330-3p. Biochem. Biophys. Res. Commun.505, 807–815. 10.1016/j.bbrc.2018.09.135

  • CrossRef
  • Google Scholar

• 12

  DengL.JiangC.ChenL.JinJ.WeiJ.ZhaoL.et al (2015). The Ubiquitination of Raga GTPase by RNF152 Negatively Regulates mTORC1 Activation. Mol. Cel. 58, 804–818. 10.1016/j.molcel.2015.03.033

  • CrossRef
  • Google Scholar

• 13

  DuJ.ZhangP.ZhaoX.HeJ.XuY.ZouQ.et al (2019). MicroRNA‐351‐5p Mediates Skeletal Myogenesis by Directly Targeting Lactamase‐β and Is Regulated by Lnc‐mg. FASEB j.33, 1911–1926. 10.1096/fj.201701394RRR

  • CrossRef
  • Google Scholar

• 14

  EisenbergI.EranA.NishinoI.MoggioM.LampertiC.AmatoA. A.et al (2007). Distinctive Patterns of microRNA Expression in Primary Muscular Disorders. Proc. Natl. Acad. Sci.104, 17016–17021. 10.1073/pnas.0708115104

  • CrossRef
  • Google Scholar

• 15

  FranciusC.ClotmanF. (2010). Dynamic Expression of the Onecut Transcription Factors HNF-6, OC-2 and OC-3 during Spinal Motor Neuron Development. Neuroscience165, 116–129. 10.1016/j.neuroscience.2009.09.076

  • CrossRef
  • Google Scholar

• 16

  GeG.YangD.TanY.ChenY.JiangD.JiangA.et al (2019). miR-10b-5p Regulates C2C12 Myoblasts Proliferation and Differentiation. Biosci. Biotechnol. Biochem.83, 291–299. 10.1080/09168451.2018.1533805

  • CrossRef
  • Google Scholar

• 17

  GengH.SongQ.ChengY.LiH.YangR.LiuS.et al (2020). MicroRNA 322 Aggravates Dexamethasone-Induced Muscle Atrophy by Targeting IGF1R and INSR. Int. J. Mol. Sci.21, 1111. 10.3390/ijms21031111

  • CrossRef
  • Google Scholar

• 18

  GermaniA.Di CarloA.MangoniA.StrainoS.GiacintiC.TurriniP.et al (2003). Vascular Endothelial Growth Factor Modulates Skeletal Myoblast Function. Am. J. Pathol.163, 1417–1428. 10.1016/S0002-9440(10)63499-2

  • CrossRef
  • Google Scholar

• 19

  GirardiF.Le GrandF. (2018). Wnt Signaling in Skeletal Muscle Development and Regeneration. Prog. Mol. Biol. Transl. Sci.153, 157–179. 10.1016/bs.pmbts.2017.11.026

  • CrossRef
  • Google Scholar

• 20

  GoG.-Y.JoA.SeoD.-W.KimW.-Y.KimY. K.SoE.-Y.et al (2020). Ginsenoside Rb1 and Rb2 Upregulate Akt/mTOR Signaling-Mediated Muscular Hypertrophy and Myoblast Differentiation. J. Ginseng Res.44, 435–441. 10.1016/j.jgr.2019.01.007

  • CrossRef
  • Google Scholar

• 21

  GrecoS.De SimoneM.ColussiC.ZaccagniniG.FasanaroP.PescatoriM.et al (2009). Common micro‐RNA Signature in Skeletal Muscle Damage and Regeneration Induced by Duchenne Muscular Dystrophy and Acute Ischemia. FASEB j.23, 3335–3346. 10.1096/fj.08-128579

  • CrossRef
  • Google Scholar

• 22

  HandschinC.RheeJ.LinJ.TarrP. T.SpiegelmanB. M. (2003). An Autoregulatory Loop Controls Peroxisome Proliferator-Activated Receptor Coactivator 1 Expression in Muscle. Proc. Natl. Acad. Sci.100, 7111–7116. 10.1073/pnas.1232352100

  • CrossRef
  • Google Scholar

• 23

  HouL.ShiJ.CaoL.XuG.HuC.WangC. (2017). Pig Has No Uncoupling Protein 1. Biochem. Biophys. Res. Commun.487, 795–800. 10.1016/j.bbrc.2017.04.118

  • CrossRef
  • Google Scholar

• 24

  HouL.XuJ.JiaoY.LiH.PanZ.DuanJ.et al (2018). MiR-27b Promotes Muscle Development by Inhibiting MDFI Expression. Cell Physiol. Biochem.46, 2271–2283. 10.1159/000489595

  • CrossRef
  • Google Scholar

• 25

  HouL.ZhuL.LiH.JiangF.CaoL.HuC. Y.et al (2019). MiR-501-3p Forms a Feedback Loop with FOS, MDFI, and MyoD to Regulate C2C12 Myogenesis. Cells8, 573. 10.3390/cells8060573

  • CrossRef
  • Google Scholar

• 26

  HuangB.JiaoY.ZhuY.NingZ.YeZ.LiQ. X.et al (2021). Mdfi Promotes C2C12 Cell Differentiation and Positively Modulates Fast-To-Slow-Twitch Muscle Fiber Transformation. Front. Cel Dev. Biol.9, 605875. 10.3389/fcell.2021.605875

  • CrossRef
  • Google Scholar

• 27

  HuangD. W.ShermanB. T.LempickiR. A. (2009). Systematic and Integrative Analysis of Large Gene Lists Using DAVID Bioinformatics Resources. Nat. Protoc.4, 44–57. 10.1038/nprot.2008.211

  • CrossRef
  • Google Scholar

• 28

  IwasakiH.ImamuraT.MorinoK.ShimosatoT.TawaM.UgiS.et al (2015). MicroRNA-494 Plays a Role in Fiber Type-specific Skeletal Myogenesis in Human Induced Pluripotent Stem Cells. Biochem. Biophys. Res. Commun.468, 208–213. 10.1016/j.bbrc.2015.10.128

  • CrossRef
  • Google Scholar

• 29

  JiangW.ChangA.SatohM.FuruichiY.TamP.ReddelR. (2000). The Distribution of Stanniocalcin 1 Protein in Fetal Mouse Tissues Suggests a Role in Bone and Muscle Development. J. Endocrinol.165, 457–466. 10.1677/joe.0.1650457

  • CrossRef
  • Google Scholar

• 30

  JiaoY.HuangB.ChenY.HongG.XuJ.HuC.et al (2018). Integrated Analyses Reveal Overexpressed Notch1 Promoting Porcine Satellite Cells' Proliferation through Regulating the Cell Cycle. Int. J. Mol. Sci.19, 271. 10.3390/ijms19010271

  • CrossRef
  • Google Scholar

• 31

  KalvariI.ArgasinskaJ.Quinones-OlveraN.NawrockiE. P.RivasE.EddyS. R.et al (2018). Rfam 13.0: Shifting to a Genome-Centric Resource for Non-coding RNA Families. Nucleic Acids Res.46, D335–D342. 10.1093/nar/gkx1038

  • CrossRef
  • Google Scholar

• 32

  KaragkouniD.ParaskevopoulouM. D.ChatzopoulosS.VlachosI. S.TastsoglouS.KanellosI.et al (2018). DIANA-TarBase V8: a Decade-Long Collection of Experimentally Supported miRNA-Gene Interactions. Nucleic Acids Res.46, D239–D245. 10.1093/nar/gkx1141

  • CrossRef
  • Google Scholar

• 33

  KayM.SoltaniB. M.AghdaeiF. H.AnsariH.BaharvandH. (2019). Hsa-miR-335 Regulates Cardiac Mesoderm and Progenitor Cell Differentiation. Stem Cel Res. Ther.10, 191. 10.1186/s13287-019-1249-2

  • CrossRef
  • Google Scholar

• 34

  KellyC. M.ManukianS.KimE.GageM. J. (2020). Differences in Stability and Calcium Sensitivity of the Ig Domains in Titin's N2A Region. Protein Sci.29, 1160–1171. 10.1002/pro.3848

  • CrossRef
  • Google Scholar

• 35

  KhanA.MathelierA. (2017). Intervene: a Tool for Intersection and Visualization of Multiple Gene or Genomic Region Sets. BMC Bioinformatics18, 287. 10.1186/s12859-017-1708-7

  • CrossRef
  • Google Scholar

• 36

  KimH. K.LeeY. S.SivaprasadU.MalhotraA.DuttaA. (2006). Muscle-specific microRNA miR-206 Promotes Muscle Differentiation. J. Cel Biol.174, 677–687. 10.1083/jcb.200603008

  • CrossRef
  • Google Scholar

• 37

  KozomaraA.BirgaoanuM.Griffiths-JonesS. (2019). miRBase: from microRNA Sequences to Function. Nucleic Acids Res.47, D155–D162. 10.1093/nar/gky1141

  • CrossRef
  • Google Scholar

• 38

  LeeC.-W.RiveraR.DubinA. E.ChunJ. (2007). LPA4/GPR23 Is a Lysophosphatidic Acid (LPA) Receptor Utilizing Gs-, Gq/Gi-Mediated Calcium Signaling and G12/13-Mediated Rho Activation. J. Biol. Chem.282, 4310–4317. 10.1074/jbc.M610826200

  • CrossRef
  • Google Scholar

• 39

  LoveM. I.HuberW.AndersS. (2014). Moderated Estimation of Fold Change and Dispersion for RNA-Seq Data with DESeq2. Genome Biol.15, 550. 10.1186/s13059-014-0550-8

  • CrossRef
  • Google Scholar

• 40

  LuJ.-J.WangQ.XieL. H.ZhangQ.SunS. H. (2017). Tumor Necrosis Factor-like Weak Inducer of Apoptosis Regulates Quadriceps Muscle Atrophy and Fiber-type Alteration in a Rat Model of Chronic Obstructive Pulmonary Disease. Tob. Induced Dis.15, 43. 10.1186/s12971-017-0148-5

  • CrossRef
  • Google Scholar

• 41

  LuoW.NieQ.ZhangX. (2013). MicroRNAs Involved in Skeletal Muscle Differentiation. J. Genet. Genomics40, 107–116. 10.1016/j.jgg.2013.02.002

  • CrossRef
  • Google Scholar

• 42

  MathenyR. W.Jr.Riddle-KottkeM. A.LeandryL. A.LynchC. M.AbdallaM. N.GeddisA. V.et al (2015). Role of Phosphoinositide 3-OH Kinase P110β in Skeletal Myogenesis. Mol. Cel. Biol.35, 1182–1196. 10.1128/MCB.00550-14

  • CrossRef
  • Google Scholar

• 43

  McFarlaneC.VajjalaA.ArigelaH.LokireddyS.GeX.BonalaS.et al (2014). Negative Auto-Regulation of Myostatin Expression Is Mediated by Smad3 and microRNA-27. PLoS One9, e87687. 10.1371/journal.pone.0087687

  • CrossRef
  • Google Scholar

• 44

  MokG. F.Lozano-VelascoE.MünsterbergA. (2017). microRNAs in Skeletal Muscle Development. Semin. Cel Develop. Biol.72, 67–76. 10.1016/j.semcdb.2017.10.032

  • CrossRef
  • Google Scholar

• 45

  MubaidS.MaJ. F.OmerA.AshourK.LianX. J.SanchezB. J.et al (2019). HuR Counteracts miR-330 to Promote STAT3 Translation during Inflammation-Induced Muscle Wasting. Proc. Natl. Acad. Sci. USA116, 17261–17270. 10.1073/pnas.1905172116

  • CrossRef
  • Google Scholar

• 46

  OberbachA.BossenzY.LehmannS.NiebauerJ.AdamsV.PaschkeR.et al (2006). Altered Fiber Distribution and Fiber-specific Glycolytic and Oxidative Enzyme Activity in Skeletal Muscle of Patients with Type 2 Diabetes. Diabetes Care29, 895–900. 10.2337/diacare.29.04.06.dc05-1854

  • CrossRef
  • Google Scholar

• 47

  OlusakinJ.MoutkineI.DumasS.PonimaskinE.PaizanisE.Soiza-ReillyM.et al (2020). Implication of 5-HT7 Receptor in Prefrontal Circuit Assembly and Detrimental Emotional Effects of SSRIs during Development. Neuropsychopharmacol45, 2267–2277. 10.1038/s41386-020-0775-z

  • CrossRef
  • Google Scholar

• 48

  PetcheyL. K.RisebroC. A.VieiraJ. M.RobertsT.BrysonJ. B.GreensmithL.et al (2014). Loss ofProx1in Striated Muscle Causes Slow to Fast Skeletal Muscle Fiber Conversion and Dilated Cardiomyopathy. Proc. Natl. Acad. Sci. USA111, 9515–9520. 10.1073/pnas.1406191111

  • CrossRef
  • Google Scholar

• 49

  ProtasiF.TakekuraH.WangY.ChenS. R. W.MeissnerG.AllenP. D.et al (2000). RYR1 and RYR3 Have Different Roles in the Assembly of Calcium Release Units of Skeletal Muscle. Biophys. J.79, 2494–2508. 10.1016/S0006-3495(00)76491-5

  • CrossRef
  • Google Scholar

• 50

  RommelC.BodineS. C.ClarkeB. A.RossmanR.NunezL.StittT. N.et al (2001). Mediation of IGF-1-Induced Skeletal Myotube Hypertrophy by PI(3)K/Akt/mTOR and PI(3)K/Akt/GSK3 Pathways. Nat. Cel Biol.3, 1009–1013. 10.1038/ncb1101-1009

  • CrossRef
  • Google Scholar

• 51

  SayedD.AbdellatifM. (2011). MicroRNAs in Development and Disease. Physiol. Rev.91, 827–887. 10.1152/physrev.00006.2010

  • CrossRef
  • Google Scholar

• 52

  SchiaffinoS.ReggianiC. (2011). Fiber Types in Mammalian Skeletal Muscles. Physiol. Rev.91, 1447–1531. 10.1152/physrev.00031.2010

  • CrossRef
  • Google Scholar

• 53

  SunL.XieH.MoriM. A.AlexanderR.YuanB.HattangadiS. M.et al (2011). Mir193b-365 Is Essential for Brown Fat Differentiation. Nat. Cel Biol.13, 958–965. 10.1038/ncb2286

  • CrossRef
  • Google Scholar

• 54

  SunQ.ZhangY.YangG.ChenX.ZhangY.CaoG.et al (2008). Transforming Growth Factor- -regulated miR-24 Promotes Skeletal Muscle Differentiation. Nucleic Acids Res.36, 2690–2699. 10.1093/nar/gkn032

  • CrossRef
  • Google Scholar

• 55

  SwoapS. J.HunterR. B.StevensonE. J.FeltonH. M.KansagraN. V.LangJ. M.et al (2000). The Calcineurin-NFAT Pathway and Muscle Fiber-type Gene Expression. Am. J. Physiol. Cell Physiol.279, C915–C924. 10.1152/ajpcell.2000.279.4.C915

  • CrossRef
  • Google Scholar

• 56

  TanB.ShaoF.YinY.XiaoH.WuM. (2015). Regulatory Mechanism of mTOR Signaling Pathway on Cell Energy Metabolism. Sci. Sin. Vitae45, 1124–1131. 10.1360/n052015-00132

  • CrossRef
  • Google Scholar

• 57

  TangZ.QiuH.LuoL.LiuN.ZhongJ.KangK.et al (2017). miR‐34b Modulates Skeletal Muscle Cell Proliferation and Differentiation. J. Cel. Biochem.118, 4285–4295. 10.1002/jcb.26079

  • CrossRef
  • Google Scholar

• 58

  TianL.CaiD.ZhuangD.WangW.WangX.BianX.et al (2020). miR-96-5p Regulates Proliferation, Migration, and Apoptosis of Vascular Smooth Muscle Cell Induced by Angiotensin II via Targeting NFAT5. J. Vasc. Res.57, 86–96. 10.1159/000505457

  • CrossRef
  • Google Scholar

• 59

  TroyA.CadwalladerA. B.FedorovY.TynerK.TanakaK. K.OlwinB. B. (2012). Coordination of Satellite Cell Activation and Self-Renewal by Par-complex-dependent Asymmetric Activation of P38α/β MAPK. Cell Stem Cell11, 541–553. 10.1016/j.stem.2012.05.025

  • CrossRef
  • Google Scholar

• 60

  Van SteenbergenA.BalteauM.GinionA.FertéL.BattaultS.RavensteinC. d. M. d.et al (2017). Sodium-myoinositol Cotransporter-1, SMIT1, Mediates the Production of Reactive Oxygen Species Induced by Hyperglycemia in the Heart. Sci. Rep.7, 41166. 10.1038/srep41166

  • CrossRef
  • Google Scholar

• 61

  VienbergS.GeigerJ.MadsenS.DalgaardL. T. (2017). MicroRNAs in Metabolism. Acta Physiol.219, 346–361. 10.1111/apha.12681

  • CrossRef
  • Google Scholar

• 62

  WalkeW.StapleJ.AdamsL.GnegyM.ChahineK.GoldmanD. (1994). Calcium-dependent Regulation of Rat and Chick Muscle Nicotinic Acetylcholine Receptor (nAChR) Gene Expression. J. Biol. Chem.269, 19447–19456. 10.1016/s0021-9258(17)32189-0

  • CrossRef
  • Google Scholar

• 63

  WangJ.ZhaoL.PengX.LiuK.ZhangC.ChenX.et al (2020). Evaluation of miR‐130 Family Members as Circulating Biomarkers for the Diagnosis of Bladder Cancer. J. Clin. Lab. Anal.34, e23517. 10.1002/jcla.23517

  • CrossRef
  • Google Scholar

• 64

  WeiW.HeH.-B.ZhangW.-Y.ZhangH.-X.BaiJ.-B.LiuH.-Z.et al (2013). miR-29 Targets Akt3 to Reduce Proliferation and Facilitate Differentiation of Myoblasts in Skeletal Muscle Development. Cell Death Dis.4, e668. 10.1038/cddis.2013.184

  • CrossRef
  • Google Scholar

• 65

  WennogleL. P.HoxieH.PengY.HendrickJ. P. (2017). Phosphodiesterase 1: A Unique Drug Target for Degenerative Diseases and Cognitive Dysfunction. Adv. Neurobiol.17, 349–384. 10.1007/978-3-319-58811-7_13

  • CrossRef
  • Google Scholar

• 66

  WuN.GuT.LuL.CaoZ.SongQ.WangZ.et al (2019). Roles of miRNA‐1 and miRNA‐133 in the Proliferation and Differentiation of Myoblasts in Duck Skeletal Muscle. J. Cel Physiol.234, 3490–3499. 10.1002/jcp.26857

  • CrossRef
  • Google Scholar

• 67

  XieC.MaoX.HuangJ.DingY.WuJ.DongS.et al (2011). KOBAS 2.0: a Web Server for Annotation and Identification of Enriched Pathways and Diseases. Nucleic Acids Res.39, W316–W322. 10.1093/nar/gkr483

  • CrossRef
  • Google Scholar

• 68

  YamamoyoH.MorinoK.NishioY.UgiS.YoshizakiT.KashiwagiA.et al (2012). MicroRNA (miRNA)-494 Regulates Mitochondrial Biogenesis in Skeletal Muscle through Mitochondrial Transcriptional Factor A (mtTFA) and Forkhead Box j3 (Foxj3). Am. J. Physiol. Endocrinol. Metab.303, E1419. 10.1152/ajpendo.00097.2012

  • CrossRef
  • Google Scholar

• 69

  YanZ.JianK.QunL.YanW.JieL. (2018). Role of miRNAs in Skeletal Muscle Aging. Clin. Interv. Aging13, 2407–2419. 10.2147/CIA.S169202

  • CrossRef
  • Google Scholar

• 70

  YinH.HeH.ShenX.ZhaoJ.CaoX.HanS.et al (2020). miR-9-5p Inhibits Skeletal Muscle Satellite Cell Proliferation and Differentiation by Targeting IGF2BP3 through the IGF2-PI3K/Akt Signaling Pathway. Int. J. Mol. Sci.21, 1655. 10.3390/ijms21051655

  • CrossRef
  • Google Scholar

• 71

  ZhangD.LiY.YaoX.WangH.ZhaoL.JiangH.et al (2016). miR-182 Regulates Metabolic Homeostasis by Modulating Glucose Utilization in Muscle. Cel Rep.16, 757–768. 10.1016/j.celrep.2016.06.040

  • CrossRef
  • Google Scholar

• 72

  ZhangM.LiD.-H.LiF.SunJ.-W.JiangR.-R.LiZ.-J.et al (2018). Integrated Analysis of MiRNA and Genes Associated with Meat Quality Reveals that Gga-MiR-140-5p Affects Intramuscular Fat Deposition in Chickens. Cel Physiol. Biochem.46, 2421–2433. 10.1159/000489649

  • CrossRef
  • Google Scholar