Abstract

Long noncoding RNA (lncRNA) is a crucial factor in the progression of insulin resistance (IR). Resveratrol (RSV) exhibits promising therapeutic potential for IR. However, there are few studies on whether RSV improves IR through lncRNA. This study aimed to determine whether RSV could influence the expression of lncRNA and to elucidate the underlying mechanism. Mice were divided into three groups: control group, high-fat diet (HFD) group, and HFD + RSV group. We conducted a high-throughput sequencing analysis to detect lncRNA and mRNA expression signatures and the ceRNA-network in the skeletal muscles of mice that were fed an HFD to induce IR. Hierarchical clustering, gene enrichment, and gene ceRNA-network analyses were subsequently conducted. Differentially expressed lncRNAs were selected and validated via reverse transcription-quantitative polymerase chain reaction (RT-qPCR). The biological functions of the selected lncRNAs were investigated by silencing the target genes via lentivirus transfection of C2C12 mouse myotube cells. RSV treatment reversed the expression of 338 mRNAs and 629 lncRNAs in the skeletal muscles of mice with HFD-induced IR. The results of the Gene Ontology and Kyoto Encyclopedia of Genes and Genomes database analyses indicated that the differentially expressed mRNAs modulated type II diabetes mellitus. After validating randomly selected lncRNAs via RT-qPCR, we identified a novel lncRNA, NONMMUT044897.2, which was upregulated in the HFD group and reversed with RSV treatment. Additionally, NONMMUT044897.2 was proven to function as a ceRNA of microRNA- (miR-) 7051-5p. Suppressor of Cytokine Signaling 1 (SOCS1) was confirmed as a target of miR-7051-5p. We further performed lentivirus transfection to knock down NONMMUT044897.2 in vitro and found that NONMMUT044897.2 silenced SOCS1 and potentiated the insulin signaling pathway. Hence, RSV mimicked the silencing effect of lentivirus transfection on NONMMUT044897.2. Our study revealed that RSV reduced IR in mouse skeletal muscles via the regulation of NONMMUT044897.2.

1. Introduction

Type 2 diabetes mellitus (T2DM) accounts for 90% of diabetes worldwide, and insulin resistance (IR) is a primary determinant of T2DM since it reduces glucose uptake and utilization [1]. Skeletal muscles play a significant role in the etiology of IR [2]. Phosphatidylinositol-3-kinase (PI3K)/protein kinase B (AKT) is the most important signaling pathway of insulin in glucose metabolism [3]. Activated AKT can promote phosphorylated GSK3β, inhibiting its activity to improve IR [4]. In recent years, researchers found that the Suppressor of Cytokine Signaling 1 (SOCS1) is an important negative regulator of the PI3K/AKT pathway [5, 6]. Overexpression of SOCS1 in the liver was found to reduce insulin sensitivity by downregulating the level of tyrosine phosphorylation of insulin receptor substrate in IR mice [5]. SOCS1 primarily inhibits the catalytic binding of insulin receptors to insulin receptor substrates to induce IR [6].

Resveratrol (3,5,4-trihydroxystilbene; RSV) is a natural polyphenol that is enriched in more than 70 kinds of plants [7]. Accumulating evidence has indicated that RSV has diverse biological activities [8, 9]. It functions as an antioxidant, antiaging, anti-inflammatory, hypoglycemic, and anti-IR agent [10, 11]. Several studies have reported that RSV possesses a significant anti-IR activity in skeletal muscles [1214]. Our earlier research has also proven that skeletal muscle IR that is caused by an HFD may be alleviated with RSV treatment [14]. RSV increases the expression of microencapsulated protein 3 (CAV-3), thereby allowing skeletal muscle cells to carry glucose when the protein GLUT4 activates the transfer from the cytoplasm to the cell membrane, which in turn increases the ability of myocytes to transport glucose and improve IR [15]. Meanwhile, RSV attenuates insulin-stimulated AKT phosphorylation by eliminating insulin-induced ROS production in skeletal muscles [16].

Long noncoding RNAs (lncRNAs), a class of RNA molecules that are more than 200 nucleotides (nt) in length, have little or no protein-coding capacity [17]. Research has shown that lncRNAs widely participate in various developmental and physiological processes [18, 19]. In addition, they are strongly correlated with the development and progression of diseases, including coronary artery diseases [20], cancers [21], and metabolic diseases [22, 23]. Recently, lncRNAs have been confirmed to inhibit miRNA activity, increase miRNA target genes, and act as competitive endogenous RNAs (ceRNA) of miRNAs [24]. Further, the effect of lncRNAs in IR has been the focus of several studies [25, 26]. Previous studies [27, 28] by our group have demonstrated that RSV can improve hepatic IR by regulating lncRNA NONMMUT058999.2 and NONMMUT008655.2 in mouse models. However, the effect of RSV on IR through the regulation of the expression of lncRNAs in skeletal muscle remains unclear. In this study, we aimed to establish in vivo and in vitro IR models in skeletal muscles, which are a different tissue than those used in previous studies, to explore whether RSV can improve skeletal muscle IR by regulating lncRNAs.

2. Materials and Methods

2.1. Animal Experiments

We equally divided 42 healthy 6-week-old C57BL/6J background mice into three groups (n = 14 in each group): the control group, the HFD group, and the HFD + RSV group. The mice weighing around 22 g were purchased from the Beijing Viton Lihua Experimental Center (China) and sustained on a standard 12 h light-dark cycle, at 20–25°C and at 40–60% humidity. HFD mice were fed D12492J (20% protein, 20% carbohydrate, and 60% fat) for 8 weeks. The feed was purchased from Beijing Huafukang Biotechnology Co., Ltd. HFD + RSV mice were intragastrically fed 100 mg/kg/day RSV solution for 6 weeks following the methods outlined in our previous study [27]. Dissolved RSV (Sigma Aldrich, St. Louis, MO, USA) with dimethyl sulfolane (Sigma Aldrich; 30 mg mL−1) was diluted with 0.9% NaCl in a ratio of 1 : 2. The control group was fed D12450J (20% protein, 70% carbohydrate, and 10% fat). Weight and food intake were measured weekly during feeding. After the feeding experiment, all mice fasting for 12 h were intraperitoneally injected with 50% glucose (1.5 g kg−1 bodyweight) to conduct the glucose tolerance tests (intraperitoneal glucose tolerance test). Blood glucose was detected on the tail vein with a glucose meter at 0, 15, 30, 60, and 120 min after injection. The IR model was validated by determining the area under the curve in accordance with the protocol presented in our previously published study [29]. Animal studies (2019E369) were approved by the ethics committee of the Hebei General Hospital, and all animal experimental procedures complied with the National Institutes of Health guide for the care and use of laboratory animals.

2.2. Serum and Tissue Samples

Six mice in each group were randomly selected and intraperitoneally injected with 1.5 U/40 g of insulin (Sigma Aldrich). All mice in three groups were anesthetized by intraperitoneal injection of 2% sodium pentobarbital after 20 min. Blood samples were gathered via cardiac puncture and centrifuged at 3000 ×  at 4°C for 10 min. The serum was then stored at −80°C for serological indicators. Collected skeletal muscles were withdrawn quickly and stored in liquid nitrogen for the follow-up study.

2.3. Serological Indicators

Detection kits for total cholesterol, triglyceride, high-density lipoprotein cholesterol, low-density lipoprotein cholesterol, and free fatty acids were acquired from Nanjing Jiancheng Institute of Biological Engineering (Jiangsu, China). Serum insulin was obtained using an ELISA kit (ALPCO Diagnostics, Salem, NH, USA). The manufacturer’s protocol was followed in all the aforementioned procedures involving the detection kits.

2.4. Western Blot

The same amount of protein with different groups was separated by sodium dodecyl sulfate-polyacrylamide gel electrophoresis, transferred to polyvinylidene fluoride (PVDF) membrane (cat. no. ISEQ00010; Merk Millipore, Billerica, MA, USA), and sealed with 5% skimmed milk for 2 h. The primary antibodies were diluted in the blocking solution at the following concentrations: β-actin (cat. no. 4970; Cell Signaling Technology, Danvers, MA, USA): rabbit antibody, 1 : 5000; GAPDH (cat. no. 10494-1-AP; Proteintech Group, Inc. 5400 Pearl Street, Suite 300 Rosemont, IL 60018, USA): rabbit antibody, 1 : 10000; AKT (cat. no. 9272; Cell Signaling Technology): rabbit antibody, 1 : 1000; p-AKT (Ser 473) (cat. no. 9271; Cell Signaling Technology): rabbit antibody, 1 : 750; GSK3β (cat. no. 12456; Cell Signaling Technology): rabbit antibody, 1 : 750; p-GSK3β (cat. no. 5558; Cell Signaling Technology): rabbit antibody, 1 : 750; GLUT4 (cat. no. 2213; Cell Signaling Technology): mouse antibody, 1 : 5000; SOCS1 (cat. no. 3950; Cell Signaling Technology): rabbit antibody, 1 : 5000. PVDF membranes and primary antibodies were incubated at 4°C for 24 h. The membrane was then treated with secondary antibodies, stored at 30°C for about 50 min and washed three times for 10 min each time. The washed membrane was fully immersed in the iPer ECL Western HRP Substrate (cat. no. MF074-01; Mei5, Beijing, China) for about 2 min and exposed using the Gel Imager System (GDS8000; UVP, California, USA) to obtain the image of the target band. Protein bands were calculated by densitometry using the ImageJ software and were normalized to β-actin or GAPDH levels.

2.5. Reverse Transcription-Quantitative Polymerase Chain Reaction (RT-qPCR)

Total RNA was extracted using the Trizol reagent (TIANGEN, Beijing, China) from mouse skeletal muscle tissues and was reverse-transcribed into cDNA using HiScript II Q RT SuperMix for qPCR. RNA was tested for purity and concentration using NanoDrop 2000 (Thermo Fisher Scientific, Wilmington, DE, USA). Amplification was performed using the SYBr® Premix ex Taq II kit (RR820A; Takara Bio, Tokyo, Japan). The Applied Biosystems 7500 real-time PCR system was used to perform RT-qPCR, with a total of 41 cycles, including 3 minutes of predenaturation at 95°C, 5 seconds at 95°C, and 32 seconds at 60°C. The melting point curve was established at 60–95°C. β-Actin and U6 were considered as internal reference controls for genes. The relative gene expression was quantified by the 2−(ΔΔCt) method [30]. The specific primers involved in this research are listed in Table 1.

2.6. cDNA Library Construction and RNA Sequencing

Total RNA was extracted using the RNeasy mini kit (Qiagen, Hilden, Germany) from skeletal muscle samples of four mice in each group. The TruSeq™ RNA Sample Preparation Kit (Illumina, San Diego, CA, USA) was used to make paired-end libraries following the manufacturer’s instructions. Ribosomal RNA was removed using the Ribo-Zero rRNA Removal Kit (Epicentre, Madison, Wisconsin, USA), and then the mRNA was fragmented into small pieces with divalent cations. First chain cDNA was synthesized using reverse transcriptase and random primers, and second chain cDNA was generated by DNA Polymerase I and RNase H. Base “A” was added at the end of the cDNAs. Purified products were PCR-amplified to create the final cDNA library. Insert sizes were confirmed using a Qubit® 2.0 Fluorometer (Life Technologies, Carlsbad, CA, USA), and mole concentrations were calculated. Clusters were created using cBot and then sequenced on the Illumina NovaSeq 6000 (Illumina) by Sinotech Genomics Co., Ltd. (Shanghai, China). Each gene fragment was counted using the Stringtie software (version: 1.3.0; Johns Hopkins University, Baltimore, MD, USA) contrast and then normalized using the TMM (trimmed mean of M values) algorithm. The Fragments Per Kilobase Million value of each gene was then calculated. The high-throughput sequencing results were uploaded to the gene expression omnibus database (accession no. GSE178415).

3. Analysis of Differentially Expressed lncRNA and mRNA

The differential expression of skeletal muscle in the three groups was analyzed based on the “edge” package in R. The threshold of the -value was confirmed by controlling the False Discovery Rate. The screening criteria for differential expression of mRNAs and lncRNAs were and absolute value of log2FC >1.0.

3.1. Functional Group Analysis

The Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways determined the potential role of the lncRNAs that were coexpressed with the differentially expressed mRNAs. The GO analysis was conducted to establish significant annotations of genes and gene products in diversified organisms using the DAVID database (https://david.abcc.ncifcrf.gov). In addition, the KEGG pathway analysis was used to identify differentially expressed mRNAs in enriched pathways. was identified as the threshold of significance.

3.2. Construction of a ceRNA Regulatory Network

Cytoscape (version 3.8.2) is a network visualization software with multiple applications for network analysis. It can be downloaded for free from https://www.cytoscape.org/.

3.3. C2C12 Cell Culture and Treatments

C2C12 mouse myotube cells were maintained in Dulbecco’s modified Eagle medium (Gibco, Waltham, MA, USA) at a density of 5 × 104 cells/cm2 containing 10% fetal bovine serum (San Diego, California, USA) and 1% penicillin/streptomycin (Wisent, Nanjing, China) at 37°C with 5% CO2. Cell differentiation was induced by incubation in Dulbecco’s modified Eagle medium containing 2% fetal bovine serum for 4 days after reaching 80% confluence. Differentiated C2C12 cells at a density of 1 × 104 cells/cm2 were incubated for 24 h with 0.25 mM palmitate (PA) (Aladdin Industries, China) [31]. At 0, 8, 16, and 24 h after the intervention of PA, the glucose concentration in the culture medium was measured by the glucose oxidase assay to determine whether IR in mice was established. Subcultured C2C12 cells were digested to prepare a cell suspension and subcultured to a 96-well culture plate. When the cells grew to about 80% confluence, RSV at concentrations of 100 μM, 50 μM, and 30 μM was added to the medium. After 24 h, 10 μL of CCK-8 was added to each well, which was protected against light. They were cultured for 20 min, and their absorbance was measured at 450 nm. The cell survival rate was then calculated. C2C12 cells (5 × 105 cells/cm2) in the logarithmic growth phase were subcultured in a 6-well plate and transfected with lentivirus. Synthesized constructs included LV3-NC (5′ to 3′ TTCTCCGAACGTGTCACGT) and LV3-NONMMUT044897.2 (5′ to 3′ GCTCTTTCAGATAAGCCTTGT), which were obtained from GenePharma Co., Ltd. China. Stable cell lines were obtained after puromycin (2 μg mL−1) selection for the PA-induced IR model and drug intervention experiments. The plated cells were grouped: the control group (CON), the PA group (PA), the PA + shRNA-NONMMUT044897.2 negative control group (PA + shRNA-NC), the PA + shRNA-NONMMUT044897.2 knockdown group (PA + shRNA-NONMMUT044897.2), and the PA + RSV 30 μM group (PA + RSV). After the IR model was established, the glucose concentration was measured for 24 h, while the NONMMUT044897.2 and miR7051-5p mRNA levels were measured via RT-qPCR. The cells at a density of 5 × 105 cells/cm2 were stimulated with 100 nM insulin, while the protein was extracted 20 min after insulin stimulation for the western blot analyses, and cells were extracted using the Trizol reagent (TIANGEN, Beijing, China).

3.4. Statistical Analyses

SPSS v23.0 was used for data analysis. The results are presented as a mean ± SD. Two-sample comparisons were analyzed using an independent sample t-test (Student’s t-test). One-way ANOVA was used for statistical analysis followed by Bonferroni’s multiple comparison test or Tamhane’s multiple comparison test. was regarded as statistically significant.

4. Results

4.1. RSV Ameliorates Body Weight, IR, and Lipid Levels in HFD-Fed Mice

After 6 weeks of RSV administration, body weight, fasting blood glucose, and insulin levels of the HFD + RSV group were greatly reduced compared with those of the HFD group (Table 2), although the daily caloric intake of the two groups was similar. The quantitative insulin sensitivity index of the HFD group was reduced compared with the indices of the CON and HFD + RSV groups (Table 2). RSV treatment reduced the upregulation of triglyceride, low-density lipoprotein cholesterol, and free fatty acids in the HFD group, while the total cholesterol was decreased but had no significance (Table 2). There was no difference in the high-density lipoprotein cholesterol of the mice (Table 2).

4.2. RSV Treatment Decreased SOCS1 and Increased the Phosphorylation of AKT and GSK3β in the HFD Group

Between the control, HFD, and HFD + RSV groups, no differences were found in the protein levels of AKT and GSK3β (Figures 1(a), 1(b), and 1(d)). The HFD group showed dramatic repression of p-AKT and p-GSK3β protein levels compared with those in the control group, while the RSV group showed a marked increase in p-AKT and p-GSK3β protein expression (Figures 1(a), 1(c), and 1(e)). Moreover, SOCS1 expression was abnormally elevated in the HFD group but decreased following RSV treatment (Figures 1(a) and 1(f)). These results suggest that RSV improves the expression of genes on the insulin signaling pathway.

4.3. RSV Systematically Modulates Skeletal Muscle Gene Expression

After standardization, 58,245 lncRNAs and 83,089 mRNAs were identified in the skeletal muscles of mice. On comparing the HFD group with the control group, we found that there were 3,276 differentially expressed lncRNAs (1,192 upregulated and 2,084 downregulated) and 2,118 differentially expressed mRNAs (314 upregulated and 1804 downregulated). Simultaneously, there were 1,640 differentially expressed lncRNAs (921 upregulated and 719 downregulated) and 604 differentially expressed mRNAs (444 upregulated and 160 downregulated) in the HFD + RSV group compared with those in the HFD group. Of the upregulated lncRNAs and mRNAs in the HFD group, 270 lncRNAs and 58 mRNAs were downregulated in the HFD + RSV group. Of the downregulated lncRNAs and mRNAs in the HFD group, 359 lncRNAs and 280 mRNAs were upregulated in the HFD + RSV group (Figures 2(a) and 2(b)). The top 30 differentially expressed lncRNAs and mRNAs are listed in Tables 3 and 4, with Fragments Per Kilobase Million = 0 eliminated. All differentially expressed lncRNAs and mRNAs in the three groups are listed in Supplementary Tables S1 and S2.

4.4. Functional Enrichment Analysis of Differentially Expressed Genes

The functions of these different mRNAs were studied via enrichment analysis. The GO analysis was conducted to classify differentially expressed mRNAs into three types: the biological process, the molecular function, and the cellular component. The most highly enriched GO terms were “myelination, ensheathment of neurons, axon ensheathment, cellular component assembly involved in morphogenesis, transition between fast and slow fiber (biological process),” “myofibril, sarcomere, contractile fiber, compact myelin, contractile fiber part (cellular component),” and “structural constituent of myelin sheath, actin binding, fatty acid synthase activity, ion gated channel activity, and gated channel activity (molecular function)” (Figure 3(a)). KEGG analysis revealed that the differentially expressed mRNAs were mostly involved in cytokine-cytokine receptor interaction, the JAK-STAT signaling pathway, hypertrophic cardiomyopathy, the prolactin signaling pathway, and type II diabetes mellitus (Figure 3(b)).

4.5. RT-qPCR Validation In Vivo

To confirm the validity of the sequencing results, we randomly picked four differentially expressed lncRNAs, two upregulated lncRNAs in HFD, a downregulated HFD + RSV (NONMMUT044897.2; NONMMUT005295.2), and two lncRNAs with a reverse trend (NONMMUT128951.1; NONMMUT145909.1). The expression levels of the selected lncRNAs were consistent with those of the sequencing results (Figures 4(a)4(d)), but there was no statistical difference in the increase of NONMMUT128951.1 in the HFD + RSV group. The SOCS1 mRNA level was increased in the HFD group and decreased with RSV treatment (Figure 4(f)). Of the verified lncRNAs, NONMMUT044897.2 had a higher expression level. In addition, the KEGG analysis indicated that SOCS1 played a vital role in the JAK-STAT signaling pathway and in type II diabetes mellitus. In the latter, SOCS1 was reportedly involved in the development of IR [5, 6]. To elucidate the interaction between NONMMUT044897.2 and SOCS1, we constructed a related lncRNA-miRNA-mRNA network diagram. The results revealed that NONMMUT044897.2 regulated SOCS1 through miR-7051-5p and miR-762 (Figure 4(g)). To further explore its potential molecular mechanism, we applied the NonCode and miRBase database analyses and found that there was base pairing in the sequence of NONMMUT044897.2 and miR-7051-5p. At the same time, according to the prediction results of TargetScan, SOCS1 was identified to be a miR-7051-5p target (Figure 4(h)). According to the above results, NONMMUT044897.2 might have regulated SOCS1 through miR-7051-5p. Therefore, this study further verified the expression of miR-7051-5p mRNA through RT-qPCR; miR-7051-5p mRNA was downregulated in the HFD group and upregulated in the HFD + RSV group (Figure 4(e)).

4.6. Establishment of a Cell Model of PA-Induced IR

The C2C12 mouse myotube cells were transferred to media with and without 0.25 mM PA. Glucose concentrations were determined at 0, 8, 16, and 24 h. There was no significant difference in control and PA groups at 0, 8, and 16 h; however, the glucose concentration of the PA group was distinctly elevated at 24 h compared with that of the control group, thereby indicating that the IR model was successfully established (Figure 5(a)). Besides, the glucose concentration was obviously decreased at 24 h in the RSV group (Figure 5(b)).

The cell survival rate of C2C12 cells 24 h after RSV administration was approximately 30–100 μm. The results showed that 30 μm of RSV had no significant influence on the survival rate of C2C12 cells (Figure 5(c)). The cell survival rate of the PA group (84%) was lower than that of the control (89.6%), while that of the 30 μM PA + RSV group (85%) was higher than that of the PA group. However, there were no statistical differences between these three groups (Figure 5(d)).

4.7. RT-qPCR Validation In Vitro

Two lncRNAs were upregulated in PA in which downregulation of PA + RSV (NONMMUT044897.2; NONMMUT139818.1) and upregulation of three lncRNAs (NONMMUT071570.2; NONMMUT065156.2; NONMMUT00000181045) were observed. The expression levels of the selected lncRNAs were consistent with those of the sequencing results (Figure 6(a)6(e)), but there was no statistical difference in the increase of NONMMUT00000181045 in the PA + RSV group. To verify the relationship between NONMMUT044897.2 and RSV in vitro, C2C12 cells were transfected with lentivirus. The results of this transfection showed that, compared with the control group, NONMMUT044897.2 expression was robustly increased in the PA and PA + shRNA-NC groups. The PA + shRNA-NONMMUT044897.2 group had a decreased expression of NONMMUT044897.2 compared with that of the PA group. RSV administration also resulted in a decrease in the expression of NONMMUT044897.2 compared with that of the PA group (Figure 6(g)).

Relative to the control group, miR-7051-5p mRNA expression in the PA and PA + shRNA-NC groups was significantly reduced, while knockdown of NONMMUT044897.2 and RSV treatment markedly increased the miR-7051-5p mRNA expression level (Figure 6(f)). The concentration of glucose in the media of the PA and PA + shRNA-NC groups was substantially increased compared with that of the control group. NONMMUT044897.2 silencing and RSV treatment strikingly overturned the glucose concentrations in the medium (Figure 6(h)). Knockdown of NONMMUT044897.2 distinctively upregulated the p-AKT, p-GSK3β, and GLUT4 protein levels and greatly reduced the SOCS1 protein level, compared with those of the PA group. RSV treatment had a similar effect in terms of NONMMUT044897.2 silencing on p-AKT, p-GSK3β, GLUT4, and SOCS1 protein levels (Figures 6(i)6(o)).

5. Discussion

A large number of experiments have verified the therapeutic effect of RSV in IR [32, 33]. In this study, we demonstrated that RSV influenced the amelioration of IR in HFD-fed mice. These findings have also been reported by several studies [11, 12]. In the HFD + RSV group, blood glucose, insulin index, blood lipids, and area under the curve were decreased. RSV treatment ameliorated HFD-induced IR in mice by restoring the insulin signaling pathway gene expression. After intervention with RSV, the protein expression levels of p-AKT, p-GSK3β, and GLUT4 increased significantly.

High-throughput sequencing showed that there were 3,276 differentially expressed lncRNAs and 2,118 differentially expressed mRNAs in the HFD + RSV group, as compared with those in the control group, which yielded 1,640 differentially expressed lncRNAs and 604 differentially expressed mRNAs. We further found 338 mRNAs and 629 lncRNAs whose expression was reversed in the HFD and the HFD + RSV groups, suggesting that RSV plays a significant role in the overall alteration of skeletal muscle gene expression. Moreover, RT-qPCR results revealed that RSV ameliorated IR by regulating the expression of lncRNAs in skeletal muscles. As mentioned above, the verified lncRNAs were consistent with those of the sequencing results. In addition, NONMMUT044897.2 was highly expressed. The KEGG analysis uncovered that the differential genes were part of T2DM. We found that NONMMUT044897.2 was associated with SOCS1, which was critically involved in T2DM. Overexpression of SOCS1 aggravated IR [5], which played an important role in T2DM. Hence, we selected this lncRNA for further study, which has not been reported before.

SOCS1 is a specific negative regulator that regulates the JAK/STAT pathway [34]. The expression of SOCS1 increases with IR and decreases with the phosphorylation of IRS-1. Overexpression of SOCS-1 could inhibit insulin-induced glycogen synthesis in L6 myotubes [6]. AKT has essential roles in many signaling pathways, such as cell survival and cell metabolism. AKT is the center of the insulin signaling pathway, which regulates glucose and lipid metabolism. Activated AKT can stimulate the translocation of insulin-sensitive GLUT4 to the cell membrane through its downstream substrate ASl60 to increase glucose uptake. It can also phosphorylate GSK3β to inhibit its activity, promote glycogen synthesis, lower blood sugar, and improve IR [35]. Furthermore, an HFD could result in the decrease of skeletal muscle IRS-1, P13K, AKT, and GLUT4 gene expression levels and reduce p-AKT (ser473) and p-GSK3β protein expression levels. Studies have reported that overexpression of SOCS1 may inhibit the phosphorylation and activation of IRS-1 [6, 34], which in turn inhibits the activation of AKT, thereby indicating that there is an important link between AKT and SOCS1. We found that, in the IR model mice, mRNA and protein expression levels of SOCS1 were significantly increased, while RSV treatment exhibited the reverse trend, thereby improving IR and decreasing blood glucose levels.

Numerous studies have demonstrated that lncRNAs may be involved in human diseases by regulating the expression of miRNAs [36, 37]. The ceRNA-network diagram uncovered the possible regulatory roles of candidate lncRNAs. NONMMUT044897.2 could regulate SOCS1 through two different miRNAs. To further clarify the relationship between NONMMUT044897.2 and SOCS1, we constructed NONMMUT044897.2 and miR-7051-5p, miR-7051-5p, and SOCS1 base-pairing maps based on NonCode, miRBase, and TargetScan databases. These results indicate that NONMMUT044897.2 might have regulated the expression of SOCS1 through miR-7051-5p. Future studies should perform luciferase assays to verify the interactions between NONMMUT044897.2, miR-7051-5p, and SOCS1.

We confirmed that knockdown of NONMMUT044897.2 increased miR-7051-5p levels and promoted the expression of genes that are involved in the insulin signaling pathways (p-AKT, p-GSK3β, and GLUT4). Meanwhile, the expression of SOCS1 was suppressed by silencing NONMMUT044897.2. Moreover, knockdown of NONMMUT044897.2 has led to reduced glucose concentration, which is similar to the phenotypes induced with RSV treatment. This indicates that RSV improves skeletal muscle IR through downregulating the lncRNA NONMMUT044897.2.

6. Conclusions

This research profiled the differential expression of lncRNAs between the IR model mice and those treated with RSV. We further revealed the potentially regulated lncRNA NONMMUT044897.2. More work remains to be done to prove the relationship between the NONMMUT044897.2/miR-7051-5p/SOCS1 and RSV in the IR model. The study also has several limitations. First, in vivo animal models are needed to further silence NONMMUT044897.2 to verify the influence on IR. Second, the overexpression of NONMMUT044897.2 should be performed in vitro. Third, knockdown of the NONMMUT044897.2 in the RSV group should be done to observe its influence on the IR model. Overall, our data indicated that RSV could promote skeletal muscle IR, at least partially, via a lncRNA NONMMUT044897.2/miR-7051-5p/SOCS1 pathway. This provides a new perspective for the RSV treatment of IR in skeletal muscles.

Data Availability

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

Disclosure

This article has been preprinted at https://www.researchsquare.com/article/rs-315719/v1.

Conflicts of Interest

The authors declare that there are no conflicts of interest.

Acknowledgments

The authors sincerely thank the teachers at the Clinical Medical Research Centre of Hebei General Hospital who helped them during the experiment. This study was supported by a grant from the Natural Science Foundation of Hebei Province (no. H201830-7071).

Supplementary Materials

All differentially expressed lncRNAs and mRNAs in the control, HFD, and HFD + RSV groups are listed in Supplementary Tables S1 and S2. (Supplementary Materials)