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

  1. Hirotaka Araki
  2. Shinjiro Hino  Is a corresponding author
  3. Kotaro Anan
  4. Kanji Kuribayashi
  5. Kan Etoh
  6. Daiki Seko
  7. Ryuta Takase
  8. Kensaku Kohrogi
  9. Yuko Hino
  10. Yusuke Ono
  11. Eiichi Araki
  12. Mitsuyoshi Nakao  Is a corresponding author
  1. Department of Medical Cell Biology, Institute of Molecular Embryology and Genetics, Kumamoto University, Japan
  2. Department of Metabolic Medicine, Faculty of Life Sciences, Kumamoto University, Japan
  3. Department of Muscle Development and Regeneration, Institute of Molecular Embryology and Genetics, Kumamoto University, Japan
  4. Department of Molecular Bone Biology, Graduate School of Biomedical Sciences, Nagasaki University, Japan

Abstract

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.

Editor's evaluation

This paper investigates the role of lysine-specific histone demethylase 1 (LSD1) in skeletal muscle responses to stresses. The findings link LSD1 to responses of skeletal muscle to glucorticoids and exercise in a fiber-type-dependent manner. The paper should be of broad interest to those studying muscle biology and physiology.

https://doi.org/10.7554/eLife.84618.sa0

Introduction

Skeletal muscle is the largest organ in humans, and the loss of muscle mass and physical activities are linked to the risks of metabolic disorders and aging-associated diseases. Skeletal muscle is composed of distinct types of myofibers. The slow-twitch fibers (type I) have high respiratory potential, and thus, a high endurance capacity, while the fast-twitch fibers (type II) have high glycolytic capacity and are specialized for dynamic movements (Schiaffino and Reggiani, 2011). In mice, type II fibers can be subdivided into three groups, type IIA, type IIX, and type IIB, with the type IIA relatively more oxidative than the types IIX and IIB. These myofibers are subject to gain or loss, with different sensitivities in response to numerous environmental cues, such as nutritional conditions, hormones, drugs, and exercise (Braun and Marks, 2015; Cartee et al., 2016; Hawley et al., 2014). For example, the fast fibers are preferentially affected in glucocorticoid (GC)-induced muscle atrophy, which is associated with certain environmental and pathological conditions, such as starvation, GC drug therapy, sepsis, cachexia, acute diabetes, and hypoxia (Braun and Marks, 2015; Schakman et al., 2013; Kuo et al., 2013; Przygodda et al., 2017). Although GCs induce atrophy by enhancing protein catabolism via the autophagy and ubiquitin-proteasome system (Sartori et al., 2021), the molecular basis for the fiber type preference is still unclear. Another environmental condition that can remarkably alter muscle phenotypes is exercise. By increasing muscle mass, resistance, and endurance, exercise confers health benefits in the young, old, and sick individuals (Cartee et al., 2016). In particular, endurance exercise promotes the formation of oxidative fibers, thus increasing the oxidative capacity in the muscle, which can effectively improve systemic metabolism in diabetic patients and prevent disuse-induced atrophy (Yan et al., 2011; Qaisar et al., 2016).

Several chromatin factors promote phenotypic plasticity in the skeletal muscle. SMYD3, a histone methyltransferase, promotes GC-induced atrophy in mice by enhancing the expression of atrophy-associated genes, such as myostatin and Met (Proserpio et al., 2013). In addition, SETD7, another histone methyltransferase, and histone deacetylases-4 and -5 have essential roles in cachexia and denervation-induced atrophy, respectively (Nayak et al., 2019; Luo et al., 2019; Moresi et al., 2010). Moreover, DNA methyltransferase 3A is essential for muscle adaptation to endurance exercise by suppressing the oxidative stress response in the muscle (Damal Villivalam et al., 2021). These previous findings in literature have established that epigenetic factors are the integral components of muscle adaptation. However, factors that define the fiber type-specific responses to environmental factors have not been identified.

Lysine-specific demethylase-1 (LSD1 or KDM1A) is a flavin-dependent amine oxidase that in general removes mono- and di-methyl groups of histone H3 at lysine 4 (H3K4) and other protein substrates (Hino et al., 2016; Shi et al., 2004; Nakao et al., 2019), suggesting its contribution to specific biological processes. Despite ubiquitous expression in various tissue types, LSD1 plays essential roles in myogenic differentiation and muscle regeneration (Choi et al., 2010; Scionti et al., 2017; Tosic et al., 2018). Notably, we recently demonstrated that LSD1 promotes fast-like myogenic differentiation by epigenetically repressing slow fiber genes in C2C12 mouse myoblasts (Anan et al., 2018). Thus, the available evidence raises the possibility that LSD1 influences the adaptive phenotypes in the muscle under environmental cues. To verify this, we generated muscle-specific/inducible LSD1-knockout mice and examined how LSD1 loss affects muscle plasticity under environmental stresses. We found that LSD1 depletion exacerbated GC-induced fast fiber atrophy while augmenting endurance exercise-induced slow fiber hypertrophy. In these adaptive processes, LSD1 regulated distinct sets of target genes, together with unique transcription factors (TFs). Collectively, our current study reveals the mechanism through which LSD1 builds the epigenetic barriers that restrict adaptive responses in the muscles.

Results

Muscle-specific LSD1 depletion exacerbates dexamethasone-induced muscle atrophy

LSD1 is also known as KDM1A. For clarity, we use the word LSD1 for the protein and Lsd1 for gene and mRNA. To investigate the role of LSD1 in muscle plasticity in vivo during postnatal life, we generated tamoxifen-inducible-, skeletal muscle-specific LSD1-knockout (LSD1-mKO) mice by crossing LSD1-floxed mice (LSD1flox/flox) and transgenic mice harboring ACTA1-creERT2 (McCarthy et al., 2012; Figure 1—figure supplement 1A). This strategy facilitated LSD1 depletion in differentiating and differentiated muscle cells. We confirmed that LSD1 was depleted after tamoxifen administration in both fast-twitch fiber-dominant (tibialis anterior [TA], gastrocnemius [Gas], and extensor digitorum longus) and slow-twitch fiber-dominant (soleus [Sol]) tissues (Figure 1—figure supplement 1B, C). The magnitude of LSD1 depletion was comparable among these tissues. Marginal amounts of LSD1 were present after Cre induction, probably due to the expression in hematopoietic and vascular cells that are resident in muscle tissues (Kerenyi et al., 2013; Yuan et al., 2022). In contrast, LSD1 expression in the heart, liver, and adipose tissues remained unaffected. Next, we depleted LSD1 in growing mice by administrating tamoxifen to pre-weaning mice (11 days of age) (Figure 1—figure supplement 2A). The expression of muscle fiber genes was largely unaffected, suggesting that LSD1 depletion in our system did not affect postnatal muscle growth (Figure 1—figure supplement 2B).

Dexamethasone (Dex), a synthetic GC, affects both the muscle mass and quality by inducing protein catabolism in fast fiber-dominant muscles (Quattrocelli et al., 2019; Figure 1—figure supplement 3A–C). To determine the involvement of LSD1 in muscle plasticity under environmental stress, we administered Dex in wild type (WT) and LSD1-mKO mice after knockout (KO) induction (Figure 1A). We found that the atrophy of TA and Gas muscles was enhanced in LSD1-mKO mice (Figure 1B), which was also evident from their reduced body weight (Figure 1—figure supplement 4A). Consistently, major hallmarks of muscle atrophy, including enhanced processing of autophagy protein LC3, increased expression of p62, attenuated phosphorylation of Akt, were pronounced in the LSD1-mKO muscle (Figure 1—figure supplement 5A–C). Another atrophy marker, 4E-BP phosphorylation, remained unaffected (Figure 1—figure supplement 5D). Of note, LSD1 depletion alone did not affect the muscle size and body weight (Figure 1B, and Figure 1—figure supplement 3D, E), indicating that the effects of Dex were enhanced by the absence of LSD1. In contrast, the Sol muscle, which was only modestly sensitive to Dex, and non-muscle tissues were unaffected by LSD1 depletion (Figure 1B and Figure 1—figure supplement 4B). Next, we evaluated the muscle performance of LSD1-mKO mice after Dex administration. LSD1 depletion led to weakened muscle strength as revealed by a whole limb grip test (Figure 1C), while LSD1 depletion alone was not sufficient to induce this phenotype (Figure 1—figure supplement 3F). In contrast, the running endurance was unaffected as assessed by a treadmill run (Figure 1—figure supplement 4C). These results suggest that LSD1 protects fast-dominant muscles from Dex-induced atrophy, thereby maintaining muscle strength.

Figure 1 with 7 supplements see all
Muscle-specific lysine-specific demethylase-1 (LSD1)-knockout exacerbates dexamethasone (Dex)-induced fast muscle atrophy.

(A) Experimental design: LSD1 deletion was induced by tamoxifen injection for 5 days (days 0–4), and Dex was administered for 7 days (days 3–9) at 1 mg/kg. (B) Muscle weight normalized by the body weight. Tibialis anterior (TA), gastrocnemius (Gas), and soleus (Sol) muscles were obtained from Dex-treated wild type (WT) (n=40) and LSD1-mKO (n=44) mice. In the control set, WT (n=22) and mKO (n=22) mice that were not subjected to Dex treatment were sacrificed 1 month after the first tamoxifen administration. Box-plot elements (center line, median; box limits, upper and lower quartiles; whiskers, min to max). (C) All-limb grip strength of Dex-treated WT (n=10) and LSD1-mKO (n=12) mice. (D) Immunofluorescence staining of individual fiber types in Gas muscles of WT and LSD1-mKO mice after Dex treatment. Left: cross section of the whole tissue (scale bars, 300 µm). Right: magnified images of dotted squares in the images on the left (scale bars, 100 µm). Type I fiber (blue), type IIA fiber (green), and laminin (red) were stained. The fibers that were negative for types I and IIA were defined as type IIB+IIX. (E) Occupancy of each fiber type based on the cross-sectional area (WT, n=6; mKO, n=6). (F) Frequency of each muscle fiber type based on the fiber number (WT, n=6, mKO, n=6). Values are mean ± standard deviation (SD). *p<0.05, **p<0.01.

Oxidative fiber ratio is increased in the fast-dominant muscle in LSD1-mKO mice after Dex administration

Because the effects of LSD1 depletion were specific to fast-dominant muscles, we examined if the fiber-type composition was affected in the KO mice under Dex treatment. Immunostaining of the Gas muscle revealed that the area and number of oxidative (type I and IIA) fibers were higher in LSD1-mKO mice than in WT (Figure 1D-F, and Figure 1—figure supplement 6A). In contrast, the area and number of glycolytic (type IIB+IIX) fibers were decreased in LSD1-mKO (Figure 1D–F, and Figure 1—figure supplement 6B). Such genotype-dependent differences were not observed in the absence of Dex (Figure 1—figure supplement 6C,D and Figure 1—figure supplement 7A,B). In both genotypes, we did not observe centralized nuclei in the muscle fiber, which is characteristic of newly regenerating myofibers (Figure 1—figure supplement 7C). This suggests that the increase in the number of slow-twitch fibers in the KO was not due to enhanced de novo fiber formation. Additionally, we analyzed the size distribution of myofibers according to cross-sectional area and found that large type I and type IIA fibers preferentially increased in the LSD1-mKO muscle, whereas large type IIB+IIX fibers decreased (Figure 1—figure supplements 6B and 7D and E). Taken together, these results suggest that oxidative fibers are preferentially formed and/or maintained in LSD1-depleted muscles under Dex treatment.

Atrophy- and slow fiber-associated genes are upregulated in fast-dominant muscles in Dex-treated LSD1-mKO mice

To gain functional insight into the role of LSD1 in muscle plasticity, we performed RNA-sequencing (RNA-seq) analysis of the Gas muscle after Dex treatment. We found 150 upregulated and 117 downregulated genes, respectively, in LSD1-mKO mice, compared with those in the WT (Figure 2A). Consistent with the phenotype, a group of genes associated with muscle atrophy were upregulated in LSD1-mKO mice (Figure 2B and Supplementary file 1). By using quantitative reverse transcription-polymerase chain reaction (qRT-PCR) analysis, we confirmed that the key atrophy mediators, including Gadd45a, Gabarapl1, and Hspb7, were significantly upregulated in the KO mice (Figure 2C; Ebert et al., 2012; Schaaf et al., 2016; Tobin et al., 2016; Bullard et al., 2016). In parallel, hypertrophy genes were downregulated in the LSD1-mKO Gas muscle (Figure 2B and Figure 2—figure supplement 1A). Furthermore, in the LSD1-mKO mice, the slow fiber genes, including Sln (encoding sarcolipin), a functional hallmark of oxidative fibers (Maurya et al., 2018), were significantly upregulated (Figure 2B and D). In contrast, fast fiber genes were downregulated in these mice. We also observed a similar expression pattern for fiber type-associated genes in the TA muscle (Figure 2—figure supplement 1B). These expression changes were not observed in the muscles that were not treated with Dex (Figure 2—figure supplement 1C and D). Intriguingly, RNA-seq analysis of the Sol muscle showed that the differentially regulated genes (WT vs. KO) were incompatible with those in the Gas muscle (Figure 2E and Figure 2—figure supplement 2A and B). In particular, a number of atrophy-associated genes were downregulated in the LSD1-mKO Sol, which was in clear contrast to the expression changes observed in the Gas (Figure 2C and Figure 2—figure supplement 2C and D). Taken together, these results indicate that LSD1 buffers Dex effects by repressing the atrophy- and slow fiber-associated transcriptome in fast-dominant muscle.

Figure 2 with 3 supplements see all
Transcriptome analysis reveals the upregulation of muscle atrophy and slow fiber genes in fast-dominant muscles in dexamethasone (Dex)-treated lysine-specific demethylase-1 (LSD1)-mKO mice.

(A) Comparison of the transcriptome data from LSD1-mKO and wild type (WT) gastrocnemius (Gas) muscles. The red and blue dots indicate significantly up- and downregulated genes, respectively, in LSD1-mKO (FDR >0.05). (B) Top 50 differentially expressed genes ranked by the FDR Z-score. In the upregulated gene list, muscle atrophy- (red) and slow fiber-associated (blue) genes are highlighted. In the downregulated list, muscle hypertrophy- (red) and fast fiber-associated genes (blue) are highlighted. The complete gene list is provided in Supplementary file 1. (C and D) Quantitative reverse transcription-polymerase chain reaction (qRT-PCR) analysis of atrophy- (C) and fiber type-associated genes (D) in the Dex-treated LSD1-mKO Gas muscles (WT, n=8; mKO, n=11). Expression levels were normalized to those of the 36B4 gene and shown as fold differences against the WT. Full descriptions of gene symbols are provided in Supplementary file 3. (E) Comparison of the transcriptome profiles of Gas and soleus (Sol) muscles from Dex-treated LSD1-mKO mice. The dark red and blue bars indicate the fraction of genes that are commonly upregulated and downregulated, respectively, in the mKO Gas and Sol muscles. The numbers of genes that were up- or downregulated in mKO Gas are also indicated. Values are mean ± SD. *p<0.05, **p<0.01.

We next tested whether LSD1 has protective function in other atrophy models. For this, we administrated a pro-diabetic agent streptozotocin (STZ), which disrupts the function of pancreatic islet leading to fast-fiber atrophy in mice (Figure 2—figure supplement 3A; O’Neill et al., 2019). LSD1-mKO did not affect the muscle weight in STZ-treated mice (Figure 2—figure supplement 3B and C). Consistently, there were no major differences in the expression of atrophy and myofiber genes in STZ-treated WT and LSD1-mKO mice (Figure 2—figure supplement 3D–E). These results suggest that the LSD1 function depends on the source of atrophy-inducing stress.

LSD1 regulates atrophy-associated genes by promoting the nuclear function of Foxk1

We previously reported that LSD1 inhibition led to an accumulation of tri-methylated histone H3 lysine 4 (H3K4me3) at the LSD1-bound promoters in differentiating myoblasts (Anan et al., 2018). To understand the mechanism through which LSD1 regulates the atrophy program, we evaluated the effects of LSD1 depletion on the H3K4 methylation status. We found that the H3K4me3 levels at the promoter regions of LSD1-mKO-suppressed genes were higher in the KO Gas muscle compared to those in the WT muscle (Figure 3A). In addition, we analyzed a previously published chromatin immunoprecipitation (ChIP)-seq dataset in mouse myoblasts (Tosic et al., 2018), and found that LSD1 bound to the promotor regions of the atrophy genes that were upregulated in LSD1-mKO Gas (Figure 3—figure supplement 1A).

Figure 3 with 2 supplements see all
Lysine-specific demethylase-1 (LSD1) regulates the expression of atrophy-associated genes through the modulation of Foxk1.

(A) Chromatin immunoprecipitation (ChIP)-qPCR analyses of H3K4me3 at Sln, Gadd45a, and Ddit4 gene promoters in LSD1-mKO gastrocnemius (Gas) muscles after dexamethasone (Dex) treatment. The data are from three independent experiments. The enrichment levels of H3K4me3 were normalized to the input DNA and presented as fold differences in the LSD1-mKO against that of wild type (WT). (B) Motif analyses of the promoter regions of genes that were upregulated in Dex-treated LSD1-mKO Gas muscles. Analysis of the atrophy-associated genes identified Foxk1 and PITX1 motifs, while significant enrichment of NFATc1 motif was found when all the upregulated genes were scanned. (C) Expression of genes encoding the transcription factors (TFs) identified by the motif analyses. Quantitative reverse transcription-polymerase chain reaction (qRT-PCR) values (n=4) are shown as the fold differences, compared with those of the Gas muscles. (D) Effects of treatment with Dex and LSD1 inhibitor (T-3775440) on the expression of Foxk1 in C2C12 myotubes. Foxk1 and 5,6-diamidino-2-phenylindole (DAPI) co-staining revealed Foxk1-positive nuclei, which are indicated by arrowheads (scale bars, 20 μm). The graph shows the percentage of Foxk1-positive nuclei (n=3). (E) Co-immunoprecipitation of LSD1 and Foxk1. C2C12 myotubes were treated with insulin for 6 hr before harvest to enhance the nuclear retention of Foxk1. Input lane contains 10% amount of the whole-cell extract. (F) Effects of Foxk1-KO on the expression of atrophy- and slow fiber-associated genes in C2C12 myotubes (n=6). qRT-PCR values are shown as the fold differences, compared with those in the control-transfected cells (n=6). (G) ChIP-qPCR analyses of LSD1 and Foxk1 at Ddit4, Gadd45a, Padi2 (negative control) gene promoters in control and Foxk1-KO C2C12 myotubes. The enrichment levels of normal rabbit IgG, LSD1, and Foxk1 were normalized to the input DNA. Values are mean ± SD. *p<0.05, **p<0.01.

To identify the transcriptional regulators that mediate the regulation of atrophy-associated genes by LSD1, we performed motif analyses on the promotor regions of upregulated genes in LSD1-mKO Gas and identified the binding motifs for Foxk1, PITX1, and NFATc1 (Figure 3B). In C2C12 cells, Foxk1 represses atrophy genes under nutrient-rich conditions and becomes sequestered from the nucleus upon starvation (Bowman et al., 2014). In addition, the repressive function of Foxk1 is achieved through its interaction with Sin3A (Bowman et al., 2014), which also forms a complex with LSD1 (Yang et al., 2018). Of note, the expression of Foxk1 in the Gas and TA muscles was higher than that in the Sol muscle (Figure 3C), raising the possibility that Foxk1 mediates the muscle site-dependent function of LSD1. Immunofluorescence and western blot analyses in C2C12 myotubes showed that Dex exposure reduced the nuclear retention of Foxk1, which was further promoted by the addition of T-3775440, an LSD1 inhibitor (Figure 3D and Figure 3—figure supplement 1B). We also performed a co-immunoprecipitation (co-IP) experiment and found a physical interaction between Foxk1 and LSD1 in C2C12 myotubes (Figure 3E). In addition, both LSD1 and Foxk1 interacted with the common partner Sin3A (Figure 3—figure supplement 1C). Furthermore, we created Foxk1-KO C2C12 cells using a CRISPR Cas9 system and demonstrated that the key LSD1-targeted atrophy genes in these cells were upregulated after differentiation (Figure 3F and Figure 3—figure supplement 2A). Foxk1-KO C2C12 myotubes also showed upregulated expression of type I and IIA fiber genes, which was consistent with the results observed in the LSD1-mKO Gas muscle (Figure 3F). In these cells, the LSD1 inhibition did not further affect the expression of atrophy and myofiber genes (Figure 3—figure supplement 2B). ChIP experiment revealed that the occupancy of LSD1 at atrophy gene loci was diminished in the absence of Foxk1 (Figure 3G and Figure 3—figure supplement 2C). These data indicate that LSD1 largely depends on Foxk1 to regulate the atrophy genes. Moreover, the upregulation of atrophy genes was also observed in the C2C12 myotubes expressing Foxk1 shRNA (Figure 3—figure supplement 2D and E). Altogether, these results suggests that LSD1 and Foxk1 cooperatively regulate the atrophy transcriptome in the presence of Dex.

Muscle-specific LSD1 depletion promotes the hypertrophy of Sol muscle and the muscle endurance after voluntary exercise

The differential effects of LSD1 depletion on the transcriptome in fast- and slow-dominant muscles under Dex administration motivated us to explore the physiological role of LSD1 in the slow-dominant muscles. For this purpose, we subjected the mice to voluntary wheel running (VWR). VWR augments muscle endurance capacity and induces transcriptomic changes in muscles in a site-dependent manner (Legerlotz et al., 2008; Kim et al., 2020; Schmitt et al., 2020). In particular, VWR has been shown to induce the hypertrophy of soleus in mice and rats (Legerlotz et al., 2008; Momken et al., 2005; Konhilas et al., 2005). To allow VWR, we kept LSD1-mKO and WT mice in a cage attached to a running wheel for over 1 month (Figure 4A). We confirmed that tamoxifen administration during the first 5 days of VWR sufficiently maintained an LSD1-deficient state until sacrifice at day 34 (Figure 4A and Figure 4—figure supplement 1A). The mice trained by VWR showed a prominent weight gain in the Sol muscle, which was further augmented by LSD1 depletion (Figure 4B and Figure 4—figure supplement 1B). Histologically, the area of type I muscle fibers significantly increased in the LSD1-KO Sol (Figure 4C and D, and Figure 4—figure supplement 1C–E). In a treadmill endurance test, LSD1-mKO mice ran a longer distance and for a longer duration than the WT mice (Figure 4E), while the grip strength was comparable between the two genotypes (Figure 4—figure supplement 1F). Consistent with the results of the augmented endurance capacity, the Sol muscle in the KO mice showed enhanced succinate dehydrogenase (SDH) staining, indicating that the OXPHOS capacity of Sol was enhanced by the loss of LSD1 (Figure 4F and Figure 4—figure supplement 1G). The total running distance during the VWR was not influenced by genotype (Figure 4—figure supplement 1H), eliminating the possibility that the observed phenotypes were related to exercise intensity. Collectively, these data show that LSD1 depletion in the muscle promoted adaptation to prolonged physical activities by enhancing hypertrophy of the slow muscle.

Figure 4 with 1 supplement see all
Lysine-specific demethylase-1 (LSD1)-mKO augments the effects of voluntary wheel running (VWR) by enhancing muscle endurance and increasing the number of oxidative fibers.

(A) Experimental design: after familiarization with a cage equipped with a running wheel, LSD1 deletion was induced by tamoxifen injection (days 0–4). Mice were subjected to VWR in the same cage for 38 days (days –3–34). (B) Muscle weight normalized to body weight in sedentary wild type (WT) (n=8) and LSD1-mKO mice (n=7) as well as in VWR-trained WT (n=23) and LSD1-mKO mice (n=22). (C) Immunofluorescence staining of individual fiber types in the whole soleus (Sol) muscles from VWR-trained WT and LSD1-mKO mice. Type I fiber (blue), type IIA fiber (green), and laminin (red) were stained. Magnified images of the tissue sections (scale bars, 100 µm). Representative images (WT and mKO, n=8) are shown. (D) Occupancy of type I fiber based on the cross-sectional area (WT and mKO, n=8). (E) Treadmill performance of WT (n=9) and LSD1-mKO (n=8) mice after VWR. Running distance and time until exhaustion are shown. (F) Succinate dehydrogenase (SDH) staining of the Sol muscles from trained WT and LSD1-mKO mice. Scale bar, 100 μm.

LSD1 depletion upregulates the expression of ERRγ and mitochondrial metabolism genes after voluntary exercise

To identify the genes regulated by LSD1 under VWR, we subjected the Sol muscle from trained WT and LSD1-mKO mice to RNA-seq analysis (Figure 5A). We performed a non-hierarchical clustering and identified a cluster of genes (cluster 6) that were upregulated in the LSD1-mKO. Consistent with the altered oxidative capacity, this cluster was associated with gene ontology terms, such as ‘mitochondrion’, ‘mitochondrion organization’, and ‘oxidative phosphorylation’ (Figure 5B). Of note, we found a significant upregulation of Esrrg gene, which encodes ERRγ, a TF that promotes oxidative metabolism in the muscle (Figure 5C and D, and Supplementary file 2; Misra et al., 2017; Fan et al., 2018; Gan et al., 2013; Narkar et al., 2011). In sedentary mice, the expression of Esrrg was not affected by LSD1 depletion (Figure 5—figure supplement 1A), suggesting a combinatorial effect of VWR and LSD1-mKO. The gene encoding PGC-1α, a pivotal transcriptional factor that promotes oxidative metabolism and endurance capacity (Handschin et al., 2007), was not affected in the LSD1-mKO mice (Figure 5D). Interestingly, the protein level of ERRγ was upregulated by VWR and further by LSD1 depletion, suggesting that the exercise-activated ERRγ function was potentiated in the LSD1-mKO muscle (Figure 5E). As we performed a motif analysis of upstream regions of the cluster 6 genes, we found a significant enrichment of the ERRγ binding motif (Figure 5F). qRT-PCR results confirmed that the known ERRγ targets as well as the genes harboring the ERRγ motif were upregulated in the LSD1-mKO Sol (Figure 5G; Cho et al., 2013; Liang et al., 2016). The expression of ERRγ and its target genes was not influenced by LSD1 depletion in the Gas and TA samples (Figure 5—figure supplement 1B).

Figure 5 with 1 supplement see all
Lysine-specific demethylase-1 (LSD1)-mKO augments the effects of voluntary wheel running (VWR) by increasing the expression of ERRγ and its target genes.

(A) RNA-seq analysis of wild type (WT) and LSD1-mKO soleus (Sol) muscles after voluntary wheel running (VWR) training. K-means clustering identified a group of genes that were upregulated in mKO (cluster 6). (B) RNA-seq analysis of WT and LSD1-mKO Sol muscles after VWR training. The genes upregulated in mKO (cluster 6) were subjected to a gene ontology analysis. (C) The expression of Esrrg in WT (n=3) and LSD1-mKO (n=3) Sol muscles after VWR training from RNA-seq data. (D) The expression of Esrrg and functionally related genes in WT (n=6) and LSD1-mKO (n=6) Sol muscles after VWR training. (E) The expression of ERRγ protein in Sol muscles from sedentary WT, trained WT, and trained LSD1-mKO mice. (F) Motif analyses of promoter regions of the cluster 6 genes. Values are mean ± SD. *p<0.05, **p<0.01. (G) Expression of ERRγ target genes in WT (n=6) and LSD1-mKO (n=6) Sol after VWR training. Quantitative reverse transcription-polymerase chain reaction (qRT-PCR) values are shown as the fold differences against the WT. Full descriptions of gene symbols are provided in the Supplementary file 3. Values are mean ± SD. *p<0.05, **p<0.01.

To further test whether Esrrg expression is downregulated by LSD1, we made use of H9c2 rat cardiomyocytes, in which ERRγ is functionally involved in the differentiation (Sakamoto et al., 2022). Inhibition of LSD1 in H9c2 cells led to significant upregulation of Esrrg and its downstream Perm1 (Figure 5—figure supplement 1C), whereas these genes were downregulated in LSD1-overexpressing H9c2 cells (Figure 5—figure supplement 1D and E). Collectively, these data suggest that LSD1 depletion enhanced oxidative metabolism and muscle endurance by promoting ERRγ activity in VWR-trained mice.

Finally, we tested whether aging affects the expression of LSD1 in the muscle. In mice, muscle weight declined with age, which was accompanied by the increased expression of atrophy genes (Figure 6A and B). In parallel, LSD1 expression was downregulated in aged mice both at RNA and protein levels (Figure 6B and C). Next, we analyzed the expression of LSD1 in individual muscle cells using single-cell RNA-seq data from the Tabula Muris Senis project (Tabula Muris Consortium, 2020). We found that Lsd1 levels decline with age in both Myod1- and Pax7-positive myogenic cells (Figure 6D). In human muscle, LSD1 expression showed a significant negative correlation with age (Figure 6E). These data, together with our observations in LSD1-mKO mice, imply that aging-associated loss of LSD1 confers the sensitivity against environmental stresses leading to qualitative and quantitative changes in the muscle.

Lysine-specific demethylase-1 (LSD1) expression is decreased in the aged muscles in mice and humans.

(A) Muscle weight in young (2.5 months of age, n=4) and old (28 months of age, n=3) male C57BL/6J mice. (B) The expression of fiber type-specific genes and atrophy-associated genes in the gastrocnemius (Gas) muscles in young and old mice. Quantitative reverse transcription-polymerase chain reaction (qRT-PCR) values are shown as fold differences against the young mice. (C) Expression of LSD1 protein in the EDL muscles from the young and old mice. Band densities were quantified by densitometry and normalized to histone H3. Values are shown as fold differences against the young. (D) Decreased expression of Lsd1 in aged muscles. Using single-cell RNA-seq data (flow cytometry) published by the Tabula Muris Senis (Tabula Muris Consortium, 2020), Lsd1 expression was analyzed in Myod1-positive and Pax7-positive muscle cells from young to old mice (3–24 months of age). (E) The expression of LSD1 in skeletal muscles of 55- to 80-year-old men and women. RNA-seq data by Tumasian et al., 2021, was analyzed. Values are mean ± SD. R2: coefficient of determination. *p<0.05, **p<0.01.

Discussion

Skeletal muscle exhibits notable plasticity in response to environmental cues, which can be beneficial or deleterious to health. In the current study, we demonstrated that the loss of LSD1 remarkably sensitized the skeletal muscle to Dex-induced atrophy and VWR-induced hypertrophy. Under these stressed conditions, LSD1 suppressed the catabolic pathways as well as the formation of type I fibers. Mechanistically, LSD1 cooperated with Foxk1 to repress protein catabolism genes, while repressing ERRγ and its downstream oxidative genes, depending on the stimulus and muscle location (Figure 7). Environmental adaptation involves transcriptional and epigenetic remodeling in response to environmental cues, which confer a phenotypic plasticity in cells (Feil and Fraga, 2012; Peters et al., 2021; Simpson et al., 2011; Gilbert and Epel, 2009; Cavalli and Heard, 2019). Epigenetic mechanisms do not only promote phenotypic alterations, but also serve as barriers that safeguard the identity and functions of cells under environmental fluctuations (Hörmanseder, 2021; Flavahan et al., 2017). While epigenetic factors that promote muscle adaptation have been reported (Proserpio et al., 2013; Nayak et al., 2019; Luo et al., 2019; Moresi et al., 2010; Damal Villivalam et al., 2021), factors that suppress environmental responses in the muscle were unknown. The current study reveals that LSD1 serves as an ‘epigenetic barrier’ that defines stress sensitivities in the muscle.

Lysine-specific demethylase-1 (LSD1) serves as an ‘epigenetic barrier’ that defines stress sensitivities in the skeletal muscle.

LSD1 attenuates glucocorticoid (GC)-induced atrophy in the fast fiber-dominant muscles, in collaboration with Foxk1, an anti-autophagic transcription factor. On the other hand, LSD1 attenuates endurance exercise-induced hypertrophy in the slow fiber-dominant muscles, by inhibiting ERRγ, a transcription factor that promotes oxidative metabolism genes. The loss of LSD1 remarkably sensitized the muscles to GC and endurance exercise.

In the current study, we generated a muscle-specific LSD1-KO mice by crossing LSD1-floxed mice with ACTA1-CreERT2 mice. This enabled us to inactivate LSD1 in differentiating and differentiated muscle cells. We did not find major phenotypic changes in these mice without environmental insults. A previous report by Tosic et al., 2018 demonstrated that the muscle-specific loss of LSD1 led to a decline in the regenerative capacity and the acquisition of brown adipocyte-like characteristics in the muscle cells. Tosic et al. made use of Pax7-cre to engineer Lsd1 gene specifically in muscle satellite cells. Thus, it is likely that the difference in the engineering method accounts for the phenotypic divergence, and that LSD1 has a maturation stage-dependent function in the muscle cell lineage.

A cascade of signaling pathways and transcriptional regulators operate under GC exposure to induce atrophic gene expression (Schakman et al., 2013; Kuo et al., 2013). We demonstrated that LSD1 inhibition combined with Dex treatment led to impaired nuclear retention of Foxk1, resulting in the de-repression of atrophy genes. Previous reports showed that, under nutrient-rich conditions, the Akt-mTORC1 pathway promotes the nuclear localization of Foxk1 to repress autophagy and other catabolic genes (Bowman et al., 2014; Sakaguchi et al., 2019). Because GCs exert their pro-atrophic effects partly through the inhibition of the Akt-mTORC1 axis (Schakman et al., 2013), this pathway may be involved in the nuclear retention of Foxk1. Thus, LSD1 counteracts GC signaling to facilitate the role of Foxk1, enabling optimal muscle adaptation to pro-atrophic stress. In addition, upon the sequestration of Foxk1 from the nucleus under starvation, Foxk1 target genes are replaced by Foxo3, which is known as a downstream effector of GR signaling (Bowman et al., 2014; Cho et al., 2010; Sandri et al., 2004). Thus, LSD1 depletion promotes the accessibility of Foxo3 or other Dex-induced TFs to atrophy genes in muscle cells.

Enhanced mitochondrial biogenesis and increased oxidative capacity are the key hallmarks of muscle adaptation to endurance exercise (Qaisar et al., 2016; Egan and Zierath, 2013). ERRγ is dominantly expressed in the slow-dominant muscles and induced in the fast muscles after exercise (Gan et al., 2013; Narkar et al., 2011). Forced expression of ERRγ in the muscle induces the formation of type I fibers and promotes endurance capacity, while the loss of ERRγ is associated with impaired oxidative metabolism and muscle endurance (Misra et al., 2017; Narkar et al., 2011; Rangwala et al., 2010). Moreover, comprehensive epigenomic analysis in the mouse muscles identified ERRγ as the causal TF for exercise-induced epigenomic reprogramming (Ramachandran et al., 2019). This evidence indicates that ERRγ is a key mediator of the metabolic and physical adaptation to exercise. Previous studies have primarily focused on the function of ERRγ in the fast-dominant muscles, and hence, its physiological relevance in the slow muscles has not been characterized. In this study, we found that LSD1 depletion promotes the adaptive response to VWR, with increased expression of ERRγ and its downstream mitochondrial metabolism genes in the Sol. We also found that VWR upregulated the protein expression of ERRγ in the Sol, and hence, our data collectively indicates that LSD1 limits the adaptive response in the slow muscle, particularly the type I fiber, by fine-tuning the exercise-induced expression of ERRγ.

Fiber-type composition in the skeletal muscle is subject to remodeling under diverse environmental stresses. Previous reports have demonstrated that Dex treatment increases the relative occupancy of type I fibers in fast-dominant muscles due to the shrinking of type IIB/X fibers but does not directly affect the growth or formation of type I fiber (Ciciliot et al., 2013; Shimizu et al., 2015). Notably, we observed an increase in the number of enlarged type I fibers in the LSD1-mKO Gas muscle. Because there were no obvious de novo fiber formation under our experimental setting (Figure 1—figure supplement 7C), it was likely that the preexisting fibers were affected by LSD1 depletion. One possible explanation is that the slow fibers in the LSD1-mKO mice became resistant to Dex. However, previous reports have shown that the type I fibers are unresponsive to GC-induced shrinking (Ciciliot et al., 2013; Shimizu et al., 2015), suggesting that this fiber type in the Gas was already resistant to Dex, regardless of the LSD1 genotype. A more plausible scenario is that type II fibers were converted to type I fibers. This possibility is backed by our previous finding that LSD1 binds to and represses a number of slow fiber-associated genes (e.g., Myh7) in C2C12 myotubes that express fast muscle-like phenotypes (Anan et al., 2018). Thus, the loss of LSD1 compromised the fast fiber identity, which in turn made the muscle cells permissive to the fiber-type conversion induced by the Dex treatment.

In our VWR study, we found that the loss of LSD1 led to an increase in the type I fiber area but not the fiber number, suggesting that these fibers were selectively enlarged. It has also been reported that VWR induces the hypertrophy of Sol with prominent effects on type IIA fibers but not on type I fibers (Momken et al., 2005; Pellegrino et al., 2005). Thus, our data indicate that LSD1 loss sensitized type I fibers to exercise-induced stress. Because LSD1-KO did not upregulate the slow fiber-associated genes, such as Myh7 and Myl2 (data not shown), LSD1 possibly suppressed the synthesis of type I fiber indirectly by negatively regulating the translation or stability of fast fiber proteins. In a previous report, the forced expression of oxidative metabolism genes in the muscle in mice enhanced type I fiber formation (Hénique et al., 2015). Thus, it is possible that the enhanced oxidative capacity in the trained LSD1-KO Sol triggered the type I fiber program. Collectively, our current study highlights the important contribution of LSD1 to the maintenance of fiber-type identity in the skeletal muscle under environmental fluctuations.

Our findings have several important implications from the biomedical standpoint. Sarcopenia is aging-associated muscle loss (Cruz-Jentoft and Sayer, 2019) and is characterized by a fast fiber-selective loss with the possible involvement of elevated GCs in aged individuals (Ciciliot et al., 2013; Clegg and Hassan-Smith, 2018). We found that LSD1 expression in mouse and human muscles declined with age and showed negative correlation with the expression of the atrophy genes (Figure 6). Revitalizing LSD1 activity using chemicals or dietary factors (e.g., riboflavin) may facilitate muscle maintenance under aging-associated stresses. Moreover, the inhibition of LSD1 in combination with aerobic exercise may be beneficial for the prevention and management of disuse-induced atrophy, in which the oxidative capacity in the slow muscles is impaired (Ciciliot et al., 2013; Ji and Yeo, 2019). The modulation of LSD1 activity in different contexts will provide new strategies for improving muscle health.

Methods

Reagents and antibodies

Dex and insulin (from bovine pancreas) were purchased from Sigma-Aldrich, and STZ was from Wako. An LSD1 inhibitor, T-3775440 hydrochloride, was from MedChemExpress. The primary antibodies used for western blot, co-IP, immunohistochemistry, and immunocytochemistry experiments were as follows: anti-LSD1 (Abcam, ab17721), anti-FOXK1 (Abcam, ab18196), anti-ERRγ (Abcam, ab12893), anti-MHC type I (Developmental Studies Hybridoma Bank [DHSB], BA-F8), anti-MHC type IIA (DHSB, SC-71), anti-Laminin (Sigma-Aldrich, L9393), anti-Akt (Cell Signaling Technology [CST], #9272), anti-Phospho-Akt (Ser473) (CST, #9271), anti-4E-BP1 (CST, #9454), anti-Phospho-4E-BP1 (Thr37/46) (CST, #2855), anti-SQSTM1(p62) (Santa Cruz, sc-28359), anti-GAPDH (Santa Cruz, sc-25778), anti-LC3(MBL, PM036), and anti-Sin3A (Santa Cruz, sc-994). The secondary antibodies used were as follows: anti-mouse IgG-horseradish peroxidase (GE Healthcare, NA931V), anti-rabbit IgG-horseradish peroxidase (GE Healthcare, NA934V or CST, #7074), Alexa Fluor 488 anti-Mouse IgG1 (Thermo Fisher, A21121), Alexa Fluor 350 anti-Mouse IgG2b (Thermo Fisher, A21140), and Cy3 anti-Rabbit IgG (Jackson ImmunoResearch, 711-165-152). The antibodies used for ChIP experiments were anti-tri-methylated histone H3K4 (Millipore, 07-473), anti-pan histone H3 (Abcam, ab1791), and normal rabbit IgG (Santa Cruz, sc-2027).

Cell culture

C2C12 mouse myoblasts were obtained from RIKEN BRC Cell Bank and maintained in Dulbecco’s Modified Eagle’s Medium (DMEM, Sigma-Aldrich) supplemented with 10% (v/v) heat-inactivated fetal bovine serum (FBS) at 37°C and 5% CO2. For myogenic induction, subconfluent C2C12 myoblasts were cultured in DMEM with 2% (v/v) horse serum, and the medium was changed every other day. For immunofluorescence study, C2C12 myotubes at differentiation day 6 were treated with 10 nM T3775440 hydrochloride or phosphate buffered saline (PBS), and then at day 8, the myotubes were treated with 1 μM Dex or ethanol for 10 hr immediately before fixation. For co-IP experiments, C2C12 myoblasts were allowed to differentiate for 8 days, followed by treatment with 1 μM insulin for 6 hr.

For CRISPR-Cas9-based genome editing, either Control Double Nickase Plasmid (Santa Cruz Biotechnology, sc-437281) or FOXK1 Double Nickase Plasmid (Santa Cruz Biotechnology, sc421682-NIC) were transfected into C2C12 cells using a Fugene 6 reagent (Promega) according to the manufacturer’s instructions. The cells were subjected to puromycin selection to establish stable lines.

H9c2 rat cardiomyocytes (ATCC) were cultured in DMEM supplemented with 10% FBS at 37°C and 5% CO2. For differentiation, cells were cultured in DMEM with 1% FBS for 7 days.

Identities of the cell lines were confirmed by gene expression profiles and morphological changes during the differentiation. In addition, C2C12 (mouse) and H9c2 (rat) cells can be distinguished by using species-specific primers. Only the cells of low passage no. were used to ensure the phenotypic consistencies. These cell lines were tested negative for mycoplasma.

Animal studies

Animal experiments were approved by the Animal Care and Use Committee of Kumamoto University (A2019-098), and conducted in accordance with its guidelines. Lsd1flox/flox; TgACTA1-CreERT2/+ conditional knockout mice (LSD1-mKO mice) were generated by crossing Lsd1-floxed mice (B6.129-Kdm1atm1.1Sho/J, Jackson Laboratories) and knock-in mice expressing a CreERT2 recombinase under the control of human ACTA1 gene promoter (Tg(ACTA1-cre/Esr1*)2Kesr/J, Jackson Laboratories). To induce Cre expression, 20 mg/mL tamoxifen in corn oil was administered intraperitoneally to 8-week-old or postnatal day 11 male LSD1-mKO or WT (Lsd1-floxed, without CreERT2) mice at a dose of 75 mg/kg for 5 consecutive days. After 1 week, RNA and protein were isolated from the tissues for qRT-PCR and western blot analyses, respectively, to evaluate the degree of LSD1 knockout (Figure 1—figure supplements 1 and 2).

To determine the effect of Dex treatment, Dex at 1 or 5 mg/kg or vehicle (PBS) was intraperitoneally administered to 8-week-old male mice for 7 consecutive days. These mice were either used for muscle function tests or sacrificed for tissue analysis after 6 hr fasting.

To evaluate the effect of insulin deficiency, STZ was dissolved in 0.1 M sodium citrate buffer (pH 4.5) to produce a concentration of 10 mg/mL. Two-hundred mg/kg STZ was administered intraperitoneally to 8-week-old male LSD1-mKO or WT mice on the fifth day of 5 consecutive days of tamoxifen administration. Insulin deficiency was confirmed 4 days later with a blood glucose reading of >350 mg/dL. After 2 weeks of STZ injection, RNA was isolated from the excised tissues for qRT-PCR.

Voluntary wheel running

For VWR test, 8-week-old control or LSD1-mKO mice were single-housed in cages with free access to a vertical running wheel (Shinano Seisakusho, Tokyo, Japan) for 38 days. After 3 days of familiarization, LSD1-KO was induced by tamoxifen injection for 5 days, consecutively, and the running distance was measured for 30 subsequent days. Before being subjected to muscle function test and tissue analysis, the mice were kept away from the running wheel and subjected to fasting for 6 hr. Sedentary controls were single-housed in similar cages but without access to a running wheel.

Muscle function tests

To assess muscle endurance, the mice were subjected to a low-intensity, run-to-exhaustion protocol on a motorized treadmill (Muromachi Kikai, Tokyo, Japan) (Fujita et al., 2018). First, the mice were allowed to familiarize with the treadmill for 10 min at 10 m/min for 2 days. The next day, the mice were placed on a treadmill moving at 10 m/min for 30 min, 11 m/min for 15 min, and 12 m/min for 15 min with a 15° incline. Finally, the speed was escalated with an increment of 1 m/min every 10 min until the mouse exhibited exhaustion. The timepoint at which the mice sat on the shock grid at the back of the treadmill for longer than 5 s was set as the endpoint.

To assess muscle strength, whole-limb grip strength was measured using a Grip Strength Meter (Columbus Instruments, Columbus, OH, USA) for mice, as previously described (Fujita et al., 2018). Peak tension (in newtons) was recorded when the mice released the grip. Two sets of 10 successive measurements were performed for each mouse and the mean maximum strength in the second set was used for data analysis.

Histological analysis

Gas and Sol muscles were dissected and immediately frozen in liquid N2-cooled isopentane. Cryosections (10 μm) were prepared using a cryostat (Leica). The sections were dried at 25°C and stored at –80°C.

For immunohistochemical analysis, the sections were fixed in acetone for 5 min at 4°C, blocked in PBS containing 3% goat serum, and 0.1% bovine serum albumin (BSA) for 2 hr at RT and incubated overnight at 4°C with a mixture of the primary antibodies (αMYH type I, αMYH type IIA, and αLaminin) in Can Get Signal immunostain Solution A (TOYOBO). After three 10 min washes in 0.1% BSA/PBS, the sections were incubated for 1 hr with a mixture of three secondary antibodies against mouse IgG2b (Alexa 350), mouse IgG1 (Alexa 488), and rabbit IgG (Cy3) in Can Get Signal immunostain Solution A, followed by three 10 min washes in 0.1% BSA/PBS. The sections were then mounted and sealed with nail polish.

The sections were imaged with a KEYENCE BZ-X800 microscope and an Olympus BX51 wide-field microscope. The cross-sectional area and numbers of type I, type IIA, and IIB+IIX fibers were determined by analyzing individual fibers that were visualized by laminin staining. Quantification was done using a BZ-X Analyzer software (KEYENCE). Example images used for detecting individual fibers are shown in Figure 1—figure supplement 6B.

SDH staining was performed by exposing the cryosections to a mixture of 130 mM succinate, 1.5 mM nitroblue tetrazolium, 0.2 M NaH2PO4, and 0.2 M Na2HPO4 for 30 min at 37°C. The sections were then washed once with 70% and 95% ethanol and twice with 99.5% ethanol. Next, after washing with xylene, the sections were dried and mounted in malinol (Muto Pure Chemicals, Tokyo, Japan).

Western blot and co-IP analyses

To prepare total cell lysates, the cells were collected and suspended in a sample buffer (0.1 M Tris-HCl, pH 6.8, 4% sodium dodecyl sulfate [SDS], 0.1 M dithiothreitol, 20% glycerol, and 0.2% bromophenol blue). To prepare tissue lysates, the frozen tissues were crashed in RIPA buffer (50 mM Tris-HCl, 150 mM NaCl, 2 mM EDTA, 1% NP-40, 0.1% SDS, and 0.5% sodium deoxycholate) by FastPrep-24 5G (MP Biomedicals, Santa Ana, CA, USA) at 6.0 m/s for 40 s with CoolPrep mode, then suspended in the sample buffer. Following sonication and centrifugation at 16,710 × g for 5 min at 4°C, the supernatant was collected and used for western blotting. Protein samples were electrophoresed on an SDS polyacrylamide gel and then transferred to a nitrocellulose membrane (Amersham Protran Premium; GE Healthcare, Waukesha, WI, USA) using a semidry method. After blocking for 1 hr using 5% skim milk or 5% BSA dissolved in PBS containing 0.3% Tween 20, the membrane was incubated overnight at 4°C with primary antibodies in Can Get Signal solution (Toyobo, Osaka, Japan). The secondary antibodies used were anti-mouse IgG and anti-rabbit IgG-horseradish peroxidase. The blots were incubated for 1 min with Western Lightning Plus-ECL solution (PerkinElmer, Waltham, MA, USA) and visualized using Image-Quant LAS4000 Mini (GE Healthcare). For quantification, the band densities were determined using ImageJ (National Institutes of Health, Bethesda, MD, USA) software.

For co-IP experiments, the cells were lysed in an IP buffer (50 mM Tris-HCl, pH 8.0, 5 mM EDTA, 150 mM NaCl, 0.5% NP-40, and 0.5% Triton X-100) and incubated with specific antibodies, followed by pulldown using Dynabeads protein A/G (Life Technologies) for subsequent western blot analysis.

Separation of nuclear and cytoplasmic fractions

Separation of nuclear and cytoplasmic fractions of C2C12 myotubes was performed as previously described (Kim et al., 2010). Briefly, cells were washed twice with cold PBS and harvested with a scraper, and then mixed with nuclear extraction buffer (20 mM HEPES, pH 7.5, containing 150 mM NaCl, 10% glycerol, 1 mM EDTA, 0.5% Triton X-100, 1 mM Na3VO4, 10 mM NaF, 1 mM PMSF, and 0.01% of protease inhibitor cocktail) and incubated for 30 min on ice. The cell suspension was passed through a 23-gauge needle 30 times and then centrifuged. The pellet and supernatant from the centrifugation are referred to as the nuclear and cytoplasmic fractions, respectively. The nuclear fraction was washed three times with nuclear extraction buffer and then suspended in the sample buffer.

Quantitative RT-PCR

Total RNA was extracted from cells or frozen mouse tissues using a TRIzol reagent (Invitrogen). cDNA samples were produced using the ReverTra Ace qRT-PCR Kit (TOYOBO). qPCR-RT was performed by the SYBR green method using a Thunderbird qPCR Mix (TOYOBO) and a StepOne Plus Sequence Detector (Applied Biosciences). Data are presented as the mean ± standard deviation (SD). Fold differences among groups were calculated by the ΔΔCt method. 36B4 (Rplp0) gene was used as an internal control. The primers are listed in Supplementary file 3.

Poly (A) RNA-seq analysis

Total RNA was extracted from the Gas and Sol muscles using a Trizol reagent. mRNA was purified using a NEBNext Poly(A) mRNA Magnetic Isolation Module (New England Biolabs [NEB], Ipswich, MA, USA), and the sequencing libraries were synthesized using a NEBNext Ultra II RNA Library Prep Kit for Illumina (NEB). Sequencing was performed using a NextSeq 500 Sequencer (Illumina, San Diego, CA, USA) with 75 bp single-end reads. The resulting reads were aligned to the mm10 reference genome by using STAR (Dobin et al., 2013). For the dataset from Dex administration experiment, Edge R was used for the quantification of gene expression changes (Robinson et al., 2010), in which a false discovery rate cut-off of 0.05 was used to determine differentially expressed genes between WT and KO mice. Heatmaps were generated by using Heatmapper (Babicki et al., 2016). For the VWR dataset, DESeq2 was used for the quantification of gene expression changes (Love et al., 2014; Takase et al., 2019). The top 1000 differentially expressed genes between WT and KO mice were extracted from the genes with mean transcripts per million (TPM) >1 by using median absolute deviation (MAD). K-means clusters and heatmaps were generated by using ComplexHeatmap (Gu et al., 2016). Gene Ontology analysis was performed by using the Molecular Signatures Database (MSigDB) v7.4 (Subramanian et al., 2005).

Motif enrichment at the promoter regions was examined using the findMotifs.pl module in HOMER (V4.9.1) (Heinz et al., 2010). The parameters used for the Dex test were ‘start –300, end +50, motif length 8, 10, 12’, whereas those for the VWR test were ‘start –2000, end +2000, motif length 10, 12’. For the Dex administration dataset, the LSD1-dependent atrophy genes and all upregulated genes in the LSD1-KO mice were subjected to motif analysis.

ChIP-qPCR

To prepare a ChIP sample, frozen Gas muscles from two mice were minced in cooled PBS by using FastPrep-24 5G (MP Biomedicals), and the lysate was crosslinked with 1% formaldehyde. Next, to isolate the nuclei, tissue lysates were homogenized in cell lysis buffer (5 mM PIPES, 85 mM KCl, and 0.5% NP-40) with a Dounce homogenizer using a loose pestle. The isolated nuclei were resuspended in SDS-lysis buffer (1% SDS, 10 mM EDTA, and 50 mM Tris-HCl) and sonicated in two steps (Kohrogi et al., 2021): first with a probe-type sonifier (Branson) and second with a water bath sonifier, Picoruptor (Diagenode). For the C2C12 ChIP sample preparation, cells were dual crosslinked with 1% formaldehyde and disuccinimidyl glutarate (Thermo Fisher Scientific) to increase the stability of protein-DNA complexes. Following cell lysis, cells were treated with 4 U of Micrococcal nuclease per 1×106 nuclei at 37°C for 20 min, then sonicated by Bioruptor (Diagenode) for 10 min. Chromatin fragments were incubated at 4°C overnight with appropriate antibodies, followed by a pulldown using protein A/G-conjugated agarose beads. Purified DNAs were subjected to qPCR using the primer sets listed in Supplementary file 4.

Immunocytochemistry

The cells were fixed with 4% paraformaldehyde in PBS for 15 min at 25°C. The cells were washed three times with PBS for 5 min and then permeabilized with PBS containing 0.5% Triton X-100 for 5 min on ice. The cells were blocked in 0.5% BSA/PBS at RT and incubated overnight at 4°C with a mixture of primary antibodies in 0.2% BSA/PBS. After three 10 min washes in 0.2% BSA/PBS, the cells were incubated for 1 hr with a mixture of secondary antibodies against mouse IgG1 (Alexa 488) and rabbit IgG (Cy3) in 0.2% BSA/PBS, followed by three washes in PBS. The cells were then mounted and sealed with nail polish. DNA was counterstained with 1 μg/mL 5,6-diamidino-2-phenylindole. Quantitative analysis was performed with BZ-X Analyzer software (KEYENCE).

Construction of lentiviral expression vectors and their transduction

The plasmids for constructing and packaging lentiviral vectors were obtained from the RIKEN Bioresource Center (https://dna.brc.riken.jp/en/rvden/cfmen). Double-stranded DNA harboring an shRNA sequence for FOXK1 was inserted into an entry vector, pENTR4-H1. Using the Gateway recombination method, an shRNA cassette was cloned into CS-RfA-EVBsd, a lentiviral plasmid. To produce lentivirus particles, packaging (pCAG-HIVgp) and envelope (pCMV-VSV-G-RSV-Rev) plasmids were co-transfected into Lenti-X 293T cells (Takara Bio) using Fugene 6 reagent (Promega). The culture supernatant was centrifuged, filtered, and stored at –80°C until use. C2C12 cells were infected with CS-RfA-EVBsd shFoxk1 (shFoxk1) or CS-RfA-EVBsd shGL3 (shGL3) (control shRNA targeting firefly luciferase gene) viruses, and then subjected to blasticidin S selection to establish stable lines. The short hairpin RNA target sequences were as follows: shFoxk1, 5’-ATCCAGTTCACATCGCTAT-3’; shGL3, 5’- CTTACGCTGAGTACTTCGA -3’.

For the forced expression of LSD1 in H9c2 cells, we constructed a lentiviral vector harboring human LSD1 cDNA under control of human elongation factor-1α (EF-1α) promoter. For this, we cloned FLAGx3-hLSD1 into a lentiviral plasmid CSII-EF-MCS (RIKEN) together with blasticidin resistance gene. After transduction, cells were selected with 10 µg/mL blasticidin S before use.

Tabula Muris Senis

Single cell RNA-seq data from mice of different ages (3, 18, and 24 m) were obtained from the Tabula Muris Senis Project (Tabula Muris Consortium, 2020). FACS-sorted single-cell data from TA muscle (designated as limb muscle) were downloaded at GitHub (https://github.com/czbiohub/tabula-muris-senis, copy archived at swh:1:rev:18219fd8d2acaf9c8d46befbdfb44cd75fadba12; Zhang, 2021) and analyzed using Scanpy (https://scanpy.readthedocs.io/en/stable/). Muscle cells were extracted according to the expression of Myod1 or Pax7.

Human muscle transcriptome data

LSD1 mRNA expression in skeletal muscles of 55- to 80-year-old men and women were analyzed using a publicly available RNA-seq data (Tumasian et al., 2021). Expression values in TPM were used for linear regression analysis by JMP software (version 10.0.2d1).

Statistical analyses

Sample size was determined based on our pilot experiments and previous studies in the related research area. Data are presented as the mean ± SD. Equality of variance was determined by using F test. If the variances were equal, two-tailed Student’s t-test was used, and if not, two-tailed Welch’s t-test was used.

Materials availability statement

Newly created materials can be provided following the material transfer agreement of Kumamoto University.

Appendix 1

Appendix 1—key resources table
Reagent type (species) or resourceDesignationSource or referenceIdentifiersAdditional information
Strain, strain background (Mus musculus)Lsd1-floxed mice (B6.129-Kdm1atm1.1Sho/J)Jackson LaboratoriesStrain #:023969
RRID:IMSR_JAX:023969
Strain, strain background (Mus musculus)Tg(ACTA1-cre/Esr1*)2Kesr/JJackson LaboratoriesStrain #:025750
RRID:IMSR_JAX:025750
Strain, strain background (Mus musculus)Lsd1flox/flox; TgACTA1-CreERT2/+ conditional knockout mice (LSD1-mKO mice)This paperN/ASee ‘Animal studies’ in Methods
Strain, strain background (Mus musculus)C57BL6/J miceCharles River LaboratoriesN/A
Cell line (Mus musculus)C2C12RIKENRCB0987
Cell line (Ruttus norvegicus)H9c2ATCCCRL-1446
AntibodyRabbit polyclonal anti-LSD1abcamab17721WB (1:500), ChIP (5 μg), Co-IP (3 μg)
AntibodyRabbit polyclonal anti-Foxk1abcamab18196WB (1:1000), ChIP (5 μg), Co-IP (3 μg)
AntibodyRabbit monoclonal anti-ERR gammaabcamab128930WB (1:1000)
AntibodyMouse monoclonal anti-MHC type IDHSBBA-F8IHC (1:50)
AntibodyMouse monoclonal anti-MHC type IIADHSBSC-71IHC (1:50)
AntibodyRabbit polyclonal anti-LamininSigma-AldrichL9393IHC (1:500)
AntibodyRabbit polyclonal anti-AktCell Signaling Technology#9272WB (1:1000)
AntibodyRabbit polyclonal anti-Phospho-AktCell Signaling Technology#9271WB (1:1000)
AntibodyRabbit polyclonal anti-LC3MBLPM036WB (1:1000)
AntibodyMouse monoclonal anti- Phospho-4E-BP1 (Thr37/46)Cell Signaling Technology#2855WB (1:1000)
AntibodyRabbit polyclonal anti-4E-BP1Cell Signaling Technology9454WB (1:1000)
AntibodyMouse monoclonal anti- SQSTM1(p62)Santa Cruzsc-28359WB (1:1000)
AntibodyRabbit polyclonal anti-GAPDHSanta Cruzsc-25778WB (1:6000)
AntibodyRabbit polyclonal anti-trimethyl-Histone H3 (Lys4)Millipore07-473ChIP (3 μg)
AntibodyRabbit polyclonal anti-Histone H3abcamab1791WB (1:5000)
Recombinant DNA reagentControl Double Nickase PlasmidSanta Cruz Biotechnologysc-437281
Recombinant DNA reagentFOXK1 Double Nickase PlasmidSanta Cruz Biotechnologysc-437282-NIC
Recombinant DNA reagentpENTR4-H1RIKEN BRC DNA BANKRDB04395
Recombinant DNA reagentpCAG-HIVgpRIKEN BRC DNA BANKRDB04394
Recombinant DNA reagentpCMV-VSV-G-RSV-RevRIKEN BRC DNA BANKRDB04393
Recombinant DNA reagentCS-RfA-EVBsdRIKEN BRC DNA BANKRDB06090
Recombinant DNA reagentCS-RfA-EVBsd shFoxk1This paperN/ASee ‘Construction of lentiviral expression vectors and their transduction’ in Methods
Recombinant DNA reagentCS-RfA-EVBsd shGL3This paperN/ASee ‘Construction of lentiviral expression vectors and their transduction’ in Methods
Recombinant DNA reagentCSII-EF-3xFLAG-hLSD1-BsdThis paperN/AA lentiviral construct used for the forced expression of human LSD1
Sequence-based reagentPrimers for qRT-PCRThis paperN/ASee Supplementary file 3
Sequence-based reagentSee Table S2 for primers for ChIP-qPCRThis paperN/ASee Supplementary file 4
Commercial assay or kitNEBNext Ultra DNA Library Prep KitNew England BiolabsE7370S
Commercial assay or kitNEBNext Poly(A) mRNA Magnetic Isolation ModuleNew England BiolabsE7490S
Chemical compound, drugTamoxifenSigma-AldrichT5648
Chemical compound, drugDexamethasoneSigma-AldrichD4902
Chemical compound, drugInsulinnacalai tesque19251-24
Chemical compound, drugStreptozotocinWako197-15153
Chemical compound, drugT-3775440 hydrochlorideMedChemExpressHY-103085
Chemical compound, drugFuGENE 6 Transfection ReagentPromegaE2691
Chemical compound, drugPuromycin dihydrochlorideSigma-AldrichP8833
Chemical compound, drugBlasticidin SFUJIFILM Wako Chemicals029-18701
Chemical compound, drugCan Get Signal immunostain Solution ATOYOBONKB-501
Chemical compound, drugmalinolMuto Pure Chemicals2009-1
Chemical compound, drugCan Get Signal Solution 1TOYOBONKB-201
Chemical compound, drugCan Get Signal Solution 2TOYOBONKB-301
Chemical compound, drugTRIzol RNA Isolation ReagentThermo Fisher Scientific15596026
Chemical compound, drugReverTra Ace qPCR RT Master MixTOYOBOFSQ-201
Chemical compound, drugTHUNDERBIRD SYBR qPCR MixTOYOBOQPS-201
Software, algorithmBZ-X Analyzer softwareKEYENCEN/A
Software, algorithmHomer v4.9.1Heinz et al., 2010http://homer.ucsd.edu/homer/
Software, algorithmSTARDobin et al., 2013https://github.com/alexdobin/STAR; Dobin, 2023
Software, algorithmEdge RRobinson et al., 2010http://bioconductor.org/packages/release/bioc/html/edgeR.html
Software, algorithmHeatmapperBabicki et al., 2016http://www.heatmapper.ca/expression/
Software, algorithmDEseq2Love et al., 2014https://bioconductor.org/packages/release/bioc/html/DESeq2.html
Software, algorithmComplexHeatmapGu et al., 2016https://github.com/jokergoo/ComplexHeatmap; Gu, 2022
Software, algorithmMolecular Signature Database v7.4Subramanian et al., 2005https://www.gsea-msigdb.org/gsea/msigdb/
Software, algorithmHOMER V4.9.1Heinz et al., 2010http://homer.ucsd.edu/homer/
Software, algorithmJMP (version 10.0.2d1)SAS Institute IncN/A
Software, algorithmGraphPad Prism version 8.4.3GraphPad SoftwareN/A
Software, algorithmIntegrative Genomics Viewer (IGV)James T et al., 2011https://software.broadinstitute.org/software/igv/

Data availability

RNA-seq data generated in this study were deposited in the Gene Expression Omnibus (accession No. GSE198911).

The following data sets were generated
    1. Araki H
    2. Hino S
    3. Nakao M
    (2022) NCBI Gene Expression Omnibus
    ID GSE198911. polyA RNA-seq of WT and LSD1-mKO mice under two different environmental stresses.
The following previously published data sets were used
    1. Tumasian et al
    (2014) NCBI Gene Expression Omnibus
    ID GSE164471. Skeletal muscle transcriptome in healthy aging.

References

  1. Book
    1. Gilbert SF
    2. Epel D
    (2009)
    Ecological Developmental Biology: Integrating Epigenetics, Medicine, and Evolution
    Sinauer Associates Inc.

Decision letter

  1. Christopher Cardozo
    Reviewing Editor; Icahn School of Medicine at Mount Sinai, United States
  2. David E James
    Senior Editor; University of Sydney, Australia

In the interests of transparency, eLife publishes the most substantive revision requests and the accompanying author responses.

[Editors' note: this paper was reviewed by Review Commons.]

https://doi.org/10.7554/eLife.84618.sa1

Author response

General Statements [optional]

We thank the reviewers for productive comments. Many of the experiments and analyses suggested have already been carried out. For the remaining experiments, we have already gathered necessary materials. Thus, we can add highly relevant revisions to the next version within a short period.

Description of the planned revisions

Insert here a point-by-point reply that explains what revisions, additional experimentations and analyses are planned to address the points raised by the referees.

Reviewer #1:

Major #1

This study primarily uses the genetic mouse model in which LSD1 gene is inactivated after tamoxifen injection in 8 weeks old mice, as shown in supplemental figure 1 B and C. 8 weeks after birth postnatal growth of muscle is not complete and the contribution of satellite cells to muscle growth is still significant. Therefore the timing of tamoxifen injection used cannot discriminate if the observed phenotype involves the function of LSD1 during the post-natal growth of the muscle or in the muscle fibers or both. One way to demonstrate the real contribution of LSD1 in the maintenance of muscle fibers plasticity under environmental stress would be to inject Tamoxifen later (around 10-12 weeks of age), in order to remove a possible bias caused by the contribution of satellite cells during the post-natal growth. At least key findings should be confirmed at this later stage.

In this study, we used ACTA1-CreERT mice to conditionally knockout LSD1 in the skeletal muscle. The ACTA1 promoter is derived from human muscle actin gene, which is not expressed in the satellite cells, and has been widely used for the transgene expression in myofibers (Stantzou et al. Development 2017). Thus, the inactivation of LSD1 occurs in the existing myofibers, and alterations in satellite cell function, if any, would be indirect effects of the loss of LSD1 in mature myocytes or differentiating myoblasts.

To test whether postnatal muscle growth was affected in our LSD1-mKO mice, we administrated tamoxifen (4OHT) to pre-weaning mice (11 days old). LSD1 depletion did not affect the expression of muscle fiber genes, when muscle tissues were isolated from mice 11 days after the start of 4OHT (Figure 1—figure supplement 2).

These evidences exclude the contribution of satellite cells in the phenotypes observed in the LSD1-mKO mice. Figure 1—figure supplement 2 has been added in the revised manuscript, and described in the Results section (P.5).

Major #2

LSD1 m-KO muscles seem to have more type I and IIA fibers than WT, even without DEX treatment. Is it possible to quantify the results in supplemental figure 4C?

As suggested, we quantitatively analyzed the fiber type compositions in Figure 1—figure supplement 6C (former Figure S4C) using the data from WT (n=4) and LSD1-mKO (n=5) mice (Figure 1—figure supplement 6D ). We did not find a significant difference between these mice, confirming our finding that the loss of LSD1 accelerates the Dex-driven phenotypic changes. We also analyzed the Soleus muscle from untreated mice, and found no obvious difference between the WT and KO (Figure 1—figure supplement 7A-B ).

Major #3

The effect on fiber type is convincing, while variations in gene expression are of quite low amplitude. However, the atrophy should be induced by other means to ensure that the effects are specific to GC/nuclear receptors pathways; Denervation? Starvation? Not all the experiments need to be repeated, just key results such as, for example, exacerbation of atrophy in LSD1 m-KO, Foxk1 increase.

We agree that testing alternative atrophy models is important for generalizing our findings. For this, we employed a model for diabetes-related muscle atrophy. A pro-diabetic agent streptozotocin (STZ) disturbs the function of pancreatic islet leading to fast-fiber atrophy (O’Neill et al. Diabetes 2019). LSD1-mKO did not affect the muscle weight in STZ-treated mice (Figure 2—figure supplement 3). Consistently, there were no major difference in the expression of atrophy genes in STZ-treated WT and LSD1-mKO mice (Figure 2—figure supplement 3 ). These results suggest that the LSD1 function depends on the source of atrophy-inducing stress, and that the loss of LSD1 sensitized the muscle to GC-mediate signaling.

Major #4

Autophagy data: the effect on the LC3I/LC3II ratio are modest. The autophagy part should be removed or completed with additional data to convincingly show that autophagy is affected. Links between LSD1 and mTOR have been published, so the mTOR pathway could be investigated in the model (S6k, S6 and 4EBP1 phosphorylation). Given AKT levels and phosphorylation are affected by the absence of LSD1 + DEX, it can be predicted that mTOR activity will change.

We have analyzed the expression of additional autophagy markers p62 and phosphorylated 4EBP1. Consistent with the upregulated expression of atrophy genes and increased LC3I/II ratio, LSD1-mKO mice had elevated levels of p62 and phosphorylated 4EBP1 (Figure 1—figure supplement 5B and D,). Altogether, the data suggest that Dex-induced muscle atrophy was exacerbated by the loss of LSD1. These data are added as Figure 1—figure supplement 5B and D, and described in the Results section (P.6).

Major #5

The ability of LSD1 to retain FOXK1 in the nucleus is an important information that should be better supported experimentally. In the absence of such information, no mechanism can be proposed for the effect of LSD1 of FOXK1. The immunofluorescence images provided are not convincing and moreover they could be interpreted by a reduction in the level of FOXK1 protein (degradation?) rather than by a nuclear exclusion in the presence of DEX. This point should be addressed, authors could include western blot of nuclear and cytoplasmic fractions to better quantify the nuclear level of FOXK1 in absence of LSD1.

As suggested, we performed WB analyses of Foxk1 using nuclear and cytoplasmic fractions. Our results show that LSD1 inhibition combined with Dex treatment led to a reduced nuclear/cytoplasmic ratio (Figure 3—figure supplement 1B ). This data is consistent with our immunofluorescence data, and altogether suggest that LSD1 promotes the nuclear retention of Foxk1. The data has been added as Figure 3—figure supplement 1B and described in the Results section (P.9).

Major #6

The absence of centralized nuclei indicates that there is no fiber regeneration but it does not exclude the possibility that satellite cells were recruited to existing fibers and thus participated to hypertrophy. To eliminate this possibility, the average nuclei/cytoplasm volume should decrease if hypertrophy results from increased protein synthesis and not myonuclei accretion. This should be checked.

We histologically analyzed the sections of Gas muscles after Dex treatment and found that there is no evidence of central nuclei in either WT or KO mice (Figure 1—figure supplement 7C).

As mentioned above (Major #1), it is unlikely that the satellite cell function was responsible for the enhanced atrophic phenotype.

Major #7

The upregulation of ERRγ in the absence of LSD1 is convincing in the VWR conditions. ERRγ level should be evaluated in the sedentary LSD1 KO mice.

We have analyzed the expression of ERRg in sedentary mice, and found no significant difference between WT and KO mice (Figure 5—figure supplement 1A ). This suggests that the loss of LSD1 in combination with VWR training led to the increased expression of ERRg. This data is included in the revised manuscript as Figure 5—figure supplement 1A and described in the Results section (P.12).

Minor #1

There is a clear difference in the number of mouse replicates between treated (Dex or VWR) and non-treated mice, regardless the genotype. Experiments with non-treated mice lack adequate numbers to make a definitive conclusion. For example, there is a huge spread in the data in Figure 1 B and 4 D. If the number of animals would have been increased, would the absence of difference hold up?

We increased the number of non-treated animals in Figures 1B and 4B as suggested. Nonetheless, we did not find any significant differences in the muscle weight (Figures 1B and 4B). These changes are reflected on original Figures 1B and 4B. For Figure 1—figure supplement 4A-B, we used violin plots to clearly visualize the differences among groups.

Minor #2

The authors claim that: "Consistent with the results of the augmented endurance capacity, the Sol muscle in the KO mice showed enhanced succinate dehydrogenase (SDH) staining, indicating that the number of oxidative fibers increased (Figure 4F and Supplemental Figure 8F)". However, supplemental figure 8 D indicates that the number of type I fibers does not change compared to WT. Authors should clarify this statement.

Indeed, we found that the area of type I fiber but not the number was increased in the LSD1-KO Sol (Figure 4D and Figure 4—figure supplement 1D). Because SDH staining reflects the OXPHOS capacity in all fiber types, it is possible that the OXPHOS capacity in the fibers other than type I had been augmented by LSD1-KO. Thus, for clarification, we will change the statement as follows: OXPHOS capacity of Sol was enhanced by the loss of LSD1 (Results section, P.11).

Reviewer #2

Methods

1. The authors used the Cre-lox system with tamoxifen to generate skeletal muscle-specific LSD1 KO mice. It is clear that both the mRNA and protein levels of LSD1 in various muscles were dramatically reduced, but there is still some LSD1 expressed in skeletal muscle, especially in Sol muscle (Supplemental Figure 1C). The author needs to think about whether it is appropriate to use the term "LSD1 knockout" or "LSD1 deficiency".

We thank the reviewer for this comment. In this study, we crossed LSD1-floxed mice with ACTA1-creERT mice. This enables the deletion of critical exons of LSD1 in mature myocytes and myogenic precursors that have initiated the differentiation program. LSD1 is a ubiquitously expressed gene, and it is known that immature myogenic cells (e.g., satellite cells, Tosic et al. Nat Commun. 2018) and other non-myogenic cells such as hematopoietic and vascular cells abundantly express LSD1 (Kerenyi et al. ELife 2013, Yuan et al. Biochem Pharmacol. 2022). Thus, it is likely that LSD1 expression by these cell types were detected in our whole muscle western blots. We added these statements in the text for clarification (Results, P.5).

Results

2. To identify the transcriptional regulators that mediate the regulation of atrophy-associated genes by LSD1, the authors performed motif analyses on the promotor regions of upregulated genes in LSD1-mKO Gas. Based on the results and other reports, they focused on Foxk1 and proved LSD1 and Foxk1 cooperatively regulate the atrophy transcriptome in the presence of Dex. However, Figure 3C showed that a transcription factor Nfatc1 is also reduced in Sol muscle similar to Foxk1. Also, other studies demonstrated that the transcription factor NFATc1 controls fiber type composition and is required for fast-to-slow fiber type switching in response to exercise in vivo. More specifically, NFATc1 inhibits MyoD-dependent fast fiber gene promoters by physically interacting with the N-terminal activation domain of MyoD and blocking recruitment of the essential transcriptional coactivator p300 (Cell Rep. 2014 Sep 25; 8(6): 1639-1648). Furthermore, it has been reported that LSD1 Controls Timely MyoD Expression via MyoD Core Enhancer Transcription (Cell Rep. 2017 Feb 21;18(8):1996-2006. doi: 10.1016/j.celrep.2017.01.078). It is unclear how the authors exclude Nfatc1 for the LSD1-mediated effects in different muscle fibers. Further experiments may be necessary to exclude Nfatc1.

We thank the reviewer for an insightful comment. In addition to Foxk1, we tested the involvement of NFATc1 in the gene regulation under LSD1-depleted state. We treated C2C12 with an LSD1 inhibitor S2101 in combination with a calcium ionophore that promotes the transcriptional function of NFATc1 by inducing its nuclear localization (Meissner et al. J Cell Physiol. 2007). While LSD1 inhibition promoted the expression of Pgc1a and Myh7, ionophore treatment had no additive effects. Because we found a physical association of Foxk1 with LSD1, we focused on the functional involvement of Foxk1 in LSD1-mediated repression of atrophy genes. We recently performed an ATAC-seq analysis in Dex-treated muscle, and found that the Foxk1 motif but not the NFATc1 motif was enriched in the LSD1-KO-specific open chromatin regions. This data further suggests the significant contribution of Foxk1 in the transcriptional regulation under LSD1 depletion.

3. In figure 3D, only merged images were colored. It would be better to show colored images for Foxk1 and DAPI.

We replaced the images with the colored ones in Figure 3D.

4. Immunofluorescence analysis in C2C12 myotubes showed that Dex exposure reduced the nuclear retention of Foxk1, which was further promoted by the addition of T-3775440, an LSD1 inhibitor (Figure 3D). The author also used Foxk1-KO C2C12 myotubes to prove LSD1 and Foxk1 cooperation to regulate the expression of type I /IIA fiber and atrophy genes in Foxk1-KO cells. Are the effects of LSD1 dependent on Foxk1 or synergistically acting with Foxk1? The treatment of LSD1 inhibitor in Foxk1-KO C2C12 may be helpful to answer this question.

As suggested, we tested the combination effect of LSD1 inhibition and Foxk1-KO. The treatment of Foxk1-KO cells with T-3775440 did not show additional effects on the expression of atrophy and slow fiber genes (Figure 3—figure supplement 2B). Next, we analyzed the chromatin association of LSD1 in Foxk1-KO cells by a ChIP assay, and found that the occupancy of LSD1 at atrophy gene loci was diminished in Foxk1-KO cells (Figure 3G and Figure 3—figure supplement 2C). These results collectively suggest that the LSD1 largely depend on Foxk1 to control atrophy genes. These data are included in the revised version as Figure 3G and Figure 3—figure supplement 2B-C, and described in the text (Results, P. 10).

5. In supplementary figure 2, body weight in the mKO+Dex group was reduced in comparison to that of WT+Dex. How about the body weight of mKO mice without Dex injection compared to that of WT? This data will be helpful to understand the effect of muscle-specific LSD1 deficiency on whole-body energy balance.

We measured the body weight of untreated mice, and found that there is no genotype effect (Figure 1—figure supplement 3D-E ). Thus, we think that LSD1-mKO alone does not influence the whole-body energy balance. We included this data as Figure 1—figure supplement 3D-E.

6. The authors analyzed the size distribution of myofibers and mentioned that large type I and type IIA fibers preferentially increased in the LSD1-mKO muscle, whereas large type IIB + IIX fibers decreased (Supplemental Figure 4, B, E, and F). It is better to show the results of statistics. If no significance were found, it should be mentioned in the result section.

We have performed statistical analyses on Figure 1—figure supplement 7D-E (Supplemental Figure 4E and F in the previous version), and found that a fraction of large type I fibers was significantly larger in KO mice. This result is indicated in the figure.

7. Page 11, To reveal the genes regulated by LSD1 under the VWR condition, the authors performed additional RNA-seq analysis using Sol muscle. The non-hierarchical clustering analysis was informative and showed signaling pathways related to ‘mitochondrion’, ‘mitochondrion organization’, and ‘oxidative phosphorylation’ were altered in the Sol muscle deficient in LSD1 under the VWR condition (Figure 5B). However, it is unclear why they focus on Err-γ to explain LSD1-KO phenotypes in Sol muscle. Is this gene also derived from RNA seq? It is better to show whether Err-γ expression is also significantly altered based on RNA seq data.

Indeed, ERRg was upregulated by LSD1-KO+VWR and was included in the Cluster 6 genes together with the OXPHOS and mitochondria-related genes (Figure 5A). These data prompted us to focus on ERRg as a potential factor that explains the LSD1-KO phenotype.

8. The authors claim that LSD1 serves as an "epigenetic barrier" that optimizes fiber type-specific responses and muscle mass under stress conditions. This claim is derived from the loss of function studies. To generalize the functions of LSD1, the gain of function studies will be also necessary. Adding the characteristics of LSD1 overexpression in C2C12 cells will further improve the quality of the manuscript.

We agree that the gain of function further strengthens the quality of our study. As suggested, we introduced LSD1 into C2C12 cells using a lentivirus vector. However, the net expression of endo- and exogenous LSD1 (mouse and human, respectively) were similar between mock- and LSD1-introduced cells (data not shown). This suggests that LSD1 protein level was maintained at a certain level in these cells.

Because Esrrg is completely silenced in C2C12 cells, it is difficult to monitor ERRg-mediated gene regulation in these cells. To overcome this, we made use of a rat cardiomyocyte cell line H9c2, in which ERRg is functionally involved in differentiation (Sakamoto et al. Nat Commun 2022). We first inhibited LSD1 using the inhibitor T-3775440 in these cells, and found that the expression of Esrrg and its downstream target Perm1 was significantly upregulated (Figure 5C, left). We next overexpressed LSD1 in these cells and found a downregulation of Esrrg and Perm1 (Figure 5C, right). These data are highly consistent with our in vivo data that LSD1 represses Esrrg in slow fiber enriched muscle. These data are added as Figure 5—figure supplement 1C-E, and described in the Results section (P. 12-13).

Discussion

9. The authors mentioned supplementary figure 10 only at the end of the manuscript of the Discussion section (page 15) without a specific explanation of the figures in the result section. The data are important in that LSD1 expression in human muscles declined with age and showed a negative correlation with the expression of the atrophy gene. It should be presented in the result section with a more detailed description.

We agree that these data are important and need further explanations. We added the description in details in the Results section (P.13) and moved the entire figure to the main figure (Figure 6).

10. There are other studies to examine LSD1 and muscle regeneration or functions (e.g. Nat Commun 9, 366 (2018). https://doi.org/10.1038/s41467-017-02740-5). More discussion to compare the current study and other studies will be necessary.

We thank the reviewer for this comment. We added a discussion about the phenotypic difference between our mice and the previously reported LSD1-KO mice (Tosic et al. 2018) (Discussion, P.14).

Description of the revisions that have already been incorporated in the transferred manuscript

Please insert a point-by-point reply describing the revisions that were already carried out and included in the transferred manuscript. If no revisions have been carried out yet, please leave this section empty.

As described above, a number of experiments and analyses requested by the reviewers have already been done. Those will be incorporated in the next version.

Description of analyses that authors prefer not to carry out

Please include a point-by-point response explaining why some of the requested data or additional analyses might not be necessary or cannot be provided within the scope of a revision. This can be due to time or resource limitations or in case of disagreement about the necessity of such additional data given the scope of the study. Please leave empty if not applicable.

As described above, for some comments, there are alternative solutions other than the ones requested by the reviewers.

Rev#1, Major#1, Major#6

In addition to the experiments using young (8wks) mice, we performed experiments using pre-weaning mice (11days), in which satellite cells actively contribute to the muscle growth. In this experiment, there were no differences in the muscle weight between WT and KO mice. In addition, our ACTA1-creERT-mediated LSD1 depletion system, the LSD1 gene was inactivated only after the onset of differentiation. From these evidences, we concluded that satellite cells are not the primary source of the LSD1-KO phenotype.

Thus, we do not think it necessary to perform experiments with older mice or reanalyze the nuclei/cytoplasm volume to test the possible contribution of satellite cells to the phenotype.

https://doi.org/10.7554/eLife.84618.sa2

Article and author information

Author details

  1. Hirotaka Araki

    1. Department of Medical Cell Biology, Institute of Molecular Embryology and Genetics, Kumamoto University, Kumamoto, Japan
    2. Department of Metabolic Medicine, Faculty of Life Sciences, Kumamoto University, Kumamoto, Japan
    Contribution
    Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Visualization, Writing - original draft
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-8679-948X
  2. Shinjiro Hino

    Department of Medical Cell Biology, Institute of Molecular Embryology and Genetics, Kumamoto University, Kumamoto, Japan
    Contribution
    Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Writing - original draft, Project administration, Writing - review and editing
    For correspondence
    s-hino@kumamoto-u.ac.jp
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-5748-016X
  3. Kotaro Anan

    Department of Medical Cell Biology, Institute of Molecular Embryology and Genetics, Kumamoto University, Kumamoto, Japan
    Contribution
    Conceptualization, Investigation, Methodology
    Competing interests
    No competing interests declared
  4. Kanji Kuribayashi

    Department of Medical Cell Biology, Institute of Molecular Embryology and Genetics, Kumamoto University, Kumamoto, Japan
    Contribution
    Formal analysis, Investigation
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-2666-1141
  5. Kan Etoh

    Department of Medical Cell Biology, Institute of Molecular Embryology and Genetics, Kumamoto University, Kumamoto, Japan
    Contribution
    Formal analysis, Visualization, Methodology
    Competing interests
    No competing interests declared
  6. Daiki Seko

    1. Department of Muscle Development and Regeneration, Institute of Molecular Embryology and Genetics, Kumamoto University, Kumamoto, Japan
    2. Department of Molecular Bone Biology, Graduate School of Biomedical Sciences, Nagasaki University, Nagasaki, Japan
    Contribution
    Methodology
    Competing interests
    No competing interests declared
  7. Ryuta Takase

    Department of Medical Cell Biology, Institute of Molecular Embryology and Genetics, Kumamoto University, Kumamoto, Japan
    Contribution
    Investigation
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-2860-2647
  8. Kensaku Kohrogi

    Department of Medical Cell Biology, Institute of Molecular Embryology and Genetics, Kumamoto University, Kumamoto, Japan
    Contribution
    Investigation
    Competing interests
    No competing interests declared
  9. Yuko Hino

    Department of Medical Cell Biology, Institute of Molecular Embryology and Genetics, Kumamoto University, Kumamoto, Japan
    Contribution
    Investigation, Methodology
    Competing interests
    No competing interests declared
  10. Yusuke Ono

    Department of Muscle Development and Regeneration, Institute of Molecular Embryology and Genetics, Kumamoto University, Kumamoto, Japan
    Contribution
    Supervision, Methodology
    Competing interests
    No competing interests declared
  11. Eiichi Araki

    Department of Metabolic Medicine, Faculty of Life Sciences, Kumamoto University, Kumamoto, Japan
    Contribution
    Supervision
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-4064-7525
  12. Mitsuyoshi Nakao

    Department of Medical Cell Biology, Institute of Molecular Embryology and Genetics, Kumamoto University, Kumamoto, Japan
    Contribution
    Conceptualization, Supervision, Funding acquisition, Investigation, Writing - original draft, Writing - review and editing
    For correspondence
    mnakao@gpo.kumamoto-u.ac.jp
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-2196-8673

Funding

Japan Society for the Promotion of Science (20H04108)

  • Shinjiro Hino

the Inter-University Research Network for High Depth Omics, IMEG, Kumamoto University

  • Shinjiro Hino

Japan Society for the Promotion of Science (21K19513)

  • Shinjiro Hino

Japan Society for the Promotion of Science (21H02686)

  • Mitsuyoshi Nakao

Japan Society for the Promotion of Science (20KK0185)

  • Mitsuyoshi Nakao

Japan Agency for Medical Research and Development (21gk0210029h0101)

  • Shinjiro Hino

Takeda Science Foundation

  • Shinjiro Hino
  • Mitsuyoshi Nakao

Nakatomi Foundation

  • Shinjiro Hino

Center for Metabolic Regulation of Healthy Aging, Kumamoto University

  • Hirotaka Araki

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

This work was supported by the following funding sources: JSPS KAKENHI (Grant Numbers: 20H04108 and 21K19513 [SH], 21H02686 and 20KK0185 [MN]), Takeda Science Foundation (SH and MN), Nakatomi Foundation (SH), the Japan Agency for Medical Research and Development (21gk0210029h0101 [SH]), Center for Metabolic Regulation of Healthy Aging, Kumamoto Univ. (HA), and the program of the Inter-University Research Network for High Depth Omics, IMEG, Kumamoto University (KE and SH). We would like to thank Editage (https://www.editage.com) for English language editing.

Ethics

Animal experiments were approved by the Animal Care and Use Committee of Kumamoto University (Permit No. A2019-098), and conducted in accordance with its guidelines.

Senior Editor

  1. David E James, University of Sydney, Australia

Reviewing Editor

  1. Christopher Cardozo, Icahn School of Medicine at Mount Sinai, United States

Publication history

  1. Preprint posted: October 9, 2022 (view preprint)
  2. Received: November 8, 2022
  3. Accepted: January 12, 2023
  4. Version of Record published: January 25, 2023 (version 1)

Copyright

© 2023, Araki et al.

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 236
    Page views
  • 27
    Downloads
  • 0
    Citations

Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Hirotaka Araki
  2. Shinjiro Hino
  3. Kotaro Anan
  4. Kanji Kuribayashi
  5. Kan Etoh
  6. Daiki Seko
  7. Ryuta Takase
  8. Kensaku Kohrogi
  9. Yuko Hino
  10. Yusuke Ono
  11. Eiichi Araki
  12. Mitsuyoshi Nakao
(2023)
LSD1 defines the fiber type-selective responsiveness to environmental stress in skeletal muscle
eLife 12:e84618.
https://doi.org/10.7554/eLife.84618

Further reading

    1. Cell Biology
    2. Chromosomes and Gene Expression
    Liangyu Zhang, Weston T Stauffer ... Abby F Dernburg
    Research Article

    Meiotic chromosome segregation relies on synapsis and crossover recombination between homologous chromosomes. These processes require multiple steps that are coordinated by the meiotic cell cycle and monitored by surveillance mechanisms. In diverse species, failures in chromosome synapsis can trigger a cell cycle delay and/or lead to apoptosis. How this key step in 'homolog engagement' is sensed and transduced by meiotic cells is unknown. Here we report that in C. elegans, recruitment of the Polo-like kinase PLK-2 to the synaptonemal complex triggers phosphorylation and inactivation of CHK-2, an early meiotic kinase required for pairing, synapsis, and double-strand break induction. Inactivation of CHK-2 terminates double-strand break formation and enables crossover designation and cell cycle progression. These findings illuminate how meiotic cells ensure crossover formation and accurate chromosome segregation.

    1. Cell Biology
    2. Physics of Living Systems
    Christa Ringers, Stephan Bialonski ... Nathalie Jurisch-Yaksi
    Research Article

    Motile cilia are hair-like cell extensions that beat periodically to generate fluid flow along various epithelial tissues within the body. In dense multiciliated carpets, cilia were shown to exhibit a remarkable coordination of their beat in the form of traveling metachronal waves, a phenomenon which supposedly enhances fluid transport. Yet, how cilia coordinate their regular beat in multiciliated epithelia to move fluids remains insufficiently understood, particularly due to lack of rigorous quantification. We combine experiments, novel analysis tools, and theory to address this knowledge gap. To investigate collective dynamics of cilia, we studied zebrafish multiciliated epithelia in the nose and the brain. We focused mainly on the zebrafish nose, due to its conserved properties with other ciliated tissues and its superior accessibility for non-invasive imaging. We revealed that cilia are synchronized only locally and that the size of local synchronization domains increases with the viscosity of the surrounding medium. Even though synchronization is local only, we observed global patterns of traveling metachronal waves across the zebrafish multiciliated epithelium. Intriguingly, these global wave direction patterns are conserved across individual fish, but different for left and right nose, unveiling a chiral asymmetry of metachronal coordination. To understand the implications of synchronization for fluid pumping, we used a computational model of a regular array of cilia. We found that local metachronal synchronization prevents steric collisions, cilia colliding with each other, and improves fluid pumping in dense cilia carpets, but hardly affects the direction of fluid flow. In conclusion, we show that local synchronization together with tissue-scale cilia alignment coincide and generate metachronal wave patterns in multiciliated epithelia, which enhance their physiological function of fluid pumping.