Research Paper Volume 11, Issue 18 pp 7587—7604
Attenuated activation of the unfolded protein response following exercise in skeletal muscle of older adults
- 1 Division of Endocrinology and Metabolism, Mayo Clinic College of Medicine, Rochester, MN 55905, USA
- 2 Division of Biostatistics and Informatics, Mayo Clinic College of Medicine, Rochester, MN 55905, USA
- 3 Department of Physical Activity and Health, Eastern Oregon University, La Grande, OR 97850, USA
Received: June 29, 2019 Accepted: September 5, 2019 Published: September 14, 2019https://doi.org/10.18632/aging.102273
How to Cite
Copyright © 2019 Hart et al. This is an open-access article distributed under the terms of the Creative Commons Attribution (CC BY 3.0) License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Sarcopenia is linked with impaired adaptive responses to exercise in aging skeletal muscle. The unfolded protein response (UPR) is an important intramyocellular molecular response pathway that is activated by exercise. The influence of age on skeletal muscle adaptive UPR in response to exercise, and the relationship to other key exercise-responsive regulatory pathways is not well-understood. We evaluated age-related changes in transcriptional markers of UPR activation following a single bout of resistance exercise in 12 young (27 ± 5yrs) and 12 older (75 ± 5yrs) healthy men and women. At baseline, there were modest differences in expression of UPR-related genes in young and older adults. Following exercise, transcriptional markers of UPR pathway activation were attenuated in older adults compared to young based on specific salient UPR-related genes and gene set enrichment analysis. The coordination of post-exercise transcriptional patterns between the UPR pathway, p53/p21 axis of autophagy, and satellite cell differentiation were less evident in older compared to young adults. In conclusion, transcriptomic analysis revealed an age-related decline in the adaptive UPR transcriptional response following a single bout of exercise that could contribute to impaired exercise responsiveness with age.
Aging is associated with the loss of skeletal muscle mass, quality, and function; decrements that have a negative influence on health span . When skeletal muscle function deteriorates to the point that activities of daily living can no longer be performed independently, quality of life is compromised, and risk of disability, morbidity, and mortality increases substantially . Resistance exercise improves muscle mass and function, but there is emerging evidence that the molecular and cellular responses to anabolic stimuli (e.g., exercise and nutrition) are attenuated in older adults; a phenomenon termed anabolic resistance . In the quest to understand the molecular drivers of age-related anabolic resistance, the physiological unfolded protein response (UPR) has emerged as a key regulatory pathway in skeletal muscle protein quality control and adaptations to exercise [3–6]. Early evidence points to altered UPR as an explanation for age and disease related changes in protein folding and accumulation and aggregation of proteins within the endoplasmic reticulum (ER) [7–10].
The UPR is a highly conserved cellular stress response within the ER. Accumulation of misfolded proteins in the ER lumen trigger ER stress and subsequent initiation of the UPR, which is mediated by three ER transmembrane sensors: RNA-dependent protein kinase-like ER eukaryotic translation initiation factor 2 alpha kinase (PERK), inositol-requiring protein 1 (IRE1), and activating transcription factor-6 (ATF6) . Detection of misfolded or unfolded proteins by the ER triggers the dissociation of glucose-regulated protein 78 (Grp78) from PERK, IRE1, and ATF6, initiating a signaling cascade to restore proteostasis within the cell . PERK phosphorylates eukaryotic translation initiation factor 2α (eIF2α) to attenuate mRNA translation [11, 13]. IRE1 phosphorylation promotes the activation and splicing of X-box-binding protein 1 (XBP1s) mRNA [11, 14]. ATF6 translocates to the Golgi, where it is processed by proteases and releases the active cytosolic fragment (ATF6f), which controls genes responsible for encoding ER associated degradation (ERAD) and XBP1 . Taken together, the UPR signaling cascade modulates chaperone and protein expression under physiological conditions to maintain proteostasis and prevent future ER stress .
The role of the UPR has been studied and characterized in several physiological and pathologic conditions. Components of the UPR play a role in feeding , exercise adaptations , and chronic diseases including diabetes, obesity, Alzheimer’s, and cancer  in tissues such as in the liver , pancreas , and adipose tissue . Under resting conditions, tonic activation of the UPR in skeletal muscle has been implicated in mechanisms that both prevent  and contribute [22, 23] to age-associated sarcopenia and anabolic resistance. Selective activation of UPR branches, sometimes referred to as the adaptive UPR, in response to both aerobic  and resistance exercise  demonstrate an important role of the UPR in mediating skeletal muscle adaptations. However, less studied is the extent to which age influences markers of the adaptive UPR in skeletal muscle in response to exercise, and the relationship to other key exercise-responsive regulatory pathways (e.g., autophagy, hypertrophy, mitochondrial biogenesis). We recently observed that healthy older adults exhibited decreased muscle strength accompanied by reduced muscle mitochondrial function and attenuated transcriptional responses to a single bout of exercise compared to young adults . Here we leveraged the availability of residual biospecimens to evaluate age-related alterations in markers associated with UPR pathway activation in response to a single bout of resistance exercise in young and older adults.
Post-exercise UPR pathway activation attenuated with age
The characteristics of the twelve young and older men and women have been previously described  and are provided in in Table 1. Most notably, there were no significant differences between the young and older adults in BMI, total body fat percentage, total skeletal muscle index (SMI), or fat and lean mass in the legs. SMI is a morphological calculation used to define sarcopenia [25, 26]. While there was no significant difference in absolute VO2 peak (L/min) between young and older adults, VO2 peak normalized to body mass, absolute leg strength (1RM), and leg strength relative to fat free mass in the leg were significantly (p≤0.05) decreased in older compared to young adults.
Table 1. Subject Characteristics.
|Young N=12||Older N=12||P Age|
|Age (years)||27 ± 5||76 ± 5|
|Height (cm)||169.6 ± 3.7||163.7 ± 8.5||0.03*|
|Weight (kg)||69.9 ± 7.0||70.7 ± 12||0.08|
|BMI (kg/m2)||24.3 ± 2.7||26.3 ± 2.8||0.10|
|SBP (mmHg)||112.4 ± 11.8||124.2 ± 11.6||0.02*|
|DBP (mmHg)||69.9 ± 9.7||64.3 ± 10.5||0.18|
|Total body Fat (%)||31.7 ± 7.7||36.6 ± 5.9||0.10|
|Fat, arms (kg)||2.0 ± 0.5||2.5 ± 0.7||0.04*|
|Fat, legs (kg)||7.6 ± 2.2||7.8 ± 2.2||0.82|
|Fat, trunk (kg)||10.6 ± 2.7||13.7 ± 5.1||0.08|
|Total lean mass (kg)||45.5 ± 7.8||42.5 ± 7.3||0.30|
|Lean, arms (kg)||5.0 ± 1.2||4.5 ± 1.1||0.31|
|Lean, legs (kg)||15.6 ± 3.1||13.9 ± 2.5||0.15|
|Lean, trunk (kg)||22 ± 3.6||21.1 ± 3.6||0.51|
|SMI (kg/m2)||15.8 ± 2.6||15.8 ± 1.4||0.67|
|1RM (AU)||84.4 ± 21.4||52.1 ± 14.6||0.001*|
|1RM (AU/kg FFM)||1.74 ± 0.2||1.15 ± 0.19||<0.0001*|
|VO2 peak (L/min)||2.2 ± 0.5||1.5 ± 0.4||0.10|
|VO2 peak (ml/kg/min)||31.4 ± 5.6||21.7 ± 4.2||<0.001*|
|M; Male, F; Female, BMI; body mass index, SBP; systolic blood pressure, DBP; diastolic blood pressure, SMI; total skeletal muscle index, 1RM; 1 repetition maximum, AU; arbitrary units. VO2 peak; peak oxygen uptake. Data are shown as mean ± SD, * p ≤ 0.05, Wilcoxon rank sum test for age group comparisons and for non-parametric distributions.|
The whole-muscle transcriptome dataset generated from RNA sequencing was used to identify age-related differences in the UPR pathway using gene set enrichment analysis (GSEA). The GSEA calculates an enrichment score (ES), which represents how frequently genes of the UPR pathway occur at the top or bottom of the overall ranked dataset , thereby providing a readout of the coordination of multiple genes within the UPR pathway. At baseline before exercise, transcripts in the UPR pathway were not remarkably different between young and older adults (Figure 1A) and the overall enrichment of the UPR pathway visualized by GSEA was not significantly enriched in older adults, although there was a notable trend (ES = 0.39, FDR = 0.0842) toward increased tonic UPR activation in older adults (Figure 1B). Age-related differences in UPR activation following exercise were assessed from the fold change (pre- to post-exercise) in the UPR-associated gene transcripts (Figure 2). The volcano plots in Figure 2 demonstrate that UPR-related genes are influenced by an acute exercise bout to a greater extent in young (Figure 2A) compared to older adults (Figure 2B). Subsequent GSEA further quantified these age-associated differences and revealed a significant increase in UPR pathway activation in young (ES = 0.57, FDR = 0.0405, Figure 2C) that was not statistically significant in older adults (ES = 0.46, FDR = 0.0619, Figure 2D). Examination of individual salient UPR-related genes reveals attenuated changes post-exercise in older adults compared to young (Figure 2E). The overall coordinated response to exercise among all detected UPR-related genes was stronger in young (Figure 3) compared to older adults (Figure 3), as evidenced by stronger correlations among UPR pathway gene transcripts from baseline to 18 hours post-exercise. There were no age-related differences in the expression of proteins involved in the UPR from baseline to 18 hours post-exercise (Supplementary Figure 2). Expression of phosphorylating proteins that are indicative of ER stress and UPR pathway activation typically peak at 48 hours post-exercise [4, 5], which could explain why no differences were detected in the present study.
Figure 1. Pathway analysis of the unfolded protein response mRNA transcripts at baseline. Baseline UPR genes transcripts are not different at baseline in young vs. older adults (A) The x-axis represents the log fold change, while y-axis represents the –log FDR-value for each gene. UPR genes are colored in blue. Dotted vertical lines represent the threshold for statistical significance (FDR < 0.05, −log [FDR value] > 1.303). Gene set enrichment analysis (GSEA) of 72 gene transcripts associated with the UPR reveals a trend towards upregulation of mRNA for the unfolded protein response in skeletal muscle under basal conditions (B) in the older compared to younger adults (FDR < 0.05 and an absolute log2 fold change R ≥ 0.5).
Figure 2. Pathway analysis of the unfolded protein response mRNA transcripts post-exercise. The unfolded protein response (UPR) in skeletal muscle volcano plots of mRNA expressed as a log fold change from baseline to 18h post-exercise in young (A) and older (B) adults. The x-axis represents the log fold change, while y-axis represents the –log FDR-value for each gene. Dotted vertical lines represent the threshold for statistical significance (FDR < 0.05, −log [FDR value] > 1.303). Gene set enrichment analysis (GSEA) was performed using Broad’s GSEA software for 72 genes related to UPR pathway and revealed a significantly greater upregulation of the UPR pathway between rest and 18h post-exercise in young (C) compared to older (D) adults (FDR < 0.05 and an absolute log2 fold change ≥ 0.5). (E) UPR up- and down-regulated mRNAs significantly altered (FDR < 0.05 and an absolute fold change of ≥ 1.2) between baseline and 18h post-exercise in young and older adults, sorted by the log2 fold change. NS = not significantly altered.
Figure 3. Pathway activation analysis demonstrated a stronger coordinated response of gene transcripts associated with the UPR post-exercise in young (bottom left) compared to older (upper right) adults. The correlation matrix heat maps demonstrate significant positive (dark red) and negative (dark blue) relationships between the fold-change from baseline to 18h post-exercise for the individual gene transcripts associated with the UPR pathway in young and older adults.
Although the small sample size of this study precludes an adequately powered assessment of sex differences or age-by-sex interactions in the post-exercise UPR transcriptional response using GSEA, a preliminary comparison of individual UPR gene transcripts in males and females within each age group suggests that sex differences may exist. For example, in the young adults the pre-to-post log change in PERK was significantly (p≤0.05) greater in females (1.24 ± 0.19) than males (0.64 ± 0.41), while in older adults the pre-to-post log change was significantly (p≤0.05) greater in the females than males for ATF3 (females: 2.54 ± 1.76, males: 0.19 ± 1.36) and Grp78 (females: 0.40 ± 0.21, males: 0.07 ± 0.25). When age groups were combined to examine sex effects irrespective of age, the pre-to-post log changes in key UPR gene transcripts were significantly (p≤0.05) greater in females compared to males (ATF3: females: 3.71 ± 1.74, males: 1.36 ± 1.94; PERK: females: 0.84 ± 0.50, males: 0.27 ± 0.40; Grp78: females: 0.61 ± 0.38, males: 0.24 ± 0.24; and XBP1: females: 0.40 ± 0.36, males: 0.12 ± .18).
The PERK arm of the UPR pathway and the relationship to autophagy and satellite cell gene transcripts
When we interrogated UPR-related transcripts specific to each arm of the UPR, we observed that genes associated with the PERK arm were most enriched in young adults post-exercise (Figure 2E). Specifically, PERK and ATF3 mRNA exhibited a significantly greater fold change from baseline to 18 hours post-exercise in young compared to older adults. Emerging evidence supports crosstalk between the UPR and autophagy to restore cellular homeostasis . In the present study, the gene transcripts associated with the p53/p21 axis of autophagy were significantly elevated in older compared to young adults at baseline (Figure 4A), but the fold change 18 hours post-exercise in these genes was attenuated in older adults (Figure 4B). As observed with the overall UPR pathway activation, young adults exhibited a more coordinated response between genes associated with the PERK arm of the UPR and the p53/p21 axis, indicated by the associations between fold change in salient gene transcripts from baseline to 18 hours post-exercise in young compared to older adults (Figure 4C and Supplementary Figure 3). Similar to UPR-related transcripts, preliminary sex differences were evident for p53/p21 gene transcripts. Young females exhibited significantly (p≤0.05) greater pre-to-post exercise log fold change in p21 (females: 5.18 ± 0.62, males: 2.88 ± 2.47), p53 (females: 1.54 ± 0.41, males: 0.74 ± 0.51), CASP3 (females: 1.56 ± 0.32, males: 0.64 ± 0.91), and GADD45A (females: 2.47 ± 1.29, males: 0.65 ± 1.31) compared to young males, and BECN1 was significantly (p≤0.05) downregulated in young females (-0.33 ± 0.10) compared to young males (-0.10 ± 0.14). When age groups were combined to examine sex effects irrespective of age, the pre-to-post log changes were significantly (p≤0.05) greater in females compared to males for p21 (females: 3.71 ± 1.74, males: 1.36 ± 1.94), p53 (females: 1.12 ± 0.52, males: 0.54 ± 0.64), SESN1 (females: 0.96 ± 0.45, males: 0.50 ± 0.47), CASP3 (females: 0.84 ± 0.50, males: 0.27 ± 0.40), GADD45A (females: 0.85 ± 0.45, males: 0.45 ± 0.45), BCL2 (females: 0.51 ± 0.34, males: 0.19 ± 0.23), BECN1 (females: 0.69 ± 0.30, males: 0.36 ± 0.29), and CYCS (females: 0.42 ± 0.41, males: 0.10 ± 0.23). Despite some intriguing evidence supporting potential sex differences in post-exercise UPR and autophagy-related transcripts, these trends should be interpreted cautiously due to the fact that this study was designed to be inclusive of both sexes, but was not powered to support a formal interrogation of sex differences or age-by-sex interactions.
Figure 4. Selected gene transcripts associated with p53/p21 axis of autophagy and the relationship to the PERK arm of the UPR. Gene transcripts associated with p53/p21 axis are elevated at baseline (A) and decreased 18h post-exercise (B) in older compared to younger adults. Correlations between gene transcripts associated with the PERK arm of the UPR and p53/p21 axis post-exercise (C) in young (bottom left) compared to older (upper right) adults. The correlation matrix heat maps demonstrate significant positive (dark red) and negative (dark blue) relationships between the fold-change from baseline to 18h post-exercise for the individual gene transcripts associated with the UPR pathway in young and older adults. *Significantly (p ≤ 0.05) different between groups. Data bars are presented as mean ± SEM.
Both p21 and PERK have recently been linked to alterations in satellite cell homeostasis and regenerative myogenesis in skeletal muscle . When the metabolic demand on skeletal muscle is low, satellite cells reside in a quiescent state and are activated during periods of disrupted homeostatic imbalance to produce myonuclei, induce hypertrophy, and repair myofibers for muscle regeneration . Skeletal muscle satellite cells express the paired box (PAX) transcription factors Pax3 and Pax7. Satellite cells express high levels of Pax7 in the quiescent state, while Pax7 expression declines in satellite cells upon differentiation during muscle regeneration and repair . Conversely, Pax3 has been documented to be expressed in satellite cells during both the quiescent and activated state . In the present study, transcripts corresponding to satellite cell markers, Pax3, Pax7, and MyoD were similar in young and older adults at baseline and 18 hours post-exercise, but there was a clear attenuation of MyoG expression post-exercise in older adults (Figure 5A). Furthermore, gene transcripts associated with the PERK arm of the UPR and p53/p21 axis exhibited stronger correlations with Pax3, Pax7, MyoD, and MyoG in young compared to older adults when expressed as a fold change from baseline to 18 hours post-exercise (Figure 5B).
Figure 5. The relationships between markers of satellite cell differentiation, the p53/p21 axis of autophagy and the PERK arm of the UPR. The fold change from baseline to 18h post-exercise in gene transcripts representing satellite cell differentiation in young and older adults (A). Correlations demonstrate significant positive (dark red) and negative (dark blue) relationships between the fold-change from baseline to 18h post-exercise for the individual gene transcripts associated with the satellite cells differentiation markers (Pax3, Pax7, MyoD, MyoG), the p53/p21 axis, and the PERK arm of the UPR pathway in young and older adults (B). Data bars are presented as mean ± SEM. *Significant correlation (p ≤ 0.05). † Significantly different (p ≤ 0.05) between young and older adults.
Stress response pathways are critical to maintaining healthspan, yet many of these pathways are impaired with aging  in conjunction with skeletal muscle impairments commonly observed after the seventh decade of life . The present study was designed to interrogate the transcriptional response of the UPR pathway to a single bout of resistance exercise in skeletal muscle of healthy young and older adults. We found that the post-exercise changes in UPR-associated gene transcripts were attenuated in older adults. Furthermore, the coordination of post-exercise gene expression patterns between the UPR pathway, p53/p21 axis of autophagy, and satellite cell differentiation were less evident in older compared to young adults. These data suggest that the transcriptional UPR pathway response to resistance exercise is attenuated with age, even in the absence of frailty or chronic disease. Further, the lower UPR response to resistance exercise observed with aging coincides with a decline in other key regulatory pathways known to be important in maintenance of skeletal muscle function .
Few studies have assessed ER stress and activation of the UPR following an acute bout of exercise in skeletal muscle [3–5, 33, 34] and far fewer have evaluated aging as a modifier of post-exercise UPR activation in skeletal muscle [4, 5]. Our data support the possibility that the transcriptional patterns of skeletal muscle UPR activation following exercise are attenuated in healthy older adults; a conclusion that contrasts two previous studies where young and older males demonstrated similar post-exercise mRNA fold changes of Grp78, IRE1, PERK, ATF6, ATF4, CHOP, eIF2α, and GADD34, and protein expression of Grp78, IRE1, and PERK [4, 5]. Evidence of age and sex influences on skeletal muscle adaptations to resistance exercise has been previously documented . Despite some initial evidence of sex-related differences in the log fold change of individual UPR and autophagy gene transcripts post-exercise between females and males in the present study, we did not feel confident in our ability to evaluate any sex by age interaction with sufficient statistical power due to the relatively low number of females within each age group. Future studies are needed to evaluate potential sex-related differences in UPR pathway activation following exercise with advancing age.
Skeletal muscle adaptations to exercise are mediated by a myriad of metabolic and molecular pathways  and the coordination among these pathways identifies patterns of coexpression that assist in the identification of potential molecular signals contributing to skeletal muscle function in the context of aging and disease. A progressive decline in the adaptive UPR, characterized by a decline ER molecular chaperones and protein folding capacity, has been implicated in the age-associated reduction of cellular functions in the brain, liver, pancreas, and spleen [37, 38]. Prior evidence suggests an early link between chronic UPR activation and anabolic resistance in skeletal muscle through the inhibition of mTORC1 [22, 23]. However, in the present study, age-related differences in the magnitude and coordination of a large number of UPR-associated gene transcripts post-exercise were present in the absence of any significant differences in the UPR pathway activation at baseline between young and older adults (Figure 1), suggesting that chronic UPR activation did not influence the age-related differences in markers of UPR pathway activation post-exercise. Alternatively, chronic low-grade inflammation and oxidative stress appear to influence the pattern of downstream UPR pathway activation in pancreatic beta cells . It is therefore interesting to note that systemic concentrations of interleukin-6 and tumor necrosis factor alpha were elevated in older adults in this study as reported previously . Based on these observations, the influence of chronic, low-grade systemic inflammation on the adaptive potential of the UPR in skeletal muscle post-exercise merits further investigation.
In the present study, correlations between the fold change of transcripts associated with the PERK arm of the UPR and the p53/p21 axis of autophagy post-exercise were stronger in young (Figure 4C) compared to older adults (Figure 4D). Unaccustomed exercise can evoke damage to cellular components which can be degraded in lysosomes through the process of autophagy . Prior evidence indicates that ER stress induces autophagy in skeletal muscle . The tumor suppressor protein, p53, is responsible for cell repair and enhancing cell survival in response to stress . Acute activation of p53 attenuates the functional decline in skeletal muscle with age through another tumor suppressor protein, p21 . Recent studies show acute induction of autophagy in skeletal muscle 48 hours following resistance exercise, as well as an age-related decline in this adaptive stress response . However, increased p53 and p21 expression under resting conditions has been associated with senescence and an aging phenotype in skeletal muscle . Indeed, in the present study, older adults exhibited elevated p53 and p21 mRNA expression at baseline (Figure 4A), but the post-exercise fold change in gene transcripts associated with the p53/p21 axis of autophagy was blunted in older compared to younger adults (Figure 4B). In support of an age-related difference in post-exercise autophagy activation in skeletal muscle, Hentilä and colleagues  recently reported an increase in lapidated LC3II, an indicator of autophagosome content, 48 hours following a single bout of unaccustomed resistance exercise, which persisted after 21 weeks of resistance training in young, but not in older men. Further studies are needed to evaluate the impact of age on the relationship between autophagy machinery induced by resistance exercise and UPR activation.
Satellite cells (SC) have been documented to serve an important role in the adaptation of skeletal muscle to resistance training [45, 46]. Upon mechanical loading or injury, SC are activated, proliferate, further differentiate, and fuse onto regenerating muscle fibers . Activation of SC is coordinated by the up- or down-regulation of paired box transcription factors 3 and 7 (Pax3 and Pax7) and the myogenic regulatory factors myogenin 5, D, and G (Myf5, MyoD, and MyoG respectively) . The sequential myogenic program of SC markers are characterized in the quiescent state as Pax7+/Myf5-/MyoD-/MyoG-, followed by Pax7+/Myf5+/MyoD-/MyoG- in the proliferative state, and Pax7-/Myf5-/MyoD-/MyoG+ in the differentiated state before returning to the quiescent state upon muscle repair or regeneration . A reduced return to quiescence accompanied by a decrease in the SC pool size with advancing age has been attributed to reduced myonuclear accretion and hypertrophy in skeletal muscle in response to exercise [48, 49]. In the present study, the fold change in post-exercise MyoG was significantly decreased in older compared to young adults, although there were no age-related differences in Pax3, Pax7, or MyoD mRNA (Figure 5A). A coordinated response between gene transcripts associated with SC, autophagy, and the PERK arm of the UPR was more apparent in young compared to older adults (Figure 5B). The significant inverse correlations between both the PERK arm of the UPR and autophagy in relationship to SC gene transcripts suggests a potential role in myogenic transitioning of SC from the quiescent to the fully differentiated state. Indeed, recent evidence supports a pivotal role of the PERK arm in the regulation of SC homeostasis during regenerative myogenesis following injury  and for the preservation of skeletal muscle mass and function , although the underlying mechanism by which the PERK arm regulates hypertrophy and expansion of the SC pool size in response to exercise is unclear at this time.
We recognize some notable limitations to the present study. UPR activation was assessed 18 hours after acute exercise, which represents an intermediate time point compared to prior studies which evaluated UPR activation at 1-3 hours or 24–48 hours post-exercise [4, 5]. We recognize our ability to definitively link the measured blunted adaptive UPR transcriptional response with robust changes phosphorylation of key UPR and p53/p21 autophagy proteins in older adults was hindered by non-optimal timing of muscle biopsies. Indeed, prior evidence suggests expression of both UPR and autophagy-associated key regulatory proteins occurs with more latency compared to gene expression [4, 5] and in the present study we were unable to detect any alterations in protein expression at 18 hours after exercise. This limitation biopsy timing was the consequence of leveraging existing biospecimens from the previous study design for which biopsy timing was optimized for the measurement of protein fractional synthesis rates. Nonetheless, we feel that the transcriptional evidence for an age-related attenuation in the adaptive UPR in the present study, in combination with our prior evidence of no age-related differences in protein fractional synthesis rates , provides insight on the regulation of muscle protein turnover and quality control in response to exercise. It is also important to consider that dietary intake has strong influence on UPR activation in skeletal muscle [23, 33]. Food intake and macronutrient composition in the present study was carefully controlled for three days prior to and throughout the study in a clinical inpatient setting, thereby decreasing the inherent variability that could potentially obscure age group differences. Finally, the average age of the older cohort in the present study was approximately 15 years older than the group studied by Hentilä and colleagues . Prior evidence supports age-related alterations in stress response pathways within skeletal muscle manifesting after the seventh decade of life . It is therefore likely that a ~15-year difference in age may influence the responsiveness of skeletal muscle to unaccustomed exercise in adults over 70 years of age.
In conclusion, older adults exhibited decreased markers of UPR activation and reduced coordination with autophagy and SC-associated gene transcripts following a single bout of unaccustomed resistance exercise. In contrast, young adults demonstrated strong coordination between UPR genes and key regulatory gene transcripts associated with autophagy and SC differentiation in skeletal muscle post-exercise. Taken together, the present findings suggest a potential age-related impairment in the post-exercise transcriptional response supporting activation of the UPR and coordination with other exercise responsive pathways (i.e., autophagy, SC differentiation) in skeletal muscle that is likely to contribute to sarcopenia and age-related attenuation of adaptive responses to exercise.
Materials and Methods
The experimental design has been previously described . In brief, twelve young (18-35 years) and 12 older (65-85 years) men and women were recruited and provided written informed consent as approved by the Mayo Foundation Institutional Review Board. The study conformed to the principles outlined in the Declaration of Helsinki. All participants completed a screening visit followed by a two-day inpatient study to evaluate acute exercise responsiveness. Participants were included in the study only if they did not participate in routine exercise training more than two days per week to limit the influence of a training effect on UPR responsiveness to a single bout of unaccustomed resistance exercise . Those who met the eligibility reported to the Clinical Research and Trials Unit (CRTU) at Mayo Clinic Hospital, St. Marys Campus following an overnight fast. The screening visit consisted of a physical examination, comprehensive blood tests, dual-energy X-ray absorptiometry (DEXA) to determine whole-body fat and fat free mass (FFM), and physical function tests. Physical function tests consisted of a whole-body peak oxygen uptake (VO2 peak) assessment performed on a stationary cycle ergometer and a maximal knee extensor strength assessment, evaluated on two separate occasions, to calculate the one repetition maximum (1RM) using methods previously described . To calculate 1RM, participants performed a warm-up of 10 repetitions on a plate-loaded knee extensor apparatus with a relatively low weight, followed by 2 sets of up to 10 repetitions with 3 min recovery intervals between sets. 1RM was calculated from the average of two sets, such that 1RM = ω·(1+r/30), where ω is the weight lifted on the plate-loaded knee extensor apparatus in arbitrary units (AU) and r is the number of repetitions per set. Leg strength was also normalized to leg FFM (AU/FFM) using results from the DEXA scan.
Prior to the first in-patient study day, participants were admitted to the CRTU at 1700 hrs and consumed nothing but water after 2200 hrs. The following morning, muscle biopsies were obtained from the right vastus lateralis at 0800 hrs. Following the biopsy, participants were given a standardized meal containing 10 kcal/kg of 20% protein, 50% carbohydrate, and 30% fat at 1200 hrs. At 1600 hrs, participants performed seated unilateral leg extension using only the left leg consisting of 8 sets of 10 repetitions at 70% of their 1 RM, determined during an outpatient visit, followed by 3 minutes of rest between each set. A second meal was given at 1800 hrs with caloric content to achieve weight maintenance, after which participants remained fasted until the second set of biopsies on the following day on the exercised leg (i.e. left vastus lateralis) at 1000 hrs, corresponding to 18 hrs post-exercise. Prior evidence suggests minimal differences in baseline gene expression between the right and left legs , such that exercise-induced changes in gene expression in one leg can be confidently compared to baseline gene expression in the opposite leg.
Total RNA was isolated and RNA libraries were prepared as previously described . In brief, total RNA was isolated from the biopsy collected before exercise and 18 hours post-exercise. Sequencing libraries were prepared with TruSeq RNA Sample Prep Kit v2. Libraries were sequenced on a HiSeq 2000 sequencer using TruSeq SBS sequencing kit version 3 and HCS version 188.8.131.52 software. The RNA-Seq data was analyzed using MAP-RSeq v.1.2.1, the Mayo Bioinformatics Core pipeline. MAP-RSeq consists of alignment with TopHat 2.0.6 against the hg19 genome build and gene counts with the HTSeq software using gene annotation files obtained from Illumina (http://cufflinks.cbcb.umd.edu/). A subset of 75 genes associated with the UPR pathway was derived from Reactome (http://www.reactome.org/) (RRID:nif-0000-03390) . Gene transcripts associated with autophagy, the UPR, or mitochondrial biogenesis that were statistically up- or downregulated with different phenotypes (young vs older adults and baseline vs exercise) were subjected to pathway analysis using WEBGESTALT software. GSEA for UPR pathway activation was performed using Broad’s GSEA software. All enrichment p values were FDR corrected (using Benjamini-Hochberg procedure). Genes with a FDR-corrected p value of ≤ 0.05 and an absolute fold change of ≥ 1.20 (where 0.00 signifies no change) were considered for further analysis.
Quantitative RT-PCR for RNA-Seq validation
Quantitative (q)RT-PCR was performed on transcripts of interest that were identified in the RNA-Seq dataset using procedures we have been previously described . Gene-specific primers are shown in Supplementary Table 1. We verified 9 mRNAs associated with the three primary transducers of the UPR (Supplementary Figure 1). Total RNA was isolated using the RNEasy fibrous tissue kit according to the manufacturer’s instructions. RNA quantity and purity were assessed by spectrophotometric analysis (Nanodrop) in which both the ratios of absorbance at 260 nm to that at 230 nm (A260/230) and the absorbance at 260 nm to that at 280 nm (A260/280) were >1.8. cDNA synthesis was performed using SuperScript III First-Strand Synthesis System for RT-PCR cDNA Synthesis Kit (Invitrogen, according to the manufacturer’s protocol. The cDNA-equivalent of 5 ng RNA was used for amplification in 384-well microtiter plates in a QuantStudio 7 cycler (Applied Biosystems, CA, USA) using SYBR green assays. Cycle threshold (CT) values for individual reactions were normalized against β2 microglobulin expression. All cDNA samples were amplified in duplicate. Relative expression was calculated using the ∆CT method. Data are presented as fold change compared with control, obtained using the ∆∆CT method.
Western blot analysis
Frozen muscle tissue was pulverized in liquid nitrogen and homogenized at 4°C in modified RIPA buffer (150 mM NaCl, 50 mM Tris-HCl pH 8.0, 1.0 mM EDTA, 1.0% (v/v) Nonidet P40) and phosphatase and protease inhibitors (Halt Protease and Phosphatase Inhibitor Cocktail, Thermo Scientific, Catalog #78446). Homogenates were incubated on ice for 5 min, followed by centrifugation at 10,000g to remove insolubilized fragments, and protein was quantitated (Pierce BCA Protein Assay Kit, Thermo Scientific, Catalog #23225). Samples were prepared for western blot in lithium dodecyl sulfate sample buffer (NuPAGE LDS Sample Buffer, Invitrogen, Catalog #NP0007) with 20mM DTT. Samples were heated at 90°C for 10 min, and 22.5 μg protein was added to each well of precast gels (NuPAGE Novex Bis-Tris Mini Gels, Invitrogen). Proteins were separated by electrophoresis and blotted on Immuno-Blot PVDF (Bio-Rad, Catalog #1620177). Membranes were then blocked with Odyssey blocking buffer (LI-COR, Catalog #927-50000) for two hours before incubating overnight with primary antibodies against CHOP (Cell Signaling, Catalog #2895), GRP78 (Santa Cruz Biotechnology, Catalog #sc-376768), p53 (Cell Signaling, Catalog #2524), PDI (Cell Signaling, Catalog #3501) or vinculin protein (EMD Millipore, Catalog #CP74). Proteins were detected using either chemiluminescent or fluorescent methods using anti-mouse or rabbit secondary antibodies. Signal intensity was determined using LI-COR 5.2.5 imaging software.
Unpaired t-tests were used to compare subject characteristics between age groups. For variables measured at baseline and post-exercise, young and older groups were compared using two-way ANOVA (age x time). When variables were not normally distributed, the Wilcoxon signed rank test was used for comparisons pre-post-exercise for both age groups. Pre- and post-exercise values were compared within age groups using paired t-tests and across age-groups using ANOVA with Tukey’s procedure to maintain 5% type 1 error rate for post-hoc testing. Exercise effect on mRNA was evaluated from pre- to post-exercise change using two-way ANOVA (age x time). Statistical analysis was performed using JMP 10 Software (SAS Institute, Cart, NC) and PRISM v7.0e (GraphPad Software Inc, La Jolla, CA).
C.R.H. wrote the manuscript and analyzed the data. Z.C.R. contributed to the protein expression measurements and analyzed data. S.D. performed bioinformatics on RNA sequencing data. K.T.P contributed to the writing of the manuscript and analysis of data. M.P. performed gene expression measurements by RT-PCR. A.Z.L contributed to the design of the study and collected the data. I.R.L. designed the study, collected data, analyzed data, and wrote the manuscript. I.R.L. is the guarantor of this work and, as such, had full access to all the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis.
We are grateful for the expert technical assistance from Roberta Soderberg, Daniel Jakaitis, Jill Schimke, Maximillian Pallauf and Melissa Aakre, and the excellent clinical support and services of the CRTU, the Dan Abraham Healthy Living Center, and the Metabolomics Resource Core.
Conflicts of Interest
C.R.H., Z.C.R., S.D., K.T.P., M.P., A.Z.L., and I.R.L. have no disclosures.
This publication was supported by the Clinical and Translational Science Award KL2 TR-000136, and the CTSA Grant Number UL TR000135 from the National Center for Advancing Translational Sciences (NCATS), a component of the National Institutes of Health (NIH). Additional funding was provided by the Strickland Career Development Award, and U24DK100469 from the National Institute of Diabetes and Digestive and Kidney Diseases which originates from the NIH Director's Common Fund. The contents of this work are solely the responsibility of the authors and do not necessarily represent the official views of the National Institutes of Health.
- 1. Liguori I, Russo G, Aran L, Bulli G, Curcio F, Della-Morte D, Gargiulo G, Testa G, Cacciatore F, Bonaduce D, Abete P. Sarcopenia: assessment of disease burden and strategies to improve outcomes. Clin Interv Aging. 2018; 13:913–27. https://doi.org/10.2147/CIA.S149232 [PubMed]
- 2. Kumar V, Selby A, Rankin D, Patel R, Atherton P, Hildebrandt W, Williams J, Smith K, Seynnes O, Hiscock N, Rennie MJ. Age-related differences in the dose-response relationship of muscle protein synthesis to resistance exercise in young and old men. J Physiol. 2009; 587:211–17. https://doi.org/10.1113/jphysiol.2008.164483 [PubMed]
- 3. Wu J, Ruas JL, Estall JL, Rasbach KA, Choi JH, Ye L, Boström P, Tyra HM, Crawford RW, Campbell KP, Rutkowski DT, Kaufman RJ, Spiegelman BM. The unfolded protein response mediates adaptation to exercise in skeletal muscle through a PGC-1α/ATF6α complex. Cell Metab. 2011; 13:160–69. https://doi.org/10.1016/j.cmet.2011.01.003 [PubMed]
- 4. Hentilä J, Ahtiainen JP, Paulsen G, Raastad T, Häkkinen K, Mero AA, Hulmi JJ. Autophagy is induced by resistance exercise in young men, but unfolded protein response is induced regardless of age. Acta Physiol (Oxf). 2018; 224:e13069. https://doi.org/10.1111/apha.13069 [PubMed]
- 5. Ogborn DI, McKay BR, Crane JD, Parise G, Tarnopolsky MA. The unfolded protein response is triggered following a single, unaccustomed resistance-exercise bout. Am J Physiol Regul Integr Comp Physiol. 2014; 307:R664–69. https://doi.org/10.1152/ajpregu.00511.2013 [PubMed]
- 6. Xiong G, Hindi SM, Mann AK, Gallot YS, Bohnert KR, Cavener DR, Whittemore SR, Kumar A. The PERK arm of the unfolded protein response regulates satellite cell-mediated skeletal muscle regeneration. eLife. 2017; 6:6. https://doi.org/10.7554/eLife.22871 [PubMed]
- 7. Gallot YS, Bohnert KR, Straughn AR, Xiong G, Hindi SM, Kumar A. PERK regulates skeletal muscle mass and contractile function in adult mice. FASEB J. 2019; 33:1946–62. https://doi.org/10.1096/fj.201800683RR [PubMed]
- 8. Bohnert KR, Gallot YS, Sato S, Xiong G, Hindi SM, Kumar A. Inhibition of ER stress and unfolding protein response pathways causes skeletal muscle wasting during cancer cachexia. FASEB J. 2016; 30:3053–68. https://doi.org/10.1096/fj.201600250RR [PubMed]
- 9. Hussain SG, Ramaiah KV. Reduced eIF2alpha phosphorylation and increased proapoptotic proteins in aging. Biochem Biophys Res Commun. 2007; 355:365–70. https://doi.org/10.1016/j.bbrc.2007.01.156 [PubMed]
- 10. Bohnert KR, Goli P, Roy A, Sharma AK, Xiong G, Gallot YS, Kumar A. The Toll-Like Receptor/MyD88/XBP1 Signaling Axis Mediates Skeletal Muscle Wasting during Cancer Cachexia. Mol Cell Biol. 2019; 39:e00184–19. https://doi.org/10.1128/MCB.00184-19 [PubMed]
- 11. Hetz C. The unfolded protein response: controlling cell fate decisions under ER stress and beyond. Nat Rev Mol Cell Biol. 2012; 13:89–102. https://doi.org/10.1038/nrm3270 [PubMed]
- 12. Pfaffenbach KT, Lee AS. The critical role of GRP78 in physiologic and pathologic stress. Curr Opin Cell Biol. 2011; 23:150–56. https://doi.org/10.1016/j.ceb.2010.09.007 [PubMed]
- 13. Bohnert KR, McMillan JD, Kumar A. Emerging roles of ER stress and unfolded protein response pathways in skeletal muscle health and disease. J Cell Physiol. 2018; 233:67–78. https://doi.org/10.1002/jcp.25852 [PubMed]
- 14. Acosta-Alvear D, Zhou Y, Blais A, Tsikitis M, Lents NH, Arias C, Lennon CJ, Kluger Y, Dynlacht BD. XBP1 controls diverse cell type- and condition-specific transcriptional regulatory networks. Mol Cell. 2007; 27:53–66. https://doi.org/10.1016/j.molcel.2007.06.011 [PubMed]
- 15. Hetz C, Chevet E, Oakes SA. Proteostasis control by the unfolded protein response. Nat Cell Biol. 2015; 17:829–38. https://doi.org/10.1038/ncb3184 [PubMed]
- 16. Pfaffenbach KT, Nivala AM, Reese L, Ellis F, Wang D, Wei Y, Pagliassotti MJ. Rapamycin inhibits postprandial-mediated X-box-binding protein-1 splicing in rat liver. J Nutr. 2010; 140:879–84. https://doi.org/10.3945/jn.109.119883 [PubMed]
- 17. Wang S, Kaufman RJ. The impact of the unfolded protein response on human disease. J Cell Biol. 2012; 197:857–67. https://doi.org/10.1083/jcb.201110131 [PubMed]
- 18. Gentile CL, Frye M, Pagliassotti MJ. Endoplasmic reticulum stress and the unfolded protein response in nonalcoholic fatty liver disease. Antioxid Redox Signal. 2011; 15:505–521. https://doi.org/10.1089/ars.2010.3790 [PubMed]
- 19. Scheuner D, Song B, McEwen E, Liu C, Laybutt R, Gillespie P, Saunders T, Bonner-Weir S, Kaufman RJ. Translational control is required for the unfolded protein response and in vivo glucose homeostasis. Mol Cell. 2001; 7:1165–76. https://doi.org/10.1016/S1097-2765(01)00265-9 [PubMed]
- 20. Sha H, He Y, Chen H, Wang C, Zenno A, Shi H, Yang X, Zhang X, Qi L. The IRE1α-XBP1 pathway of the unfolded protein response is required for adipogenesis. Cell Metab. 2009; 9:556–64. https://doi.org/10.1016/j.cmet.2009.04.009 [PubMed]
- 21. Tezze C, Romanello V, Desbats MA, Fadini GP, Albiero M, Favaro G, Ciciliot S, Soriano ME, Morbidoni V, Cerqua C, Loefler S, Kern H, Franceschi C, et al. Age-Associated Loss of OPA1 in Muscle Impacts Muscle Mass, Metabolic Homeostasis, Systemic Inflammation, and Epithelial Senescence. Cell Metab. 2017; 25:1374–1389.e6. https://doi.org/10.1016/j.cmet.2017.04.021 [PubMed]
- 22. Deldicque L, Bertrand L, Patton A, Francaux M, Baar K. ER stress induces anabolic resistance in muscle cells through PKB-induced blockade of mTORC1. PLoS One. 2011; 6:e20993. https://doi.org/10.1371/journal.pone.0020993 [PubMed]
- 23. Deldicque L, Cani PD, Philp A, Raymackers JM, Meakin PJ, Ashford ML, Delzenne NM, Francaux M, Baar K. The unfolded protein response is activated in skeletal muscle by high-fat feeding: potential role in the downregulation of protein synthesis. Am J Physiol Endocrinol Metab. 2010; 299:E695–705. https://doi.org/10.1152/ajpendo.00038.2010 [PubMed]
- 24. Lalia AZ, Dasari S, Robinson MM, Abid H, Morse DM, Klaus KA, Lanza IR. Influence of omega-3 fatty acids on skeletal muscle protein metabolism and mitochondrial bioenergetics in older adults. Aging (Albany NY). 2017; 9:1096–129. https://doi.org/10.18632/aging.101210 [PubMed]
- 25. Nishigori T, Tsunoda S, Obama K, Hisamori S, Hashimoto K, Itatani Y, Okada K, Sakai Y. Optimal Cutoff Values of Skeletal Muscle Index to Define Sarcopenia for Prediction of Survival in Patients with Advanced Gastric Cancer. Ann Surg Oncol. 2018; 25:3596–603. https://doi.org/10.1245/s10434-018-6728-7 [PubMed]
- 26. Teranaka H, Fujisawa S, Sakuma T, Teshigawara H, Matsumura A, Ando T, Suzuki T, Negoro A, Ogusa E, Aoki J, Miyashita K, Ishii Y, Takahashi H, et al. Clinical Impact of Sarcopenia and Skeletal Muscle Mass Change during Chemotherapy on Outcomes of Diffuse Large B-Cell Lymphoma. Blood. 2017 (Suppl 1); 130:5221.
- 27. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, Mesirov JP. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA. 2005; 102:15545–50. https://doi.org/10.1073/pnas.0506580102 [PubMed]
- 28. Senft D, Ronai ZA. UPR, autophagy, and mitochondria crosstalk underlies the ER stress response. Trends Biochem Sci. 2015; 40:141–48. https://doi.org/10.1016/j.tibs.2015.01.002 [PubMed]
- 29. Zammit PS. All muscle satellite cells are equal, but are some more equal than others? J Cell Sci. 2008; 121:2975–82. https://doi.org/10.1242/jcs.019661 [PubMed]
- 30. Collins CA, Gnocchi VF, White RB, Boldrin L, Perez-Ruiz A, Relaix F, Morgan JE, Zammit PS. Integrated functions of Pax3 and Pax7 in the regulation of proliferation, cell size and myogenic differentiation. PLoS One. 2009; 4:e4475. https://doi.org/10.1371/journal.pone.0004475 [PubMed]
- 31. Kourtis N, Tavernarakis N. Cellular stress response pathways and ageing: intricate molecular relationships. EMBO J. 2011; 30:2520–31. https://doi.org/10.1038/emboj.2011.162 [PubMed]
- 32. Welle S, Brooks AI, Delehanty JM, Needler N, Thornton CA. Gene expression profile of aging in human muscle. Physiol Genomics. 2003; 14:149–59. https://doi.org/10.1152/physiolgenomics.00049.2003 [PubMed]
- 33. Deldicque L, Cani PD, Delzenne NM, Baar K, Francaux M. Endurance training in mice increases the unfolded protein response induced by a high-fat diet. J Physiol Biochem. 2013; 69:215–25. https://doi.org/10.1007/s13105-012-0204-9 [PubMed]
- 34. Kim HJ, Jamart C, Deldicque L, An GL, Lee YH, Kim CK, Raymackers JM, Francaux M. Endoplasmic reticulum stress markers and ubiquitin–proteasome pathway activity in response to a 200-km run. Med Sci Sports Exerc. 2011; 43:18–25. https://doi.org/10.1249/MSS.0b013e3181e4c5d1 [PubMed]
- 35. Martel GF, Roth SM, Ivey FM, Lemmer JT, Tracy BL, Hurlbut DE, Metter EJ, Hurley BF, Rogers MA. Age and sex affect human muscle fibre adaptations to heavy-resistance strength training. Exp Physiol. 2006; 91:457–64. https://doi.org/10.1113/expphysiol.2005.032771 [PubMed]
- 36. Egan B, Zierath JR. Exercise metabolism and the molecular regulation of skeletal muscle adaptation. Cell Metab. 2013; 17:162–84. https://doi.org/10.1016/j.cmet.2012.12.012 [PubMed]
- 37. Martínez G, Duran-Aniotz C, Cabral-Miranda F, Vivar JP, Hetz C. Endoplasmic reticulum proteostasis impairment in aging. Aging Cell. 2017; 16:615–23. https://doi.org/10.1111/acel.12599 [PubMed]
- 38. Naidoo N. ER and aging-Protein folding and the ER stress response. Ageing Res Rev. 2009; 8:150–59. https://doi.org/10.1016/j.arr.2009.03.001 [PubMed]
- 39. Chan JY, Cooney GJ, Biden TJ, Laybutt DR. Differential regulation of adaptive and apoptotic unfolded protein response signalling by cytokine-induced nitric oxide production in mouse pancreatic beta cells. Diabetologia. 2011; 54:1766–76. https://doi.org/10.1007/s00125-011-2139-z [PubMed]
- 40. Hubbard VM, Valdor R, Macian F, Cuervo AM. Selective autophagy in the maintenance of cellular homeostasis in aging organisms. Biogerontology. 2012; 13:21–35. https://doi.org/10.1007/s10522-011-9331-x [PubMed]
- 41. Yorimitsu T, Nair U, Yang Z, Klionsky DJ. Endoplasmic reticulum stress triggers autophagy. J Biol Chem. 2006; 281:30299–304. https://doi.org/10.1074/jbc.M607007200 [PubMed]
- 42. Beyfuss K, Hood DA. A systematic review of p53 regulation of oxidative stress in skeletal muscle. Redox Rep. 2018; 23:100–117. https://doi.org/10.1080/13510002.2017.1416773 [PubMed]
- 43. Baker DJ, Weaver RL, van Deursen JM. p21 both attenuates and drives senescence and aging in BubR1 progeroid mice. Cell Rep. 2013; 3:1164–74. https://doi.org/10.1016/j.celrep.2013.03.028 [PubMed]
- 44. Edwards MG, Anderson RM, Yuan M, Kendziorski CM, Weindruch R, Prolla TA. Gene expression profiling of aging reveals activation of a p53-mediated transcriptional program. BMC Genomics. 2007; 8:80–80. https://doi.org/10.1186/1471-2164-8-80 [PubMed]
- 45. Nederveen JP, Snijders T, Joanisse S, Wavell CG, Mitchell CJ, Johnston LM, Baker SK, Phillips SM, Parise G. Altered muscle satellite cell activation following 16 wk of resistance training in young men. Am J Physiol Regul Integr Comp Physiol. 2017; 312:R85–92. https://doi.org/10.1152/ajpregu.00221.2016 [PubMed]
- 46. Kadi F, Schjerling P, Andersen LL, Charifi N, Madsen JL, Christensen LR, Andersen JL. The effects of heavy resistance training and detraining on satellite cells in human skeletal muscles. J Physiol. 2004; 558:1005–12. https://doi.org/10.1113/jphysiol.2004.065904 [PubMed]
- 47. Rudnicki MA, Le Grand F, McKinnell I, Kuang S. The molecular regulation of muscle stem cell function. Cold Spring Harb Symp Quant Biol. 2008; 73:323–31. https://doi.org/10.1101/sqb.2008.73.064 [PubMed]
- 48. Snijders T, Nederveen JP, McKay BR, Joanisse S, Verdijk LB, van Loon LJ, Parise G. Satellite cells in human skeletal muscle plasticity. Front Physiol. 2015; 6:283. https://doi.org/10.3389/fphys.2015.00283 [PubMed]
- 49. Snijders T, Verdijk LB, van Loon LJ. The impact of sarcopenia and exercise training on skeletal muscle satellite cells. Ageing Res Rev. 2009; 8:328–38. https://doi.org/10.1016/j.arr.2009.05.003 [PubMed]
- 50. Catoire M, Mensink M, Boekschoten MV, Hangelbroek R, Müller M, Schrauwen P, Kersten S. Pronounced effects of acute endurance exercise on gene expression in resting and exercising human skeletal muscle. PLoS One. 2012; 7:e51066–51066. https://doi.org/10.1371/journal.pone.0051066 [PubMed]
- 51. Croft D, Mundo AF, Haw R, Milacic M, Weiser J, Wu G, Caudy M, Garapati P, Gillespie M, Kamdar MR, Jassal B, Jupe S, Matthews L, et al. The Reactome pathway knowledgebase. Nucleic Acids Res. 2014; 42:D472–77. https://doi.org/10.1093/nar/gkt1102 [PubMed]