Skip to main content

miR-1/206 downregulates splicing factor Srsf9 to promote C2C12 differentiation

Abstract

Background

Myogenesis is driven by specific changes in the transcriptome that occur during the different stages of muscle differentiation. In addition to controlled transcriptional transitions, several other post-transcriptional mechanisms direct muscle differentiation. Both alternative splicing and miRNA activity regulate gene expression and production of specialized protein isoforms. Importantly, disruption of either process often results in severe phenotypes as reported for several muscle diseases. Thus, broadening our understanding of the post-transcriptional pathways that operate in muscles will lay the foundation for future therapeutic interventions.

Methods

We employed bioinformatics analysis in concert with the well-established C2C12 cell system for predicting and validating novel miR-1 and miR-206 targets engaged in muscle differentiation. We used reporter gene assays to test direct miRNA targeting and studied C2C12 cells stably expressing one of the cDNA candidates fused to a heterologous, miRNA-resistant 3′ UTR. We monitored effects on differentiation by measuring fusion index, myotube area, and myogenic gene expression during time course differentiation experiments.

Results

Gene ontology analysis revealed a strongly enriched set of putative miR-1 and miR-206 targets associated with RNA metabolism. Notably, the expression levels of several candidates decreased during C2C12 differentiation. We discovered that the splicing factor Srsf9 is a direct target of both miRNAs during myogenesis. Persistent Srsf9 expression during differentiation impaired myotube formation and blunted induction of the early pro-differentiation factor myogenin as well as the late differentiation marker sarcomeric myosin, Myh8.

Conclusions

Our data uncover novel miR-1 and miR-206 cellular targets and establish a functional link between the splicing factor Srsf9 and myoblast differentiation. The finding that miRNA-mediated clearance of Srsf9 is a key myogenic event illustrates the coordinated and sophisticated interplay between the diverse components of the gene regulatory network.

Background

A complex network of integrated transcriptional and post-transcriptional regulatory mechanisms controls skeletal muscle gene expression. Overall, skeletal muscle displays one of the most tissue-specific splicing profiles, and changes in alternative splicing trigger proper temporal gene expression patterns during myogenesis [1,2,3,4,5,6,7]. Moreover, biochemical and biophysical properties of many components of the contractile machinery can be fine-tuned by selectively expressing specialized isoforms. Regulated splicing is orchestrated by the combinatorial interaction between cis-regulatory elements and trans-acting factors. This process efficiently creates functionally diverse proteins from a single gene; for instance, the fast troponin T primary transcript can be alternatively spliced into 64 functionally distinct isoforms [8]. The critical importance of alternative splicing in skeletal muscle is even more apparent when this process is misregulated, as occurs in a number of muscle disorders [1]. For instance, perturbations in alternative splicing ultimately result in an inability to transition from a fetal to an adult splicing pattern in myotonic dystrophy 1 [9, 10].

Members of both hnRNP and SR protein families play important roles during muscle development [11,12,13]. SR proteins constitute a highly conserved group of splicing factors required for both constitutive and alternative splicing [14, 15]. They are structurally characterized by the presence of an N-terminal RNA recognition motif (RRM) and a C-terminal arginine/serine-rich domain (RS domain), which interacts with other cellular factors during spliceosome assembly and splice site selection. A number of reports underscore the vital importance of certain SR protein in vivo. For example, global SR gene knockouts exhibit embryonic lethality [13, 16,17,18,19], and Srsf10 (Srp38) knockout embryos die mid- to late-gestation of multiple defects affecting both cardiac and skeletal muscle [13, 19]. Finally, while conditional cardiac Srsf1 knockout animals die postnatally from defective juvenile-to-adult heart remodeling [18], conditional cardiac knockouts of Srsf2 (SC35) and Srsf4 (Srp75) are viable but develop dilated and hypertrophic cardiomyopathies, respectively [20, 21].

While alternative splicing modulates the expression of different gene isoforms, miRNAs can control their cellular levels by base-pairing with miRNA response elements (MREs) present in target mRNAs. Formation of this duplex can subsequently block translation or trigger degradation of the mRNA. Importantly, disrupted miRNA processing in muscle leads to a perinatal lethal phenotype characterized by muscle hypoplasia, abnormal myofiber organization, and increased cell death [22]. Several specific miRNAs are enriched in skeletal muscle, including miR-1 and miR-206, which share a common seed sequence. These miRNAs suppress myoblast proliferation and promote muscle differentiation in both animal and cell culture models [23,24,25,26,27]. Both miR-1 and miR-206 are minimally expressed in proliferating myoblasts, but ectopic expression or inhibition can force or prevent differentiation, respectively [24,25,26,27]. Notably, miR-1 may also act as a key myogenic fate determinant as forced expression in HeLa cells shifts their transcriptome to a typical muscle profile [28].

miRNAs also control myogenesis by modulating the cellular concentration of several splicing factors. For example, miR-133a regulates nPTB, miR-30-5p targets MBNL, and miR-222 regulates Rbm24 [12, 29, 30]. In addition, the RNA helicases Ddx5 and Ddx17, which cooperate with the splicing factors hnRNP H/F and the transcription factor MyoD to regulate transcription of myogenic genes in C2C12 mouse myoblasts, are direct miR-1 and miR-206 targets [31].

Through a combination of in silico predictions and correlated myogenesis functional analyses in C2C12 cells, we report herein that miR-1 and miR-206 targeting of the transcript encoding the SR protein Srsf9 is a significant myogenic event. The data presented highlight the cellular relationship between miRNA expression and levels of splicing factors that control the temporal production of specialized isoforms during muscle differentiation.

Methods

RNA isolation, cDNA synthesis, and qPCR

Total RNA was isolated with TRI Reagent (MRC TR 118) according to the manufacturer’s protocol. Tissue samples were homogenized directly into TRI reagent using a dispersion tool (IKA T10 Basic S1) while cultured cells were washed with PBS and scraped directly into TRI reagent. All qPCR was performed on a Bio-Rad CFX96 attached to a C1000 thermocycler. miRNA expression was assessed by TaqMan-based qPCR (ThermoFisher, Part Number 4427975; Assay IDs: sno202 = 001232, miR-206 = 000510, miR-1 = 002222). miRNAs were reverse transcribed with a TaqMan MicroRNA Reverse Transcription Kit (ThermoFisher, 4366596) and qPCR was performed with TaqMan Universal PCR Master Mix, No AmpErase UNG (ThermoFisher, 4324018), all according to the manufacturer’s instructions. Relative expression was analyzed with the ΔΔCT method. mRNA expression was assessed with SYBR Green-based qPCR. All primer sequences are listed in Additional file 6: Table S1. Total RNA was reverse transcribed with random hexamer primers using the Superscript II Reverse Transcriptase kit (ThermoFisher, 18064-022). qPCR was performed with SYBR Green PCR Master Mix (ThermoFisher, 4312704) according to the manufacturer’s instructions. Relative expression was analyzed with the Pfaffl standard curve method.

RNA sequencing

Total RNA from 3 independent growing myoblast cultures before differentiation (BD1, BD2, BD3) and 3 independent day 6 myotube cultures after differentiation (AD1, AD2, AD3) was oligo dT-selected and reverse transcribed with random hexamers, followed by second strand synthesis, adapter ligation, and 5′ end selection. Libraries were sequenced with an Illumina HiSEQ. A total of 210 million paired end reads were aligned to the mm9 version of the mouse genome downloaded from UCSC Genome Browser. Additional file 7: Table S2 presents reads per million aligned. BAM files are available upon request.

miR-1/206 target selection and gene ontology analysis

TargetScan 7.2 was queried to retrieve the set of predicted miR-1/206 mouse target genes [32]. This was crossed with a list of all mRNAs in the RNA-Seq dataset that decreased by 1.3-fold or more during differentiation, regardless of statistical significance. These permissive criteria were chosen to be as inclusive as possible for potential miR-1/206 targets in C2C12. This list was analyzed with DAVID (Database for Annotation, Visualization and Integrated Discovery version 6.8 [33,34,35]) to uncover enriched gene ontology (GO) clusters of related gene sets. The functional clustering tool was used with default settings and the mouse genome as background. Clusters were ranked by descending enrichment score. Enrichment score is the –log10 of the geometric mean of the p values of all individual terms in the cluster. A cutoff of 0.05 or below for the geometric mean of p values, which corresponds to an enrichment score of greater than 1.3, was chosen.

Cloning and mutagenesis

All primer sequences are listed in Additional file 6: Table S1. miRNA overexpression constructs were generated by cloning either the stem loop region (miR-450a-1 and miR-1) or a larger 1 kilobase region encompassing the stem loop (miR-206; the stem loop sequence alone was insufficient for proper processing and subsequent targeting; data not shown) into pcDNA3.1(-) through standard molecular cloning techniques. miR-450a-1 and miR-1 stem loops were assembled with primers through overlap extension PCR and cloned into the EcoRI/BamHI sites of the plasmid while the 1 kb miR-206 locus was PCR-amplified from mouse genomic DNA and inserted between the HindIII and XhoI sites of the plasmid. 3′ UTR reporter constructs were generated in psiCheck2, a dual-luciferase plasmid where UTRs are cloned downstream of Renilla luciferase and firefly luciferase serves as an internal control. All sequences were cloned between the XhoI and NotI sites. Positive controls for miR-450a-1, miR-1, and miR-206 (2x450a-1, 2x1, and 2x206) were constructed by introducing 2 repeats (spaced by 2 A residues) of the mature miRNA antisense sequence downstream of Renilla luciferase using annealed complementary oligonucleotides. 3′ UTRs were PCR-amplified from mouse genomic DNA and inserted into the plasmid using standard molecular cloning techniques. The Srsf9 MRE was mutated by replacing the natural MRE with the reverse sequence through inverse PCR of the entire psiCheck2-Srsf9 wild-type 3′ UTR plasmid. iProof high fidelity DNA polymerase (Bio-Rad, 172-5302) was used according to the manufacturer’s instructions with phosphorylated primers. pEGFP-Srsf9 (a C-terminally GFP-tagged Srsf9 expression construct) was generated by amplifying Srsf9 cDNA from C2C12 myoblast cDNA and inserting it between the EcoRI and BamHI sites of pEGFP-N1 with standard molecular cloning techniques. A C-terminal GFP tag in a related vector, pEGFP-N3, has been shown to be tolerated by other SR proteins [36]. All clones were verified by Sanger sequencing.

Cell culture, transfection, and stable cell line generation

C2C12 cells were grown as myoblasts in growth medium (GM): high glucose DMEM (Invitrogen 11960069) supplemented with 20% fetal bovine serum, 2 mM l-glutamine, 100 U/mL penicillin and 100 μg/mL streptomycin, and 1 mM sodium pyruvate. Myotubes were differentiated by changing media to differentiation medium (DM): high glucose DMEM supplemented with 5% adult horse serum, 2 mM l-glutamine, 100 U/mL penicillin and 100 μg/mL streptomycin, and 1 mM sodium pyruvate. When differentiating, DM was refreshed every day to prevent media acidification. For luciferase assays, cells were plated in triplicate in 6-well dishes at a density of 50,000 cells/well 24 h before transfection. Cells were transfected 24 h after plating using the transfection reagent TransIT-LT1 (Mirus, MIR 2305) according to the manufacturer’s instructions. When ectopically expressing a miRNA along with a 3′ UTR-linked reporter gene, cells were harvested 24 h post-transfection. When differentiating cells for endogenous induction of miRNAs, day 0 time points were collected 24 h post-transfection and differentiation was triggered at the same time for later time points. At harvest, cells were washed twice with phosphate-buffered saline (PBS) solution and stored in an ultralow freezer in order to process all time points together at the end of the experiment. For stable cell line generation, C2C12 cells were transfected with plasmids (either pcDNA3.1(-), pEGFP-N1, or pEGFP-N1-Srsf9) using TransIT-LT1. Twenty-four hours post-transfection, G418 antibiotic selection was initiated by adding GM supplemented with 250 μg/mL G418. Selection was continued until control untransfected cells all died and cells were split whenever necessary to keep density below 60% to prevent spontaneous differentiation. Cell lines were maintained as pools to minimize potential genomic locus insertional effects.

Tissue collection

Animal work was reviewed by the University of Colorado Boulder Institutional Animal Care and Use Committee and approved under protocols 1002.07 and 1002.08. Wild-type C57Bl/6 mice were bred and housed at the University of Colorado Boulder under standard conditions. Adult animals were anesthetized by inhaled isoflurane and sacrificed by cervical dislocation followed by pneumothorax. Embryos were anesthetized on ice and sacrificed by decapitation. Whole hindlimbs were collected for embryonic time points and males and females were equally represented (sex was determined by PCR for the Sry gene using genomic DNA from the tail as template and primers against β-glucuronidase on chromosome 5 as positive control; see Additional file 6: Table S1 for primer sequences). Adult soleus was collected from 6-month-old male mice. Tissue samples were flash frozen in liquid nitrogen and stored in an ultralow freezer. Animal numbers for 15.5 dpc, 17.5 dpc, 19.5 dpc, and adults were 4, 6, 6, and 4, respectively.

Luciferase assays

A Dual-Luciferase Reporter Assay System (Promega, E1960) was used according to the manufacturer’s instructions with a Turner Designs TD-20/20 luminometer. Briefly, frozen cells were lysed in 1X Passive Lysis Buffer and reporter gene activities were measured first in LARII reagent (firefly luciferase internal control) and second in Stop & Glo reagent (Renilla luciferase experimental reporter gene). Read times for both were 10 s. Renilla/firefly ratios were compared.

Western blotting and in-gel fluorescent imaging

Cell lysates were prepared in ice cold RIPA buffer (50 mM Tris pH 8.0, 1 mM EDTA, 0.5 mM EGTA, 1% Triton X-100, 0.5% deoxycholate, 0.1% SDS, 140 mM NaCl) with 1X cOmplete EDTA-free Protease Inhibitor Cocktail (Roche, 11873580001). Total protein concentration was measured with a BCA protein assay kit (Pierce, 23250) and 20 μg total protein was resolved by 10% PAGE and transferred to nitrocellulose membrane. Standard western blotting techniques were employed with all antibodies diluted in TBS/0.1% tween-20/4% non-fat dry milk. 1:2000 α-GFP (Santa Cruz Biotechnology, sc-8334, rabbit polyclonal), 1:1000 α-Srsf9 (MBL, RN081PW, rabbit polyclonal), and 1:5000 α-Gapdh (Cell Signaling Technology, 2118S, rabbit monoclonal 14C10) were used as primary antibodies and horseradish peroxidase-linked goat anti-rabbit (Jackson Immunoresearch, 111-035-144) was used as secondary antibody. In-gel GFP autofluorescence was measured by directly imaging the gel with a Typhoon FLA 9500 (GE Healthcare). For differentiation time courses, equal loading was assessed by Ponceau S staining to reveal total protein.

Immunostaining and microscopy

Cells were seeded at a density of 100,000 cells per well in 6-well dishes on glass coverslips that had been silanized, gelatinized, and glutaraldehyde cross-linked to promote adhesion of myotubes. Cells were harvested by washing once in 37 °C PBS + 0.4% glucose and fixed 5′ at room temperature in 3% paraformaldehyde diluted in PBS. Paraformaldehyde was quenched with 50 mM ammonium chloride. Fixed cells were permeabilized 5′ at room temperature in PBS/0.1% triton X-100 (Permeabilization Solution (PS)) and blocked in Blocking Solution (BS: PS/1% BSA/1% normal goat serum) while rocking. Myosin heavy chain was probed with F59 as primary antibody (Developmental Studies Hybridoma Bank; custom production from hybridoma cells) and AlexaFluor568-linked goat anti-mouse IgG1 as secondary antibody (Invitrogen, A-21124) diluted 1:200. DNA was stained with 300 nM 4′,6-diamidino-2-phenylindole (DAPI; Sigma D9542). All imaging (for GFP, AlexaFluor568, and DAPI) was performed on a Nikon TE2000 inverted fluorescent microscope with a 20X objective connected to a Nikon DS-QiMc-U3 camera controlled through the NIS-Elements AR software version 4.00.03. Exposure times for a given channel were kept constant for all slides. Raw images were processed with the Image5D plugin in ImageJ using the same settings for a given channel across all images.

Nuclear fusion index and myotube area calculation

Fusion indices were calculated from the percent of total nuclei residing in syncytia (syncytium defined as a cell with 2 or more nuclei). Eight to 10 non-overlapping 20X fields of view were analyzed and averaged for each cell line. The number of nuclei in syncytia was tabulated and divided by the total nuclei in each field. Myotube area per field of view was calculated as follows. Red (myosin) and blue (DNA) channels for a given field of view were individually opened in ImageJ. Each was inverted and thresholded (the same threshold for a given channel was applied to all images) to create a binary black and white mask. The Analyze Particles function was run on each myosin and DNA image and the area of the myosin image was divided by the area of the DNA image (to normalize for potential differences in cell density across cell lines and across different fields of view). A representative set of thresholded image masks is displayed in Additional file 1: Figure S1.

Graphing and statistical analysis

All graphing was done with GraphPad Prism and statistical significance was assessed by one-way ANOVA with Dunnett’s post-test to compare with the control condition or Tukey’s post-test to compare all pairs as noted for multiple comparisons. Data are presented as means with error bars representing standard error of the mean. To assess potential differences in myogenin and perinatal myosin heavy chain expression over a differentiation time course, non-linear regression analysis was used to fit quadratic curves to the data points for GFP control and Srsf9-GFP cell lines. Curve fits were compared to determine whether they were statistically different. A p value cutoff of 0.05 was the minimum for significance. Other cutoffs are as noted in the figure legends.

Results

Splicing factors are enriched amongst predicted miR-1/206 targets and they decrease in expression during myotube formation

In order to identify functional categories of putative miR-1/206 targets in myogenesis, we took a combined in silico and in vitro cell culture model approach. We reasoned that miR-1/206 target genes would be expressed in myoblasts as they are poised to differentiate but that their concentrations would be reduced during the differentiation process by the increasing levels of miR-1/206. To this end, in both C2C12 myoblasts and day 6 myotubes, we initially measured the amounts of miR-1 and miR-206 by qPCR as well as the global changes in their mRNA transcriptomes by RNA-Seq (Additional file 7: Table S2). These analyses showed that while miR-1 and miR-206 levels are very low in myoblasts, both miRNAs are robustly expressed in differentiated myotubes, which have both the molecular and morphological characteristics of mature muscle (Additional file 2: Figure S2, panels A and B). This also confirms previous measurements of miR-1 and miR-206 in differentiating C2C12 cells [26]. Since miRNAs often exert mild effects on individual targets, we used the RNA-Seq dataset to assemble a list of genes downregulated by 1.3-fold or more during differentiation as this is the average fold change in mRNA targets observed for multiple miRNAs [37, 38]. We then retrieved the set of 896 predicted miR-1/206 targets from TargetScan and found that the levels of 354 decreased during myogenesis (Additional file 8: Table S3).

In order to evaluate whether particular functional categories are enriched in this gene set, we employed the DAVID functional annotation clustering tool, which groups functional categories with similar gene sets from multiple different annotation sources to capture biological themes. Of the 354 candidate C2C12 miR-1/206 targets, 351 had associated DAVID IDs. We found 21 clusters (Additional file 9: Table S4) with an enrichment score greater than 1.30 (indicating an average p value of below 0.05 for all terms in the cluster). While the top two clusters were related to DNA binding and transcriptional regulation and included known targets such as E2f5 [39], Pax3 [40], and Sox9 [41], both the third and fifth clusters included genes related to RNA metabolism. Notably, one quarter of the SR family of splicing factors appeared in these clusters. Since regulated splicing plays a central role in promoting myogenesis [4, 6, 7], we investigated whether miR-1/206 controls muscle differentiation by specifically targeting SR splicing factors.

mRNAs encoding SR protein family members are miR-1/206 targets during myogenesis

Our bioinformatics analysis predicted mRNAs encoding four of these proteins, Srsf1 (ASF/SF2), Srsf3 (Srp20), Srsf9 (Srp30c), and the SR-related protein Tra2b, as miR-1/206 targets. Since Srsf1 has an established critical role in heart development [18], we focused our attention on the remaining candidates. We assayed miRNA targeting in C2C12 myoblasts, which, as already noted, have extremely low levels of both miR-1 and miR-206. We co-transfected cells with constructs expressing a reporter Renilla luciferase gene fused to the mouse 3′ UTRs of our RNA candidates and either miR-1 or miR-206 overexpression plasmids. We determined the specificity of miR-1 or miR-206 by measuring the targeting activity of miR-450a-1, a non-myogenic miRNA which does not have overlapping predicted targets with miR-1/206. Luciferase quantification showed that this miRNA did not cross-react with the selected SR 3′ UTRs or with 3′ UTRs carrying tandem copies of the complete reverse complement of either miR-206 (2x206) or miR-1 (2x1) (data not shown). However, it efficiently downregulated the positive control construct 2x450a-1, constructed similarly to 2x206 and 2x1 (Additional file 3: Figure S3). To determine the sensitivity of the assay, we next assessed the activity of miR-1/206 on the 2x206 and 2x1 constructs as well as on the 3′ UTR of Ccnd1, previously identified as a miR-206 target in C2C12 cells [42]. We also monitored promoter interference causing potential transcriptional repression of the reporter construct using a miRless construct, which expresses the Renilla luciferase RNA without miRNA target sites. Figure 1a shows the result of this analysis with each luciferase signal measured in the presence of miR-1 or miR-206 normalized to the same constructs co-transfected with the control miRNA, miR-450a-1 (see “Methods”). The data, presented as fold change versus miR-450a-1 control, show that miRless expression did not change appreciably in the presence of miR-1 or miR-206 and that miR-1 and miR-206 efficiently targeted 2x1 and 2x206 along with the positive control Ccnd1 (Fig. 1a). Furthermore, both miR-1 and miR-206 also reduced luciferase activity by targeting the 3′ UTRs of Srsf9 and Tra2b fused to the reporter gene while they did not have any effect on the Srsf3 construct (Fig. 1b). While statistically significant, the negative regulation of the Tra2b 3′ UTR was modest in magnitude. Thus, we focused our next set of experiments only on Srsf9 activity and its potential role in muscle differentiation. We substantiated the specific miR-1/206 targeting of the Srsf9 3′ UTR by MRE mutagenesis. When we reversed the orientation of the predicted MRE in the Srsf9 3′ UTR, which preserves positioning of any unrecognized flanking elements, this mutant construct (termed Srsf9 MRE Rev) restored reporter gene activity to levels measured in the presence of the control miR-450a-1 (Fig. 1c). Moreover, endogenous Srsf9 mRNA expression decreased in myoblasts when we transfected miR-1 or miR-206 overexpression plasmids (Additional file 4: Figure S4 A). Taken together, these results establish that expression of the miR-1/206 family can directly modulate the level of Srsf9 in C2C12 myoblasts.

Fig. 1
figure 1

Srsf9 is a direct target of miR-1 and miR-206. a miR-206 and miR-1 target 3′ UTRs with cognate MREs in C2C12 myoblasts. Cells were transfected with reporter constructs and expression constructs for miR-206, miR-1, or the control miRNA, miR-450a-1. Renilla luciferase linked to two tandem copies of the reverse complement of miR-206 or miR-1 (2x206 and 2x1, respectively) or to the 3′ UTR of known miR-1/206 target Ccnd1 were positive controls. A 3′ UTR with no MREs in the empty reporter construct (miRless) was the negative control. Cells were harvested 24 h after transfection. The graph presents fold changes (FC) in firefly-normalized Renilla luciferase activity for miR-206- or miR-1-transfected vs. miR-450a-1-transfected cells. The dashed line at y = 1 represents miR-450a-1-transfected levels. **p ≤ 0.01; ***p ≤ 0.001 vs. miR-450a-1. b 3′ UTRs from Srsf9 and Tra2b but not Srsf3 are sensitive to miR-206 and miR-1 activity. C2C12 myoblasts were transfected with reporter constructs containing the Srsf3, Srsf9, or Tra2b 3′ UTRs along with miRNA expression constructs. Data were collected and analyzed as in (a). c Srsf9 is a direct target of miR-1/206. The predicted MRE in the Srsf9 3′ UTR was reversed in the reporter construct to abolish miR-1/206 binding. C2C12 myoblasts were transfected and analyzed as in panel (a). d The Srsf9 3′ UTRs are sensitive to endogenous myogenic cues. 3′ UTR reporter constructs were transfected into C2C12 myoblasts which were then differentiated for 0, 1, 2, or 4 days. Firefly-normalized Renilla luciferase activity was measured. Fold changes for differentiating vs. day 0 myoblasts were calculated and normalized to the miRless control. The dashed line at y = 1 represents day 0 levels. *p ≤ 0.05; **p ≤ 0.01 vs. D0 For all panels, N = 3 independent cultures

To broaden the data obtained in myoblasts, we next surveyed the activity of the reporter construct during C2C12 differentiation which exposes the Srsf9 3′ UTR to physiological concentrations of both miR-1 and miR-206. Time course analysis showed an inversely proportional relationship between the decreased reporter luciferase activity and the concomitant increase in miR-1 and miR-206 expression occurring during the activation of the C2C12 differentiation program (Fig. 1d; Additional file 2: Figure S2). Interestingly, we observed a similar pattern for the Ccnd1 reporter construct (Fig. 1d). Consistently, the mRNA expression levels of Srsf9 and Ccnd1 measured 6 days after differentiation were significantly reduced (Fig. 2a), which is also in accordance with the RNA-Seq data. As expected, Srsf9 protein levels correspondingly decreased during differentiation (Additional file 4: Figure S4 B). We also queried a quantitative pSILAC mass spectrometry dataset that compiles protein level changes in HeLa cells 8 h after miR-1 overexpression. The median fold change for Srsf9 peptides was − 3.71, which supports direct miR-1 targeting [37, 43].

Fig. 2
figure 2

Srsf9 mRNA levels are inversely proportional to miR-1/206 levels. a Srsf9 mRNA levels decrease in differentiating C2C12 cells. Expression was assessed by qPCR in growing myoblasts (day 0) and in day 6 myotubes. Gapdh was the reference gene. Expression levels in the differentiating samples were normalized to myoblast levels. Myogenin (MyoG) was the positive control as its expression increases during differentiation. *p ≤ 0.05; **p ≤ 0.01 vs. D0 N = 3 independent cultures for all genes and time points. b miR-206 and miR-1 levels increase during mouse hindlimb skeletal muscle development. Expression was assessed by qPCR and normalized to sno202 in whole hindlimb for developmental time points and in the soleus from 6-month-old male mice. Fold change relative to 15.5 dpc is presented. N for 15.5 dpc, 17.5 dpc, 19.5 dpc, and adult was 4, 6, 6, and 4, respectively. ***p ≤ 0.001 vs. 15.5 dpc. c Srsf9 and Ccnd1 mRNAs decrease during mouse hindlimb skeletal muscle development. Expression was assessed by qPCR with 18S rRNA as the reference gene. Expression levels were normalized to 15.5 dpc. *p ≤ 0.05; **p ≤ 0.01 vs 15.5 dpc Ns were the same as in panel (b)

We extended this analysis to embryonic limb formation in the developing mouse. As shown in Fig.2b, we found that both miR-1 and miR-206 expression, which increased between 3- and 5-fold as embryogenesis progressed, peaked during the transition to mature adult muscle. In agreement with the inverse correlation observed in the C2C12 analysis, the mRNA levels of Srsf9 and Ccnd1 showed a modest decline in expression from embryonic day 15.5 to day 19.5 but a marked decrease from the embryonic to adult mouse transition (Fig. 2c).

The inability to downregulate Srsf9 expression during myogenesis results in impaired differentiation

Finally, we investigated whether the downregulation of Srsf9 expression plays an important role in controlling the sequential myogenic pathway. To this end, we generated C2C12 stable cell lines expressing a GFP-tagged version of Srsf9 fused to a heterologous 3′ UTR containing the SV40 polyA signal as well as GFP only or empty vector controls. After selection, we maintained pools of stable transfectants and specifically measured transgene RNA and protein expression during differentiation. qPCR and western blot analysis carried out with GFP-specific primers and antibody revealed that both Srsf9-GFP and GFP were stably expressed in our cell pools (Additional file 5: Figure S5 A, C). qPCR with Srsf9 gene-specific primers indicated that the composite expression level of endogenous Srsf9 and Srsf9-GFP was only 1.7-fold higher than in control myoblasts (Additional file 5: Figure S5 B). More importantly, time course analysis showed that the presence of the heterologous 3′ UTR stabilized the level of Srsf9-GFP mRNA and protein throughout cell differentiation (Additional file 5: Figure S5 D, E). Moreover, we found that, as revealed by co-localization with DAPI stain, Srsf9-GFP was properly localized in the nucleus of the majority of the cells, while the GFP control showed diffuse fluorescence signal (Fig. 3a). Interestingly, Srsf9-GFP was occasionally diffuse in the cytoplasm of myotubes; other SR proteins are known to shuttle and play extranuclear roles, particularly in the absence of their nuclear RNA targets [44]. To determine potential functional repercussions caused by sustained Srsf9 expression, we first calculated the fusion index, a method frequently used to quantify the extent of C2C12 differentiation. Scoring of the number of nuclei residing in a fused syncytium (defined as a cell with 2 or more nuclei) at day 6 of differentiation showed that the average fusion index of cells expressing Srsf9-GFP was significantly lower than the ones expressing either GFP or the empty vector control (Fig. 3b). Since Srsf9-GFP differentiated myotubes imaged at day 6 appeared to have reduced width compared with controls, we measured whether the myotube-positive area of each field of view was statistically smaller when normalized to DNA-positive area, a control for potential differences in cell density. After imaging the different cell pools stained with DAPI and a myosin heavy chain antibody (considered a terminal marker of muscle differentiation), we created cell masks by imposing equal intensity thresholds to the acquired images and then divided the myosin area by the DAPI area derived from the binary derivatives. As shown in Fig. 3c, the stable expression of the Srsf9-GFP construct, which is refractory to miR-1/206 targeting activity, resulted in a significant reduction of myotube area, a clear indication of impaired differentiation. Supporting this observation, induction of the perinatal myosin isoform (Myh8), the most expressed myosin in mammalian skeletal muscle during the early perinatal period [45] and one of the most abundant myosins expressed in differentiated C2C12 (Additional file 7: Table S2), was significantly blunted in Srsf9-GFP cells (Fig. 3d). This is indicated not only by decreased expression at individual time points but also by statistically different curves fit to each series when comparing GFP control and Srsf9-GFP cells. Finally, as the myogenesis program is fundamentally orchestrated at the transcriptional level, we measured induction of the muscle regulatory factor myogenin which is robustly expressed at early myogenic time points and then decreases as myotubes approach terminal differentiation. Although both the GFP controls and the Srsf9-GFP cells induced myogenin relative to D0 growing myoblast expression levels, this induction was much weaker in the Srsf9-GFP cells (Fig. 3e). Since the series could still be fit to a quadratic curve, with myogenin expression in Srsf9-GFP cells peaking mid-time course and then declining at day 6, we believe the data do not support a simple delay of differentiation. Taken together, these data strongly suggest that proper progression throughout the steps of muscle differentiation requires temporal Srsf9 downregulation.

Fig. 3
figure 3

mir-1/206-resistant Srsf9-GFP expression blunts myogenic differentiation. a C2C12 cell lines stably expressing GFP alone (GFP control) or GFP-tagged Srsf9 linked to a heterologous miR-1/206-resistant 3′ UTR (Srsf9-GFP) were differentiated for 6 days, immunostained, and imaged (20X objective). Myosin heavy chain, a marker of terminal differentiation, was visualized with the F59 antibody and an AlexaFluor568-linked secondary antibody. DNA was visualized by DAPI staining and GFP or Srsf9-GFP was visualized by GFP autofluorescence. Scale bars denote 6.5 μm. b Srsf9-GFP persistence reduces fusion. The fusion index of Srsf9-GFP cells was significantly lower than either pC-Empty cells (an empty expression plasmid) or GFP control cells. Eight to 10 non-overlapping fields of view were analyzed for each cell line. *p < 0.05 vs. pC-Empty; ###p < 0.001 vs. GFP control. c Srsf9-GFP persistence reduces myotube size. Myotube area was assessed by calculating the myosin-positive area of D6 myotubes. The myosin-positive area/DAPI-positive area was calculated for 8–10 non-overlapping fields of view. Myosin-positive area was significantly lower than both pC-Empty and GFP control cell lines. *p ≤ 0.05 vs. pC-Empty; #p < 0.05 vs. GFP control. d Srsf9-GFP persistence reduces myosin mRNA expression. mRNA levels of perinatal myosin (Myh8) were assessed by qPCR at days 0, 1, 2, 4, and 6 of differentiation. Statistically different quadratic curves were fit to GFP control and Srsf9-GFP series. e Srsf9-GFP persistence reduces MyoG induction. MyoG mRNA levels were measured by qPCR as in panel (d). Statistically different quadratic curves were fit to GFP control and Srsf9-GFP. For d and e, asterisks next to the Srsf9-GFP legend denote p < 0.0001 for the curve fit comparison

Discussion

In this report, we used a computational analysis of differentiating C2C12 myoblasts to identify miR-1/206 targets relevant to myogenesis. We found that the 3′ UTRs of several RNA-binding proteins, which belong to the SR splicing factor family, were highly enriched in our bioinformatics assessment. Accordingly, we discovered that downregulation of one of its members, Srsf9, controls myoblast differentiation. One limitation of our computational method is that it predicts miR-1/206 targets based on the decrease of their RNA levels; targets that are translationally inhibited but not regulated at the mRNA level, such as Hdac4 and Igf1 [25, 46, 47], escaped our analysis. However, several genome-wide studies on miRNAs, which include miR-1, have shown good correlation between RNA and protein levels [28, 38].

miR-1 and miR-206 are highly conserved members of the myomiR family that also includes miR-133a/b, miR-208a/b, miR-486, and miR-499. This ensemble of miRNAs, specifically expressed in cardiac and skeletal muscle, governs muscle differentiation, maintenance, and plasticity [48]. miR-206 global knockout (KO) mice have no reported gross muscle defects at baseline, but an apparent protective role during regeneration from acute or chronic muscle damage has emerged from studies where the miR-206KO mouse is (i) injured with cardiotoxin, (ii) crossed with the mdx muscular dystrophy mouse, or (iii) crossed with a mouse amyotrophic lateral sclerosis model [23, 46]. Interestingly, a different genetic model of miR-206 loss of function that deleted not only miR-206 but also miR-133b and a portion of the long non-coding RNA linc-MD1 (whose expression is normally mutually exclusive to miR-133b and can function as a miR-133 sponge) has no exaggerated response in the face of muscle damage [49]. While this contradicts the previously described studies, it is worth noting that miR-1/206 and miR-133 have been reported to have opposing roles: miR-1/206 promotes differentiation while miR-133 promotes proliferation [25, 26]. miR-1-1/1-2 double KO mice die shortly after birth from severe cardiac defects and so it is unclear whether they have a skeletal muscle phenotype [50].

To date, several reports have shown that miR-1/206 controls the expression levels of genes implicated in transcription (Pax3, Pax7, Hdac4 [24, 40, 46, 51, 52]) and cellular proliferation (Pola1 [26]). In this study, we found that downregulation of the splicing factor Srsf9 by miR-1/206 targeting is also essential for proper C2C12 differentiation. Srsf9 cooperates with several other RNA-binding proteins to repress or activate regulated splicing. It modulates the inclusion of SMN exon 7 [53], the alternative splicing of glucocorticoid receptor beta and gonadotropin-releasing hormone [54, 55], inclusion/exclusion of tau exons 2 and 10 [56, 57], and generation of the hnRNP A1B isoform [58, 59]. Notably, the main hnRNP A1 isoform regulates the alternative splicing of pyruvate kinase-M (PK-M). While the PK-M1 isoform is highly enriched in differentiated myotubes, myoblasts express the PK-M2 isoform which gives them a proliferative advantage [60]. Thus, it is tempting to speculate that, by controlling the shift in pyruvate kinase isoform expression, the combined activity of Srsf9 and miR-1/206 during differentiation lays the groundwork for metabolic adaptation. Moreover, it has been proposed that Srsf9 can interact cooperatively or compete with other splicing regulators for binding to high-affinity sites present in alternatively spliced pre-mRNAs [59]. Therefore, even a small decrease in Srsf9 cellular concentration could lead to profound changes in global splicing patterns. The SR protein cellular function is still not completely defined. Global mapping of RNA targets determined by uv cross-linking and immunoprecipitation (CLIP/iCLIP) has been carried out for Srsf1, 2, 3, and 4, and the eight drosophila SR homologs [36, 61, 62] but not for Srsf9. These studies performed in mouse embryo fibroblasts (MEFs) and P19 embryonal carcinoma cells revealed intriguing findings. For example, (i) specific alternative splicing patterns apparently controlled by a single SR protein actually hinge upon an intricate network of relationships with many other SR proteins; (ii) SR proteins can control gene expression by binding intronless transcripts as well as ncRNAs. Since no such map is available for Srsf9, it is difficult to predict through which molecular pathways Srsf9 exerts control over myoblast differentiation. Srsf9 high-affinity binding sites have been identified from a randomized pool of RNA sequences by a SELEX approach [63]. However, a computational survey of SELEX-derived consensus sequences, which are short and degenerate, revealed frequent overrepresentation in RNAs [64].

miR-1 and miR-206 have both been proposed as tumor suppressors. Interestingly, Srsf9 depletion reduces viability of prostate cancer cells [65]. Moreover, in two human bladder cancer cell lines, Srsf9 downregulation by miR-1 overexpression resulted in significant reduction in cell proliferation, migration, and invasion [66]. Accordingly, siRNA-mediated Srsf9 knockdown promoted apoptosis [66, 67]. These data strongly suggest that in bladder cancer cells, the tumor suppressive activity of miR-1 triggers apoptosis through direct Srsf9 inhibition. However, we did not observe any impact on cell survival when Srsf9 expression was stabilized during C2C12 differentiation.

Conclusions

We report here that miR-1/206 targets the Srsf9 3′ UTR during C2C12 myoblast differentiation. This is the first report showing a direct correlation between this member of the SR protein family and muscle maturation. We demonstrate that Srsf9 downregulation is necessary to achieve robust cell fusion and muscle-specific gene expression. Based on the complex web of functional interactions occurring amongst the RNA splicing proteins, we suggest that the persistence of Srsf9 alters several alternative splicing events and impairs proper production of specific protein isoforms driving myoblast maturation. The data presented also emphasize the importance of maintaining appropriate miRNA-mediated gene repression in adult muscle to avoid re-expression of fetal genes that could cause muscle disease.

Availability of data and materials

All data generated or analyzed during this study are included in this published article and its supplementary files.

Abbreviations

DAPI:

4′,6-Diamidino-2-phenylindole

DAVID:

Database for Annotation, Visualization, and Integrated Discovery

dpc:

Days post-coitum

FC:

Fold change

GFP:

Green fluorescent protein

MRE:

miRNA response element

PBS:

Phosphate-buffered saline

qPCR:

Quantitative PCR

RNA-Seq:

RNA sequencing

SR protein:

Serine/arginine-rich protein

TBS:

Tris-buffered saline

UTR:

Untranslated region

References

  1. Nakka K, Ghigna C, Gabellini D, Dilworth FJ. Diversification of the muscle proteome through alternative splicing. Skelet Muscle. 2018;8:8.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  2. Castle JC, Zhang C, Shah JK, Kulkarni AV, Kalsotra A, Cooper TA, et al. Expression of 24,426 human alternative splicing events and predicted cis regulation in 48 tissues and cell lines. Nat Genet. 2008;40:1416–25.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Llorian M, Smith CW. Decoding muscle alternative splicing. Curr Opin Genet Dev. 2011;21:380–7.

    Article  CAS  PubMed  Google Scholar 

  4. Bland CS, Wang ET, Vu A, David MP, Castle JC, Johnson JM, et al. Global regulation of alternative splicing during myogenic differentiation. Nucleic Acids Res. 2010;38:7651–64.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Tapial J, Ha KCH, Sterne-Weiler T, Gohr A, Braunschweig U, Hermoso-Pulido A, et al. An atlas of alternative splicing profiles and functional associations reveals new regulatory programs and genes that simultaneously express multiple major isoforms. Genome Res. 2017;27:1759–68.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, Van Baren MJ, et al. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010;28:511–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Brinegar AE, Xia Z, Loehr JA, Li W, Rodney GG, Cooper TA. Extensive alternative splicing transitions during postnatal skeletal muscle development are required for calcium handling functions. Elife. 2017;6:e27192.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Breitbart RE, Nguyen HT, Medford RM, Destree AT, Mahdavi V, Nadal-Ginard B. Intricate combinatorial patterns of exon splicing generate multiple regulated troponin T isoforms from a single gene. Cell. 1985;41:67–82.

    Article  CAS  PubMed  Google Scholar 

  9. Batra R, Charizanis K, Manchanda M, Mohan A, Li M, Finn DJ, et al. Loss of MBNL leads to disruption of developmentally regulated alternative polyadenylation in RNA-mediated disease. Mol Cell. 2014;56:311–22.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Han H, Irimia M, Ross PJ, Sung HK, Alipanahi B, David L, et al. MBNL proteins repress ES-cell-specific alternative splicing and reprogramming. Nature. 2013;498:241–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Hall MP, Nagel RJ, Fagg WS, Shiue L, Cline MS, Perriman RJ, et al. Quaking and PTB control overlapping splicing regulatory networks during muscle cell differentiation. RNA. 2013;19:627–38.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Boutz PL, Chawla G, Stoilov P, Black DL. MicroRNAs regulate the expression of the alternative splicing factor nPTB during muscle development. Genes Dev. 2007;21:71–84.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Wei N, Cheng Y, Wang Z, Liu Y, Luo C, Liu L, et al. SRSF10 plays a role in myoblast differentiation and glucose production via regulation of alternative splicing. Cell Rep. 2015;13:1647–57.

    Article  CAS  PubMed  Google Scholar 

  14. Long JC, Caceres JF. The SR protein family of splicing factors: master regulators of gene expression. Biochem J. 2009;417:15–27.

    Article  CAS  PubMed  Google Scholar 

  15. Zhou Z, Fu X-D. Regulation of splicing by SR proteins and SR protein-specific kinases. Chromosoma. 2013;122:191–207.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Jumaa H, Wei G, Nielsen PJ. Blastocyst formation is blocked in mouse embryos lacking the splicing factor SRp20. Curr Biol. 1999;9:899–902.

    Article  CAS  PubMed  Google Scholar 

  17. Wang HY, Xu X, Ding JH, Bermingham JR, Fu XD. SC35 plays a role in T cell development and alternative splicing of CD45. Mol Cell. 2001;7:331–42.

    Article  CAS  PubMed  Google Scholar 

  18. Xu X, Yang D, Ding JH, Wang W, Chu PH, Dalton ND, et al. ASF/SF2-regulated CaMKIIδ alternative splicing temporally reprograms excitation-contraction coupling in cardiac muscle. Cell. 2005;120:59–72.

    Article  CAS  PubMed  Google Scholar 

  19. Feng Y, Valley MT, Lazar J, Yang AL, Bronson RT, Firestein S, et al. SRp38 Regulates alternative splicing and is required for Ca(2+) handling in the embryonic heart. Dev Cell. 2009;16:528–38.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Ding JH, Xu X, Yang D, Chu PH, Dalton ND, Ye Z, et al. Dilated cardiomyopathy caused by tissue-specific ablation of SC35 in the heart. EMBO J. 2004;23:885–96.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Larrasa-Alonso J, Villalba M, Sanchez-Cabo F, Marti-Gomez C, Ortiz-Sanchez P, Lopez-Olaneta MM, et al. Loss of SRSF4 in cardiomyocytes induces hypertrophy, diastolic dysfunction and risk of sudden death. Eur Heart J. 2017;38:411.

    Article  Google Scholar 

  22. O’Rourke JR, Georges SA, Seay HR, Tapscott SJ, McManus MT, Goldhamer DJ, et al. Essential role for Dicer during skeletal muscle development. Dev Biol. 2007;311:359–68.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  23. Liu N, Williams AH, Maxeiner JM, Bezprozvannaya S, Shelton JM, Richardson JA, et al. microRNA-206 promotes skeletal muscle regeneration and delays progression of Duchenne muscular dystrophy in mice. J Clin Invest. 2012;122:2054–65.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Chen JF, Tao Y, Li J, Deng Z, Yan Z, Xiao X, et al. microRNA-1 and microRNA-206 regulate skeletal muscle satellite cell proliferation and differentiation by repressing Pax7. J Cell Biol. 2010;190:867–79.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Chen JF, Mandel EM, Thomson JM, Wu Q, Callis TE, Hammond SM, et al. The role of microRNA-1 and microRNA-133 in skeletal muscle proliferation and differentiation. Nat Genet. 2006;38:228–33.

    Article  CAS  PubMed  Google Scholar 

  26. Kim HK, Lee YS, Sivaprasad U, Malhotra A, Dutta A. Muscle-specific microRNA miR-206 promotes muscle differentiation. J Cell Biol. 2006;174:677–87.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Goljanek-Whysall K, Pais H, Rathjen T, Sweetman D, Dalmay T, Munsterberg A. Regulation of multiple target genes by miR-1 and miR-206 is pivotal for C2C12 myoblast differentiation. J Cell Sci. 2012;125:3590–600.

    Article  CAS  PubMed  Google Scholar 

  28. Lim LP, Lau NC, Garrett-Engele P, Grimson A, Schelter JM, Castle J, et al. Microarray analysis shows that some microRNAs downregulate large numbers of target mRNAs. Nature. 2005;433:769–73.

    Article  CAS  PubMed  Google Scholar 

  29. Zhang B-W, Cai H-F, Wei X-F, Sun J-J, Lan X-Y, Lei C-Z, et al. mir-30-5p regulates muscle differentiation and alternative splicing of muscle-related genes by targeting MBNL. Int J Mol Sci. 2016;17:182.

    Article  PubMed Central  CAS  Google Scholar 

  30. Cardinali B, Cappella M, Provenzano C, Garcia-Manteiga JM, Lazarevic D, Cittaro D, et al. MicroRNA-222 regulates muscle alternative splicing through Rbm24 during differentiation of skeletal muscle cells. Cell Death Dis. 2016;7:e2086.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Dardenne E, Polay Espinoza M, Fattet L, Germann S, Lambert MP, Neil H, et al. RNA helicases DDX5 and DDX17 dynamically orchestrate transcription, miRNA, and splicing programs in cell differentiation. Cell Rep. 2014;7:1900–13.

    Article  CAS  PubMed  Google Scholar 

  32. TargetScan. www.targetscan.org/vert_72/. Accessed 20 Mar 2019.

  33. DAVID Bioinformatics Resources 6.8 https://david.ncifcrf.gov/home.jsp. Accessed 22 Mar 2019.

  34. Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4:44–57.

    Article  CAS  Google Scholar 

  35. Huang DW, Sherman BT, Lempicki RA. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009;37:1–13.

    Article  CAS  Google Scholar 

  36. Änkö ML, Müller-McNicoll M, Brandl H, Curk T, Gorup C, Henry I, et al. The RNA-binding landscapes of two SR proteins reveal unique functions and binding to diverse RNA classes. Genome Biol. 2012;13:R17.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  37. Selbach M, Schwanhäusser B, Thierfelder N, Fang Z, Khanin R, Rajewsky N. Widespread changes in protein synthesis induced by microRNAs. Nature. 2008;455:58–63.

    Article  CAS  PubMed  Google Scholar 

  38. Hendrickson DG, Hogan DJ, McCullough HL, Myers JW, Herschlag D, Ferrell JE, et al. Concordant regulation of translation and mRNA abundance for hundreds of targets of a human microRNA. PLoS Biol. 2009;7:e1000238.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  39. Zhang X, Zhang E, Ma Z, Pei R, Jiang M, Schlaak JF, et al. Modulation of hepatitis B virus replication and hepatocyte differentiation by MicroRNA-1. Hepatology. 2011;53:1476–85.

    Article  CAS  PubMed  Google Scholar 

  40. Boutet SC, Cheung TH, Quach NL, Liu L, Prescott SL, Edalati A, et al. Alternative polyadenylation mediates microRNA regulation of muscle stem cell function. Cell Stem Cell. 2012;10:327–36.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Zhang Y-J, Xu F, Zhang Y-J, Li H-B, Han J-C, Li L. miR-206 inhibits non small cell lung cancer cell proliferation and invasion by targeting SOX9. Int J Clin Exp Med. 2015;8:9107–13.

    CAS  PubMed  PubMed Central  Google Scholar 

  42. Alteri A, De Vito F, Messina G, Pompili M, Calconi A, Visca P, et al. Cyclin D1 is a major target of miR-206 in cell differentiation and transformation. Cell Cycle. 2013;12:3781–90.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. pSILAC. https://psilac.mdc-berlin.de/proteingroups/906/3/. Accessed 23 Apr 2019.

  44. Twyffels L, Gueydan C, Kruys V. Shuttling SR proteins: more than splicing factors. FEBS J. 2011;278:3246–55.

    Article  CAS  PubMed  Google Scholar 

  45. Gokhin DS, Ward SR, Bremner SN, Lieber RL. Quantitative analysis of neonatal skeletal muscle functional improvement in the mouse. J Exp Biol. 2008;211:837–43.

    Article  CAS  PubMed  Google Scholar 

  46. Williams AH, Valdez G, Moresi V, Qi X, McAnally J, Elliott JL, et al. MicroRNA-206 delays ALS progression and promotes regeneration of neuromuscular synapses in mice. Science (80- ). 2009;326:1549–54.

    Article  CAS  Google Scholar 

  47. Shan ZX, Lin QX, Fu YH, Deng CY, Zhou ZL, Zhu JN, et al. Upregulated expression of miR-1/miR-206 in a rat model of myocardial infarction. Biochem Biophys Res Commun. 2009;381:597–601.

    Article  CAS  PubMed  Google Scholar 

  48. Horak M, Novak J, Bienertova-Vasku J. Muscle-specific microRNAs in skeletal muscle development. Dev Biol. 2016;410:1–13.

    Article  CAS  PubMed  Google Scholar 

  49. Boettger T, Wüst S, Nolte H, Braun T. The miR-206/133b cluster is dispensable for development, survival and regeneration of skeletal muscle. Skelet Muscle. 2014;4:23.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  50. Heidersbach A, Saxby C, Carver-Moore K, Huang Y, Ang YS, de Jong PJ, et al. MicroRNA-1 regulates sarcomere formation and suppresses smooth muscle gene expression in the mammalian heart. Elife. 2013;2:e01323.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  51. Hirai H, Verma M, Watanabe S, Tastad C, Asakura Y, Asakura A. MyoD regulates apoptosis of myoblasts through microRNA-mediated down-regulation of Pax3. J Cell Biol. 2010;191:347–65.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Dey BK, Gagan J, Dutta A. miR-206 and -486 induce myoblast differentiation by downregulating Pax7. Mol Cell Biol. 2011;31:203–14.

    Article  CAS  PubMed  Google Scholar 

  53. Young PJ, DiDonato CJ, Hu D, Kothary R, Androphy EJ, Lorson CL. SRp30c-dependent stimulation of survival motor neuron (SMN) exon 7 inclusion is facilitated by a direct interaction with hTra2 beta 1. Hum Mol Genet. 2002;11:577–87.

    Article  CAS  PubMed  Google Scholar 

  54. Xu Q, Leung DYM, Kisich KO. Serine-arginine-rich protein p30 directs alternative splicing of glucocorticoid receptor pre-mRNA to glucocorticoid receptor β in neutrophils. J Biol Chem. 2003;278:27112–8.

    Article  CAS  PubMed  Google Scholar 

  55. Park E, Han J, Son GH, Lee MS, Chung S, Park SH, et al. Cooperative actions of Tra2α with 9G8 and SRp30c in the RNA splicing of the gonadotropin-releasing hormone gene transcript. J Biol Chem. 2006;281:401–9.

    Article  CAS  PubMed  Google Scholar 

  56. Kondo S, Yamamoto N, Murakami T, Okumura M, Mayeda A, Imaizumi K. Tra2β, SF2/ASF and SRp30c modulate the function of an exonic splicing enhancer in exon 10 of tau pre-mRNA. Genes to Cells. 2004;9:121–30.

    Article  CAS  PubMed  Google Scholar 

  57. Wang Y, Wang J, Gao L, Lafyatis R, Stamm S, Andreadis A. Tau exons 2 and 10, which are misregulated in neurodegenerative diseases, are partly regulated by silencers which bind a SRp30c·SRp55 complex that either recruits or antagonizes htra2β1. J Biol Chem. 2005;280:14230–9.

    Article  CAS  PubMed  Google Scholar 

  58. Simard MJ, Chabot B. SRp30c is a repressor of 3’ splice site utilization. Mol Cell Biol. 2002;22:4001–10.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. Paradis C, Cloutier P, Shkreta L, Toutant J, Klarskov K, Chabot B. hnRNP I/PTB can antagonize the splicing repressor activity of SRp30c. RNA. 2007;13:1287–300.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Clower CV, Chatterjee D, Wang Z, Cantley LC, Vander Heiden MG, Krainer AR. The alternative splicing repressors hnRNP A1/A2 and PTB influence pyruvate kinase isoform expression and cell metabolism. Proc Natl Acad Sci. 2010;107:1894–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  61. Pandit S, Zhou Y, Shiue L, Coutinho-Mansfield G, Li H, Qiu J, et al. Genome-wide analysis reveals SR protein cooperation and competition in regulated splicing. Mol Cell. 2013;50:223–35.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  62. Bradley T, Cook ME, Blanchette M. SR proteins control a complex network of RNA-processing events. RNA. 2015;21:75–92.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  63. Tian H, Kole R. Strong RNA splicing enhancers identified by a modified method of cycled selection interact with SR protein. J Biol Chem. 2001;276:33833–9.

    Article  CAS  PubMed  Google Scholar 

  64. Wang J, Smith PJ, Krainer AR, Zhang MQ. Distribution of SR protein exonic splicing enhancer motifs in human protein-coding genes. Nucleic Acids Res. 2005;33:5053–62.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  65. Zhu J, Gong JY, Goodman OB, Cartegni L, Nanus DM, Shen R. Bombesin attenuates pre-mRNA splicing of glucocorticoid receptor by regulating the expression of serine-arginine protein p30c (SRp30c) in prostate cancer cells. Biochim Biophys Acta - Mol Cell Res. 2007;1773:1087–94.

    Article  CAS  Google Scholar 

  66. Yoshino H, Enokida H, Chiyomaru T, Tatarano S, Hidaka H, Yamasaki T, et al. Tumor suppressive microRNA-1 mediated novel apoptosis pathways through direct inhibition of splicing factor serine/arginine-rich 9 (SRSF9/SRp30c) in bladder cancer. Biochem Biophys Res Commun. 2012;417:588–93.

    Article  CAS  PubMed  Google Scholar 

  67. Yoshino H, Chiyomaru T, Enokida H, Kawakami K, Tatarano S, Nishiyama K, et al. The tumour-suppressive function of miR-1 and miR-133a targeting TAGLN2 in bladder cancer. Br J Cancer. 2011;104:808–18.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We thank Drs. John and Christine Seidman for technical assistance with RNA-Seq data collection and analysis. We thank Dr. Joe Dragavon and the BioFrontiers Advanced Light Microscopy Core facility for technical support during image collection. We thank Dr. Eunhee Chung for assistance with embryonic tissue collection.

Funding

This work was supported by NIH GM29090. KKB was supported by National Institutes of Health post-doctoral training grant 2T32HL007822-11A2.

Author information

Authors and Affiliations

Authors

Contributions

KKB and LAL conceived and designed the study. KKB performed bioinformatics analyses, luciferase assays, and stable cell line generation and accompanying analysis, analyzed data, and co-wrote the manuscript. MB was a major contributor to data analysis and co-wrote the manuscript. EKP performed mutagenesis on reporter gene constructs and performed luciferase assays. MMP assisted with cloning of reporter gene constructs. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Leslie A. Leinwand.

Ethics declarations

Ethics approval and consent to participate

All animal work was approved by the University of Colorado Boulder Institutional Animal Care and Use Committee under protocols 1002.07 and 1002.08.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1: Figure S1.

Representative masked images used for D6 myotube area calculation. (A) Myosin heavy chain, a marker of terminal differentiation, was visualized with the F59 antibody and an AlexaFluor568-linked secondary antibody. (B) DNA was visualized by DAPI staining. (C) The binary myosin mask corresponding to Panel A. (D) The binary DNA mask corresponding to Panel B.

Additional file 2: Figure S2.

miR-1 and miR-206 levels increase during C2C12 differentiation. Expression was assessed by qPCR and normalized to sno202 in growing myoblasts (D0) and days 1, 2, 4, and 6 of differentiation. Fold changes vs D0 are presented. * = p < 0.05; *** = p < 0.001

Additional file 3: Figure S3.

Non-myogenic miR-450a-1 is processed and active when ectopically expressed in myoblasts. C2C12 myoblasts were co-transfected with a reporter construct containing two tandem copies of the reverse complement of miR-450a-1 (2x450a-1) in the 3’ UTR along with expression constructs for miR-206, miR-1, or miR-450a-1. Cells were harvested 24 hours later. Firefly-normalized Renilla luciferase activity was equivalent between miR-206 and miR-1 expression but significantly down-regulated in the presence of miR-450a-1. *** = p < 0.001 vs miR-206; ### = p < 0.001 vs miR-1 N = 3 independent cultures for each.

Additional file 4: Figure S4

. Endogenous Srsf9 mRNA and protein decrease with increased miR-1/206 expression. (A) Srsf9 mRNA levels decrease in C2C12 myoblasts transfected with miR-1 or miR-206 overexpression plasmids. Expression was assessed by qPCR and 18S rRNA was the reference gene. Expression levels were normalized to control myoblast levels (indicated by the dashed line). * = p ≤ 0.05; ** = p ≤ 0.01 vs. Control myoblasts. N = 3 independent cultures for each condition. (B) Srsf9 protein levels decrease in differentiating C2C12. Srsf9 was detected by western blot with a Srsf9-specific antibody (top panel). Ponceau S staining indicates equal protein loading (bottom panel).

Additional file 5: Figure S5.

miR-1/206-resistant Srsf9-GFP expression does not change during differentiation of stable cell lines. (A) GFP or Srsf9-GFP mRNA levels in the corresponding stable C2C12 cell lines were compared to a negative control C2C12 cell line (pC-Empty) incorporating the empty expression vector pCDNA3.1(-). Expression was assessed by qPCR with GFP-specific primers and normalized to 18S rRNA levels. Fold change relative to pC-Empty is presented. N = 3 independent cultures for each. (B) Total Srsf9 (endogenous plus Srsf9-GFP) in Srsf9-GFP myoblasts is only 1.7-fold higher than endogenous Srsf9 levels in GFP control myoblasts. Expression was assessed by qPCR with Srsf9-specific primers and normalized to 18S rRNA levels. Fold change relative to GFP Control is presented. N = 3 independent cultures for each. * = p < 0.05 (C) Srsf9-GFP protein is expressed in the stable cell line. Expression in pC-Empty, GFP, and Srsf9-GFP cells was assessed by western blot with a GFP-specific antibody. (D) Srsf9-GFP mRNA levels are stable during differentiation of the Srsf9-GFP cell line. Expression was assessed by qPCR as in A. There is no statistical difference amongst the time points. (E) Srsf9-GFP protein levels are stable during differentiation of the Srsf9-GFP cell line. Expression was assessed by in-gel GFP autofluorescence (top panel). The gel was subsequently transferred to nitrocellulose which was then Ponceau S stained to reveal equal total protein loading (bottom panel).

Additional file 6: Table S1.

Primer sequences. All primer sequences are presented 5’ ➔ 3’.

Additional file 7: Table S2.

C2C12 RNA–Seq data. mRNA from C2C12 day 6 myotubes (AD1, AD2, AD3) and proliferating myoblasts (BD1, BD2, BD3) was sequenced with a paired end protocol on an Illumina HiSEQ. The table presents reads per million aligned to the mouse mm9 genome. Column A is the UCSC Gene Symbol and Columns B-E list chromosomal coordinates and coding strand designation.

Additional file 8: Table S3.

TargetScan-predicted miR-1/206 targets that decrease during C2C12 differentiation. Predicted mouse miR-1/206 targets from TargetScan v. 7.2 were crossed with all mRNAs that decreased by 1.3-fold or greater during C2C12 differentiation. This filtered the list of all 896 candidate targets to 354 that were down-regulated at the mRNA level. This is the list that was used for functional clustering gene ontology analysis with the DAVID database. Fold change from D0 to D6 is presented along with MRE scores from TargetScan.

Additional file 9: Table S4.

Functional Clustering output for DAVID analysis of 354 candidate miR-1/206 targets. Candidate miR-1/206 targets were uploaded to DAVID and the functional clustering tool was applied to the gene set. 21 GO term clusters of related gene sets had an enrichment score greater than 1.3, corresponding to an average p-value of less than 0.05 for all GO categories in a cluster. Clusters are ranked by descending Enrichment Score and include all output columns from DAVID.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Bjorkman, K.K., Buvoli, M., Pugach, E.K. et al. miR-1/206 downregulates splicing factor Srsf9 to promote C2C12 differentiation. Skeletal Muscle 9, 31 (2019). https://doi.org/10.1186/s13395-019-0211-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13395-019-0211-4

Keywords