Next Article in Journal
Illegitimate and Repeated Genomic Integration of Cell-Free Chromatin in the Aetiology of Somatic Mosaicism, Ageing, Chronic Diseases and Cancer
Previous Article in Journal
The Influence of Diet and Obesity on Gene Expression in SLE
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 10
Issue 6
10.3390/genes10060406
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Transcriptome Changes of Skeletal Muscle RNA-Seq Speculates the Mechanism of Postprandial Hyperglycemia in Diabetic Goto-Kakizaki Rats During the Early Stage of T2D

by
Wenlu Zhang
,
Yuhuan Meng
,
Shuying Fu
,
Xingsong Li
,
Zixi Chen
,
Lizhen Huang
and
Hongli Du
*
School of Biology and Biological Engineering, South China University of Technology, Guangzhou 510000, China
*
Author to whom correspondence should be addressed.
Genes 2019, 10(6), 406; https://doi.org/10.3390/genes10060406
Submission received: 12 March 2019 / Revised: 20 May 2019 / Accepted: 23 May 2019 / Published: 28 May 2019
(This article belongs to the Section Animal Genetics and Genomics)
Browse Figures
Versions Notes

Abstract

:
To address how skeletal muscle contributes to postprandial hyperglycemia, we performed skeletal muscle transcriptome analysis of diabetic Goto-Kakizaki (GK) and control Wistar rats by RNA sequencing (RNA-Seq). We obtained 600 and 1785 differentially expressed genes in GK rats compared to those Wistar rats at three and four weeks of age, respectively. Specifically, Tbc1d4, involved in glucose uptake, was significantly downregulated in the skeletal muscle of GK aged both three and four weeks compared to those of age-matched Wistar rats. Pdk4, related to glucose uptake and oxidation, was significantly upregulated in the skeletal muscle of GK aged both three and four weeks compared to that of age-matched Wistar rats. Genes (Acadl, Acsl1 and Fabp4) implicated in fatty acid oxidation were significantly upregulated in the skeletal muscle of GK aged four weeks compared to those of age-matched Wistar rats. The overexpression or knockout of Tbc1d4, Pdk4, Acadl, Acsl1 and Fabp4 has been reported to change glucose uptake and fatty acid oxidation directly in rodents. By taking the results of previous studies into consideration, we speculated that dysregulation of key dysregulated genes (Tbc1d4, Pdk4, Acadl, Acsl1 and Fabp4) may lead to a decrease in glucose uptake and oxidation, and an increase in fatty acid oxidation in GK skeletal muscle at three and four weeks, which may, in turn, contribute to postprandial hyperglycemia. Our research revealed transcriptome changes in GK skeletal muscle at three and four weeks. Tbc1d4, Acadl, Acsl1 and Fabp4 were found to be associated with early diabetes in GK rats for the first time, which may provide a new scope for pathogenesis of postprandial hyperglycemia.

1. Introduction

Type 2 diabetes (T2D), the most prevalent form of diabetes, has affected more than 380 million people worldwide in 2017. Over a quarter of patients suffering from T2D are Chinese (109 million) [1]. Isolated postprandial hyperglycemia is one of the characteristics of diabetes, with the diagnosis criteria requiring a 2-hour plasma glucose concentration ≥11.1 mmol/L and a fasting plasma glucose concentration <7.0 mmol/L, according to a World Health Organization Consultation [2]. However, due to its inconvenient detection method, diagnoses of T2D using isolated postprandial hyperglycemia was not recommended for routine clinical use, with fasting plasma glucose used more commonly in practical diagnoses. Isolated postprandial hyperglycemia is common in Asian diabetes subjects [3], and 46.6% of undiagnosed diabetes subjects in China have isolated postprandial hyperglycemia [4]. The data from a Diabetes Epidemiology Collaborative Analysis of Diagnosis Criteria in Asia (DECODA) study showed that more than 50% of Asian diabetic patients had isolated postprandial hyperglycemia [5]. These findings indicate that isolated postprandial hyperglycemia is one of the main characteristics of Asian or Chinese T2D at the early stage. Therefore, researching the postprandial hyperglycemia mechanism is important for prevention and treatment of T2D. However, there are few researches concerning the molecular mechanisms of postprandial hyperglycemia in T2D patients at the early stage.
Goto-Kakizaki (GK) rats are a non-obese animal model for T2D, generated by selective breeding from Wistar rats with impaired glucose tolerance (IGT) [6,7]. GK rats display postprandial hyperglycemia and glucose intolerance as early as the first week after birth, and show fasting hyperglycemia after four weeks [8], suggesting that postprandial hyperglycemia appears earlier than fasting hyperglycemia in GK rats at early stages. Asian or Chinese T2D patients also have the characteristic of isolated postprandial hyperglycemia at the early stage. Insulin resistance has not been observed in GK rats from one to four weeks of age, while increased insulin sensitivity in GK rats from ages one to three weeks has been reported [8,9]. Thus, postprandial hyperglycemia in GK rats during the early stage of three and four weeks may not result from peripheral insulin resistance. As skeletal muscle is responsible for 70–80% postprandial glucose disposal [10,11], it may play a key role in postprandial hyperglycemia. However, the molecular mechanisms of how skeletal muscle contributes to postprandial hyperglycemia remain unclear. Taken together, GK rats are a suitable model for investigating the molecular mechanisms of skeletal muscle in postprandial hyperglycemia.
To investigate the role of skeletal muscle on postprandial hyperglycemia, diabetic GK and control Wistar rats aged three and four weeks were used in our study. We performed transcriptome and quantitative RT-qPCR analysis of the skeletal muscle of GK and Wistar rats to reveal the molecular mechanisms of postprandial hyperglycemia of T2D. Our study may provide a new scope for dysfunction research of postprandial hyperglycemia in the future.

2. Materials and Methods

2.1. Animal Breeding and Tissues Collection

In the present study, 10 male rats in each of four groups (GK rats aged three weeks, Wistar rats aged three weeks, GK rats aged four weeks and Wistar rats aged four weeks) were used, totaling 40 subjects. Rats were raised under a 12-hour light-12-hour dark cycle, at a temperature of 20 to 25 °C and 60 ± 5% humidity at the SLAC Laboratory Animal Co., Ltd. (Shanghai, China). All protocols were carried out in accordance with the approved guidelines of the Institutional Animal Care and Use Committee (IACUC) with Ethics certificate IACUC2014029. The rats were not fasted overnight. Blood was collected from the orbital plexus behind the eyeball using EDTA (4 mM final concentration) as an anticoagulant [12,13]. Plasma was obtained from the blood samples by centrifugation (2000×g, 4 °C, 15 min), divided into aliquots, and stored at −80 °C. All rats were killed by cervical dislocation within 2 hours. We harvested the red part of the gastrocnemius muscle from all 10 rats in each group, placed the muscle portion into liquid nitrogen immediately, and then put the portions into a −80 °C freezer for future studies. The rats were weaned at the end of the second week. Six samples from each group were used for RNA sequencing (RNA-Seq).

2.2. Measurement of Plasma Glucose and Insulin Concentration

Random plasma glucose concentrations were measured using an automatic Dry Biochemical Analyzer FUJIFILM DRI-CHEM 7000i (Saitama, Japan). Plasma insulin concentrations were measured using a Thermo scientific Rat Insulin ELISA Kit (Thermo Scientific, Waltham, MA, USA), in accordance with the manufacturer’s instructions.

2.3. RNA Extraction, Library Construction and Sequencing

The total RNA of the red part of the gastrocnemius muscle was extracted with TRIzol reagent (Cat#15596-018, Life Technologies, Carlsbad, CA, USA), following the manufacturer’s instructions. The kaiaoK5500® Spectrophotometer (Kaiao, Beijing, China) was used to measure RNA purity. The RNA integrity and concentration were measured using the RNA Nano 6000 Assay Kit (Agilent Technologies, Santa Clara, CA, USA) on the Bioanalyzer 2100 system (Agilent Technologies) (the RNA mass is shown in Supplementary Table S1). Ribosomal RNA was removed using Epicentre Ribo-Zero™ Gold Kits (Epicentre, Madison, WI, USA). Sequencing libraries were subsequently generated following the manufacturer’s recommendations with a varied index label by NEBNext® Ultra™ Directional RNA Library Prep Kit for Illumina (NEB, Beverly, MA, USA). After the preparation of the library, the insert size was assessed by the Agilent Bioanalyzer 2100 system (Agilent Technologies), and qualified using the Taqman fluorescence probe of the AB Step One Plus Real-Time PCR system. Then, RNA-Seq was performed on the Illumina HiSeq X Ten system (Illumina, San Diego, CA, USA) following the HiSeq X Ten User Guide and 150 bp paired-end reads were generated.

2.4. Differentially Expressed Genes

After filtering and quality control steps, we used STAR [14] version 020201 to align the clean reads to the reference genome (Rattus norvegicus: Rnor 6.0 version 86, Ensembl). Subsequently, we used Stringtie [15] version 1.3.0 to assemble the transcripts. The fragments per kilobase of exon per million fragments mapped  (FPKM) of all annotated genes were obtained by Ballgown [16] version 2.10.0. Finally, a Bayes-regularized t-test, with a false discovery rate (FDR) correction using Cyber-T bayesreg.R [17], was used to obtained the differentially expressed genes (DEGs). FDR < 0.05 was regarded as statistically significant.

2.5. Enrichment Analysis of Kyoto Encyclopedia of Genes and Genomes (KEGG) Pathways

The Database for Annotation, Visualization, and Integrated Discovery (DAVID) was used to obtain the Kyoto Encyclopedia of Genes and Genomes (KEGG). A Benjamini correction of P value < 0.05 was set as the threshold of significance.

2.6. Prediction of the Transcriptional Factors of Key Dysregulated Genes

To investigate the potential upstream regulators that affect gene expression in muscle of GK rats, we used the Match tool in TRANSFANC [18] and JASPAR [19] to predict the transcriptional factors (TFs) of key dysregulated genes. The searching length was 2 kb upstream of each gene locus. The TFs with Relative score < 0.9 predicted by JASPAR were removed. The TFs with core matches <0.9 and matrix matches <0.9 predicted by TRANSFAC were filtered. After filtering, the overlap of TFs in TRANSFAC and JASPAR were regarded as predicted TFs.

2.7. Validation of RNA-Seq Results with RT-qPCR

RT-qPCR analyses were conducted in the LightCycler 96 system (Roche, Basel, Switzerland) to verify the DEGs by SYBR® Premix Ex Taq™ II reagent (TaKaRa, Dalian, China). The methods of 2−ΔΔCt was used to quantify the relative mRNA levels, in which Wistar rats were used as the controls and β-actin (Actb) [20,21] was used as an internal reference.

2.8. Statistical Analysis

All data were expressed as mean ± SEM unless otherwise stated. Significant differences were assessed using a two-tailed student t-test, with P values less than 0.05 considered statistically significant.

3. Results

3.1. The Characterization of Rats

GK random plasma glucose concentration was significantly increased compared to the control at both three and four weeks (Figure 1, P < 0.001). Moreover, GK plasma glucose concentration at four weeks was significantly higher than that of GK at three weeks (Figure 1, P < 0.001). There was no significant difference in plasma insulin concentration between GK and Wistar rats at both three and four weeks (Supplementary Table S2). The weight of GK rats was similar with Wistar rats at three and four weeks (Supplementary Table S2).

3.2. The Differentially Expressed Genes and Enriched Pathways

In total, 12,300 genes were detected with expression (mean FPKM > 1) in GK and Wistar skeletal muscle at three and four weeks. This detected gene number was similar to that of a previous study [22]. Among the detected genes, and in comparison to the Wistar rats, 600 and 1785 differentially expressed genes (DEGs) were obtained from GK skeletal muscle at three and four weeks, respectively (FDR < 0.05, Supplementary Table S3). Of these DEGs, there were 380 overlap DEGs at both three and four weeks. We found 220 DEGs were specific at three weeks and 1405 DEGs were specific at four weeks (Figure 2A). Additionally, 273 and 1027 DEGs were upregulated at three and four weeks, respectively, while 327 and 758 DEGs were downregulated at three and four weeks, respectively (Supplementary Table S3). A KEGG signaling pathway enrichment analysis with DEGs was conducted using DAVID. There was no enriched pathway at three weeks (Benjamini correction of P < 0.05). In contrast, six KEGG pathways were enriched at four weeks (Figure 2B, Benjamini correction of P < 0.05) and the PI3K-Akt signaling pathway was ranked first with the lowest Benjamini corrected P value. These results indicate the dynamic diabetic processes of GK rats from three to four weeks. It has been reported that mice with skeletal muscle which specifically overexpresses Nr4a1 (orphan nuclear receptor 4A1, also known as Nur77) have enhanced fatty acid oxidation [23]. This suggests that upregulation of the gene Nr4a1 in GK skeletal muscle at four weeks—which was enriched in the PI3K-Akt signaling pathway—may facilitate fatty acid oxidation in GK skeletal muscle compared to that of age-matched Wistar rats.

3.3. The Differentially Expressed Genes Related to Glucose Uptake and Oxidation

Since glucose transport is the rate-limiting step for glucose uptake and utilization in muscle [24,25,26], and glucose transport, uptake and oxidation are potentially implicated in the regulation of glucose homeostasis, we focused on the genes related to glucose transport, uptake and oxidation.
There was no significant difference in the expression of Slc2a4 (glucose transporter 4, also termed as Glut4) between GK and Wistar skeletal muscle at three or four weeks (Figure 3A,B). The FPKM of gene Tbc1d4 (TBC1 domain family, member 4, also termed AS160) was significantly reduced in GK skeletal muscle compared to that of age-matched Wistar rats at both three and four weeks (FDR < 0.05 and FDR < 0.001, respectively) (Figure 3A,B). The relative expression of Tbc1d4 was verified as downregulated without significant differences by RT-qPCR (Figure 3C).
The expression of Pdk4 (pyruvate dehydrogenase kinase 4) was significantly increased in GK skeletal muscle compared to that of age-matched Wistar rats at three and four weeks (FDR < 0.001) (Figure 3A and 3B), which was verified as upregulated by RT-qPCR (Figure 3C) with significant difference (Figure 3C). The expression of two confirmed Pdk4 transcription factors, Cebpb (CCAAT/enhancer-binding protein β, C/EBPβ) and Foxo1 (forkhead box O1), were increased in GK skeletal muscle in comparison to age-matched Wistar rats at three and four weeks (Figure 3), indicating that the two Pdk4 transcription factors, Cebpb and Foxo1, may be responsible for Pdk4 upregulation.
As glycolysis is an important pathway for glucose oxidation, we also focused on the genes related to glycolysis. The expression of genes related to glycolysis, including Hk2 (hexokinase 2), Pfkm (phosphofructokinase, muscle), and Pkm (pyruvate kinase) showed no significant difference between GK and Wistar rats (Supplementary Table S4).

3.4. The Differentially Expressed Genes Related to Fatty Acid Oxidation

As Randle et al. pointed out, increased oxidation of fatty acids could repress glucose oxidation [27,28]. Genes (Acadl, Acsl1, Fabp4 and Cpt1a) related to fatty acid oxidation may be implicated in the regulation of glucose oxidation. Among those genes, the expression of Acadl (acyl-CoA dehydrogenase, long chain), Acsl1 (acyl-CoA synthetase long-chain family member 1) and Fabp4 (fatty acid binding protein 4) showed no significant difference in GK skeletal muscle at three weeks compared to that of age-matched Wistar rats (Figure 4A,B), and were significantly increased at four weeks, as verified by RT-qPCR (Figure 4C). Cpt1a (carnitine palmitoyltransferase 1A), encoded by Cpt1a and responsible for transportation of long-chain fatty acids into mitochondria, was significantly increased in GK skeletal muscle in comparison to age-matched Wistar rats at three and four weeks (mean FPKM < 6, Supplementary Table S4). Other genes involved in the oxidation of fatty acids did not show significant differences between the GK and Wistar rat groups (Supplementary Table S4).

3.5. The Predicted Transcriptional Factors of Key Dysregulated Genes

To investigate the factors affecting the expression of key dysregulated genes, we predicted the transcriptional factor (TFs) of genes Acadl, Acsl1, Fabp4, Pdk4 and Tbcld4. The TFs predicted by TRANSFAC and JASPAR are shown in Supplementary Table S7. There were 12, 9, 17, 16 and 16 overlapping TFs of Acadl, Acsl1, Fabp4, Pdk4 and Tbcld4, respectively, according to TRANSFAC and JASPAR (Supplementary Table S7). The overlap TFs with dysregulated expression at three weeks and four weeks were Jun (proto-oncogene, AP-1 transcription factor subunit), Myc (proto-oncogene, bHLH transcription factor), Sox9 (SRY box 9), Srebf1 (sterol regulatory element binding transcription factor 1) and Junb (JunB proto-oncogene, AP-1 transcription factor subunit). Jun, which was the predicted TF of Acsl1, Fabp4 and Tbc1d4, was significantly downregulated in GK skeletal muscle compared to that of age-matched Wistar rats at three weeks, but was significantly upregulated in GK skeletal muscle at four weeks compared to age-matched Wistar rats (Figure 5). Myc was the predicted TF of Acsl1 and was significantly upregulated in GK skeletal muscle compared to that of age-matched Wistar rats at four weeks, but with no significant difference at three weeks. Sox9, which was the predicted TF of Fabp4 and Pdk4, was significantly upregulated in GK skeletal muscle compared to Wistar rats at four weeks, but was without significant difference at three weeks. Srebf1 was significantly downregulated in GK skeletal muscle compared to that of age-matched Wistar rats at three and four weeks, and was the predicted TF of Pdk4. Junb was the predicted TF of Tbc1d4 and was significantly upregulated in GK skeletal muscle compared to that of age-matched Wistar rats at three and four weeks. Therefore, dysregulation of the predicted TFs may lead to the dysregulated expression of Acadl, Acsl1, Fabp4, Pdk4 and Tbcld4, although further study is needed to verify this.

4. Discussion

It is known that postprandial hyperglycemia is important in whole-body glycemic control in T2D patients [29,30]. In addition, postprandial hyperglycemia is one of the early events in T2D patients [31]. As such, it is essential to reduce postprandial hyperglycemia of T2D patients [32]. Since GK rats are a suitable model for investigating the mechanism of postprandial hyperglycemia in skeletal muscle, we performed transcriptome analysis on the skeletal muscle of GK and Wistar (control) rats by RNA-Seq. From this analysis, 600 and 1785 DEGs were obtained from the skeletal muscle of GK rats at three and four weeks, respectively, compared to the Wistar rats. In considering previous studies, we predicted that dysregulated expression of key dysregulated genes (Acadl, Acsl1, Fabp4, Pdk4 and Tbcld4) could lead to decreased glucose uptake and oxidation, as well as enhanced fatty acid oxidation, in GK skeletal muscle at three and four weeks. These dysregulated metabolic processes may lead to the accumulation of glucose in the blood of GK individuals, which in turn may result in postprandial hyperglycemia in GK rats during an early stage (three and four weeks).
Glucose uptake is a complex process that involves a number of genes. Among these genes, Slc2a4 encodes the transporter protein Glut4. Glut4 mediates postprandial glucose transport in skeletal muscle [33,34], and is also necessary for glucose homeostasis in muscles [35,36,37,38]. In fact, muscle-specific Slc2a4 knockout mice have been shown to exhibit significantly reduced glucose uptake in their muscles [39,40,41]. It has also been reported that Tbc1d4, encoded by the gene Tbc1d4, is responsible for regulating glucose transport [42,43]. Whole-body Tbc1d4 knockout mice have exhibited significantly decreased insulin-stimulated glucose uptake in their skeletal muscle [44,45]. The overexpression of Pdk4 could decrease glucose uptake, and knockdown of Pdk4 has been shown to increase glucose uptake in neonatal myocytes [46]. Glycogen synthesis is also implicated in glucose disposal in skeletal muscle in the postprandial state [47], which may influence the homeostasis of glucose. Gsk3b (glycogen synthase kinase 3 beta) and Gys1 (glycogen synthase 1) play key roles in glycogen synthesis [48,49]. It has been reported that insulin-stimulated glucose uptake is markedly decreased in the skeletal muscle of GK rats weighing 280–300 g [50,51,52]. In addition, a study on Pdk4 expression has shown that it is significantly higher in GK skeletal muscle from 4–12 weeks [53], which is consistent with our results. In the present study, there was no expression difference in Slc2a4 (Figure 3), Gsk3b or Gys1 (Supplementary Table S4) in GK and Wistar skeletal muscle at three and four weeks. Therefore, the significantly decreased Tbc1d4 expression and the significantly increased Pdk4 expression in GK skeletal muscle—compared to that of age-matched Wistar rats—may have led to the observed decrease in glucose uptake of GK skeletal muscle at three and four weeks. Further study is needed to confirm the reduced glucose uptake of GK skeletal muscle at three and four weeks.
It is well known that there is metabolic competition between fatty acid oxidation and glucose oxidation, wherein increased fatty acid oxidation can inhibit glucose oxidation in heart and skeletal muscle [27,28]. The genes (Acadl, Acsl1 and Fabp4) related to fatty acid oxidation are also involved in regulating glucose oxidation, which would therefore affect the glucose content of blood [12,54]. Herein, the expression of genes (Acadl, Acsl1, Fabp4 and Cpt1a) related to fatty acid oxidation have exhibited an increase in GK skeletal muscle compared to that of age-matched Wistar rats at three weeks, and even more notably at four weeks. It has been reported that muscle-specific Acadl knockout mice exhibit reduced fatty acid oxidation [55]. Furthermore, skeletal muscle-specific deficiency Acsl1 mice have been shown to exhibit a significant decrease in fatty acid oxidation and lower plasma glucose concentration [54]. The gene product of Fabp4 has been shown to band and transport intracellular long-chain fatty acids. Whole-body Fabp4 knockout mice exhibit a reduction in free fatty acid release and a reduced plasma glucose level [12,13]. It has been suggested that overexpression of Cpt1a leads to enhanced fatty acid oxidation in hepatocytes, β-cells and muscle cells [56,57,58,59]. Besides regulating glucose uptake, Pdk4 is also implicated in regulating oxidation of glucose and fatty acids, with upregulated expression of Pdk4 reducing glucose oxidation and enhancing fatty acid oxidation in myocardium and skeletal muscle [60]. It has been demonstrated that skeletal muscle-specific overexpression of Nr4a1 facilitates fatty acid oxidation in the skeletal muscle of mice [23]. Reduced glucose oxidation in the skeletal muscle of GK rats aged 8–10 weeks has also been reported [51,52]. Moreover, it has been reported that the muscle lipid contents in GK rats is lower than that of Wistar rats aged four weeks and seven months [9,61]. Thus, the increased expression of key dysregulated genes (Acadl, Acsl1, Fabp4 and Pdk4) in GK skeletal muscle—in comparison to the age-matched Wistar rats—may have resulted in the enhanced fatty acid oxidation and reduced glucose oxidation in the GK skeletal muscle at three and four weeks. Additionally, the increased Nr4a1 expression in GK skeletal muscle compared to that of Wistar rats may be associated with the enhanced fatty acid oxidation and reduced glucose oxidation in GK skeletal muscle at four weeks. Further study is needed to confirm the increase in fatty acid oxidation and reduction in glucose oxidation in the skeletal muscle of GK rats aged three and four weeks.
By integrating previous studies and the results of our analysis, a putative mechanism of postprandial hyperglycemia in GK skeletal muscle during the early stage of three and four weeks is proposed in Figure 6. The putative mechanism suggests that the dysregulated expression of certain genes (Tbcld4, Pdk4, Acadl, Acsl1, and Fabp4) may result in increased fatty acid oxidation and reduced glucose uptake and oxidation in GK skeletal muscle during the early stage of three and four weeks. This may result in the accumulation of glucose in the blood of GK individuals and, in turn, postprandial hyperglycemia in those rats.

5. Conclusions

In the present study, we revealed transcriptome changes in GK skeletal muscle during the early stages of three and four weeks and found that the expression of Tbc1d4, Pdk4, Acadl, Acsl1 and Fabp4 were dysregulated. This is the first study that suggests Tbc1d4, Acadl, Acsl1 and Fabp4 are associated with early diabetes in GK rats, which may provide a new scope for pathogenesis of postprandial hyperglycemia and controlling postprandial glucose levels in T2D patients at early stages.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4425/10/6/406/s1, Table S1 The information of RNA mass, Table S2 The insulin concentration and body weight of rats. (Data are shown as mean ± SEM. n = 10), Table S3 The differentially expressed genes in GK rats compared with Wistar rats at three and four weeks, Table S4 The details of genes related to glycolysis, fatty acid oxidation and glycogen synthesis at three and four weeks, Table S5 The TFs predicted by TRANSFAC and JASPAR.

Author Contributions

H.D. conceived and designed the study; W.Z. and Y.M. carried out the data analysis and wrote the manuscript; W.Z. and S.F. performed the experiments; Y.M., S.F., X.L., Z.C., L.H., and H.D. revised the manuscript. All authors approved the final manuscript.

Funding

This work was supported by the National Key R&D Program of China (2018YFC0910201, 2017YFD0400300), the Key R&D Program of Guangdong Province (2019B020226001), and the Science and the Technology Planning Project of Guangzhou (201704020176).

Conflicts of Interest

The authors declare that they have no competing interests.

Data Availability

The RNA-Seq data used to support this study is uploaded to NCBI, the SRA accession is PRJNA540098.

References

  1. Cho, N.; Shaw, J.; Karuranga, S.; Huang, Y.; Fernandes, J.D.R.; Ohlrogge, A.; Malanda, B. IDF Diabetes Atlas: Global estimates of diabetes prevalence for 2017 and projections for 2045. Diabetes Res. Clin. Pract. 2018, 138, 271–281. [Google Scholar]
  2. Report of the Expert Committee on the Diagnosis and Classification of Diabetes Mellitus. Diabetes. Care 1997, 20, 1183–1197. [CrossRef]
  3. Qiao, Q.; Nakagami, T.; Tuomilehto, J.; Borch-Johnsen, K.; Balkau, B.; Iwamoto, Y.; Tajima, N.; Group International Diabetes Epidemiology; Decoda Study Group. Comparison of the fasting and the 2-h glucose criteria for diabetes in different Asian cohorts. Diabetologia 2000, 43, 1470–1475. [Google Scholar] [CrossRef] [Green Version]
  4. Yang, W.; Lu, J.; Weng, J.; Jia, W.; Ji, L.; Xiao, J.; Shan, Z.; Liu, J.; Tian, H.; Ji, Q. Prevalence of diabetes amongmen and women in China. N. Engl. J. Med. 2010, 362, 1090–1101. [Google Scholar] [CrossRef]
  5. Qiao, Q.; Hu, G.; Tuomilehto, J.; Nakagami, T.; Balkau, B.; Borch-Johnsen, K.; Ramachandran, A.; Mohan, V.; Iyer, S.R.; Tominaga, M.; et al. Age- and sex-specific prevalence of diabetes and impaired glucose regulation in 11 Asian cohorts. Diabetes Care 2003, 26, 1770–1780. [Google Scholar]
  6. GoTo, Y.; Kakizaki, M.; Masaki, N. Spontaneous diabetes produced by selective breeding of normal Wistar rats. Proc. Japan Acad. 1975, 51, 80–85. [Google Scholar] [CrossRef]
  7. Goto, Y.; Kakizaki, M.; Masaki, N. Production of Spontaneous Diabetic Rats by Repetition of Selective Breeding. Tohoku J. Exp. Med. 1976, 119, 85–90. [Google Scholar] [CrossRef] [Green Version]
  8. Movassat, J.; Bailbe, D.; Lubrano-Berthelier, C.; Picarel-Blanchot, F.; Bertin, E.; Mourot, J.; Portha, B. Follow-up of GK rats during prediabetes highlights increased insulin action and fat deposition despite low insulin secretion. Am. J. Physiol. Metab. 2008, 294, E168–E175. [Google Scholar] [CrossRef]
  9. Nie, J.; Xue, B.; Sukumaran, S.; Jusko, W.J.; Dubois, D.C.; Almon, R.R. Differential muscle gene expression as a function of disease progression in Goto-Kakizaki diabetic rats. Mol. Cell. Endocrinol. 2011, 338, 10–17. [Google Scholar] [CrossRef] [Green Version]
  10. DeFronzo, R.A.; Jacot, E.; Jéquier, E.; Maeder, E.; Wahren, J.; Felber, J.P. The Effect of Insulin on the Disposal of Intravenous Glucose: Results from Indirect Calorimetry and Hepatic and Femoral Venous Catheterization. Diabetes 1981, 30, 1000–1007. [Google Scholar] [CrossRef]
  11. Baron, A.D.; Brechtel, G.; Wallace, P.; Edelman, S.V. Rates and tissue sites of non-insulin- and insulin-mediated glucose uptake in humans. Am. J. Physiol. Metab. 1988, 255, 769–774. [Google Scholar] [CrossRef]
  12. Uysal, K.T.; Scheja, L.; Wiesbrock, S.M.; Bonner-Weir, S.; Hotamisligil, G.S. Improved Glucose and Lipid Metabolism in Genetically Obese Mice Lacking aP2. Endocrinology 2000, 141, 3388–3396. [Google Scholar] [CrossRef]
  13. Scheja, L.; Makowski, L.; Uysal, K.T.; Wiesbrock, S.M.; Shimshek, D.R.; Meyers, D.S.; Morgan, M.; Parker, R.A.; Hotamisligil, G.S. Altered insulin secretion associated with reduced lipolytic efficiency in aP2-/- mice. Diabetes 1999, 48, 1987–1994. [Google Scholar] [CrossRef]
  14. Dobin, A.; Gingeras, T.R. Mapping RNA-seq Reads with STAR. Curr. Protoc. Bioinformatics 2015, 51, 11–19. [Google Scholar] [PubMed]
  15. Pertea, M.; Pertea, G.M.; Antonescu, C.M.; Chang, T.-C.; Mendell, J.T.; Salzberg, S.L. StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat. Biotechnol. 2015, 33, 290–295. [Google Scholar] [CrossRef] [Green Version]
  16. Frazee, A.C.; Pertea, G.; E Jaffe, A.; Langmead, B.; Salzberg, S.L.; Leek, J.T. Ballgown bridges the gap between transcriptome assembly and expression analysis. Nat. Biotechnol. 2015, 33, 243–246. [Google Scholar] [CrossRef] [Green Version]
  17. Kayala, M.A.; Baldi, P. Cyber-T web server: differential analysis of high-throughput data. Nucleic Acids Res. 2012, 40, W553–W559. [Google Scholar] [CrossRef]
  18. Matys, V. TRANSFAC®: transcriptional regulation, from patterns to profiles. Nucleic Acids Res. 2003, 31, 374–378. [Google Scholar] [CrossRef]
  19. Khan, A.; Fornes, O.; Stigliani, A.; Gheorghe, M.; Castro-Mondragon, J.A.; van der Lee, R.; Bessy, A.; Cheneby, J.; Kulkarni, S.R.; Tan, G. JASPAR 2018: Update of the open-access database of transcription factor binding profiles and its web framework. Nucleic Acids Res. 2018, 46, D260–D266. [Google Scholar] [CrossRef]
  20. Qi, Z.; He, J.; Zhang, Y.; Shao, Y.; Ding, S. Exercise training attenuates oxidative stress and decreases p53 protein content in skeletal muscle of type 2 diabetic Goto-Kakizaki rats. Radic. Boil. Med. 2011, 50, 794–800. [Google Scholar] [CrossRef]
  21. Tachibana, N.; Yamashita, Y.; Nagata, M.; Wanezaki, S.; Ashida, H.; Horio, F.; Kohno, M. Soy β-conglycinin improves glucose uptake in skeletal muscle and ameliorates hepatic insulin resistance in Goto-Kakizaki rats. Nutr. Res. 2014, 34, 160–167. [Google Scholar] [CrossRef]
  22. Rubin, A.; Salzberg, A.C.; Imamura, Y.; Grivitishvilli, A.; Tombran-Tink, J. Identification of novel targets of diabetic nephropathy and PEDF peptide treatment using RNA-seq. BMC Genom. 2016, 17, 1047. [Google Scholar] [CrossRef]
  23. Chao, L.C.; Wroblewski, K.; Ilkayeva, O.R.; Stevens, R.D.; Bain, J.; Meyer, G.A.; Schenk, S.; Martinez, L.; Vergnes, L.; Narkar, V.A.; et al. Skeletal muscle Nur77 expression enhances oxidative metabolism and substrate utilization. J. Lipid Res. 2012, 53, 2610–2619. [Google Scholar] [CrossRef] [Green Version]
  24. Fink, R.; Wallace, P.; Brechtel, G.; Olefsky, J. Evidence that glucose transport is rate-limiting for in vivo glucose uptake. Metab. Clin. Exp. 1992, 41, 897–902. [Google Scholar] [CrossRef]
  25. Ziel, F.H.; Venkatesan, N.; Davidson, M.B. Glucose Transport Is Rate Limiting for Skeletal Muscle Glucose Metabolism in Normal and STZ-Induced Diabetic Rats. Diabetes 1988, 37, 885–890. [Google Scholar] [CrossRef]
  26. Ren, J.M.; A Marshall, B.; A Gulve, E.; Gao, J.; Johnson, D.W.; O Holloszy, J.; Mueckler, M. Evidence from transgenic mice that glucose transport is rate-limiting for glycogen deposition and glycolysis in skeletal muscle. J. Boil. Chem. 1993, 268, 16113–16115. [Google Scholar]
  27. Randle, P.; Garland, P.; Hales, C.; Newsholme, E. The glucose fatty-acid cycle. Its role in insulin sensitivity and the metabolic disturbances of diabetes mellitus. Lancet 1963, 281, 785–789. [Google Scholar] [CrossRef]
  28. Hue, L.; Taegtmeyer, H. The Randle cycle revisited: a new head for an old hat. Am. J. Physiol. Metab. 2009, 297, E578–E591. [Google Scholar] [CrossRef] [Green Version]
  29. Monnier, L.; Lapinski, H.; Colette, C. Contributions of Fasting and Postprandial Plasma Glucose Increments to the Overall Diurnal Hyperglycemia of Type 2 Diabetic Patients: Variations with increasing levels of HbA1c. Diabetes Care 2003, 26, 881–885. [Google Scholar] [CrossRef]
  30. Woerle, H.J.; Pimenta, W.P.; Meyer, C.; Gosmanov, N.R.; Szoke, E.; Szombathy, T.; Mitrakou, A.; Gerich, J.E. Diagnostic and therapeutic implications of relationships between fasting, 2-hour postchallenge plasma glucose and hemoglobin a1c values. Arch. Intern. Med. 2004, 164, 1627–1632. [Google Scholar] [CrossRef]
  31. Monnier, L.; Colette, C.; Dunseath, G.J.; Owens, D.R. The Loss of Postprandial Glycemic Control Precedes Stepwise Deterioration of Fasting With Worsening Diabetes. Diabetes Care 2007, 30, 263–269. [Google Scholar] [CrossRef] [Green Version]
  32. Woerle, H.J.; Neumann, C.; Zschau, S.; Tenner, S.; Irsigler, A.; Schirra, J.; Gerich, J.E.; Göke, B. Impact of fasting and postprandial glycemia on overall glycemic control in type 2 diabetes Importance of postprandial glycemia to achieve target HbA1c levels. Diabetes Res. Clin. Pract. 2007, 77, 280–285. [Google Scholar] [CrossRef]
  33. Birnbaum, M.J. Identification of a novel gene encoding an insulin-responsive glucose transporter protein. Cell 1989, 57, 305–315. [Google Scholar] [CrossRef]
  34. Kaestner, K.H.; Christy, R.J.; McLenithan, J.C.; Braiterman, L.T.; Cornelius, P.; Pekala, P.H.; Lane, M.D. Sequence, tissue distribution, and differential expression of mRNA for a putative insulin-responsive glucose transporter in mouse 3T3-L1 adipocytes. Proc. Acad. Sci. 1989, 86, 3150–3154. [Google Scholar] [CrossRef]
  35. Abel, E.D.; Peroni, O.; Kim, J.K.; Kim, Y.-B.; Boss, O.; Hadro, E.; Minnemann, T.; Shulman, G.I.; Kahn, B.B. Adipose-selective targeting of the GLUT4 gene impairs insulin action in muscle and liver. Nat. Cell Boil. 2001, 409, 729–733. [Google Scholar] [CrossRef]
  36. Dugani, C.B.; Klip, A. Glucose transporter 4: cycling, compartments and controversies. EMBO Rep. 2005, 6, 1137–1142. [Google Scholar] [CrossRef]
  37. Graham, T.; Kahn, B. Tissue-specific Alterations of Glucose Transport and Molecular Mechanisms of Intertissue Communication in Obesity and Type 2 Diabetes. Horm. Metab. 2007, 39, 717–721. [Google Scholar] [CrossRef] [Green Version]
  38. Huang, S.; Lifshitz, L.M.; Jones, C.; Bellve, K.D.; Standley, C.; Fonseca, S.; Corvera, S.; Fogarty, K.E.; Czech, M.P. Insulin Stimulates Membrane Fusion and GLUT4 Accumulation in Clathrin Coats on Adipocyte Plasma Membranes. Mol. Cell. Boil. 2007, 27, 3456–3469. [Google Scholar] [CrossRef] [Green Version]
  39. Kim, J.K.; Zisman, A.; Fillmore, J.J.; Peroni, O.D.; Kotani, K.; Perret, P.; Zong, H.; Dong, J.; Kahn, C.R.; Kahn, B.B.; et al. Glucose toxicity and the development of diabetes in mice with muscle-specific inactivation of GLUT4. J. Clin. Investig. 2001, 108, 153–160. [Google Scholar] [CrossRef]
  40. Zisman, A.; Peroni, O.D.; Abel, E.D.; Michael, M.D.; Mauvais-Jarvis, F.; Lowell, B.B.; Wojtaszewski, J.F.; Hirshman, M.F.; Virkamaki, A.; Goodyear, L.J.; et al. Targeted disruption of the glucose transporter 4 selectively in muscle causes insulin resistance and glucose intolerance. Nat. Med. 2000, 6, 924–928. [Google Scholar] [CrossRef]
  41. Kim, Y.-B.; Peroni, O.D.; Aschenbach, W.G.; Minokoshi, Y.; Kotani, K.; Zisman, A.; Kahn, C.R.; Goodyear, L.J.; Kahn, B.B. Muscle-Specific Deletion of the Glut4 Glucose Transporter Alters Multiple Regulatory Steps in Glycogen Metabolism. Mol. Cell. Boil. 2005, 25, 9713–9723. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Sano, H.; Kane, S.; Mı̂inea, C.P.; Asara, J.M.; Lane, W.S.; Garner, C.W.; Sano, E.; Miinea, C.P.; Lienhard, G.E. Insulin-stimulated Phosphorylation of a Rab GTPase-activating Protein Regulates GLUT4 Translocation. J. Boil. Chem. 2003, 278, 14599–14602. [Google Scholar] [CrossRef] [Green Version]
  43. Sakamoto, K.; Holman, G.D. Emerging role for AS160/TBC1D4 and TBC1D1 in the regulation of GLUT4 traffic. Am. J. Physiol. Metab. 2008, 295, E29–37. [Google Scholar] [CrossRef] [PubMed]
  44. Wang, H.Y.; Ducommun, S.; Quan, C.; Xie, B.; Li, M.; Wasserman, D.H.; Sakamoto, K.; Mackintosh, C.; Chen, S. AS160 deficiency causes whole-body insulin resistance via composite effects in multiple tissues. Biochem. J. 2013, 449, 479–489. [Google Scholar] [CrossRef] [PubMed]
  45. Lansey, M.N.; Walker, N.N.; Hargett, S.R.; Stevens, J.R.; Keller, S.R. Deletion of Rab GAP AS160 modifies glucose uptake and GLUT4 translocation in primary skeletal muscles and adipocytes and impairs glucose homeostasis. Am. J. Physiol. Metab. 2012, 303, E1273–E1286. [Google Scholar] [CrossRef] [Green Version]
  46. Qin, X.; Hou, Z.; Guo, Y.; Wu, J.; Gao, F.; Li, T.; Xu, J.; Liu, Z.; Zheng, H.; Zhang, X. Glucose oxidation positively regulates glucose uptake and improves cardiac function recovery after myocardial reperfusion. Am. J. Physiol. Metab. 2017, 313, E577–E585. [Google Scholar] [Green Version]
  47. Jue, T.; Stein, P.; Shulman, G.I.; Rothman, D.L.; DeFronzo, R.A.; Shulman, R.G. Quantitation of Muscle Glycogen Synthesis in Normal Subjects and Subjects with Non-Insulin-Dependent Diabetes by13C Nuclear Magnetic Resonance Spectroscopy. New Engl. J. Med. 1990, 322, 223–228. [Google Scholar]
  48. Macaulay, K.; Blair, A.S.; Hajduch, E.; Hundal, H.S.; Terashima, T.; Baba, O.; Sutherland, C. Constitutive Activation of GSK3 Down-regulates Glycogen Synthase Abundance and Glycogen Deposition in Rat Skeletal Muscle Cells. J. Boil. Chem. 2005, 280, 9509–9518. [Google Scholar] [CrossRef] [Green Version]
  49. Wang, Y.; Wang, Y.; Zhong, T.; Guo, J.; Li, L.; Zhang, H.; Wang, L. Transcriptional regulation of pig GYS1 gene by glycogen synthase kinase 3β (GSK3β). Mol. Cell. Biochem. 2017, 424, 203–208. [Google Scholar] [CrossRef]
  50. Song, X.M.; Kawano, Y.; Krook, A.; Ryder, J.W.; Efendic, S.; Roth, R.A.; Wallberg-Henriksson, H.; Zierath, J. Muscle fiber type-specific defects in insulin signal transduction to glucose transport in diabetic GK rats. Diabetes 1999, 48, 664–670. [Google Scholar] [CrossRef]
  51. Bitar, M.S.; Al-Saleh, E.; Al-Mulla, F. Oxidative stress—mediated alterations in glucose dynamics in a genetic animal model of type II diabetes. Life Sci. 2005, 77, 2552–2573. [Google Scholar] [CrossRef] [PubMed]
  52. Bitar, M.S.; Wahid, S.; Pilcher, C.W.; Al-Saleh, E.; Al-Mulla, F. α-lipoic acid mitigates insulin resistance in Goto-Kakizaki rats. Horm. Metab. Res. 2004, 36, 542–549. [Google Scholar] [CrossRef]
  53. Nie, J.; DuBois, D.C.; Jusko, W.J.; Almonm, R.R. Mechanistic population modeling of diabetes disease progression in Goto-Kakizaki rat muscle. Biopharm. Drug Dispos. 2011, 32, 50–63. [Google Scholar] [CrossRef]
  54. Li, L.O.; Grevengoed, T.J.; Paul, D.S.; Ilkayeva, O.; Koves, T.R.; Pascual, F.; Newgard, C.B.; Muoio, D.M.; Coleman, R.A. Compartmentalized acyl-CoA metabolism in skeletal muscle regulates systemic glucose homeostasis. Diabetes 2015, 64, 23–35. [Google Scholar] [CrossRef]
  55. Chegary, M.; te Brinke, H.; Ruiter, J.P.; Wijburg, F.A.; Stoll, M.S.; Minkler, P.E.; van Weeghel, M.; Schulz, H.; Hoppel, C.L.; Wanders, R.J. Mitochondrial long chain fatty acid β-oxidation in man and mouse. Biochim. Biophys. Acta 2009, 1791, 806–815. [Google Scholar] [CrossRef] [PubMed]
  56. Perdomo, G.; Commerford, S.R.; Richard, A.M.; Adams, S.H.; Corkey, B.E.; O’Doherty, R.M.; Brown, N.F. Increased β-oxidation in muscle cells enhances insulin-stimulated glucose metabolism and protects against fatty acid-induced insulin resistance despite intramyocellular lipid accumulation. J. Biol. Chem. 2004, 279, 27177–27186. [Google Scholar] [CrossRef] [PubMed]
  57. Herrero, L.; Rubí, B.; Sebastián, D.; Serra, D.; Asins, G.; Maechler, P.; Prentki, M.; Hegardt, F.G. Alteration of the malonyl-CoA/carnitine palmitoyltransferase I interaction in the β-cell impairs glucose-induced insulin secretion. Diabetes 2005, 54, 462–471. [Google Scholar] [CrossRef]
  58. Stefanovic-Racic, M.; Perdomo, G.; Mantell, B.S.; Sipula, I.J.; Brown, N.F.; O’Doherty, R.M. A moderate increase in carnitine palmitoyltransferase 1a activity is sufficient to substantially reduce hepatic triglyceride levels. Am. J. Physiol. Metab. 2008, 294, E969–E977. [Google Scholar] [CrossRef] [Green Version]
  59. Akkaoui, M.; Cohen, I.; Esnous, C.; Lenoir, V.; Sournac, M.; Girard, J.; Prip-Buus, C. Modulation of the hepatic malonyl-CoA-carnitine palmitoyltransferase 1A partnership creates a metabolic switch allowing oxidation of de novo fatty acids. Biochem. J. 2009, 420, 429–438. [Google Scholar] [CrossRef]
  60. Zhao, G.; Burgess, S.C.; Rosaaen-Stowe, K.A.; Inagaki, T.; Latif, S.; Shelton, J.M.; McAnally, J.; Bassel-Duby, R.; Harris, R.A.; Kliewer, S.A.; et al. Overexpression of pyruvate dehydrogenase kinase 4 in heart perturbs metabolism and exacerbates calcineurin-induced cardiomyopathy. Am. J. Physiol. Circ. Physiol. 2008, 294, 936–943. [Google Scholar] [CrossRef]
  61. Macia, M.; Pecchi, É.; Vilmen, C.; Desrois, M.; Lan, C.; Portha, B.; Bernard, M.; Bendahan, D.; Giannesini, B. Insulin Resistance Is Not Associated with an Impaired Mitochondrial Function in Contracting Gastrocnemius Muscle of Goto-Kakizaki Diabetic Rats In Vivo. PLOS ONE 2015, 10, e0129579. [Google Scholar]
Genes 10 00406 g001 550
Figure 1. The random plasma glucose concentration of Goto-Kakizaki (GK) and Wistar rats (data are shown as mean ± SEM. n = 10, *** P < 0.001).
Figure 1. The random plasma glucose concentration of Goto-Kakizaki (GK) and Wistar rats (data are shown as mean ± SEM. n = 10, *** P < 0.001).
Genes 10 00406 g001
Genes 10 00406 g002 550
Figure 2. The differentially expressed genes (DEGs) and enriched pathways in GK skeletal muscle. (A) Venn diagram of differentially expressed genes in GK skeletal muscle compared to age-matched Wistar rats at three and four weeks of age. (B) The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways enriched in GK skeletal muscle at four weeks (Benjamini correction of P < 0.05). Adjp represents the Benjamini correction of P values.
Figure 2. The differentially expressed genes (DEGs) and enriched pathways in GK skeletal muscle. (A) Venn diagram of differentially expressed genes in GK skeletal muscle compared to age-matched Wistar rats at three and four weeks of age. (B) The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways enriched in GK skeletal muscle at four weeks (Benjamini correction of P < 0.05). Adjp represents the Benjamini correction of P values.
Genes 10 00406 g002
Genes 10 00406 g003 550
Figure 3. The differentially expressed genes related to glucose uptake and oxidation. (A) The mean fragments per kilobase of exon per million fragments mapped  (FPKM) of Pdk4, Slc2a4, Tbc1d4, Cebpb and Foxo1 in skeletal muscle of GK and Wistar rats. (B) The fold change (meanE/meanC) and corresponding false discovery rate (FDR) values of Pdk4, Slc2a4, Tbc1d4, Cebpb and Foxo1 in GK rats in comparison with those of Wistar rats. (C) The fold change (2−ΔΔCt) of GK genes Pdk4, Tbc1d4, Cebpb and Foxo1 expression compared with Wistar rats by RT-qPCR (data are shown as mean ± SEM. n = 6, * P < 0.05, *** P < 0.001). MeanC represents the mean FPKM of Wistar rats, while meanE represents the mean FPKM of GK rats.
Figure 3. The differentially expressed genes related to glucose uptake and oxidation. (A) The mean fragments per kilobase of exon per million fragments mapped  (FPKM) of Pdk4, Slc2a4, Tbc1d4, Cebpb and Foxo1 in skeletal muscle of GK and Wistar rats. (B) The fold change (meanE/meanC) and corresponding false discovery rate (FDR) values of Pdk4, Slc2a4, Tbc1d4, Cebpb and Foxo1 in GK rats in comparison with those of Wistar rats. (C) The fold change (2−ΔΔCt) of GK genes Pdk4, Tbc1d4, Cebpb and Foxo1 expression compared with Wistar rats by RT-qPCR (data are shown as mean ± SEM. n = 6, * P < 0.05, *** P < 0.001). MeanC represents the mean FPKM of Wistar rats, while meanE represents the mean FPKM of GK rats.
Genes 10 00406 g003
Genes 10 00406 g004 550
Figure 4. The differentially expressed genes related to fatty acid oxidation. (A) The mean FPKM of Acadl, Acsl1 and Fabp4 in skeletal muscle of GK and Wistar rats. (B) The fold change (meanE /meanC) and corresponding FDR values of Acadl, Acsl1 and Fabp4 in GK rats in comparison with those of Wistar rats. (C) The fold change (2−ΔΔCt) of GK genes Acadl, Acsl1 and Fabp4 expression compared with Wistar rats by RT-qPCR (data are shown as mean ± SEM. n = 6, * P < 0.05, *** P < 0.001). MeanC represents the mean FPKM of Wistar rats, while meanE represents the mean FPKM of GK rats.
Figure 4. The differentially expressed genes related to fatty acid oxidation. (A) The mean FPKM of Acadl, Acsl1 and Fabp4 in skeletal muscle of GK and Wistar rats. (B) The fold change (meanE /meanC) and corresponding FDR values of Acadl, Acsl1 and Fabp4 in GK rats in comparison with those of Wistar rats. (C) The fold change (2−ΔΔCt) of GK genes Acadl, Acsl1 and Fabp4 expression compared with Wistar rats by RT-qPCR (data are shown as mean ± SEM. n = 6, * P < 0.05, *** P < 0.001). MeanC represents the mean FPKM of Wistar rats, while meanE represents the mean FPKM of GK rats.
Genes 10 00406 g004
Genes 10 00406 g005 550
Figure 5. The network of Acadl, Acsl1, Fabp4, Pdk4 and Tbc1d4 and their dysregulated transcription factors (TFs). (A) The network of Acadl, Acsl1, Fabp4, Pdk4 and Tbc1d4 and their dysregulated TFs at three weeks. (B) The network of Acadl, Acsl1, Fabp4, Pdk4 and Tbc1d4 and their dysregulated TFs at four weeks. Grey represents the genes expressed without significant difference; cyan represents the downregulated expression of genes; and red represents the upregulated expression of genes. The circled genes represent the target genes of the TFs, while the diamonds represent the individual TFs.
Figure 5. The network of Acadl, Acsl1, Fabp4, Pdk4 and Tbc1d4 and their dysregulated transcription factors (TFs). (A) The network of Acadl, Acsl1, Fabp4, Pdk4 and Tbc1d4 and their dysregulated TFs at three weeks. (B) The network of Acadl, Acsl1, Fabp4, Pdk4 and Tbc1d4 and their dysregulated TFs at four weeks. Grey represents the genes expressed without significant difference; cyan represents the downregulated expression of genes; and red represents the upregulated expression of genes. The circled genes represent the target genes of the TFs, while the diamonds represent the individual TFs.
Genes 10 00406 g005
Genes 10 00406 g006 550
Figure 6. The putative mechanism of postprandial hyperglycemia in GK skeletal muscle during the early stage of three and four weeks. The downregulated expression of Tbc1d4 would reduce the uptake of glucose. The upregulated expression of Acadl, Acsl1 and Fabp4 could contribute to enhanced fatty acid oxidation, which could repress the oxidation of glucose. Moreover, increased Pdk4 expression may not only reduce glucose uptake, but could also reduce glucose oxidation and enhance fatty acid oxidation. Reduced glucose uptake and oxidation, as well as enhanced fatty acid oxidation, may result in postprandial hyperglycemia. The red boxes represent upregulated expression, while the blue box represents downregulated expression. The upward red arrow represents an increase, while the downward blue arrow represents a reduction. The dotted lines indicate that the precise role is unclear.
Figure 6. The putative mechanism of postprandial hyperglycemia in GK skeletal muscle during the early stage of three and four weeks. The downregulated expression of Tbc1d4 would reduce the uptake of glucose. The upregulated expression of Acadl, Acsl1 and Fabp4 could contribute to enhanced fatty acid oxidation, which could repress the oxidation of glucose. Moreover, increased Pdk4 expression may not only reduce glucose uptake, but could also reduce glucose oxidation and enhance fatty acid oxidation. Reduced glucose uptake and oxidation, as well as enhanced fatty acid oxidation, may result in postprandial hyperglycemia. The red boxes represent upregulated expression, while the blue box represents downregulated expression. The upward red arrow represents an increase, while the downward blue arrow represents a reduction. The dotted lines indicate that the precise role is unclear.
Genes 10 00406 g006

Share and Cite

MDPI and ACS Style

Zhang, W.; Meng, Y.; Fu, S.; Li, X.; Chen, Z.; Huang, L.; Du, H. Transcriptome Changes of Skeletal Muscle RNA-Seq Speculates the Mechanism of Postprandial Hyperglycemia in Diabetic Goto-Kakizaki Rats During the Early Stage of T2D. Genes 2019, 10, 406. https://doi.org/10.3390/genes10060406

AMA Style

Zhang W, Meng Y, Fu S, Li X, Chen Z, Huang L, Du H. Transcriptome Changes of Skeletal Muscle RNA-Seq Speculates the Mechanism of Postprandial Hyperglycemia in Diabetic Goto-Kakizaki Rats During the Early Stage of T2D. Genes. 2019; 10(6):406. https://doi.org/10.3390/genes10060406

Chicago/Turabian Style

Zhang, Wenlu, Yuhuan Meng, Shuying Fu, Xingsong Li, Zixi Chen, Lizhen Huang, and Hongli Du. 2019. "Transcriptome Changes of Skeletal Muscle RNA-Seq Speculates the Mechanism of Postprandial Hyperglycemia in Diabetic Goto-Kakizaki Rats During the Early Stage of T2D" Genes 10, no. 6: 406. https://doi.org/10.3390/genes10060406

APA Style

Zhang, W., Meng, Y., Fu, S., Li, X., Chen, Z., Huang, L., & Du, H. (2019). Transcriptome Changes of Skeletal Muscle RNA-Seq Speculates the Mechanism of Postprandial Hyperglycemia in Diabetic Goto-Kakizaki Rats During the Early Stage of T2D. Genes, 10(6), 406. https://doi.org/10.3390/genes10060406

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop