- Open Access
A role for RNA post-transcriptional regulation in satellite cell activation
Skeletal Muscle volume 2, Article number: 21 (2012)
Satellite cells are resident skeletal muscle stem cells responsible for muscle maintenance and repair. In resting muscle, satellite cells are maintained in a quiescent state. Satellite cell activation induces the myogenic commitment factor, MyoD, and cell cycle entry to facilitate transition to a population of proliferating myoblasts that eventually exit the cycle and regenerate muscle tissue. The molecular mechanism involved in the transition of a quiescent satellite cell to a transit-amplifying myoblast is poorly understood.
Satellite cells isolated by FACS from uninjured skeletal muscle and 12 h post-muscle injury from wild type and Syndecan-4 null mice were probed using Affymetrix 430v2 gene chips and analyzed by Spotfiretm and Ingenuity Pathway analysis to identify gene expression changes and networks associated with satellite cell activation, respectively. Additional analyses of target genes identify miRNAs exhibiting dynamic changes in expression during satellite cell activation. The function of the miRNAs was assessed using miRIDIAN hairpin inhibitors.
An unbiased gene expression screen identified over 4,000 genes differentially expressed in satellite cells in vivo within 12 h following muscle damage and more than 50% of these decrease dramatically. RNA binding proteins and genes involved in post-transcriptional regulation were significantly over-represented whereas splicing factors were preferentially downregulated and mRNA stability genes preferentially upregulated. Furthermore, six computationally identified miRNAs demonstrated novel expression through muscle regeneration and in satellite cells. Three of the six miRNAs were found to regulate satellite cell fate.
The quiescent satellite cell is actively maintained in a state poised to activate in response to external signals. Satellite cell activation appears to be regulated by post-transcriptional gene regulation.
Skeletal muscle is terminally differentiated and thus, requires a population of resident adult stem cells, satellite cells, for maintenance and repair [1–3]. Satellite cells are typically mitotically quiescent in resting muscle and activate to prepare for cell cycle entry by HGF [4, 5], nitric oxide , and TNFα , upon a muscle injury. Intracellular p38α/β MAPK and downstream signaling is stimulated upon satellite cell activation, permitting MyoD induction (Troy et al.)a, S-phase entry [8, 9], and subsequent proliferation. A subset of satellite cells self-renew to maintain the satellite cell pool (Troy et al.)a[10, 11] and generate a rapidly proliferating transit-amplifying myoblast population (Troy et al.)a.
The transition from a quiescent satellite cell to a proliferating, transit amplifying myoblast was thought to require extensive transcriptional induction as quiescent satellite cells have a low ratio of cytoplasmic volume to nuclear volume, few cellular organelles, tightly packed heterochromatin, and are believed to be metabolically inactive [12, 13]. However, recent evidence suggests that satellite cell quiescence is ‘active’ and satellite cells are poised to react to external stimuli after muscle damage . Moreover, quiescent fibroblasts exhibit high metabolic activity  in agreement with a quiescent state that is far from ‘quiet’. Interestingly, a growing pool of data demonstrates that cell fate determination is reliant on post-transcriptional gene regulation [16–20] and may provide mechanisms to maintain quiescent satellite cells in a ready state.
One such RNA post-transcriptional mechanism, microRNA-mediated gene silencing, regulates skeletal muscle specification and myogenic differentiation [21–23]. MicroRNAs (miRNA) are a class of small non-coding RNAs that bind to target mRNA in a sequence specific manner to mediate gene silencing [24–27] and can target and silence protein expression from tens to hundreds of mRNAs [26, 27]. Furthermore, miRNAs modulate stem cell fate decisions [28–31] and may have similar functions in satellite cells. Recent studies identify miR-489 and miR-206 expression in quiescent satellite cells [32, 33], however, it is likely that many uncharacterized miRNAs play roles in the transition of a quiescent satellite cell to transit-amplifying myoblast.
To understand the mechanisms involved in satellite cell activation, we previously screened a number of candidate genes for changes in expression from freshly isolated satellite cells and from satellite cells isolated at either 12 h post-muscle injury or 48 h post-muscle injury to represent quiescent, activated, and proliferating satellite cells, respectively. Although unbiased gene expression screens have been performed on satellite cells, these studies have either compared freshly isolated satellite cells to satellite cells expanded in culture [14, 34] or to satellite cells in diseased skeletal muscle . Neither of these studies directly compared satellite cells prior to and following induced muscle injury in vivo and thus, the reported gene expression changes specific to cell culture or specific to diseased muscle may not reliably identify gene expression changes associated with satellite cell activation in vivo. Here, we report global gene expression profiles and candidate miRNAs associated with quiescent and activated satellite cells as well as identify a novel function for miR-16, miR-106b, and miR-124 in satellite cell fate determination. From these analyses, we posit that satellite cell activation is primarily regulated by post-transcriptional gene regulation as opposed to transcriptional induction.
All animal procedures were performed according to protocol number 1012.01 approved by Institutional Animal Care and Use Committee at the University of Colorado at Boulder. Mice were housed in a pathogen-free environment at the University of Colorado at Boulder. All mice sacrificed were female and between 3 and 6 months of age. Wild type mice were C57Bl/6xDBA2 (B6D2F1/J, Jackson Labs) and syndecan-4−/− mice carry homozygous deletion of syndecan-4 in the C57Bl/6 background .
Fluorescence-activated cell sorting of satellite cells
The tibialis anterior muscles of 3-month-old female B6D2F1/J or syndecan-4−/− mice were injured by injection with 50 μL 1.2% BaCl2 in saline prior to harvest or harvested from uninjured hind limbs. The tibialis anterior muscles were dissected from the hind limb, minced, and digested in 400 U/mL collagenase in Ham’s F-12C at 37°C for 1 h, vortexing frequently. Collagenase was inactivated by the addition of horse serum and debris was removed by sequential straining through 70 μm and 40 μm cell strainers (BD Falcon). Cells were gently centrifuged and the cell pellets were incubated at 4°C with 1:100 rabbit anti-syndecan-3 antibody in Ham’s F-12C with 15% horse serum followed by an incubation on ice with Cy5 conjugated anti-rabbit-IgG (Molecular Probes). Satellite cells were sorted based on syndecan-3 immunoreactivity on a MoFlo Legacy cell sorter (Dako Cytomation) directly into RNA lysis buffer (PicoPure RNA Isolation kit, Arcturus).
Myofiber explant culture and immunostaining
All hind limb muscles were dissected, connective tissue removed, and individual muscle groups isolated followed by digestion in 400 U/mL collagenase in Ham’s F-12C at 37°C. Single myofibers were isolated and grown in Ham’s F-12C supplemented with 15% horse serum and 0.5 nM FGF-2 prior to fixation in 4% PFA. Fibers were blocked in 10% normal goat serum in phosphate buffered saline followed by antibody staining. Primary antibodies were rabbit anti-cmet (Santa Cruz) at 1:100, mouse anti-MyoD (Novocastra) at 1:10, mouse anti-Pax7 at 1:5 (Developmental Studies Hybridoma Bank), and rabbit anti-MyoD C-20 at 1:500 (Santa Cruz Biotechnology). Secondary antibodies were Alexa-488 conjugated anti-mouse IgG, Alexa-594 conjugated anti-rabbit IgG, Alexa-555 conjugated anti-mouse IgG, and Alexa-647 conjugated anti-rabbit IgG (Molecular Probes). All images taken on a Nikon Eclipse E800 microscope with a Nikon 40x/0.75 differential interference contrast M lens and analyzed with Slidebook (Intelligent Imaging Innovations, Inc.).
RNA was isolated from satellite cells using the PicoPure RNA Isolation kit (Arcturus) followed by two rounds of linear T7-based amplification (RiboAmp HA kit: Arcturus). The RNA equivalent of 5,000 cells was hybridized to Affymetrix mouse 430v2 GeneGhips (MOE430v2) according to manufacturer’s instructions. GeneChips were scanned at the University of Colorado at Boulder on an Affymetrix GeneChip Scanner 3000 and spot intensities were recovered in the GeneChip Operating System (Affymetrix).
Microarray data processing and analysis
All analysis was performed using Spotfire™ DecisionSite 2 for Microarray Analysis. The raw CEL data files were normalized using GC Robust Multi-array Analysis (GCRMA). The raw CEL data files, microarray metadata, and GCRMA normalized expression values were deposited in GEO datasets (GSE38870). One wild type freshly isolated satellite cell replicate consistently clustered with the wild type satellite cells 12 h post-injury replicates (via hierarchical, Self-Organizing Map (SOM), k-means) indicating myogenic commitment and was removed from our analysis. The hierarchical cluster and associated dendrogram were generated using the log2-value for relative probe intensity using the Unweighted Pair Group Method with Arithmetic Mean (UPGMA) with Euclidean distance as the similarity measure. The significance between genotypes and time points was determined using the multifactor analysis of variance (ANOVA) with a false discovery rate (FDR) ≤ 0.05 and Bonferroni adjustment. Fold change was calculated as 2abs(difference), where difference is the log2 difference between samples compared. Venn diagrams were generated in Spotfire™ using the list comparison function.
Gene ontology and biological pathway analysis
Unique gene identifiers (gene symbol, entrez gene ID, or Affymetrix probeset ID) were uploaded to the Database for Annotation, Visualization and Integrated Discovery (DAVID http://david.abcc.ncifcrf.gov/), FunNet (http://www.funnet.info/), and ProfCom (http://webclu.bio.wzw.tum.de/profcom/). The mouse genome reference dataset for each algorithm was used as background. No further settings were required for DAVID. Analysis setting for FunNet were conventional functional analysis with the specificity enrichment computation for GO and false discovery rate ≤ 5%. Analysis settings for ProfCom were up degree 1 and exclude. The default settings for each algorithm were used to identify enriched GO terms. Affymetrix probeset IDs and the log2 difference between wild type freshly isolated satellite cells and satellite cells isolated 12 h post-injury were uploaded and analyzed using IPA v9.0 (Ingenuity® Systems http://www.ingenuity.com). Analysis settings were to consider both all direct and indirect molecules and/or relationships using the Mouse Genome 430 2.0 Array as a reference dataset.
Computational prediction of miRNAs
The Srivastava lab algorithm assessed mRNA sequences in human, mouse, and rat for miRNA seed matches and required base-paring of miRNA nucleotides 2 to 7 with binding energy ≤ −14 kcal/mol and flanking energy ≥ −7 kcal/mol. Secondary structure was used to eliminate false-positives by removing those seed matches with secondary elements that stabilize mRNA. Priority 1 calls had a destabilizing mRNA element while Priority 2 calls did not contain a destabilizing element in at least one species. GeneAct used the miRanda algorithm to identify miRNA target sites across three mammalian species. False-positives were eliminated with the differential binding site search against genes that were constitutively expressed in satellite cells.
mRNA was extracted from satellite cells using the RNeasy Kit according to the manufacturer’s protocol (Qiagen). miRNA was extracted using both the RNAqueous-micro kit the mirVana miRNA isolation kit according to the manufacturer’s protocols (Ambion) with the following modifications. For satellite cells, the mirVana manufacturer’s protocol was followed using volumes and columns for the RNAqueous-micro kit. For MM14 cells, the mirVana manufacturer’s protocol was followed. RNA concentration was determined using a NanoDrop 2000 spectrophotometer (Thermo Scientific).
The Superscript III First Strand cDNA Synthesis kit was used to generate cDNA from mRNA according to manufacturer’s instructions (Invitrogen). The Ncode miRNA qRT-PCR system (Invitrogen) was used to generate cDNA from miRNA according to manufacturer’s instructions. Briefly, RNA was poly-adenylated with Poly-A polymerase followed by cDNA transcription with Superscript III reverse transcriptase using a primer similar to Oligo-dT with a unique 5’ end. Quantitative RT-PCR was performed using SYBR-Green (Applied Biosystems) or SYBR-GreenER (Invitrogen) on either an ABI 7500 Fast or ABI 7900 Real-Time PCR machine (Applied Biosystems). Primer sequences are listed in Table 1.
RNase protection assay
Candidate miRNAs were screened using the mirVana miRNA detection kit according to the manufacturer’s protocol (Ambion). All probes were radio-labeled with 32P-UTP using the mirVana miRNA probe construction kit (Ambion) according to manufacturer’s protocol.
Isolation of quiescent satellite cells
Quiescent satellite cells were isolated following IP injection with 75 mg/kg of SB203580 (Alexis Corporation) and kept in 25 μM SB203580 through the isolation as described (Hausburg et al., Submitted). All satellite cells were isolated as follows. Hind limb muscles of 3- to 6-month-old female B6D2F1/J mice were dissected and digested in 400 U/mL collagenase in Ham’s F-12C at 37°C for 1 h with periodic vortexing. The collagenase was inactivated with horse serum and debris was removed with sequential straining through 70 μm and 40 μm cell strainers (BD Falcon). Satellite cells were either isolated at the interface of a 40%/70% percoll gradient (GE Healthcare) or plated in Ham’s F-12C supplemented with 15% horse serum and 0.5nM FGF-2 for various times before RNA isolation.
Myofiber explant cultures were transfected using Lipofectamine 2000 (Invitrogen) according to manufacturer’s protocol with a 2.5:1 ratio of lipofectamine:nucleic acid. 200nM miRIDIAN hairpin inhibitors (Dharmacon) against miR-16, miR-93, miR-106b, and miR-124 were co-transfected with pEGFP-C1-H2B. Scrambled control was the miRIDIAN hairpin inhibitor negative control 1 (Dharmacon).
Identification of gene expression changes associated with satellite cell activation in vivo
An unbiased global gene expression analysis using Affymetrix GeneChips was performed to identify changes occurring during the transition of satellite cells from quiescence to a population of proliferating myoblasts in vivo. To accomplish this, we identified changes in gene expression profiles between freshly isolated satellite cells and satellite cells isolated 12 h or 48 h following BaCl2-induced muscle injury. The time points chosen correspond to activated satellite cells which do not express MyoD protein (freshly isolated), committed myoblasts marked by MyoD expression (12 h post-muscle injury), and proliferating myoblasts (48 h post-muscle injury; Figure 1A) [36, 37]. Syndecan-4 null satellite cells fail to activate, express MyoD, or enter the cell cycle appropriately within 48 h post-injury, and are incapable of skeletal muscle repair (Figure 1B) [36, 37]. Therefore, we eliminated the genes whose expression changes following a muscle injury in Sdc4−/− satellite cells from our analyses as these genes were unlikely to be involved in satellite cell activation. Wild type and Sdc4−/− satellite cells were isolated by fluorescence activated cell sorting (FACS) using anti-syndecan-3 antibodies as a marker for quiescent and proliferating satellite cells  from uninjured tibialis anterior (TA) muscle (Figure 1C) and TA muscles 12 h post-injury (Figure 1D).
RNA was extracted from the isolated cells and processed for hybridization to Affymetrix v2.0 mouse GeneChips. The GeneChip data (Additional file 1, Additional file 2, Additional file 3, Additional file 4, Additional file 5, Additional file 6) were analyzed with Spotfire™ DecisionSite for Microarry Analysis software and an initial hierarchical cluster dendrogram generated. Visualization of the relationships between gene expression profiles show Sdc4−/− samples, regardless of time post-injury, cluster with freshly isolated wild type satellite cells (Figure 1E) supporting our observations that Sdc4−/− satellite cells do not activate appropriately within 48 h following an induced muscle injury (Figure 1A, B). Within the wild type dataset, we observed that committed myoblasts isolated 12 h post-injury exhibited the most divergent gene expression profiles, suggesting that these committed satellite cells differ substantially from either quiescent satellite cells or proliferating myoblasts (Figure 1E). To focus on genes that may be involved in satellite cell activation, we chose to further compare gene expression changes occurring within the first 12 h post-muscle injury.
Over 4,000 genes are specifically regulated during satellite cell activation
A comparative analysis of gene expression profiles from wild type and Sdc4−/− satellite cells within the first 12 h following satellite cell activation identified a cohort of genes unique to satellite cell activation. In wild type satellite cells, 5,162 genes change significantly between satellite cells isolated from uninjured muscle and those isolated 12 h post-injury as defined by a ≥ 2-fold change with an ANOVA P ≤ 0.01 (Figure 2A). In contrast, 2,236 genes similarly changed expression in Sdc4−/− satellite cells isolated from uninjured TA muscle and TA muscles 12 h post-injury (Figure 2A). Eighty percent (4,093) of the genes differentially expressed in WT satellite cells do not significantly change in Sdc4−/− satellite cells as identified by Venn analysis (Figure 2A, B; Additional file 7). We reasoned that the metabolic changes occurring during satellite cell activation as well as the induction of the transcription factor MyoD and cell cycle entry would result in a large cohort of induced genes. Surprisingly, more than half (56%) of the genes differentially expressed in satellite cells by 12 h post-injury decrease in relative expression (Figure 2B). Moreover, the magnitude of change for genes that decrease is on average three-fold greater than the magnitude of change for genes that increase following muscle injury. Less than 10% of genes whose expression is increased change more than four-fold (22), while 70% of downregulated genes change more than four-fold (22) and 3% decrease more than 64-fold (26) (Figure 2C). These observations suggest that quiescent satellite cells express a cohort of genes that maintains and regulates the quiescent state, are likely critical for interaction with the satellite cell niche, and are necessary for interpreting signals for exit from quiescence. Furthermore, our results support the idea that satellite cell quiescence is actively maintained, consistent with a prior report examining freshly isolated satellite cells and satellite cells isolated from dystrophic muscle .
To further test the idea that satellite cell quiescence is actively maintained, we analyzed genes in the cohort that significantly change 12 h post-injury involved in cell cycle progression. We would expect cell cycle progression genes to be induced during activation and found that genes modulating the G1/S phase transition are among those that increase in relative expression (Figure 2D, green). In contrast, cell cycle inhibitors decrease in relative expression (Figure 2D, red) as expected for the transition of satellite cells from mitotic quiescence to an activated state in preparation for cell cycle entry. Moreover, these genes do not change expression significantly in Sdc4−/− satellite cells 12 h post-injury (Figure 2A, B; Additional file 7) consistent with their impaired cell cycle activation and MyoD induction.
Genes involved in RNA post-transcriptional regulation are significantly enriched during satellite cell activation
Gene ontology (GO) classifications were used to aid in identifying potential mechanisms regulating satellite cell activation. The Database for Annotation, Visualization and Integrated Discovery (DAVID) was used to identify enrichment of general molecular function categories [39, 40] during activation of satellite cells. The GO category of Molecular Function: Binding is the most significantly over-represented GO category during satellite cell activation (P value = 7.03 × 10-43 compared to the mouse genome), where a 7% increase in the total percentage of genes classified as binding occurs in satellite cells within the first 12 h post-muscle injury (Figure 3A). Further refinement of ontological categories reveals that GO terms unfolded protein binding, actin binding, and RNA binding were enriched an average of 1.5-fold over three independent gene annotation algorithms: FunNet [41, 42], ProfCom , and DAVID [39, 40] (Figure 3B, Table 2). Thus, in the first 12 h post-muscle injury, major changes occur in genes involved in RNA binding, the unfolded protein response, and in actin binding. The changes in RNA binding proteins and unfolded protein response may be involved in the down-regulation of genes necessary to maintain a quiescent satellite cell, while changes in actin binding are likely to reflect changes in satellite cell motility  that accompany repair of skeletal muscle tissue.
A fourth independent analysis focused on biological networks (IPA-Ingenuity Pathway Analysis http://www.ingenuity.com) ranked RNA Post-Transcriptional Modification in the top biological network (Table 3; Table 4). Thus, from four independent methods of gene expression analysis, a much higher proportion of genes involved in post-transcriptional RNA regulation change expression in the transition from a quiescent satellite cell to a committed myoblast, suggesting a role for post-transcriptional regulation of RNA in this transition. Therefore, we further analyzed genes involved in RNA post-transcriptional modification to further characterize individual genes and to develop hypotheses regarding the function of these genes in the transition of satellite cells from mitotic quiescence to cell cycle entry.
Splicing factors are preferentially downregulated during satellite cell activation
The top ranking biological network identified by IPA has associated biological functions of RNA Post-Transcriptional Modification, Developmental Disorder, and Genetic Disorder (Table 3), where RNA Post-Transcriptional Modification includes mRNA decay and stabilization, mRNA splicing, and miRNA-mediated gene silencing. When we examined genes within this network, nearly half of the genes whose expression changes during satellite cell activation are involved in RNA processing and splicing (Figure 4). Moreover, many genes in this interaction network are implicated in either muscle function, muscle disease, or cell fate decisions (Table 4). Consistent with our prior observations that downregulated genes exhibit greater fold change in relative expression occurring during the first 12 h post-muscle injury, the majority of RNA processing and splicing genes identified in the interaction network (87% and 88%, respectively) decrease post-muscle injury (Figure 4). Furthermore, of the 154 RNA binding proteins identified by GO analysis (Figure 3B, Additional file 8), 69% decrease in relative expression (Figure 5A; Additional file 8). These data show that RNA binding proteins are highly over-represented in quiescent satellite cells and suggest that regulation of RNA plays an important role in maintaining the quiescent state and in the transition to a cycling myoblast.
An examination of the top quartile of RNA binding proteins with the most significant changes in expression reveal that the majority (45%) of upregulated genes are positive mRNA regulators (Figure 5B), while downregulated genes include a similar distribution of all RNA protein functions with splicing factors exhibiting a slight majority (37%: Figure 5C). Furthermore, 24% of the previously identified RNA processing of splicing factors (Figure 4) are present in this top quartile. To confirm the microarray expression data, we performed qPCR validation finding that approximately 80% of the tested RNA binding proteins and both tested cell cycle regulators have consistent gene expression profiles (Table 5). Thus, genes involved in mRNA regulation may play diverse roles in satellite cell activation and promote the conversion of a quiescent satellite cell to a proliferating myoblast.
Dynamic regulation of miRNAs during muscle regeneration
MicroRNA-mediated gene silencing regulates alternative splicing [52, 53] as well as mRNA stability factors [54, 55]. Moreover, miRNAs regulate stem cell fate determination [28, 30, 31, 56] suggesting a potential role for miRNAs in the transition of quiescent satellite cells to proliferating myoblasts. The low levels of RNA present in quiescent satellite cells combined with infrequent satellite cell abundance in uninjured muscle prevented successful unbiased screen for miRNAs. Therefore, we assessed whether genes involved in miRNA biogenesis and gene silencing including argonautes1-4 (Eif2c1-4), Dgcr8, and Dicer, as well as other genes associated with miRNA function, are expressed in quiescent satellite cell (Additional file 1, Additional file 2, Additional file 3, Additional file 4, Additional file 5, Additional file 6). Although present, Argonautes1-4 (Eif2c1-4), Dgcr8, and Dicer are not in the cohort of differentially expressed genes (Additional file 7), suggesting that major changes in miRNA processing do not occur during satellite cell activation.
To identify potential miRNAs involved in the transition of quiescent satellite cells to proliferating myoblasts, we applied miRNA target prediction algorithms to identify putative miRNAs regulating genes whose expression changes rapidly during the first 48 h post-muscle injury. Initially, we established a minimum gene expression value in freshly isolated satellite cells to reduce the cohort to 641 genes (ANOVA P ≤ 0.01) that either increased or decreased by ≥ two-fold during the first 48 h following muscle injury. As gene expression levels did not change for 47 of these genes in Sdc4−/− satellite cells (ANOVA P ≥ 0.9), we chose these 47 genes to pursue as potential miRNA targets involved in satellite cell activation (Table 6). Candidate miRNAs were then computationally identified using two independent algorithms, one developed by the group of Deepak Srivastava (unpublished) and GeneAct (http://promoter.colorado.edu/geneact/)  (Figure 6A). The union of both algorithms identified 12 miRNA candidates with six, miR-16, miR-93, miR-106b, miR-107, miR-124, and miR-200b, being detected in cultured primary satellite cells or proliferating satellite cell derived MM14 cells by ribonuclease protection assay (Table 7). All six miRNAs present in primary myoblasts and MM14 cells were detectable in uninjured tibialis anterior muscle (Figure 6B).
We observed dramatic regulation of these six miRNAs following a muscle injury when compared to control, uninjured tibialis anterior muscle. Four of the identified miRNAs (miR-93, miR-107, miR-124, and miR-200b) changed expression levels by more than two-fold during the first 5 days following induced muscle injury (Figure 6B-F; Table 8). In the uninjured TA muscle, miR-200b decreased three-fold by 12 h post-injury, while miR-93 and miR-124 increased significantly 12 h post-muscle injury (Figure 6B, C, F; Table 8). Within 48 h post-injury, the relative levels of miR-93, miR-107, and miR-124 had decreased levels well below those present in uninjured muscle and remained low at 5 days post-injury (Figure 6C-E). In contrast, miR-106b remained elevated following injury while miR-16 trended to slightly lower expression (Figure 6B-E). The rapid changes in miRNA relative expression and their presence in skeletal muscle suggest that these miRNAs may play important roles in the regeneration of skeletal muscle and validates our approach to identify such miRNAs.
Relative expression of miRNAs in satellite cells following muscle injury
We asked whether the six miRNAs that change expression during muscle regeneration are present in satellite cells, muscle tissue, or both. The relative expression levels for each miRNA in the uninjured tibialis anterior muscle was normalized to 1 and the relative levels in satellite cells isolated from uninjured muscle and in proliferating satellite cells isolated 48 h post-injury examined. Surprisingly, all six miRNAs were expressed at low to undetectable levels in freshly isolated satellite cells (Figure 6G). In proliferating satellite cells, all but miR-106b were expressed at levels substantially lower than that found in the tibialis anterior muscle (Figure 6H). Although MyoD protein is not detectable in freshly isolated satellite cells and they have not yet entered S-phase, freshly isolated satellite cells are not quiescent since the p38α/β MAPK is activated . Therefore, to identify miRNAs present in quiescent satellite cells, mice were injected with SB203580, a p38α/β MAPK inhibitor, 1.5 h prior to satellite cell isolation (Hausburg et al., Submitted), and relative miRNA levels examined. Remarkably, we found that miR-124 was expressed 35-fold higher in quiescent satellite cells than in freshly isolated or proliferating satellite cells (Figure 6I, J; Table 8). Moreover, miR-124 was expressed at levels 10-fold greater in quiescent satellite cells than in uninjured skeletal muscle (compare Figure 6I and 6B, Figure 6J; Table 8), suggesting that the primary source of miR-124 in uninjured muscle is the satellite cell population. In contrast to miR-124, miR-16 and miR-93 are present at low to undetectable levels in quiescent satellite cells and are induced in freshly isolated satellite cells (Figure 6G, I). The expression level of miR-16 is maintained in proliferating cells, while miR-93 declines and miR-106b is dramatically induced in proliferating satellite cells as compared to freshly isolated and quiescent satellite cells (Figure 6G-I).
miR-16, miR-106b, and miR-124 regulate satellite cell fate
The changes in relative levels of miR-16, miR-93, miR-106b, and miR-124 in satellite cells following a muscle injury suggests that these four miRNAs may play a role in the transition from a quiescent satellite cell to a proliferating myoblast. To test this idea, inhibitors for each miRNA were transfected into myofiber-associated satellite cells immediately following isolation and the cultures fixed and assayed at 3 days post-isolation and 5 days post-isolation. The total number of Pax7+ cells transfected with the scrambled RNA control inhibitor declined two-fold between 3 and 5 days in culture, indicative of differentiation (Figure 7A, B). Between 3 and 5 days in culture, the Pax7+/MyoD + decreased three-fold accompanied by the appearance of the Pax7+/MyoD- ‘reserve’ population (Figure 7C-F). In contrast, inhibition of miR-124 increased the percentage of Pax7+/MyoD- ‘reserve’ cells at 3 and 5 days of culture, as did inhibition of miR-106b (Figure 7A-F). Inhibition of miR-16 elevated the total number of Pax7+ cells at 5 days of culture and inhibition of miR-93 did not have any detectable effect (Figure 7-F). We further examined the role of miRNAs in satellite cell activation using Ingenuity® System’s IPA and identified PTEN signaling and Cell Cycle Regulation by BTG Family Proteins as the top canonical pathway regulated by miR-16, miR-93, miR-106b, and miR-124 in the transition of a quiescent satellite cell to a proliferating myoblasts (Figure 8). Many predicted targets of miR-16 and the miR-93/106b family inhibit cell cycle progression and cell growth. These targets are downregulated during satellite cell activation, consistent with increased expression of miR-16, miR-93, and miR-106b in proliferating satellite cells as compared to quiescent satellite cells (Figure 8). Conversely, predicted target genes of miR-124 promote cell cycle progression and are upregulated during satellite cell activation when miR-124 is downregulated (Figure 8). These data demonstrate that a number of miRNAs regulate satellite cell fate following a muscle injury and support the idea that post-transcriptional regulation of RNA plays a critical role in satellite cell activation and maintenance of satellite cell quiescence.
The low cytoplasmic to nuclear ratio, low organelle number, and mitotic quiescence of resident satellite cells [12, 13] lead to the speculation that metabolic activity in these cells is low. Indeed, quiescent satellite cells with high levels of Pax7 express reduced levels of mitochondrial genes . Moreover, the delay to the first cell cycle division (Troy et al.)a coupled with the dramatic increase in cell size and mobility [36, 44] suggests that satellite cell activation and cell cycle entry would require transcriptional induction of a large cohort of genes similar to that observed in serum stimulated fibroblasts . However, we and others have postulated that quiescent satellite cells are poised for activation awaiting a critical signaling event [8, 14]. The data presented here further support this hypothesis and provide the first direct comparison of quiescent satellite cells with activated satellite cells and proliferating myoblasts derived from uninjured and injured skeletal muscle, respectively. The prior analyses performed compared freshly isolated satellite cells with satellite cells isolated from dystrophic mice  or cultured cells [14, 34] and are expected to identify gene expression changes associated with a diseased environment or a culture environment, respectively. Since the most significant reductions in gene expression occur within the first 12 h post-muscle injury, the metabolic and signaling events in a quiescent satellite cell are thus predicted to be divergent from those of a proliferating myoblast. Moreover, these data suggest that comparisons of freshly isolated satellite cells with proliferating myoblasts may not identify critical regulatory mechanisms involved in satellite cell activation [14, 34].
Here, we used computational methods to initially identify that RNA post-transcriptional mechanisms are likely to maintain the quiescent satellite cell phenotype and to promote the conversion of the quiescent satellite cell to a transit-amplifying myoblast. Recent studies indicate that RNA post-transcriptional mechanisms, specifically alternative splicing, mRNA stability, and miRNA-mediated gene silencing, regulate stem cell pluripotency and progression through differentiation [17–20]. As satellite cells are an adult stem cell population, similar mechanisms may mediate the transition from quiescence to a population of proliferating myoblasts. We found that splicing factors may play roles in regulating the transition of satellite cells from quiescence to proliferating myoblasts. Consistent with these observations, the relative molar ratios of splicing factors guide alternative splicing  and these factors often function combinatorially to direct the expression of different mRNA splice variants . Together with published data, our observations of dynamic splicing factor expression in satellite cells following muscle injury suggests that unique cohorts of mRNA species regulate the conversion of quiescent adult stem cells to the committed proliferating myoblast.
In addition to splicing factors, we found that mRNA binding proteins regulating mRNA stabilization and mRNA decay are preferentially upregulated in satellite cells following muscle injury. Indeed, the AU-rich element (ARE) binding protein HuR (Elavl1) is reported to stabilize MyoD and myogenin mRNA in skeletal muscle cell lines derived from satellite cells [62, 63] potentially participating in satellite cell activation and commitment to myogenesis. Elavl1 is one of the most significantly upregulated RNA binding genes in satellite cells within 12 h of muscle injury (See Figure 5C; Additional file 8) and corresponds to MyoD expression. Interestingly, the ARE-binding proteins, tristetraprolin (Zfp36) and family members Zfp36l1 and Zfp36l2, decrease dramatically during the conversion from quiescence to activated satellite cells within the first 12 h post-muscle injury (See Figure 5D; Additional file 8). However, tristetraprolin and HuR have opposing functions as they counter-regulate expression of the same mRNAs [64, 65] and may act as an agonist–antagonist pair for many genes that promote commitment to myogenesis. In agreement with these data, we have demonstrated that the Zfp36 family directly targets MyoD mRNA and functions to regulate satellite cell fate during satellite cell activation and self-renewal (Hausburg et al., Submitted).
A recent report identifies miR-489 as an important miRNA maintaining satellite cell quiescence  suggesting that miRNA-mediated gene silencing functions in the transition of a quiescent satellite cell to a proliferating myoblast. We identified a cohort of genes that significantly change expression in satellite cells within the first 48 h following muscle injury to computationally predict cognate miRNAs that may regulate these targets with two independent prediction algorithms. Six miRNAs not previously reported in skeletal muscle were selected for further analysis and all six were observed to be dynamically regulated in relative levels during induced muscle injury. Moreover, four of the six miRNAs were expressed in satellite cells (miR-16, miR-93, miR-106b, and miR-124), while two were likely present only in differentiated muscle (miR-107 and miR-200b). Comparing relative levels in muscle tissue and satellite cells revealed that miR-124 is likely only expressed in satellite cells, while miR-16, miR-93, and miR-106b are most likely expressed in satellite cells and in differentiated muscle fibers. Pathways predicted to be targeted by these miRNAs include cell cycle progression as well as PTEN signaling, which is involved in stem cell self-renewal  and muscle regeneration .
To directly test whether these four miRNAs regulate satellite cell behavior, we transfected inhibitors for each miRNA into myofiber-associated satellite cells shortly after isolation and examined the effects on satellite cell fate at 3 and 5 days post-myofiber isolation. Inhibition of two miRNAs, miR-106b and miR-124, increased the relative number of progenitor or ‘reserve’ satellite cells (Pax7+/MyoD-) relative to a control, suggesting that these miRNAs participate in the regulation of satellite cell fate and satellite cell self-renewal. In contrast, miR-16 enhanced the relative numbers of Pax7+ cells but did not appear to alter the percentage of Pax7+/MyoD + myoblasts or Pax7+/MyoD- ‘reserve’ cells relative to a scrambled control. Of the four miRNAs tested, the loss of miR-93 elicited no detectable changes in the numbers of ‘reserve’ Pax7+/MyoD- satellite cells and Pax7+/MyoD + satellite cells suggesting functions in cellular processes other than cell fate determination.
We believe that RNA post-transcriptional regulation plays a critical role in the transition of a quiescent satellite cell to a transit-amplifying myoblast. At each time point we examined, including quiescent satellite cells (freshly isolated in the presence of a p38α/β MAPK inhibitor), activated satellite cells (12 h post-muscle injury), and proliferating myoblasts (48 h post-muscle injury), we found extensive changes in genes involved in post-transcriptional RNA regulation, including mRNA splicing, mRNA stability, and miRNA-mediated gene silencing. We conclude that satellite cell quiescence is actively maintained via combinatorial contributions primarily mediated through post transcriptional mRNA regulation and identified four miRNAs that likely play a role in the conversion of quiescent satellite cells to proliferating myoblasts.
Following submission of our manuscript, the following was published demonstrating post-transcriptional regulation of myf-5 during satellite cells activation. Crist CG, Montarras D, Buckingham M (2012) Muscle Satellite Cells Are Primed for Myogenesis but Maintain Quiescence with Sequestration of Myf5 mRNA Targeted by microRNA-31 in mRNP Granules. Cell Stem Cell 11: 118-126.
Lepper C, Partridge TA, Fan C-M: An absolute requirement for Pax7-positive satellite cells in acute injury-induced skeletal muscle regeneration. Development. 2011, 138: 3639-3646. 10.1242/dev.067595.
Murphy MM, Lawson JA, Mathew SJ, Hutcheson DA, Kardon G: Satellite cells, connective tissue fibroblasts and their interactions are crucial for muscle regeneration. Development. 2011, 138: 3625-3637. 10.1242/dev.064162.
Sambasivan R, Yao R, Kissenpfennig A, Van Wittenberghe L, Paldi A, Gayraud-Morel B, Guenou H, Malissen B, Tajbakhsh S, Galy A: Pax7-expressing satellite cells are indispensable for adult skeletal muscle regeneration. Development. 2011, 138: 3647-3656. 10.1242/dev.067587.
Allen RE, Sheehan SM, Taylor RG, Kendall TL, Rice GM: Hepatocyte growth factor activates quiescent skeletal muscle satellite cells in vitro. J Cell Physiol. 1995, 165: 307-312. 10.1002/jcp.1041650211.
Tatsumi R, Anderson JE, Nevoret CJ, Halevy O, Allen RE: HGF/SF is present in normal adult skeletal muscle and is capable of activating satellite cells. Dev Biol. 1998, 194: 114-128. 10.1006/dbio.1997.8803.
Anderson JE: A role for nitric oxide in muscle repair: nitric oxide-mediated activation of muscle satellite cells. Mol Biol Cell. 2000, 11: 1859-1874.
Li YP: TNF-alpha is a mitogen in skeletal muscle. Am J Physiol Cell Physiol. 2003, 285: C370-C376.
Jones NC, Tyner KJ, Nibarger L, Stanley HM, Cornelison DD, Fedorov YV, Olwin BB: The p38alpha/beta MAPK functions as a molecular switch to activate the quiescent satellite cell. J Cell Biol. 2005, 169: 105-116. 10.1083/jcb.200408066.
Troy A, Cadwallader AB, Fedorov Y, Tyner K, Tanaka KK, Olwin BB: Coordination of satellite cell activation and self-renewal by Par-Complex-Dependent Asymmetric Activation of P38a/β MAPK. Cell Stem Cell. 2012, 11: 541-553. 10.1016/j.stem.2012.05.025.
Zhang K, Sha J, Harter ML: Activation of Cdc6 by MyoD is associated with the expansion of quiescent myogenic satellite cells. J Cell Biol. 2010, 188: 39-48. 10.1083/jcb.200904144.
Kuang S, Kuroda K, Le Grand F, Rudnicki MA: Asymmetric self-renewal and commitment of satellite stem cells in muscle. Cell. 2007, 129: 999-1010. 10.1016/j.cell.2007.03.044.
Le Grand F, Jones AE, Seale V, Scimè A, Rudnicki MA: Wnt7a activates the planar cell polarity pathway to drive the symmetric expansion of satellite stem cells. Cell Stem Cell. 2009, 4: 535-547. 10.1016/j.stem.2009.03.013.
Mauro A: Satellite cell of skeletal muscle fibers. J Biophys Biochem Cytol. 1961, 9: 493-495. 10.1083/jcb.9.2.493.
Ishikawa H: Electron microscopic observations of satellite cells with special reference to the development of mammalian skeletal muscles. Z Anat Entwicklungsgesch. 1966, 125: 43-63. 10.1007/BF00521974.
Pallafacchina G, François S, Regnault B, Czarny B, Dive V, Cumano A, Montarras D, Buckingham M: An adult tissue-specific stem cell in its niche: a gene profiling analysis of in vivo quiescent and activated muscle satellite cells. Stem Cell Res. 2010, 4: 77-91. 10.1016/j.scr.2009.10.003.
Lemons JM, Feng XJ, Bennett BD, Legesse-Miller A, Johnson EL, Raitman I, Pollina EA, Rabitz HA, Rabinowitz JD, Coller HA: Quiescent fibroblasts exhibit high metabolic activity. PLoS Biol. 2010, 8: e1000514-10.1371/journal.pbio.1000514.
Pritsker M, Doniger TT, Kramer LC, Westcot SE, Lemischka IR: Diversification of stem cell molecular repertoire by alternative splicing. Proc Natl Acad Sci U S A. 2005, 102: 14290-14295. 10.1073/pnas.0502132102.
Atlasi Y, Mowla SJ, Ziaee SA, Gokhale PJ, Andrews PW: OCT4 spliced variants are differentially expressed in human pluripotent and nonpluripotent cells. Stem Cells. 2008, 26: 3068-3074. 10.1634/stemcells.2008-0530.
Calvanese V, Lara E, Suárez-Alvarez B, Abu Dawud R, Vázquez-Chantada M, Martínez-Chantar ML, Embade N, López-Nieva P, Horrillo A, Hmadcha A, Soria B, Piazzolla D, Herranz D, Serrano M, Mato JM, Andrews PW, López-Larrea C, Esteller M, Fraga MF: Sirtuin 1 regulation of developmental genes during differentiation of stem cells. Proc Natl Acad Sci U S A. 2010, 107: 13736-13741. 10.1073/pnas.1001399107.
Rao S, Zhen S, Roumiantsev S, McDonald LT, Yuan GC, Orkin SH: Differential roles of Sall4 isoforms in embryonic stem cell pluripotency. Mol Cell Biol. 2010, 30: 5364-5380. 10.1128/MCB.00419-10.
Gabut M, Samavarchi-Tehrani P, Wang X, Slobodeniuc V, O’Hanlon D, Sung HK, Alvarez M, Talukder S, Pan Q, Mazzoni EO, Nedelec S, Wichterle H, Woltjen K, Hughes TR, Zandstra PW, Nagy A, Wrana JL, Blencowe BJ: An alternative splicing switch regulates embryonic stem cell pluripotency and reprogramming. Cell. 2011, 147: 132-146. 10.1016/j.cell.2011.08.023.
Chen JF, Mandel EM, Thomson JM, Wu Q, Callis TE, Hammond SM, Conlon FL, Wang DZ: The role of microRNA-1 and microRNA-133 in skeletal muscle proliferation and differentiation. Nat Genet. 2006, 38: 228-233. 10.1038/ng1725.
Kim HK, Lee YS, Sivaprasad U, Malhotra A, Dutta A: Muscle-specific microRNA miR-206 promotes muscle differentiation. J Cell Biol. 2006, 174: 677-687. 10.1083/jcb.200603008.
Flynt AS, Li N, Thatcher EJ, Solnica-Krezel L, Patton JG: Zebrafish miR-214 modulates Hedgehog signaling to specify muscle cell fate. Nat Genet. 2007, 39: 259-263. 10.1038/ng1953.
Lee RC, Feinbaum RL, Ambros V: The C. elegans heterochronic gene lin-4 encodes small RNAs with antisense complementarity to lin-14. Cell. 1993, 75: 843-854. 10.1016/0092-8674(93)90529-Y.
Lee RC, Ambros V: An extensive class of small RNAs in Caenorhabditis elegans. Science. 2001, 294: 862-864. 10.1126/science.1065329.
Lewis BP, Shih IH, Jones-Rhoades MW, Bartel DP, Burge CB: Prediction of mammalian microRNA targets. Cell. 2003, 115: 787-798. 10.1016/S0092-8674(03)01018-3.
Lewis BP, Burge CB, Bartel DP: Conserved seed pairing, often flanked by adenosines, indicates that thousands of human genes are microRNA targets. Cell. 2005, 120: 15-20. 10.1016/j.cell.2004.12.035.
Chen CZ, Li L, Lodish HF, Bartel DP: MicroRNAs modulate hematopoietic lineage differentiation. Science. 2004, 303: 83-86. 10.1126/science.1091903.
Johnston RJ, Chang S, Etchberger JF, Ortiz CO, Hobert O: MicroRNAs acting in a double-negative feedback loop to control a neuronal cell fate decision. Proc Natl Acad Sci U S A. 2005, 102: 12449-12454. 10.1073/pnas.0505530102.
Georgantas RW, Hildreth R, Morisot S, Alder J, Liu CG, Heimfeld S, Calin GA, Croce CM, Civin CI: CD34+ hematopoietic stem-progenitor cell microRNA expression and function: a circuit diagram of differentiation control. Proc Natl Acad Sci U S A. 2007, 104: 2750-2755. 10.1073/pnas.0610983104.
Ivey KN, Muth A, Arnold J, King FW, Yeh RF, Fish JE, Hsiao EC, Schwartz RJ, Conklin BR, Bernstein HS, Srivastava D: MicroRNA Regulation of Cell Lineages in Mouse and Human Embryonic Stem Cells. Cell Stem Cell. 2008, 2: 219-229. 10.1016/j.stem.2008.01.016.
Cheung TH, Quach NL, Charville GW, Liu L, Park L, Edalati A, Yoo B, Hoang P, Rando TA: Maintenance of muscle stem-cell quiescence by microRNA-489. Nature. 2012, 482: 524-528. 10.1038/nature10834.
Boutet SC, Cheung TH, Quach NL, Liu L, Prescott SL, Edalati A, Iori K, Rando TA: Alternative polyadenylation mediates microRNA regulation of muscle stem cell function. Cell Stem Cell. 2012, 10: 327-336. 10.1016/j.stem.2012.01.017.
Fukada S, Uezumi A, Ikemoto M, Masuda S, Segawa M, Tanimura N, Yamamoto H, Miyagoe-Suzuki Y, Takeda S: Molecular signature of quiescent satellite cells in adult skeletal muscle. Stem Cells. 2007, 25: 2448-2459. 10.1634/stemcells.2007-0019.
Echtermeyer F, Streit M, Wilcox-Adelman S, Saoncella S, Denhez F, Detmar M, Goetinck P: Delayed wound repair and impaired angiogenesis in mice lacking syndecan-4. J Clin Invest. 2001, 107: R9-R14. 10.1172/JCI10559.
Cornelison DD, Wilcox-Adelman SA, Goetinck PF, Rauvala H, Rapraeger AC, Olwin BB: Essential and separable roles for Syndecan-3 and Syndecan-4 in skeletal muscle development and regeneration. Genes Dev. 2004, 18: 2231-2236. 10.1101/gad.1214204.
Hall JK, Banks GB, Chamberlain JS, Olwin BB: Prevention of muscle aging by myofiber-associated satellite cell transplantation. Sci Transl Med. 2010, 2: 57-83.
Cornelison DD, Filla MS, Stanley HM, Rapraeger AC, Olwin BB: Syndecan-3 and syndecan-4 specifically mark skeletal muscle satellite cells and are implicated in satellite cell maintenance and muscle regeneration. Dev Biol. 2001, 239: 79-94. 10.1006/dbio.2001.0416.
da Huang W, Sherman BT, Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009, 4: 44-57.
da Huang W, Sherman BT, Lempicki RA: Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009, 37: 1-13. 10.1093/nar/gkn923.
Prifti E, Zucker JD, Clement K, Henegar C: FunNet: an integrative tool for exploring transcriptional interactions. Bioinformatics. 2008, 24: 2636-2638. 10.1093/bioinformatics/btn492.
Prifti E, Zucker JD, Clément K, Henegar C: Interactional and functional centrality in transcriptional co-expression networks. Bioinformatics. 2010, 26: 3083-3089. 10.1093/bioinformatics/btq591.
Antonov AV, Schmidt T, Wang Y, Mewes HW: ProfCom: a web tool for profiling the complex functionality of gene groups identified from high-throughput data. Nucleic Acids Res. 2008, 36: W347-W351. 10.1093/nar/gkn239.
Stark DA, Karvas RM, Siegel AL, Cornelison DD: Eph/ephrin interactions modulate muscle satellite cell motility and patterning. Development. 2011, 138: 5279-5289. 10.1242/dev.068411.
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. 10.1101/gad.1500707.
Smith RA, Meade K, Pickford CE, Holley RJ, Merry CL: Glycosaminoglycans as regulators of stem cell differentiation. Biochem Soc Trans. 2011, 39: 383-387. 10.1042/BST0390383.
Mourelatos Z, Dostie J, Paushkin S, Sharma A, Charroux B, Abel L, Rappsilber J, Mann M: Dreyfuss G: miRNPs: a novel class of ribonucleoproteins containing numerous microRNAs. Genes Dev. 2002, 16: 720-728. 10.1101/gad.974702.
Heinrichs A: MicroRNAs get a boost. Nat Rev Mol Cell Biol. 2009, 10: 302-
Johnson SM, Grosshans H, Shingara J, Byrom M, Jarvis R, Cheng A, Labourier E, Reinert KL, Brown D, Slack FJ: RAS is regulated by the let-7 microRNA family. Cell. 2005, 120: 635-647. 10.1016/j.cell.2005.01.014.
Cheung TH, Kwan YL, Hamady M, Liu X: Unraveling transcriptional control and cis-regulatory codes using the software suite GeneACT. Genome Biol. 2006, 7: R97-10.1186/gb-2006-7-10-r97.
Rocheteau P, Gayraud-Morel B, Siegl-Cachedenier I, Blasco MA, Tajbakhsh S: A subpopulation of adult skeletal muscle stem cells retains all template DNA strands after cell division. Cell. 2012, 148: 112-125. 10.1016/j.cell.2011.11.049.
Kops GJ, Medema RH, Glassford J, Essers MA, Dijkers PF, Coffer PJ, Lam EW, Burgering BM: Control of cell cycle exit and entry by protein kinase B-regulated forkhead transcription factors. Mol Cell Biol. 2002, 22: 2025-2036. 10.1128/MCB.22.7.2025-2036.2002.
Cáceres JF, Stamm S, Helfman DM, Krainer AR: Regulation of alternative splicing in vivo by overexpression of antagonistic splicing factors. Science. 1994, 265: 1706-1709. 10.1126/science.8085156.
Barberan-Soler S, Medina P, Estella J, Williams J, Zahler AM: Co-regulation of alternative splicing by diverse splicing factors in Caenorhabditis elegans. Nucleic Acids Res. 2011, 39: 666-674. 10.1093/nar/gkq767.
Figueroa A, Cuadrado A, Fan J, Atasoy U, Muscat GE, Muñoz-Canoves P, Gorospe M, Muñoz A: Role of HuR in skeletal myogenesis through coordinate regulation of muscle differentiation genes. Mol Cell Biol. 2003, 23: 4991-5004. 10.1128/MCB.23.14.4991-5004.2003.
van der Giessen K, Di-Marco S, Clair E, Gallouzi IE: RNAi-mediated HuR depletion leads to the inhibition of muscle cell differentiation. J Biol Chem. 2003, 278: 47119-47128. 10.1074/jbc.M308889200.
Young LE, Sanduja S, Bemis-Standoli K, Pena EA, Price RL, Dixon DA: The mRNA binding proteins HuR and tristetraprolin regulate cyclooxygenase 2 expression during colon carcinogenesis. Gastroenterology. 2009, 136: 1669-1679. 10.1053/j.gastro.2009.01.010.
Chen F, Shyu AB, Shneider BL: Hu antigen R and tristetraprolin: counter-regulators of rat apical sodium-dependent bile acid transporter by way of effects on messenger RNA stability. Hepatology. 2011, 54: 1371-1378. 10.1002/hep.24496.
Groszer M, Erickson R, Scripture-Adams DD, Dougherty JD, Le Belle J, Zack JA, Geschwind DH, Liu X, Kornblum HI, Wu H: PTEN negatively regulates neural stem cell self-renewal by modulating G0-G1 cell cycle entry. Proc Natl Acad Sci U S A. 2006, 103: 111-116. 10.1073/pnas.0509939103.
Hu Z, Wang H, Lee IH, Modi S, Wang X, Du J, Mitch WE: PTEN inhibition improves muscle regeneration in mice fed a high-fat diet. Diabetes. 2010, 59: 1312-1320. 10.2337/db09-1155.
Kimura E, Hidaka K, Kida Y, Morisaki H, Shirai M, Araki K, Suzuki M, Yamamura KI, Morisaki T: Serine-arginine-rich nuclear protein Luc7l regulates myogenesis in mice. Gene. 2004, 341: 41-47.
Yang T, Adamson TE, Resnick JL, Leff S, Wevrick R, Francke U, Jenkins NA, Copeland NG, Brannan CI: A mouse model for Prader-Willi syndrome imprinting-centre mutations. Nat Genet. 1998, 19: 25-31. 10.1038/ng0598-25.
Ito T, Saso K, Arimitsu N, Sekimizu K: Defective FESTA/EAF2-mediated transcriptional activation in S-II-deficient embryonic stem cells. Biochem Biophys Res Commun. 2007, 363: 603-609. 10.1016/j.bbrc.2007.09.009.
Yokoyama S, Ito Y, Ueno-Kudoh H, Shimizu H, Uchibe K, Albini S, Mitsuoka K, Miyaki S, Kiso M, Nagai A, Hikata T, Osada T, Fukuda N, Yamashita S, Harada D, Mezzano V, Kasai M, Puri PL, Hayashizaki Y, Okado H, Hashimoto M, Asahara H: A systems approach reveals that the myogenesis genome network is regulated by the transcriptional repressor RP58. Dev Cell. 2009, 17: 836-848. 10.1016/j.devcel.2009.10.011.
Qi Y, Hoshino M, Wada Y, Marubuchi S, Yoshimura N, Kanazawa I, Shinomiya K, Okazawa H: PQBP-1 is expressed predominantly in the central nervous system during development. Eur J Neurosci. 2005, 22: 1277-1286. 10.1111/j.1460-9568.2005.04339.x.
Bäumer D, Lee S, Nicholson G, Davies JL, Parkinson NJ, Murray LM, Gillingwater TH, Ansorge O, Davies KE, Talbot K: Alternative splicing events are a late feature of pathology in a mouse model of spinal muscular atrophy. PLoS Genet. 2009, 5: e1000773-10.1371/journal.pgen.1000773.
Baguma-Nibasheka M, Kablar B: Altered retinal cell differentiation in the AP-3 delta mutant (Mocha) mouse. Int J Dev Neurosci. 2009, 27: 701-708. 10.1016/j.ijdevneu.2009.07.005.
We thank Tiffany Antwine for technical assistance and members of the Olwin laboratory for critical discussion of this manuscript.
The project described was supported by the National Institutes of Health grant AR049446 to B. B. Olwin, a grant from the Muscular Dystrophy Association to B.B. Olwin and a Muscular Dystrophy Association Development Grant to D.D.W. Cornelison.
The authors do not have any competing interests.
NHF carried out bioinformatic analyses of the microarray data, designed and carried out the miRNA studies, and drafted the manuscript. MH carried out qRT-PCR for mRNA expression and many of the bioinformatic analyses of the microarray data. NDB carried out qRT-PCR for miRNA expression. CP carried out qRT-PCR for mRNA expression. DS carried out miRNA prediction. DDWC carried out satellite cell isolation and hybridization for the microarray. BBO conceived of the study, participated in its design and coordination, and helped to draft the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 7: Satellite cell activation gene expression profile. Excel spreadsheet (*.xlsx) with columns for probe set ID, log 2 expression data, ANOVA P value, fold change, and gene ID for unique genes that significantly change (ANOVA P ≤ 0.01, ≥ two-fold change) in WT and not in Sdc4−/−satellite cells 12 h after BaCl 2 -induced injury. The ANOVA P value and fold change was calculated for changes occurring in the first 12 h post-muscle injury according to genotype. (XLSX 870 KB)
Additional file 8: Molecular Function: RNA Binding genes regulated during satellite cell activation. Excel spreadsheet (*.xlsx) of RNA binding proteins from Additional file 7 with gene identifiers of Probe set ID, representative gene symbol, and entrez gene ID. The relative expression data for genes that significantly change (ANOVA P ≤ 0.01, ≥ two-fold change) in wild type and in Sdc4−/− satellite cells isolated from uninjured TA muscle or from the TA 12 h post-muscle injury is represented as a log2. (XLSX 68 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.