miR-1/206 downregulates splicing factor Srsf9 to promote C2C12 differentiation
Skeletal Muscle volume 9, Article number: 31 (2019)
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.
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.
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.
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.
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 . 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 . 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 , 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 . 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 .
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 .
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.
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.
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 . 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 . 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.
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.
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.
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 . 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 , Pax3 , and Sox9 , 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 , 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 . 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.
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].
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 . 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  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.
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 . 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 . 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 .
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 ). 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 , 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 . 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 . 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 . However, a computational survey of SELEX-derived consensus sequences, which are short and degenerate, revealed frequent overrepresentation in RNAs .
miR-1 and miR-206 have both been proposed as tumor suppressors. Interestingly, Srsf9 depletion reduces viability of prostate cancer cells . Moreover, in two human bladder cancer cell lines, Srsf9 downregulation by miR-1 overexpression resulted in significant reduction in cell proliferation, migration, and invasion . 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.
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.
Database for Annotation, Visualization, and Integrated Discovery
Green fluorescent protein
miRNA response element
- SR protein:
Nakka K, Ghigna C, Gabellini D, Dilworth FJ. Diversification of the muscle proteome through alternative splicing. Skelet Muscle. 2018;8:8.
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.
Llorian M, Smith CW. Decoding muscle alternative splicing. Curr Opin Genet Dev. 2011;21:380–7.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Long JC, Caceres JF. The SR protein family of splicing factors: master regulators of gene expression. Biochem J. 2009;417:15–27.
Zhou Z, Fu X-D. Regulation of splicing by SR proteins and SR protein-specific kinases. Chromosoma. 2013;122:191–207.
Jumaa H, Wei G, Nielsen PJ. Blastocyst formation is blocked in mouse embryos lacking the splicing factor SRp20. Curr Biol. 1999;9:899–902.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
TargetScan. www.targetscan.org/vert_72/. Accessed 20 Mar 2019.
DAVID Bioinformatics Resources 6.8 https://david.ncifcrf.gov/home.jsp. Accessed 22 Mar 2019.
Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4:44–57.
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.
Ä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.
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.
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.
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.
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.
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.
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.
pSILAC. https://psilac.mdc-berlin.de/proteingroups/906/3/. Accessed 23 Apr 2019.
Twyffels L, Gueydan C, Kruys V. Shuttling SR proteins: more than splicing factors. FEBS J. 2011;278:3246–55.
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.
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.
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.
Horak M, Novak J, Bienertova-Vasku J. Muscle-specific microRNAs in skeletal muscle development. Dev Biol. 2016;410:1–13.
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.
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.
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.
Dey BK, Gagan J, Dutta A. miR-206 and -486 induce myoblast differentiation by downregulating Pax7. Mol Cell Biol. 2011;31:203–14.
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.
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.
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.
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.
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.
Simard MJ, Chabot B. SRp30c is a repressor of 3’ splice site utilization. Mol Cell Biol. 2002;22:4001–10.
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.
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.
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.
Bradley T, Cook ME, Blanchette M. SR proteins control a complex network of RNA-processing events. RNA. 2015;21:75–92.
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.
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.
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.
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.
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.
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.
This work was supported by NIH GM29090. KKB was supported by National Institutes of Health post-doctoral training grant 2T32HL007822-11A2.
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.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
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.
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
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.
. 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).
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).
Primer sequences. All primer sequences are presented 5’ ➔ 3’.
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.
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.
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.
About this article
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