Sex differences in muscle protein expression and DNA methylation in response to exercise training.

BACKGROUND: Exercise training elicits changes in muscle physiology, epigenomics, transcriptomics, and proteomics, with males and females exhibiting differing physiological responses to exercise training. However, the molecular mechanisms contributing to the differing adaptations between the sexes are poorly understood. METHODS: We performed a meta-analysis for sex differences in skeletal muscle DNA methylation following an endurance training intervention (Gene SMART cohort and E-MTAB-11282 cohort). We investigated for sex differences in the skeletal muscle proteome following an endurance training intervention (Gene SMART cohort). Lastly, we investigated whether the methylome and proteome are associated with baseline cardiorespiratory fitness (maximal oxygen consumption; VO2max) in a sex-specific manner. RESULTS: Here, we investigated for the first time, DNA methylome and proteome sex differences in response to exercise training in human skeletal muscle (n = 78; 50 males, 28 females). We identified 92 DNA methylation sites (CpGs) associated with exercise training; however, no CpGs changed in a sex-dependent manner. In contrast, we identified 189 proteins that are differentially expressed between the sexes following training, with 82 proteins differentially expressed between the sexes at baseline. Proteins showing the most robust sex-specific response to exercise include SIRT3, MRPL41, and MBP. Irrespective of sex, cardiorespiratory fitness was associated with robust methylome changes (19,257 CpGs) and no proteomic changes. We did not observe sex differences in the association between cardiorespiratory fitness and the DNA methylome. Integrative multi-omic analysis identified sex-specific mitochondrial metabolism pathways associated with exercise responses. Lastly, exercise training and cardiorespiratory fitness shifted the DNA methylomes to be more similar between the sexes. CONCLUSIONS: We identified sex differences in protein expression changes, but not DNA methylation changes, following an endurance exercise training intervention; whereas we identified no sex differences in the DNA methylome or proteome response to lifelong training. Given the delicate interaction between sex and training as well as the limitations of the current study, more studies are required to elucidate whether there is a sex-specific training effect on the DNA methylome. We found that genes involved in mitochondrial metabolism pathways are differentially modulated between the sexes following endurance exercise training. These results shed light on sex differences in molecular adaptations to exercise training in skeletal muscle.

Genetic variants within the COL5A1 gene are associated with ligament injuries in physically active populations from Australia, South Africa, and Japan.

Previous small-scale studies have shown an association between the COL5A1 gene and anterior cruciate ligament (ACL) injury risk. In this larger study, the genotype and allele frequency distributions of the COL5A1 rs12722 C/T and rs10628678 AGGG/deletion (AGGG/-) indel variants were compared between participants: (i) with ACL injury in independent and combined cohorts from South-Africa (SA) and Australia (AUS) vs controls (CON), and (ii) with any ligament (ALL) or only ACL injury in a Japanese (JPN) cohort vs CON. Samples were collected from SA (235 cases; 232 controls), AUS (362 cases; 80 controls) and JPN (500 cases; 1,403 controls). Genomic DNA was extracted and genotyped. Distributions were compared, and inferred haplotype analyses performed. No independent associations were noted for rs12722 or rs10628678 when the combined SA + AUS cohort was analysed. However, the C-deletion (rs12722-rs10628678) inferred haplotype was under-represented (p = 0.040, OR = 0.15, CI = 0.04-0.56), while the T-deletion inferred haplotype was over-represented in the female SA + AUS ACL participants versus controls (p < 0.001, OR = 4.74, CI = 1.66-13.55). Additionally, the rs12722 C/C genotype was under-represented in JPN CON vs ACL (p = 0.039, OR = 0.52, 0.27-1.00), while the rs10628678 -/- genotype was associated with increased risk of any ligament injuries (p = 0.035, OR = 1.31, CI = 1.02-1.68) in the JPN cohort. Collectively, these results highlight that a region within the COL5A1 3'-UTR is associated with ligament injury risk. This must be evaluated in larger cohorts and its functional relevance to the structure and capacity of ligaments and joint biomechanics be explored.Highlights The COL5A1 T-deletion inferred haplotype (rs12722-rs10628678) was associated with an increased risk of ACL rupture in the combined SA and AUS female participants.The COL5A1 C-deletion inferred haplotype (rs12722-rs10628678) was associated with a decreased risk of ACL rupture in the combined SA and AUS female participants.The COL5A1 rs12722 C/C and rs10628678 -/- genotypes were associated with increased risk of ACL rupture and of ligament injuries in JPN, respectively.A region within the COL5A1 3'-UTR is associated with risk of ligament injury, including ACL rupture, and therefore the functional significance of this region on ligament capacity and joint biomechanics requires further exploration.

Uncovering the effects of gender affirming hormone therapy on skeletal muscle and epigenetics: protocol for a prospective matched cohort study in transgender individuals (the GAME study).

INTRODUCTION: Gender affirming hormone therapy (GAHT) is increasingly used by transgender individuals and leads to shifts in sex hormone levels. Skeletal muscle is highly responsive to hormone activity, with limited data on the effects of GAHT on different human tissues. Here, we present the protocol for the GAME study (the effects of Gender Affirming hormone therapy on skeletal Muscle training and Epigenetics), which aims to uncover the effects of GAHT on skeletal muscle 'omic' profiles (methylomics, transcriptomics, proteomics, metabolomics) and markers of skeletal muscle health and fitness. METHODS AND ANALYSIS: This study is a prospective age-matched cohort study in transgender adults commencing GAHT (n=80) and age-matched individuals not commencing GAHT (n=80), conducted at Austin Health and Victoria University in Victoria, Australia. Assessments will take place prior to beginning GAHT and 6 and 12 months into therapies in adults commencing GAHT. Age-matched individuals will be assessed at the same time points. Assessments will be divided over three examination days, involving (1) aerobic fitness tests, (2) muscle strength assessments and (3) collection of blood and muscle samples, as well as body composition measurements. Standardised diets, fitness watches and questionnaires will be used to control for key confounders in analyses. Primary outcomes are changes in aerobic fitness and muscle strength, as well as changes in skeletal muscle DNA methylation and gene expression profiles. Secondary outcomes include changes in skeletal muscle characteristics, proteomics, body composition and blood markers. Linear mixed models will be used to assess changes in outcomes, while accounting for repeated measures within participants and adjusting for known confounders. ETHICS AND DISSEMINATION: The Austin Health Human Research Ethics Committee (HREC) and Victoria University HREC granted approval for this study (HREC/77146/Austin-2021). Findings from this project will be published in open-access, peer-reviewed journals and presented to scientific and public audiences. TRIAL REGISTRATION NUMBER: ACTRN12621001415897; Pre-results.

Skeletal muscle methylome and transcriptome integration reveals profound sex differences related to muscle function and substrate metabolism.

Nearly all human complex traits and diseases exhibit some degree of sex differences, with epigenetics being one of the main contributing factors. Various tissues display sex differences in DNA methylation; however, this has not yet been explored in skeletal muscle, despite skeletal muscle being among the tissues with the most transcriptomic sex differences. For the first time, we investigated the effect of sex on autosomal DNA methylation in human skeletal muscle across three independent cohorts (Gene SMART, FUSION, and GSE38291) using a meta-analysis approach, totalling 369 human muscle samples (222 males and 147 females), and integrated this with known sex-biased transcriptomics. We found 10,240 differentially methylated regions (DMRs) at FDR < 0.005, 94% of which were hypomethylated in males, and gene set enrichment analysis revealed that differentially methylated genes were involved in muscle contraction and substrate metabolism. We then investigated biological factors underlying DNA methylation sex differences and found that circulating hormones were not associated with differential methylation at sex-biased DNA methylation loci; however, these sex-specific loci were enriched for binding sites of hormone-related transcription factors (with top TFs including androgen (AR), estrogen (ESR1), and glucocorticoid (NR3C1) receptors). Fibre type proportions were associated with differential methylation across the genome, as well as across 16% of sex-biased DNA methylation loci (FDR < 0.005). Integration of DNA methylomic results with transcriptomic data from the GTEx database and the FUSION cohort revealed 326 autosomal genes that display sex differences at both the epigenome and transcriptome levels. Importantly, transcriptional sex-biased genes were overrepresented among epigenetic sex-biased genes (p value = 4.6e-13), suggesting differential DNA methylation and gene expression between male and female muscle are functionally linked. Finally, we validated expression of three genes with large effect sizes (FOXO3A, ALDH1A1, and GGT7) in the Gene SMART cohort with qPCR. GGT7, involved in antioxidant metabolism, displays male-biased expression as well as lower methylation in males across the three cohorts. In conclusion, we uncovered 8420 genes that exhibit DNA methylation differences between males and females in human skeletal muscle that may modulate mechanisms controlling muscle metabolism and health.

Individual physiological and mitochondrial responses during 12 weeks of intensified exercise.

AIM: Observed effects of exercise are highly variable between individuals, and subject-by-training interaction (i.e., individual response variability) is often not estimated. Here, we measured mitochondrial (citrate synthetase, cytochrome-c oxidase, succinate dehydrogenase, and mitochondrial copy-number), performance markers (Wpeak , lactate threshold [LT], and VO2peak ), and fiber type proportions/expression (type I, type IIa, and type IIx) in multiple time points during 12-week of high-intensity interval training (HIIT) to investigate effects of exercise at the individual level. METHODS: Sixteen young (age: 33.1 ± 9.0 years), healthy men (VO2peak 35-60 ml/min/kg and BMI: 26.4 ± 4.2) from the Gene SMART study completed 12-week of progressive HIIT. Performance markers and muscle biopsies were collected every 4 weeks. We used mixed-models and bivariate growth models to quantify individual response and to estimate correlations between variables. RESULTS: All performance markers exhibited significant (Wpeak 0.56 ± 0.33 p = 0.003, LT 0.37 ± 0.35 p = 0.007, VO2peak 3.81 ± 6.13 p = 0.02) increases overtime, with subject-by-training interaction being present (95% CI: Wpeak 0.09-0.24, LT 0.06-0.18, VO2peak 0.27-2.32). All other measurements did not exhibit significant changes. Fiber type IIa proportions at baseline was significantly associated with all physiological variables (p < 0.05), and citrate synthetase and cytochrome-c oxidase levels at baseline and overtime (i.e., intercept and slope) presented significant covariance (p < 0.05). Finally, low correlations between performance and mitochondrial markers were observed. CONCLUSION: We identified a significant subject-by-training interaction for the performance markers. While for all other measures within-subject variability was too large and interindividual differences in training efficacy could not be verified. Changes in measurements in response to exercise were not correlated, and such disconnection should be further investigated by future studies.

Mapping Robust Genetic Variants Associated with Exercise Responses.

This review summarised robust and consistent genetic variants associated with aerobic-related and resistance-related phenotypes. In total we highlight 12 SNPs and 7 SNPs that are robustly associated with variance in aerobic-related and resistance-related phenotypes respectively. To date, there is very little literature ascribed to understanding the interplay between genes and environmental factors and the development of physiological traits. We discuss future directions, including large-scale exercise studies to elucidate the functional relevance of the discovered genomic markers. This approach will allow more rigour and reproducible research in the field of exercise genomics.

Mitochondrial respiration variability and simulations in human skeletal muscle: The Gene SMART study.

Mitochondrial respiration using the oxygraph-2k respirometer (Oroboros) is widely used to estimate mitochondrial capacity in human skeletal muscle. Here, we measured mitochondrial respiration variability, in a relatively large sample, and for the first time, using statistical simulations, we provide the sample size required to detect meaningful respiration changes following lifestyle intervention. Muscle biopsies were taken from healthy, young men from the Gene SMART cohort, at multiple time points. We utilized samples for each measurement with two technical repeats using two respirometer chambers (n = 160 pairs of same muscle after removal of low-quality samples). We measured the Technical Error of measurement (TEM ) and the coefficient of variation (CV) for each mitochondrial complex. There was a high correlation between measurements from the two chambers (R > 0.7 P < .001) for all complexes, but the TEM was large (7.9-27 pmol s-1  mg-1 ; complex dependent), and the CV was >15% for all complexes. We performed statistical simulations of a range of effect sizes at 80% power and found that 75 participants (with duplicate measurements) are required to detect a 6% change in mitochondrial respiration after an intervention, while for interventions with 11% effect size, ~24 participants are sufficient. The high variability in respiration suggests that the typical sample sizes in exercise studies may not be sufficient to capture exercise-induced changes.

Aerobic capacity and telomere length in human skeletal muscle and leukocytes across the lifespan.

A reduction in aerobic capacity and the shortening of telomeres are hallmarks of the ageing process. We examined whether a lower aerobic capacity is associated with shorter TL in skeletal muscle and/or leukocytes, across a wide age range of individuals. We also tested whether TL in human skeletal muscle (MTL) correlates with TL in leukocytes (LTL). Eighty-two recreationally active, healthy men from the Gene SMART cohort (31.4±8.2 years; body mass index (BMI)=25.3±3.3kg/m2), and 11 community dwelling older men (74.2±7.5years-old; BMI=28.7±2.8kg/m2) participated in the study. Leukocytes and skeletal muscle samples were collected at rest. Relative telomere length (T/S ratio) was measured by RT-PCR. Associations between TL, aerobic capacity (VO2 peak and peak power) and age were assessed with robust linear models. Older age was associated with shorter LTL (45% variance explained, P<0.001), but not MTL (P= 0.7). Aerobic capacity was not associated with MTL (P=0.5), nor LTL (P=0.3). MTL and LTL were correlated across the lifespan (rs=0.26, P=0.03). In healthy individuals, age explain most of the variability of LTL and this appears to be independent of individual aerobic capacity. Individuals with longer LTL also have a longer MTL, suggesting that there might be a shared molecular mechanism regulating telomere length.