Glucocorticoid Receptor (NR3C1) Variants Associate with the Muscle Strength and Size Response to Resistance Training

Methylome and proteome integration in human skeletal muscle uncover group and individual responses to high-intensity interval training

Exercise is a major beneficial contributor to muscle metabolism, and health benefits acquired by exercise are a result of molecular shifts occurring across multiple molecular layers (i.e., epigenome, transcriptome, and proteome). Identifying robust, across-molecular level targets associated with exercise response, at both group and individual levels, is paramount to develop health guidelines and targeted health interventions. Sixteen, apparently healthy, moderately trained (VO2 max = 51.0 ± 10.6 mL min-1  kg-1 ) males (age range = 18-45 years) from the Gene SMART (Skeletal Muscle Adaptive Responses to Training) study completed a longitudinal study composed of 12-week high-intensity interval training (HIIT) intervention. Vastus lateralis muscle biopsies were collected at baseline and after 4, 8, and 12 weeks of HIIT. DNA methylation (~850 CpG sites) and proteomic (~3000 proteins) analyses were conducted at all time points. Mixed models were applied to estimate group and individual changes, and methylome and proteome integration was conducted using a holistic multilevel approach with the mixOmics package. A total of 461 proteins significantly changed over time (at 4, 8, and 12 weeks), whilst methylome overall shifted with training only one differentially methylated position (DMP) was significant (adj.p-value < .05). K-means analysis revealed cumulative protein changes by clusters of proteins that presented similar changes over time. Individual responses to training were observed in 101 proteins. Seven proteins had large effect-sizes >0.5, among them are two novel exercise-related proteins, LYRM7 and EPN1. Integration analysis showed bidirectional relationships between the methylome and proteome. We showed a significant influence of HIIT on the epigenome and more so on the proteome in human muscle, and uncovered groups of proteins clustering according to similar patterns across the exercise intervention. Individual responses to exercise were observed in the proteome with novel mitochondrial and metabolic proteins consistently changed across individuals. Future work is required to elucidate the role of these proteins in response to exercise.

The Glucocorticoid Receptor Polymorphism Landscape in Patients With Diamond Blackfan Anemia Reveals an Association Between Two Clinically Relevant Single Nucleotide Polymorphisms and Time to Diagnosis

NR3C1, the gene encoding the glucocorticoid receptor, is polymorphic presenting numerous single nucleotide polymorphisms (SNPs) some of which are emerging as leading cause in the variability of manifestation and/or response to glucocorticoids in human diseases. Since 60–80% of patients with Diamond Blackfan anemia (DBA), an inherited pure red cell aplasia induced by mutations in ribosomal protein genes became transfusion independent upon treatment with glucocorticoids, we investigated whether clinically relevant NR3C1 SNPs are associated with disease manifestation in DBA. The eight SNPs rs10482605, rs10482616, rs7701443, rs6189/rs6190, rs860457, rs6198, rs6196, and rs33388/rs33389 were investigated in a cohort of 91 European DBA patients. Results were compared with those observed in healthy volunteers (n=37) or present in public genome databases of Italian and European populations. Although, cases vs. control analyses suggest that the frequency of some of the minor alleles is significantly altered in DBA patients with respect to healthy controls or to the Italian or other European registries, lack of consistency among the associations across different sets suggests that overall the frequency of these SNPs in DBA is not different from that of the general population. Demographic data (47 females and 31 males) and driver mutations (44 S and 29 L genes and eight no-known mutation) are known for 81 patients while glucocorticoid response is known, respectively, for 81 (36 responsive and 45 non-responsive) and age of disease onsets for 79 (55 before and 24 after 4months of age) patients. Neither gender nor leading mutations were associated with the minor alleles or with disease manifestation. In addition, none of the SNPs met the threshold in the response vs. non-responsive groups. However, two SNPs (rs6196 and rs860457) were enriched in patients manifesting the disease before 4months of age. Although the exact biomechanistical consequences of these SNPs are unknown, the fact that their configuration is consistent with that of regulatory regions suggests that they regulate changes in glucocorticoid response during ontogeny. This hypothesis was supported by phosphoproteomic profiling of erythroid cells expanded ex vivo indicating that glucocorticoids activate a ribosomal signature in cells from cord blood but not in those from adult blood, possibly providing a compensatory mechanism to the driving mutations observed in DBA before birth.

An intron variant of the GLI family zinc finger 3 (GLI3) gene differentiates resistance training-induced muscle fiber hypertrophy in younger men

We examined the association between genotype and resistance training-induced changes (12 wk) in dual x-ray energy absorptiometry (DXA)-derived lean soft tissue mass (LSTM) as well as muscle fiber cross-sectional area (fCSA; vastus lateralis; n = 109; age = 22 ± 2 y, BMI = 24.7 ± 3.1 kg/m² ). Over 315 000 genetic polymorphisms were interrogated from muscle using DNA microarrays. First, a targeted investigation was performed where single nucleotide polymorphisms (SNP) identified from a systematic literature review were related to changes in LSTM and fCSA. Next, genome-wide association (GWA) studies were performed to reveal associations between novel SNP targets with pre- to post-training change scores in mean fCSA and LSTM. Our targeted investigation revealed no genotype-by-time interactions for 12 common polymorphisms regarding the change in mean fCSA or change in LSTM. Our first GWA study indicated no SNP were associated with the change in LSTM. However, the second GWA study indicated two SNP exceeded the significance level with the change in mean fCSA (P = 6.9 × 10^-7 for rs4675569, 1.7 × 10^-6 for rs10263647). While the former target is not annotated (chr2:205936846 (GRCh38.p12)), the latter target (chr7:41971865 (GRCh38.p12)) is an intron variant of the GLI Family Zinc Finger 3 (GLI3) gene. Follow-up analyses indicated fCSA increases were greater in the T/C and C/C GLI3 genotypes than the T/T GLI3 genotype (P < .05). Data from the Auburn cohort also revealed participants with the T/C and C/C genotypes exhibited increases in satellite cell number with training (P < .05), whereas T/T participants did not. Additionally, those with the T/C and C/C genotypes achieved myonuclear addition in response to training (P < .05), whereas the T/T participants did not. In summary, this is the first GWA study to examine how polymorphisms associate with the change in hypertrophy measures following resistance training. Future studies are needed to determine if the GLI3 variant differentiates hypertrophic responses to resistance training given the potential link between this gene and satellite cell physiology.

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.

Relationship of Anabolic Hormones With Motor Unit Characteristics in Quadriceps Muscle in Healthy and Frail Aging Men

CONTEXT

Anabolic hormones are important factors in maintaining muscle mass for aging men, but their role in overall motor unit structure and function is unclear.

OBJECTIVE

The objective of this work is to determine associations of anabolic and reproductive hormone levels with motor unit characteristics in quadriceps muscle in older healthy and frail men.

DESIGN

This work is an observational cohort study of community-dwelling men.

PARTICIPANTS

Participants included healthy and frail men younger than 65 years.

INTERVENTION

No intervention was performed.

OUTCOME MEASURE

Quantitative assessments of electromyography-derived motor unit potential size (MUP) and compound muscle action potential size (CMAP) of the vastus lateralis muscle.

RESULTS

We studied 98 men (mean ± SD: age 73 ± 6 years; body mass index [BMI] 25.7 ± 4.0 kg/m2; diabetes 11%) of whom 45% were prefrail and 18% frail. After adjusting for age, BMI, and prevalent diabetes, higher total and free testosterone levels were significantly related to larger CMAP (total testosterone: β [95% CI]: 0.3 [0.08-0.53]; free testosterone: 0.34 [0.13-0.56]). Exploratory analysis showed the relationship between free testosterone and CMAP was stronger in frail rather than robust men. In univariate analyses, estradiol was associated with CMAP size (0.37 [0.16-0.57]); and vitamin D was associated with MUP size (0.22 [0.01-0.43]) but these relationships were no longer significant after adjusting for potential confounders.

CONCLUSION

Our data highlight the associations between androgen levels and the electrophysiological characteristics of older men, particularly in the frail. Clinical trials involving administration of androgens will help to elucidate the potential benefits of intervention on neuromuscular function and/or frailty status.

The Pathway to Cancer Cachexia: MicroRNA-Regulated Networks in Muscle Wasting Based on Integrative Meta-Analysis

Cancer cachexia is a multifactorial syndrome that leads to significant weight loss. Cachexia affects 50%-80% of cancer patients, depending on the tumor type, and is associated with 20%-40% of cancer patient deaths. Besides the efforts to identify molecular mechanisms of skeletal muscle atrophy-a key feature in cancer cachexia-no effective therapy for the syndrome is currently available. MicroRNAs are regulators of gene expression, with therapeutic potential in several muscle wasting disorders. We performed a meta-analysis of previously published gene expression data to reveal new potential microRNA-mRNA networks associated with muscle atrophy in cancer cachexia. We retrieved 52 differentially expressed genes in nine studies of muscle tissue from patients and rodent models of cancer cachexia. Next, we predicted microRNAs targeting these differentially expressed genes. We also include global microRNA expression data surveyed in atrophying skeletal muscles from previous studies as background information. We identified deregulated genes involved in the regulation of apoptosis, muscle hypertrophy, catabolism, and acute phase response. We further predicted new microRNA-mRNA interactions, such as miR-27a/Foxo1, miR-27a/Mef2c, miR-27b/Cxcl12, miR-27b/Mef2c, miR-140/Cxcl12, miR-199a/Cav1, and miR-199a/Junb, which may contribute to muscle wasting in cancer cachexia. Finally, we found drugs targeting MSTN, CXCL12, and CAMK2B, which may be considered for the development of novel therapeutic strategies for cancer cachexia. Our study has broadened the knowledge of microRNA-regulated networks that are likely associated with muscle atrophy in cancer cachexia, pointing to their involvement as potential targets for novel therapeutic strategies.

Genetic predisposition score predicts the increases of knee strength and muscle mass after one-year exercise in healthy elderly

This study aims to identify a genetic predisposition score from a set of candidate gene variants that predicts the response to a one-year exercise intervention. 200 participants (aged 60-83 years) were randomly assigned to a fitness (FIT), whole-body vibration (WBV) and control group. Participants in the exercise (FIT and WBV) groups performed a one-year intervention program. Whole-body skeletal muscle mass (SMM) and isometric knee extension strength (PT_IM60) were measured before and after the intervention. A set of 170 muscle-related single nucleotide polymorphisms (SNPs) were genotyped. Stepwise regression analysis was applied to select significantly contributing SNPs for baseline and relative change parameters. A data-driven genetic predisposition score (GPS) was calculated by adding up predisposing alleles for each of the phenotypes. GPS was calculated based on 4 to 8 SNPs which were significantly related to the corresponding phenotypes. These SNPs belong to genes that are involved in myoblast differentiation, muscle and bone growth, myofiber contraction, cytokines and DNA methylation. GPS was related to baseline PT_IM60 and relative changes of SMM and PT_IM60 in the exercise groups, explaining the variance of the corresponding parameter by 3.2%, 14% and 27%, respectively. Adding one increasing allele in the GPS increased baseline PTIM₆₀ by 4.73 Nm, and exercise-induced relative changes of SMM and PT_IM60 by 1.78% and 3.86% respectively. The identified genetic predisposition scores were positively related to baseline knee extension strength and muscle adaptations to exercise in healthy elderly. These findings provide supportive genetic explanations for high and low responders in exercise-induced muscle adaptations.

Tensor decomposition-based unsupervised feature extraction identifies candidate genes that induce post-traumatic stress disorder-mediated heart diseases

Background

Although post-traumatic stress disorder (PTSD) is primarily a mental disorder, it can cause additional symptoms that do not seem to be directly related to the central nervous system, which PTSD is assumed to directly affect. PTSD-mediated heart diseases are some of such secondary disorders. In spite of the significant correlations between PTSD and heart diseases, spatial separation between the heart and brain (where PTSD is primarily active) prevents researchers from elucidating the mechanisms that bridge the two disorders. Our purpose was to identify genes linking PTSD and heart diseases.

Methods

In this study, gene expression profiles of various murine tissues observed under various types of stress or without stress were analyzed in an integrated manner using tensor decomposition (TD).

Results

Based upon the obtained features, ∼ 400 genes were identified as candidate genes that may mediate heart diseases associated with PTSD. Various gene enrichment analyses supported biological reliability of the identified genes. Ten genes encoding protein-, DNA-, or mRNA-interacting proteins—ILF2, ILF3, ESR1, ESR2, RAD21, HTT, ATF2, NR3C1, TP53, and TP63—were found to be likely to regulate expression of most of these ∼ 400 genes and therefore are candidate primary genes that cause PTSD-mediated heart diseases. Approximately 400 genes in the heart were also found to be strongly affected by various drugs whose known adverse effects are related to heart diseases and/or fear memory conditioning; these data support the reliability of our findings.

Conclusions

TD-based unsupervised feature extraction turned out to be a useful method for gene selection and successfully identified possible genes causing PTSD-mediated heart diseases.

Electronic supplementary material

The online version of this article (doi:10.1186/s12920-017-0302-1) contains supplementary material, which is available to authorized users.