Skip to main content
Download PDF

Identifying key genes for diabetic kidney disease by bioinformatics analysis

BMC Nephrology volume 24, Article number: 305 (2023) Cite this article

Abstract

Background

There are no reliable molecular targets for early diagnosis and effective treatment in the clinical management of diabetic kidney disease (DKD). To identify novel gene factors underlying the progression of DKD.

Methods

The public transcriptomic datasets of the alloxan-induced DKD model and the streptozotocin-induced DKD model were retrieved to perform an integrative bioinformatic analysis of differentially expressed genes (DEGs) shared by two experimental animal models. The dominant biological processes and pathways associated with DEGs were identified through enrichment analysis. The expression changes of the key DEGs were validated in the classic db/db DKD mouse model.

Results

The downregulated and upregulated genes in DKD models were uncovered from GSE139317 and GSE131221 microarray datasets. Enrichment analysis revealed that metabolic process, extracellular exosomes, and hydrolase activity are shared biological processes and molecular activity is altered in the DEGs. Importantly, Hmgcs2, angptl4, and Slco1a1 displayed a consistent expression pattern across the two DKD models. In the classic db/db DKD mice, Hmgcs2 and angptl4 were also found to be upregulated while Slco1a1 was downregulated in comparison to the control animals.

Conclusions

In summary, we identified the common biological processes and molecular activity being altered in two DKD experimental models, as well as the novel gene factors (Hmgcs2, Angptl4, and Slco1a1) which may be implicated in DKD. Future works are warranted to decipher the biological role of these genes in the pathogenesis of DKD.

Peer Review reports

Introduction

Diabetic kidney disease (DKD) is one of the most prevalent complications of diabetes and the main cause of end-stage renal disease (ESRD), accounting for ~ 50% of the total cases of ESRD in the United States and affecting approximately 30% of the global population [1,2,3]. Currently, there are no reliable molecular targets for early diagnosis and effective treatment in the clinical management of DKD [4, 5]. Therefore, identifying novel pathogenic factors closely associated with DKD development and progression not only facilitates the understanding of the molecular mechanisms of DKD, but also provides novel targets for therapeutic opportunities.

The etiology and pathogenesis of DKD is complex, which involves hierarchical physiological and pathological processes [3, 6, 7]. Different risk factors have been recognized as contributors to DKD, including lipid metabolic process [8], inflammation [9], oxidation–reduction process and ischemia [10] and oxidative stress [11]. Recent molecular and cellular researches explored new fields of pathophysiology of DKD, such as mitochondria dysfunction [12], podocyte autophagy [13], and genetic and epigenetic regulation [14]. Various types of animal models of DKD have been developed as valuable tools to investigate the pathogenesis of DKD and uncover important molecular events in the progression of DKD [15,16,17]. Nevertheless, different chemical inducers may incur a distinct physiological and pathological process, which may not be the direct driving factors or consequences of DKD. The identification of common gene factors helps define the biological changes shared by different models, and these genes are more likely to be disease-relevant and serve as a targetable factor for therapeutic interventions.

The rapid development of bioinformatic analysis has provided a platform to systematically survey the transcriptomic changes associated with disease conditions [18]. In this study, we retrieved and reanalyzed microarray datasets from the Gene Expression Omnibus (GEO), which are derived from an alloxan- or streptozotocin-induced DKD mouse models. We reported three previously uncharacterized genes, Hmgcs2, Angptl4, and Slco1a1, which showed consistent changes in the two animal models, and their expression changes were also verified in the classical db/db DKD mouse. Future functional characterization of these genes in DKD animal models will shed light on the unappreciated mechanisms of DKD progression.

Materials and methods

Gene expression omnibus dataset

Microarray data were downloaded from the GEO database (http://www.ncbi.nlm.nih.gov/geo). The GSE139317 dataset generated on the GPL21163 Agilent-074809 Sure-Print G3 Mouse GE v2 8 × 60 K Microarray platform comprises of 6 control and 9 DKD samples from alloxan-treated mice. The GSE131221 dataset derived from the GPL22740 Agilent-074036 Sure-Print G3 Rat GE v2 8 × 60 K Microarray G4858A platform contains 5 control and 7 DKD samples from streptozotocin-induced rats.

Identification of differentially expressed genes (DEGs)

The Bioconductor R package was utilized for the differential gene expression analysis of the microarray data. Genes with a p value less than 0.05 and |log2-fold change (FC)| greater than 1.5 were considered the DEGs with statistical significance.

Functional enrichment analysis

All identified DEGs in GSE139317 and GSE131221 were subjected to functional enrichment analysis by using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway database (https://www.kegg.jp/kegg/kegg1.html, KEGG is developed by Kanehisa Laboratories) [19]. For data annotation, DAVID bioinformatics resources (https://david.ncifcrf.gov/) was used for gene id conversion. Gene ontology enrichment analysis was performed using G:Profiler: a web server for functional enrichment analysis, with the threshold for significance at an adjusted p value < 0.05. In order to build a protein–protein interaction (PPI) network, we used the String database (https://string-db.org/) to retrieve the relationship among the candidates. The PPI network was constructed by Cytoscape software using a minimum Required Interaction score at 0.9 and with the parameter "hide unconnected nodes in the network.”

Animal model of DKD

Male C57BL/6 mice (15 weeks; n = 5) and male db/db mice (15 weeks; n = 5) were purchased from Cavens Laboratory Animal Co., Ltd (Changzhou, China). The animals were fed with a standard rodent diet and water ad libitum. Diabetes mellitus in mice was defined as a blood glucose level > 16.7 mmol/L for 3 consecutive days, and the animals were euthanized after 16 weeks using a euthanizing chamber connected to a CO2 cylinder. HE staining was performed to verify the disrupted glomerular capillary structures, and PAS staining showed that the glomerular basement membrane was thickened and obvious glycogen deposition, which is the sign of DKD in the db/db mouse model [19].

Quantitative real-time PCR

Total RNA extracted from the renal tissues using RNAiso Plus reagent (9109, TaKaRa, Kyoto, Japan). 1 μg of RNA sample was reverse-transcribed into complementary DNA using the first strand cDNA synthesis kit (RR047A; TaKaRa, Kyoto, Japan). Quantitative real-time PCR was conducted using SYBR green-based assay (11201ES03; Yeasen Biotech, Shanghai, China) on the CFX96 Real-Time system (Bio-Rad, Hercules, CA). β-actin gene was used as the internal reference for relative gene expression analysis based on the 2−ΔΔCt method. The sequences of primers for real-time PCR are in the Table 1.

Table 1 Primers for qPCR in db/db mouse
Full size table

Statistics

Statistical analysis was performed using the GraphPad Prism software (USA). Data were presented as means ± standard deviation (SD). The significant differences were analyzed using the unpaired Student’s -t-test between the two groups. Data with p values of 0.05 or less were considered statistically significant.

Results

Project workflow

In the current study, we formulated an integrative analysis plan (Fig. 1) to identify the potential pathogenic factors involved in DKD. A comparative analysis of microarray datasets derived from two different DKD models was performed to profile the DEGs, the associated biological processes, and KEGG pathway. The common DEGs of the two datasets were validated in the classical DKD model, which may be employed as candidates for diagnosis or treatment.

Fig. 1
figure 1

Study workflow. DKD, Diabetic kidney disease; STZ, streptozotocin; DEGs, differentially expressed genes; KEGG, Kyoto Encyclopedia of Genes and Genomes

Full size image

Differential expression analysis of genes in alloxan-induced DKD model

Alloxan-induced diabetic kidney malfunction was widely used to elucidate the pathophysiology of DKD. Gene expression profiles of 9 DKD and 6 sham samples were quality-checked and all data showed a median centered distribution (Fig. 2A), indicating they are statistically comparable. The volcano map in Fig. 2B, C displays the general distribution of gene expression changes, indicating an equivalent number of genes being upregulated or downregulated in DKD samples. In order to better screen the key genes with dramatic changes in DKD, we tightened the criteria of DEGs by using a |log2 FC| greater than 1.5 and an adjusted P value less than 0.05. As a result, 13 upregulated DEGs and 27 downregulated DEGs were identified between the DKD group and the sham group (shown in the Table 2), and the top 10 DEGs (Hmgcs2, Angptl7, Anxa13, Gm10639, Gsta2, Cyp4a10, Il34, Lpl, Inmt, Bhmt) were shown in Fig. 2D.

Fig. 2
figure 2

Differential expression analysis of GSE139317 dataset of alloxan-induced DKD model. A Boxplot presentation of gene expression distribution of 9 DKD and 6 sham samples. B, C Volcano map of DEGs. The DEGs were identified between the DKD group and the sham group with the criteria of |log2 FC| greater than 1 and the adjusted P value less than 0.05. D Top 10 DEGs identified between the DKD group and the sham group with the criteria of |log2 FC| greater than 1.5 and the adjusted P value less than 0.05

Full size image
Table 2 Upregulated and downregulated genes in alloxan-induced DKD
Full size table

To gain more biological insights into the DEGs, functional enrichment analysis of the DEGs was performed to explore the Gene Ontology (GO) annotations and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. The significantly enriched biological process mostly included metabolic process, oxidation–reduction process, lipid metabolic process, extracellular exosomes, extracellular region, hydrolase activity and heme binding (Fig. 3A-C). Metabolic pathways, retinol metabolism, PPAR signaling pathway and fatty acid degradation were the top KEGG pathways associated with DEGs (Fig. 3D). Furthermore, the protein–protein interaction (PPI) network analysis of the DGEs revealed a key module comprising of multiple members of the cytochrome P450 superfamily of enzymes (Cyp2b10, Cyp4a10, Cyp4a12a, Cyp4a12b), which are involved in the oxidation–reduction metabolism of fatty acids (Fig. 3E). Combining the above data, we conclude that the oxidation–reduction, metabolism, extracellular exosome, extracellular region, and hydrolase activity are the major biological process and cellular components altered in the alloxan-induced DKD model. The lipid-related metabolism suggests that DKD can be improved through diet and exercise, while the extracellular matrix suggests that fibrosis is also a point of concern in DKD.

Fig. 3
figure 3

Functional enrichment analysis of the DEGs. A-C GO terms significantly enriched in the DGEs: BP, biological processes; CC, cellular component; MF, molecular function. D KEGG pathway enrichment analysis of DGEs. E Protein‐protein interaction network of DGEs. The above analyses were based on DEGs with the criteria of |log2 FC| greater than 1.5 and the adjusted P value less than 0.05

Full size image

Profiling of DEGs in streptozotocin (STZ)-induced DKD model

Next, we profiled the DEGs in streptozotocin (STZ)-induced DKD model. The dataset comprises of 7 DKD and 5 sham samples (Fig. 4A), in which the DEGs were initially identified with the criteria of |log2 FC| greater than 1 and the adjusted P value less than 0.05 (Fig. 4B, C). Likewise, the criteria of |log2 FC| greater than 1.5 and the adjusted P value less than 0.05 were further used to identify the DEGs with dramatic changes. 27 upregulated genes and 34 downregulated genes between the DKD and the sham groups were identified (shown in the Table 3), and the top 10 DEGs (Car3, Tff3, Abcb1b, Slco1a1, Kap, Prima1, Rgn, LOC100911353, Slc7a12, Angptl4) were shown in Fig. 4D. In this regard, the GO and KEGG enrichment analysis showed a broad spectrum of biological processes and pathways were associated with the DEGs, including the response to drugs, inflammatory response, metabolic process, negative regulation of apoptotic process, ion transport, extracellular exosomes, extracellular space, hydrolase activity, complement and coagulation cascade, bile secretion, glutathione metabolism and histidine metabolism (Fig. 5A-D). Consistent with these data, PPI network indicated a wide range of DEGs involved in different cellular processes were implicated in STZ-induced DKD (Fig. 5E). Together, these data indicate that except for the extracellular exosome, extracellular space, and hydrolase activity, there are other biological processes such as complement and coagulation cascade, bile secretion, glutathione metabolism and histidine metabolism that are affected in the pathophysiological process of STZ-induced DKD model.

Fig. 4
figure 4

Differential expression analysis of GSE131221 dataset of streptozotocin (STZ)-induced DKD model. A Boxplot presentation of gene expression profiles of 7 DKD and 5 sham samples. B, C Volcano map of DEGs. The DEGs were identified between the DKD group and the sham group with the criteria of |log2 FC| greater than 1 and adjusted P value less than 0.05. D Top 10 DEG lists identified between the DKD and the sham group with the criteria of |log2 FC| greater than 1.5 and the adjusted P value less than 0.05

Full size image
Table 3 Upregulated and downregulated DEGs in STZ-induced DKD
Full size table
Fig. 5
figure 5

Enrichment analysis. A-C GO terms significantly enriched in the DGEs: BP, biological processes; CC, cellular component; MF, molecular function. D KEGG pathway enrichment analysis related to DGEs. E Protein‐protein interaction networks related to DGEs

Full size image

Identifying the common DEGs in two DKD models

The above analyses demonstrated that although there are common biological processes affected in both models of DKD, there are specific processes exclusively enriched in each model, some changes incurred may be unique to the chemical inducer in alloxan and STZ-induced DKD. These changes are unlikely to be the true signals signifying the pathogenic process of DKD. To find the gene factors shared by the two models, we analyzed the intersection of the downregulated and upregulated genes with the criteria of |log2 FC| greater than 1.5 and the adjusted P value less than 0.05. (Fig. 6A, B). Three genes, Hmgcs2, Angptl4, and Slco1a1 displayed a consistent change in both alloxan and STZ-induced DKD models (Fig. 6C). As shown in Fig. 6D, E, Hmgcs2 and Angptl4 showed dramatic upregulation while Slco1a1 expression was suppressed in both DKD models. Relevant to our analyses, the activity of Hmgcs2 in the kidneys and 24-h urinary excretion of the ketone body β-hydroxybutyrate (β-HB) were reported to be increased in db/db mice [20]. The upregulation of Angptl4 was found in the development of diabetic nephropathy [21], indicating a potential role for Angptl4 in DN for the detection of a diabetic kidney disease [22] and as a therapeutic target [23]. Consistently, a previous study reported that TNF-α and TGF-β1 downregulated Slco1a1 in primary hepatocytes, suggesting a proinflammatory role [24] and hypermethylation of Slco1a1 locus in diabetic mice was observed [25]. Taken together, our analyses suggest the potential contribution of Hmgcs2, Angptl4, and Slco1a1 to the pathogenesis of DKD, and the functional characterization of these genes is conducive to the understanding of the underlying mechanisms of DKD.

Fig. 6
figure 6

The identification of common DEGs affected in the two DKD models. A, B The number of DEGs was significantly upregulated or downregulated in GSE139317 or GSE131221, with the criteria of |log2 FC| greater than 1.5 and the adjusted P value less than 0.05. C The Venn diagram shows the common DEGs of GSE139317 and GSE131221 datasets. D, E The relative expression of Hmgcs2, Angptl4, and Slco1a1 in alloxan-induced model (GSE139317) and STZ-induced DKD model (GSE131221). *P < 0.05, **P < 0.01, ***P < 0.001

Full size image

Validation of the expression change of Hmgcs2, Angptl4 and Slco1a1 in db/ db DKD mice

To further confirm the changes of these genes, we used the db/db mice as the well-established animal model of DKD study [26], to detect the expression pattern of Hmgcs2, Angptl4 and Slco1a1 (Fig. 7A). To confirm the onset of diabetic symptoms, we measured the blood glucose level (Fig. 7B) and body weight (Fig. 7C) of the wildtype controls (WT) and db/db mice at the end of week 16, which showed an elevated blood glucose level and gain in body weight. H&E and PAS staining of renal tissue sections showed the disrupted glomerular capillary structures, verifying the onset of DKD in the db/db mouse model (Fig. 7D). Compared with the WT mice, the expression of Hmgcs2 (Fig. 7E) and Angptl4 (Fig. 7F) were significantly increased in the renal tissues of db/db mice, while the expression of Slco1a1 (Fig. 7G) was suppressed in db/db mice. Taken together, these data further support the uncharacterized roles of Hmgcs2, Angptl4, and Slco1a1 in the pathogenesis of DKD.

Fig. 7
figure 7

The expression of Hmgcs2, Angptl4, and Slco1a1 in db/db mice. A Schematic illustration of the experimental design. B The blood glucose levels in db/db mice and WT mice were measured at week 16 (n = 5 per group). C The body weights in db/db mouse, and b/db mouse, and WT mice were measured at week 16 (n = 5 per group). D-F The relative expressions of Hmgcs2, Angptl4, and Slco1a1 in db/db mice and b/db mice and WT mice were measured at week 16. *P < 0.05, **P < 0.01, ***P < 0.001

Full size image

Discussion

Despite intense efforts to understand the signaling pathways and molecular players driving DKD, the pathogenesis of DKD remains as an intricate process. Because of the fact that a single DKD model may not comprehensively reflect the true signals underlying DKD progression, the integrative analysis of different models could provide an unbiased view of the biological processes and key genes contributing to the progression of DKD. In the current study, we conducted in-depth analysis of the transcriptome data of two widely used DKD models (alloxan-induced model and streptozotocin-induced model). We revealed the common and unique DEGs and biological processes altered in these two DKD models. Most importantly, three genes, Slco1a1, Angptl4, and Hmgcs2 displayed a consistent change across the two models, which was further verified in db/db mice, implying an important role of Slco1a1, Angptl4 and Hmgcs2 as candidate targets for the diagnosis and intervention of DKD.

The pathogenesis of DKD is complex, which is driven by both genetic and environmental factors [27, 28]. Several genes have been proposed to influence the initiation and development of DKD. The C-X-C motif chemokine ligand 1 (CXCL1) was indicated to serve as a proinflammatory mediatorin the DKD [29]. In addition, von willebrand factor (VWF) involved in intrarenal thrombosis was suggested to mediate the deterioration of renal function [30]. Importantly, Spleen tyrosine kinase (SYK) is a protein kinase mediating the secretion of high glucose-induced IL-1β [31] and TGF-β [32]. Metabolic changes associated with obesity, diabetes, aging, and nutritional status also play critical roles in DKD development [33]. Furthermore, studies reported that effective control of smoking, blood pressure, blood glucose, and lipid intake could significantly improve the DKD [34]. However, given the complexity of DKD and the limited strategies of treatment in human patients [35, 36], it is imperative to uncover novel molecular players and develop targeted therapies for DKD management.

Alloxan and streptozotocin (STZ) are the most widely used chemicals to induce DKD in animal models [15, 37]. However, alloxan and streptozotocin can each trigger gene expression changes which are irrelevant to the DKD, thus the biased conclusion may not reflect the true signals of DKD in a single animal model. For example, our analysis demonstrated that the GO annotations and KEGG pathway of DEGs in alloxan-induced DKD are mainly involved in metabolic process, oxidation–reduction process, lipid metabolic process, response to stilbenoid, ethanol catabolic process, triglyceride biosynthetic process, metabolic pathways, retinol metabolism, PPAR signaling pathway and fatty acid degradation. While, in the STZ-induced DKD, the mainly significant enriched biological process includes transport, response to drug, inflammatory response, metabolic process, extracellular exosomes, extracellular space, hydrolase activity, negative regulation of apoptotic process, ion transport, complement and coagulation cascades, bile secretion, glutathione metabolism, histidine metabolism. Therefore, the common processes affected in both models are the metabolic processes, extracellular exosome, extracellular space and hydrolase activity in both models.

Through our PPI network analysis, several members of the network were previously linked to DKD progression, such as Slc22a1 [38], Cyp24a1 [39], Gstm1 [40] and Vwf [41]. Among them, Hmgcs2, Angptl4, and Slco1a1 were the genes of interest that are commonly affected in alloxan and STZ-induced DKD. The activity of Hmgcs2 in the kidneys and 24-h urinary excretion of the ketone body β-hydroxybutyrate (β-HB) were found to be increased in db/db mice [19, 20]. Upregulation of Angptl4 in diabetic nephropathy was proposed to contribute to the development of diabetic nephropathy [21], indicating a potential role of Angptl4 for the detection of a diabetic kidney disease [22] and as a potential therapeutic target [23]. Furthermore, a previous study reported that TNF-α and TGF-β1 downregulate Slco1a1 in primary hepatocytes to modulate hepatic fibrosis and inflammation, suggesting an proinflammatory role of Slco1a1 [24], and the hypermethylation of Slco1a1 locus in diabetic mice was observed [25]. Given the fact that these genes are commonly altered in two different DKD models and their potential functional engagement in DKD progression, future works are warranted to functionally characterize their roles in DKD to provide insights into the molecular mechanisms of DKD and help formulate novel intervention strategies.

In summary, the current study identified Hmgcs2, Angptl4, and Slco1a1 as critical genes in the pathogenesis of two DKD animal models, which are of significance in the understanding of the underlying mechanisms of DKD and may serve as promising targets for DKD clinical management. However, further verification of these genes in clinical samples and functional characterization are needed to confirm their implications in DKD.

Availability of data and materials

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

References

  1. Ogurtsova K, et al. IDF Diabetes Atlas: Global estimates for the prevalence of diabetes for 2015 and 2040. Diabetes Res Clin Pract. 2017;128:40–50.

    Article  CAS  PubMed  Google Scholar 

  2. Saeedi P, et al. Global and regional diabetes prevalence estimates for 2019 and projections for 2030 and 2045: Results from the International Diabetes Federation Diabetes Atlas, 9(th) edition. Diabetes Res Clin Pract. 2019;157: 107843.

    Article  PubMed  Google Scholar 

  3. Lv JC, Zhang LX. Prevalence and disease burden of chronic kidney disease. Adv Exp Med Biol. 2019;1165:3–15.

    Article  CAS  PubMed  Google Scholar 

  4. Wang C, et al. New urinary biomarkers for diabetic kidney disease. Biomark Res. 2013;1(1):9.

    Article  PubMed  PubMed Central  Google Scholar 

  5. National Kidney Foundation. KDOQI Clinical Practice Guideline for Diabetes and CKD: 2012 Update. Am J Kidney Dis. 2012;60(5):850–86.

  6. Romagnani P, et al. Chronic kidney disease. Nat Rev Dis Primers. 2017;3:17088.

    Article  PubMed  Google Scholar 

  7. Wang YN, et al. Chronic kidney disease: Biomarker diagnosis to therapeutic targets. Clin Chim Acta. 2019;499:54–63.

    Article  CAS  PubMed  Google Scholar 

  8. Gai Z, et al. Lipid accumulation and chronic kidney disease. Nutrients. 2019;11(4):722.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Mihai S, Codrici E, Popescu ID. Inflammation-related mechanisms in chronic kidney disease prediction, progression, and outcome. J Immunol Res. 2018;2018:2180373.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Pefanis A, et al. Regulated necrosis in kidney ischemia-reperfusion injury. Kidney Int. 2019;96(2):291–301.

    Article  PubMed  Google Scholar 

  11. Daenen K, et al. Oxidative stress in chronic kidney disease. Pediatr Nephrol. 2019;34(6):975–91.

    Article  PubMed  Google Scholar 

  12. Galvan DL, Green NH, Danesh FR. The hallmarks of mitochondrial dysfunction in chronic kidney disease. Kidney Int. 2017;92(5):1051–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Liu WJ, et al. Lysosome restoration to activate podocyte autophagy: a new therapeutic strategy for diabetic kidney disease. Cell Death Dis. 2019;10(11):806.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Nowak N. Protective factors as biomarkers and targets for prevention and treatment of diabetic nephropathy: From current human evidence to future possibilities. J Diabetes Investig. 2020;11(5):1085–96.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Ighodaro OM, Adeosun AM, Akinloye OA. Alloxan-induced diabetes, a common model for evaluating the glycemic-control potential of therapeutic compounds and plants extracts in experimental studies. Medicina (Kaunas). 2017;53(6):365–74.

    Article  PubMed  Google Scholar 

  16. Pérez-López L, et al. Animal models and renal biomarkers of diabetic nephropathy. Adv Exp Med Biol. 2021;1307:521–51.

    Article  PubMed  Google Scholar 

  17. Brosius FC 3rd, et al. Mouse models of diabetic nephropathy. J Am Soc Nephrol. 2009;20(12):2503–12.

    Article  PubMed  Google Scholar 

  18. Petryszak R, et al. Expression Atlas update–a database of gene and transcript expression from microarray- and sequencing-based functional genomics experiments. Nucleic Acids Res. 2014;42:D926-32.

    Article  CAS  PubMed  Google Scholar 

  19. Kanehisa M, Furumichi M, Sato Y, Kawashima M, Ishiguro-Watanabe M. KEGG for taxonomy-based analysis of pathways and genomes. Nucleic Acids Res. 2023;51(D1):D587–92.

    Article  CAS  PubMed  Google Scholar 

  20. Tomita I, et al. SGLT2 inhibition mediates protection from diabetic kidney disease by promoting ketone body-induced mTORC1 inhibition. Cell Metab. 2020;32(3):404-419.e6.

    Article  CAS  PubMed  Google Scholar 

  21. Ma J, et al. Upregulation of podocyte-secreted angiopoietin-like-4 in diabetic nephropathy. Endocrine. 2015;49(2):373–84.

    Article  CAS  PubMed  Google Scholar 

  22. Al Shawaf E, Abu-Farha M. ANGPTL4: A Predictive Marker for Diabetic Nephropathy. J Diabetes Res. 2019;2019:4943191.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Chugh SS, et al. Angiopoietin-like 4 based therapeutics for proteinuria and kidney disease. Front Pharmacol. 2014;5:23.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Tanaka N, et al. Disruption of phospholipid and bile acid homeostasis in mice with nonalcoholic steatohepatitis. Hepatology. 2012;56(1):118–29.

    Article  CAS  PubMed  Google Scholar 

  25. Marumo T, et al. Diabetes induces aberrant DNA methylation in the proximal tubules of the kidney. J Am Soc Nephrol. 2015;26(10):2388–97.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Sharma K, McCue P, Dunn SR. Diabetic kidney disease in the db/db mouse. Am J Physiol Renal Physiol. 2003;284(6):F1138–44.

    Article  CAS  PubMed  Google Scholar 

  27. Reidy K, et al. Molecular mechanisms of diabetic kidney disease. J Clin Invest. 2014;124(6):2333–40.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Lin YC, et al. Update of pathophysiology and management of diabetic kidney disease. J Formos Med Assoc. 2018;117(8):662–75.

    Article  CAS  PubMed  Google Scholar 

  29. Niu S, Bian Z, Tremblay A. Broad Infiltration of Macrophages Leads to a Proinflammatory State in Streptozotocin-Induced Hyperglycemic Mice. J Immunol. 2016;197(8):3293–301.

    Article  CAS  PubMed  Google Scholar 

  30. Dhanesha N, et al. ADAMTS13 retards progression of diabetic nephropathy by inhibiting intrarenal thrombosis in mice. Arterioscler Thromb Vasc Biol. 2017;37(7):1332–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Qiao Y, et al. Spleen tyrosine kinase promotes NLR family pyrin domain containing 3 inflammasome-mediated IL-1β secretion via c-Jun N-terminal kinase activation and cell apoptosis during diabetic nephropathy. Mol Med Rep. 2018;18(2):1995–2008.

    CAS  PubMed  PubMed Central  Google Scholar 

  32. Yang WS, et al. Spleen tyrosine kinase mediates high glucose-induced transforming growth factor-β1 up-regulation in proximal tubular epithelial cells. Exp Cell Res. 2012;318(15):1867–76.

    Article  CAS  PubMed  Google Scholar 

  33. Winnicki E, et al. Use of the kidney failure risk equation to determine the risk of progression to end-stage renal disease in children with chronic kidney disease. JAMA Pediatr. 2018;172(2):174–80.

    Article  PubMed  Google Scholar 

  34. Thomas MC, et al. Diabetic kidney disease. Nat Rev Dis Primers. 2015;1:15018.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Doshi SM, Friedman AN. Diagnosis and management of type 2 diabetic kidney disease. Clin J Am Soc Nephrol. 2017;12(8):1366–73.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Alicic RZ, Rooney MT, Tuttle KR. Diabetic kidney disease: challenges, progress, and possibilities. Clin J Am Soc Nephrol. 2017;12(12):2032–45.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Lenzen S. The mechanisms of alloxan- and streptozotocin-induced diabetes. Diabetologia. 2008;51(2):216–26.

    Article  CAS  PubMed  Google Scholar 

  38. Sallinen R, et al. Association of the SLC22A1, SLC22A2, and SLC22A3 genes encoding organic cation transporters with diabetic nephropathy and hypertension. Ann Med. 2010;42(4):296–304.

    Article  CAS  PubMed  Google Scholar 

  39. Tourigny A, et al. CYP24A1 exacerbated activity during diabetes contributes to kidney tubular apoptosis via caspase-3 increased expression and activation. PLoS ONE. 2012;7(10): e48652.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Yang Y, et al. Glutathione S-transferase T1 deletion is a risk factor for developing end-stage renal disease in diabetic patients. Int J Mol Med. 2004;14(5):855–9.

    CAS  PubMed  Google Scholar 

  41. Hirano T, et al. Vascular endothelial markers, von Willebrand factor and thrombomodulin index, are specifically elevated in type 2 diabetic patients with nephropathy: comparison of primary renal disease. Clin Chim Acta. 2000;299(1–2):65–75.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

Authors appreciate all participants for their continuous participation and contributions in this research project.

Funding

This work was supported by Grants from the National Natural Science Foundation of China (81960156), High Level Health Technology Personnel Training Project of Yunnan Province(D-2018044), Health Science and Technology Program of Yunnan Province[202201AY070001-061],Yunnan Province Clinical Research Center for Metabolic diseases (202102AA100056) and Yunnan Province Clinical Research Center for kidney diseases (202102AA100060.

Author information

Author notes

  1. Yushan Xu and Lan Li contributed equally to this work.

Authors and Affiliations

  1. Department of Endocrinology, The First Affiliated Hospital of Kunming Medical University, Kunming, 650031, China

    Yushan Xu

  2. Department of Diabetes, The First Affiliated Hospital of Kunming Medical University, Kunming, 650031, China

    Lan Li

  3. Department of General Practice, The First Affiliated Hospital of Kunming Medical University, Kunming, 650031, China

    Ping Tang, Jingrong Zhang, Ruxian Zhong, Jingmei Luo, Jie Lin & Lihua Zhang

Authors
  1. Yushan Xu
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Lan Li
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Ping Tang
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Jingrong Zhang
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Ruxian Zhong
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Jingmei Luo
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Jie Lin
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Lihua Zhang
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

LHZ designed and supervised the study. YSX and LL explored literature. PT, LL and RXZ collected the data. YSX, PT, RXZ, and LL analyzed the data. JML, PT and LHZ wrote the paper. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Lihua Zhang.

Ethics declarations

Ethics approval and consent to participate

Animal protocols were approved by the Ethics Committee of the Kunming Medical University and were conducted in accordance with the guidelines outlined by the Guide for the Care and Use of Laboratory Animals published by the National Institutes of Health. The study was carried out in compliance with the ARRIVE guidelines.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1:

Supporting Table S1. Primers for qPCR in db/db mouse.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Xu, Y., Li, L., Tang, P. et al. Identifying key genes for diabetic kidney disease by bioinformatics analysis. BMC Nephrol 24, 305 (2023). https://doi.org/10.1186/s12882-023-03362-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12882-023-03362-4

Keywords

BMC Nephrology

ISSN: 1471-2369

Contact us