Figures
Abstract
Background
Berberine is used to treat diabetes and dyslipidemia. However, the effect of berberine on specific diabetes treatment targets is unknown. In the current study, we investigated the effect of berberine on the random plasma glucose, glycated hemoglobin (HbA1C), AST, ALT, BUN and CREA levels of Zucker diabetic fatty (ZDF) rats, and we identified and verified the importance of potential therapeutic target genes to provide molecular information for further investigation of the mechanisms underlying the anti-diabetic effects of berberine.
Methods
ZDF rats were randomly divided into control (Con), diabetic (DM) and berberine-treated (300 mg⋅kg−1, BBR) groups. After the ZDF rats were treated with BBR for 12 weeks, its effect on the random plasma glucose and HbA1C levels was evaluated. Aspartate aminotransferase (AST), alanine aminotransferase (ALT), blood urea nitrogen (BUN), CREA and OGTT were measured from blood, respectively. The levels of gene expression in liver samples were analyzed using an Agilent rat gene expression 4x44K microarray. The differentially expressed genes (DEGs) were screened as those with log2 (Con vs DM) ≥ 1 and log2 (BBR vs DM) ≥ 1 expression levels, which were the genes with up-regulated expression, and those with log2 (Con vs DM) ≤ -1 and log2 (BBR vs DM) ≤ -1 expression levels, which were the genes with down-regulated expression; the changes in gene expression were considered significant at P<0.05. The functions of the DEGs were determined using gene ontology (GO) and pathway analysis. Furthermore, a protein-protein interaction (PPI) network was constructed using STRING and Cytoscape software. The expression levels of the key node genes in the livers of the ZDF rats were also analyzed using qRT-PCR.
Results
We found that 12 weeks of berberine treatment significantly decreased the random plasma glucose, HbA1C levels and improved glucose tolerance. There was a tendency for berberine to reduce AST, ALT, BUN except increase CREA levels. In the livers of the BBR group, we found 154 DEGs, including 91 genes with up-regulated expression and 63 genes with down-regulated expression. In addition, GO enrichment analysis showed significant enrichment of the DEGs in the following categories: metabolic process, localization, cellular process, biological regulation and response to stimulus process. After the gene screening, KEGG pathway analysis showed that the target genes are involved in multiple pathways, including the lysine degradation, glycosaminoglycan biosynthesis-chondroitin sulfate/dermatan sulfate and pyruvate metabolism pathways. By combining the results of PPI network and KEGG pathway analyses, we identified seven key node genes. The qRT-PCR results confirmed that the expression of the RHOA, MAPK4 and DLAT genes was significantly down-regulated compared with the levels in DM group, whereas the expression of the SgK494, DOT1L, SETD2 and ME3 genes was significantly up-regulated in the BBR group.
Conclusion
Berberine can significantly improve glucose metabolism and has a protective effects of liver and kidney function in ZDF rats. The qRT-PCR results for the crucial DEGs validated the microarray results. These results suggested that the RHOA, MAPK4, SGK494, DOT1L, SETD2, ME3 and DLAT genes are potential therapeutic target genes for the treatment of diabetes.
Citation: Wu YS, Chen Y-T, Bao Y-T, Li Z-M, Zhou X-J, He J-N, et al. (2016) Identification and Verification of Potential Therapeutic Target Genes in Berberine-Treated Zucker Diabetic Fatty Rats through Bioinformatics Analysis. PLoS ONE 11(11): e0166378. https://doi.org/10.1371/journal.pone.0166378
Editor: Michael Bader, Max Delbruck Centrum fur Molekulare Medizin Berlin Buch, GERMANY
Received: September 5, 2016; Accepted: October 27, 2016; Published: November 15, 2016
Copyright: © 2016 Wu et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: All microarray data files are available from ArrayExpress database (accession number: E-MTAB-5178).
Funding: This work was supported by the Project of National Great New Drug Research and Development (No.2012ZX09503001-001), and Traditional Chinese Medicine Open funds of Zhejiang Chinese Medical University (No. 752223A00201/005/019). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Diabetes mellitus (DM), which is characterized by disordered lipid metabolism and a high circulating blood glucose level, is a metabolic syndrome [1]. The main form of diabetes found in 90~95% of diabetic patients globally is type 2 diabetes mellitus (T2DM) [2]. Diabetes can cause many complications, including neuropathy, nephropathy and retinopathy [3]. Diabetes seriously affects the quality of life of patients.
The anti-diabetes activity of Rhizoma Coptidis was documented in the book “Notes of Elite Physicians” written by Hongjing Tao 1500 years ago [3, 4]. Berberine (BBR) is the major isoquinoline alkaloid constituent of the Chinese herb Rhizoma Coptidis [1,5,6,7], which has the beneficial characteristic of regulating glucose and lipid metabolism and has been extensively used in the treatment of obesity, diabetes and hypercholesterolemia [3]. Zucker diabetic fatty (ZDF) rats are an animal model of type 2 diabetes mellitus due to a leptin receptor gene mutation. [8]. However, the exact anti-diabetic gene targets of berberine in diabetic model rats have not yet been fully elucidated.
Thus, in this study, we examined the hypoglycemic activity of berberine in ZDF rats and showed that berberine treatment changed the levels of gene expression in the liver to improve the diabetic metabolism. To identify the exact gene targets, we screened the genes, constructed a protein-protein interaction network, and validated the expression levels of the crucially differentially expressed genes (DEGs) that were identified by protein-protein interaction network analysis using qRT-PCR.
Methods
Groups and treatment
Six-week-old male ZDF rats (housed two per cage) were obtained from the Animal Center of Zhejiang Chinese Medicine University. To induce diabetes symptoms, ZDFfa/fa rats were fed high-sugar/high-fat feed for 7 weeks, whereas ZDFfa/+ rats were fed the basic feed during this period. Afterward, the ZDFfa/fa rats were randomly divided into two groups that had a similar degree of hyperglycemia: the model group (DM) and the BBR group (n = 6 per group). The ZDFfa/+ rats served as the control (Con) group. The start and end of gavage, the mental and dietary status of animals were observed to monitor animal health in every day during the experimental procedure. The rats in the control and model groups were administered 0.5% gum tragacanth (CAS: 9000-65-1, Shanghai Macklin Biochemical Co., Ltd. Shanghai, China) in normal saline via gavage, whereas the rats in the BBR group were administered berberine (K151212, Xi’an Kai Lai Biological Engineering Co., Ltd. Xi’an, China) at 300 mg/kg dissolved in 0.5% gum tragacanth in normal saline via gavage. All of the animals were housed at 25°C with 50~70% humidity and a 12 h light/12 h dark cycle with free access to regular food and water. Before treatment, the first measurements (0 weeks) of the random plasma glucose, HbA1C, AST, ALT, BUN and CREA levels were taken. The random plasma glucose level was measured at the 1st, 3rd, 5th, 8th, 10th, and 12th weeks and the HbA1C, AST, ALT, BUN and CREA levels were measured at the 12th week. Throughout the experimental period, no animals died prior to the experimental endpoint. Blood samples were collected from heart after pentobarbital sodium (45mg/kg, i.p.) anesthesia, and then rapidly excised liver tissues, which were quickly frozen in liquid nitrogen, and stored at -80°C for use in the microarray and quantitative real-time PCR assays. At the end of the experiments, euthanasia occurred under sodium pentobarbital anesthesia followed by cardiac puncture/livers removal for all animals. This study was approved by the animal care and welfare committee of Zhejiang Chinese Medical University (permit number: ZSLL-20130-106). All of the procedures involving animals were performed under sodium pentobarbital anesthesia to minimize their suffering.
Measurement of Blood Biochemical Parameters
Random plasma glucose levels in ZDF rats after berberine administration.
Sera were prepared from the blood collected at 0 weeks and at weeks 1, 3, 5, 8, 10, and 12 of the treatment via tail vein puncture. The random plasma glucose levels were measured using accu-chek® Performa (Roche Co., Ltd. Shanghai. China).
Glycosylated hemoglobin (HbA1C) in ZDF rats after berberine administration.
Blood was collected at 0 weeks and at the 12th week of treatment via tail vein puncture. The levels of glycosylated hemoglobin (HbA1C) were determined using a glycosylated hemoglobin analyzer (Quotient Diagnostics Ltd.).
AST, ALT, BUN and CREA levels in ZDF rats after berberine administration.
Blood samples were obtained by puncturing retro-orbital plexus at 0 weeks and collected from heart under anesthesia at 12th week. The serum samples were obtained by centrifugation (1000g × 10 min, 4°C) for the measurement of blood biochemical parameters. AST, ALT, BUN and CREA levels were measured by a 7020 full automatic biochemical analyzer (Hitachi, Japan).
Glucose tolerance test.
Upon completion of 12 h fasting, before the start as 0 weeks and end as 12th week of the study. Blood glucose was measured by withdrawing blood from tail vein. All rats received 2.5g/kg glucose by oral gavage, and blood glucose level was measured at four time points (0, 30, 60, 120 min after gavage). Respectively, the area under the curve of oral glucose tolerance test (AUCOGTT) was measured through 1/4*(a+b)+1/4*(b+c)+1/2*(c+d), a,b,c,d stand for 0, 30, 60, and 120 min’s blood glucose level. The blood glucose level observed between 30, 60, and 120 min, among various groups compared to 0 min, was used for assessing glucose tolerance. The glucometer used was the accu-chek® Performa (Roche Co., Ltd. Shanghai. China).
RNA preparation and microarray assay
The liver samples acquired from the control, DM, and BBR groups (n = 3 per group) were utilized in a microarray assay. Agilent rat gene expression 4x44K microarrays, which cover approximately 26,930 genes, were chosen for this study. The total RNA was purified using a QIAGEN RNeasy Mini Kit according to the manufacturer’s instructions (Qiagen, 217004). The quantity of the total RNA in each sample was determined using a NanoDrop instrument, and the RNA integrity was assessed using an Agilent Bioanalyzer 2100 instrument (Agilent Technology, Santa Clara, CA, USA). The microarray assay was independently repeated in triplicate. After hybridization and washing, the microarray chips were scanned using a DNA chip scanner (Agilent Technologies). Feature Extraction software was used to determine the signal intensities.
GO term enrichment and pathway analysis
The Gene Ontology (GO) platform (http://www.geneontology.org) is a community-based bioinformatics resource in which terms describing gene-product attributes are organized in three structured networks [9]. GO enrichment analysis was performed for the 91 genes with significantly increased expression and the 63 genes with significantly decreased expression using the Gene Ontology database. We defined the significance of GO term enrichment for the DEGs according to the P-value, with a cut-off value of 0.05 [9]. We assessed the liver DEGs according to their association with the biological process, molecular function and cellular component categories.
The pathways corresponding to the DEGs were investigated using the KEGG Orthology Based Annotation System 2.0 database (http://kobas.cbi.pku.edu.cn/), which is the most recent version of KOBAS. The function of KOBAS is to identify the human diseases and significantly enriched pathways that correspond to a given set of genes. The significance and selection of the associated pathways was based on the P-value, with a cut-off value of 0.05 [9].
Protein-protein interaction (PPI) network construction
The Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) database (http://string-db.org/) identifies the interactions of gene products, including not only the direct physical interactions of proteins but also their functional interactions. To evaluate the interactions among the DEGs, we uploaded these genes and drew a color-coded protein-protein interaction network graph. We then imported the PPI data in text format into the Cytoscape program (http://www.cytoscape.org/) to visualize the relationships and used its network analyzer plug-in to analyze the PPI network.
Quantitative real-time PCR analysis
Quantitative real-time PCR (qRT-PCR) analyses were performed using SYBR Green to validate the microarray and PPI network data. Total RNA was isolated from the whole livers of rats (n = 6 per group) using the TaKaRa MiniBEST Universal RNA Extraction kit (TaKaRa, Clontech). The amount and quality of the total RNA were determined using spectrophotometry at 260/280 nm (Thermo Scientific, USA). The RNA was reverse-transcribed using PrimerScriptTM RT Master Mix (Perfect Real Time). Previously published primers were used to amplify the RHOA, MAPK4, DOT1L, SETD2, ME3 and DLAT genes [10,11,12,13,14], and primers obtained from PrimerBank (https://pga.mgh.harvard.edu/primerbank/) were used to amplify the SGK494 gene (Table 1). GAPDH was used as the internal control. The RHOA, MAPK4, SGK494, DOT1L, SETD2, ME3 and DLAT amplification conditions were as follows: pre-denaturation at 95°C for 30 sec, followed by 40 cycles of denaturation at 95°C for 5 sec and annealing at 60°C for 30 sec. The dissolution curve conditions were 65°C for 0.05 sec and 95°C for 0.5 sec. The relative mRNA expression levels in the Con, DM and BBR groups were calculated using the comparative Ct method.
Results
Berberine decreased the random plasma glucose level of DM rats
The random plasma glucose level of the DM rats was significantly higher than that of the control rats at 0 weeks (P<0.01) and at the 1st week (P<0.01), 3rd week (P<0.01), 5th week (P<0.01), 8th week (P<0.01), 10th week (P<0.01), and 12th week of the experimental period (P<0.01). The random plasma glucose level in the BBR group was significantly lower at the 1st week (P<0.01), 3rd week (P<0.01), 5th week (P<0.01), 10th week (P<0.05), and 12th week (P<0.05) of treatment than the levels in the DM group (Fig 1).
▲P<0.05,▲▲P<0.01 the DM group versus the control group; ★P<0.05, ★★P<0.01 the BBR group versus the DM group. The data presented are the mean values ± SD, n = 6.
Berberine reduced the glycosylated hemoglobin (HbA1C) level of DM rats
At 0 weeks and at the 12th week of the experimental period, the HbA1C levels were higher in the DM group than in the control group (P<0.01) (Fig 2). The HbA1C level of the DM rats treated with BBR for 12 weeks was significantly reduced compared with that of the untreated DM rats (P<0.05) [6].
▲▲P<0.01 versus the control group; ★P<0.05, ★★P<0.01 versus the DM group. The data presented are the mean values ± SD, number of rats per group n = 6.
Effect on AST, ALT, BUN and CREA
The biochemical results showed that the DM group with a significant increase (P<0.01 or P<0.05) in the serum marker enzymes aspartate, alanine transaminases and blood urea nitrogen (AST, ALT and BUN) at 0 weeks. Under our experimental conditions, there were no statistical differences in AST, ALT and CREA between the BBR group and DM group (P>0.05), but expect the BUN at the 12th week (P<0.01; Fig 3).
▲P<0.05, ▲▲P<0.01 the DM group versus the control group; ★★P<0.01 the BBR group versus the DM group. Data are expressed as mean ± SD, number of rats per group n = 6. AST: Aspartate aminotransferase; ALT: alanine aminotransferase; BUN: blood urea nitrogen; CREA: creatinine.
Effect on glucose tolerance
As shown in Table 2, the blood glucose level and the area under the curve (AUC) were significantly elevated in the DM group when compared with the normal control group (P<0.01) at 0 weeks and at the 12th week. After 12 weeks of berberine treatment, the area under the curve (AUC) level was significantly reduced in BBR group (P<0.05) compared with DM group.
Microarray analysis of gene expression
The gene expression data obtained using the Agilent rat gene expression 4x44K microarray were extracted using GenePix pro 6.1 software (Agilent Technologies). Quantile normalization and subsequent data processing were performed using the Agilent_Analyze_V1.0 program (Agilent Technologies). Quantile normalization and PCA analysis of the raw data (Fig 4A and Fig 4B) showed that our data were symmetrical and representative. Among the DEGs (S1 Table) detected using microarray analysis, 661 genes had differentially up-regulated expression and 546 genes had differentially down-regulated expression (S2 Table). The genes with log2 (Con vs DM) ≥ 1 and log2 (BBR vs DM) ≥ 1 expression levels and the genes with log2 (Con vs DM) ≤ -1 and log2 (BBR vs DM) ≤ -1 expression levels at P<0.05 were combined, leading to the identification of 154 genes with significantly different expression levels in the livers of the BBR group. Of these genes, 91 had significantly increased expression levels and 63 had significantly decreased expression levels (S3 Table).
A. Box plot of the genes that were differentially expressed in the Con, DM and BBR groups; B. PCA plot of the genes that were differentially expressed in the Con, DM and BBR groups. The plot was constructed based on the fold-change in expression level and the significance of the differences, as determined from the P-values.
GO enrichment and pathway analysis
Both up-regulated and down-regulated genes were observed in the liver samples of the ZDF rats. GO analysis revealed that the functions of these genes are related to many biological processes that are important for the occurrence and development of diabetes, such as the processes in the cellular component organization or biogenesis, metabolic process, biological regulation and response to stimulus categories. The molecular functions associated with these genes include transporter activity, catalytic activity, enzyme regulator activity and antioxidant activity, and the cellular components corresponding to these genes include organelle, cell part, macromolecular complex and extracellular region (Fig 5).
(A) Biological process enrichment; (B) molecular function enrichment; (C) cell component enrichment.
The main pathways for both the up-regulated genes and the down-regulated genes are the lysine degradation, glycosaminoglycan biosynthesis-chondroitin sulfate/dermatan sulfate, thyroid hormone synthesis, hematopoietic cell lineage, and pyruvate metabolism pathways (Table 3).
Interaction network construction and network analysis
We evaluated 154 DEGs using the STRING version 10.0 database to identify the interactions among their gene products; this analysis identified 145 proteins. Then, a PPI network was constructed using Cytoscape software (Fig 6) to provide insight into the molecular mechanisms underlying the actions of the DEGs in the liver However, the network contained many isolated nodes, which did not provide useful information [15]. Therefore, we deleted the DEGs corresponding to isolated nodes and pairs of linked nodes, including Ppm1b, Atp2c1, Sod3, Lyc2, Hrh4, Drd2, Rpl3, RGD1565713, Rps2, Gar1, Hnrnpa3, Clk1, and Srsf7, among others (S3 Table). The resulting network was composed of 20 nodes and 21 edges, with 5 as the maximum degree of connectivity of a node and 1 as the minimum. The average degree of connectivity of the nodes in the network was 2.1. When the degree of connectivity of a node was much greater than average, the node was considered significantly more important than the other nodes and was called the central node. The central node of a PPI network is the core protein in the network. In this study, the nodes were analyzed using a network analyzer, with the degree of connectivity of a node set as ≥ the mean value for all nodes to identify candidate therapeutic target genes. RHOA, MAPK4, Sgk494, Setd2 and DLAT were identified as such genes. Based on the results of the KEGG pathway and protein-protein network analyses, Dot1l and Me3 were also identified as candidate target genes.
Proteins encoded by up-regulated (red) and down-regulated (green) genes are indicated by color.
The seven genes with significantly different expression levels in nine liver samples (three each from the Con, DM, and BBR groups) were used to generate a heat map based on hierarchical clustering analysis of the normalized microarray expression data. The liver samples were collected from randomly selected members of the control, DM, and BBR groups (Fig 7).
Heat map showing the expression levels of the seven genes that were differentially expressed in the liver samples taken from the Con (n = 3), DM (n = 3) and BBR (n = 3) groups. Each row represents one gene, and each column represents a liver sample. The gene symbols are shown on the left side of the rows. The relative gene expression levels are depicted according to the color scale shown at the right. Red indicates high-level expression; blue indicates low-level expression.
Validation of the microarray data using real-time qRT-PCR
The validity of the microarray-based identification of genes with significantly different expression in the liver was confirmed using qRT-PCR. This method identified seven genes with significantly different expression. As shown in Fig 5, the expression levels of DOT1L, SGK494, ME3 and SETD2 were significantly higher in the BBR group than in the DM group (Fig 8A, Fig 8B, Fig 8C and Fig 8D), whereas the expression levels of MAPK4, DLAT and RHOA were significantly lower in the BBR group than in the DM group (Fig 8E, Fig 8F and Fig 8G). The expression level of each of these seven DEGs was significantly different at P<0.05, demonstrating that the qRT-PCR results regarding the seven DEGs were consistent with the microarray results (Fig 8H).
Relative expression levels of the seven DEGs (RHOA, MAPK4, SGK494, DOT1L, SETD2, ME3, and DLAT) determined using qRT-PCR. The P-values were calculated using a one-way ANOVA (**P<0.01;*P<0.05). The qRT-PCR results were consistent with the microarray data.
Discussion
The results of our study showed that berberine treatment of ZDF rats, which exhibit significant hyperglycemia as well as insulin resistance, obesity and hyperlipidemia, significantly reduced their random plasma glucose and HbA1C levels. These results suggest that berberine exerts its anti-diabetic effect in ZDF rats by moderating the diabetic metabolism. These results are similar to those obtained in previous studies [7,16,17,18,19]. The administered dosage (300mg/kg) of berberine was determined according to the equation used for the adult maximal dosage in the clinic [7].The toxicity effect of berberine was evaluated by assaying liver and kidney function markers, the aspartate aminotransferase (AST), alanine aminotransferase (ALT), blood urea nitrogen (BUN), creatinine (CREA). Although the biochemical results showed no marked statistical differences, but we still observed a tendency, the hepatoprotective and renoprotective effect of berberine by regulating liver and kidney function markers under our experimental conditions. The pathogenesis and pathological features of diabetes are not yet clear, and the vast majority of researchers consider diabetes to be a multifactorial disease. GO analysis of the genes with differentially down-regulated and up-regulated expression in the BBR group indicated that these genes were associated with the categories of metabolic process, cellular process (ontology: biological process), catalytic activity, binding (ontology: molecular functions), cell part, and organelle (ontology: cellular functions). Interestingly, in studying the long noncoding RNAs that were differentially expressed in sperm samples collected from diabetic and non-diabetic mice, Jiang [9] found that these RNAs were associated with the GO categories of metabolic process, cell part, binding and catalytic activity. The results of the present study are consistent with these findings. KEGG pathway analysis revealed that some of the genes that were differentially expressed in the BBR group are involved in lysine degradation and pyruvate metabolism. Lysine acetylation is vital for metabolic pathway activity and may regulate the balance between energy storage and expenditure [20]. Pyruvate is a central biochemical node connecting amino acid, carbohydrate, and fatty acid metabolism [21]. We found that the expression levels of DOT1L and SETD2, which are involved in lysine degradation, were up-regulated in the BBR group. In addition, the expression levels of two genes involved in pyruvate metabolism were significantly altered in the BBR group, with down-regulation of DLAT and up-regulation of ME3. Lysine and pyruvate both play important roles in glycometabolism. Consequently, the results of our study provide a basis for research on the pathogenic mechanisms of diabetes at the molecular level.
The results of the PPI network and pathway analyses of the genes that were significantly differentially expressed in the livers of the ZDF rats suggest that seven genes, RHOA, MAPK4, SGK494, DOT1L, SETD2, ME3 and DLAT, play important roles in diabetes and are candidate targets for diabetes treatment.
Ras homolog gene family, member A (RHOA) is an important member of the Rho small G protein superfamily, which regulates cytoskeletal architecture and interacts with the downstream critical target ROCK to affect cell shrinkage, proliferation, migration and adhesion. The RHOA/ROCK signaling pathway plays critical roles in the heart, liver, lungs and kidneys [22]. Excessive activation of the RHOA/ROCK signaling pathway plays an important role in the development of the chronic complications of diabetes [23]. RHOA also involved in renal integrity in diabetic mice, and its inhibitors attenuate diabetic nephropathy in different models of diabetes [24]. RHOA/ROCK signaling activated by TGF-β to regulate Nox4 activity, which might be associated with abnormalities in renal structure [25]. Berberine ameliorates the microvascular complications of diabetes, the renoprotective and anti-oxidative stress effect of BBR by inhibiting the activation of RHOA/ROCK signaling [26]. RHOA was significantly down-regulated in the BBR group in this study, which is consistent with these findings.
Mitogen-activated protein kinase 4 (MAPK4) is a protein serine/threonine kinase that specifically regulates the activity of mitogen-activated protein kinase. The MAPK-mediated signaling pathway plays a major role in cell growth, proliferation, differentiation and apoptosis. In addition, this pathway is involved in insulin resistance and protein phosphorylation. Up-regulated MAPK4 expression can lead to insulin resistance, and MAPK4 overexpression inhibits insulin-stimulated glucose uptake by negatively regulating the insulin signal-transduction pathway [27]. The ZDF rat diabetic model used in this study is a typical model of insulin resistance. The microarray and PCR results showed that the MAPK4 expression was significantly increased in the DM group compared with the control group. In contrast, the MAPK4 expression was significantly decreased in the BBR group compared with the DM group.
Uncharacterized serine/threonine-protein kinase (SGK494) is an important member of the AGC protein kinase family. AGC kinases not only are involved in diverse cellular functions but also may be associated with human diseases, such as diabetes, obesity, inflammation and viral infections [28,29]. To the best of our knowledge, no studies further elucidating the role of this kinase in diabetes have been conducted, and there are no reports of the effect of BBR treatment on SGK494. Our study is the first to show the increased expression of SGK494 in BBR-treated ZDF rats.
DOT1L, a nucleosomal histone 3 (H3)-specific methyltransferase, participates in a variety of cellular processes, including reprogramming, development, differentiation, and proliferation [30]. Zhou et al. found that DOT1L encodes histone H3K79 methyltransferase Dot1a, which contributes to the recovery of renal injury in streptozotocin-induced diabetic rats [31].DOT1L deficiency results in the upregulation of endothelin-1 (ET-1) expression at both mRNA and protein levels in diabetic kidney, and suggest that associated with decreased expression of Dot1a [31]. DOT1L also regulates water homeostasis, and identify AQP5 as a potential novel target of Dot1a, an AQP2 binging partner and regulator. Moreover, depletion of DOT1L led to endogenous AQP5 upregulation, but the upregulated AQP5 contribute to polyuria in patients with diabetic nephropathy, suggesting that DOT1L is critical in renal water homeostasis with DN [32]. In the current study, the PCR results showed higher expression of DOT1L in the BBR group than in the DM group, suggesting that the DOT1L is a potential target about the renoprotective effect of BBR.
SET domain-containing 2 (SETD2) is a histone methyltransferase that trimethylates lysine 36 of H3 (H3K36) and is involved in transcriptional elongation, RNA processing, and DNA repair [33]. The H3K36 site can be modified by methylation, and its modified state is controlled by the dynamic regulation of the methyl transferase and the demethylase that are specific to H3K36. Histone demethylation plays a role in the transcriptional regulation of gluconeogenesis, which has important implications for the treatment of diabetes [34].
Malic enzyme (ME) is a lipogenic enzyme whose activity level is correlated with the rate of fatty acid synthesis [13]. The mitochondrial NADP malic enzyme (ME3) is one isoform of ME. The function of mitochondrial MEs is to supply adequate amounts of pyruvate for the operation of the Krebs cycle [35], and their expression is synergistically regulated by insulin and glucose [36]. Hasan et al. found that ME3 participates in the insulin secretion by pancreatic β-cells, and knockdown of ME3 inhibited insulin release stimulated by glucose, pyruvate or 2-aminobicyclo [2,2,1] heptane-2-carboxylic acid-plus-glutamine, suggesting that ME3 far more than ME1 or ME2 is necessary for insulin release [37]. The knockdown of ME3 led to the inhibition of insulin release might cause by lowering its malic enzyme catalytic activity. In this study, we found that ME3 was expressed at a higher level in the BBR group than in the DM group might suggest that BBR improves insulin release through increase the expression of ME3 and its catalytic activity.
Dihydrolipoamide S-acetyltransferase (DLAT) is the E2 component of the large mitochondrial pyruvate dehydrogenase-containing complex that catalyzes the conversion of pyruvate to acetyl coenzyme A [38], and a responsible gene for pyruvate dehydrogenase complex defect causes of mitochondrial disorders [39]. Moreover, Goguet-Rubio et al revel a central role for DLAT in the pyruvate metabolism program regulated through E4 transcription factor 1 [40]. Lee et al. used cDNA microarray analysis to study the placental gene expression related to glucose metabolism and found nine genes with altered expression levels, including DLAT, which related to the tricarboxylic acid cycle pathway. Through increased mRNA and protein expression of DLAT to increase glycolysis and gluconeogenesis [41]. In this study, we found that the DLAT expression was significantly down-regulated in the BBR-treated ZDF rats. Exact mechanisms whereby BBR and DLAT interact and how this impacts on pyruvate metabolism warrant further investigation.
In this study, we did not focus on the DEGs represented as pairs of linked nodes in the PPI network. However, these genes may play roles in diabetes; one of these genes encodes Ppm1b (protein phosphatase, Mg2+/Mn2+ dependent, 1B, Ppm1b), which increases the expression of PPARγ. In adipocytes, PPARγ not only promotes the removal of lipids and free fatty acids (FFAs) but also inhibits the production of proinflammatory cytokines, ameliorating insulin resistance [42]; this finding suggests that Ppm1b and other genes may be worth further study under certain conditions. Therefore, we will investigate the simply linked genes in a follow-up study.
In conclusion, based on current experimental evidence, our study showed that BBR treatment improved the glucose metabolism and has protective effects against hepatorenal impairment of diabetic rats. We found that the RHOA, MAPK4, SGK494, DOT1L, SETD2, ME3 and DLAT genes were significantly differentially expressed in the livers of the BBR and DM groups and thus are potential targets for the treatment of T2DM. These results provide molecular information that will facilitate elucidating the mechanisms through which berberine alleviates T2DM, and further studies are required to confirm the results of this study.
Supporting Information
S1 Table. The differentially expressed genes identified in the control, DM and BBR rat liver samples.
https://doi.org/10.1371/journal.pone.0166378.s001
(XLSX)
S2 Table. The 661 up-regulated and 546 down-regulated genes identified in the control, DM and BBR rat liver samples.
https://doi.org/10.1371/journal.pone.0166378.s002
(XLSX)
S3 Table. The 91 up-regulated and 63 down-regulated genes identified in the control, DM and BBR rat liver samples.
https://doi.org/10.1371/journal.pone.0166378.s003
(XLSX)
Acknowledgments
This work was supported by the Project of National Great New Drug Research and Development (No. 2012ZX09503001-001), and Traditional Chinese Medicine Open funds of Zhejiang Chinese Medical University (No. 752223A00201/005/019). We are grateful to Ms Yuanxiao Yang, Ms Ying He for constructive discussion, Ms Ling Xuan, Ms Jie Ying Zhang and Mr Xue Ming Chen for literature search.
Author Contributions
- Conceptualization: YSW YTC.
- Formal analysis: YSW YTB ZML.
- Funding acquisition: CYL.
- Investigation: YSW YTB ZML.
- Methodology: YSW YTC.
- Project administration: YSW YTC CYL.
- Resources: CYL.
- Software: YSW.
- Supervision: CYL.
- Validation: YSW YTC ZML.
- Visualization: YSW.
- Writing – original draft: YSW XJZ JNH SJD.
- Writing – review & editing: YSW YTC.
References
- 1. Wei S, Zhang M, Yu Y, Lan X, Yao F, Yan X, et al. Berberine attenuates development of the hepatic gluconeogenesis and lipid metabolism disorder in type 2 diabetic mice and in palmitate-incubated HepG2 cells through suppression of the HNF-4α miR122 pathway. PloS one. 2016; 11(3): e0152097. pmid:27011261
- 2. Pang B, Zhao LH, Zhou Q, Zhao TY, Wang H, Gu CJ, et al. Application of berberine on treating type 2 diabetes mellitus. Int J Endocrinol. 2015;2015: 905749. pmid:25861268
- 3. Zhang Q, Xiao X, Li M, Li W, Yu M, Zhang H, et al. Berberine moderates glucose metabolism through the GnRH-GLP-1 and MAPK pathways in the intestine. BMC Complement Alternat Med. 2014;14: 188. pmid:24912407
- 4. Xu M, Xiao Y, Yin J, Hou W, Yu X, Shen L, et al. Berberine promotes glucose consumption independently of AMP-activated protein kinase activation. PLOS ONE. 2014;9: e103702. pmid:25072399
- 5. Li F, Zhao YB, Wang DK, Zou X, Fang K, Wang KF. Berberine relieves insulin resistance via the cholinergic anti-inflammatory pathway in HepG2 cells. J Huazhong Univ Sci Technolog Med Sci 2016;36: 64–69. pmid:26838742
- 6. Liu C, Wang Z, Song Y, Wu D, Zheng X, Li P, et al. Effects of berberine on amelioration of hyperglycemia and oxidative stress in high glucose and high fat diet-induced diabetic hamsters in vivo. BioMed Res Int. 2015;2015: 313808. pmid:25705654
- 7. Dong Y, Chen YT, Yang YX, Zhou XJ, Dai SJ, Tong JF, et al. Metabolomics study of type 2 diabetes mellitus and the antidiabetic effect of berberine in Zucker diabetic fatty rats using Uplc-ESI-Hdms. Phytother Res. 2016;30: 823–828. pmid:26888689
- 8. Scarda A, Franzin C, Milan G, Sanna M, Dal Prà C, Pagano C, et al. Increased adipogenic conversion of muscle satellite cells in obese Zucker rats. Int J Obes (Lond). 2010;34: 1319–1327. pmid:20231840
- 9. Jiang GJ, Zhang T, An T, Zhao DD, Yang XY, Zhang DW, et al. Differential expression of long noncoding RNAs between sperm samples from diabetic and non-diabetic mice. PloS One. 2016; 11(4):
- 10. Chen J, Li Q, Dong R, Gao H, Peng H, Wu Y. The effect of the Ras homolog gene family (Rho), member A/Rho associated coiled‑coil forming protein kinase pathway in atrial fibrosis of type 2 diabetes in rats. Exp Ther Med. 2014;8: 836–840. pmid:25120610
- 11. Rousseau J, Klinger S, Rachalski A, Turgeon B, Déléris P, Vigneault E, et al. Targeted inactivation of Mapk4 in mice reveals specific nonredundant functions of Erk3/Erk4 subfamily mitogen-activated protein kinases. Mol Cell Biol. 2010;30: 5752–5763. pmid:20956558.
- 12. Suzuki MG, Ito H, Aoki F. Effects of RNAi-mediated knockdown of histone methyltransferases on the sex-specific mRNA expression of imp in the silkworm Bombyx mori. Int J Mol Sci. 2014;15: 6772–6796. pmid:24758924.
- 13. Bragoszewski P, Habior A, Walewska-Zielecka B, Ostrowski J. Expression of genes encoding mitochondrial proteins can distinguish nonalcoholic steatosis from steatohepatitis. Acta Biochim Pol. 2007;54: 341–348. pmid:17502923.
- 14. Lee MH, Jeon YJ, Lee SM, Park MH, Jung SC, Kim YJ. Placental gene expression is related to glucose metabolism and fetal cord blood levels of insulin and insulin-like growth factors in intrauterine growth restriction. Early Hum Dev. 2010;86: 45–50. pmid:20106611
- 15. Guo WG, Zhang Y, Ge D, Zhang YX, Lu CL, Wang Q, et al. Bioinformatics analyses combined microarray identify the desregulated microRNAs in lung cancer. Eur Rev Med Pharmacol Sci. 2013;17: 1509–1516 pmid:23771539
- 16. Tang YH, Sun ZL, Fan MS, Li ZX, Huang CG. Anti-diabetic effects of TongGuanWan, a Chinese traditional herbal formula in C57BL/KsJ-db/db mice. Planta Med 2012;78: 18–23. pmid:22002851
- 17. Martínez-Abundis E, Méndez-Del Villar M, Pérez-Rubio KG, Zuñiga LY, Cortez-Navarrete M, Ramírez-Rodriguez A, et al. Novel nutraceutic therapies for the treatment of metabolic syndrome. World J Diabetes. 2016;7: 142–152. pmid:27076875
- 18. Lan J, Zhao Y, Dong F, Yan Z, Zheng W, Fan J, et al. Meta-analysis of the effect and safety of berberine in the treatment of type 2 diabetes mellitus, hyperlipemia and hypertension. J Ethnopharmacol. 2015;161: 69–81. pmid:25498346
- 19. Di Pierro F, Putignano P, Villanova N, Montesi L, Moscatiello S, Marchesini G. Preliminary study about the possible glycemic clinical advantage in using a fixed combination of berberis aristata and Silybum marianum standardized extracts versus only berberis aristata in patients with type 2 diabetes. Clin Pharmacol. 2013;5: 167–174. pmid:24277991
- 20. Iyer A, Fairlie DP, Brown L. Lysine acetylation in obesity, diabetes and metabolic disease. Immunol Cell Biol. 2012;90: 39–46. pmid:22083525
- 21. Vacanti NM, Divakaruni AS, Green CR, Parker SJ, Henry RR, Ciaraldi TP, et al. Regulation of substrate utilization by the mitochondrial pyruvate carrier. Mol Cell. 2014;56: 425–435. pmid:25458843
- 22. Knipe RS, Tager AM, Liao JK. The Rho kinases: critical mediators of multiple profibrotic processes and rational targets for new therapies for pulmonary fibrosis. Pharmacol Rev. 2015;67: 103–117. pmid:25395505
- 23. Zhou H, Li Y. Long-term diabetic complications may be ameliorated by targeting Rho kinase. Diabetes Metab Res Rev. 2011;27: 318–330. pmid:21309060
- 24. Bruder-Nascimento T, Callera G, Montezano A C, Antunes TT, He Y, Cat AN et al. Renoprotective Effects of Atorvastatin in Diabetic Mice: Downregulation of RhoA and Upregulation of Akt/GSK3. PLoS One. 2016; 11: e0162731. pmid:27649495
- 25. Lan A, Du J. Potential role of Akt signaling in chronic kidney disease. Nephrol Dial Transp. 2015; 30:385–394. pmid:24891436
- 26. Xie X, Chang X, Chen L, Huang K, Huang J, Wang S, et al. Berberine ameliorates experimental diabetes-induced renal inflammation and fibronectin by inhibiting the activation of RhoA/ROCK signaling. Mol Cell Endocrinol. 2013;381: 56–65. pmid:23896433
- 27. Xu H, Dembski M, Yang Q, Yang D, Moriarty A, Tayber O, et al. Dual specificity mitogen-activated protein (MAP) kinase phosphatase-4 plays a potential role in insulin resistance. J Biol Chem. 2003;278: 30187–30192. pmid:12777378
- 28. Arencibia JM, Pastor-Flores D, Bauer AF, Schulze JO, Biondi RM. AGC protein kinases: from structural mechanism of regulation to allosteric drug development for the treatment of human diseases. Biochim Biophys Acta. 2013;1834: 1302–1321. pmid:23524293
- 29. Pearce LR, Komander D, Alessi DR. The nuts and bolts of AGC protein kinases. Nat Rev Mol Cell Biol. 2010;11: 9–22. pmid:20027184
- 30. Kim W, Choi M, Kim JE. The histone methyltransferase Dot1/DOT1L as a critical regulator of the cell cycle. Cell Cycle. 2014;13: 726–738. pmid:24526115
- 31. Zhou Q, Liu K, Wu H, Chen L, Pouranan V, Yuan M, et al. Spironolactone rescues Dot1a-Af9-mediated repression of endothelin-1 and improves kidney injury in streptozotocin-induced diabetic rats. PLoS One. 2012;7: e47360. pmid:23077601
- 32. Wu H, Chen L, Zhang X, Wu H, Chen L, Zhang X, Zhou Q, Li JM, Berger S et al. Aqp5 is a new transcriptional target of Dot1a and a regulator of Aqp2. PLoS One. 2013; 8: e53342. pmid:23326416
- 33. Al Sarakbi W, Sasi W, Jiang WG, Roberts T, Newbold RF, Mokbel K. The mRNA expression of SETD2 in human breast cancer: correlation with clinico-pathological parameters. BMC Cancer. 2009;9: 290. pmid:19698110
- 34. Pan D, Mao C, Zou T, Yao AY, Cooper MP, Boyartchuk V, et al. The histone demethylase Jhdm1a regulates hepatic gluconeogenesis. PLoS Genet. 2012;8: e1002761. pmid:22719268
- 35. Pongratz RL, Kibbey RG, Shulman GI, Cline GW. Cytosolic and mitochondrial malic enzyme isoforms differentially control insulin secretion. J Biol Chem. 2007;282: 200–207. pmid:17102138
- 36. Dentin R, Girard J, Postic C. Carbohydrate responsive element binding protein (ChREBP) and sterol regulatory element binding protein-1c (SREBP-1c): two key regulators of glucose metabolism and lipid synthesis in liver. Biochimie. 2005;87: 81–86. pmid:15733741
- 37. Hasan NM, Longacre MJ, Stoker SW, Kendrick MA, MacDonald MJ. Mitochondrial Malic Enzyme 3 is important for insulin secretion in pancreatic β-cells. Mol Endocrinol. 2015;29: 396–410. pmid:25594249.
- 38. Head RA, Brown RM, Zolkipli Z, Shahdadpuri R, King MD, Clayton PT, et al. Clinical and genetic spectrum of pyruvate dehydrogenase deficiency: dihydrolipoamide acetyltransferase (E2) deficiency. Ann Neurol. 2005;58: 234–241. pmid:16049940
- 39. Ciara E,Rokicki D,Halat P, Karkucińska-Więckowska A,Piekutowska-Abramczuk D, Mayr J et al. Difficulties in recognition of pyruvate dehydrogenase complex deficiency on the basis of clinical and biochemical features. The role of next-generation sequencing. Mol Genet Metab Rep. 2016; 7: 70–76. pmid:27144126
- 40. Goguet-Rubio P, Seyran B, Gayte L, Bernex F, Sutter A, Delpech H, et al. E4F1-mediated control of pyruvate dehydrogenase activity is essential for skin homeostasis. Proc Natl Acad Sci U S A.2016; 11: 11004–11009. doi: https://doi.org/10.1073/pnas.1602751113
- 41. Lee MH, Jeon YJ, Lee SM, Park MH, Jung SC, Kim YJ. Placental gene expression is related to glucose metabolism and fetal cord blood levels of insulin and insulin-like growth factors in intrauterine growth restriction. Early Hum Dev. 2010;86: 45–50. pmid:20106611
- 42. Tasdelen I, van Beekum O, Gorbenko O, Fleskens V, van den Broek NJ, Koppen A, et al. The serine/threonine phosphatase PPM1B (PP2Cβ) selectively modulates PPARγ activity. Biochem J. 2013; 451: 45–53. pmid:23320500