Exercise is associated with younger methylome and transcriptome profiles in human skeletal muscle.

Exercise training prevents age-related decline in muscle function. Targeting epigenetic aging is a promising actionable mechanism and late-life exercise mitigates epigenetic aging in rodent muscle. Whether exercise training can decelerate, or reverse epigenetic aging in humans is unknown. Here, we performed a powerful meta-analysis of the methylome and transcriptome of an unprecedented number of human skeletal muscle samples (n = 3176). We show that: (1) individuals with higher baseline aerobic fitness have younger epigenetic and transcriptomic profiles, (2) exercise training leads to significant shifts of epigenetic and transcriptomic patterns toward a younger profile, and (3) muscle disuse "ages" the transcriptome. Higher fitness levels were associated with attenuated differential methylation and transcription during aging. Furthermore, both epigenetic and transcriptomic profiles shifted toward a younger state after exercise training interventions, while the transcriptome shifted toward an older state after forced muscle disuse. We demonstrate that exercise training targets many of the age-related transcripts and DNA methylation loci to maintain younger methylome and transcriptome profiles, specifically in genes related to muscle structure, metabolism, and mitochondrial function. Our comprehensive analysis will inform future studies aiming to identify the best combination of therapeutics and exercise regimes to optimize longevity.

Early cardiac aging linked to impaired stress-resistance and transcriptional control of stress response, quality control and mitochondrial pathways.

Phenotypic and transcriptomic evidence of early cardiac aging, and associated mechanisms, were investigated in young to middle-aged male mice (C57Bl/6; ages 8, 16, 32, 48 wks). Left ventricular gene expression (profiled via Illumina MouseWG-6 BeadChips), contractile and coronary function, and stress-resistance were assessed in Langendorff perfused hearts under normoxic conditions and following ischemic insult (20 min global ischemia-45 min reperfusion; I-R). Baseline or normoxic contractile function was unaltered by age, while cardiac and coronary 'reserves' (during ß-adrenoceptor stimulation; 1 µM isoproterenol) declined by 48 wks. Resistance to I-R injury fell from 16 to 32 wks. Age-dependent transcriptional changes In un-stressed hearts were limited to 104 genes (>1.3-fold; 0.05 FDR), supporting: up-regulated innate defenses (glutathione and xenobiotic metabolism, chemotaxis, interleukins) and catecholamine secretion; and down-regulated extracellular matrix (ECM), growth factor and survival (PI3K/Akt) signaling. In stressed (post-ischemic) myocardium, ∼15-times as many genes (1528) were age-dependent, grouped into 6 clusters (>1.3-fold change; 0.05 FDR): most changing from 16 wks (45 % up/44 % down), a further 5 % declining from 32 wks. Major age-dependent Biological Processes in I-R hearts reveal: declining ATP metabolism, oxidative phosphorylation, cardiac contraction and morphogenesis, phospholipid metabolism and calcineurin signaling; increasing proteolysis and negative control of MAPK; and mixed changes in nuclear transport and angiogenic genes. Pathway analysis supports reductions in: autophagy, stress response, ER protein processing, mRNA surveillance and ribosome/translation genes; with later falls in mitochondrial biogenesis, oxidative phosphorylation and proteasome genes in I-R hearts. Summarizing, early cardiac aging is evident from 16 to 32 wks in male mice, characterized by: declining cardiovascular reserve and stress-resistance, transcriptomic evidence of constitutive stress and altered catecholamine and survival/growth signaling in healthy hearts; and declining stress response, quality control, mitochondrial energy metabolism and cardiac modeling processes in stressed hearts. These very early changes, potentially key substrate for advanced aging, may inform approaches to healthy aging and cardioprotection in the adult heart.

The muscle proteome reflects changes in mitochondrial function, cellular stress and proteolysis after 14 days of unilateral lower limb immobilization in active young men.

Skeletal muscle unloading due to joint immobilization induces muscle atrophy, which has primarily been attributed to reductions in protein synthesis in humans. However, no study has evaluated the skeletal muscle proteome response to limb immobilization using SWATH proteomic methods. This study characterized the shifts in individual muscle protein abundance and corresponding gene sets after 3 and 14 d of unilateral lower limb immobilization in otherwise healthy young men. Eighteen male participants (25.4 ±5.5 y, 81.2 ±11.6 kg) underwent 14 d of unilateral knee-brace immobilization with dietary provision and following four-weeks of training to standardise acute training history. Participant phenotype was characterized before and after 14 days of immobilization, and muscle biopsies were obtained from the vastus lateralis at baseline (pre-immobilization) and at 3 and 14 d of immobilization for analysis by SWATH-MS and subsequent gene-set enrichment analysis (GSEA). Immobilization reduced vastus group cross sectional area (-9.6 ±4.6%, P <0.0001), immobilized leg lean mass (-3.3 ±3.9%, P = 0.002), unilateral 3-repetition maximum leg press (-15.6 ±9.2%, P <0.0001), and maximal oxygen uptake (-2.9 ±5.2%, P = 0.044). SWATH analyses consistently identified 2281 proteins. Compared to baseline, two and 99 proteins were differentially expressed (FDR <0.05) after 3 and 14 d of immobilization, respectively. After 14 d of immobilization, 322 biological processes were different to baseline (FDR <0.05, P <0.001). Most (77%) biological processes were positively enriched and characterized by cellular stress, targeted proteolysis, and protein-DNA complex modifications. In contrast, mitochondrial organization and energy metabolism were negatively enriched processes. This study is the first to use data independent proteomics and GSEA to show that unilateral lower limb immobilization evokes mitochondrial dysfunction, cellular stress, and proteolysis. Through GSEA and network mapping, we identify 27 hub proteins as potential protein/gene candidates for further exploration.

Effect of short-term hindlimb immobilization on skeletal muscle atrophy and the transcriptome in a low compared with high responder to endurance training model.

Skeletal muscle atrophy is a physiological response to disuse, aging, and disease. We compared changes in muscle mass and the transcriptome profile after short-term immobilization in a divergent model of high and low responders to endurance training to identify biological processes associated with the early atrophy response. Female rats selectively bred for high response to endurance training (HRT) and low response to endurance training (LRT; n = 6/group; generation 19) underwent 3 day hindlimb cast immobilization to compare atrophy of plantaris and soleus muscles with line-matched controls (n = 6/group). RNA sequencing was utilized to identify Gene Ontology Biological Processes with differential gene set enrichment. Aerobic training performed prior to the intervention showed HRT improved running distance (+60.6 ± 29.6%), while LRT were unchanged (-0.3 ± 13.3%). Soleus atrophy was greater in LRT vs. HRT (-9.0 ±8.8 vs. 6.2 ±8.2%; P<0.05) and there was a similar trend in plantaris (-16.4 ±5.6% vs. -8.5 ±7.4%; P = 0.064). A total of 140 and 118 biological processes were differentially enriched in plantaris and soleus muscles, respectively. Soleus muscle exhibited divergent LRT and HRT responses in processes including autophagy and immune response. In plantaris, processes associated with protein ubiquitination, as well as the atrogenes (Trim63 and Fbxo32), were more positively enriched in LRT. Overall, LRT demonstrate exacerbated atrophy compared to HRT, associated with differential gene enrichments of biological processes. This indicates that genetic factors that result in divergent adaptations to endurance exercise, may also regulate biological processes associated with short-term muscle unloading.

Meta-analysis of genome-wide DNA methylation and integrative omics of age in human skeletal muscle.

BACKGROUND: Knowledge of age-related DNA methylation changes in skeletal muscle is limited, yet this tissue is severely affected by ageing in humans. METHODS: We conducted a large-scale epigenome-wide association study meta-analysis of age in human skeletal muscle from 10 studies (total n = 908 muscle methylomes from men and women aged 18-89 years old). We explored the genomic context of age-related DNA methylation changes in chromatin states, CpG islands, and transcription factor binding sites and performed gene set enrichment analysis. We then integrated the DNA methylation data with known transcriptomic and proteomic age-related changes in skeletal muscle. Finally, we updated our recently developed muscle epigenetic clock (https://bioconductor.org/packages/release/bioc/html/MEAT.html). RESULTS: We identified 6710 differentially methylated regions at a stringent false discovery rate <0.005, spanning 6367 unique genes, many of which related to skeletal muscle structure and development. We found a strong increase in DNA methylation at Polycomb target genes and bivalent chromatin domains and a concomitant decrease in DNA methylation at enhancers. Most differentially methylated genes were not altered at the mRNA or protein level, but they were nonetheless strongly enriched for genes showing age-related differential mRNA and protein expression. After adding a substantial number of samples from five datasets (+371), the updated version of the muscle clock (MEAT 2.0, total n = 1053 samples) performed similarly to the original version of the muscle clock (median of 4.4 vs. 4.6 years in age prediction error), suggesting that the original version of the muscle clock was very accurate. CONCLUSIONS: We provide here the most comprehensive picture of DNA methylation ageing in human skeletal muscle and reveal widespread alterations of genes involved in skeletal muscle structure, development, and differentiation. We have made our results available as an open-access, user-friendly, web-based tool called MetaMeth (https://sarah-voisin.shinyapps.io/MetaMeth/).

Low responders to endurance training exhibit impaired hypertrophy and divergent biological process responses in rat skeletal muscle.

NEW FINDINGS: What is the central question of this study? The extent to which genetics determines adaptation to endurance versus resistance exercise is unclear. Previously, a divergent selective breeding rat model showed that genetic factors play a major role in the response to aerobic training. Here, we asked: do genetic factors that underpin poor adaptation to endurance training affect adaptation to functional overload? What is the main finding and its importance? Our data show that heritable factors in low responders to endurance training generated differential gene expression that was associated with impaired skeletal muscle hypertrophy. A maladaptive genotype to endurance exercise appears to dysregulate biological processes responsible for mediating exercise adaptation, irrespective of the mode of contraction stimulus. ABSTRACT: Divergent skeletal muscle phenotypes result from chronic resistance-type versus endurance-type contraction, reflecting the principle of training specificity. Our aim was to determine whether there is a common set of genetic factors that influence skeletal muscle adaptation to divergent contractile stimuli. Female rats were obtained from a genetically heterogeneous rat population and were selectively bred from high responders to endurance training (HRT) or low responders to endurance training (LRT; n = 6/group; generation 19). Both groups underwent 14 days of synergist ablation to induce functional overload of the plantaris muscle before comparison to non-overloaded controls of the same phenotype. RNA sequencing was performed to identify Gene Ontology biological processes with differential (LRT vs. HRT) gene set enrichment. We found that running distance, determined in advance of synergist ablation, increased in response to aerobic training in HRT but not LRT (65 ± 26 vs. -6 ± 18%, mean ± SD, P < 0.0001). The hypertrophy response to functional overload was attenuated in LRT versus HRT (20.1 ± 5.6 vs. 41.6 ± 16.1%, P = 0.015). Between-group differences were observed in the magnitude of response of 96 upregulated and 101 downregulated pathways. A further 27 pathways showed contrasting upregulation or downregulation in LRT versus HRT in response to functional overload. In conclusion, low responders to aerobic endurance training were also low responders for compensatory hypertrophy, and attenuated hypertrophy was associated with differential gene set regulation. Our findings suggest that genetic factors that underpin aerobic training maladaptation might also dysregulate the transcriptional regulation of biological processes that contribute to adaptation to mechanical overload.

An epigenetic clock for human skeletal muscle.

BACKGROUND: Ageing is associated with DNA methylation changes in all human tissues, and epigenetic markers can estimate chronological age based on DNA methylation patterns across tissues. However, the construction of the original pan-tissue epigenetic clock did not include skeletal muscle samples and hence exhibited a strong deviation between DNA methylation and chronological age in this tissue. METHODS: To address this, we developed a more accurate, muscle-specific epigenetic clock based on the genome-wide DNA methylation data of 682 skeletal muscle samples from 12 independent datasets (18-89 years old, 22% women, 99% Caucasian), all generated with Illumina HumanMethylation (HM) arrays (HM27, HM450, or HMEPIC). We also took advantage of the large number of samples to conduct an epigenome-wide association study of age-associated DNA methylation patterns in skeletal muscle. RESULTS: The newly developed clock uses 200 cytosine-phosphate-guanine dinucleotides to estimate chronological age in skeletal muscle, 16 of which are in common with the 353 cytosine-phosphate-guanine dinucleotides of the pan-tissue clock. The muscle clock outperformed the pan-tissue clock, with a median error of only 4.6 years across datasets (vs. 13.1 years for the pan-tissue clock, P < 0.0001) and an average correlation of ρ = 0.62 between actual and predicted age across datasets (vs. ρ = 0.51 for the pan-tissue clock). Lastly, we identified 180 differentially methylated regions with age in skeletal muscle at a false discovery rate < 0.005. However, gene set enrichment analysis did not reveal any enrichment for gene ontologies. CONCLUSIONS: We have developed a muscle-specific epigenetic clock that predicts age with better accuracy than the pan-tissue clock. We implemented the muscle clock in an r package called Muscle Epigenetic Age Test available on Bioconductor to estimate epigenetic age in skeletal muscle samples. This clock may prove valuable in assessing the impact of environmental factors, such as exercise and diet, on muscle-specific biological ageing processes.