DNA methylation assessment from human slow- and fast-twitch skeletal muscle fibers

DNA methylation of exercise-responsive genes differs between trained and untrained men

BACKGROUND

Physical activity is well known for its multiple health benefits and although the knowledge of the underlying molecular mechanisms is increasing, our understanding of the role of epigenetics in long-term training adaptation remains incomplete. In this intervention study, we included individuals with a history of > 15 years of regular endurance or resistance training compared to age-matched untrained controls performing endurance or resistance exercise. We examined skeletal muscle DNA methylation of genes involved in key adaptation processes, including myogenesis, gene regulation, angiogenesis and metabolism.

RESULTS

A greater number of differentially methylated regions and differentially expressed genes were identified when comparing the endurance group with the control group than in the comparison between the strength group and the control group at baseline. Although the cellular composition of skeletal muscle samples was generally consistent across groups, variations were observed in the distribution of muscle fiber types. Slow-twitch fiber type genes MYH7 and MYL3 exhibited lower promoter methylation and elevated expression in endurance-trained athletes, while the same group showed higher methylation in transcription factors such as FOXO3, CREB5, and PGC-1α. The baseline DNA methylation state of those genes was associated with the transcriptional response to an acute bout of exercise. Acute exercise altered very few of the investigated CpG sites.

CONCLUSIONS

Endurance- compared to resistance-trained athletes and untrained individuals demonstrated a different DNA methylation signature of selected skeletal muscle genes, which may influence transcriptional dynamics following a bout of acute exercise. Skeletal muscle fiber type distribution is associated with methylation of fiber type specific genes. Our results suggest that the baseline DNA methylation landscape in skeletal muscle influences the transcription of regulatory genes in response to an acute exercise bout.

DNA Methylation in the Adaptive Response to Exercise

Emerging evidence published over the past decade has highlighted the role of DNA methylation in skeletal muscle function and health, including as an epigenetic transducer of the adaptive response to exercise. In this review, we aim to synthesize the latest findings in this field to highlight: (1) the shifting understanding of the genomic localization of altered DNA methylation in response to acute and chronic aerobic and resistance exercise in skeletal muscle (e.g., promoter, gene bodies, enhancers, intergenic regions, un-annotated regions, and genome-wide methylation); (2) how these global/regional methylation changes relate to transcriptional activity following exercise; and (3) the factors (e.g., individual demographic or genetic features, dietary, training history, exercise parameters, local epigenetic characteristics, circulating hormones) demonstrated to alter both the pattern of DNA methylation after exercise, and the relationship between DNA methylation and gene expression. Finally, we discuss the changes in non-CpG methylation and 5-hydroxymethylation after exercise, as well as the importance of emerging single-cell analyses to future studies-areas of increasing focus in the field of epigenetics. We anticipate that this review will help generate a framework for clinicians and researchers to begin developing and testing exercise interventions designed to generate targeted changes in DNA methylation as part of a personalized exercise regimen.

Fast and slow muscle fiber transcriptome dynamics with lifelong endurance exercise

We investigated fast and slow muscle fiber transcriptome exercise dynamics among three groups of men: lifelong exercisers (LLE, n = 8, 74 ± 1 yr), old healthy nonexercisers (OH, n = 9, 75 ± 1 yr), and young exercisers (YE, n = 8, 25 ± 1 yr). On average, LLE had exercised ∼4 day·wk-1 for ∼8 h·wk-1 over 53 ± 2 years. Muscle biopsies were obtained pre- and 4 h postresistance exercise (3 × 10 knee extensions at 70% 1-RM). Fast and slow fiber size and function were assessed preexercise with fast and slow RNA-seq profiles examined pre- and postexercise. LLE fast fiber size was similar to OH, which was ∼30% smaller than YE (P < 0.05) with contractile function variables among groups, resulting in lower power in LLE (P < 0.05). LLE slow fibers were ∼30% larger and more powerful compared with YE and OH (P < 0.05). At the transcriptome level, fast fibers were more responsive to resistance exercise compared with slow fibers among all three cohorts (P < 0.05). Exercise induced a comprehensive biological response in fast fibers (P < 0.05) including transcription, signaling, skeletal muscle cell differentiation, and metabolism with vast differences among the groups. Fast fibers from YE exhibited a growth and metabolic signature, with LLE being primarily metabolic, and OH showing a strong stress-related response. In slow fibers, only LLE exhibited a biological response to exercise (P < 0.05), which was related to ketone and lipid metabolism. The divergent exercise transcriptome signatures provide novel insight into the molecular regulation in fast and slow fibers with age and exercise and suggest that the ∼5% weekly exercise time commitment of the lifelong exercisers provided a powerful investment for fast and slow muscle fiber metabolic health at the molecular level.NEW & NOTEWORTHY This study provides the first insights into fast and slow muscle fiber transcriptome dynamics with lifelong endurance exercise. The fast fibers were more responsive to exercise with divergent transcriptome signatures among young exercisers (growth and metabolic), lifelong exercisers (metabolic), and old healthy nonexercisers (stress). Only lifelong exercisers had a biological response in slow fibers (metabolic). These data provide novel insights into fast and slow muscle fiber health at the molecular level with age and exercise.

Exertional heat stroke causes long-term skeletal muscle epigenetic reprogramming, altered gene expression, and impaired satellite cell function in mice

The effect of exertional heat stroke (EHS) exposure on skeletal muscles is incompletely understood. Muscle weakness is an early symptom of EHS but is not considered a major target of multiorgan injury. Previously, in a preclinical mouse model of EHS, we observed the vulnerability of limb muscles to a second EHS exposure, suggesting hidden processes contributing to declines in muscle resilience. Here, we evaluated the possible molecular origins of EHS-induced declines in muscle resilience. Female C57BL/6 mice [total n = 56; 28/condition, i.e., EHS and exercise control (EXC)] underwent forced wheel running at 37.5°C/40% relative humidity until symptom limitation (unconsciousness). EXC mice exercised identically at room temperature (22-23°C). After 1 mo of recovery, the following were assessed: 1) specific force and caffeine-induced contracture in soleus (SOL) and extensor digitorum longus (EDL) muscles; 2) transcriptome and DNA methylome responses in gastrocnemius (GAST); and 3) primary satellite cell function (proliferation and differentiation). There were no differences in specific force in either SOL or EDL from EXC. Only EHS solei exhibited lower caffeine sensitivity. EHS GAST exhibited higher RNA expression of genes encoding structural proteins of slow fibers, heat shock proteins, and myogenesis. A total of ∼2,500 differentially methylated regions of DNA that could potentially affect many cell functions were identified. Primary satellite cells exhibited suppressed proliferation rates but normal differentiation responses. Results demonstrate long-term changes in skeletal muscles 1 mo after EHS that could contribute to declines in muscle resilience. Skeletal muscle may join other, more recognized tissues considered vulnerable to long-term effects of EHS.NEW & NOTEWORTHY Exertional heat stroke (EHS) in mice induces long-term molecular and functional changes in limb muscle that could reflect a loss of "resilience" to further stress. The phenotype was characterized by altered caffeine sensitivity and suppressed satellite cell proliferative potential. This was accompanied by changes in gene expression and DNA methylation consistent with ongoing muscle remodeling and stress adaptation. We propose that EHS may induce a prolonged vulnerability of skeletal muscle to further stress or injury.

Tropomyosin 3 (TPM3) function in skeletal muscle and in myopathy

The tropomyosin genes (TPM1-4) contribute to the functional diversity of skeletal muscle fibers. Since its discovery in 1988, the TPM3 gene has been recognized as an indispensable regulator of muscle contraction in slow muscle fibers. Recent advances suggest that TPM3 isoforms hold more extensive functions during skeletal muscle development and in postnatal muscle. Additionally, mutations in the TPM3 gene have been associated with the features of congenital myopathies. The use of different in vitro and in vivo model systems has leveraged the discovery of several disease mechanisms associated with TPM3-related myopathy. Yet, the precise mechanisms by which TPM3 mutations lead to muscle dysfunction remain unclear. This review consolidates over three decades of research about the role of TPM3 in skeletal muscle. Overall, the progress made has led to a better understanding of the phenotypic spectrum in patients affected by mutations in this gene. The comprehensive body of work generated over these decades has also laid robust groundwork for capturing the multiple functions this protein plays in muscle fibers.

DNA methylation alterations in muscle of critically ill patients

BACKGROUND

Intensive care unit (ICU)-acquired weakness can persist beyond ICU stay and has been associated with long-term functional impairment of ICU survivors. Recently, DNA methylation alterations were found in the blood of ICU patients, partially explaining long-term developmental impairment of critically ill children. As illness-induced aberrant DNA methylation theoretically could also be involved in long-term weakness, we investigated whether the DNA methylation signature in muscle of adult critically ill patients differs from that in muscle of healthy controls.

METHODS

Genome-wide methylation was determined (Infinium® HumanMethylationEPIC BeadChips) in DNA extracted from skeletal muscle biopsies that had been collected on Day 8 ± 1 in ICU from 172 EPaNIC-trial patients [66% male sex, median age 62.7 years, median body mass index (BMI) 25.9 kg/m² ] and 20 matched healthy controls (70% male sex, median age 58.0 years, median BMI 24.4 kg/m² ). Methylation status of individual cytosine-phosphate-guanine (CpG) sites of patients and controls was compared with F-tests, using the Benjamini-Hochberg false discovery rate to correct for multiple comparisons. Differential methylation of DNA regions was assessed with bump hunting, with 1000 permutations assessing uncertainty, expressed as family-wise error rate. Gene expression was investigated for 10 representative affected genes.

RESULTS

In DNA from ICU patients, 565 CpG sites, associated with 400 unique genes, were differentially methylated as compared with controls (average difference 3.2 ± 0.1% ranging up to 16.9%, P < 0.00005). Many of the associated genes appeared highly relevant for muscle structure and function/weakness, including genes involved in myogenesis, muscle regeneration, nerve/muscle membrane excitability, muscle denervation/re-innervation, axon guidance/myelination/degeneration/regeneration, synapse function, ion channelling with especially calcium signalling, metabolism (glucose, protein, and fat), insulin signalling, neuroendocrine hormone regulation, mitochondrial function, autophagy, apoptosis, oxidative stress, Wnt signalling, transcription regulation, muscle fat infiltration during regeneration, and fibrosis. In patients as compared with controls, we also identified two hypomethylated regions, spanning 18 and 3 CpG sites in the promoters of the HIC1 and NADK2 genes, respectively (average differences 5.8 ± 0.01% and 12.1 ± 0.04%, family-wise error rate <0.05). HIC1 and NADK2 play important roles in muscle regeneration and postsynaptic acetylcholine receptors and in mitochondrial processes, respectively. Nine of 10 investigated genes containing DNA methylation alterations were differentially expressed in patients as compared with controls (P ≤ 0.03).

CONCLUSIONS

Critically ill patients present with a different DNA methylation signature in skeletal muscle as compared with healthy controls, which in theory could provide a biological basis for long-term persistence of weakness in ICU survivors.

TRIAL REGISTRATION

ClinicalTrials.gov: NCT00512122, registered on 31 July 2007.

Muscle-Specific Cellular and Molecular Adaptations to Late-Life Voluntary Concurrent Exercise

Murine exercise models can provide information on factors that influence muscle adaptability with aging, but few translatable solutions exist. Progressive weighted wheel running (PoWeR) is a simple, voluntary, low-cost, high-volume endurance/resistance exercise approach for training young mice. In the current investigation, aged mice (22-mo-old) underwent a modified version of PoWeR for 8 wk. Muscle functional, cellular, biochemical, transcriptional, and myonuclear DNA methylation analyses provide an encompassing picture of how muscle from aged mice responds to high-volume combined training. Mice run 6-8 km/d, and relative to sedentary mice, PoWeR increases plantarflexor muscle strength. The oxidative soleus of aged mice responds to PoWeR similarly to young mice in every parameter measured in previous work; this includes muscle mass, glycolytic-to-oxidative fiber type transitioning, fiber size, satellite cell frequency, and myonuclear number. The oxidative/glycolytic plantaris adapts according to fiber type, but with modest overall changes in muscle mass. Capillarity increases markedly with PoWeR in both muscles, which may be permissive for adaptability in advanced age. Comparison to published PoWeR RNA-sequencing data in young mice identified conserved regulators of adaptability across age and muscles; this includes Aldh1l1 which associates with muscle vasculature. Agrn and Samd1 gene expression is upregulated after PoWeR simultaneous with a hypomethylated promoter CpG in myonuclear DNA, which could have implications for innervation and capillarization. A promoter CpG in Rbm10 is hypomethylated by late-life exercise in myonuclei, consistent with findings in muscle tissue. PoWeR and the data herein are a resource for uncovering cellular and molecular regulators of muscle adaptation with aging.

Daily running enhances molecular and physiological circadian rhythms in skeletal muscle

Objective

Exercise is a critical component of a healthy lifestyle and a key strategy for the prevention and management of metabolic disease. Identifying molecular mechanisms underlying adaptation in response to chronic physical activity is of critical interest in metabolic physiology. Circadian rhythms broadly modulate metabolism, including muscle substrate utilization and exercise capacity. Here, we define the molecular and physiological changes induced across the daily cycle by voluntary low intensity daily exercise.

Methods

Wildtype C57BL6/J male and female mice were housed with or without access to a running wheel for six weeks. Maximum running speed was measured at four different zeitgeber times (ZTs, hours after lights on) using either electrical or manual stimulation to motivate continued running on a motorized treadmill. RNA isolated from plantaris muscles at six ZTs was sequenced to establish the impact of daily activity on genome-wide transcription. Patterns of gene expression were analyzed using Gene Set Enrichment Analysis (GSEA) and Detection of Differential Rhythmicity (DODR). Blood glucose, lactate, and ketones, and muscle and liver glycogen were measured before and after exercise.

Results

We demonstrate that the use of mild electrical shocks to motivate running negatively impacts maximum running speed in mice, and describe a manual method to motivate running in rodent exercise studies. Using this method, we show that time of day influences the increase in exercise capacity afforded by six weeks of voluntary wheel running: when maximum running speed is measured at the beginning of the nighttime active period in mice, there is no measurable benefit from a history of daily voluntary running, while maximum increase in performance occurs at the end of the night. We show that daily voluntary exercise dramatically remodels the murine muscle circadian transcriptome. Finally, we describe daily rhythms in carbohydrate metabolism associated with the time-dependent response to moderate daily exercise in mice.

Conclusions

Collectively, these data indicate that chronic nighttime physical activity dramatically remodels daily rhythms of murine muscle gene expression, which in turn support daily fluctuations in exercise performance.

•

Daily voluntary running dramatically remodels the mouse muscle circadian transcriptome.

•

Daily voluntary running maximally increases mouse running speed in the late active period.

•

Muscle and liver glycogen content exhibit robust daily rhythms in laboratory mice.

•

Use of mild electric shocks to motivate running in mice impairs maximum running speed.

The Aging Athlete: Paradigm of Healthy Aging

The Exercise Boom of the 1970's resulted in the adoption of habitual exercise in a significant portion of the population. Many of these individuals are defying the cultural norms by remaining physically active and competing at a high level in their later years. The juxtaposition between masters athletes and non-exercisers demonstrate the importance of remaining physically active throughout the lifespan on physiological systems related to healthspan (years of healthy living). This includes ~50% improved maximal aerobic capacity (VO₂max) and enhanced skeletal muscle health (size, function, as well as metabolic and communicative properties) compared to non-exercisers at a similar age. By taking a reductionist approach to VO₂max and skeletal muscle health, we can gain insight into how aging and habitual exercise affects the aging process. Collectively, this review provides a physiological basis for the elite performances seen in masters athletes, as well as the health implications of lifelong exercise with a focus on VO₂max, skeletal muscle metabolic fitness, whole muscle size and function, single muscle fiber physiology, and communicative properties of skeletal muscle. This review has significant public health implications due to the potent health benefits of habitual exercise across the lifespan.

Innovation in culture systems to study muscle complexity

Endogenous skeletal muscle development, regeneration, and pathology are extremely complex processes, influenced by local and systemic factors. Unpinning how these mechanisms function is crucial for fundamental biology and to develop therapeutic interventions for genetic disorders, but also conditions like sarcopenia and volumetric muscle loss. Ex vivo skeletal muscle models range from two- and three-dimensional primary cultures of satellite stem cell-derived myoblasts grown alone or in co-culture, to single muscle myofibers, myobundles, and whole tissues. Together, these systems provide the opportunity to gain mechanistic insights of stem cell behavior, cell-cell interactions, and mature muscle function in simplified systems, without confounding variables. Here, we highlight recent advances (published in the last 5 years) using in vitro primary cells and ex vivo skeletal muscle models, and summarize the new insights, tools, datasets, and screening methods they have provided. Finally, we highlight the opportunity for exponential advance of skeletal muscle knowledge, with spatiotemporal resolution, that is offered by guiding the study of muscle biology and physiology with in silico modelling and implementing high-content cell biology systems and ex vivo physiology platforms.

Methylome of skeletal muscle tissue in patients with hypertension and diabetes undergoing cardiopulmonary bypass

Background: Epigenomic changes occurring during surgery have been neglected in research; diabetes and hypertension can affect the epigenome but little is known about the epigenetics of skeletal muscle (SKM). Methods: DNA methylation was profiled via Illumina MethylationEPIC arrays in SKM samples obtained at the beginning and end of heart surgery with cardiopulmonary bypass. Results: Methylation in patients with hypertension and diabetes was significantly different, more so for uncontrolled diabetes; hypertension alone produced minimal effect. The affected pathways involved IL-1, IL-12, IL-18, TNF-α, IFN-γ, VEGF, NF-κB and Wnt signaling, apoptosis and DNA damage response. Significant changes occurred during surgery and included loci in the Hippo-YAP/TAZ pathway. Conclusion: Cardiopulmonary bypass surgery affects the SKM methylome, and the combination of hypertension and diabetes induces changes in the SKM epigenome in contrast to hypertension alone.

Genome-Wide Analysis of H3K27me3 in Porcine Embryonic Muscle Development

The trimethylation of histone H3 lysine 27 (H3K27me3) is one of the most important chromatin modifications, which is generally presented as a repressive mark in various biological processes. However, the dynamic and global-scale distribution of H3K27me3 during porcine embryonic muscle development remains unclear. Here, our study provided a comprehensive genome-wide view of H3K27me3 and analyzed the matching transcriptome in the skeletal muscles on days 33, 65, and 90 post-coitus from Duroc fetuses. Transcriptome analysis identified 4,124 differentially expressed genes (DEGs) and revealed the key transcriptional properties in three stages. We found that the global H3K27me3 levels continually increased during embryonic development, and the H3K27me3 level was negatively correlated with gene expression. The loss of H3K27me3 in the promoter was associated with the transcriptional activation of 856 DEGs in various processes, including skeletal muscle development, calcium signaling, and multiple metabolic pathways. We also identified for the first time that H3K27me3 could enrich in the promoter of genes, such as DES, MYL1, TNNC1, and KLF5, to negatively regulate gene expression in porcine satellite cells (PSCs). The loss of H3K27me3 could promote muscle cell differentiation. Taken together, this study provided the first genome-wide landscape of H3K27me3 in porcine embryonic muscle development. It revealed the complex and broad function of H3K27me3 in the regulation of embryonic muscle development from skeletal muscle morphogenesis to myofiber maturation.

Comparing the epigenetic landscape in myonuclei purified with a PCM1 antibody from a fast/glycolytic and a slow/oxidative muscle

Muscle cells have different phenotypes adapted to different usage, and can be grossly divided into fast/glycolytic and slow/oxidative types. While most muscles contain a mixture of such fiber types, we aimed at providing a genome-wide analysis of the epigenetic landscape by ChIP-Seq in two muscle extremes, the fast/glycolytic extensor digitorum longus (EDL) and slow/oxidative soleus muscles. Muscle is a heterogeneous tissue where up to 60% of the nuclei can be of a different origin. Since cellular homogeneity is critical in epigenome-wide association studies we developed a new method for purifying skeletal muscle nuclei from whole tissue, based on the nuclear envelope protein Pericentriolar material 1 (PCM1) being a specific marker for myonuclei. Using antibody labelling and a magnetic-assisted sorting approach, we were able to sort out myonuclei with 95% purity in muscles from mice, rats and humans. The sorting eliminated influence from the other cell types in the tissue and improved the myo-specific signal. A genome-wide comparison of the epigenetic landscape in EDL and soleus reflected the differences in the functional properties of the two muscles, and revealed distinct regulatory programs involving distal enhancers, including a glycolytic super-enhancer in the EDL. The two muscles were also regulated by different sets of transcription factors; e.g. in soleus, binding sites for MEF2C, NFATC2 and PPARA were enriched, while in EDL MYOD1 and SIX1 binding sites were found to be overrepresented. In addition, more novel transcription factors for muscle regulation such as members of the MAF family, ZFX and ZBTB14 were identified.

Complex tissues like skeletal muscle contain a variety of cells which confound the analysis of each cell type when based on homogenates, thus only about half of the cell nuclei in muscles reside inside the muscle cells. We here describe a labelling and sorting technique that allowed us to study the epigenetic landscape in purified muscle cell nuclei leaving the other cell types out. Differences between a fast/glycolytic and a slow/oxidative muscle were studied. While all skeletal muscle fibers have a similar make up and basic function, they differ in their physiology and the way they are used. Thus, some fibers are fast contracting but fatigable, and are used for short lasting explosive tasks such as sprinting. Other fibers are slow and are used for more prolonged tasks such as standing or long distance running. Since fiber type correlate with metabolic profile these features can also be related to metabolic diseases. We here show that the epigenetic landscape differed in gene loci corresponding to the differences in functional properties, and revealed that the two types are enriched in different gene regulatory networks. Exercise can alter muscle phenotype, and the epigenetic landscape might be related to how plastic different properties are.

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.

Molecular Profiling of DNA Methylation and Alternative Splicing of Genes in Skeletal Muscle of Obese Rabbits

DNA methylation and the alternative splicing of precursor messenger RNAs (pre-mRNAs) are two important genetic modification mechanisms. However, both are currently uncharacterized in the muscle metabolism of rabbits. Thus, we constructed the Tianfu black rabbit obesity model (obese rabbits fed with a 10% high-fat diet and control rabbits from 35 days to 70 days) and collected the skeletal muscle samples from the two groups for Genome methylation sequencing and RNA sequencing. DNA methylation data showed that the promoter regions of 599 genes and gene body region of 2522 genes had significantly differential methylation rates between the two groups, of which 288 genes had differential methylation rates in promoter and gene body regions. Analysis of alternative splicing showed 555 genes involved in exon skipping (ES) patterns, and 15 genes existed in differential methylation regions. Network analysis showed that 20 hub genes were associated with ubiquitinated protein degradation, muscle development pathways, and skeletal muscle energy metabolism. Our findings suggest that the two types of genetic modification have potential regulatory effects on skeletal muscle development and provide a basis for further mechanistic studies in the rabbit.

Human adipose and skeletal muscle tissue DNA, RNA, and protein content

The purpose of this project was to provide a profile of DNA, RNA, and protein content in adipose tissue, which is relatively understudied in humans, to gain more insight into the amount of tissue that may be required for various analyses. Skeletal muscle tissue was also investigated to provide a direct comparison into potential differences between these two highly metabolically active tissues. Basal adipose and skeletal muscle tissue samples were obtained from 10 (7 M, 3 W) recreationally active participants [25 ± 1 yr; 84 ± 3 kg, maximal oxygen consumption (V̇o_2max): 3.5 ± 0.2 L/min, body fat: 29 ± 2%]. DNA, RNA, and protein were extracted and subsequently analyzed for quantity and quality. DNA content of adipose and skeletal muscle tissue was 52 ± 14 and 189 ± 44 ng DNA·mg tissue^-1, respectively (P < 0.05). RNA content of adipose and skeletal muscle tissue was 46 ± 14 and 537 ± 72 ng RNA·mg tissue^-1, respectively (P < 0.05). Protein content of adipose and skeletal muscle tissue was 4 ± 1 and 177 ± 10 µg protein·mg tissue^-1, respectively (P < 0.05). In summary, human adipose had 28% of the DNA, 9% of the RNA, and 2% of the protein found in skeletal muscle per mg of tissue. This information should be useful across a wide range of human clinical investigation designs and various laboratory analyses.NEW & NOTEWORTHY This investigation studied DNA, RNA, and protein contents of adipose and skeletal muscle tissues from young active individuals. A series of optimization steps were investigated to aid in determining the optimal approach to extract high-yield and high-quality biomolecules. These findings contribute to the knowledge gap in adipose tissue requirements for molecular biology assays, which is of increasing importance due to the growing interest in adipose tissue research involving human exercise physiology research.

Genome-wide DNA methylation profiles provide insight into epigenetic regulation of red and white muscle development in Chinese perch Siniperca chuatsi

Fish skeletal muscles are composed of spatially well-separated fiber types, namely, red and white muscles with different physiological functions and metabolism. To compare the DNA methylation profiles of the two types of muscle tissues and identify potential candidate genes for the muscle growth and development under epigenetic regulation, genome-wide DNA methylation of the red and white muscle in Chinese perch Siniperca chuatsi were comparatively analyzed using bisulfate sequencing methods. An average of 0.9 billion 150-bp paired-end reads were obtained, of which 86% were uniquely mapped to the genome. Methylation mostly occurred at CG sites at a ratio of 94.43% in the red muscle and 93.16% in the white muscle. The mean methylation levels at C-sites were 5.95% in red muscle and 5.83% in white muscle, whereas the mean methylation levels of CG, CHG, and CHH were 73.23%, 0.62%, and 0.67% in red muscle, and 71.01%, 0.62%, and 0.67% in white muscle, respectively. A total of 4192 differentially methylated genes (DMGs) were identified significantly enriched in cell signaling pathways related to skeletal muscle differentiation and growth. Various muscle-related genes, including myosin gene isoforms and regulatory factors, are differentially methylated in the promoter region between the red and white muscles. Further analysis of the transcriptional expression of these genes showed that the muscle regulatory factors (myf5, myog, pax3, pax7, and twitst2) and myosin genes (myh10, myh16, myo18a, myo7a, myo9a, and myl3) were differentially expressed between the two kinds of muscles, consistent with the DNA methylation analysis results. ELISA assays confirmed that the level of 5mC in red muscle was significantly higher than in white muscle (P < 0.05). The RT-qPCR assays revealed that the expression levels of the three DNA methylation transferase (dnmt) subtypes, dnmt1, dnmt3ab, and dnmt3bb1, were significantly higher in red muscle than in white muscle. The higher DNA methylation levels in the red muscle may result from higher DNA methylation transferase expression in the red muscles. Thus, this study might provide a theoretical foundation to better understand epigenetic regulation in the growth and development of red and white muscles in animals, at least in Chinese perch fish.

Genetic and epigenetic regulation of skeletal muscle ribosome biogenesis with exercise

KEY POINTS

Ribosome biogenesis and MYC transcription are associated with acute resistance exercise (RE) and are distinct from endurance exercise in human skeletal muscle throughout a 24 h time course of recovery. A PCR-based method for relative ribosomal DNA (rDNA) copy number estimation was validated by whole genome sequencing and revealed that rDNA dosage is positively correlated with ribosome biogenesis in response to RE. Acute RE modifies rDNA methylation patterns in enhancer, intergenic spacer and non-canonical MYC-associated regions, but not the promoter. Myonuclear-specific rDNA methylation patterns with acute mechanical overload in mice corroborate and expand on rDNA findings with RE in humans. A genetic predisposition for hypertrophic responsiveness may exist based on rDNA gene dosage.

ABSTRACT

Ribosomes are the macromolecular engines of protein synthesis. Skeletal muscle ribosome biogenesis is stimulated by exercise, although the contribution of ribosomal DNA (rDNA) copy number and methylation to exercise-induced rDNA transcription is unclear. To investigate the genetic and epigenetic regulation of ribosome biogenesis with exercise, a time course of skeletal muscle biopsies was obtained from 30 participants (18 men and 12 women; 31 ± 8 years, 25 ± 4 kg m^-2 ) at rest and 30 min, 3 h, 8 h and 24 h after acute endurance (n = 10, 45 min cycling, 70% V ̇ O 2 max ) or resistance exercise (n = 10, 4 × 7 × 2 exercises); 10 control participants underwent biopsies without exercise. rDNA transcription and dosage were assessed using quantitative PCR and whole genome sequencing. rDNA promoter methylation was investigated using massARRAY EpiTYPER and global rDNA CpG methylation was assessed using reduced-representation bisulphite sequencing. Ribosome biogenesis and MYC transcription were associated primarily with resistance but not endurance exercise, indicating preferential up-regulation during hypertrophic processes. With resistance exercise, ribosome biogenesis was associated with rDNA gene dosage, as well as epigenetic changes in enhancer and non-canonical MYC-associated areas in rDNA, but not the promoter. A mouse model of in vivo metabolic RNA labelling and genetic myonuclear fluorescence labelling validated the effects of an acute hypertrophic stimulus on ribosome biogenesis and Myc transcription, and also corroborated rDNA enhancer and Myc-associated methylation alterations specifically in myonuclei. The present study provides the first information on skeletal muscle genetic and rDNA gene-wide epigenetic regulation of ribosome biogenesis in response to exercise, revealing novel roles for rDNA dosage and CpG methylation.

Difference in potential DNA methylation impact on gene expression between fast- and slow-type myofibers

Skeletal muscles are comprised of two major types of myofibers, fast and slow. It is hypothesized that once myofiber type is determined, muscle fiber-type specificity is maintained by an epigenetic mechanism, however, this remains poorly understood. To address this, we conducted a comprehensive CpG methylation analysis with a reduced representation of bisulfite sequencing (RRBS). Using GFP-myh7 mouse, we visually distinguished and separately pooled slow-type and myh7-negative fast-type fibers for analyses. A total of 31,967 and 26,274 CpGs were hypermethylated by ≥10% difference in the fast- and slow-type fibers, respectively. Notably, the number of promoter-hypermethylated genes with downregulated expression in the slow-type fibers was 3.5 times higher than that in the fast-type fibers. Gene bodies of the fast-type-specific myofibrillar genes Actn3, Tnnt3, Tnni2, Tnnc2, and Tpm1 were hypermethylated in the slow-type fibers, whereas those of the slow-type-specific genes Myh7, Tnnt1, and Tpm3 were hypermethylated in the fast-type fibers. Each of the instances of gene hypermethylation was associated with the respective downregulated expression. In particular, a relationship between CpG methylation sites and the transcription variant distribution of Tpm1 was observed, suggesting a regulation of Tpm1 alternative promoter usage by gene body CpG methylation. An association of hypermethylation with the regulation of gene expression was also observed in the transcription factors Sim2 and Tbx1. These results suggest not only a myofiber type-specific regulation of gene expression and alternative promoter usage by gene body CpG methylation but also a dominant effect of promoter-hypermethylation on the gene expressions in slow myofibers.

The myonuclear DNA methylome in response to an acute hypertrophic stimulus

In addition to multi-nucleated muscle fibres, numerous resident and infiltrating mononuclear cells populate the muscle compartment. As most epigenetic assays in skeletal muscle are conducted on whole tissue homogenates, essentially nothing is known about regulatory processes exclusively within muscle fibres in vivo. Utilizing a novel genetically modified mouse model developed by our laboratory, we (1) outline a simple and rapid workflow for isolating pure myonuclei from small tissue samples via fluorescent activated cell sorting and extracting high-quality large-fragment DNA for downstream analyses, and (2) provide information on myonuclear and interstitial cell nuclear CpG DNA methylation via reduced representation bisulphite sequencing (RRBS) using mice that were subjected to an acute mechanical overload of the plantaris muscle. In 3-month-old mice, myonuclei are ~50% of total nuclei in sham and ~30% in 3-d overloaded muscle, the difference being attributable to mononuclear cell infiltration and proliferation with overload. In purified myonuclei, pathway analysis of hypomethylated promoter regions following overload was distinct from interstitial nuclei and revealed marked regulation of factors that converge on the master regulator of muscle growth mTOR, and on autophagy. Specifically, acute hypomethylation of Rheb, Rictor, Hdac1, and Hdac2, in addition to a major driver of ribosome biogenesis Myc, reveals the epigenetic regulation of hypertrophic signalling within muscle fibres that may underpin the long-term growth response to loading. This study provides foundational information on global myonuclear epigenetics in vivo using RRBS, and demonstrates the importance of isolating specific nuclear populations to study the epigenetic regulation of skeletal muscle fibre adaptation.

Epigenetic Responses to Acute Resistance Exercise in Trained vs. Sedentary Men

Bagley, JR, Burghardt, KJ, McManus, R, Howlett, B, Costa, PB, Coburn, JW, Arevalo, JA, Malek, MH, and Galpin, AJ. Epigenetic responses to acute resistance exercise in trained vs. sedentary men. J Strength Cond Res 34(6): 1574-1580, 2020-Acute resistance exercise (RE) alters DNA methylation, an epigenetic process that influences gene expression and regulates skeletal muscle adaptation. This aspect of cellular remodeling is poorly understood, especially in resistance-trained (RT) individuals. The study purpose was to examine DNA methylation in response to acute RE in RT and sedentary (SED) young men, specifically targeting genes responsible for metabolic, inflammatory, and hypertrophic muscle adaptations. Vastus lateralis biopsies were performed before (baseline), 30 minutes after, and 4 hours after an acute RE bout (3 × 10 repetitions at 70% 1 repetition maximum [1RM] leg press and leg extension) in 11 RT (mean ± SEM: age = 26.1 ± 1.0 years; body mass = 84.3 ± 0.2 kg; leg press 1RM = 412.6 ± 25.9 kg) and 8 SED (age = 22.9 ± 1.1 years; body mass = 75.6 ± 0.3 kg; leg press 1RM = 164.8 ± 22.5 kg) men. DNA methylation was analyzed through methylation sensitive high-resolution melting using real-time polymerase chain reaction. Separate 2 (group) × 3 (time) repeated-measures analyses of variance and analyses of covariance were performed to examine changes in DNA methylation for each target gene. Results showed that acute RE (a) hypomethylated LINE-1 (measure of global methylation) in RT but not SED, (b) hypermethylated metabolic genes (GPAM and SREBF2) in RT, while lowering SREBF2 methylation in SED, and (c) did not affect methylation of genes associated with inflammation (IL-6 and TNF-α) or hypertrophy (mTOR and AKT1). However, basal IL-6 and TNF-α were lower in SED compared with RT. These findings indicate the same RE stimulus can illicit different epigenetic responses in RT vs. SED men and provides a molecular mechanism underpinning the need for differential training stimuli based on subject training backgrounds.

Exertional heat stroke leads to concurrent long-term epigenetic memory, immunosuppression and altered heat shock response in female mice

KEY POINTS

Exposure to exertional heat stroke (EHS) has been linked to increased long-term decrements of health. Epigenetic reprogramming is involved in the response to heat acclimation; however, whether the long-term effects of EHS are mediated by epigenetic reprogramming is unknown. In female mice, we observed DNA methylation reprogramming in bone marrow-derived (BMD) monocytes as early as 4 days of recovery from EHS and as late as 30 days compared with sham exercise controls. Whole blood, collected after 30 days of recovery from EHS, exhibited an immunosuppressive phenotype when challenged in vitro by lipopolysaccharide. After 30 days of recovery from EHS, BMD monocytes exhibited an altered in vitro heat shock response. The location of differentially methylated CpGs are predictive of both the immunosuppressive phenotype and altered heat shock responses.

ABSTRACT

Exposure to exertional heat stroke (EHS) has been linked to increased susceptibility to a second heat stroke, infection and cardiovascular disease. Whether these clinical outcomes are mediated by an epigenetic memory is unknown. Using a preclinical mouse model of EHS, we investigated whether EHS exposure produces a lasting epigenetic memory in monocytes and whether there are phenotypic alterations that may be consistent with these epigenetic changes. Female mice underwent forced wheel running at 37.5°C/40% relative humidity until symptom limitation, characterized by CNS dysfunction. Results were compared with matched exercise controls at 22.5°C. Monocytes were isolated from bone marrow after 4 or 30 days of recovery to extract DNA and analyse methylation. Broad-ranging alterations to the DNA methylome were observed at both time points. At 30 days, very specific alterations were observed to the promoter regions of genes involved with immune responsiveness. To test whether these changes might be related to phenotype, whole blood at 30 days was challenged with lipopolysaccharide (LPS) to measure cytokine secretion; monocytes were also challenged with heat shock to quantify mRNA expression. Whole blood collected from EHS mice showed markedly attenuated inflammatory responses to LPS challenge. Furthermore, monocyte mRNA from EHS mice showed significantly altered responses to heat shock challenge. These results demonstrate that EHS leads to a unique DNA methylation pattern in monocytes and altered immune and heat shock responsiveness after 30 days. These data support the hypothesis that EHS exposure can induce long-term physiological changes that may be linked to altered epigenetic profiles.

Transcriptome and epigenome diversity and plasticity of muscle stem cells following transplantation

Adult skeletal muscles are maintained during homeostasis and regenerated upon injury by muscle stem cells (MuSCs). A heterogeneity in self-renewal, differentiation and regeneration properties has been reported for MuSCs based on their anatomical location. Although MuSCs derived from extraocular muscles (EOM) have a higher regenerative capacity than those derived from limb muscles, the molecular determinants that govern these differences remain undefined. Here we show that EOM and limb MuSCs have distinct DNA methylation signatures associated with enhancers of location-specific genes, and that the EOM transcriptome is reprogrammed following transplantation into a limb muscle environment. Notably, EOM MuSCs expressed host-site specific positional Hox codes after engraftment and self-renewal within the host muscle. However, about 10% of EOM-specific genes showed engraftment-resistant expression, pointing to cell-intrinsic molecular determinants of the higher engraftment potential of EOM MuSCs. Our results underscore the molecular diversity of distinct MuSC populations and molecularly define their plasticity in response to microenvironmental cues. These findings provide insights into strategies designed to improve the functional capacity of MuSCs in the context of regenerative medicine.

Selection signatures of Fuzhong Buffalo based on whole-genome sequences

BACKGROUND

Fuzhong buffalo, a native breed of Guangxi Zhuang Autonomous Region, is traditionally used as a draft animal to provide farm power in the rice cultivation. In addition, the Fuzhong buffalo also prepared for the bullfighting festival organized by the locals. The detection of the selective signatures in its genome can help in elucidating the selection mechanisms in its stamina and muscle development of a draft animal.

RESULTS

In this study, we analyzed 27 whole genomes of buffalo (including 15 Fuzhong buffalo genomes and 12 published buffalo genomes from Upper Yangtze region). The ZHp, ZFst, π-Ratio, and XP-EHH statistics were used to identify the candidate signatures of positive selection in Fuzhong buffalo. Our results detected a set of candidate genes involving in the pathways and GO terms associated with the response to exercise (e.g., ALDOA, STAT3, AKT2, EIF4E2, CACNA2D2, TCF4, CDH2), immunity (e.g., PTPN22, NKX2-3, PIK3R1, ITK, TMEM173), nervous system (e.g., PTPN21, ROBO1, HOMER1, MAGI2, SLC1A3, NRG3, SNAP47, CTNNA2, ADGRL3). In addition, we also identified several genes related to production and growth traits (e.g., PHLPP1, PRKN, MACF1, UCN3, RALGAPA1, PHKB, PKD1L). Our results depicted several pathways, GO terms, and candidate genes to be associated with response to exercise, immunity, nervous system, and growth traits.

CONCLUSIONS

The selective sweep analysis of the Fuzhong buffalo demonstrated positive selection pressure on potential target genes involved in behavior, immunity, and growth traits, etc. Our findings provided a valuable resource for future research on buffalo breeding and an insight into the mechanisms of artificial selection.

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.

Adaptive responses of histone modifications to resistance exercise in human skeletal muscle

Exercise training causes epigenetic changes in skeletal muscle, although it is unclear how resistance exercise (RE) affects histone modifications. The present study was carried out to investigate the effects of acute RE and RE training on gene expression profiles and histone modifications in human skeletal muscle. Healthy male adults were assigned to acute RE (n = 9, age = 20.5±4.3yr, BMI = 28.0±6.8kg/m2) or RE training (n = 21, age = 23.7±2.5yr, BMI = 24.2±2.7kg/m2) groups. Biopsy samples were obtained from the vastus lateralis muscle before and three hours after a single bout of acute RE, or 3-days after 10 weeks of RE training. RNA sequencing analysis revealed that 153 genes with GO terms including muscle development, stress response, metabolism, cell death, and transcription factor were significantly up-regulated (+291% vs. pre-acute RE) upon acute RE. Expressions of these genes were also greater (+9.6% vs. pre-RE training, p<0.05) in RE trained subjects. Significant up-regulation of acetylated histone 3 (H3) (+235%) and H3 mono-methylated at lysine 4 (+290%) and tri-methylated at lysine 27 (+849%), whereas down-regulation of H3.3 variant (-39%) distributions relative to total H3 were observed at transcriptionally activated loci after acute RE compared to pre-acute RE levels. Interestingly, the distribution of acetylated H3 was found to be up-regulated as compared to the level of total H3 after RE training (+40%, p<0.05). These results indicate that a single bout of RE drastically alters both gene expressions and histone modifications in human skeletal muscle. It is also suggested that enhanced histone acetylation is closely related to up-regulation of gene expressions after RE training.

Epigenetics of Skeletal Muscle-Associated Genes in the ASB, LRRC, TMEM, and OSBPL Gene Families

Much remains to be discovered about the intersection of tissue-specific transcription control and the epigenetics of skeletal muscle (SkM), a very complex and dynamic organ. From four gene families, Leucine-Rich Repeat Containing (LRRC), Oxysterol Binding Protein Like (OSBPL), Ankyrin Repeat and Socs Box (ASB), and Transmembrane Protein (TMEM), we chose 21 genes that are preferentially expressed in human SkM relative to 52 other tissue types and analyzed relationships between their tissue-specific epigenetics and expression. We also compared their genetics, proteomics, and descriptions in the literature. For this study, we identified genes with little or no previous descriptions of SkM functionality (ASB4, ASB8, ASB10, ASB12, ASB16, LRRC14B, LRRC20, LRRC30, TMEM52, TMEM233, OSBPL6/ORP6, and OSBPL11/ORP11) and included genes whose SkM functions had been previously addressed (ASB2, ASB5, ASB11, ASB15, LRRC2, LRRC38, LRRC39, TMEM38A/TRIC-A, and TMEM38B/TRIC-B). Some of these genes have associations with SkM or heart disease, cancer, bone disease, or other diseases. Among the transcription-related SkM epigenetic features that we identified were: super-enhancers, promoter DNA hypomethylation, lengthening of constitutive low-methylated promoter regions, and SkM-related enhancers for one gene embedded in a neighboring gene (e.g., ASB8-PFKM, LRRC39-DBT, and LRRC14B-PLEKHG4B gene-pairs). In addition, highly or lowly co-expressed long non-coding RNA (lncRNA) genes probably regulate several of these genes. Our findings give insights into tissue-specific epigenetic patterns and functionality of related genes in a gene family and can elucidate normal and disease-related regulation of gene expression in SkM.

Transcriptional memory in skeletal muscle. Don't forget (to) exercise

Transcriptional memory describes an ancient and highly conserved form of cellular learning that enables cells to benefit from recent experience by retaining a mitotically inheritable but reversible memory of the initial transcriptional response when encountering an environmental or physiological stimulus. Herein, we will review recent progress made in the understanding of how cells can make use of diverse constituents of the epigenetic toolbox to retain a transcriptional memory of past states and perturbations. Specifically, we will outline how these mechanisms will help to improve our understanding of skeletal muscle plasticity in health and disease. We describe the epigenetic road map that allows skeletal muscle fibers to navigate through training-induced adaptation processes, and how epigenetic memory marks can preserve an autobiographical history of lifestyle behavior changes, pathological challenges, and aging. We will further consider some key findings in the field of exercise epigenomics to emphasize major challenges when interpreting dynamic changes in the chromatin landscape in response to acute exercise and training.

Single-cell transcriptional profiles in human skeletal muscle

Skeletal muscle is a heterogeneous tissue comprised of muscle fiber and mononuclear cell types that, in addition to movement, influences immunity, metabolism and cognition. We investigated the gene expression patterns of skeletal muscle cells using RNA-seq of subtype-pooled single human muscle fibers and single cell RNA-seq of mononuclear cells from human vastus lateralis, mouse quadriceps, and mouse diaphragm. We identified 11 human skeletal muscle mononuclear cell types, including two fibro-adipogenic progenitor (FAP) cell subtypes. The human FBN1+ FAP cell subtype is novel and a corresponding FBN1+ FAP cell type was also found in single cell RNA-seq analysis in mouse. Transcriptome exercise studies using bulk tissue analysis do not resolve changes in individual cell-type proportion or gene expression. The cell-type gene signatures provide the means to use computational methods to identify cell-type level changes in bulk studies. As an example, we analyzed public transcriptome data from an exercise training study and revealed significant changes in specific mononuclear cell-type proportions related to age, sex, acute exercise and training. Our single-cell expression map of skeletal muscle cell types will further the understanding of the diverse effects of exercise and the pathophysiology of muscle disease.

Demethylation and derepression of genomic retroelements in the skeletal muscles of aged mice

Changes in DNA methylation influence the aging process and contribute to aging phenotypes, but few studies have been conducted on DNA methylation changes in conjunction with skeletal muscle aging. We explored the DNA methylation changes in a variety of retroelement families throughout aging (at 2, 20, and 28 months of age) in murine skeletal muscles by methyl‐binding domain sequencing (MBD‐seq). The two following contrasting patterns were observed among the members of each repeat family in superaged mice: (a) hypermethylation in weakly methylated retroelement copies and (b) hypomethylation in copies with relatively stronger methylation levels, representing a pattern of “regression toward the mean” within a single retroelement family. Interestingly, these patterns depended on the sizes of the copies. While the majority of the elements showed a slight increase in methylation, the larger copies (>5 kb) displayed evident demethylation. All these changes were not observed in T cells. RNA sequencing revealed a global derepression of retroelements during the late phase of aging (between 20 and 28 months of age), which temporally coincided with retroelement demethylation. Following this methylation drift trend of “regression toward the mean,” aging tended to progressively lose the preexisting methylation differences and local patterns in the genomic regions that had been elaborately established during the early period of development.

One week, but not 12 hours, of cast immobilization alters promotor DNA methylation patterns in the nNOS gene in mouse skeletal muscle

KEY POINTS

DNA methylation may play an important role in regulating gene expression in skeletal muscle to adapt to physical activity and inactivity. Neuronal nitric oxide synthase (nNOS) in skeletal muscle is a key regulator of skeletal muscle mass; however, it is unclear whether nNOS expression is regulated by DNA methylation. We found that 1 week of cast immobilization increased nNOS DNA methylation levels and downregulated nNOS gene expression in atrophic slow-twitch soleus muscle from the mouse leg. These changes were not detected in non-atrophic fast-twitch extensor digitorum longus muscle. Twelve hours of cast immobilization decreased nNOS gene expression, whereas nNOS DNA methylation levels were unchanged, suggesting that downregulation of nNOS gene expression by short-term muscle inactivity is independent of the DNA methylation pattern. These findings contribute to a better understanding of the maintenance of skeletal muscle mass and prevention of muscle atrophy by epigenetic mechanisms via the nNOS/NO pathway.

ABSTRACT

DNA methylation is a mechanism that controls gene expression in skeletal muscle under various environmental stimuli, such as physical activity and inactivity. Neuronal nitric oxide synthase (nNOS) regulates muscle atrophy in skeletal muscle. However, the mechanisms regulating nNOS expression in atrophic muscle remain unclear. We hypothesized that nNOS expression in atrophic muscle is regulated by DNA methylation of the nNOS promotor in soleus (Sol; slow-twitch fibre dominant) and extensor digitorum longus (EDL; fast-twitch fibre dominant) muscles. One week of cast immobilization induced significant muscle atrophy in Sol but not in EDL. We showed that 1 week of cast immobilization increased nNOS DNA methylation levels in Sol, although only a minor change was detected in EDL. Consistent with the increased DNA methylation levels in atrophic Sol, the gene expression levels of total nNOS and nNOSµ (i.e. the major splicing variant of nNOS in skeletal muscle) decreased. The abundance of the nNOS protein and cell membrane (especially type IIa fibre) immunoreactivity also decreased in atrophic Sol. These changes were not observed in EDL after 1 week of cast immobilization. Furthermore, despite the lack of significant atrophy, 12 h of cast immobilization decreased gene expression levels of total nNOS and nNOSµ in Sol. However, no association was detected between nNOS DNA methylation and gene expression. The expression of the nNOSβ gene, another splicing variant of nNOS, in EDL was unchanged by cast immobilization, whereas its expression was not detected in Sol. We concluded that chronic adaptation of nNOS gene expression in cast immobilized muscle may involve nNOS DNA methylation.

Transcriptional regulatory region and DNA methylation analysis of TNNI1 gene promoters in Gaoyou duck skeletal muscle (Anas platyrhynchos domestica)

1. The slow skeletal muscle troponin I (TNNI1) gene has been found to be specifically expressed in slow muscle fibres and plays an important role in muscle development. The aim of this study was to determine the active control area of duck TNNI1 and identify the potential cis-regulatory elements in the promoter. 2. In this study, the TNNI1 promoter was first cloned by genome walking and the sequences were analysed using bioinformatics software. Firefly luciferase reporter gene vectors, driven by a series of constructs with progressive deletions, were used to identify the core transcriptional regulatory region of the duck TNNI1 gene. The methylation status of the CpG island in the TNNI1 promoter was detected in skeletal muscle on embryonic days 21 and 27, by bisulphite sequencing PCR (BSP). 3. The results showed two CpG islands presented in the promoter region, with one of the CpG islands located in the core transcriptional regulatory region (-2078/-885 bp). The total methylation levels of the 14 CpG sites were not altered between breast and leg muscles on embryonic days 21 and 27. However, four CpG sites (loci of positions 4, 11, 13, and 14) showed dramatically different methylation levels between breast and leg muscles at embryonic days 21 and 27. Analysis showed that multiple CpG sites had a significant correlation between the methylation levels of the CpG sites and mRNA expressions in skeletal muscle. Multiple transcription factor binding sites including Sp1, c-Myc, Oct-1 and NF-kB motifs were identified and might be responsible for transcriptional regulation of the TNNI1 gene. 4. These findings contribute to further understanding of the fundamental mechanism for transcriptional regulation of the TNNI1 gene in ducks.

Genome-wide methylation analysis reveals differentially methylated loci that are associated with an age-dependent increase in bovine fibroblast response to LPS

BACKGROUND

Differences in DNA methylation are known to contribute to the development of immune-related disorders in humans but relatively little is known about how methylation regulates immune function in cattle. Utilizing whole-transcriptome analyses of bovine dermal fibroblasts, we have previously identified an age and breed-dependent up-regulation of genes within the toll-like receptor 4 (TLR4) pathway that correlates with enhanced fibroblast production of IL-8 in response to lipopolysaccharide (LPS). Age-dependent differences in IL-8 production are abolished by treatment with 5-aza-2-deoxycytidine and Trichostatin A (AZA-TSA), suggesting epigenetic regulation of the innate response to LPS. In the current study, we performed reduced representation bisulfite sequencing (RRBS) on fibroblast cultures isolated from the same animals at 5- and 16-months of age to identify genes that exhibit variable methylation with age. To validate the role of methylation in gene expression, six innate response genes that were hyper-methylated in young animals were assessed by RT-qPCR in fibroblasts from animals at different ages and from different breeds.

RESULTS

We identified 14,094 differentially methylated CpGs (DMCs) that differed between fibroblast cultures at 5- versus 16-months of age. Of the 5065 DMCs that fell within gene regions, 1117 were located within promoters, 1057 were within gene exons and 2891 were within gene introns and 67% were more methylated in young cultures. Transcription factor enrichment of the promoter regions hyper-methylated in young cultures revealed significant regulation by the key pro-inflammatory regulator, NF-κB. Additionally, five out of six chosen genes (PIK3R1, FES, NFATC1, TNFSF13 and RORA) that were more methylated in young cultures showed a significant reduction in expression post-LPS treatment in comparison with older cultures. Two of these genes, FES and NFATC1, were similarly down-regulated in Angus cultures that also exhibit a low LPS response phenotype.

CONCLUSIONS

Our study has identified immune-related loci regulated by DNA methylation in cattle that may contribute to differential cellular response to LPS, two of which exhibit an identical expression profile in both low-responding age and breed phenotypes. Methylation biomarkers of differential immunity may prove useful in developing selection strategies for replacement cows that are less susceptible to severe infections, such as coliform mastitis.