Dynamic enhancers control skeletal muscle identity and reprogramming

Sucrose overconsumption impairs AgRP neuron dynamics and promotes palatable food intake

Rapid gut-brain communication is critical to maintain energy balance and is disrupted in diet-induced obesity. In particular, the role of carbohydrate overconsumption in the regulation of interoceptive circuits in vivo requires further investigation. Here, we report that an obesogenic high-sucrose diet (HSD) selectively blunts silencing of hunger-promoting agouti-related protein (AgRP) neurons following intragastric delivery of glucose, whereas we previously showed that overconsumption of a high-fat diet (HFD) selectively attenuates lipid-induced neural silencing. By contrast, both HSD and HFD reversibly dampen rapid AgRP neuron inhibition following chow presentation and promote intake of more palatable foods. Our findings reveal that excess sugar and fat pathologically modulate feeding circuit activity in both macronutrient-dependent and -independent ways and thus may additively exacerbate obesity.

Transcriptional programming of translation by BCL6 controls skeletal muscle proteostasis

Skeletal muscle is dynamically controlled by the balance of protein synthesis and degradation. Here we discover an unexpected function for the transcriptional repressor B cell lymphoma 6 (BCL6) in muscle proteostasis and strength in mice. Skeletal muscle-specific Bcl6 ablation in utero or in adult mice results in over 30% decreased muscle mass and force production due to reduced protein synthesis and increased autophagy, while it promotes a shift to a slower myosin heavy chain fibre profile. Ribosome profiling reveals reduced overall translation efficiency in Bcl6-ablated muscles. Mechanistically, tandem chromatin immunoprecipitation, transcriptomic and translational analyses identify direct BCL6 repression of eukaryotic translation initiation factor 4E-binding protein 1 (Eif4ebp1) and activation of insulin-like growth factor 1 (Igf1) and androgen receptor (Ar). Together, these results uncover a bifunctional role for BCL6 in the transcriptional and translational control of muscle proteostasis.

Molecular control of endurance training adaptation in male mouse skeletal muscle

Skeletal muscle has an enormous plastic potential to adapt to various external and internal perturbations. Although morphological changes in endurance-trained muscles are well described, the molecular underpinnings of training adaptation are poorly understood. We therefore aimed to elucidate the molecular signature of muscles of trained male mice and unravel the training status-dependent responses to an acute bout of exercise. Our results reveal that, even though at baseline an unexpectedly low number of genes define the trained muscle, training status substantially affects the transcriptional response to an acute challenge, both quantitatively and qualitatively, in part associated with epigenetic modifications. Finally, transiently activated factors such as the peroxisome proliferator-activated receptor-γ coactivator 1α are indispensable for normal training adaptation. Together, these results provide a molecular framework of the temporal and training status-dependent exercise response that underpins muscle plasticity in training.

Multi-omic identification of key transcriptional regulatory programs during endurance exercise training

Transcription factors (TFs) play a key role in regulating gene expression and responses to stimuli. We conducted an integrated analysis of chromatin accessibility, DNA methylation, and RNA expression across eight rat tissues following endurance exercise training (EET) to map epigenomic changes to transcriptional changes and determine key TFs involved. We uncovered tissue-specific changes and TF motif enrichment across all omic layers, differentially accessible regions (DARs), differentially methylated regions (DMRs), and differentially expressed genes (DEGs). We discovered distinct routes of EET-induced regulation through either epigenomic alterations providing better access for TFs to affect target genes, or via changes in TF expression or activity enabling target gene response. We identified TF motifs enriched among correlated epigenomic and transcriptomic alterations, DEGs correlated with exercise-related phenotypic changes, and EET-induced activity changes of TFs enriched for DEGs among their gene targets. This analysis elucidates the unique transcriptional regulatory mechanisms mediating diverse organ effects of EET.

Reprogramming of cis-regulatory networks during skeletal muscle atrophy in male mice

A comprehensive atlas of cis-regulatory elements and their dynamic activity is necessary to understand the transcriptional basis of cellular structure maintenance, metabolism, and responses to the environment. Here we show, using matched single-nucleus chromatin accessibility and RNA-sequencing from juvenile male C57BL6 mice, an atlas of accessible chromatin regions in both normal and denervated skeletal muscles. We identified cell-type-specific cis-regulatory networks, highlighting the dynamic regulatory circuits mediating transitions between myonuclear types. Through comparison of normal and perturbed muscle, we delineated the reprogramming of cis-regulatory networks in response to denervation, described the interplay of promoters/enhancers and target genes. We further unveil a hierarchical structure of transcription factors that delineate a regulatory network in atrophic muscle, identifying ELK4 as a key atrophy-related transcription factor that instigates muscle atrophy through TGF-β1 regulation. This study furnishes a rich genomic resource, essential for decoding the regulatory dynamics of skeletal muscle in both physiological and pathological states.

Integrated single-cell multiome analysis reveals muscle fiber-type gene regulatory circuitry modulated by endurance exercise

Endurance exercise is an important health modifier. We studied cell-type specific adaptations of human skeletal muscle to acute endurance exercise using single-nucleus (sn) multiome sequencing in human vastus lateralis samples collected before and 3.5 hours after 40 min exercise at 70% VO2max in four subjects, as well as in matched time of day samples from two supine resting circadian controls. High quality same-cell RNA-seq and ATAC-seq data were obtained from 37,154 nuclei comprising 14 cell types. Among muscle fiber types, both shared and fiber-type specific regulatory programs were identified. Single-cell circuit analysis identified distinct adaptations in fast, slow and intermediate fibers as well as LUM-expressing FAP cells, involving a total of 328 transcription factors (TFs) acting at altered accessibility sites regulating 2,025 genes. These data and circuit mapping provide single-cell insight into the processes underlying tissue and metabolic remodeling responses to exercise.

Nuclear receptors as potential therapeutic targets in peripheral arterial disease and related myopathy

Peripheral arterial disease (PAD) is a prevalent cardiovascular complication of limb vascular insufficiency, causing ischemic injury, mitochondrial metabolic damage and functional impairment in the skeletal muscle, and ultimately leading to immobility and mortality. While potential therapies have been mostly focussed on revascularization, none of the currently available pharmacological treatments are fully effective in PAD, often leading to amputations, particularly in chronic metabolic diseases. One major limitation of focussed angiogenesis and revascularization as a therapeutic strategy is a limited effect on metabolic restoration and muscle regeneration in the affected limb. Therefore, additional preclinical investigations are needed to discover novel treatment options for PAD preferably targeting multiple aspects of muscle recovery. In this review, we propose nuclear receptors expressed in the skeletal muscle as potential candidates for ischemic muscle repair in PAD. We review classic steroid and orphan receptors that have been reported to be involved in the regulation of paracrine muscle angiogenesis, oxidative metabolism, mitochondrial biogenesis and muscle regeneration, and discuss how these receptors could be critical for recovery from ischemic muscle damage. Furthermore, we identify existing gaps in our understanding of nuclear receptor signalling in the skeletal muscle and propose future areas of research that could be instrumental in exploring nuclear receptors as therapeutic candidates for treating PAD.

Exosome-derived circKIF20B suppresses gefitinib resistance and cell proliferation in non-small cell lung cancer

BACKGROUND

The gefitinib resistance mechanism in non-small cell lung cancer (NSCLC) remains unclear, albeit exosomal circular RNA (circRNA) is known to possibly play a vital role in it.

METHODS

We employed high-throughput sequencing techniques to detect the expressions of exosomal circRNA both in gefitinib-resistant and gefitinib-sensitive cells in this study. The circKIF20B expression was determined in serum exosomes and tissues of patients by qRT-PCR. The structure, stability, and intracellular localization of circKIF20B were verified by Sanger sequencing, Ribonuclease R (RNase R)/actinomycin D (ACTD) treatments, and Fluorescence in situ hybridization (FISH). The functions of circKIF20B were investigated by 5-Ethynyl-20-deoxyuridine (EdU), flow cytometry, Cell Counting Kit-8 (CCK-8), oxygen consumption rate (OCR), and xenograft model. Co-culture experiments were performed to explore the potential ability of exosomal circKIF20B in treating gefitinib resistance. The downstream targets of circKIF20B were determined by luciferase assay, RNA pulldown, and RNA immunoprecipitation (RIP).

RESULTS

We found that circKIF20B was poorly expressed in the serum exosomes of gefitinib-resistant patients (n = 24) and the tumor tissues of patients with NSCLC (n = 85). CircKIF20B was negatively correlated with tumor size and tumor stage. Decreasing circKIF20B was found to promote gefitinib resistance by accelerating the cell cycle, inhibiting apoptosis, and enhancing mitochondrial oxidative phosphorylation (OXPHOS), whereas increasing circKIF20B was found to restore gefitinib sensitivity. Mechanistically, circKIF20B is bound to miR-615-3p for regulating the MEF2A and then altering the cell cycle, apoptosis, and mitochondrial OXPHOS. Overexpressing circKIF20B parental cells can restore sensitivity to gefitinib in the recipient cells by upregulating the exosomal circKIF20B expression.

CONCLUSIONS

This study revealed a novel mechanism of circKIF20B/miR-615-3p/MEF2A signaling axis involving progression of gefitinib resistance in NSCLC. Exosomal circKIF20B is expected to be an easily accessible and alternative liquid biopsy candidate and potential therapeutic target in gefitinib-resistant NSCLC. The schematic diagram of mechanism in this study. Exosomal circKIF20B inhibits gefitinib resistance and cell proliferation by arresting the cell cycle, promoting apoptosis, and reducing OXPHOS via circKIF20B/miR-615-3p/MEF2A axis in NSCLC.

Dynamic chromatin architectures provide insights into the genetics of cattle myogenesis

BACKGROUND

Sharply increased beef consumption is propelling the genetic improvement projects of beef cattle in China. Three-dimensional genome structure is confirmed to be an important layer of transcription regulation. Although genome-wide interaction data of several livestock species have already been produced, the genome structure states and its regulatory rules in cattle muscle are still limited.

RESULTS

Here we present the first 3D genome data in Longissimus dorsi muscle of fetal and adult cattle (Bos taurus). We showed that compartments, topologically associating domains (TADs), and loop undergo re-organization and the structure dynamics were consistent with transcriptomic divergence during muscle development. Furthermore, we annotated cis-regulatory elements in cattle genome during myogenesis and demonstrated the enrichments of promoter and enhancer in selection sweeps. We further validated the regulatory function of one HMGA2 intronic enhancer near a strong sweep region on primary bovine myoblast proliferation.

CONCLUSIONS

Our data provide key insights of the regulatory function of high order chromatin structure and cattle myogenic biology, which will benefit the progress of genetic improvement of beef cattle.

Effect of 11-Deoxycorticosterone in the Transcriptomic Response to Stress in Rainbow Trout Skeletal Muscle

In aquaculture, many stressors can negatively affect growth in teleosts. It is believed that cortisol performs glucocorticoid and mineralocorticoid functions because teleosts do not synthesize aldosterone. However, recent data suggest that 11-deoxycorticosterone (DOC) released during stress events may be relevant to modulate the compensatory response. To understand how DOC modifies the skeletal muscle molecular response, we carried out a transcriptomic analysis. Rainbow trout (Oncorhynchus mykiss) were intraperitoneally treated with physiological doses of DOC in individuals pretreated with mifepristone (glucocorticoid receptor antagonist) or eplerenone (mineralocorticoid receptor antagonist). RNA was extracted from the skeletal muscles, and cDNA libraries were constructed from vehicle, DOC, mifepristone, mifepristone plus DOC, eplerenone, and eplerenone plus DOC groups. The RNA-seq analysis revealed 131 differentially expressed transcripts (DETs) induced by DOC with respect to the vehicle group, mainly associated with muscle contraction, sarcomere organization, and cell adhesion. In addition, a DOC versus mifepristone plus DOC analysis revealed 122 DETs related to muscle contraction, sarcomere organization, and skeletal muscle cell differentiation. In a DOC versus eplerenone plus DOC analysis, 133 DETs were associated with autophagosome assembly, circadian regulation of gene expression, and regulation of transcription from RNA pol II promoter. These analyses indicate that DOC has a relevant function in the stress response of skeletal muscles, whose action is differentially modulated by GR and MR and is complementary to cortisol.

LSD1 defines the fiber type-selective responsiveness to environmental stress in skeletal muscle

Skeletal muscle exhibits remarkable plasticity in response to environmental cues, with stress-dependent effects on the fast-twitch and slow-twitch fibers. Although stress-induced gene expression underlies environmental adaptation, it is unclear how transcriptional and epigenetic factors regulate fiber type-specific responses in the muscle. Here, we show that flavin-dependent lysine-specific demethylase-1 (LSD1) differentially controls responses to glucocorticoid and exercise in postnatal skeletal muscle. Using skeletal muscle-specific LSD1-knockout mice and in vitro approaches, we found that LSD1 loss exacerbated glucocorticoid-induced atrophy in the fast fiber-dominant muscles, with reduced nuclear retention of Foxk1, an anti-autophagic transcription factor. Furthermore, LSD1 depletion enhanced endurance exercise-induced hypertrophy in the slow fiber-dominant muscles, by induced expression of ERRγ, a transcription factor that promotes oxidative metabolism genes. Thus, LSD1 serves as an 'epigenetic barrier' that optimizes fiber type-specific responses and muscle mass under the stress conditions. Our results uncover that LSD1 modulators provide emerging therapeutic and preventive strategies against stress-induced myopathies such as sarcopenia, cachexia, and disuse atrophy.

Identification of enhancers responsible for the coordinated expression of myosin heavy chain isoforms in skeletal muscle

Background

Skeletal muscles consist of fibers of differing contractility and metabolic properties, which are primarily determined by the content of myosin heavy chain (MYH) isoforms (MYH7, MYH2, MYH1, and MYH4). The regulation of Myh genes transcription depends on three-dimensional chromatin conformation interaction, but the mechanistic details remain to be determined.

Results

In this study, we characterized the interaction profiles of Myh genes using 4C-seq (circular chromosome conformation capture coupled to high-throughput sequencing). The interaction profile of Myh genes changed between fast quadriceps and slow soleus muscles. Combining chromatin immunoprecipitation-sequencing (ChIP-seq) and transposase accessible chromatin with high-throughput sequencing (ATAC-seq), we found that a 38 kb intergenic region interacting simultaneously with fast Myh genes promoters controlled the coordinated expression of fast Myh genes. We also identified four active enhancers of Myh7, and revealed that binding of MYOG and MYOD increased the activity of Myh7 enhancers.

Conclusions

This study provides new insight into the chromatin interactions that regulate Myh genes expression.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08737-9.

Ranking reprogramming factors for cell differentiation

Transcription factor over-expression is a proven method for reprogramming cells to a desired cell type for regenerative medicine and therapeutic discovery. However, a general method for the identification of reprogramming factors to create an arbitrary cell type is an open problem. Here we examine the success rate of methods and data for differentiation by testing the ability of nine computational methods (CellNet, GarNet, EBseq, AME, DREME, HOMER, KMAC, diffTF and DeepAccess) to discover and rank candidate factors for eight target cell types with known reprogramming solutions. We compare methods that use gene expression, biological networks and chromatin accessibility data, and comprehensively test parameter and preprocessing of input data to optimize performance. We find the best factor identification methods can identify an average of 50-60% of reprogramming factors within the top ten candidates, and methods that use chromatin accessibility perform the best. Among the chromatin accessibility methods, complex methods DeepAccess and diffTF have higher correlation with the ranked significance of transcription factor candidates within reprogramming protocols for differentiation. We provide evidence that AME and diffTF are optimal methods for transcription factor recovery that will allow for systematic prioritization of transcription factor candidates to aid in the design of new reprogramming protocols.

ATAC-Seq of a Single Myofiber from Mus musculus

Chromatin accessibility is a key determinant of gene expression that can be altered under different physiological and disease conditions. Skeletal muscle is made up of myofibers that are highly plastic and adaptive. Therefore, assessing the genome-wide chromatin state of myofibers under various conditions is very important to gain insight into the epigenetic state of myonuclei. The rigid nature of myofibers, as well as the low number of myonuclei that they contain, have rendered genome-wide studies with myofibers challenging. In recent years, ATAC-Seq from whole muscle and single nucleus ATAC-Seq have been performed. However, these techniques cannot distinguish between different fiber and cell types present in the muscle. In addition, due to the limited depth capacity obtained from single nucleus ATAC-Seq, an extensive comparative analysis cannot be performed. Here, we introduce a protocol where we combine the isolation of a single myofiber with OMNI ATAC-Seq. This protocol allows for genome-wide analysis of accessible chromatin regions of a selected single myofiber at a sufficient depth for comparative analysis under various physiological and disease conditions. This protocol can also allow for a specific myofiber to be selected, such as a regenerating myofiber. In the future, this protocol can help identify global changes in chromatin state under various conditions, as well as between different types of myofibers. Graphical abstract.

Transcriptional states and chromatin accessibility during bovine myoblasts proliferation and myogenic differentiation

Objectives

Although major advances have been made in bovine epigenome study, the epigenetic basis for fetal skeletal muscle development still remains poorly understood. The aim is to recapitulated the time course of fetal skeletal muscle development in vitro, and explore the dynamic changes of chromatin accessibility and gene expression during bovine myoblasts proliferation and differentiation.

Methods

PDGFR‐ cells were isolated from bovine fetal skeletal muscle, then cultured and induced myogenic differentiation in vitro in a time‐course study (P, D0, D2,and D4). The assay for transposase‐accessible chromatin sequencing (ATAC‐seq) and RNA sequencing (RNA‐seq) were performed.

Results

Among the enriched transcriptional factors with high variability, we determined the effects of MAFF, ZNF384, and KLF6 in myogenesis using RNA interference (RNAi). In addition, we identified both stage‐specific genes and chromatin accessibility regions to reveal the sequential order of gene expression, transcriptional regulatory, and signal pathways involved in bovine skeletal muscle development. Further investigation integrating chromatin accessibility and transcriptome data was conducted to explore cis‐regulatory regions in line with gene expression. Moreover, we combined bovine GWAS results of growth traits with regulatory regions defined by chromatin accessibility, providing a suggestive means for a more precise annotation of genetic variants of bovine growth traits.

Conclusion

Overall, these findings provide valuable information for understanding the stepwise regulatory mechanisms in skeletal muscle development and conducting beef cattle genetic improvement programs.

State of Knowledge on Molecular Adaptations to Exercise in Humans: Historical Perspectives and Future Directions

For centuries, regular exercise has been acknowledged as a potent stimulus to promote, maintain, and restore healthy functioning of nearly every physiological system of the human body. With advancing understanding of the complexity of human physiology, continually evolving methodological possibilities, and an increasingly dire public health situation, the study of exercise as a preventative or therapeutic treatment has never been more interdisciplinary, or more impactful. During the early stages of the NIH Common Fund Molecular Transducers of Physical Activity Consortium (MoTrPAC) Initiative, the field is well-positioned to build substantially upon the existing understanding of the mechanisms underlying benefits associated with exercise. Thus, we present a comprehensive body of the knowledge detailing the current literature basis surrounding the molecular adaptations to exercise in humans to provide a view of the state of the field at this critical juncture, as well as a resource for scientists bringing external expertise to the field of exercise physiology. In reviewing current literature related to molecular and cellular processes underlying exercise-induced benefits and adaptations, we also draw attention to existing knowledge gaps warranting continued research effort. © 2021 American Physiological Society. Compr Physiol 12:3193-3279, 2022.

A fast Myosin super enhancer dictates muscle fiber phenotype through competitive interactions with Myosin genes

The contractile properties of adult myofibers are shaped by their Myosin heavy chain isoform content. Here, we identify by snATAC-seq a 42 kb super-enhancer at the locus regrouping the fast Myosin genes. By 4C-seq we show that active fast Myosin promoters interact with this super-enhancer by DNA looping, leading to the activation of a single promoter per nucleus. A rainbow mouse transgenic model of the locus including the super-enhancer recapitulates the endogenous spatio-temporal expression of adult fast Myosin genes. In situ deletion of the super-enhancer by CRISPR/Cas9 editing demonstrates its major role in the control of associated fast Myosin genes, and deletion of two fast Myosin genes at the locus reveals an active competition of the promoters for the shared super-enhancer. Last, by disrupting the organization of fast Myosin, we uncover positional heterogeneity within limb skeletal muscles that may underlie selective muscle susceptibility to damage in certain myopathies.

The contractile properties of adult myofibers are shaped by their Myosin heavy chain isoform content. Here the authors show that a super enhancer controls the spatiotemporal expression of the genes at the fast myosin heavy chain locus by DNA looping and that this expression profile is recapitulated in a rainbow transgenic mouse model of the locus.

Application of ATAC-Seq for genome-wide analysis of the chromatin state at single myofiber resolution

Myofibers are the main components of skeletal muscle, which is the largest tissue in the body. Myofibers are highly adaptive and can be altered under different biological and disease conditions. Therefore, transcriptional and epigenetic studies on myofibers are crucial to discover how chromatin alterations occur in the skeletal muscle under different conditions. However, due to the heterogenous nature of skeletal muscle, studying myofibers in isolation proves to be a challenging task. Single-cell sequencing has permitted the study of the epigenome of isolated myonuclei. While this provides sequencing with high dimensionality, the sequencing depth is lacking, which makes comparisons between different biological conditions difficult. Here, we report the first implementation of single myofiber ATAC-Seq, which allows for the sequencing of an individual myofiber at a depth sufficient for peak calling and for comparative analysis of chromatin accessibility under various physiological and disease conditions. Application of this technique revealed significant differences in chromatin accessibility between resting and regenerating myofibers, as well as between myofibers from a mouse model of Duchenne Muscular Dystrophy (mdx) and wild-type (WT) counterparts. This technique can lead to a wide application in the identification of chromatin regulatory elements and epigenetic mechanisms in muscle fibers during development and in muscle-wasting diseases.

Six1 promotes skeletal muscle thyroid hormone response through regulation of the MCT10 transporter

BACKGROUND

The Six1 transcription factor is implicated in controlling the development of several tissue types, notably skeletal muscle. Six1 also contributes to muscle metabolism and its activity is associated with the fast-twitch, glycolytic phenotype. Six1 regulates the expression of certain genes of the fast muscle program by directly stimulating their transcription or indirectly acting through a long non-coding RNA. We hypothesized that additional mechanisms of action of Six1 might be at play.

METHODS

A combined analysis of gene expression profiling and genome-wide location analysis data was performed. Results were validated using in vivo RNA interference loss-of-function assays followed by measurement of gene expression by RT-PCR and transcriptional reporter assays.

RESULTS

The Slc16a10 gene, encoding the thyroid hormone transmembrane transporter MCT10, was identified as a gene with a transcriptional enhancer directly bound by Six1 and requiring Six1 activity for full expression in adult mouse tibialis anterior, a predominantly fast-twitch muscle. Of the various thyroid hormone transporters, MCT10 mRNA was found to be the most abundant in skeletal muscle, and to have a stronger expression in fast-twitch compared to slow-twitch muscle groups. Loss-of-function of MCT10 in the tibialis anterior recapitulated the effect of Six1 on the expression of fast-twitch muscle genes and led to lower activity of a thyroid hormone receptor-dependent reporter gene.

CONCLUSIONS

These results shed light on the molecular mechanisms controlling the tissue expression profile of MCT10 and identify modulation of the thyroid hormone signaling pathway as an additional mechanism by which Six1 influences skeletal muscle metabolism.

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.

Epigenetic rewiring of skeletal muscle enhancers after exercise training supports a role in whole-body function and human health

Objectives

Regular physical exercise improves health by reducing the risk of a plethora of chronic disorders. We hypothesized that endurance exercise training remodels the activity of gene enhancers in skeletal muscle and that this remodeling contributes to the beneficial effects of exercise on human health.

Methods and results

By studying changes in histone modifications, we mapped the genome-wide positions and activities of enhancers in skeletal muscle biopsies collected from young sedentary men before and after 6 weeks of endurance exercise. We identified extensive remodeling of enhancer activities after exercise training, with a large subset of the remodeled enhancers located in the proximity of genes transcriptionally regulated after exercise. By overlapping the position of enhancers with genetic variants, we identified an enrichment of disease-associated genetic variants within the exercise-remodeled enhancers.

Conclusion

Our data provide evidence of a functional link between epigenetic rewiring of enhancers to control their activity after exercise training and the modulation of disease risk in humans.

•

Exercise training changes in skeletal muscle gene expression is enriched for secreted factors.

•

The activity of skeletal muscle enhancers undergoes substantial remodeling after exercise training.

•

Skeletal muscle enhancer activity and gene transcription are strongly associated.

•

Exercise training-remodeled enhancer regions are enriched for GWAS SNPs associated with human traits and diseases.

Myod1 and GR coordinate myofiber-specific transcriptional enhancers

Skeletal muscle is a dynamic tissue the size of which can be remodeled through the concerted actions of various cues. Here, we investigated the skeletal muscle transcriptional program and identified key tissue-specific regulatory genetic elements. Our results show that Myod1 is bound to numerous skeletal muscle enhancers in collaboration with the glucocorticoid receptor (GR) to control gene expression. Remarkably, transcriptional activation controlled by these factors occurs through direct contacts with the promoter region of target genes, via the CpG-bound transcription factor Nrf1, and the formation of Ctcf-anchored chromatin loops, in a myofiber-specific manner. Moreover, we demonstrate that GR negatively controls muscle mass and strength in mice by down-regulating anabolic pathways. Taken together, our data establish Myod1, GR and Nrf1 as key players of muscle-specific enhancer-promoter communication that orchestrate myofiber size regulation.

Single-nucleus RNA-seq and FISH identify coordinated transcriptional activity in mammalian myofibers

Skeletal muscle fibers are large syncytia but it is currently unknown whether gene expression is coordinately regulated in their numerous nuclei. Here we show by snRNA-seq and snATAC-seq that slow, fast, myotendinous and neuromuscular junction myonuclei each have different transcriptional programs, associated with distinct chromatin states and combinations of transcription factors. In adult mice, identified myofiber types predominantly express either a slow or one of the three fast isoforms of Myosin heavy chain (MYH) proteins, while a small number of hybrid fibers can express more than one MYH. By snRNA-seq and FISH, we show that the majority of myonuclei within a myofiber are synchronized, coordinately expressing only one fast Myh isoform with a preferential panel of muscle-specific genes. Importantly, this coordination of expression occurs early during post-natal development and depends on innervation. These findings highlight a previously undefined mechanism of coordination of gene expression in a syncytium.

Whether skeletal muscle fibre gene expression is coordinated as a whole in different nuclei in the fibre is unclear. Here, the authors use single nucleus RNAseq and ATACseq to show the transcriptome heterogeneity of muscle nuclei in the adult mouse fibre, with correlations between the two datasets.

Changes in Northern Elephant Seal Skeletal Muscle Following Thirty Days of Fasting and Reduced Activity

Northern elephant seals (NES, Mirounga angustirostris) undergo an annual molt during which they spend ∼40 days fasting on land with reduced activity and lose approximately one-quarter of their body mass. Reduced activity and muscle load in stereotypic terrestrial mammalian models results in decreased muscle mass and capacity for force production and aerobic metabolism. However, the majority of lost mass in fasting female NES is from fat while muscle mass is largely preserved. Although muscle mass is preserved, potential changes to the metabolic and contractile capacity are unknown. To assess potential changes in NES skeletal muscle during molt, we collected muscle biopsies from 6 adult female NES before the molt and after ∼30 days at the end of the molt. Skeletal muscle was assessed for respiratory capacity using high resolution respirometry, and RNA was extracted to assess changes in gene expression. Despite a month of reduced activity, fasting, and weight loss, skeletal muscle respiratory capacity was preserved with no change in OXPHOS respiratory capacity. Molt was associated with 162 upregulated genes including those favoring lipid metabolism. We identified 172 downregulated genes including those coding for ribosomal proteins and genes associated with skeletal muscle force transduction and glucose metabolism. Following ∼30 days of molt, NES skeletal muscle metabolic capacity is preserved although mechanotransduction may be compromised. In the absence of exercise stimulus, fasting-induced shifts in muscle metabolism may stimulate pathways associated with preserving the mass and metabolic capacity of slow oxidative muscle.

Myogenesis control by SIX transcriptional complexes

SIX homeoproteins were first described in Drosophila, where they participate in the Pax-Six-Eya-Dach (PSED) network with eyeless, eyes absent and dachsund to drive synergistically eye development through genetic and biochemical interactions. The role of the PSED network and SIX proteins in muscle formation in vertebrates was subsequently identified. Evolutionary conserved interactions with EYA and DACH proteins underlie the activity of SIX transcriptional complexes (STC) both during embryogenesis and in adult myofibers. Six genes are expressed throughout muscle development, in embryonic and adult proliferating myogenic stem cells and in fetal and adult post-mitotic myofibers, where SIX proteins regulate the expression of various categories of genes. In vivo, SIX proteins control many steps of muscle development, acting through feedforward mechanisms: in the embryo for myogenic fate acquisition through the direct control of Myogenic Regulatory Factors; in adult myofibers for their contraction/relaxation and fatigability properties through the control of genes involved in metabolism, sarcomeric organization and calcium homeostasis. Furthermore, during development and in the adult, SIX homeoproteins participate in the genesis and the maintenance of myofibers diversity.