Abstract
Keisuke Hitachi1, Kunihiro Tsuchida1
1Division for Therapies Against Intractable Diseases, Institute for Comprehensive Medical Science (ICMS), Fujita Health University, Toyoake, Aichi 470-1192, Japan
Correspondence to:
Kunihiro Tsuchida, email: [email protected]
Keywords: myostatin, miRNA, callipyge, IG-DMR ncRNA, skeletal muscle
Received: July 20, 2016 Accepted: December 08, 2016 Published: December 15, 2016
ABSTRACT
Myostatin, a member of the transforming growth factor-beta superfamily, is a negative regulator of skeletal muscle growth and development. Myostatin inhibition leads to increased skeletal muscle mass in mammals; hence, myostatin is considered a potential therapeutic target for skeletal muscle wasting. However, downstream molecules of myostatin in the skeletal muscle have not been fully elucidated. Here, we identified the Dlk1-Dio3 locus at the mouse chromosome 12qF1, also called as the callipyge locus in sheep, as a novel downstream target of myostatin. In skeletal muscle of myostatin knockout mice, the expression of mature miRNAs at the Dlk1-Dio3 locus was significantly increased. The increased miRNA levels are caused by the transcriptional activation of the Dlk1-Dio3 locus, because a significant increase in the primary miRNA transcript was observed in myostatin knockout mice. In addition, we found increased expression of coding and non-coding genes (Dlk1, Gtl2, Rtl1/Rtl1as, and Rian) at the Dlk1-Dio3 locus in myostatin-deficient skeletal muscle. Moreover, epigenetic changes, associated with the regulation of the Dlk1-Dio3 locus, were observed in myostatin knockout mice. Taken together, this is the first report demonstrating the role of myostatin in regulating the Dlk1-Dio3 (the callipyge) locus in the skeletal muscle.
INTRODUCTION
Myostatin, a cytokine belonging to the transforming growth factor-beta superfamily, negatively regulates skeletal muscle growth and development through the regulation of anabolic and catabolic pathways in skeletal muscles [1]. Systemic administration of myostatin results in decreased skeletal muscle mass in mice [2], whereas myostatin-null mice exhibit 200–300% increase in skeletal muscle weight [3]. In human, a mutation that may be related to mis-splicing in the myostatin gene was found in a child with muscle hypertrophy [4]. The loss of skeletal muscle mass is a clinically important problem in diseases (cachexia) and aging (sarcopenia). Thus, myostatin is a promising therapeutic target for skeletal muscle wasting due to several diseases and aging; however, its downstream molecules in the skeletal muscle have not been fully identified.
The callipyge (CLPG) phenotype in sheep shows postnatal skeletal muscle hypertrophy of the lower body [5]. A point mutation (A to G transition) on the regulatory element located between Dlk1 and Gtl2 (also called Meg3) genes is associated with the callipyge phenotype [6]. The paternally inherited CLPG mutation can induce skeletal muscle hypertrophy with an increase in Dlk1 and Rtl1 (also known as PEG11) expression [7, 8]. Contrary to the paternally inherited CLPG mutation, the maternally inherited CLPG mutation decreases the expression levels of the paternally expressed genes with an increase in Gtl2, antiRtl1 (Rtl1as), and Rian/Meg8 expression [7, 8]. In addition, the maternally inherited CLPG mutation upregulates the expression of a miRNA cluster containing more than 100 miRNAs, located at the Dlk1-Dio3 locus [9]. Previous reports have indicated the roles of these miRNAs in skeletal muscle regeneration and mitochondrial biogenesis [10, 11]. Although the transcription factor MEF2a was shown to be involved in regulating the expression of this miRNA cluster in skeletal muscle regeneration [10], the precise mechanism that controls the transcription of the Dlk1-Dio3 locus in skeletal muscle remains unclear.
We have previously examined the global miRNA expression profile in the skeletal muscle of myostatin knockout mice by microarray analysis and found 11 miRNAs whose expression was upregulated or downregulated in myostatin-deficient skeletal muscle [12]. We showed that one of these miRNAs, miR-486, is a direct molecular target of myostatin in the regulation of skeletal muscle mass [12]. However, the role of other miRNAs whose expression was altered in myostatin knockout mice remains to be determined. To characterize whether these miRNAs are associated with the muscular phenotype in myostatin knockout mice, here we validated the miRNA expression and found the increased expression of the miRNA cluster at the Dlk1-Dio3 locus in myostatin-deficient skeletal muscle. In addition to miRNAs, myostatin deficiency also increased the expression of a large number of genes located within the Dlk1-Dio3 locus in the skeletal muscle. Moreover, we found epigenetic changes, decreased DNA methylation at the Dlk1-Dio3 locus and increased IG-DMR ncRNA expression, in myostatin knockout mice. Taken together, the Dlk1-Dio3 locus was shown to be a new downstream target of myostatin signaling, suggesting that myostatin deficiency leads to activation of transcription at the Dlk1-Dio3 locus via the regulation of epigenetic modulation.
RESULTS
Increased expression of miRNAs encoded by the Dlk1-Dio3 locus in myostatin-deficient skeletal muscle
We previously identified 11 miRNAs which were upregulated or downregulated in myostatin knockout mice by microarray analysis, and validated the expression of miR-486 and miR-206 [12]. However, the expression of the remaining 9 miRNAs (miR-411, miR-434-3p, miR-299*, miR-193, miR-146b, miR-379, miR-193b, miR-22, and miR-223) has not been verified. To reveal the role of these miRNAs in the muscular phenotype of myostatin knockout mice, we quantified their expression by qRT-PCR, and observed a significant increase in the expression of miR-411, miR-434-3p, miR-379, and miR-193b in myostatin knockout mice (Figure 1A). Among these, the increased expressions of miR-411, miR-434-3p, and miR-379 were of interest because these miRNAs are expressed from the mouse chromosome 12qF1 (chr12qF1), called as the Dlk1-Dio3 locus (Figure 1B). The Dlk1-Dio3 locus is highly conserved in mammals and contains both paternally and maternally expressed genes [13]. In addition, the miRNA cluster (more than 100 miRNAs) is maternally transcribed from the Dlk1-Dio3 locus [14]. We randomly picked 9 miRNAs (miR-337, miR-540-3p, miR-127, miR-434-5p, miR-329, miR-543-3p, miR-376a, miR-300, and miR-381) expressed from the Dlk1-Dio3 locus and validated their expression by qRT-PCR. The results of this quantification showed that the expression of all miRNAs tested was significantly increased in myostatin knockout mice (Figure 1C). These results indicate that myostatin deficiency increases the expression of the miRNA cluster located in the Dlk1-Dio3 locus in the skeletal muscle.
Figure 1: The expression of miRNAs derived from the Dlk1-Dio3 locus is upregulated in myostatin-deficient skeletal muscle. (A) A significant increase in miR-411, miR-434-3p, miR-379, and miR-193b expression in myostatin-deficient skeletal muscle at 13 weeks of age was determined by quantitative RT-PCR. (B) Schematic diagram of the Dlk1-Dio3 locus structure. This locus contains imprinting genes; Dlk1, Rtl1, and Dio3 (shown as blue boxes) are paternally expressed, whereas Gtl2, Rtl1as, Rian, Mirg (shown as white boxes), and miRNAs (arrowheads) are maternally expressed. Red arrowheads indicate the miRNAs whose expression was validated by qRT-PCR in Figure 1C. (C) Quantitative RT-PCR analysis shows a significant increase in the expression of miR-337, miR-540-3p, miR-127, miR-434-5p, miR-329, miR-543-3p, miR-376a, miR-300, and miR-381 in myostatin-deficient skeletal muscle at 13 weeks of age. The wild-type mice are the same age as the myostatin knockout mice. The results were normalized to U6 small RNA expression. Data are the mean ± SD (n = 5). *P < 0.05. **P < 0.01. ***P < 0.001.
Myostatin deficiency leads to remarkable skeletal muscle hypertrophy in mice (Figure 2A). To clarify the role of these miRNAs in skeletal muscle hypertrophy induced by myostatin deficiency, we transfected miRNA mimics (miR-127, miR-300, miR-329, miR-337-3p, miR-376a, miR-379, miR-381, miR-411, miR-434-5p, and miR-540-3p) into C2C12 myotubes and measured the myotube diameter. Immunofluorescence staining of myotubes with an antibody specific to the myosin heavy chain, which is a terminally differentiated marker of skeletal muscle cells, demonstrated that overexpression of miR-411 and miR-540-3p significantly increased the C2C12 myotube diameter (miR-411; 100 ± 4.89 vs. 112 ± 2.94, P = 0.022, and miR-540-3p; 100 ± 4.89 vs. 113 ± 3.02, P = 0.017) (Figure 2B and 2C). Overexpression of miR-434-5p tended to increase the myotube diameter (100 ± 4.89 vs. 108 ± 4.23, P = 0.091), whereas overexpression of the remaining 7 miRNAs (miR-127, miR-300, miR-329, miR-337-3p, miR-376a, miR-379, and miR-381) showed no significant effect on myotube diameter (Figure 2D).
Figure 2: miR-411 and miR-540-3p increase the C2C12 myotube diameters. (A) Myostatin knockout mice show a significant increase in the skeletal muscle mass. (B) Confluent C2C12 cells were differentiated into myotubes and transfected with 25 nM of miRNA mimics. Myotubes were transfected with miR-411, miR-434-5p, and miR-540-3p mimics and immunostained with an anti-myosin heavy chain antibody and DAPI. (C and D) Myotube diameter was measured as described in the materials and methods section and was shown as the percent of the control. Transfection of miR-411 and miR-540-3p mimics, but not other miRNAs, significantly increased the myotube diameter. Bar, 250 μm. Error bars indicate the SD of 3 independent experiments. *P < 0.05.
A previous study has revealed that the miRNA cluster at the Dlk1-Dio3 locus targets SFRP genes, which encode secreted inhibitors of the Wnt signaling pathway, in murine skeletal muscle [10]. In addition, myostatin was shown to negatively regulate skeletal muscle growth and development through the downregulation of Wnt signaling pathway [15–17]. We therefore evaluated the expression levels of SFRP genes (Sfrp1, Sfrp2, Sfrp4, and Sfrp5) in myostatin-deficient skeletal muscle by qRT-PCR. A significant decrease in Sfrp1 expression was observed in both 5- and 13-week-old myostatin knockout mice (Figure 3A and 3B). At 5 weeks of age, the expression level of Sfrp5 was also significantly reduced in myostatin knockout mice. Taken together, these results suggest that increased expression of the miRNA cluster on chr12qF1 would contribute to decrease the expression of Sfrp family genes in myostatin knockout mice. This could explain one of the mechanisms for the muscular phenotype in myostatin knockout mice.
Figure 3: Decreased expression of SFRP family genes in myostatin-deficient skeletal muscle. Quantification of expression levels of SFRP (Secreted frizzled-related protein) genes, which are secreted inhibitors of Wnt signaling pathway, was performed by qRT-PCR. A significant decrease in the expression of Sfrp1 was observed in myostatin-deficient skeletal muscle at 5 (A) and 13 (B) weeks of age. In addition, a significant decrease in Sfrp5 expression was also detected in myostatin knockout mice at 5 weeks of age (A). The wild-type mice are the same age as the myostatin knockout mice. The results were normalized to the expression of Rpl26. Data are mean ± SD (n = 5). *P < 0.05, ***P < 0.001.
Myostatin deficiency increased transcription at the Dlk1-Dio3 locus in skeletal muscle
To determine whether myostatin deficiency activates the transcription of chr12qF1 miRNAs, we further quantified the expression levels of the primary transcripts of miR-127 (pri-miR-127) and miR-411 (pri-miR-411) by qRT-PCR. The quantification results showed a significant increase in pri-miR-127 expression in myostatin knockout mice (Figure 4A). The expression level of pri-miR-411 also tended to increase (Figure 4A). Besides miRNAs, the Dlk1-Dio3 locus includes paternally expressed genes (Dlk1, Rtl1, and Dio3) and maternally expressed genes (Gtl2, Rtl1as, Rian, and Mirg). This prompted us to quantify the expression levels of these chr12qF1 genes in myostatin-deficient skeletal muscle. Quantitative RT-PCR analysis showed that the expression levels of Dlk1, Gtl2, and Rian, but not Mirg, were significantly increased in myostatin knockout mice (Figure 4B). Strand-specific semi-quantitative RT-PCR analysis also showed a significant increase in Rtl1 and Rtl1as expression (Figure 4C and 4D). These results indicate that myostatin deficiency increased the expression of chr12qF1 transcripts including miRNAs, regardless of the paternal or maternal chromosome.
Figure 4: Myostatin deficiency increases transcription at the Dlk1-Dio3 locus. (A) Quantitative RT-PCR analysis shows a significant increase in the expression of the primary transcript of miR-127 in myostatin-deficient skeletal muscle at 13 weeks of age. (B) Myostatin deficiency increases the expression levels of the maternally expressed Gtl2 and Rian, as well as the paternally expressed Dlk1. The results were normalized to Rpl26 expression. Data are the mean ± SD (n = 5). *P < 0.05. (C) Strand-specific RT-PCR shows increased expression of Rtl1 and Rtl1as in myostatin-deficient skeletal muscle. (D) Quantitative analyses of Rtl1 and Rtl1as expression using Image J software. The wild-type mice are the same age as the myostatin knockout mice. Data are the mean ± SD (n = 3). **P < 0.01, ***P < 0.001.
We have previously shown that skeletal muscle hypertrophy in myostatin knockout mice is initially observed at 5 weeks of age [12]. Consequently, we next evaluated the expression levels of chr12qF1 transcripts at 5, 9, and 13 weeks of age by qRT-PCR. Although the expression of chr12qF1 transcripts was higher in myostatin knockout mice compared to that in wild-type mice, the expression level of these transcripts was drastically decreased with increasing age (Figure 5A–5D). The same trend was observed for another imprinting gene, Igf2, but not in the non-imprinting miRNA, miR-486 (Figure 5E and 5F), suggesting that the expression of imprinting genes is gradually decreased in the skeletal muscle, as mice get older.
Figure 5: Mice growth decreases the expression levels of imprinting genes in the skeletal muscle. Quantitative RT-PCR showing the decreased expression of miR-411 (A) miR-540-3p (B) pri-miR-411 (C) and Gtl2 (D) in both wild-type and myostatin knockout mice with age. (E) The expression level of Igf2, the growth-related imprinting gene, is also decreased in the skeletal muscle as mice grow. (F) The expression of miR-486 that is increased in myostatin knockout mice is not changed with age. The expression level of either U6 or Rpl26 was used to normalize the qRT-PCR results. Data are the mean ± SD (n = 5). *P < 0.05, **P < 0.01, ***P < 0.001.
Altered DNA methylation status and non-coding RNA expression in myostatin-deficient skeletal muscle
The DNA methylation levels of two differentially methylated regions (IG-DMR and Gtl2-DMR) located between the Dlk1 and Gtl2 genes affect the expression of chr12qF1 transcripts [18–20]. We therefore examined whether the increased expression of genes at the Dlk1-Dio3 locus in myostatin knockout mice depends on the changes in DNA methylation status using bisulfite sequencing. Although the DNA methylation status of IG-DMR was not changed (Figure 6A and 6B), that of Gtl2-DMR was found to be decreased in myostatin knockout mice (Figure 6C and 6D). We further investigated whether myostatin deficiency affects the expression of de novo DNA methyltransferases (Dnmt3a1 and Dnmt3a2) in the skeletal muscle by qRT-PCR and found that the expression of Dnmt3a2 but not Dnmt3a1 was significantly decreased in myostatin knockout mice (Figure 6E).
Figure 6: Methylation status of the IG-DMR and Gtl2-DMR domain in myostatin-deficient skeletal muscle. Bisulfite sequencing analysis shows the methylation maps around the IG-DMR (A) and Gtl2-DMR (C) domain. The CpG sites are indicated by circles (black circles indicate methylated and white circles indicate unmethylated cytosine sites, respectively) and each row represents an individual clone. The methylation ratio of the IG-DMR (B) and Gtl2-DMR (D) domain in the myostatin knockout mice is shown as the percentage of that in wild-type mice. (E) Quantitative RT-PCR showing decreased Dnmt3a2 expression in myostatin-deficient skeletal muscle. The wild-type mice are the same age as the myostatin knockout mice. The results were normalized to the expression of Rpl26. Data are the mean ± SD (n = 5). ***P < 0.001.
Recently, IG-DMR and IPW RNAs were identified to regulate the expression of almost all the transcripts at the Dlk1-Dio3 locus [21, 22]. IG-DMR RNA is an enhancer RNA-like ncRNA expressed from the IG-DMR domain at the Dlk1-Dio3 locus. IPW RNA is a paternally expressed long ncRNA (lncRNA) located in the Prader-Willi syndrome-associated locus. Notably, loss-of-function of IG-DMR ncRNA and gain-of-function of IPW lncRNA resulted in a drastic decrease in chr12qF1 transcripts in mouse ESCs and human iPSCs, respectively. We thus examined the expressions of these ncRNAs in myostatin knockout mice. Intriguingly, qRT-PCR analysis showed a drastic increase in IG-DMR ncRNA expression in myostatin knockout mice (Figure 7A), whereas the expression level of IPW lncRNA was not changed (Figure 7B). From these results, we concluded that myostatin deficiency affected the expression of epigenetic modulators including DNA methyltransferase and ncRNA in the skeletal muscle.
Figure 7: Increased expression of IG-DMR ncRNA in myostatin-deficient skeletal muscle. (A) Quantitative RT-PCR analysis shows the increased expression of IG-DMR ncRNA in myostatin-deficient skeletal muscle at 13 weeks of age. (B) Myostatin deficiency does not affect the expression level of IPW lncRNA. The wild-type mice are the same age as the myostatin knockout mice. The results were normalized to the expression of Rpl26. Data are the mean ± SD (n = 5). **P < 0.01.
DISCUSSION
Myostatin is a powerful negative regulator of skeletal muscle growth and development. Although myostatin inhibition in the skeletal muscle is valuable for medical and agricultural applications, its downstream molecules in the skeletal muscle have not been fully identified. In this study, we found higher expression of mature miRNAs derived from the mouse chromosome 12qF1, called as the Dlk1-Dio3 locus, in the tibialis anterior muscle of myostatin knockout mice compared to that of wild-type mice. Our results are consistent with a recent study by Wu et al., wherein the expression of miR-431 derived from the Dlk1-Dio3 locus is increased in myostatin knockout mice [23]. In addition to the miRNAs, we also found increased expression levels of both coding and non-coding genes at the Dlk1-Dio3 locus, regardless of paternal or maternal origin, in myostatin-deficient skeletal muscle. These observations indicate that myostatin deficiency leads to activate the transcription at the Dlk1-Dio3 locus. In fact, we showed that an increase in mature miR-127 expression is the result of activated transcription of the primary miR-127 transcript in myostatin knockout mice. Notably, the genes at the Dlk1-Dio3 locus are transcribed as one long polycistronic RNA in the skeletal muscle [24]. Taken together, our study is the first to show that myostatin deficiency increases the global gene expression at the Dlk1-Dio3 locus in murine skeletal muscle.
In sheep, the Dlk1-Dio3 locus is associated with the double-muscle phenotype called callipyge, which shows a 30–40% postnatal increase in the mass of hindlimb muscles [5]. A paternally inherited callipyge mutation is responsible for an increase in Dlk1 and Rtl1 expression [6]. Transgenic overexpression of Dlk1 or Rtl1 in mice resulted in increased skeletal muscle mass [25, 26], whereas Dlk1 ablation resulted in reduced body weight and skeletal muscle mass [27], indicating that Dlk1 and Rtl1 are responsible for muscular hypertrophy phenotype in the callipyge sheep. Intriguingly, bi-allelic mutation of this locus is not characterized by muscle hypertrophy, because increased expression of maternally expressed ncRNAs might inhibit the Dlk1 expression [25]. Most recently, Gao et al. showed that deletion of the maternally expressed miR-379/miR-544 cluster located between the Rian and Mirg genes resulted in callipyge-like muscular hypertrophy with increased Dlk1 expression in mice [28]. Thus, increased expression of miRNAs at this locus is likely to attenuate muscular hypertrophy phenotype caused by Dlk1 overexpression. We did find a small but a significant increase in Dlk1 expression and a remarkable increase in Rtl1 expression in the myostatin knockout mice. Concomitantly, a significant increase in maternally expressed miRNAs at the Dlk1-Dio3 locus was also observed. Although we demonstrated the potency of the Dlk1-Dio3 locus derived miRNAs (miR-411 and miR-540-3p) to induce myotube hypertrophy using an in vitro model, muscular hypertrophy phenotype observed in the callipyge mutation and myostatin deficiency would be caused by slightly different mechanism. In fact, a remarkable increase in skeletal muscle mass in myostatin knockout mice was observed after 5 weeks of age [12], whereas muscle hypertrophy was not observed in mice with the deletion of the miR-379/miR-544 cluster after 3 weeks of age [28]. In addition to the Dlk1-Dio3 locus, myostatin signaling was shown to regulate the expression of miR-486 [12], miR-29b/c [29], and E3 ubiquitin ligase Atrogin-1/MAFbx [30] in skeletal muscle. Thus, in spite of an increase in miRNAs on the Dlk1-Dio3 locus, the increased expression of Dlk1 and Rtl1 might be sufficient to induce the muscular hypertrophy in myostatin knockout mice. Alternatively, increased expression of miRNAs at the Dlk1-Dio3 locus is likely to contribute to the promotion of muscle regeneration rather than the skeletal muscle hypertrophy in myostatin knockout mice. Snyder et al. showed that the expression of miRNAs from the Dlk1-Dio3 locus is downregulated in Mef2a knockout mice and a shortage of these miRNAs results in upregulation of the SFRP genes that impair skeletal muscle regeneration by attenuating Wnt signaling [10]. Consistent with a previous report [15], we observed significantly decreased expression of SFRP genes in myostatin-deficient skeletal muscle in parallel with increased expression of miRNAs at the Dlk1-Dio3 locus (Figures 1 and 3). Considering the fact that myostatin inhibition facilitates muscle regeneration [31], increased expression of miRNAs at the Dlk1-Dio3 locus in myostatin-deficient skeletal muscle might contribute to the promotion of muscle regeneration in addition to inducing muscular hypertrophy.
In this study, myostatin deficiency was shown to affect the expression of epigenetic modulators; Dnmt3a2 and IG-DMR ncRNA. Dnmt3a2 is one of the de novo DNA methyltransferases in humans and mice [32], and reduced Dnmt3a2 expression causes hypomethylation of imprinted genes in mouse ESCs [33]. Transcription from the Dlk1-Dio3 locus is controlled by Gtl2-DMR that is hypomethylated in myostatin knockout mice. Thus, increased expression of the Dlk1-Dio3 locus in myostatin-deficient skeletal muscle would be a result of a decrease in Dnmt3a2 expression. It is also worthwhile to note that we used purchased wild-type mice as a control, therefore, it is possible that the differences in diet and environment would contribute to the different DNA methylation status at the Gtl2-DMR. In addition to Dnmt3a2, we observed greater than 10-fold increase in IG-DMR ncRNA expression in myostatin-deficient skeletal muscle. IG-DMR ncRNA is produced from the IG-DMR domain at the Dlk1-Dio3 locus and the knockdown of the IG-DMR ncRNA drastically reduces the expression of both paternally and maternally expressed genes at the Dlk1-Dio3 locus in ESCs [21]. Although whether IG-DMR ncRNA is required for activation of the Dlk1-Dio3 locus in the skeletal muscle remains to be determined, it seems likely that increased expression of IG-DMR ncRNA would be one of the factors that activates the expression at the Dlk1-Dio3 locus in myostatin knockout mice.
Clinical trials of myostatin blockade against muscular dystrophy and sarcopenia are now ongoing [34]. Thus, it is important to elucidate the role of the Dlk1-Dio3 locus in the muscle hypertrophy induced by myostatin blockade in postnatal and adult skeletal muscle.
MATERIALS AND METHODS
Animals
All animal experiments were conducted with approval from the Institutional Animal Care and Use Committee of Fujita Health University. Myostatin knockout mice previously described [3] were a kind gift from Dr. S.-J. Lee. As a control, age-matched C57BL/6 mice purchased from the Japan SLC, Inc. were used.
Cell culture, immunofluorescence, and quantification of myotube diameters
Mouse C2C12 myoblasts were grown in Dulbecco’s modified Eagle medium (DMEM) supplemented with 10% fetal bovine serum at 37°C and 5% CO2. Confluent C2C12 myoblasts were differentiated into myotubes using DMEM supplemented with 2% horse serum. Myotubes were transfected with 25 nM of miRNA mimics (miCENTURY OX miNaturals, Cosmo Bio) using Lipofectamine 2000 (Thermo Fischer Scientific) at 5 days after induction of differentiation. As a control, 25 nM of miRNA control (Non-target RNA, Cosmo Bio) was transfected. Two days after transfection, the myotubes were fixed and stained with an anti-myosin heavy chain antibody (MF20, Developmental Studies Hybridoma Bank) and DAPI. Cell images were visualized and myotube diameters were measured as described previously [12].
RNA isolation, reverse transcriptase reaction, and quantitative-PCR
RNA samples were prepared from skeletal muscle as described previously [12]. In brief, the tibialis anterior muscles were surgically dissected from 5-, 9-, and 13-week-old mice under anesthesia and were quickly frozen in liquid nitrogen. Total RNAs were prepared using a miRNeasy Mini kit (Qiagen) according to the manufacturer’s instructions.
The reverse transcriptase reaction was performed using a PrimeScript RT-PCR kit (TaKaRa) and a miScript Reverse Transcription Kit (Qiagen) for the detection of mRNAs and mature miRNAs, respectively. A SuperScript III First-Strand Synthesis System with Random primers (Thermo Fisher Scientific) was used for the detection of IG-DMR and IPW ncRNA according to the manufacturer’s instructions.
The expression of mRNAs and ncRNAs was estimated by quantitative real-time PCR using a SYBR Premix Ex Taq (TaKaRa), and that of mature miRNAs were quantified by a miScript SYBR Green PCR Kit (Qiagen).
Strand-specific RT-PCR
A strand-specific reverse transcriptase reaction was performed for detecting the paternally expressed Rtl1 transcript and the maternally expressed Rtl1as transcript using a SuperScript III First-Strand Synthesis System using the previously reported RT-primer [35], Rtl1 primer 5′- GGAGCCACTTCATGCCTAAGACGA-3′ and Rtl1as primer is 5′-GTGGAGAACTTCGCTGTCATCGC-3′. Semi-quantitative PCR for both strand-derived transcripts was conducted using TaKaRa Ex Taq (TaKaRa). The PCR results were quantified with the Image J 1.49v software, and the results were normalized to Rpl26 expression. The primers used are listed in Supplementary Table S1.
Genomic DNA isolation and bisulfite sequencing
Genomic DNA was isolated from the tibialis anterior muscles of 13-week-old mice using PureLink Genomic DNA Kits (Thermo Fisher Scientific). Genomic DNA (200 ng) was bisulfite treated using a Cells-to-CpG Bisulfite Conversion Kit (Thermo Fisher Scientific) according to the manufacturer’s instructions. The PCR reaction was performed with an EpiTaq HS (TaKaRa) under the following conditions: denaturation at 94°C for 60 s and 35 cycles each of 94°C for 30 s, 55°C for 30 s, and 72°C for 60 s. A second nested PCR reaction was performed using the following conditions: denaturation at 94°C for 60 s and 25 cycles each of 94°C for 30 s, 55°C for 30 s, and 72°C for 60 s. The primers used are listed in Supplementary Table S1. The PCR products were subcloned into the pGEM T-easy vector (Promega), and clones derived from each independent subject (n = 3) were sequenced. The methylation status of the region (higher than 97% sequence identity) was determined and analyzed with QUMA (http://quma.cdb.riken.jp/top/quma_main_j.html).
ACKNOWLEDGMENTS
We thank Drs. H. Ageta, and M. Nakatani, and A. Uezumi for discussion and technical advice.
CONFLICTS OF INTEREST
The authors declare no conflicts of interest.
GRANT SUPPORT
This work was supported in part by JSPS KAKENHI (25860151 and 16K08599), Intramural Research Grant (26-8) for Neurological and Psychiatric Disorders of NCNP, and a Grant-in-Aid from the NAKATOMI Foundation.
REFERENCES
1. Lee SJ. Regulation of muscle mass by myostatin. Annu Rev Cell Dev Biol. 2004; 20:61–86. doi: 10.1146/annurev.cellbio.20.012103.135836.
2. Zimmers TA, Davies MV, Koniaris LG, Haynes P, Esquela AF, Tomkinson KN, McPherron AC, Wolfman NM, Lee SJ. Induction of cachexia in mice by systemically administered myostatin. Science. 2002; 296:1486–8. doi: 10.1126/science.1069525.
3. McPherron AC, Lawler AM, Lee SJ. Regulation of skeletal muscle mass in mice by a new TGF-beta superfamily member. Nature. 1997; 387:83–90. doi: 10.1038/387083a0.
4. Schuelke M, Wagner KR, Stolz LE, Hübner C, Riebel T, Kömen W, Braun T, Tobin JF, Lee SJ. Myostatin mutation associated with gross muscle hypertrophy in a child. N Engl J Med. 2004; 350:2682–8. doi: 10.1056/NEJMoa040933.
5. Koohmaraie M, Shackelford SD, Wheeler TL, Lonergan SM, Doumit ME. A muscle hypertrophy condition in lamb (callipyge): characterization of effects on muscle growth and meat quality traits. J Anim Sci. 1995; 73:3596–607. doi: 10.2527/1995.73123596x.
6. Freking BA, Murphy SK, Wylie AA, Rhodes SJ, Keele JW, Leymaster KA, Jirtle RL, Smith TPL. Identification of the single base change causing the callipyge muscle hypertrophy phenotype, the only known example of polar overdominance in mammals. Genome Res. 2002; 12:1496–506. doi: 10.1101/gr.571002.
7. Charlier C, Segers K, Karim L, Shay T, Gyapay G, Cockett N, Georges M. The callipyge mutation enhances the expression of coregulated imprinted genes in cis without affecting their imprinting status. Nat Genet. 2001; 27:367–9. doi: 10.1038/86856.
8. Bidwell CA, Kramer LN, Perkins AC, Hadfield TS, Moody DE, Cockett NE. Expression of PEG11 and PEG11AS transcripts in normal and callipyge sheep. BMC Biol. 2004; 2:17. doi: 10.1186/1741-7007-2-17.
9. Davis E, Caiment F, Tordoir X, Cavaillé J, Ferguson-Smith A, Cockett N, Georges M, Charlier C. RNAi-mediated allelic trans-interaction at the imprinted Rtl1/Peg11 locus. Curr Biol. 2005; 15:743–9. doi: 10.1016/j.cub.2005.02.060.
10. Snyder CM, Rice AL, Estrella NL, Held A, Kandarian SC, Naya FJ. MEF2A regulates the Gtl2-Dio3 microRNA mega-cluster to modulate WNT signaling in skeletal muscle regeneration. Development. 2013; 140:31–42. doi: 10.1242/dev.081851.
11. Yamamoto H, Morino K, Nishio Y, Ugi S, Yoshizaki T, Kashiwagi A, Maegawa H. MicroRNA-494 regulates mitochondrial biogenesis in skeletal muscle through mitochondrial transcription factor A and Forkhead box j3. Am J Physiol Endocrinol Metab. 2012; 303: E1419–27. doi: 10.1152/ajpendo.00097.2012.
12. Hitachi K, Nakatani M, Tsuchida K. Myostatin signaling regulates Akt activity via the regulation of miR-486 expression. Int J Biochem Cell Biol. 2014; 47:93–103. doi: 10.1016/j.biocel.2013.12.003.
13. Paulsen M, Takada S, Youngson NA, Benchaib M, Charlier C, Segers K, Georges M, Ferguson-Smith AC. Comparative sequence analysis of the imprinted Dlk1-Gtl2 locus in three mammalian species reveals highly conserved genomic elements and refines comparison with the Igf2-H19 region. Genome Res. 2001; 11:2085–94. doi: 10.1101/gr.206901.
14. Seitz H, Youngson N, Lin SP, Dalbert S, Paulsen M, Bachellerie JP, Ferguson-Smith AC, Cavaillé J. Imprinted microRNA genes transcribed antisense to a reciprocally imprinted retrotransposon-like gene. Nat Genet. 2003; 34:261–2. doi: 10.1038/ng1171.
15. Steelman CA, Recknor JC, Nettleton D, Reecy JM. Transcriptional profiling of myostatin-knockout mice implicates Wnt signaling in postnatal skeletal muscle growth and hypertrophy. FASEB J. 2006; 20:580–2. doi: 10.1096/fj.05-5125fje.
16. Takata H, Terada K, Oka H, Sunada Y, Moriguchi T, Nohno T. Involvement of Wnt4 signaling during myogenic proliferation and differentiation of skeletal muscle. Dev Dyn. 2007; 236:2800–7. doi: 10.1002/dvdy.21327.
17. Bernardi H, Gay S, Fedon Y, Vernus B, Bonnieu A, Bacou F. Wnt4 activates the canonical β-catenin pathway and regulates negatively myostatin: functional implication in myogenesis. Am J Physiol, Cell Physiol. 2011; 300: C1122–38. doi: 10.1152/ajpcell.00214.2010.
18. Lin SP, Youngson N, Takada S, Seitz H, Reik W, Paulsen M, Cavaillé J, Ferguson-Smith AC. Asymmetric regulation of imprinting on the maternal and paternal chromosomes at the Dlk1-Gtl2 imprinted cluster on mouse chromosome 12. Nat Genet. 2003; 35:97–102. doi: 10.1038/ng1233.
19. Takahashi N, Okamoto A, Kobayashi R, Shirai M, Obata Y, Ogawa H, Sotomaru Y, Kono T. Deletion of Gtl2, imprinted non-coding RNA, with its differentially methylated region induces lethal parent-origin-dependent defects in mice. Hum Mol Genet. 2009; 18:1879–88. doi: 10.1093/hmg/ddp108.
20. Zhou Y, Cheunsuchon P, Nakayama Y, Lawlor MW, Zhong Y, Rice KA, Zhang L, Zhang X, Gordon FE, Lidov HGW, Bronson RT, Klibanski A. Activation of paternally expressed genes and perinatal death caused by deletion of the Gtl2 gene. Development. 2010; 137:2643–52. doi: 10.1242/dev.045724.
21. Kota SK, Llères D, Bouschet T, Hirasawa R, Marchand A, Begon-Pescia C, Sanli I, Arnaud P, Journot L, Girardot M, Feil R. ICR noncoding RNA expression controls imprinting and DNA replication at the Dlk1-Dio3 domain. Dev Cell. 2014; 31:19–33. doi: 10.1016/j.devcel.2014.08.009.
22. Stelzer Y, Sagi I, Yanuka O, Eiges R, Benvenisty N. The noncoding RNA IPW regulates the imprinted DLK1-DIO3 locus in an induced pluripotent stem cell model of Prader-Willi syndrome. Nat Genet. 2014; 46:551–7. doi: 10.1038/ng.2968.
23. Wu R, Li H, Li T, Zhang Y, Zhu D. Myostatin regulates miR-431 expression via the Ras-Mek-Erk signaling pathway. Biochem Biophys Res Commun. 2015; 461:224–9. doi: 10.1016/j.bbrc.2015.03.150.
24. Tierling S, Dalbert S, Schoppenhorst S, Tsai CE, Oliger S, Ferguson-Smith AC, Paulsen M, Walter J. High-resolution map and imprinting analysis of the Gtl2-Dnchc1 domain on mouse chromosome 12. Genomics. 2006; 87:225–35. doi: 10.1016/j.ygeno.2005.09.018.
25. Davis E, Jensen CH, Schroder HD, Farnir F, Shay-Hadfield T, Kliem A, Cockett N, Georges M, Charlier C. Ectopic expression of DLK1 protein in skeletal muscle of padumnal heterozygotes causes the callipyge phenotype. Curr Biol. 2004; 14:1858–62. doi: 10.1016/j.cub.2004.09.079.
26. Xu X, Ectors F, Davis EE, Pirottin D, Cheng H, Farnir F, Hadfield T, Cockett N, Charlier C, Georges M, Takeda H. Ectopic expression of retrotransposon-derived PEG11/RTL1 contributes to the callipyge muscular hypertrophy. PLoS ONE. 2015; 10:e0140594. doi: 10.1371/journal.pone.0140594.
27. Waddell JN, Zhang P, Wen Y, Gupta SK, Yevtodiyenko A, Schmidt JV, Bidwell CA, Kumar A, Kuang S. Dlk1 is necessary for proper skeletal muscle development and regeneration. PLoS ONE. 2010; 5:e15055. doi: 10.1371/journal.pone.0015055.
28. Gao YQ, Chen X, Wang P, Lu L, Zhao W, Chen C, Chen CP, Tao T, Sun J, Zheng YY, Du J, Li CJ, Gan ZJ, et al. Regulation of DLK1 by the maternally expressed miR-379/miR-544 cluster may underlie callipyge polar overdominance inheritance. Proc Natl Acad Sci USA. 2015; 112:13627–32. doi: 10.1073/pnas.1511448112.
29. Goodman CA, McNally RM, Hoffmann FM, Hornberger TA. Smad3 induces atrogin-1, inhibits mTOR and protein synthesis, and promotes muscle atrophy in vivo. Mol Endocrinol. 2013; 27:1946–57. doi: 10.1210/me.2013-1194.
30. Sartori R, Milan G, Patron M, Mammucari C, Blaauw B, Abraham R, Sandri M. Smad2 and 3 transcription factors control muscle mass in adulthood. Am J Physiol, Cell Physiol. 2009; 296:C1248–57. doi: 10.1152/ajpcell.00104.2009.
31. Yaden BC, Croy JE, Wang Y, Wilson JM, Datta-Mannan A, Shetler P, Milner A, Bryant HU, Andrews J, Dai G, Krishnan V. Follistatin: a novel therapeutic for the improvement of muscle regeneration. J Pharmacol Exp Ther. 2014; 349:355–71. doi: 10.1124/jpet.113.211169.
32. Weisenberger DJ, Velicescu M, Preciado-Lopez MA, Gonzales FA, Tsai YC, Liang G, Jones PA. Identification and characterization of alternatively spliced variants of DNA methyltransferase 3a in mammalian cells. Gene. 2002; 298:91–9. doi: 10.1016/S0378-1119(02)00976-9.
33. Popkie AP, Zeidner LC, Albrecht AM, D’Ippolito A, Eckardt S, Newsom DE, Groden J, Doble BW, Aronow B, McLaughlin KJ, White P, Phiel CJ. Phosphatidylinositol 3-kinase (PI3K) signaling via glycogen synthase kinase-3 (Gsk-3) regulates DNA methylation of imprinted loci. Journal of Biological Chemistry. 2010; 285:41337–47. doi: 10.1074/jbc.M110.170704.
34. Garber K. No longer going to waste. Nat Biotechnol. 2016; 34:458–61. doi: 10.1038/nbt.3557.
35. Riordan JD, Keng VW, Tschida BR, Scheetz TE, Bell JB, Podetz-Pedersen KM, Moser CD, Copeland NG, Jenkins NA, Roberts LR, Largaespada DA, Dupuy AJ. Identification of rtl1, a retrotransposon-derived imprinted gene, as a novel driver of hepatocarcinogenesis. PLoS Genetics. 2013; 9:e1003441. doi: 10.1371/journal.pgen.1003441.
