Bmi1 enhances skeletal muscle regeneration through MT1-mediated oxidative stress protection in a mouse model of dystrophinopathy

Enhanced polycomb complex protein Bmi1 expression in adult stem cells of the skeletal muscle leads to improved muscle function in a model of Duchenne Muscular Dystrophy via metallothionein1-mediated protection from oxidative stress.

Skeletal muscle is characterized by a remarkable capacity to regenerate after injury, mainly due to the function of satellite cells, the main skeletal muscle stem cells (Brack and Rando, 2012;Wang and Rudnicki, 2012). Polycomb group (PcG) proteins are essential regulators of stem cell function during normal development and in adult organs. They form multi-protein chromatin-associated complexes that play an essential role in the genome-wide epigeneticmediated remodeling of gene expression during myogenic differentiation of satellite cells, mainly through posttranslational modifications of histones (Asp et al., 2011). Ezh2 and Bmi1 are required for adult satellite cell homeostasis and proliferation in response to muscle injury, an effect mediated at least in part by repression of the ink4a locus (Juan et al., 2011;Robson et al., 2011). Importantly, although Bmi1 is expressed in several types of cancer and its mechanism of action may be similar in a non-neoplastic and neoplastic context, its overexpression does not initiate tumorigenesis (He et al., 2009;Yadirgi et al., 2011).
An emerging role for PcG proteins is their involvement in DNA repair (Liu et al., 2009;Facchino et al., 2010;Ismail et al., 2010;Ginjala et al., 2011;Pan et al., 2011). Bmi1 / -derived cells show significant mitochondrial dysfunction accompanied by sustained increase in reactive oxygen species (ROS) production that are sufficient to engage the DNA repair pathway (Liu et al., 2009), which is in turn impaired, thus leading to a magnified cellular damage. The balance between intracellular ROS and antioxidant molecules is vital in determining the rate of oxidative damage accumulation and the impaired function of satellite cells in aging and in myopathies, in which decreased anti-oxidative capacity has been documented (Fulle et al., 2005;Whitehead et al., 2006;Tidball and Wehling-Henricks, 2007). X-linked Duchenne muscular dystrophy (DMD) is the most common primary myopathy caused by the loss of the dystrophin protein from the plasma membrane, which causes loss of its integrity and fiber damage during repeated cycles of muscle degeneration and regeneration (Duncan, 1989). The proliferative capacity of myogenic cells was reported to be rapidly exhausted in dystrophin-deficient muscle, also because they are more sensitive to oxidative stress injury, leading to reduced and defective regeneration of the muscle as the disease progresses (Blau et al., 1983(Blau et al., , 1985Disatnik et al., 1998). Moreover, enzymatic adaptations to exercise-induced production of ROS and free radical damage are significantly decreased in dystrophic compared with normal muscles (Faist et al., 1998(Faist et al., , 2001. Overall, an impaired protection against ROS in dystrophic muscle appears to contribute to disease progression as also indicated by the beneficial, albeit transient, effect of antioxidants in ameliorating the skeletal muscle pathophysiology of DMD patients (Whitehead et al., 2008).
Metallothionein 1 (MT1) and MT2 are ubiquitously expressed (Kägi and Hunziker, 1989) low molecular weight, cysteine-rich zinc binding proteins. Although the role of MT1 in promoting cell proliferation is controversial (Smith et al., 2008), studies on MT-null liver cells showed their failure to regenerate after oxidative stress injury (Oliver et al., 2006).
Here, we show that overexpression of Bmi1 in the satellite cells significantly improves muscle strength through enhanced MT1-mediated protection of these cells from oxidative stress in a mouse model of dystrophinopathies but not after acute traumatic injury.

Bmi1 expression in mouse models of acute traumatic and chronic degenerative skeletal muscle injuries
To understand the potential impact of fine tuning Bmi1 expression in muscle injury, we characterized its expression profile in satellite cells at representative time points (3 and 10 d after i. and 10 d.a.i. (mean ± SEM from three independent experiments with n = 3 uninjured, n = 3 3 d.a.i., and n = 2 10 d.a.i.; ***, P < 0.001). (c) Immunofluorescence for Pax7 and Bmi1 on muscle sections from uninjured mice, injured mice (10 d.a.i.), and mdx mice (representative results from n = 4 uninjured, n = 2 3 d.a.i., n = 2 10 d.a.i., and n = 4 mdx). Arrows indicate Bmi1 staining intensity (expression level) in Pax7 +ve satellite cells in three different conditions (uninjured, 10 d.a.i., and mdx). (d) Quantification of Bmi1 staining of Pax7 +ve cells with pixel intensity above 400 at 3 and 10 d.a.i. compared with uninjured controls (mean ± SEM from two independent experiments with n = 4 uninjured, n = 2 3 d.a.i., n = 2 10 d.a.i., and n = 4 mdx; *, P < 0.05). (e) Bmi1 expression in satellite cells obtained from mdx single fibers after switching from growth to differentiation-inducing culturing conditions (48 h in differentiation conditions) compared with control satellite cells (mean ± SEM from two independent experiments with n = 3 for each condition; *, P < 0.05). Bar, 50 µm. (Fig. 1,c and d [quantification]). Interestingly, although Bmi1 expression increased significantly in control satellite cell cultures obtained from single fibers after switching from growth to differentiation-inducing conditions, no significant change was observed in mdx satellite cells (Fig. 1 e). These results show that the levels of Bmi1 expression are differentially regulated during the regenerative process in a pathological specific context.

Conditional Bmi1 expression increases the satellite cell pool and does not prevent myogenic differentiation in vitro
To explore the biological effect of enhancing Bmi1 expression in the satellite cells, we used a mouse line where Bmi1 expression can be activated in a cell-and time-specific fashion upon Cre-mediated recombination (Yadirgi et al., 2011). Transduction of satellite cell cultures derived from STOP-FloxBmi1 muscles with Adeno-Cre (A-Cre) virus induced efficient recombination and increased Bmi1 expression at protein injury [d.a.i.]) in a well-established model of acute traumatic muscle injury: the freeze injury model (Gayraud-Morel et al., 2007). Satellite cells were isolated 3 and 10 d.a.i. by magnetic activated cell sorting using SM/C-2.6 antibody (Fukada et al., 2004; Fig. 1 a). qRT-PCR analysis revealed significant increase in the expression of Bmi1 in these cells at both time points (Fig.1 b), a result which was confirmed by an increased percentage of Pax7 +ve cells expressing high level of Bmi1 in the muscle tissue, as assessed by double immunolabeling for Bmi1 and Pax7 and quantification of staining intensity (van der Laak et al., 2000;Gavet and Pines, 2010;Fig. 1, c and d [quantification]).
The mdx mouse is a well characterized model of chronic degenerative muscle pathology (Bulfield et al., 1984), which is genetically and biochemically similar to the human DMD (Collins and Morgan, 2003). We observed down-regulation of Bmi1 expression in Pax7 +ve satellite cells in the muscle tissue STOPFloxBmi1 single fibers, infected with A-GFP or A-Cre viruses and immunostained for Bmi1 (representative results from three independent experiments with n = 3 for each condition). (c and d) Levels of expression of Bmi1 are analyzed with WB (c; representative blot from 2 independent experiments with n = 4 biological samples) and qRT-PCR on the same samples (d; mean ± SEM from three independent experiments with n = 3 for each condition; ***, P < 0.001).
(e and f) Representative images (e; 6 independent experiments with n = 8 for each condition) of satellite cell cultures infected with either A-GFP or A-Cre viruses (Bmi1 Over cells; f) and stained for Pax7, MyoD, and MyHC, 2 and 5 d.a.i. in differentiation conditions. (g and h) Quantification of the findings in normal (g; mean ± SEM from six independent experiments with n = 8 for each condition; *, P < 0.05) and dystrophic context (h; mean ± SEM from 4 independent experiments with n = 4 for each condition; **, P < 0.01). (i and j) Representative images of the BrdU staining on A-GFP (i) and A-Cre-infected satellite cells (j) in four different conditions (six independent experiments with n = 8 for each condition): growth, 2 d differentiation, 5 d differentiation, and restimulation (i and j). (k) Quantification of BrdU +ve cells in growth, differentiation, and upon high serum restimulation conditions compared with A-GFPinfected satellite cells (mean ± SEM from six independent experiments with n = 8 for each condition; *, P < 0.05; ***, P < 0.001). (l) Quantification of Pax7 +ve satellite cells in Bmi1 Over satellite cells compared with A-GFP-infected satellite cells in growth, differentiation, and upon high serum restimulation conditions (mean ± SEM from four independent experiments with n = 4 for each condition; *, P < 0.05; **, P < 0.01). Bars, 62.5 µm. Bmi1 Over : Bmi1-overexpressing cells.

2620
Bmi1 overexpression in skeletal muscle regeneration | Di Foggia et al.
To investigate whether these results could be explained by modulation of satellite cell proliferation by Bmi1, satellite cells were treated with BrdU in four different culturing conditions: growth (proliferation medium), early differentiation (2 d in differentiation medium), late differentiation (5 d in differentiation medium), and restimulation (2 d in proliferation medium after 5 d of differentiation, Fig. 2, i and j). Under growth conditions, 59% of Bmi1 Over satellite cells had incorporated BrdU compared with 44% of the controls, and in early and late differentiation, 40% as compared with 5% and 8% as compared with 2%, respectively (Fig. 2,i,j,and k [quantification]). These results indicate that Bmi1 overexpression prolongs proliferation even in differentiation-inducing conditions. When satellite cell cultures were returned to high serum conditions after 5 d of differentiation (restimulation), Next, satellite cell cultures derived from single fibers isolated from the SOL muscle of postnatal day 60 (P60) STOPFloxBmi1 mice were infected with A-GFP or A-Cre virus and induced to differentiate by switching to low serum conditions for 2 d (Fig. 2, e-g). The number of Pax7 +ve cells was increased by 14% in cultures overexpressing Bmi1 (Bmi1 Over ), whereas no alteration in the proportion of MyoD +ve cells (29% in both conditions) was found (Fig. 2,e,f,and g [quantification]). The higher number of Pax7 +ve ;MyoD -ve cells in Bmi1 Over cultures suggested that the stem cell pool of satellite cells was increased (Fig. 2 g). Importantly, Bmi1 overexpression in satellite cells derived from P60 mice did not inhibit differentiation and the percentage of MyHC +ve cells remained constant as compared with controls ( Fig. 2, e, f, and g [quantification]). Similar results were obtained also when (a) Schematic temporal summary of the experiment. (b and c) Quantification of Bmi1 expression upon Tamoxifen injections in mdx;Pax7 Bmi1 mice compared with controls at RNA (b, mean ± SEM from four independent experiments with n = 12 mdx;Pax7 Bmi1 and n = 12 mdx; ***, P < 0.001) and at protein level (c; arrowheads indicate a Pax7 +ve cell then included in the insert, representative results from two independent experiments with n = 4 mdx, n = 3 mdx;Pax7 Bmi1 , n = 12 fields for mdx, and n = 11 fields for mdx;Pax7 Bmi1 ). (d) Quantification of Bmi1 pixel intensity in mdx;Pax7 Bmi1 mice compared with controls (mean ± SEM from two independent experiments with n = 4 mdx, n = 3 mdx;Pax7 Bmi1 , n = 12 fields for mdx, and n = 11 fields for mdx;Pax7 Bmi1 ; ***, P < 0.001). (e) The progression of the muscle pathology is followed by the forelimb grip strength test before each treadmill session. The table reports the statistical significant differences between groups and the coefficient representing the improving (if positive) or worsening (if negative) of the phenotype (four independent experiments with n = 5 controls, n = 9 mdx, n = 11 mdx;Pax7 Bmi1 ; mdx vs. mdx;Pax7 Bmi1 : **, P < 0.01; mdx vs. control: ***, P < 0.001; and control vs. mdx;Pax7 Bmi1 : *, P < 0.05). (f and g) Quantification of Pax7 +ve cells in mdx;Pax7 Bmi1 and mdx control (f) and quantification of double-positive cells Pax7 +ve /Ki67 +ve in mdx;Pax7 Bmi1 and mdx control (g; mean ± SEM from n = 3 mdx, n = 7 mdx;Pax7 Bmi1 ; *, P < 0.05). (h) Representative images of muscle sections stained for Pax7 and Ki67 in mdx;Pax7 Bmi1 as compared with mdx mice (n = 3 mdx, n = 7 mdx;Pax7 Bmi1 ; *, P < 0.05). Mdx;Pax7CreER T2 ;STOPFloxBmi1: mdx;Pax7 Bmi1 (Bmi1 overexpression in satellite cells in a dystrophic environment). Bars: (c) 125 µm; (h) 62.5 µm.
staining, revealed activation of the construct in 75% of the Pax7CreER T2 ;R26R injected mice (unpublished data). LacZ/Laminin double labeling confirmed the specificity of the recombination in satellite cells (unpublished data). FACS analysis revealed that at least 65% of the satellite cells isolated from Pax7CreER T2 ;R26R muscles by means of M-Cadherin labeling (Cooper et al., 1999) were -galactosidase positive, hence confirming high recombination efficiency (unpublished data).
To test the hypothesis that Bmi1 overexpression in satellite cells could enhance regeneration in a chronic injury model, triple transgenic mice mdx;Pax7CreER T2 ;STOPFloxBmi1, here called mdx;Pax7 Bmi1 , were generated. Mdx mice are more susceptible to exercise-induced injury, and a cycle of forced exercise on a treadmill has often been applied to accelerate muscle pathology and to test the efficacy of novel treatments in preclinical trials (De Luca et al., 2003;Whitehead et al., 2006). incorporation of BrdU was found in 12% of the cells, whereas almost 70% of the Bmi1 Over satellite cells incorporated BrdU after 2 d of growth restimulation (Fig. 2,i,j,and k [quantification]). Importantly, at this time point a higher number of Pax7 +ve cells was observed, in keeping with an increase in the pool of satellite cells even after a prolonged induction of differentiation ( Fig. 2 l).

Conditional activation of Bmi1 expression in satellite cells improves muscle strength in dystrophic muscles but not after traumatic injury in vivo
To assess the impact of Bmi1 overexpression on muscle regeneration in vivo, we used an inducible CreER line targeting satellite cells (Pax7CreER T2 ; Mourikis et al., 2012). Crosses with the R26R reporter mouse (Soriano, 1999), followed by Tamoxifen-mediated activation of Cre expression in double transgenic mice and whole mount LacZ mice compared with controls at RNA (b; mean ± SEM from two independent experiments with n = 5 controls and n = 3 Pax7 Bmi1 ; *, P < 0.05) and at protein level (via double immunolabeling for Pax7 and Bmi1; c; arrows indicate a Pax7 +ve cell that is included in the insert, representative results from one experiment with n = 3 controls and n = 3 Pax7 Bmi1 ); quantification of the average pixel intensity of Bmi1 staining in Pax7 +ve cells (d; mean ± SEM with n = 3 controls and n = 3 Pax7 Bmi1 ; ***, P < 0.001). (e) H&E staining of hind limbs 10 d.a.i. in injured controls and injured Pax7 Bmi1 mice (top, representative results from two independent experiments with n = 6 controls and n = 3 Pax7 Bmi1 ). Gomori Trichrome staining for endomysial fibrosis in the same samples of H&E (injured controls and injured Pax7 Bmi1 mice; bottom; representative pictures from one experiment with n = 3 controls and n = 3 Pax7 Bmi1 ). (f) Quantification of the ratio of regeneration in Pax7 Bmi1 mice compared with controls (f, mean ± SEM from two independent experiments with n = 6 controls and n = 3 Pax7 Bmi1 ). (g) Quantification of the CSA distribution of CNF in injured Pax7 Bmi1 and injured controls: the percentage of CNF of control mice (gray), Pax7 Bmi1 mice (blue), and the distribution of mature fibers (striped gray) is shown (mean ± SEM from two independent experiments and five independent experiments with n = 8 injured controls, n = 3 Pax7 Bmi1 , and n = 15 uninjured controls; *, P < 0.05; **, P < 0.01; ***, P < 0.001). (h) Treadmill performance in terms of the time run (h, mean ± SEM from three independent experiments with n = 12 uninjured controls, n = 7 injured controls and n = 3 Pax7 Bmi1 . **, P < 0.01) in injured control mice (gray) as compared with uninjured mice (striped), as well as in the distances run (i; mean ± SEM from three independent experiments with n = 7 injured controls, n = 3 Pax7 Bmi1 , and n = 12 uninjured controls; *, P < 0.05; **, P < 0.01). Bars, 125 µm.
was up-regulated (5.73 ± 2.03-fold) also in Bmi1 Over cultures upon restimulation (Fig. 6 a, last two bars on the right). IPA system biology analysis on all differentially expressed genes Triple transgenic mice, mdx, and nontransgenic littermates were injected with Tamoxifen, followed by a 6-wk period of running on a treadmill for 30 min twice weekly (Fig. 3 a). Bmi1 overexpression was confirmed both at RNA level in whole muscle tissue (Fig. 3 b) and at protein level specifically in the targeted satellite cells by double immunolabeling for Pax7/Bmi1 (Fig. 3, c and d [quantification]). Before each treadmill session, a forelimb grip test was performed to follow the progression of the disease. The forelimb strength of mdx mice was clearly compromised compared with controls ( Fig. 3 e) but the strength of both groups (mdx and control littermates) was stable throughout the duration of the experiment. Strikingly, mdx;Pax7 Bmi1 mice showed a significant improvement in muscle strength from the fourth week onwards (Fig. 3 e).
Histological examination of the muscle did not reveal morphological differences between mdx and mdx;Pax7 Bmi1 , in particular an impact on the degree of endomysial fibrosis or on adipose tissue dedifferentiation was excluded (unpublished data). Immunolabeling for Pax7 showed an increased percentage of Pax7 +ve cells in mdx;Pax7 Bmi1 mice as compared with mdx controls (Fig. 3, h and f [quantification]). Double labeling for Pax7 and the proliferation marker Ki67 revealed that it was the Ki67 ve component of the Pax7 +ve population which was increased in mdx;Pax7 Bmi1 (Fig. 3, h and g [quantification]), suggesting that at this time point overexpression of Bmi1 in Pax7 +ve cells had increased the pool of quiescent satellite cells in a dystrophic environment. Interestingly, when a muscle freeze injury (FZ) was performed on the TA of hind limbs of Pax7CreER T2 ;STOPFloxBmi1 compound mutant mice-here called Pax7 Bmi1 -after Tamoxifen injection ( Fig. 4 a) to induce Bmi1 overexpression ( Fig. 4, b-d), no significant differences were found morphologically (Fig. 4, e-g) or functionally (Fig. 4, h and i) in injured Pax7 Bmi1 as compared with injured controls.

Transcriptome analysis reveals increased expression of MT1 in satellite cells upon activation of Bmi1 expression
To gain insights into the downstream mechanisms mediating the observed phenotype, we first analyzed the canonical Bmi1 downstream target genes p16 ink4a , p19 Arf , and p21 waf1/cip1 . No changes were found in the expression level of these genes in Bmi1 over satellite cells cultured under differentiation-inducing conditions or upon restimulation (unpublished data). Alternative ink4a/p21-independent mechanisms were then explored by whole-genome transcriptome analysis on the same cultures (Illumina platform mouse v2). A set of 254 genes was identified (Table S1), the expression of which was deregulated in satellite cell cultures overexpressing Bmi1 compared with controls (P < 0.05 and absolute fold change > 1.9). Out of 10 genes (Ahnak, Bmi1, DNAjb6, Ehd1, Hspd1, Itgb1, Klf6, MT1, p21 waf1/cip1 , and Tpm4) selected for validation analysis, all but one (p21 waf1/cip1 ) were confirmed to be differentially expressed (Figs. 5 a and 6 a; and not depicted). MT1 was the most significantly up-regulated gene (5.27 ± 0.99-fold) in these cultures (Fig. 6, a and b) and importantly, its expression cell cultures compared with A-GFP-infected controls (mean ± SEM from two independent experiments with n = 3 for each condition; *, P < 0.05; **, P < 0.01; ***, P < 0.001). (b and c) System biology analysis using the Ingenuity platform reporting genes in the top canonical pathways (b) and the top tox lists (c; n = 3 for each condition). (d and e) qRT-PCR analysis of expression of redox homeostasis genes in Bmi1 Over satellite cell cultures compared with A-GFP-infected controls (mean ± SEM from three independent experiments with n = 5 for each condition; *, P < 0.05; **, P < 0.01) and in Bmi1 / satellite cell cultures as compared with controls (e; mean ± SEM from two independent experiments with n = 4 for each condition; *, P < 0.05; **, P < 0.01; ***, P < 0.001).

MT1 is a crucial mediator of the enhanced satellite cell pool elicited by Bmi1 overexpression in vitro
Because the DMD pathology is linked to accumulation of ROS in satellite cells with a consequent increase in the oxidative stress induced cellular damage Zhuang et al., 2001;Wozniak and Anderson, 2007), we analyzed whether Bmi1 overexpression would exert a protective role against oxidative stress-induced cellular damage in satellite cells. We observed a decrease in ROS production in Bmi1 Over satellite cell cultures incubated with the fluorescent dye H 2 DCFDA, followed by flow cytometry (Fig. 6 c, white and black bars), which led to decreased oxidative damage to the cells as assessed by 8-hydroxydeoxyguanosine (8-OHdG) staining, a marker of oxidative stress-dependent DNA damage (Won et al., 1999;Fig. 6, e and d [quantification]).
showed hypoxia signaling in the cardiovascular system, p53 signaling, and protein ubiquitination pathway (Fig. 5 b) to be the most significantly deregulated top canonical pathways. Hypoxia-inducible factor signaling and NRF2-mediated oxidative stress response were most significantly affected among the top tox lists (Fig. 5 c). To further extend and validate these findings, the expression of a selection of genes involved in Bmi1-mediated control of redox homeostasis through generation of ROS (Liu et al., 2009) was analyzed. We found a significant repression of Alox15, Bnip3L, Cdo1, and Duox2 (Fig. 5 d). Down-regulation of MT1 expression (Fig. 6 a, first two bars; and Fig. 7 a) and increased expression of Alox5, Bnip3L, Cdo1, Duox1, and Duox2 in Bmi1 / satellite cells confirmed the specificity of these associations (Fig. 5 e). (a) Relative expression level of MT1 in Bmi1 / , Bmi1 Over , and control satellite cells (mean ± SEM from two independent experiments with n = 3 for each condition; *, P < 0.05; **, P < 0.01). (b) Expression levels of MT1 are also analyzed at the protein level in Bmi1 Over satellite cells both in growth and differentiation-inducing conditions (left, three independent experiments with n = 3 for each condition). Efficient knock down of MT1 is also shown in proliferation (P) and differentiation (D) conditions (b, right; representative results from three independent experiments with n = 3 for each condition). WB analysis of the levels of 3-NT modification in Bmi1 Over satellite cells both in proliferation and differentiation-inducing conditions (b, left; three independent experiments with n = 3 for each condition). (c) The amount of ROS production is measured via DCFDA relative fluorescence intensity in Bmi1 Over satellite cells compared with controls and treated with either sh Scramble or knockdown of MT1 (mean ± SEM from two independent experiments with n = 4 for each condition; *, P < 0.05; ***, P < 0.001). (d and e) Quantification of the percentage of strongly 8OHdG +ve satellite cells in Bmi1 Over cultures with and without knockdown of MT1 (mean ± SEM from two independent experiments with n = 4 for each condition; *, P < 0.05) and representative examples of immunolabeling (e) are shown. (f) Representative examples of immunolabeling for Pax7 and Edu in Bmi1 Over satellite cell cultures upon MT1 knockdown (two independent experiments with n = 4 for each condition). (g) Quantification of the percentage of Edu +ve or Pax7 +ve or Edu +ve /Pax7 +ve cells is shown (mean ± SEM from two independent experiments with n = 4 for each condition; *, P < 0.05). (h) Total concentration of nitrosylated proteins (ng/µl) in Bmi1 Over mdx;Bmi1 Over satellite cells compared with A-GFP-infected controls, with and without Mt1 knockdown (mean ± SEM from two independent experiments with n = 4; ***, P < 0.001; *, P < 0.05). (i and j) The level of SNO-Mef2c is shown in mdx;Bmi1 Over satellite cells compared with controls A-GFP-infected satellite cells as assessed by a modified biotin switch assay, followed by immunoblot (i) and quantification (j; mean ± SEM from two independent experiments with n = 4; *, P < 0.05). cells isolated from Bmi1 / mice showed increased ROS production ( Fig. 7 b), increased number of 8-OHdG +ve cells (Fig. 7, c and d [quantification]), and increased levels of 3-NT in Bmi1 / muscles and satellite cell cultures (Fig. 7 a), hence supporting the interpretation that the observed protection from oxidative stress is a Bmi1-dependent effect.

Bmi1-driven improved muscle strength in dystrophic mice is mediated by MT1 up-regulation
Next, we assessed whether up-regulation of MT1 mediates the improved muscle strength observed in the mdx model upon Bmi1 overexpression in satellite cells. MT1 up-regulation was detected in 33% of mdx;Pax7 Bmi1 (Fig. 8 a). Importantly, no MT1 up-regulation was ever detected in mdx control mice (Fig. 8 a), in keeping with a Bmi1-induced effect. Reevaluation of the muscle strength data with focus on MT1 expression revealed that mdx;Pax7 Bmi1 mice with MT1 upregulation were the animals with the best improvement in forelimb grip strength (Fig. 8 b) as compared with mdx control mice. mdx;Pax7 Bmi1 mice up-regulating MT1 showed a behavior similar to that of unaffected control mice, whereas mdx;Pax7 Bmi1 with no MT1 up-regulation showed a performance comparable to mdx control mice (Fig. 8 b). These findings show that Bmi1-induced MT1 up-regulation molecularly mediates the impact of Bmi1 overexpression on muscle strength. Furthermore, centrally nucleated fiber (CNF) percentage was significantly reduced in mdx;Pax7 Bmi1 with MT1 up-regulation compared with mdx control mice and to mdx;Pax7 Bmi1 mice without MT1 up-regulation (Fig. 8,c and d [quantification]), suggesting that muscle fibers were more resistant to exercise-induced damage in mdx;Pax7 Bmi1 with MT1 up-regulation. The reduced number of H2AX +ve (a marker of DNA double-strand break) satellite cells detected in these mice (double immunofluorescence for Pax7 and H2AX; Fig. 8, e and f [quantification]) is in keeping with the protective effect having been exerted specifically on the satellite cells. We did not observe differences in the number of apoptotic These findings suggest that the amount of oxidative damage incurred by satellite cells in vitro depends on Bmi1 expression levels and raise the possibility that this effect is mediated by MT1. To test this hypothesis functionally, we set out to knock down MT1 expression (MT1kd) in both Bmi1 Over and mdx;Bmi1 Over cultures and we analyzed ROS production, the percentage of 8-OHdG +ve cells, and the number of proliferating satellite cells and/or Pax7 +ve cells. Western blot analysis confirmed efficient knockdown of MT1 in mdx;Bmi1 Over cultures in growth and in differentiation conditions (Fig. 6 b). A significant increase in ROS production was observed upon MT1kd in Bmi1 Over and in mdx;Bmi1 Over satellite cells (Fig. 6 c and not depicted) with concomitant increase in the percentage of 8-OHdG +ve cells (Fig. 6, d and e) and significantly fewer EdU +ve , Pax7 +ve and EdU +ve ;Pax7 +ve cells (Fig. 6, f and g [quantification]), compared with controls treated with scrambled sequences.

2626
Bmi1 overexpression in skeletal muscle regeneration | Di Foggia et al.
satellite cells or in the total number of apoptotic nuclei (double immunolabeling for Pax7 and cCaspase3; Fig. 8 h and not depicted), or in the percentage of necrotic fibers (not depicted), but a generalized decrease in mitochondrial oxidative enzyme activities was detected by means of in vitro biochemical analysis in mdx;Pax7 Bmi1 with MT1 up-regulation (Fig. 8 g). To functionally validate these findings, we transplanted satellite cells isolated from mdx;STOPFloxBmi1 mice induced to overexpress Bmi1 in vitro and either up-regulating MT1 (Bmi1-induced) or not (ShRNA-MT1 mediated knockdown) into the TA muscle of young mdx mice. We show that the contribution of the cells to the regenerative process was more efficient when MT1 was up-regulated, as demonstrated by the reduced number of donor satellite cell-derived GFP-positive fibers upon MT1 knockdown (Fig. 9, a and b [quantification]).

Bmi1 expression is reduced in patients affected by DMD
To translate the significance of these findings to human pathology, we analyzed the expression of Bmi1 in DMD muscles and in morphologically normal controls. Four DMD patients aged 2.2, 3.2, 6.1, and 9.5 yr and six age-matched controls were used. Reduction of Bmi1 expression was found in DMD muscle biopsies at all ages (Fig. 10 a). Immunolabeling for Bmi1 performed on adjacent sections confirmed a reduction in the number of Bmi1 +ve cells in DMD cases (Fig. 10, a and b) and double labeling for Bmi1 and Pax7 revealed lower levels of Bmi1 expression in Pax7 +ve satellite cells in Figure 9. MT1 up-regulation is essential for the efficient contribution of Bmi1 Over satellite cells to muscle regeneration upon transplantation in the mdx mouse. (a) Satellite cells are isolated from mdx;STOPFloxBmi1 mice, cultured for 14 d, infected with A-Cre virus to activate Bmi1 expression, and treated either with Sh scramble and Mt1 knockdown, both viruses expressing GFP in infected cells. 6 d after treatment, the cells are transplanted in the TA muscles of mdx recipient mice. The transplanted muscles are collected after 27 d. Representative pictures of mdx muscles transplanted with A-Cre; mdx;STOPFloxBmi1 + Sh Scr satellite cells (a, left) or A-Cre; mdx;STOPFloxBmi1 + Sh MT1 satellite cells (a, right) are shown (a, representative results from n = 4 for each condition). (b) Quantification of GFP +ve fibers in A-Cre; mdx;STOPFloxBmi1 + Sh MT1 transplanted satellite as compared with A-Cre; mdx;STOPFloxBmi1 + Sh MT1 transplanted satellite (mean ± SEM from one experiment with n = 4 A-Cre;Mdx;STOPFloxBmi1+Sh Scr and n = 4 Mt1 A-Cre; Mdx;STOPFloxBmi1+ShMt1; *, P < 0.05). Bar, 125 µm.

Figure 10. Bmi1 expression is reduced in the muscle of human DMD cases. (a) H&E
staining of a DMD case as compared with an age-matched morphologically normal control. The same muscles are double stained for Bmi1 and Pax7 (representative results from two different cases). (b) Quantification of Bmi1 +ve cells in DMD muscle biopsies and age matched controls (mean ± SEM from two independent experiments with n = 2 control and DMD in both age groups, five high power fields/case; *, P < 0.05; **, P < 0.01). (c) Levels of expression of Bmi1 are analyzed by qRT-PCR in the muscle biopsies of DMD patients as compared with age matched unaffected muscles (mean ± SEM from two independent experiments, n = 3 control and n = 2 DMD in both age groups, two technical replicas for each case; *, P < 0.05; **, P < 0.01). (d) Expression levels of MT1 are measured by qRT-PCR in the muscle biopsies of DMD patients >5 yr compared with age-matched unaffected muscles (mean ± SEM from two independent experiments with n = 3 control and n = 2 DMD in both age groups, two technical replicas for each case; *, P < 0.05; **, P < 0.01). Bar: (H&E) 62.5 µm; (immunofluorescence) 31 µm. pool seen in dystrophic muscle conditions after repeated degeneration/regeneration cycles. Genome-wide transcriptome analysis of satellite cells overexpressing Bmi1 revealed significantly increased expression of MT1, which has been shown to be necessary to maintain proliferation and normal cellular functions under stress conditions, including when cells reenter the cell cycle from a quiescent state (Zbinden et al., 2010). In keeping with this notion, improvement of muscle strength is seen only in those mdx mice in which up-regulation of MT1 occurred concomitantly with overexpression of Bmi1 in satellite cells. The reduced contribution to regeneration after transplantation in a dystrophic environment of Bmi1 Over satellite cells upon MT1 knockdown lends additional support to this conclusion. At present it is unclear why up-regulation of MT1 occurs only in a proportion of mdx mice overexpressing Bmi1 and it will be important to investigate this fully in future studies. However, the correlation of muscle strength and MT1 up-regulation suggests that Bmi1 exerts its effect through protection against oxidative stress-induced cellular damage in vitro and in vivo, hence delaying muscle wasting through reduction of ROS-induced oxidative stress. The reduction of H2AX +ve satellite cells observed in the mdx;Pax7 Bmi1 with concomitant MT1 up-regulation is in keeping with this conclusion, as it is the MT1-dependent reduction of 3-NT levels in mdx;Pax7 Bmi1 satellite cell cultures and the reduced levels of nitrosylated Mef2c, a transcription factor with an essential role in muscle regeneration (Liu et al., 2014;Panda et al., 2014). We also observed a decrease in metabolic and mitochondrial oxidative enzyme activities in these mice, raising the possibility that because satellite cells are better protected against exercise-induced oxidative stress in a dystrophic environment, metabolic adaptation is less prominent. Although additional studies will be required to fully elucidate the impact of Bmi1 overexpression on metabolism, it is intriguing that decreased mitochondrial oxidative phosphorylation activities have been shown to increase resistance to metabolic diseases through prevention of ROS accumulation (Pospisilik et al., 2007) and to be associated with better tolerance to oxidative stress (Patel et al., 2013). It is currently unclear what mediates Bmi1-induced up-regulation of MT1 expression in the skeletal muscle. As the main mode of function of PcG genes is through transcriptional repression, it is possible to speculate that a MT1 repressor may be involved in this process.
Finally, we show here an important correspondence of Bmi1 expression between mice and humans, including its depletion in human DMD muscles. Interestingly there is evidence that MT1 expression increases in human skeletal muscle 24 h after intense exercise (Urso et al., 2006); here, we show decreased expression of MT1 in human DMD muscle, raising the possibility that manipulation of Bmi1 expression could be a useful approach to enhance satellite cell-driven regeneration also in human primary myopathies. Identifying novel pharmacological agents to enhance and sustain muscle regeneration in DMD is highly desirable as they could have a tremendous impact on DMD therapy either alone or, most DMD patients (Fig. 10 a). Importantly, reduced expression of Bmi1 (Fig. 10 c) and MT1 (Fig. 10 d) was detected also in DMD patients at the RNA level.

DISCUSSION
Here, we show for the first time that mild overexpression of a chromatin modifier gene has a functionally relevant impact on muscle regeneration. Bmi1 overexpression in satellite cells confers enhanced protection from exercise-induced injury and improves muscle strength in a mouse model of DMD. We also demonstrate that the mechanism underlying the effect of enhanced physiological function of Bmi1 in the satellite cells is mediated by MT1-driven modulation of resistance against oxidative stress.
Skeletal muscle is an important and highly relevant platform to assess how manipulation of stem cell properties may enhance a functionally significant regenerative process in vivo. We show that the expression level of Bmi1 is modulated in response to pathological conditions in the satellite cells. Indeed, Bmi1 is up-regulated upon freeze injury in wild-type mice, whereas reduction of the proportion of Pax7 +ve cells with high Bmi1 expression is found in the mdx mouse, with the latter finding corroborated by the lack of Bmi1 up-regulation in mdx satellite cell cultures upon induction of differentiation. In keeping with these findings, overexpression of Bmi1 in satellite cells has a beneficial effect solely in the mdx mouse, where downregulation of its expression was observed.
The best characterized molecular mechanisms mediating Bmi1 function relies on transcriptional repression of the ink4a locus, encoding for p16 ink4a and p19 arf ( Jacobs et al., 1999), as well as of p21 waf1/cip1 . However, we found no change in the expression levels of these genes upon activation of Bmi1 expression raising the intriguing possibility that a novel mechanism may mediate the observed phenotype. The generation of free radicals and ROS/reactive nitrogen species (RNS) is a physiological, continuous process in aerobic organisms. When the production of ROS/RNS exceeds the endogenous antioxidant buffering capacity, oxidative stress-induced cellular damage occurs. There is increasing evidence that oxidative stress level is elevated in diseases such as DMD and in normal aging (Fulle et al., 2004;Whitehead et al., 2006;Jackson, 2009;Cozzoli et al., 2011;Spassov et al., 2011). Bmi1 has an important role in mitochondrial function and redox homeostasis in thymocytes and in neurons (Liu et al., 2009;Abdouh et al., 2012); however, the mechanisms mediating this role of Bmi1 are currently unclear. Up-regulation of Bmi1 expression in Pax7 +ve satellite cells increases the pool of quiescent Ki67 ve cells in vivo after 6 wk of training on a treadmill, indicating that Bmi1 is involved in the regulation of self-renewal of these cells. To explain these results, two scenarios can be envisaged: it is possible that a phase of increased proliferation has also occurred in vivo, perhaps immediately after induction of Bmi1 overexpression in a similar fashion as we observed in vitro. Alternatively, satellite cells overexpressing Bmi1 are better able to reenter quiescence, hence preventing the typical exhaustion of the satellite cell 2628 Bmi1 overexpression in skeletal muscle regeneration | Di Foggia et al.

Microscopy and quantification.
Fluorescent and bright field image capture was performed using an epifluorescent (Leica) or confocal (Meta 510 LSM or LSM 710; Zeiss) microscope. In myoblast cultures, all DAPIpositive cells and antigen-positive cells were counted and the percentage of Antigen +ve cells was calculated. Data from three or more cultures were pooled to give a population mean ± SEM.
For the quantification of the intensity staining, all samples were stained simultaneously and the pictures were acquired with an LSM 710 confocal microscope at the same ratio of emissions. For the analysis with ImageJ software (National Institutes of Health), the shape of each Pax7 +ve cell was drawn around the nucleus and the intensity of Bmi1 staining calculated for each cell as mean of integrated density (mean intensity × area). The mean integrated density was calculated among all cells per each field acquired. The percentage of cells above or below certain pixel intensity (determined by the mean integrated density between all measured cells) was calculated out of the total number of measured cells per each muscle section analyzed. The method was also validated using Definiens Tissue Studio (IF) 2.1.1, action library 3.6.1, and Developer XD 2.1.1 (Definiens AG). In particular for the validation on sections stained for Pax7 and Bmi1, an initial analysis was performed using Tissue Studio to identify the nuclei based on DAPI stain and Pax7. Nuclei were identified using a stain threshold of 1.62 (RGB space converted to HSD model to give a range of 0-2; van der Laak et al., 2000) and typical size of 18 µm 2 . Pax7 was processed twice using thresholds of 0 and 25 (RGB space) to classify all nonnuclear tissue (Marker 1) and areas with significant Pax7 staining (Marker 2), respectively. These are overlaid to give classifications of marker 1 and markers 1 and 2. Developer was used to perform post-processing on the classification resultant from Tissue Studio. Before classification of nuclei as positive or negative, nuclei with area <5 µm 2 were merged with connecting nuclei. In cases where these small nuclei did not border another nucleus, no change occurred. In cases where a small nucleus borders multiple nuclei, the merge resulting in the lowest elliptical fit was selected. Several dynamic thresholds were calculated for each image: the mean Pax7 stain intensity for each area classified as marker 1 was used to give Pax7 Background; areas of marker 1 <10 µm 2 were excluded from this calculation to prevent skewing of resultant threshold toward the upper limit of 24 because these are generally isolated areas within the intensely stained regions (markers 1 and 2) just below the threshold of 25. The threshold for classifying nuclei as positive Pax7 threshold was calculated as: maximum Pax7 pixel value  (2 × Pax7 Background). A second set of thresholds were calculated based on the Bmi1 stain to segment positive nuclei staining as below mean, low intensity, and high intensity. The mean Bmi1 stain intensity in all pixels with intensity <50 (RGB space) gives mean Bmi1, and the range between mean Bmi1 and the maximum Bmi1 value was divided by two to give low Bmi1 and high Bmi1. Nuclei with more than five pixels with Pax7 greater than Pax7 threshold were classified as positive nuclei. Positive nuclei were then segmented based on BMI1 stain as below mean, low intensity, and high intensity. To compare this analysis with the manual analysis with ImageJ, the mean Bmi1 intensity for each single slice of the stack where Bmi1 was on focus was multiplied to the area (in pixel) so to obtain integrated density values. The percentage of cells above or below certain pixel intensity (determined by the mean integrated importantly, in combination with genetic therapies aiming at reconstituting dystrophin expression.

Satellite cell culture.
Single fibers from the soleus (SOL) and the extensor digitorum longus (EDL) muscles of adult Bmi1 / , STOPFloxBmi1, mdx;STOPFloxBmi1, and nontransgenic littermates were isolated and cultured to obtain pure satellite cell cultures according to standard protocols. In brief, the muscles were dissected from P60 mice being careful not to stretch the muscles. Tendon and other connective tissues were carefully removed and the muscles incubated in 0.2% Collagenase type 1 (Worthington Biochemical Corporation) at 37°C for 2 h. Fibers were liberated by gentle trituration in DMEM. Single fibers were placed in a 24-well plate (BD) coated with 1mg/ml Geltrex Reduced Growth Factor Basement Membrane Matrix in DMEM (Invitrogen) and allowed to adhere for 5 min at 37°C before 500 µl of plating medium was added (DMEM supplemented with 10% horse serum [PAA] and 0.5% chick embryo extract [Gibco] with antibiotics) and incubated at 37°C in 5% CO 2 . After 3 d, the media was changed to proliferation media (DMEM supplemented with 20% FCS [PAA] and 2% chick embryo extract [Gibco] plus antibiotics). To induce differentiation, satellite cells were plated at 2,000 cells per well in an 8-well chambered slides (BD) and cultured in differentiation media (DMEM supplemented with 2% horse serum and antibiotics) for 2 or 5 d. Satellite cells derived from the same fiber of STOPFloxBmi1 and mdx;STOPFloxBmi1 mice were set up in parallel wells and treated with either A-GFP or A-Cre (Akagi et al., 1997) at the time of plating in the differentiating media. To control for nonspecific effects of viral treatment and overexpression of any exogenous proteins, infection with A-GFP was used as control for all experiments where A-Cre infection was performed.
Magnetic activated satellite cells sorting. TA muscles were isolated from C57BL/6J mice at time zero (uninjured controls), and 3 and 10 d after freeze injury. The muscles were dissected and finely chopped with scissors and incubated with 0.2% Collagenase in DMEM for 30 min at 37°C while gently shaking (250 rpm). This step of the procedure was performed twice. The cells were then filtered with 30-µm filters (Miltenyi Biotec). 1 ml of serum was added and the suspension was centrifuged for 20 min at 1,200 rpm at 4°C. The supernatant was decanted and the staining was performed with the SM/C-2.6 antibody (Fukada et al., 2004;1:200 for 40 min on ice) and with anti-Biotin MicroBeads (Miltenyi Biotec) as secondary antibody according to manufacturer's protocol. After magnetic sorting the cells were centrifuged 20 min at 2,500 rpm and the pellet was snap frozen on dry ice.
When satellite cells were analyzed, cell pellets were collected from different culture conditions and lysed with RIPA buffer as mentioned above. 40 µg of protein extracted from cell samples was loaded per lane on the Nu-Page 4-12% Bis-Tris gel.
Production of shRNA-Mt1 lentiviral vector. A pGIPZ-lentiviral shRNAmir vector containing a hairpin sequence targeting MT1 is commercially available from Open Biosystems. The shRNA-containing lentiviral vector was cotransfected with lentiviral packaging, using the calcium phosphate method into 293T cells. The medium was refreshed and viruses were harvested after 48 h of transfection, passed through 0.45-µm filters, concentrated by PEG precipitation, and stored at 80°C. The infectious titer was determined by FACS analysis of GFP-positive 293T cells.
Lentiviral-mediated Mt1 knockdown. Satellite cells were seeded in multiple-well plates and transduced with concentrated lentivirus containing shRNA that targets MT1, along with A-Cre (to induce Bmi1 overexpression) in proliferation medium. The medium was replaced on the second day. 96 h after transduction, cells were switched to differentiation medium for 48 h and the transduction efficiency was analyzed. Cells were collected for Western blot, immunocytochemistry, or ROS assay.

ROS assay.
Satellite cells were treated with 40 µM H 2 DCFDA in differentiation medium for 1 h at 37°C before being harvested and resuspended in 500 µl of PBS. Cells were then immediately processed for flow cytometry analysis using an LSRII instrument (BD) and the relative fluorescence intensity was calculated.

ELISA for 3-NT detection. Satellite cells prepared from mdx;STOPFloxBmi1
EDL muscles were transfected with A-Cre to overexpress Bmi1 and a proportion then treated with LV-shMt1. Cells infected with A-GFP were used as controls. Proteins were extracted from the cells and concentration was measured by BCA protein assay kit (Thermo Fisher Scientific). The amount of 3-NT protein in the lysates was determined by the Nitrotyrosine ELISA kit (Abcam) according to the manufacturer's protocol. The OD was recorded at 450 nm and the data were analyzed with the program XLfit 5.3.1.3 (IBDS). density between all measured cells) was calculated out of the total number of measured cells per each muscle section analyzed.
For the analysis of satellite cells in vitro, satellite cells were stained for Bmi1 and counterstained with DAPI. Immunofluorescent images taken at 16-bit were analyzed using ImageJ and Developer XD (2.1.1; Definiens AG). The analysis in ImageJ was conducted using the same settings as above. For the analysis in Definiens, pixels with intensity above 30 (RGB space) were classified as nucleus, those below 30 remained unclassified; and small unclassified areas bordering only nucleus were reclassified as nucleus and merged in. The distance of each pixel classified as nucleus to unclassified was then calculated to facilitate the splitting of connected nuclei. Nuclei with irregular shapes (roundness >0.5; roundness is based on the difference between the largest enclosing ellipse and the smallest enclosed ellipse) were then segmented based on these distances; local minima provide seeds for nuclei. These seeds were then grown to fill the space originally classified as nuclei, with separate nuclei growing to share a border rather than joining, thus splitting connected nuclei. Nuclei with area < 100 µm 2 were then removed. For each nucleus, the sum of the Bmi1 staining for each pixel was exported and this value was multiplied for the area of each single nucleus. For both the analyses, the fold change of increased Bmi1 staining was calculated using the average of Bmi1 staining of A-GFP-infected cells as one per each coverslip analyzed.

RNA extraction and qPCR analysis.
Satellite cell cultures were harvested after 48 h of either differentiation or restimulation. The cell pellets were lysed in RLT buffer from RNeasy Micro purification kit (QIAGEN). Total RNA was isolated using the same kit. Digestion with DNase (on column) was performed to remove genomic DNA. cDNA synthesis was performed with SuperScript III Reverse transcription kit (Invitrogen). For the sorted satellite cells, the retrotranscription was performed with QuantiTect Whole Transcriptome kit (QIAGEN) according to manufacturer's protocols. Analysis of gene expression was performed with a Rotor-Gene apparatus (Corbett Life Science) using SYBR Green JumpStart TaqReadyMix (Sigma-Aldrich) according to the manufacturer's instruction. A list of the primers used is reported. The Ct values of all the genes analyzed were normalized with GAPDH and fold changes were calculated. qRT-PCR analysis for Bmi1 was performed with TaqMan assays (Applied Biosystems) with FAM-labeled probes in 96-well plates using the 7500 RT-PCR machine (Applied Biosystems) according to the manufacturer's instructions (Bmi1: Mm 00776122_gH). The cDNA content was normalized against the expression of the housekeeping gene -actin (Mm 00607939_s1). Technical duplicates for each culture or muscle tissue and three or more independently derived cultures or biological samples were analyzed. The following primers were used in qPCR experiments: Alox5 sense, 5-CATGTTACCGC-

2630
Bmi1 overexpression in skeletal muscle regeneration | Di Foggia et al.
2008a). The researcher who performed all the analysis was blinded to the animal's genotype.
Enforced treadmill exercise and histological analysis. The dystrophic animals mdx;Pax7 Bmi1 , mdx control littermates and C57BL/6J-BL/10 mice (2-4 mo) were used in this experiment. Mice were injected with Tamoxifen. Starting from 3 d after the last injection, one mouse was put in each lane of the treadmill and forced to run at the speed of 12 m/min for 30 min, twice a week for 6 wk. This protocol was not exhausting and generally all the mice were able to perform the exercise. Whenever an avoiding behavior was observed, the mice were stimulated to run by a gentle touch of the tail or with the presence of a soft barrier at the end of the treadmill lane (De Luca, 2008b). The mice were sacrificed 6 d after the last treadmill session. One forelimb was embedded in OCT and frozen in isopentane cooled at the liquid nitrogen temperature. The other forelimb was dissected and the muscles snap-frozen in liquid nitrogen to perform RNA extraction and to check Bmi1 expression levels.
The frozen muscles were serially sectioned at 8 µm and processed for H&E staining and/or immunohistochemistry. The percentage of CNF was counted in the injured area of three different sections along the length of the leg for the same animal. The researcher who performed all the analysis was blinded to the animal's genotype.
Assessment of muscle strength. The dystrophic animals mdx;Pax7 Bmi1 , mdx control littermates and C57BL/6J-BL/10 mice (2-4 mo) performing the enforced treadmill exercise were used in this experiment. Before each treadmill session, a forelimb grip test was performed (De Luca, 2008c) to assess the course of the disease. The test was performed using a horizontal grip strength meter (Linton). In brief, each mouse was lifted by the tail toward the metal pull bar. Once the mouse had grasped the bar with the forelimbs, the mouse was pulled back in a horizontal plane. The protocol was repeated three consecutive times within 30 s. The grip strength was normalized against the body weight.
Cell transplantation experiments. Satellite cells isolated from mdx;STOPFloxBmi1 were treated as for the MT1 kd experiments. 6 d later, cells were harvested and 20,000 cells were transplanted into surgically exposed TA muscles of 3-wk-old mdx mice. After 27 d, the muscle tissues were fixed in 4% PFA for 6 h, left on in 30% sucrose, and embedded in OCT.
Serial sections (8 µm thick) were cut, mounted, and analyzed at the microscope for GFP positivity (two sections/mouse; three fields/section). The same settings were applied to take pictures of MT1 kd and MT1 SCr transplanted mice. For the analysis of the strong GFP positivity a defined threshold was applied to all pictures and the GFP +ve fibers were counted out of the total number of fibers/microscopic field. The percentage of strong GFP +ve fibers was calculated out of the total number of GFP +ve fibers/field. The researcher who performed all the analysis was blinded to the different treatments of the transplanted satellite cells.

Microarray analysis.
Three independently prepared satellite cell cultures isolated from STOPFloxBmi1 mice, treated with A-Cre infection and induced to differentiate for 2 d were analyzed by whole-genome Illumina platform mouse v2 as compared with the same cells treated with A-GFP. Integrity of starting total RNA was evaluated using capillary electrophoresis (Bioanalyzer 2100; Agilent Technologies) and quantified using a NanoDrop 1000 (NanoDrop). Aliquots of RNA (250 ng) samples were amplified according to the specifications of the Illumina TotalPrep RNA Amplification kit (Ambion). The cRNA samples obtained were applied to the MouseWG-6 v2.0 Expression BeadChip arrays (Illumina) and hybridized according to the manufacturer's specification. Each array on the BeadChip targets >45,200 transcripts using 3-µm beads bearing covalently attached 50-base oligonucleotide probes. The BeadChips were scanned with the BeadSrray system scanner (Illumina). The hybridization images signal intensity was extracted and background subtracted using Bead-Studio software (version 3.3.7 Illumina). The data discussed in this publication have been deposited in NCBI's Gene Expression Omnibus and are accessible through GEO Series accession no. GSE41878.
Modified biotin switch assay to detect SNO-MEF2c. The biotinswitch assay (BST) is a modified immunoblot analysis used to detect the levels of S-nitrosylated proteins. Here, we used a modified version of this assay. In brief, cell lysates were prepared from A-Cre-or A-GFP-treated mdx;STOPFloxBmi1 satellite cell cultures in HENS Buffer (Pierce S-Nitrosylation Western Blot kit; Thermo Fisher Scientific) and the BCA Protein Assay was used (Thermo Fisher Scientific) to determine the protein concentration. 300 µg of protein from each sample was mixed with 20 mM methyl methane thiosulfonate (MMTS) for 30 min at room temperature to block free thiol groups. After removing excess MMTS by acetone precipitation (1 h at 20°C), S-nitrosothiols were reduced to thiols with 20 mM sodium ascorbate. In the negative control, ultrapure water was used instead of sodium ascorbate. Labeling reagent containing tandem mass tag (TMT) was added to label the newly formed thiols. After acetone precipitation, samples were resuspended into HENS buffer. TMT labeled proteins were loaded onto spin columns with pre-added immobilized anti-TMT resin (Thermo Fisher Scientific) and eluted. The eluted proteins were concentrated using Pierce Concentrator (Thermo Fisher Scientific) and further analyzed by immunoblotting. Input was used for the control blot.
Pax7CreERT2 T2 tamoxifen injection. Tg:Pax7CreER T2 mice were crossed with R26R to optimize the injections protocol. 2-mo-old littermate Pax7CreER T2 ;R26R and R26R (as controls) were injected with Tamoxifen (Sigma-Aldrich) 20 mg/ml suspended in corn oil (Sigma-Aldrich; 4.5 mg/40 g of body weight) twice a day for five consecutive days. 5 d after the last injection the mice were culled and the efficiency of recombination was assessed. Pax7CreER T2 mice were bred with the STOPFloxBmi1 mice to overexpress Bmi1 in a cell type-specific fashion and the mdx mice for the experiments in a dystrophic model. In particular, the following compound mutant mice were generated: Pax7CreER T2 ;STOPFloxBmi1 (Pax7 Bmi1 ) and mdx;Pax7CreER T2 ;STOPFloxBmi1 (mdx;Pax7 Bmi1 ).
Acute muscle injury and histological analysis. Pax7 Bmi1 and littermate control mice (2-4 mo) were used in this experiment. Mice were injected with Tamoxifen (as described above) to induce Bmi1 overexpression in satellite cells. 10 d after the last injection, muscle freeze injury (FZ) was performed on these mice. 2% isoflurane inhalation was used to anesthetize the mice. The left hind limb was shaved and then wet with 70% ethanol to remove hair. To induce the damage to the Tibialis Anterior (TA) muscle, an incision of 5 mm was done in the skin above the muscle; then, a dry ice precooled probe was applied on the muscle for 10 s. The incision was closed by veterinary glue (VETBOND Tissue Adhesive). The mice were either sacrificed or used in a test for treadmill performance 10 d.a.i. Both the injured and contralateral hind limbs were collected. The injured leg was embedded in OCT and frozen in liquid nitrogen-precooled isopentane, whereas the contralateral muscles were snap-frozen in liquid nitrogen and processed for RNA extraction to check Bmi1 expression levels.
The frozen muscles were serially sectioned at 8 µm (with a step of 120 µm) and processed for H&E staining and/or immunohistochemistry. The percentage of CNF was counted in the injured area of three different sections along the length of the leg for each animal. In the same sections, the crosssectional areas (CSAs) of the fibers in regenerating and uninjured areas (as control) were measured using the MetaMorph software (Life Sciences Research Imaging Systems). The researcher who performed all the analysis was blinded to the animal's genotype.
Treadmill performance test after FZ injury. Tamoxifen-injected Pax7 Bmi1 and littermate controls mice (2-4 mo) were used in this experiment. 10 d after the muscle freeze injury mice performed a functional test running on a treadmill (LE8710M, Panlab Harvard Apparatus). 3 d before the test, mice were trained for each day for 5 min, at the speed of 10 m/min, with an inclination of 5°. The day of the test, the mice started running for 4 min (speed 10 m/min) and then every 4 min the speed was increased of 2 m/min until the mice were tired. Tiredness was defined as the inability of the mouse to go back to the treadmill belt despite five mechanical stimulations (De Luca,