Molecular Medicine Reports Molecular Medicine Reports 1791-2997 1791-3004 D.A. Spandidos 10.3892/mmr.2019.10661 mmr-20-05-4175 Articles Expression patterns of regulatory lncRNAs and miRNAs in muscular atrophy models induced by starvation in vitro and in vivo LeiSi 1 * SheYanling 1 * ZengJie 2 ChenRui 1 ZhouShanyao 1 ShiHuacai 1 Guangdong Traditional Medical and Sports Injury Rehabilitation Research Institute, Guangdong Second Provincial General Hospital, Guangzhou, Guangdong 510317, P.R. China Department of Medical Ultrasonics, The Third Affiliated Hospital, Sun Yat-Sen University, Guangzhou, Guangdong 510630, P.R. China Correspondence to: Dr Rui Chen, Guangdong Traditional Medical and Sports Injury Rehabilitation Research Institute, Guangdong Second Provincial General Hospital, 466 Xin Gang Zhong Road, Guangzhou, Guangdong 510317, P.R. China, E-mail: rui.c.med@163.com

Contributed equally

 112019 10092019 20 5 4175 4185 13032019 30072019 Copyright: © Lei et al. 2019 This is an open access article distributed under the terms of the Creative Commons Attribution-NonCommercial-NoDerivs License, which permits use and distribution in any medium, provided the original work is properly cited, the use is non-commercial and no modifications or adaptations are made.

Starvation or severe deprivation of nutrients, which is commonly seen in surgical patients, can result in catabolic changes in skeletal muscles, such as muscle atrophy. Therefore, it is important to elucidate the underlying molecular regulatory mechanisms during skeletal muscle atrophy. In the present study, muscular atrophy was induced by starvation and the results demonstrated that myosin heavy chain was decreased, whereas muscle RING finger protein 1 and atrogin-1 were increased, both in vitro and in vivo. The impact of starvation on the expression patterns of long non-coding RNAs (lncRNAs) and microRNAs (miRNAs) was next determined. The expression patterns of miR-23a, miR-206 and miR-27b in the starved mice exhibited similar trends as those in starved C2C12 cells in vitro, whereas the expression patterns of six other miRNAs (miR-18a, miR-133a, miR-133b, miR-186, miR-1a and miR-29b) differed between the in vivo and the in vitro starvation models. The present study indicated that in vitro expression of the selected miRNAs was not completely consistent with that in vivo. By contrast, lncRNAs showed excellent consistency in their expression patterns in both the in vitro and in vivo starvation models; six of the lncRNAs (Atrolnc-1, long intergenic non-protein coding RNA of muscle differentiation 1, Myolinc, lncRNA myogenic differentiation 1, Dum and muscle anabolic regulator 1) were significantly elevated in starved tissues and cells, while lnc-mg was significantly decreased, compared with the control groups. Thus, lncRNAs involved in muscle atrophy have the potential to be developed as diagnostic tools.

 long non-coding RNAs microRNAs atrophy starvation skeletal muscle Introduction

Starvation or severe deprivation of nutrients, which is commonly seen in surgical patients, can result in various adaptive and catabolic changes in skeletal muscles, such as muscle atrophy (1,2). Except starvation and fasting, multiple physiological and pathological conditions, including immobilization, mechanical unloading, cancer and chronic disease, can also trigger muscle atrophy (3,4). Our previous studies demonstrated that decreases in actin fiber size resulted in reduced expression of myogenin, and elevated expression of atrogin-1 induced by cobalt dichloride mimicking hypoxia, which subsequently led to skeletal muscle atrophy (5,6).

Muscle atrophy is usually accompanied by a decrease in skeletal muscle mass, an important prognostic indicator of many diseases. Patients with muscular atrophy often have increased morbidity and mortality (7). Decreased muscle protein synthesis accounts predominantly for skeletal muscle atrophy. There are several protective and therapeutic measures against muscular atrophy, such as intake of essential amino acids (8), insulin-like growth factor (IGF)-I treatment (9) and anabolic androgenic steroids (10); however, these measures can cause increased drug resistance and cardiac events, limiting their therapeutic utility for the treatment of muscle atrophy (4,11). Although major advances in our understanding of protein loss in muscle atrophy have been made recently (12–15), it is important to elucidate the underlying molecular regulatory mechanisms of dysfunctional muscle anabolism during skeletal muscle atrophy, in order to develop promising strategies to prevent and treat muscular atrophy.

Long non-coding RNAs (lncRNAs) are a class of non-coding RNAs >200 nucleotides in length that make up the majority of the transcriptome (16,17). lncRNAs have emerged as significant regulators in multiple physiological and pathological processes. lncRNAs have important roles in skeletal muscle differentiation, development and muscular atrophy (3,18–24) through various mechanisms, such as competitive endogenous RNA and cis- and trans-regulatory mechanisms (4,17,18,20,23,25). Long intergenic non-protein coding RNA muscle differentiation 1 (lincMD1) regulates the expression of mastermind like transcriptional coactivator 1 and myocyte enhancer factor 2C in the muscle differentiation program (18). Myogenesis-associated lncRNA (lnc-mg) regulates IGF-2, myosin heavy chain (MYHC), lactamase β and myogenic differentiation factor to control muscle differentiation and development (24,26). MYHC is a crucial component of myosin, which participates in the formation of the cytoskeleton and provides a force for muscular contraction. Atrolnc-1 causes increased expression of muscle RING finger protein 1 (MURF1), leading to myofiber atrophy in mice with chronic kidney disease (3). The abnormal expression of lncRNAs can result in various pathologies, including cancer, and cardiac and muscle diseases (27,28). Although mechanisms involving the obvious functions of lncRNAs continue to be identified, few studies have provided a comprehensive perspective of lncRNAs in terms of their functions in skeletal muscle atrophy.

MicroRNAs (miRNAs, ~22 nucleotides in length) constitute a class of highly conserved, small endogenous non-coding RNA molecules that negatively regulate gene expression at the post-transcriptional level (29–32). Most miRNAs are ubiquitously expressed. Some miRNAs are mainly expressed in muscle but may also be observed at low levels in other tissues. These miRNAs have important roles in muscle differentiation, development and atrophy (33–35), and are predicted to participate in the regulatory networks of myogenesis, muscle fiber type composition, muscle growth and homeostasis (33). The miRNAs associated with muscle differentiation, development and muscular atrophy include miR-133a (36,37), miR-133b (38), miR-206 (39,40), miR-186 (41), miR-23a (42,43), miR-27b (44), miR-29b (32), miR-1a (39) and miR-18a (45,46). miR-1 promotes myoblast differentiation and regeneration (36,37,39,47). miR-29b is required for loss of muscle mass in cases of dexamethasone, tumor necrosis factor-α and H2O2 treatment-induced muscle atrophy (32). Given the crucial roles of miRNAs in muscle proliferation and differentiation, changes in their regulation, and their potential involvement in muscle atrophy during starvation, could be expected. Changes in the expression of muscle development-associated miRNAs in a disease state may induce muscle atrophy, however it is unclear whether these miRNAs are causally involved in adaptive or compensatory responses to muscle atrophy (33).

lncRNAs and miRNAs are both non-coding RNAs, of which miRNAs have been extensively investigated in skeletal muscle development (32,33,35–47). An increasing number of studies has demonstrated that lncRNAs are also involved in skeletal muscle development (3,18,20–23). lncRNAs and miRNAs can mutually restrict muscle development. lncRNA H19 modulates let-7 availability by acting as a molecular sponge. H19 depletion leads to accelerated muscle differentiation, indicating that this lncRNA counteracts muscle differentiation (48). lncRNA Sirt1 antisense fully bound to sirtuin 1 mRNA forms an RNA duplex that promotes myoblast proliferation and inhibits differentiation by competing with miR-34a (49). Few studies have provided a comprehensive perspective of lncRNAs and miRNAs in terms of their regulation of skeletal muscle atrophy. Therefore, the present study aimed to investigate the expression levels of regulatory lncRNAs and miRNAs in mice with muscular atrophy induced by starvation in vitro and in vivo.

In the present study, a serum starvation C2C12 cell model and a starved muscular atrophy mouse model were established, and analyzed by reverse transcription-quantitative PCR (RT-qPCR), western blot analysis, immunofluorescence staining and hematoxylin and eosin (H&E) staining. Subsequently, lncRNAs and miRNAs associated with muscle differentiation or atrophy, whose expression trends in the muscle differentiation or atrophy models were consistent with the findings of previous studies (3,18–24,32,33,35–47), were selected to identify differences of regulatory non-coding RNA expression. The current findings indicated that lncRNAs showed similar expression trends in vitro and in vivo and may have great potential as a diagnostic tool.

 Materials and methods <sec> <title>Cell culture

The mouse myoblast C2C12 cell line was purchased from the Stem Cell Bank of the Type Culture Collection of the Chinese Academy of Sciences and cultured at 37°C in 5% CO2 and high-glucose Dulbecco's modified Eagle's medium (cat. no. 12100046; Gibco; Thermo Fisher Scientific, Inc.) supplemented with 10% (v/v) fetal bovine serum (cat. no. SH30070.03; HyClone; GE Healthcare Life Sciences). For myotube differentiation, C2C12 myoblasts were incubated in 12-well plates with 2% horse serum (cat. no. SH30074.03; HyClone; GE Healthcare Life Sciences) for 72 h, according to previous studies (17,50). For starvation studies, the medium of the differentiated myotubes was replaced with serum-free medium for an additional 48 h (51).

 Mouse strains and starvation-induced muscular atrophy model

Male C57BL/6 mice (age, 6 weeks; weight, 20–22 g) were acquired from the Animal Laboratory of Sun Yat-sen University (Guangzhou, China). Mice were maintained under standard conditions of 22±2°C at a relative humidity of 50–60% and housed under alternate 12 h dark/light cycle conditions. Mice were allowed free access to food and drinking water. All animal procedures were in accordance with the Suggestions for the Care and Use of Laboratory Animals of the Ministry of Science and Technology of the People's Republic of China (52). The study was approved by the Ethics Committee of Guangdong Second Provincial General Hospital (no. 2019-YJSWZ-001).

Following a seven-day acclimation period, the mice were randomly divided into two groups (n=6 mice per group): The starvation/fasting model group and the control group. The starvation/fasting-induced muscular atrophy model group was generated by giving mice access to drinking water for 48 h only, whereas the control group mice had access to food and drinking water. The humane endpoint used in the present study was body weight loss no more than 20%. The mice exhibited normal grooming and the total body weight loss was <20% when the mice were sacrificed. Mice were weighed and sacrificed by cervical dislocation under terminal anesthesia by inhalation of isoflurane, and tibialis anterior muscles were dissected, weighed and placed in 4% neutral formalin fixative solution for 24 h at 25°C, or immediately snap frozen in liquid nitrogen before storing at −80°C.

 RNA extraction and RT-qPCR Isolation of RNA from cultured C2C12 cells

For in vitro analyses, cell samples were washed with PBS before lysis in TRIzol reagent (cat. no. 15596018; Thermo Fisher Scientific, Inc.). RNA extraction was performed according to the manufacturer's instructions.

 Isolation of RNA from tibialis anterior muscles

Frozen tibialis anterior muscles were lysed in TRIzol reagent (Thermo Fisher Scientific, Inc.) and homogenized with a tissue homogenizer (Shanghai Jingxin Industrial Development Co., Ltd.). RNA extraction was performed as described above.

 RT-qPCR

For mRNA and lncRNA detection, total RNA was reverse-transcribed into cDNA using the PrimeScript RT Master mix (cat. no. RR036A; Takara Bio, Inc.). For miRNA detection, total RNA was reverse-transcribed into cDNA using the Mir-X™ miRNA First-Strand Synthesis kit (cat. no. 638313; Takara Bio, Inc.). The abundance of mRNAs, lncRNAs and miRNAs was measured with SYBR-Green Mix (cat. no. RR820A; Takara Bio, Inc.) on a StepOne Plus Real-Time PCR System (Applied Biosystems; Thermo Fisher Scientific, Inc.). PCR for detecting mRNAs and lncRNAs was performed as follows: Denaturation at 95°C for 30 sec, followed by 40 cycles of denaturation at 95°C for 5 sec, annealing at 55°C for 30 sec, and extension at 72°C for 30 s. PCR for detecting miRNAs was performed as follows: Denaturation at 95°C for 10 sec, followed by 40 cycles of denaturation at 95°C for 5 sec, annealing and extension at 60°C for 30 sec. mRNA and lncRNA expression was normalized to that of 18S RNA, and miRNA expression was normalized to that of U6 using the 2−ΔΔCq method (53). A universal reverse primer, mRQ 3′ Primer (cat. no. 638313; Takara Bio, Inc.), was used for the all the miRNAs. Primer sequences are listed in Tables I and II.

 Immunofluorescence staining and myotube diameter measurements

Myotube diameter was analyzed by MYHC staining. Harvested cells were fixed with a 4% paraformaldehyde solution three times for 10 min at 25°C, blocked with goat serum (cat. no. C0265; Beyotime Institute of Biotechnology) for 1 h at 25°C, and then incubated with mouse anti-MYHC primary antibody (1:80; cat. no. MAB4470; R&D Systems, Inc.) overnight at 4°C. After removing the solution containing the primary antibody, cells were incubated with Alexa Fluor 594-conjugated AffiniPure Goat Anti-Mouse IgG (H+L; 1:100; cat. no. AS077; ABclonal Biotech Co., Ltd.) for 1 h. DAPI (cat. no. KGA215; KeyGen Biotech Co., Ltd.) was used to stain nuclei. Images were acquired at magnification, ×400 using a fluorescent Leica DMIL LED microscope (Leica Microsystems GmbH). Fiber diameters were measured in randomly selected fields from three different wells of the control and three different wells of the serum-starved myotubes using ImageJ software (v1.44P; National Institutes of Health). Three diameters were measured per myotube, and 80 myotubes were measured per well.

 H&E staining and determining the cross-sectional area (CSA) of muscle fibers

Muscles were placed in neutral formalin fixative solution for 24 h at 25°C. To assess tissue morphology, tissues were embedded in paraffin, then 10 µm thick transverse sections of tibialis anterior muscles were stained with hematoxylin for 5 min and eosin for 1 min (both at 25°C) and sealed with neutral balsam. After H&E staining, six pairs of samples were placed under a microscope at magnification, ×200 to capture light microscope images (Leica Microsystems GmbH). The obtained images were analyzed by ImageJ software (v1.44P; National Institutes of Health) to measure muscle fiber CSA. The measured data were used to calculate the mean ± standard deviation.

 Western blot analysis

Serum-starved C2C12 myotubes were used to extract total protein. Frozen tibialis anterior muscles were homogenized with a tissue homogenizer (Shanghai Jingxin Industrial Development Co., Ltd.) to extract total protein. RIPA lysis buffer (100 µl; cat. no. P0013D; Beyotime Institute of Biotechnology) was used to extract protein. Equal quantities of protein (20 µg), determined using an enhanced Bicinchoninic Acid Protein Assay Kit (cat. no. P0010S; Beyotime Institute of Biotechnology) were separated via 10–12% SDS-PAGE and transferred onto polyvinylidene fluoride membranes. Membranes were blocked with 5% non-fat milk for 1 h at 25°C and incubated with primary antibodies targeting MURF1 (1:1,000; cat. no. ab77577; Abcam), atrogin-1 (1:2,000; cat. no. ab168372; Abcam), MYHC (1:1,000; cat. no. MAB4470; R&D Systems, Inc.) or tubulin (1:5,000; cat. no. AC021; ABclonal Biotech Co., Ltd.) overnight at 4°C. Tubulin was used as the control for protein quantification. Membranes were then incubated with horseradish peroxide-conjugated goat anti-mouse (1:5,000; cat. no. AS003; ABclonal Biotech Co., Ltd.) or anti-rabbit (1:10,000; cat. no. AS014; ABclonal Biotech Co., Ltd.) secondary antibodies for 1 h at 25°C. Band intensity was determined using a chemiluminescent imaging system (Tanon-V8 Pro; Tanon Sciences and Technology, Co., Ltd.). Densitometry was performed to determine protein expression using TanonImage (version 1.00; Tanon Sciences and Technology, Co., Ltd.)

 Statistical analysis

GraphPad Prism 5.0 (GraphPad Software, Inc.) was used to perform statistical analyses. Student's t-test was applied for the comparison of two groups. P<0.05 was considered to indicate a statistically significant difference.

 Results <sec> <title>Establishment of the serum starvation C2C12 cell model in vitro

Immunofluorescence staining results indicated that, under serum-starvation conditions, the diameter of starved C2C12 myotubes was thinner compared with control (regularly cultured) myotubes (Fig. 1A and B). The protein expression levels of MYHC were significantly decreased in starved myotubes compared with control myotubes (Fig. 1C). An increase in the transcript and protein expression levels of MURF1 and atrogin-1 were also detected in starved myotubes, by RT-qPCR and western blot analysis respectively (Fig. 1D and E). These results suggested that starvation-induced muscle atrophy was successfully established in C2C12 myotubes.

 lncRNA expression patterns in atrophic C2C12 cells

The levels of various lncRNAs associated with muscle development were detected in order to identify differences in regulatory lncRNA expression between control and serum-starved C2C12 cells by RT-qPCR. lncRNAs whose altered expression in the starvation model was consistent with previous studies were selected for subsequent investigation (3,18–24). It was revealed that seven lncRNAs matched these criteria. Six of these lncRNAs [Atrolnc-1, lincMD1, Myolinc, lncMyoD, Dum and muscle anabolic regulator 1 (MAR1)] exhibited significantly increased levels in starved C2C12 cells compared with control C2C12 cells (Fig. 2A-F). By contrast, significantly lower lnc-mg levels were observed in starved C2C12 myotubes (Fig. 2G). In the differentiated starved model, the expression levels of lincMD1, Myolinc, LncMyoD, Dum, MAR1, lnc-mg and Atrolnc-1 were significantly different compared with control C2C12, suggesting that these lncRNAs were involved in skeletal muscle differentiation (lincMD1, Myolinc, LncMyoD, Dum, MAR1 and lnc-mg) and atrophy (Atrolnc-1) under starvation conditions.

 miRNA expression patterns in atrophic C2C12 cells

As with the lncRNAs, the levels of various miRNAs associated with muscle development were detected in order to identify differences in regulatory miRNA expression between control and serum-starved C2C12 cells by RT-qPCR; miRNAs with expression trends consistent with previous literature (32,33,35–47) were selected for further experiments. It was revealed that nine miRNAs matched these criteria. In starved C2C12 cells, five of the miRNAs (miR-133a, miR-133b, miR-206, miR-186 and miR-23a) exhibited no significant change in expression levels (Fig. 3A-E). Four miRNAs showed different levels of expression between control and serum-starved C2C12 cells. miR-27b and miR-29b expression (Fig. 3F and G) was significantly decreased in starved C2C12 cells; however, miR-1a and miR-18a expression levels were significantly higher (Fig. 3H and I).

 Establishment of a muscular atrophic mouse model induced by starvation in vivo

The ratio of tibial anterior muscle to body weight was decreased significantly in starved mice compared with normally fed control mice (Fig. 4A). H&E staining revealed that the CSA of mouse tibialis anterior muscles was significantly decreased following food deprivation for 48 h (Fig. 4B and C). The protein expression levels of MYHC were decreased in starved mouse tibialis anterior muscles compared with control mouse muscles (Fig. 4D). Additionally, the transcript and protein expression levels of MURF1 and atrogin-1 were increased in starved tibialis anterior muscles compared with control mouse muscles (Fig. 4E and F). These results indicated that starvation induced muscle atrophy in mice.

 lncRNA expression patterns during mouse starvation

The expression levels of the seven lncRNAs that were selected in the present study were detected in the tibialis anterior muscles of the starved and control mice by RT-qPCR. Six of the lncRNAs (Atrolnc-1, lincMD1, Myolinc, LncMyoD, Dum and MAR1) exhibited significantly higher levels in starved tibialis anterior muscles compared with control muscles (Fig. 5A-F). Significantly lower lnc-mg levels were observed in starved tibialis anterior muscles (Fig. 5G). The differences in the expression levels of the seven lncRNAs in starved mice reflected the patterns of differentiated serum-starved C2C12 cells in vitro. The significantly differential expression of lncRNAs between control and starved mice indicated that these lncRNAs are associated with skeletal muscle proliferation, differentiation and atrophy under starvation conditions.

 miRNA expression patterns during mouse starvation

In starved mouse tibialis anterior muscles, three miRNAs (miR-206, miR-23a and miR-18a) showed no significant change in expression levels (Fig. 6C, E and I). Six miRNAs (miR-133a, miR-133b, miR-186, miR-27b and miR-1a) showed significantly decreased levels of expression, following starvation in mice (Fig. 6A, B, D, F and H). By contrast, the miR-29b expression levels (Fig. 6G) were significantly higher in starved mice compared with control. The expression patterns of miR-23a, miR-206 and miR-27b in starved mice exhibited the same changes as those of the differentiated serum-starved C2C12 cells in vitro, whereas the expression patterns of the other six miRNAs differed between the two starvation models.

 Discussion

Muscular atrophy is beginning to be considered in the pathogenesis of many physiological and pathological conditions, such as starvation and severe deprivation of nutrients, muscle disuse, Duchenne progressive muscular dystrophy, cancer and chronic disease. Therefore, it is important to elucidate the molecular regulatory mechanisms underlying dysfunctional muscle anabolism during skeletal muscle atrophy, and to develop strategies to prevent and treat muscular atrophy.

Isolation of muscles from mice revealed few vascular endothelial cells, lymphatic vessels and nerves among the muscle fibers. However, it is difficult at present to purify the skeletal muscle of mice from contamination with endothelial cells or other cell types due to the limitation of current technical tools. Based on other studies on lncRNAs and miRNAs in the skeletal muscle (19,24,32), muscle samples were harvested directly without considering the effects of endothelial cells or any other cell types. Therefore, the entire tibialis anterior muscle was used to explore gene expression and the underlying molecular regulatory mechanisms of muscle dysfunctional metabolism.

Myotube diameter and MYHC intensity are closely associated with skeletal muscle function (54). MYHC participates in the formation of the cytoskeleton and provides the force for muscles to contract normally. In the cell starvation model of the present study, cell myotube diameter was dramatically thinner and MYHC protein expression levels were significantly decreased. Changes in myotube diameter and MYHC expression following starvation may indicate changes in muscle functions.

The present study analyzed the expression patterns of lncRNAs and miRNAs during starvation-induced muscular atrophy in vitro and in vivo. Several lncRNAs, including lnc-mg, Atrolnc-1, lincMD1, Myolinc, lncMyoD, Dum and MAR1 differed in their expression levels between control and atrophic muscles. Lnc-mg (24,26), lincMD1 (18) and Dum (22) regulate the expression of MYHC and have important roles in skeletal muscle cell differentiation and development. Lnc-mg regulates IGF-2 to promote muscle differentiation and development, by functioning as a competitive endogenous RNA for miR-125b (24). In the muscular atrophic models (in vitro and in vivo) of the present study, lnc-mg and MYHC were both decreased. Thus, it is speculated that downregulated lnc-mg may serve an important role in starvation-induced muscular atrophy and could be positively associated with downregulated MYHC. The increased expression of lincMD1 and Dum may be compensatory responses to atrophy. Atrolnc-1 expression was significantly higher in the starvation-induced atrophic models (in vitro and in vivo). This result was consistent with the fact that Atrolnc-1 is remarkably enhanced in atrophic muscles in atrophic mouse models of chronic kidney disease, starvation and cancer (3). Atrolnc-1 interacts with the A20 binding inhibitor of nuclear factor κB-1 (NF-κB-1), promoting its activation and resulting in increased MURF1 (3). The high expression of Atrolnc-1 observed in the present study further indicates that muscular atrophy occurred in starved cells and mice.

LincMD1 regulates the expression of mastermind like transcriptional coactivator 1 and myocyte enhancer factor 2C to promote muscle differentiation by acting as a sponge for miR-133 and miR-135 (18). The present demonstrated that lincMD1 expression levels significantly increased during muscular atrophy, during which myoblasts are required for differentiation and fusion in order to replace damaged myofibers. In addition, lincMD1 promotes the expression of human antigen R protein to consolidate the muscle regulatory factor mRNAs (20,55). The expression levels of the other four lncRNAs, Myolinc, LncMyoD, Dum and MAR1, were all elevated under atrophic conditions, which may be explained by the fact that myoblasts are required for differentiation and fusion, in order to replace damaged myofibers under atrophic conditions.

The expression patterns of miR-206, miR-23a and miR-27b in starved mice exhibited the same changes as those in differentiated serum-starved C2C12 cells in vitro, whereas the expression patterns of the other six miRNAs differed between the two starvation models. miR-23a and miR-206 exhibited no significant change in expression levels, indicating that their roles in starvation-induced skeletal muscle atrophy may be minimal. miR-27b promotes muscle differentiation and inhibits proliferation by inhibiting MyoD family inhibitor (MDFI) expression (44). MDFI functions as a negative regulator of muscle growth by inhibiting muscle hypertrophy (56). However, miR-27b expression was significantly decreased in the starvation-induced atrophic models, possibly because differentiation was decreased under atrophic conditions. miR-29b expression in the starvation-induced atrophic mouse model was significantly higher than in control mice, consistent with a study that demonstrated that miR-29b was elevated in an in vivo atrophy model (32). However, the expression of miR-29b in serum-starved C2C12 myotubes in vitro was significantly reduced compared to control myotubes, which was not consistent with our in vivo atrophy model; the role of miR-29b in primary myoblasts therefore remains unclear, which was a limitation of the present study, and the reasons for this phenomenon remain to be determined.

The seven lncRNAs selected exhibited the same expression trends in the in vivo and in vitro starved models. Among the nine miRNAs selected, three miRNAs showed consistent change in trends in the two models, whereas six miRNAs showed inconsistent change. Therefore, the consistency of lncRNA expression was superior; the inconsistent changes in the trends of miRNAs in the two starved models require further investigation.

In conclusion, the present study analyzed the expression patterns of lncRNAs and miRNAs in starvation/fasting-induced muscular atrophy models. The current findings indicated that lncRNAs had similar expression changes in vitro and in vivo. The present results provided a novel insight into the lncRNAs involved in muscle atrophy and suggested that they may serve as potential diagnostic tools. A follow-up investigation is required to better understand the roles of lncRNAs in muscle atrophy.

 Acknowledgements

Not applicable.

 Funding

The current study was supported by grants from the Natural Science Foundation of Guangdong Province (grant no. 2018A030313591), the Science and Technology Planning Project of Guangdong Province (grant no. 2017ZC0333), the Administration of Traditional Chinese Medicine of Guangdong Province (grant no. 20181014), the Science and Technology Planning Project of Haizhu District (grant no. 2018-87), the Science Foundation of Guangdong Second Provincial General Hospital (grant nos. YQ2017-012, YN2017-002, YN2017-003 and YQ2018-005) and the National Natural Science Foundation of China (grant no. 81601503).

 Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

 Authors' contributions

RC conceived and designed the experiments. SL, YS and RC performed the experiments in vivo and analyzed the data. SL, YS, JZ, SZ and HS performed the experiments in vitro and analyzed the data. SL and RC wrote the manuscript. All authors read and approved the final manuscript.

 Ethics approval and consent to participate

The present study was approved by the Ethics Committee of Guangdong Second Provincial General Hospital (approval no. 2019-YJSWZ-001).

 Patient consent for publication

Not applicable.

 Competing interests

The authors declare that they have no competing interests.

 References WallimannTDolderMSchlattnerUEderMHornemannTKraftTStolzMCreatine kinase: An enzyme with a central role in cellular energy metabolismMAGMA6116119199810.1007/BF026609279803379 ZhaoCShangLWangWJacobsDOMyocellular creatine and creatine transporter serine phosphorylation after starvationJ Surg Res1051016200210.1006/jsre.2002.643112069495 SunLSiMLiuXChoiJMWangYThomasSSPengHHuZLong-noncoding RNA Atrolnc-1 promotes muscle wasting in mice with chronic kidney diseaseJ Cachexia Sarcopenia Muscle9962974201810.1002/jcsm.12321300434446204593 ZhangZKLiJGuanDLiangCZhuoZLiuJLuAZhangGZhangBTLong noncoding RNA lncMUMA reverses established skeletal muscle atrophy following mechanical unloadingMol Ther2626692680201810.1016/j.ymthe.2018.09.01430415659 ChenRJiangTSheYXuJLiCZhouSShenHShiHLiuSEffects of cobalt chloride, a hypoxia-mimetic agent, on autophagy and atrophy in skeletal C2C12 myotubesBiomed Res Int201770975802017287069505494548 ChenRSheYFuQChenXShiHLeiSZhouSOuJLiuYDifferentially expressed coding and noncoding RNAs in CoCl2-induced cytotoxicity of C2C12 cellsEpigenomics11423438201910.2217/epi-2018-008730785338 MakRHIkizlerATKovesdyCPRajDSStenvinkelPKalantar-ZadehKWasting in chronic kidney diseaseJ Cachexia Sarcopenia Muscle2925201110.1007/s13539-011-0026-6214756753063874 ZwartSRDavis-StreetJEPaddon-JonesDFerrandoAAWolfeRRSmithSMAmino acid supplementation alters bone metabolism during simulated weightlessnessJ Appl Physiol (1985)99134140200510.1152/japplphysiol.01406.200415691900 AlzghoulMbGerrardDWatkinsBAHannonKEctopic expression of IGF-I and shh by skeletal muscle inhibits disuse-mediated skeletal muscle atrophy and bone osteopenia in vivoFASEB J18221223200410.1096/fj.03-0293fje14597562 GregoryCMVandenborneKHuangHFOttenwellerJEDudleyGAEffects of testosterone replacement therapy on skeletal muscle after spinal cord injurySpinal Cord412328200310.1038/sj.sc.310137012494317 FinkleWDGreenlandSRidgewayGKAdamsJLFrascoMACookMBFraumeniJFJrHooverRNIncreased risk of non-fatal myocardial infarction following testosterone therapy prescription in menPLoS One9e85805201410.1371/journal.pone.0085805244896733905977 MilanGRomanelloVPescatoreFArmaniAPaikJHFrassonLSeydelAZhaoJAbrahamRGoldbergALRegulation of autophagy and the ubiquitin-proteasome system by the FoxO transcriptional network during muscle atrophyNat Commun66670201510.1038/ncomms7670258588074403316 Cid-DíazTSantos-ZasIGonzález-SánchezJGurriarán-RodríguezUMosteiroCSCasabiellXGarcía-CaballeroTMoulyVPazosYCamiñaJPObestatin controls the ubiquitin-proteasome and autophagy-lysosome systems in glucocorticoid-induced muscle cell atrophyJ Cachexia Sarcopenia Muscle8974990201710.1002/jcsm.12222286756645700440 MillerBFBaehrLMMusciRVReidJJPeelorFFIIIHamiltonKLBodineSCMuscle-specific changes in protein synthesis with aging and reloading after disuse atrophyJ Cachexia Sarcopenia MuscleJul162019(Epub ahead of print). doi: 10.1002/jcsm.1247010.1002/jcsm.1247031313502 TanakaMSugimotoKFujimotoTXieKTakahashiTAkasakaHKurinamiHYasunobeYMatsumotoTFujinoHRakugiHPreventive effects of low-intensity exercise on cancer cachexia-induced muscle atrophyFASEB J3378527862201910.1096/fj.201802430R30916585 BoltañaSValenzuela-MirandaDAguilarAMackenzieSGallardo-EscárateCLong noncoding RNAs (lncRNAs) dynamics evidence immunomodulation during ISAV–Infected Atlantic salmon (Salmo salar)Sci Rep622698201610.1038/srep22698269397524778034 ChenRJiangTSheYXieSZhouSLiCOuJLiuYComprehensive analysis of lncRNAs and mRNAs with associated co-expression and ceRNA networks in C2C12 myoblasts and myotubesGene647164173201810.1016/j.gene.2018.01.03929331478 CesanaMCacchiarelliDLegniniISantiniTSthandierOChinappiMTramontanoABozzoniIA long noncoding RNA controls muscle differentiation by functioning as a competing endogenous RNACell147358369201110.1016/j.cell.2011.10.031220000143234495 GongCLiZRamanujanKClayIZhangYLemire-BrachatSGlassDJA long non-coding RNA, LncMyoD, regulates skeletal muscle differentiation by blocking IMP2-mediated mRNA translationDev Cell34181191201510.1016/j.devcel.2015.05.00926143994 LegniniIMorlandoMMangiavacchiAFaticaABozzoniIA feedforward regulatory loop between HuR and the long noncoding RNA linc-MD1 controls early phases of myogenesisMol Cell6506514201410.1016/j.molcel.2013.12.012 MilitelloGHosenMRPonomarevaYGellertPWeirickTJohnDHindiSMMamchaouiKMoulyVDöringCA novel long non-coding RNA myolinc regulates myogenesis through TDP-43 and Filip1J Mol Cell Biol10102117201810.1093/jmcb/mjy02529618024 WangLZhaoYBaoXZhuXKwokYKSunKChenXHuangYJauchREstebanMALncRNA Dum interacts with Dnmts to regulate Dppa2 expression during myogenic differentiation and muscle regenerationCell Res25335350201510.1038/cr.2015.21256866994349245 ZhangZLiJGuanDLiangCZhuoZLiuJLuAZhangGZhangBTA newly identified lncRNA MAR1 acts as a miR-487b sponge to promote skeletal muscle differentiation and regenerationJ Cachexia Sarcopenia Muscle9613626201810.1002/jcsm.12281295123575989759 ZhuMLiuJXiaoJYangLCaiMShenHChenXMaYHuSWangZLnc-mg is a long non-coding RNA that promotes myogenesisNat Commun814718201710.1038/ncomms14718282815285353601 XiongWJiangYAiYLiuSWuXRCuiJGQinJYLiuYXiaYXJuYHMicroarray analysis of long non-coding RNA expression profile associated with 5-fluorouracil-based chemoradiation resistance in colorectal cancer cellsAsian Pac J Cancer Prev1633953402201510.7314/APJCP.2015.16.8.339525921151 DuJZhangPZhaoXHeJXuYZouQLuoJShenLGuHTangQMicroRNA-351-5p mediates skeletal myogenesis by directly targeting lactamase-β and is regulated by lnc-mgFASEB J3319111926201910.1096/fj.201701394RRR30216112 NeguemborMVJothiMGabelliniDLong noncoding RNAs, emerging players in muscle differentiation and diseaseSkelet Muscle48201410.1186/2044-5040-4-8246850023973619 Simionescu-BankstonAKumarANoncoding RNAs in the regulation of skeletal muscle biology in health and diseaseJ Mol Med (Berl)94853866201610.1007/s00109-016-1443-y273774064957971 BartelDPMicroRNAs: Genomics, biogenesis, mechanism, and functionCell116281297200410.1016/S0092-8674(04)00045-514744438 IveyKNSrivastavaDMicroRNAs as developmental regulatorsCold Spring Harb Perspect biol7a008144201510.1101/cshperspect.a008144261343124484971 LeiZSluijterJPvan MilAMicroRNA therapeutics for cardiac regenerationMini Rev Med Chem15441451201510.2174/138955751566615032412391325807943 LiJChanMCYuYBeiYChenPZhouQChengLChenLZieglerORoweGCmiR-29b contributes to multiple types of muscle atrophyNat Commun815201201710.1038/ncomms15201285412895458521 HorakMNovakJBienertova-VaskuJMuscle-specific microRNAs in skeletal muscle developmentDev Biol410113201610.1016/j.ydbio.2015.12.01326708096 Lagos-QuintanaMRauhutRYalcinAMeyerJLendeckelWTuschlTIdentification of tissue-specific microRNAs from mouseCurr Biol12735739200210.1016/S0960-9822(02)00809-612007417 WaldenTBTimmonsJAKellerPNedergaardJCannonBDistinct expression of muscle-specific MicroRNAs (myomirs) in brown adipocytesJ Cell Physiol218444449200910.1002/jcp.2162118937285 ChenJMandelEMThomsonJMWuQCallisTEHammondSMConlonFLWangDZThe role of microRNA-1 and microRNA-133 in skeletal muscle proliferation and differentiationNat Genet38228233200610.1038/ng172516380711 McCarthyJJEsserKAMicroRNA-1 and microRNA-133a expression are decreased during skeletal muscle hypertrophyJ Appl Physiol (1985)102306313200710.1152/japplphysiol.00932.200617008435 BoutzPLChawlaGStoilovPBlackDLMicroRNAs regulate the expression of the alternative splicing factor nPTB during muscle developmentGene Dev217184200710.1101/gad.150070717210790 AndersonCCatoeHWernerRMIR-206 regulates connexin43 expression during skeletal muscle developmentNucleic Acids Res3458635871200610.1093/nar/gkl743170626251635318 KimHKLeeYSSivaprasadUMalhotraADuttaAMuscle-specific microRNA miR-206 promotes muscle differentiationJ Cell Biol174677687200610.1083/jcb.200603008169238282064311 AntoniouAMastroyiannopoulosNPUneyJBPhylactouLAMiR-186 inhibits muscle cell differentiation through myogenin regulationJ Biol Chem28939233935201410.1074/jbc.M113.507343243854283924261 GuanLHuXLiuLXingYZhouZLiangXYangQJinSBaoJGaoHbta-miR-23a involves in adipogenesis of progenitor cells derived from fetal bovine skeletal muscleSci Rep743716201710.1038/srep43716282551765334644 MercatelliNFittipaldiSDe PaolaEDimauroIParonettoMPJacksonMJCaporossiDMiR-23-TrxR1 as a novel molecular axis in skeletal muscle differentiationSci Rep77219201710.1038/s41598-017-07575-0287753215543121 HouLXuJJiaoYLiHPanZDuanJGuTHuCWangCMiR-27b promotes muscle development by inhibiting MDFI expressionCell Physiol Biochem4622712283201810.1159/00048959529734192 LiuCWangMChenMZhangKGuLLiQYuZLiNMengQMiR-18a induces myotubes atrophy by down-regulating IgfIInt J Biochem Cell Biol90145154201710.1016/j.biocel.2017.07.02028782600 LiuCChenMWangMPiWLiNMengQMiR-18a regulates myoblasts proliferation by targeting Fgf1PLoS One13e2015512018 ChenJTaoYLiJDengZYanZXiaoXWangDZMicroRNA-1 and microRNA-206 regulate skeletal muscle satellite cell proliferation and differentiation by repressing Pax7J Cell Biol190867879201010.1083/jcb.200911036208199392935565 KallenANZhouXBXuJQiaoCMaJYanLLuLLiuCYiJSZhangHThe imprinted H19 lncRNA antagonizes let-7 microRNAsMol Cell52101112201310.1016/j.molcel.2013.08.02724055342 WangGWangYXiongYChenXCMaMLCaiRGaoYSunYMYangGSPangWJSirt1 AS lncRNA interacts with its mRNA to inhibit muscle formation by attenuating function of miR-34aSci Rep621865201610.1038/srep21865269026204763196 ChenRJiangTLeiSSheYShiHZhouSOuJLiuYExpression of circular RNAs during C2C12 myoblast differentiation and prediction of coding potential based on the number of open reading frames and N6-methyladenosine motifsCell Cycle1718321845201810.1080/15384101.2018.1502575300804266133337 LiFLiXPengXSunLJiaSWangPMaSZhaHYuQHuoHGinsenoside Rg1 prevents starvation-induced muscle protein degradation via regulation of AKT/mTOR/FoxO signaling in C2C12 myotubesExp Ther Med1412411247201710.3892/etm.2017.4615287816215526193 The Ministry of Science and Technology of the People's Republic of ChinaGuidance suggestions for the care and use of laboratory animalsSep302006 LivakKJSchmittgenTDAnalysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) methodMethods25402408200110.1006/meth.2001.126211846609 KimWKimJParkHJeonJDevelopment of microfluidic stretch system for studying recovery of damaged skeletal muscle cellsMicromachines (Basel)9E671201810.3390/mi912067130567359 von RoretzCBeauchampPDi MarcoSGallouziIHuR and myogenesis: Being in the right place at the right timeBiochim Biophys Acta181316631667201110.1016/j.bbamcr.2011.01.03621315776 ChenCMKrautNGroudineMWeintraubHI-mf, a novel myogenic repressor, interacts with members of the MyoD familyCell86731741199610.1016/S0092-8674(00)80148-88797820

Establishment of the serum starvation-induced C2C12 cell atrophic model in vitro. (A and B) MYHC and DAPI staining showed different myotube diameters between control and serum-starved C2C12 cells (n=3). (C) MYHC protein expression levels were detected by western blotting (n=3). (D) mRNA (n=6) and (E) protein (n=3) expression levels of MURF1 and atrogin-1 were observed in control and serum-starved C2C12 cells. *P<0.05, **P<0.01 and ****P<0.0001. MYHC, myosin heavy chain; MURF1, muscle RING finger protein 1.

lncRNA expression patterns in serum-starved C2C12 cells. The expression levels of lncRNAs (A) Atrolnc-1, (B) lincMD1, (C) Myolinc, (D) LncMyoD, (E) Dum, (F) MAR1 and (G) lnc-mg were detected in control C2C12 cells and differentiated serum-starved C2C12 cells. All data are presented as the mean ± standard deviation. *P<0.05, **P<0.01 and ***P<0.001. (n=6). lncRNA, long non-coding RNA; lincMD1, long intergenic non-protein coding RNA, muscle differentiation 1; lnc-mg, myogenesis-associated lncRNA; MAR1, muscle anabolic regulator 1.

miRNA expression patterns in atrophic C2C12 cells. The expression levels of (A) miR-133a, (B) miR-133b, (C) miR-206, (D) mR-186, (E) miR-23a, (F) miR-27b, (G) miR-29b, (H) miR-1a and (I) miR-18a were detected in control C2C12 cells and serum-starved C2C12 cells. All data are presented as the mean ± standard deviation. *P<0.05 and **P<0.01 (n=6). miRNA/miR, microRNA; ns, no significance.

Establishment of the muscular atrophic mouse model induced by starvation in vivo. (A) The ratio of tibial anterior muscle to total body weight decreased in the starved mice. (B and C) Hematoxylin-eosin staining indicates decreased CSA between control and starved tibialis anterior muscles. (D) MYHC protein expression levels were detected by western lot analysis in starved and control tibialis anterior muscles. (E) mRNA and (F) protein expression levels of MURF1 and atrogin-1 in starved and control tibialis anterior muscles. All data are presented as the mean ± standard deviation. *P<0.05, **P<0.01 and ***P<0.001 (n=6). MYHC, myosin heavy chain; CSA, cross-sectional area; MURF1, muscle RING finger protein 1.

lncRNA expression patterns during mouse starvation. The expression levels of lncRNAs (A) Atrolnc-1, (B) lincMD1, (C) Myolinc, (D) LncMyoD, (E) Dum, (F) MAR1 and (G) lnc-mg were detected in tibialis anterior muscles from normally fed and starved mice. All data are presented as the mean ± standard deviation. *P<0.05, **P<0.01 and ***P<0.001 (n=6 mice per group). lncRNA, long non-coding RNA; lincMD1, long intergenic non-protein coding RNA muscle differentiation 1; lnc-mg, myogenesis-associated lncRNA; MAR1, muscle anabolic regulator 1.

miRNA expression patterns during mouse starvation. The expression levels of (A) miR-133a, (B) miR-133b, (C) miR-206, (D) mR-186, (E) miR-23a, (F) miR-27b, (G) miR-29b, (H) miR-1a and (I) miR-18a were detected in normally fed and starved mice. All data are presented as the mean ± standard deviation. *P<0.05, **P<0.01 and ****P<0.0001 (n=6 mice per group). miRNA/miR, microRNA; ns, no significance.

Sequences of mRNA and long non-coding RNA primers.

Gene Forward primer (5′-3′) Reverse primer (5′-3′)
18S GTAACCCGTTGAACCCCATT CCATCCAATCGGTAGTAGCG
MURF1 GTGTGAGGTGCCTACTTGCTC GCTCAGTCTTCTGTCCTTGGA
Atrogin-1 GAGTGGCATCGCCCAAAAGA TCTGGAGAAGTTCCCGTATAAGT
Atrolnc-1 CAGCTGCCTACACCTGAAGA AGGGCTCGCAGATTACACC
lincMD1 AGTGATTGAGGTGGACAGAAGG CCCATTGAGGAGCATAGAACC
Myolinc CGGTGCTATGGTTCTGATCG TATGTGGGAAATACAGGGACA
LncMyoD ACCCAAGGCAAGAAAAGTAGCA ACTCACGAGTCAGCGGCAGAAC
Dum CACAAAGACAGGCAGACAGAC TACCAAGCAGGTTCCTACGG
MAR1 CCAAAGGACTGTCTTGGAACA AACAGCACTGAGCAGGGAC
lnc-mg CTGCATCACGGAAGGAGATA AACAATCCATCCTCATTGGC

LincMD1, long intergenic non-protein coding RNA muscle differentiation 1; lnc-mg, myogenesis-associated lncRNA; MAR1, muscle anabolic regulator 1.

Sequences of miRNA primers.

miRNA Sequence (5′-3′)
U6 F: GGAACGATACAGAGAAGATTAGC
R: TGGAACGCTTCACGAATTTGCG
miR-133a TTTGGTCCCCTTCAACCAGC
miR-133b GGTCCCCTTCAACCAGCTA
miR-206-3p GGAATGTAAGGAAGTGTGTGG
miR-186 GAATTCTCCTTTTGGGCTAAAA
miR-23a-3p ATCACATTGCCAGGGATTTCC
miR-27b TTCACAGTGGCTAAGTTCTGC
miR-29b-3p TAGCACCATTTGAAATCAGTGTT
miR-1a-3p TGGAATGTAAAGAAGUATGTAT
miR-18a GCCATCTAGTGCAGATAGAAAA

miRNA/miR, microRNA; F, forward; R, reverse.