Skip to main content

Exercise promotes satellite cell contribution to myofibers in a load-dependent manner



Satellite cells (SCs) are required for muscle repair following injury and are involved in muscle remodeling upon muscular contractions. Exercise stimulates SC accumulation and myonuclear accretion. To what extent exercise training at different mechanical loads drive SC contribution to myonuclei however is unknown.


By performing SC fate tracing experiments, we show that 8 weeks of voluntary wheel running increased SC contribution to myofibers in mouse plantar flexor muscles in a load-dependent, but fiber type-independent manner. Increased SC fusion however was not exclusively linked to muscle hypertrophy as wheel running without external load substantially increased SC fusion in the absence of fiber hypertrophy. Due to nuclear propagation, nuclear fluorescent fate tracing mouse models were inadequate to quantify SC contribution to myonuclei. Ultimately, by performing fate tracing at the DNA level, we show that SC contribution mirrors myonuclear accretion during exercise.


Collectively, mechanical load during exercise independently promotes SC contribution to existing myofibers. Also, due to propagation of nuclear fluorescent reporter proteins, our data warrant caution for the use of existing reporter mouse models for the quantitative evaluation of satellite cell contribution to myonuclei.


Muscle tissue has a remarkable ability to adapt its structure and microenvironment in response to diverse conditions such as injury, contractile activity, mechanical overload, and nutrients [1,2,3]. Satellite cells (SCs) are a unique population of stem cells that remain quiescent for the majority of time but upon activation can proliferate asymmetrically and fuse with myofibers [4]. It has been hypothesized that SCs play a supportive role in muscular remodeling by exercise, because SC accumulation as well as myonuclear accretion have been observed during episodes of hypertrophy [5,6,7,8]. Yet, the quantitative extent of SC contribution to myonuclei was not demonstrated by these studies. Genetic lineage experiments, in which Pax7+ cells and their derived fibers are fluorescently labeled, have recently shown that there is a considerable (up to 30%) contribution of SCs to uninjured adult mouse skeletal muscle fibers [9, 10]. To what extent SCs fuse with existing myofibers during exercise and whether SC fusion mirrors myonuclear accretion under different exercise intensities (with different load) remains to be elucidated.

While there is accumulating data showing SC are required for muscle hypertrophy [6, 8] the functional role of SCs in muscle endurance adaptations is much less understood [11]. Running exercise increases the number of myofibers containing an elevated number of SCs [12, 13], and SC number is directly proportional to running performance [14]. Nevertheless, SC number but not myonuclear number increased in several rat hind limb muscles after a period of free wheel running [15]. Moreover, no defects in intrinsic adaptations to endurance exercise (e.g., muscle vascularization, fiber type shift, or succinate dehydrogenase activity) were reported when SCs were ablated using the Pax7CreERT2-DTA model or X-ray irridation, arguing against the necessity of SCs for adaptations to endurance exercise [16, 17]. Consequently, it is currently unknown whether exercise training aimed to evoke endurance adaptions, but not muscle hypertrophy, stimulates SC contribution to myonuclei.

An excellent model to study the involvement of SCs in muscle adaptations to exercise is voluntary wheel running. Mice have a strong intrinsic drive to run large distances in their active phase [18], which keeps stress and mouse handling at a minimum. Furthermore, by applying resistance to the running wheel, muscle hypertrophy can be induced [19,20,21,22,23,24,25], and reducing wheel resistance (free running) leads to favorable endurance adaptations without pronounced effects on muscle hypertrophy [16, 26]. Thus, using this exercise model in combination with genetic lineage tracing experiments of SCs, we aimed to evaluate whether SC contribution to myonuclei is mediated by increasing load during exercise training and whether this manifests in SC-mediated myonuclear accretion.



All mice were housed in individually ventilated cages (3–4 littermates per cage) at standard housing conditions (22 °C, 12 h light/dark cycle), with ad libitum access to chow (KlibaNafag, diet #3436 and diet #3437) and water. Health status of all mouse lines was regularly monitored according to FELASA guidelines. All mouse lines were maintained on a C57BL/6J background: Pax7CreERT2 (Jackson labs stock no. 017763), RosamTmG (Jackson labs stock no. 007576), and RosaCAG-LSL-ntdTomato (Ai75D, Jackson labs stock no. 025106). RosamTmG and RosaCAG-LSL-ntdTomato mice were intercrossed with Pax7CreERT2 mice to generate Pax7CreERT2/+, RosamTmG/+ (Pax7mTmG/+) and Pax7CreERT2/+, and RosaCAG-LSL-ntdTomato/+ (Pax7nTom/+). All mice used in this study were female. At the age of 15 weeks, mice received tamoxifen via intraperitoneal injection at a dose of 1.0 mg/day for 5 days to induce Cre-mediated recombination, and a washout period of 7 days was allowed before experiments were initiated.

Experimental procedures

All animal procedures were approved by the veterinary office of the Canton of Zürich (license nr. ZH254-16). Mice were individually housed in open cages equipped with a running wheel device (TSE Systems, Bad Homburg vor der Höhe, Germany) for the duration of the intervention. Non-running control mice were also singly housed in cages of equal dimensions but without running wheels. The running wheel device continuously records wheel movements out of which total distance (km), speed (m.s−1), number of running bouts, and resistance on the wheel (N) were extracted. Additionally, to increase the force needed to rotate the wheel, resistance (0–100%) can be added. To calculate total external work, we used the equation W = Pt where W is work, P is power, and t is time. To calculate the power of the wheel at each braking resistance, we used the equation P = 2π × f × M where f is the angular frequency of the wheel, and M is the torque at a given braking resistance. TSE Systems provided a torque braking resistance curve. Mice were randomized in three groups: sedentary controls (no running), voluntary low load endurance wheel running where there is no extra resistance on the wheel (VRun), and voluntary high load resistance wheel running where extra resistance is placed on the wheel (VResRun). A representation of the study design is shown in Fig. 1a. Prior to the 8 weeks of voluntary wheel running, both groups VRun and VResRun were familiarized for 5 days to the running wheel without extra resistance. After familiarization, the VResRun group was subjected to a progressive weekly increase in resistance on the wheel: The load on the wheel was 50% in week 1, 60% in week 2, 70% in week 3, 72% in week 4, and 74% from weeks 5 to 8. The minimal distance which mice needed to run was set at 2 km per night. In case a mouse did not reach 2 km for 2 consecutive nights, resistance was slightly reduced for 3 nights after which it was raised again. All mice ran at 74% resistance during the last 2 weeks.

Fig. 1
figure 1

Running pattern and muscle hypertrophy are load-dependent. a Experimental set-up. b Average running distances and c average external work per night throughout the exercise protocol (n = 16 mice per group). d Soleus weight normalized to tibia length. e Representative cross-sections stained for laminin (red) and f quantification (n = 7–8 mice per group) of fiber cross-sectional area distribution in soleus muscle. g Representative pictures of single myofibers isolated from soleus and (h quantification of fiber volume of 11–16 single myofibers per mouse (n = 4 mice per group). i Muscle force production by grip strength and j time-to-fatigue in ex vivo stimulated muscle. Statistics, one-way ANOVA test with Tukey correction for multiple comparisons (d, h, i, j) or two-way ANOVA test with a Bonferroni post hoc test (b, c, f, k). (*p < 0.05; **p < 0.01; ***p < 0.001). Each dot represents a single mouse (d, i, j) or muscle fiber (h). Bar graphs (d, i, j) and line graphs (b, c, f, k) represent mean ± SEM (error bars). Violin plot (h). Scale bars, 50 μm (e) and 25 μm (g). Pax7mTmG/+ female mice on a C57BL/6J background were used for all analyses

Sample collection

Tissue collection was performed 24 h after the last exercise session after a 4-h fast. Mice were anaesthetized using Ketamine (80–100 mg/kg), Xylazine (10–15 mg/kg), and Acepromazine (2–5 mg/kg) via intraperitoneal injection 5 min before sacrifice. The depth of anesthesia was confirmed by testing pedal withdrawal reflex. Subsequently, the m. soleus (SOL), m. plantaris (PLT), m. gastrocnemius (GAS), m. tibialis anterior (TA), and m. extensor digitorum longus (EDL) were harvested, weighed, and frozen in OCT embedding matrix (CellPath) in liquid nitrogen-cooled isopentane for histochemical analysis or fixed in 4% PFA for single myofiber isolation. Tibia length was assessed by a digital caliper. After sample collection, animals were euthanized, and major bleeding was induced to confirm death.


Frozen sections (10 μm) of muscle embedded in OCT of mid-belly level were made using a cryostat (Leica CM 1950) and collected on Superfrost Ultra Plus slides (Thermo Scientific, Zug, Switzerland). The following antibodies were used: anti-MyHC-I (1:50 dilution, BA-F8 from hybridoma, Iowa City, IA, USA), anti-MyHC-IIa (1:200 dilution, SC-71 from hybridoma), anti-MyHC-IIb (1:100 dilution, BF-F3 from hybridoma), anti-laminin (1:200 dilution, PA1-16730, ThermoFisher), anti-PCM1 (1:1000 dilution, HPA023370, Sigma, Buchs, Switzerland), eMHC (1:20 dilution, F1.652 from hybridoma), Alexa Fluor 488 goat anti-mouse IgG2B (1:250 dilution, ThermoFisher), Alexa Fluor 350 goat anti-mouse IgG1 (1:250 dilution, A-21120, ThermoFisher), Alexa Fluor 568 goat anti-mouse IgM (1:250 dilution, A-21043, ThermoFisher), Alexa Fluor 647 goat anti-rabbit IgG (1:250 dilution, A-21244, ThermoFisher), Alexa Fluor 568 goat anti-mouse IgG1 (1:400 dilution, A-21124, ThermoFisher), and Alexa Fluor 647 conjugated wheat germ agglutinin (WGA, 1:50 dilution, W32466, ThermoFisher).

Immunohistochemical muscle fiber staining was conducted as previously described [27]. In short, sections were dried and washed for 5 min in PBS supplemented with 0.05% Triton-X-100 (PBST) and subsequently blocked for 60 min in PBST + 10% goat serum (16210064, ThermoFisher). Afterwards a primary antibody cocktail diluted in PBST + 10% goat serum was applied for 120 min against MyHC-I, MyHC-IIa, MyHC-IIb, and laminin. After washing 3 times for 5 min, a secondary antibody cocktail, diluted in PBST + 10% goat serum, was applied for 60 min. Slides were mounted after a 3 x 5 min wash and sealed with glass cover slips. A similar protocol was used for eMHC. Sections were dried and washed for 5 min in PBST and subsequently blocked for 60 min in PBST + 10% goat serum (16210064, ThermoFisher). Sections were stained overnight at 4 °C with eMHC (F1.652) and laminin diluted in PBST + 10% goat serum. Next day after washing 3 times for 5 min, anti-mouse and anti-rabbit secondary antibodies were diluted in PBST + 10% goat serum for 1 h. Slides were mounted. For myonuclear staining, samples were stained with anti-PCM1 as previously described [28]. Slides were pre-incubated in 2% BSA in PBS. Sections were stained with a rabbit polyclonal antibody against PCM1 in staining solution (5% BSA in PBS, 0.2% Triton-X-100) overnight at 4 °C. The next day, the sections were washed three times for 5 min in PBS and stained with an anti-rabbit secondary antibody in 2% BSA in PBS for 1 h. Sections were washed three times for 10 min with PBS. For detection of the fiber boundaries, the sections were stained with WGA in PBS for 30 min. Nuclei were co-stained using Hoechst dye 33342 (dilution 1:5000, H3570, ThermoFisher).

Images from sections were captured at × 10 using an epifluorescent microscope (Zeiss Axio observer Z.1, Zeiss, Oberkochen, Germany). Composite images were stitched together using the tiles module in the ZEN 2011 imaging software (Zeiss). All images were captured at the same exposure time. Approximately 100 (SOL, PLT) and 200 (GAS) muscle fibers, divided in different regions over the whole muscle section, were analyzed for PCM1+/Hoechst+ nuclei with the Zen software (ZEN 2011 imaging software, Zeiss).

Satellite cell contribution to myofibers in Pax7mTmG/+ mice was quantified based on a previously described method of myofiber masking [29]. Tilescan images were processed using an Image-J plug-in to define myofiber cross-sectional area (CSA) [30], followed by a manual refinement of the CSA outlines. An overlay of those outlines on the corresponding mGFP image was performed, after which the data was compiled in Excel, and a threshold was set for the mean gray value of mGFP. A frozen section isolated from a vehicle-treated mouse (negative control sample) was used to determine the threshold and to ensure the elimination of autofluorescence. mGFP images were overlaid with the muscle fiber type staining to identify mGFP+ oxidative (MyHCI) and glycolytic (MyHCIIa, MyHCIIb, and MyHCIIx) fibers.

Single myofiber isolation

For single myofiber analysis, the SOL from the Pax7nTom/+ mice was harvested, fixed in 4% paraformaldehyde/PBS for 1 h at 4 °C and transferred to PBS containing 2% horse serum. To release individual myofibers, single fibers were mechanically teased apart, strained, and washed with PBS. Isolated single myofibers were stained with Hoechst for nuclear visualization and gently mounted on a glass coverslip. Images were captured at × 10 using a confocal microscope (Olympus olympus fluoview FV 3000, Olympus, Hamburg, Germany). Eleven to sixteen fibers per mouse from random areas from the muscle were analyzed. Total number of nTom+ nuclei and Hoechst+ nuclei along 500 μm fiber length were quantified using a spot detection algorithm developed in Imaris (Bitplane). Myofiber volume was approximated as the volume of a cylinder using the average radius measured along the fiber length in Imaris. Fibers isolated from vehicle-treated mice were used to set a threshold to ensure the elimination of autofluorescence.

Cell culture


Freshly isolated primary myoblasts were cultured in growth medium, which contained a 1:1 ratio DMEM (ThermoFisher Scientific, 12320032) and Ham’s F-10 nutrient mix (ThermoFisher Scientific, 22390058) supplemented 20% FBS (Thermo Scientific, 10270106) and 10 ng/ml basic-FGF (ThermoFisher Scientific, PHG0266). Cells were differentiated in low-serum differentiation medium (Thermo Scientific, 41965039), supplemented with 2% horse serum (Thermo Scientific, 16050-122). Cells were routinely cultured in 21% O2 and 5% CO2. Every other day, differentiation medium was changed until myotubes reached full differentiation. All media were supplemented with 100 units/ml penicillin and 100 μg/ml streptomycin.

Isolation of SCs/myoblasts

Isolation of SCs/myoblasts was done as previously described [31]. Briefly, muscle tissue was digested in Hank’s Balanced Salt Solution (HBSS, Thermo Fisher, 24020117) supplemented with 1.5% bovine serum albumin (BSA) and 2 mg/ml collagenase type II (ThermoFisher Scientific, 17101015) for 1 h at 37 °C. After centrifugation, the cell pellet was then filtered using 100- and 40-μm cell strainers, and a heterogeneous cell population was purified by FACS sorting or by serial preplating. For FACS, SCs were sorted based on positive alpha 7-integrin (1:100, R&D Systems FAB3518P) and absence of Sca1 (1:1000, Bio-Legend 122511), CD31 (1:1000, Bio-Legend, 102413), and CD45 (1:1000, Bio-Legend 103121).

Generation of nuclear H2B-GFP (H2B-nGFP) myoblasts

Lentivirus was generated with PEG-it Virus Precipitation Solution (System Biosciences, Palo Alto, CA, USA no LV810A-1) using a LV H2B-GFP plasmid (Addgene no 25999) in HEK 293 cells. Myoblasts were transfected with 20 μl of concentrated lentivirus for 24 h using growth medium supplemented with 4 μg/ml polybrene transfection reagent (TR1003-G, Sigma-Aldrich).

Co-culturing experiments

nTom+ myoblasts, together with committed WT myoblasts, were co-cultured at 1:4 ratio (nTom+/WT). nTom+ myoblasts were co-cultured with myoblasts stably expressing H2B-nGFP at a ratio of 1:1 ratio (nTom+/H2B-nGFP+). Two hours after seeding, growth medium was changed to differentiation medium. Every other day, differentiation medium was changed until myotubes reached full differentiation, and immunofluorescent images were taken every 24 h via an epifluorescent microscope, equipped with a cell incubator (Nikon eclipse Ti2, Nikon, Amstelveel, Netherlands). Fluorescent pictures were captured every 15 min to generate a time-lapse video (Windows movie maker, Windows, Redmond, WA, US).

SC fate tracing at the DNA level using RT-PCR for Cre-mediated recombination in bulk muscle and sorted myonuclei

To perform SC fate tracing at the DNA level, genomic DNA was isolated from muscle or myonuclei using QIAamp DNA Micro and Mini Kit (Qiagen, Hilden, Germany). Myonuclei were isolated as previously described [32]. Nuclei were isolated with EZ PREP buffer (Sigma, NUC-101). Frozen soleus muscles (10 mg) were homogenized using a glass dounce tissue grinder in 2 ml of ice-cold EZ PREP and incubated on ice for 5 min, with an additional 2 ml of ice-cold EZ PREP and filtered through a 100-μm cell strainer. Nuclei were centrifuged at 500×g for 5 min at 4 °C, washed with 4 ml ice-cold EZ PREP, and incubated on ice for 5 min. After centrifugation, the nuclei were stained for 45 min with an antibody against PCM-1 (1:1500, HPA023370, Sigma), followed by an Alexa647-anti-rabbit secondary antibody for 30 min (1:1500) and Hoechst 33342 (1:5000), and filtered through a 35-μm cell strainer. Hoechst+ and Alexa647+ myonuclei were directly sorted into RLTplus buffer (Qiagen). A final amount of 10,000–20,000 nuclei was used for genomic DNA isolation. A SYBR Green-based master mix (ThermoFisher Scientific, A25778) was applied for RT-PCR analysis. Recombination rates were calculated from the relative expression of recombined levels normalized to an internal control. The delta–delta CT method was used to normalize the data. The percentage recombination was calculated relative to positive controls. To generate a positive control, we crossed HSACreERT2 mice [33] with RosamTmG mice to generate HSACreERT2/+; RosamTmG/+ mice. In these mice, tamoxifen injection leads to the excision of mTomato from the DNA in myonuclei and the subsequent expression of mGFP. Additionally, we sorted SCs from Pax7nTom/+ mice. The following primers were used to detect recombination (see also Fig. 5a, b): RosamTmG, recombined amplicon: Fw5'- GGGCTCGACACTAGTGAACC-3' Rv 5'-GGTGATGATCCGGAACCCTT-3', Internal control: Fw 5'-AGCGAACGGACAGGAGAATG-3', Rv 5'-ACTTGTGGCCGTTTACGTCG-3', RosaCAG-LSL-ntdTomato, recombined amplicon: Fw 5'-TTATTGTGCTGTCTCATCATTTTGG-3' Rv 5'-CTTCGCGAGACCGGTAATCT-3', Internal control: Fw5'- CTGTTCCTGTACGGCATGG-3' Rv 5'-GGCATTAAAGCAGCGTATCC-3'.

Grip strength

Grip strength of the limbs was measured using a force tension apparatus (grip strength meter, Bioseb, Vitrolles, France) after 8 weeks of voluntary wheel running or no running. Mice were lifted by the tail and were made to hold a metal grid with all limbs. Total limb maximal grip strength was registered in gram during three consecutive attempts, and the result was set as the average of the attempts.

Ex vivo muscle force assay

The SOL was carefully dissected and placed in Krebs–Henseleit buffer (120 mM NaCl, 4.8 mM KCl, 25 mM NaHCO3, 2.5 mM CaCl2, 1.2 mM KH2PO4, 2 mM MgSO4) supplemented with 25 mM glucose at 37 °C and bubbled with 95% O2–5% CO2 to stabilize the pH at 7.4. The distal tendon was securely connected to a fixed bottom plate and the proximal tendon to the arm of a servomotor (800A in vitro muscle apparatus, Aurora Scientific). Muscle length was adjusted until a single stimulus pulse elicited maximum force during a twitch (optimal length, Lo) under isometric conditions. After 15 min of equilibration, a fatigue protocol was started. To fatigue the muscle, tetanic contractions of 60 Hz were produced, once every 2 s, for a total of 300 s. Time-to-fatigue at 75%, an assessment for the time it takes for a muscle to fall to 75% of its initial force, was calculated for every mouse.

Statistical analysis

Statistical significance was determined by a one-way ANOVA with Tukey correction for multiple comparisons or two-way ANOVA with Bonferroni correction with Graphpad Prism 8.2. Data are presented as means ± SEM, and values of p < 0.05 were considered statistically significant. The degree of significance is depicted as follows: *p < 0.05, **p < 0.01, and ***p < 0.001.


Running pattern, hypertrophy, and muscle performance in VRun vs. VResRun

To study how increased load during exercise affects contribution of SCs to myofibers, we used 8 weeks of voluntary low load–endurance type–running (VRun) and voluntary high load–resistance type–running (VResRun) as exercise models (Fig. 1a). A progressive increase in load reduced running distance (Fig. 1b), running speed, and total daily running time without affecting the number of running bouts from week 3 onwards (Figure S1A-C). In contrast, external work was higher over the 8-week time period (Fig. 1c). In agreement with previous observations [23, 24, 34], extra resistance on the wheel caused an increase in the weight of m. soleus (SOL), while unloaded running did not change SOL weight (Fig. 1d). Body weight, lean mass, and fat mass were not affected by exercise (Table S1). Further analysis on cross-sectional area (CSA) confirmed that VResRun, but not VRun, induced fiber hypertrophy in SOL (Fig. 1e, f). This was exclusively caused by an increase in the CSA of MyHCIIa muscle fibers, while MyHCI fiber CSA remained unaffected (Figure S2A-E). Finally, in accordance with the data on whole muscle, VResRun increased the volume of isolated single myofibers from SOL, while VRun did not affect myofiber volume (Fig. 1g, h).

Despite unchanged muscle weight in PLT (Table S1), we measured increased fiber CSA in PLT upon VResRun as well as VRun. CSA of MyHCIIa fibers increased both after VRun and VResRun compared to control, whereas MyHCIIb-IIx fiber CSA only increased after VResRun (Figure S2F-J). The discrepancy between muscle mass and CSA in the PLT could be due to difficulties in anatomically dissecting the muscle [35]. Muscle mass of m. gastrocnemius (GAS), m. tibialis anterior (TA), and m. extensor digitorum longus (EDL) remained unaffected by VRun and VResRun (Table S1). In addition, no change in fiber CSA was found for GAS, TA, and EDL (Figure S2K,L,M). We also performed immunofluorescent stainings for embryonic myosin heavy chain (eMHC) on VRun and VResRun SOL muscle samples, using glycerol injected muscle as a positive control. We however could not detect eMHC+ fibers in VRun nor VResRun whereas extensive eMHC+ fibers were detected in glycerol damaged muscle (Figure S3A). This shows that at the time of muscle sampling (8 weeks of running), no active regeneration was ongoing. Subsequently, we also evaluated the number of centrally nucleated fibers. We observed similar fractions of nucleated fibers in both VRun as well as VResRun, whereas negligible centrally nucleated fibers were observed under control conditions (Figure S3B). A discrepancy between eMHC+ and centrally nucleated fibers has been observed before upon synergist ablation and tenotomy [36]. Potentially, this can be the sign of ongoing fiber splitting or other myofiber remodeling processes which do not necessarily requires satellite cells [36].

Finally, to analyze whether VResRun resulted in different functional adaptations when compared to VRun, we tested functional strength (grip strength) as well as muscle endurance capacity. After 8 weeks, grip strength was increased in VResRun (Fig. 1i), showing that VResRun increased functional strength. In addition, time-to-fatigue at 75%, an assessment for the time it takes for an ex vivo stimulated muscle to fall to 75% of its initial force, was higher after VRun and VResRun compared to control (Fig. 1j), indicating that free running improved intrinsic muscle endurance capacity. Collectively, these data confirm that VResRun, but not VRun, induces an increase in muscle mass in both SOL and PLT, which was induced by hypertrophy of the glycolytic muscle fibers. On the other hand, VRun leads to more pronounced endurance adaptations.

Exercise promotes SC contribution to myofibers in a load-dependent manner

To investigate whether voluntary wheel running affects SC contribution to various hind limb muscles, and whether this is load-dependent, we genetically labeled SCs and performed fate-tracing experiments into myonuclei. To do so, we used a mouse model which was previously used to study SC involvement in muscle regeneration [37] and muscle homeostasis [9]. Briefly, we intercrossed Pax7CreERT2/CreERT2 mice, which efficiently induces recombination in SCs, with mice carrying a RosamTmG/mTmG reporter [38]. The resulting Pax7CreERT2/+; RosamTmG/+ mice (Pax7mTmG/+) ubiquitously express membrane-bound Tomato (mTomato), but following tamoxifen injection, Cre-mediated recombination results in the excision of mTomato and expression of membrane-bound GFP (mGFP) in SCs. In these mice, the presence of GFP+ myofibers indicates that at least one SC fused with those myofibers. In our experiments, we labeled SCs by 5 consecutive daily IP tamoxifen injections and initiated the exercise protocol 1 week later. Subsequently, hind limb muscles were harvested after 8 weeks of VRun or VResRun or no running (control) (Fig. 2a), and the contribution of SCs to myofibers was evaluated. Importantly, we opted to use 15 weeks old mice for all our experiments since this is an age where postnatal muscle growth has ceased and SC contribution to myofibers is restricted to maintenance of muscle homeostasis [10, 36, 39,40,41]. Indeed, in agreement with previous work [9, 10], we found that in the non-exercising control mice, 5–20% of all myofibers were GFP+, with the highest contribution of SCs in the SOL, confirming that SC fusion occurs in hind limb muscles of sedentary mice (Fig. 2b–d). Interestingly, running increased the SC contribution to SOL, PLT, and GAS, whereas we did not observe a clear increase in SC contribution to fibers upon exercise in TA, nor EDL (Fig. 2b–d), confirming previous observations from our group that wheel running specifically recruits the hind limb plantar flexors [42]. Moreover, SC contribution to the SOL was load-dependent, as VResRun further increased the number of GFP+ myofibers compared to VRun alone, resulting in 80% GFP+ myofibers vs. 60% in VRun condition. In PLT and GAS, we also observed a 38% increase in GFP+ myofibers upon VResRun, while the increase by VRun did not reach statistical significance (Fig. 2b–d).

Fig. 2
figure 2

Exercise promotes SC contribution to myofibers in a load-dependent manner. aPax7CreERT2 mice were crossed with RosamTmG mice to genetically trace SC fusion to myofibers. Upon tamoxifen treatment SCs will express mGFP. mGFP+ myofibers indicate fusion of SCs. b A cryosection from a vehicle-treated mouse (negative control sample) was used to determine the threshold and to ensure the elimination of autofluorescence. c Representative cross-sections of soleus (SOL), plantaris (PLT), gastrocnemius (GAS), tibialis anterior (TA), and extensor digitorum longus (EDL) muscle showing GFP+ myofibers to which SCs contributed. d Quantification of GFP+ myofibers. Statistics, one-way ANOVA test with Tukey correction for multiple comparisons (*p < 0.05; **p < 0.01; ***p < 0.001). Each dot represents a single mouse. Bar graphs represent mean ± SEM (error bars). Scale bar, 100 μm

To investigate whether the increase in SC contribution to myofibers by exercise was different between fiber types, we made sequential cryosections of SOL, overlaid the GFP signal with fiber type immunostainings, and assessed the relative amount of GFP+ fibers to MyHC type (Figure S4A). VResRun induced the largest increase in GFP+ fibers resulting in ~ 80% of all MyHCI and MyHCIIa fibers becoming GFP positive, while this was only ~ 60% in VRun, indicating that 8 weeks of increased load promoted SC contribution to myofibers but that there was no fiber type preference (Figure S4B,C). Moreover, in agreement with previous data [43], both running modalities led to a shift towards a more oxidative fiber type (MyHCI) with the largest shift in the VRun group (Figure S4D). Thus, because we observed similar contribution of SCs to MyHCI versus MyHCIIa, relative GFP contribution after exercise mirrored the shift in fiber type (Figure S4E). Altogether, we show that exercise drives SC contribution to muscle fibers in a load-dependent manner but independent of muscle fiber type.

Nuclear propagation renders nuclear reporter mice inadequate to assess SC contribution to myonuclei

Our data from the Pax7mTmG/+ mice indicated that exercise strongly increases SC contribution to myofibers. Nevertheless, this model does not allow us to assess the exact number of SCs that fuse with the myofiber since the GFP signal diffuses throughout the sarcolemma, nor does it provide information on myonuclear accretion upon exercise training. Therefore, we generated a nuclear reporter mouse by crossing the aforementioned Pax7CreERT2 with a RosaCAG-LSL-ntdTomato mouse. The RosaCAG-LSL-ntdTomato mouse expresses a loxP-flanked stop codon before the nuclear tdTomato fluorescent red (nTom) reporter. Upon tamoxifen injection, the nuclear stop codon is excised to indelibly label Pax7+ SCs and their derivative nuclei with nTom (Fig. 3a). We used the Pax7CreERT2/+; RosaCAG-LSL-ntdTomato/+ (Pax7nTom/+) mice to assess myonuclear accretion and SC contribution to myonuclei.

Fig. 3
figure 3

Nuclear reporter mice show high SC contribution during exercise. aPax7CreERT2 mice were crossed with RosaCAG-LSL-ntdTomato mice to genetically trace SC-derived nuclei. Upon tamoxifen treatment, SCs and SC-derived myonuclei will express nTomato (nTom). b Quantification of Hoechst+ and c nTom+ nuclei per 100 μm fiber length (11–16 single myofibers per mouse, n = 4 mice per group). d Represented pictures of single myofibers isolated from soleus. Statistics, one-way ANOVA test with Tukey correction for multiple comparisons (*p < 0.05; **p < 0.01; ***p < 0.001). Violin plot is shown (b, c). Each dot represents a single mouse. Scale bar, 25 μm

First, by isolating single muscle fibers from SOL (Fig. 3a) after 8 weeks of exercise training, we found that VResRun, and to a lesser extent VRun, increased the number of Hoechst+ nuclei on single fibers compared to control by respectively 16 and 14% (p < 0.05 and p = 0.05, Fig. 3b). Subsequently, we also analyzed myonuclear content on cryosections using a myonuclear-specific PCM1 staining. Stainings on SOL sections corroborated these findings (Figure S5A-B). Additionally, we also found a load-dependent increase in myonuclear content in PLT and GAS (Figure S5A-B). Intriguingly, compared to myonuclear (Hoechst+) accretion, we observed a much larger increase in nTom+ nuclei after both VRun (+ 1.8 Hoechst+ nuclei vs. + 7.0 nTom+ nuclei per 100 μm fiber length) and VResRun (+ 2.0 Hoechst+ nuclei vs. + 8.3 nTom+ nuclei per 100 μm fiber length) on single fibers (Fig. 3b–d), suggesting (1) high turnover of myonuclei during exercise or (2) nuclear reporter mice are inadequate to measure SC fusion in myocytes due to nuclear travelling of nTom, a process recently termed “propagation” by Taylor-Weiner et al. [44]. To investigate the latter hypothesis, we performed a series of in vitro co-culturing experiments. In a first experiment, we cultured nTom+ myoblasts together with myoblasts which stably expressed GFP in the nucleus (H2B-nGFP) at a 1:1 ratio and differentiated them to myotubes for 144 h (Fig. 4a). In this experiment, protein propagation between nuclei in multinucleated myotubes would lead to nTom+/H2B-nGFP+ myonuclei. Time-course imaging revealed that once H2B-nGFP+ and nTom+ myoblast started to fuse, double positive nuclei (yellow) became apparent, resulting in 45% of all nuclei that turned double positive for H2B-nGFP and nTom at full differentiation (Fig. 4b, c), showing that there is reporter protein propagation. Since only limited number of SCs contribute to myonuclei in myofibers in vivo, we performed a second experiment, where we cultured nTom+ myoblasts together with committed myocytes which do not express a reporter protein (WT) at a ratio of 1:4 (nTom+/WT) and further differentiated them to myotubes for 144 h (Fig. 4d). Even though myoblasts completely stop proliferating 48 h after initiation of differentiation (data not shown), sequential imaging showed that an increasing fraction of nuclei became positive for nTom over the 144-h follow-up period (Fig. 4e–g). In fact, while less than 30% of nuclei were nTom+ 48 h after the start of the experiment, almost 60% of all nuclei were nTom+ at 144 h (Fig. 4e). Altogether, this data shows that as soon as SC fuse to existing fibers and contribute to the myonuclear pool, during in vitro differentiation or in vivo following exercise, proteins from the newly acquired myonuclei propagate to old residing myonuclei within the same myofiber. In addition, we also observed that the original nTom+ myoblasts expressed higher levels of nTom when compared to the nTom levels in the differentiated cells, suggesting that nTom protein was diluted between nuclei (Fig. 4f). Finally, time-lapse microscopy imaging showed that propagation occurs at very fast rates, as nTom was transported towards neighboring nTom myonuclei within 15 min of myoblast fusion (Video S1). In conclusion, it is not possible to determine the exact number of SCs that contribute to myonuclei using common lineage tracing reporter models.

Fig. 4
figure 4

Fluorescent nuclear protein propagation renders reporter mice inadequate to assess SC fusion in myofibers. a Nuclear H2B-nGFP expressing and nTom+ myoblast were seeded at 1:1 ratio and differentiated for 144 h. b Double positive nuclei for H2B-nGFP+ and nTom+ (marked in yellow) were assessed over time in relation to total nuclei (Hoechst+). c Representative pictures of H2B-nGFP+ and nTom+ myoblast co-cultures over time. d nTom+ and committed WT myoblast were seeded at a 1:4 ratio and differentiated for 144 h. e Double positive nuclei for nTom+ and Hoechst+ were assessed over time in relation to total nuclei (Hoechst+). f nTom intensity was measured 1 h after seeding (T0) and upon full differentiation (T144). Data is represented as percentage of nTom+ nuclei per intensity range in relation to the total nTom+ nuclei. g Representative pictures of WT and nTom+ myoblast co-cultures over time. Statistics, one-way ANOVA test with Tukey correction for multiple comparisons (*p < 0.05; **p < 0.01; ***p < 0.001). Bar graphs represent mean ± SEM (error bars). Scale bar, 50 μm

Additional file 6: Video S1. In vitro nuclear propagation. Time-lapse epifluorescence microscopy of a fusing nTom+ myoblast with WT myoblasts (nTom-). White arrow shows propagation of nTom. Scale bar, 50 μm.

SC fate tracing at the DNA level confirms load-dependent SC fusion during exercise

Since our data shows that counting the number of nTom+ nuclei in Pax7nTom/+ mice does not allow to evaluate the exact contribution of SCs to myonuclei during exercise, we used a RT-PCR on genomic DNA to measure Cre-mediated recombination. This allows to trace the fate of SC at the DNA level. To this end, we sorted Hoechst+/PCM1+ myonuclei from Pax7nTom/+ control, VRun, and VResRun mice (SOL muscle) (Fig. 5b) and isolated DNA from these myonuclei. We designed primers specific for the recombined sequence and an internal control (housekeeping gene) (Fig. 5a). We found almost 20% recombination after 8 weeks of VResRun, indicating that 20% of all myonuclei were recombined, while only 12% and 4% recombination was observed in VRun and control respectively (Fig. 5c). To further substantiate these data, we extracted genomic DNA from bulk muscle (SOL) from control, VRun, and VResRun Pax7mTmG/+ mice and observed a similar load-dependent increase in SC contribution to nuclei upon exercise (Fig. 5d, e). Collectively, these observations show that mechanical load during exercise independently promotes SC contribution to existing myofibers which manifests in myonuclear accretion.

Fig. 5
figure 5

Myonuclear-specific Cre-mediated recombination. a Schematic diagram of the RosaCAG-lsl-ntdTomato construct before and after Cre-mediated recombination. Representation of primers used for quantitative RT-PCR analysis showing amplicon of recombinant DNA and amplicon of internal control for normalization. b Representative flow cytometry plots depicting the gating strategy for sorting of myonuclei. Isolated nuclei were analyzed by side scatter (SSC) and Hoechst. Hoechst+ PCM1+ single nuclei were sorted for further analysis. Purity of sorted myonuclei is shown in the right panel. c RT-PCR performed on sorted Hoechst+ PCM1+ nuclei. DNA isolated from satellite cells from Pax7CreERT2/+; RosaCAG-lsl-ntdTomato/+ mice with the presence of one recombined allele in each satellite cell served as a 100% reference (positive control). d Schematic diagram of the RosamTmG construct before and after Cre-mediated recombination. Representation of primers used for quantitative RT-PCR analysis showing amplicon of recombinant DNA and amplicon of internal control for normalization. e RT-PCR performed on genomic DNA isolated from bulk muscle. DNA isolated from tamoxifen injected HSAiCre/+; RosamTmG/+ mice with the presence of one recombined allele in each myonucleus, served as a 100% reference (positive control). Recombination rates were calculated from the relative expression of recombined levels normalized for internal control. Statistics, one-way ANOVA test with Tukey correction for multiple comparisons (*p < 0.05; **p < 0.01; ***p < 0.001). Each dot represents a single mouse. Bar graphs represent mean ± SEM (error bars)


Satellite cells are required for muscle regeneration and in addition contribute to muscle homeostasis under sedentary conditions [9, 10], but how they contribute to training adaptations is still poorly understood. Using two SC fate tracing mouse models, we show that exercise promotes the contribution of SCs to myofibers in a load-dependant manner. Moreover, we found that SC contribution to myofibers is present to a similar extent in oxidative and glycolytic muscle fibers.

Mice which ran against an external load (VResRun) showed higher SC contribution when compared to free runners (VRun). Notwithstanding, VRun extensively increased SC fusion without stimulating fiber hypertrophy. In agreement with other studies [34], we found clear hypertrophy after resistance running but no consistent increases in fiber size/volume after voluntary running: SOL wet weight did not increase, fiber cross sectional area was not changed, and we also did not observe an increase in SOL single fiber volume after VRun. Some studies have reported hypertrophy during short term (4–8 weeks) voluntary free wheel running [17, 23] while others failed to do so [34]. Likely, the effects of short-term voluntary free wheel running are minor and potentially only reach significance when corrected for body weight [23]. Against this background, we found significant SC contribution to muscle where we did not observe hypertrophy. More specifically, we report SC contribution (in both MHCI as well as MHCIIa fibers) of SOL upon VRun. We also show SC contribution to MHCI fibers in SOL upon VResRun, which did not exhibit enhanced fiber cross sectional area. Finally, we also found increased myonuclear density in GAS, despite lack of increased muscle weight. Thus, our data shows that SC fusion can be uncoupled of muscle hypertrophy under specific conditions.

The potential signals promoting exercise-induced SC-to-myonuclear contribution remain elusive and have been hypothesized to be both local and systemic. When rodents walk or run, most of the work is performed by the plantar flexors. The dorsi flexors undergo much less mechanical load since they only keep the foot up in the swing phase [45]. In this respect, recent work from our group has shown that acute VResRun increases mTORC1 signaling exclusively in SOL, PLT, and the oxidative part of the GAS but not in TA [42]. Moreover, muscle hypertrophy is predominantly observed in plantar flexors upon resistance running, as opposed to dorsi flexors [24, 34, 46]. Using Pax7CreERT2/+; RosamTmG/+ mice, we show that voluntary running strongly increases SC fusion in the plantar flexors SOL and PLT, partially in the GAS, and not in EDL or TA dorsiflexors. Furthermore, adding resistance to the wheel further increased SC fusion, indicating that the extra muscular work independently promotes SC fusion. Detailed analyses of myonuclear accretion via immunostainings of PCM1 and nuclei count on single fibers when compared to satellite cell tracing at the DNA level confirm that load-dependent augmentation of SC fusion corroborates very well with increases in myonuclear number. Thus, our data suggests that local rather than systemic signals regulate exercise-induced SC fusion and subsequent myonuclear accretion.

Our study did not focus on elucidating whether or not satellite cell contribution to myonuclei is driven by muscle damage or a consequence of exercise-dependent signals. We could not detect eMHC+ fibers in VRun nor VResRun. This shows that at the time of muscle sampling (8 weeks of running), no active regeneration was ongoing. Studies have shown that high intensity exercise (predominantly eccentric exercise) leads to muscle damage which is associated with satellite cell activation [47,48,49]. However, recent work has also shown that local injury occurs in exercising muscle but that this might not be responsible for widespread SC activation upon resistance type exercise [50]. Thus, the mechanisms that promote SC activation, proliferation, and fusion upon exercise/load are likely to be different when compared to injury-induced regeneration. Potentially, even temporal differences between different stimuli exist [8]. Since myonuclear accumulation occurs gradually during training [8], it is unlikely that it is exclusively driven by damage. Future studies will need to reveal the exact contribution of both mechanisms to exercise-induced myonuclear accumulation.

Whether SCs are required for muscle hypertrophy is a heavily debated topic [51] and likely depends on age, genetic mutation, models used to induce hypertrophy, and time of sampling after the applied stimulus [6, 8, 19,20,21, 25, 52, 53]. Moreover, impairing SC contribution by the frequently used Pax7CreER-DTA model might evoke undesirable adaptations such as impairments of proprioception and reduced voluntary wheel running, which do not reflect physiological circumstances [16]. Although we could not causally link SC fusion with hypertrophy in this study since we used reporter models to trace SC contributions under physiological circumstances, our tracing approach allows us to evaluate whether SC fusion occurs upon exercise and whether concomitant myonuclear accretion occurs in the absence or presence of hypertrophy. VResRun induced a 23% increase in SOL wet weight, which was exclusively caused by the increase in CSA of MyHCIIa fibers, while no hypertrophy was observed in VRun in any fiber type. This is remarkable given that even though SC fusion was higher in VResRun, we observed similar numbers of GFP positive MyHCI versus MyHCIIa fibers in both VRun and VResrun. Thus, the lack of hypertrophy in the heavily SC fused MyHCI fibers promotes the hypothesis that myonuclear accretion might also play a role in endurance adaptations to interval training [11, 54] and that myonuclear accretion is not merely a way to withstand to expansion of the myonuclear domain [55].

The Pax7mTmG/+ mice only provide a rough estimation of SC contribution to myonuclei under exercise conditions since the mGFP diffuses throughout the sarcolemma. Therefore, to get a quantitative insight into the extent of SC contribution to myonuclei, we used a nuclear reporter mouse in which the fluorescent reporter (nTom) is retained exclusively within the Pax7+ cells and all SC-derived myonuclei. Strikingly, we found that nearly all myonuclei became nTom+ after 8 weeks of VRun and VResRun, despite only a modest 14% and 16% increase in myonuclear number by VRun and VResRun respectively. During the preparation of this manuscript, an extensive report provided direct evidence of fluorescent protein traveling between nuclei of primary differentiated myotubes, a phenomenon termed nuclear propagation [44]. In concordance, spatial distribution of specific transcription factors, which may exit the myonucleus and activate receptors in other regions of the myofiber has been reported [56], suggesting nuclear transport of various proteins between nuclei in myofibers. We confirmed propagation in our model by a series of co-culturing experiments using H2B-nGFP and nTom expressing myoblasts. After 6 days of differentiation, when myotubes reached full differentiation, nearly 60% of all nuclei became nTom+, while WT and nTom were initially seeded at a 1 to 4 ratio. Similar results were obtained when nTom+, and H2B-nGFP+ myoblasts were differentiated in the same culture. Thus, the high number of nTom+ myonuclei, which we observed after voluntary running, is a vast overestimation of SC fusion. Importantly, these data warrant against the use of nuclear reporter mice in multinucleated cells such as myofibers and suggest a careful re-interpretation of previous data published with Pax7 nuclear reporter mice [57,58,59].

To circumvent nuclear propagation artifacts in our Pax7nTom/+ mice, we assessed recombination at the DNA level from both sorted myonuclei as well as bulk muscle derived from two Pax7 reporter models. Intriguingly, the increase in the fraction of recombined myonuclei from sedentary control to VRun (+ 8%) and (VResRun) (+ 12%) was almost identical to the increases in myonuclear accretion measured on single fiber (Hoechst+) and on cryosections (PCM1+). Myonuclei have a half-life of at least 15 years [60] and are remarkably resistant to apoptosis, even when the myofiber is challenged by unloading or atrophy [61]. This suggests that myonuclear turnover during exercise training is low and that the ~ 16% increase in myonuclei upon exercise as observed in this study is the direct result of SC fusion to myofibers, without active nuclear apoptosis. Thus, based on our data, a major question that arises is whether these new myonuclei are different in terms of transcriptional or epigenetic signature from the old residing ones. In addition, whether they contribute to training adaptations and whether there is a difference between endurance and resistance-based training is yet to be explored. Future research using single nuclei transcriptomics may elucidate potential reprogramming of newly derived myonuclei during exercise.


In summary, this study is the first to assess SC fusion and concomitant myonuclear accretion in a physiological setting of increased mechanical load. Using a Pax7 membrane-marker reporter, we show increased load-dependent SC fusion in the plantar flexors of the mouse hind limb after 8 weeks of voluntary running. Predominantly, the glycolytic type MyHCIIa fibers in SOL underwent hypertrophy, but SC contribution occurred to the same extent in both in MyHCI and MyHCIIa fibers. Thus, the lack of hypertrophy in the heavily SC fused MyHCI fibers strongly suggests that myonuclear accretion might also play a role in endurance adaptations to interval training. Furthermore, in an attempt to objectively quantify SC-derived myonuclear accretion after voluntary running, we used a Pax7 nuclear reporter mouse and unexpectedly observed massive increase in nTom positive nuclei. Further in vitro co-culture experiments showed that this was due to propagation of nTom between nuclei sharing the same cytoplasm. Hence, future investigations should warrant caution in using nuclear reporter constructs in multinucleated cells. Finally, assessment of Cre-mediated recombination in Pax7CreERT2 reporter models indicated that SC contribution to myofibers mirrors myonuclear accretion during exercise.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.



Satellite cell


Tibialis anterior


Extensor digitorum longus








Paired Box Protein 7


Membrane-bound Tomato


Nuclear tdTomato


Membrane-bound Green Fluorescent Protein


Pericentriolar material 1


Myosin heavy chain


Mammalian target of rapamycin


Cross-sectional area


  1. Caiozzo VJ. Plasticity of skeletal muscle phenotype: mechanical consequences. Muscle Nerve. 2002;26:740–68.

    PubMed  Google Scholar 

  2. Morton RW, Murphy KT, McKellar SR, Schoenfeld BJ, Henselmans M, Helms E, et al. A systematic review, meta-analysis and meta-regression of the effect of protein supplementation on resistance training-induced gains in muscle mass and strength in healthy adults. Br J Sports Med. 2017:bjsports-2017-097608 Available from:

  3. Chargé SBP, Rudnicki MA. Cellular and molecular regulation of muscle regeneration. Physiol Rev. 2004;84:209–38.

    PubMed  Google Scholar 

  4. Tierney MT, Sacco A. Satellite cell heterogeneity in skeletal muscle homeostasis. Trends Cell Biol [Internet]. Elsevier Ltd. 2016;26:434–44 Available from:

    CAS  Google Scholar 

  5. Dungan CM, Murach KA, Frick KK, Jones SR, Crow SE, Englund DA, et al. Elevated myonuclear density during skeletal muscle hypertrophy in response to training is reversed during detraining. Am J Physiol Physiol. 2019;316:C649–54.

    CAS  Google Scholar 

  6. Goh Q, Millay DP. Requirement of myomaker-mediated stem cell fusion for skeletal muscle hypertrophy. Elife. 2017;6:1–19.

    Google Scholar 

  7. Fry CS, Lee JD, Jackson JR, Kirby TJ, Stasko SA, Liu H, et al. Regulation of the muscle fiber microenvironment by activated satellite cells during hypertrophy. FASEB J. 2014;28:1654–65.

    CAS  PubMed  PubMed Central  Google Scholar 

  8. Goh Q, Song T, Petrany MJ, Cramer AA, Sun C, Sadayappan S, et al. Myonuclear accretion is a determinant of exercise-induced remodeling in skeletal muscle. Elife. 2019;8:1–19.

    Google Scholar 

  9. Keefe AC, Lawson JA, Flygare SD, Fox ZD, Colasanto MP, Mathew SJ, et al. Muscle stem cells contribute to myofibres in sedentary adult mice. Nat Commun Nature Publishing Group. 2015;6:1–11.

    Google Scholar 

  10. Pawlikowski B, Pulliam C, Betta ND, Kardon G, Olwin BB. Pervasive satellite cell contribution to uninjured adult muscle fibers. Skelet Muscle. 2015;5:1–13 Available from:

    Google Scholar 

  11. Abreu P, Mendes SVD, Ceccatto VM, Hirabara SM. Satellite cell activation induced by aerobic muscle adaptation in response to endurance exercise in humans and rodents. Life Sci. 2017;170:33–40 Available from:

    CAS  PubMed  Google Scholar 

  12. Cisterna B, Giagnacovo M, Costanzo M, Fattoretti P, Zancanaro C, Pellicciari C, et al. Adapted physical exercise enhances activation and differentiation potential of satellite cells in the skeletal muscle of old mice. J Anat. 2016;228:771–83.

    PubMed  PubMed Central  Google Scholar 

  13. Shefer G, Rauner G, Yablonka-Reuveni Z, Benayahu D. Reduced satellite cell numbers and myogenic capacity in aging can be alleviated by endurance exercise. PLoS One. 2010;5.

  14. Kurosaka M, Naito H, Ogura Y, Kojima A, Goto K, Katamoto S. Effects of voluntary wheel running on satellite cells in the rat plantaris muscle. J Sports Sci Med. 2009;8:51–7.

    PubMed  PubMed Central  Google Scholar 

  15. Smith HK, Merry TL. Voluntary resistance wheel exercise during post-natal growth in rats enhances skeletal muscle satellite cell and myonuclear content at adulthood. Acta Physiol. 2012;204:393–402.

    CAS  Google Scholar 

  16. Jackson JR, Kirby TJ, Fry CS, Cooper RL, Mccarthy JJ, Peterson CA, et al. Reduced voluntary running performance is associated with impaired coordination as a result of muscle satellite cell depletion in adult mice. Skelet Muscle [Internet]. Skelet Muscle; 2015;1–17. Available from:

  17. Li P, Akimoto T, Zhang M, Williams RS, Yan Z. Resident stem cells are not required for exercise-induced fiber-type switching and angiogenesis but are necessary for activity-dependent muscle growth. Am J Phys Cell Phys. 2006;290:1461–8.

    Google Scholar 

  18. De Bono JP, Adlam D, Paterson DJ, Channon KM. Novel quantitative phenotypes of exercise training in mouse models. Am J Physiol Integr Comp Physiol. 2005;290:R926–34.

    Google Scholar 

  19. O’Connor RS, Pavlath GK, McCarthy JJ, Esser KA. Last word on point:counterpoint: satellite cell addition is/is not obligatory for skeletal muscle hypertrophy. J Appl Physiol. 2007;103:1107.

    PubMed  Google Scholar 

  20. Murach KA, Fry CS, Kirby TJ, Jackson JR, Lee JD, White SH, et al. Starring or supporting role? Satellite cells and skeletal muscle fiber size regulation. Physiology (Bethesda). 2018;33:26–38 Available from:

    CAS  Google Scholar 

  21. Cornelison DDW. “Known Unknowns”: current questions in muscle satellite cell biology. 1st Curr. Top. Dev. Biol. Elsevier Inc.; 2018. Available from:

  22. Soffe Z, Radley-Crabb HG, McMahon C, Grounds MD, Shavlakadze T. Effects of loaded voluntary wheel exercise on performance and muscle hypertrophy in young and old male C57Bl/6J mice. Scand J Med Sci Sports. 2016;26:172–88.

    CAS  Google Scholar 

  23. Legerlotz K, Elliott B, Guillemin B, Smith HK. Voluntary resistance running wheel activity pattern and skeletal muscle growth in rats. Exp Physiol. 2008;93:754–62.

    PubMed  Google Scholar 

  24. White Z, Terrill J, White RB, McMahon C, Sheard P, Grounds MD, et al. Voluntary resistance wheel exercise from mid-life prevents sarcopenia and increases markers of mitochondrial function and autophagy in muscles of old male and female C57BL/6J mice. Skelet Muscle. 2016;6:1–21 Available from:

    CAS  Google Scholar 

  25. Mccarthy JJ, Dupont-versteegden EE, Fry CS, Murach KA, Peterson CA. Methodological issues limit interpretation of negative effects of satellite cell depletion on adult muscle hypertrophy. Development. 2017;144:1363–7.

    CAS  PubMed  Google Scholar 

  26. Manzanares G, Brito-Da-Silva G, Gandra PG. Voluntary wheel running: patterns and physiological effects in mice. Braz J Med Biol Res. 2019;52:1–9.

    Google Scholar 

  27. Bergmeister KD, Gröger M, Aman M, Willensdorfer A, Manzano-Szalai K, Salminger S, et al. Automated muscle fiber type population analysis with ImageJ of whole rat muscles using rapid myosin heavy chain immunohistochemistry. Muscle Nerve. 2016;54:292–9.

    CAS  PubMed  Google Scholar 

  28. Winje IM, Bengtsen M, Eftestøl E, Juvkam I, Bruusgaard JC, Gundersen K. Specific labelling of myonuclei by an antibody against pericentriolar material 1 (PCM1) on skeletal muscle tissue sections. Acta Physiol. 2018;12:e13034.

    Google Scholar 

  29. Tierney MT, Stec MJ, Sacco A. Assessing muscle stem cell clonal complexity during aging. In: Turksen K, editor. Methods Mol Biol vol 2045. New York: Humana; 2018.

    Google Scholar 

  30. Noirez P, Torres S, Cebrian J, Agbulut O, Peltzer J, Butler-Browne G, et al. TGF-β1 favors the development of fast type identity during soleus muscle regeneration. J Muscle Res Cell Motil. 2006:1–8.

  31. D’Hulst G, Soro-Arnaiz I, Masschelein E, Veys K, Fitzgerald G, Smeuninx B, et al. PHD1 controls muscle mTORC1 in a hydroxylation-independent manner by stabilizing leucyl tRNA synthetase. Nat Commun Springer US. 2020;11:1–15.

    Google Scholar 

  32. Habib N, Avraham-Davidi I, Basu A, Burks T, Shekhar K, Hofree M, et al. Massively parallel single-nucleus RNA-seq with DroNc-seq. Nat Methods. 2017;14:955–8.

    CAS  PubMed  PubMed Central  Google Scholar 

  33. McCarthy JJ, Srikuea R, Kirby TJ, Peterson CA, Esser KA. Inducible Cre transgenic mouse strain for skeletal muscle-specific gene targeting. Skelet Muscle. 2012:1–7.

  34. Konhilas JP, Widegren U, Allen DL, Paul AC, Cleary A, Leinwand LA. Loaded wheel running and muscle adaptation in the mouse. Am J Physiol Circ Physiol. 2005;289:H455–65.

    CAS  Google Scholar 

  35. Egner IM, Bruusgaard JC. Gundersen K. Development: Satellite cell depletion prevents fiber hypertrophy in skeletal muscle; 2016.

    Google Scholar 

  36. Murach KA, White SH, Wen Y, Ho A, Dupont-Versteegden EE, McCarthy JJ, et al. Differential requirement for satellite cells during overload-induced muscle hypertrophy in growing versus mature mice. Skelet Muscle. 2017;7:1–13.

    Google Scholar 

  37. Murphy MM, Keefe AC, Lawson JA, Flygare SD, Yandell M, Kardon G. Transiently active wnt/β-catenin signaling is not required but must be silenced for stem cell function during muscle regeneration. Stem Cell Reports. 2014;3:475–88 Available from:

    CAS  PubMed  PubMed Central  Google Scholar 

  38. Muzumdar MD, Tasic B, Miyamichi K, Li L, Luo L. A global double-fluorescent Cre reporter mouse. Genesis. 2007;45:593–605.

    CAS  PubMed  Google Scholar 

  39. Pallafacchina G, Blaauw B, Schiaffino S. Role of satellite cells in muscle growth and maintenance of muscle mass. Nutr Metab Cardiovasc Dis. 2013;23:12–8 Available from:

    Google Scholar 

  40. White RB, Biérinx A, Gnocchi VF, Zammit PS. Dynamics of muscle fibre growth during postnatal mouse development. BMC Dev Biol. 2010;10:1–11.

    Google Scholar 

  41. Bruusgaard JC, Liestø K, Gundersen K. Distribution of myonuclei and microtubules in live muscle fibers of young, middle-aged, and old mice. J Appl Physiol. 2006;100:2024–30.

    CAS  PubMed  Google Scholar 

  42. D’Hulst G, Palmer SP, Masschelein E, Bar-Nur O, De Bock K. Voluntary resistance running as a model to induce mTOR activation in mouse skeletal muscle. Front Physiol. 2019;10:1–9.

    Google Scholar 

  43. Yan Z, Okutsu M, Akhtar YN, Lira VA. Regulation of exercise-induced fiber type transformation, mitochondrial biogenesis, and angiogenesis in skeletal muscle. J Appl Physiol. 2011;110:264–74.

    CAS  PubMed  Google Scholar 

  44. Taylor-weiner H, Grigsby CL, Ferreira DMS, Dias JM, Stevens MM. Modeling the transport of nuclear proteins along single skeletal muscle cells. Proc Natl Acad Sci. 2020:1–9.

  45. Carlson-Kuhta P, Trank TV, Smith JL. Forms of forward quadrupedal locomotion. III. A comparison of posture, hindlimb kinematics, and motor patterns for downslope and level walking. J Neurophysiol. 1998;79:1702–16.

    PubMed  Google Scholar 

  46. Call JA, Mckeehen JN, Novotny SA, Lowe DA. Progressive resistance voluntary wheel running in the mdx mouse. Muscle Nerve. 2010;42:871–80.

    CAS  PubMed  PubMed Central  Google Scholar 

  47. Dreyer HC, Blanco CE, Sattler FR, Schroeder ET, Wiswell RA. Satellite cell numbers in young and older men 24 hours after eccentric exercise. Muscle Nerve. 2006;33:242–53.

    PubMed  Google Scholar 

  48. Crameri RM, Langberg H, Magnusson P, Jensen CH, Schrøoder HD, Olesen JL, et al. Changes in satellite cells in human skeletal muscle after a single bout of high intensity exercise. J Physiol. 2004;558:333–40.

    CAS  PubMed  PubMed Central  Google Scholar 

  49. Parise G, McKinnell IW, Rudnicki MA. Muscle satellite cell and atypical myogenic progenitor response following exercise. Muscle Nerve. 2008;37:611–9.

    CAS  PubMed  Google Scholar 

  50. Fukuda S, Kaneshige A, Kaji T, Noguchi YT, Takemoto Y, Zhang L, et al. Sustained expression of HeyL is critical for the proliferation of muscle stem cells in overloaded muscle. Elife. 2019;8:1–21.

    Google Scholar 

  51. Egner IM, Bruusgaard JC, Gundersen K. An apparent lack of effect of satellite cell depletion on hypertrophy could be due to methodological limitations. Response to ‘Methodological issues limit interpretation of negative effects of satellite cell depletion on adult muscle hypertrophy.’. Development. 2017;144:1365–7.

    CAS  PubMed  Google Scholar 

  52. Blaauw B, Canato M, Agatea L, Toniolo L, Mammucari C, Masiero E, et al. Inducible activation of Akt increases skeletal muscle mass and force without satellite cell activation. FASEB J. 2009;23:3896–905.

    CAS  PubMed  Google Scholar 

  53. Raffaello A, Milan G, Masiero E, Carnio S, Lee D, Lanfranchi G, et al. JunB transcription factor maintains skeletal muscle mass and promotes hypertrophy. J Cell Biol. 2010;191:101–13.

    CAS  PubMed  PubMed Central  Google Scholar 

  54. Hawley JA, Hargreaves M, Joyner MJ, Zierath JR. Integrative biology of exercise. Cell. 2014;159:738–49 Available from:

    CAS  PubMed  Google Scholar 

  55. Conceiçao M, Vechin FC, Lixandrão ME, Damas F, Libardi CA, Tricoli V, et al. Muscle fiber hypertrophy and myonuclei addition. Med Sci Sports Exerc. 2018;50:1385–93 Available from:

    PubMed  Google Scholar 

  56. Fontaine B, Sassoon D, Buckingham M, Changeux JP. Detection of the nicotinic acetylcholine receptor alpha-subunit mRNA by in situ hybridization at neuromuscular junctions of 15-day-old chick striated muscles. EMBO J. 1988;7:603–9.

    CAS  PubMed  PubMed Central  Google Scholar 

  57. Bachman JF, Klose A, Liu W, Paris ND, Blanc RS, Schmalz M, et al. Prepubertal skeletal muscle growth requires Pax7-expressing satellite cell-derived myonuclear contribution. Development. 2018;145:dev167197.

    PubMed  PubMed Central  Google Scholar 

  58. Tierney MT, Stec MJ, Rulands S, Simons BD, Sacco A. Muscle stem cells exhibit distinct clonal dynamics in response to tissue repair and homeostatic aging. Cell Stem Cell. 2017;22:119–27 Available from:

    PubMed  PubMed Central  Google Scholar 

  59. Liu W, Klose A, Forman S, Paris ND, Wei-LaPierre L, Cortés-Lopéz M, et al. Loss of adult skeletal muscle stem cells drives age-related neuromuscular junction degeneration. Elife. 2017;6:1–15.

    Google Scholar 

  60. Spalding KL, Bhardwaj RD, Buchholz BA, Druid H, Frisén J. Retrospective birth dating of cells in humans. Cell. 2005;122:133–43.

    CAS  PubMed  Google Scholar 

  61. Bruusgaard JC, Gundersen K. In vivo time-lapse microscopy reveals no loss of murine myonuclei during weeks of muscle atrophy. J Clin Invest. 2008;118:1450–7.

    CAS  PubMed  PubMed Central  Google Scholar 

Download references


We gratefully acknowledge Nicola Bundschuh and Paola Gilardoni for their skilled technical assistance and Andrew Palmer for his thoughtful discussions about the experiments.


ETH Zurich internal funding

Author information

Authors and Affiliations



Evi Masschelein, conceptualization, data curation, formal analysis, investigation, methodology, writing original draft, and review and editing; Gommaar D’Hulst, conceptualization, data curation, formal analysis, investigation, methodology, writing original draft, and review and editing; Laura Hinte, formal analysis, methodology, and manuscript editing; Joel Zvick, formal analysis, methodology, and manuscript editing; Inés Soro-Arnaiz, formal analysis, methodology, and manuscript editing; Tatiane Gorski, formal analysis, methodology, and manuscript editing; Ferdinand von Meyenn, conceptualization, investigation, and manuscript editing; Ori Bar-Nur, conceptualization, investigation, and manuscript editing; and Katrien De Bock, conceptualization, formal analysis, funding acquisition, investigation, methodology, writing original draft, project administration, and editing. The author(s) read and approved the final manuscript.

Corresponding author

Correspondence to Katrien De Bock.

Ethics declarations

Ethics approval

All animal procedures were approved by the Veterinary office of the Canton of Zürich (license nr. ZH254-16).

Consent for publication

Not applicable

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

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

Supplementary information

Additional file 1: Figure S1.

Running pattern. (A) Average running speed, (B) running time and (C) number of bouts per night throughout the exercise protocol (n = 16 mice per group). Statistics: two-way ANOVA test with a Bonferroni post hoc test. Group x time interaction is indicated in the figure. (*p <0.05; **p < 0.01; ***p < 0.001). Line graphs represent mean ± SEM (error bars).

Additional file 2: Figure S2.

Muscle hypertrophy is muscle dependent, fiber type dependent and load-dependent. (A) Representative cross-sections for m. soleus (SOL) stained for MyHCI (green), MyHCIIa (blue) and Laminin (red). (B,C) Quantification of fiber cross-sectional area distribution and (D,E) average fiber cross-sectional area. (F) Representative cross-sections for plantaris (PLT) stained for MyHCIIa (blue), MyHCIIb (yellow) and laminin. (G,H) quantification of fiber cross-sectional area distribution and (I,J) average fiber cross-sectional area. (K) Fiber cross-sectional area distribution for gastrocnemius (GAS), (L) tibialis anterior (TA) and (M) extensor digitorum longus (EDL). Statistics: one-way ANOVA test with Tukey correction for multiple comparisons (*p < 0.05; **p < 0.01; ***p < 0.001) (a, p < 0.05 VResRun compared to Control; b, p < 0.05 VResRun compared to VRun; c, p < 0.05 VRun compared to Control). Bar graphs and line graphs represent mean ± SEM (error bars). Scale bar, 100 μm.

Additional file 3: Figure S3.

Fiber remodeling and central nuclei. (A) Myofibers stained for embryonic myosin heavy chain (eMHC), laminin and hoechst. No eMHC+ fibers were detected in Control, VRun, VResRun. Glycerol injected muscle was used as a positive control. (B) Quantification of fibers containing one or more centrally located nucleus. Arrows indicate central nuclei. Bar graph represents mean ± SEM (error bars). Statistics: one-way ANOVA test with Tukey correction for multiple comparisons (*p < 0.05). Each dot represents a single mouse. Scale bar, 50 μm.

Additional file 4: Figure S4.

SC fusion is not fiber type dependent. (A) Representative sequential cross-sections for soleus for mGFP and mTomato (left panel) and stained for MyHCI (yellow) and MyHCIIa (blue) (right panel). (B) Quantification of GFP+ fibers for MyHCI and (C) MyHCIIa. (D) Fiber type distribution and (E) fiber type distribution of GFP+ fibers. Statistics: one-way ANOVA test with Tukey correction for multiple comparisons (*p < 0.05; **p < 0.01; ***p < 0.001). Each dot represents a single mouse. Bar graphs and line graphs represent mean ± SEM (error bars). White arrows indicate GFP+ MyHCI+ fibers, purple arrows indicate GFP+ MyHCIIa+ fibers. Scale bar, 100 μm.

Additional file 5: Figure S5.

Myonuclear accretion is load-dependent. (A) Representative cross-sections of soleus, plantaris and gastrocnemius stained for PCM1 (yellow) and Hoechst (blue). (B) Quantification of PCM1+/Hoechst+ nuclei per fiber. Statistics: one-way ANOVA test with Tukey correction for multiple comparisons (*p < 0.05; **p < 0.01; ***p < 0.001). Bar graphs and line graphs represent mean ± SEM (error bars). White arrows indicate PCM1+ nuclei. Scale bar, 50 μm.

Additional file 7: Table S1.

Phenotypic characterization and muscle weights normalized to tibia length. Statistics: one-way ANOVA test with Tukey correction for multiple comparisons (*p < 0.05 VResRun vs. Control). Values represent mean ± SEM. n = 9-12 mice per group.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Masschelein, E., D’Hulst, G., Zvick, J. et al. Exercise promotes satellite cell contribution to myofibers in a load-dependent manner. Skeletal Muscle 10, 21 (2020).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: