- Research
- Open access
- Published:
Identification and validation of immune-related biomarkers and potential regulators and therapeutic targets for diabetic kidney disease
BMC Medical Genomics volume 16, Article number: 90 (2023)
Abstract
Background
Diabetic kidney disease (DKD) is a major complication of diabetes and the leading cause of end-stage renal disease worldwide. Renal inflammation and infiltration of immune cells contribute to the development and progression of DKD. Thus, the aim of the present study was to identify and validate immune-related biomarkers and analyze potential regulators including transcription factors (TFs), microRNAs (miRNAs), and drugs for DKD.
Methods
Immune-related genes from the ImmPort database and glomeruli samples from GSE1009 and GSE30528 were used to identify differentially expressed immune-related genes (DEIRGs) of DKD. The expression level and clinical correlation analyses of DEIRGs were verified in the Nephroseq database. Murine podocytes were cultured to construct the high glucose-induced podocyte injury model. The reliability of the bioinformatics analysis was experimentally validated by RT-qPCR in podocytes. Networks among DEIRGs, regulators, and drugs were constructed to predict potential regulatory mechanisms for DKD.
Results
DKD-associated DEIRGs were identified. CCL19 and IL7R were significantly upregulated in the DKD group and negatively correlated with glomerular filtration rate (GFR). GHR, FGF1, FYN, VEGFA, F2R, TGFBR3, PTGDS, FGF9, and SEMA5A were significantly decreased in the DKD group and positively correlated with GFR. RT-qPCR showed that the relative mRNA expression levels of GHR, FGF1, FYN, TGFBR3, PTGDS, FGF9, and SEMA5A were significantly down-regulated in the high glucose-induced podocyte injury group. The enriched regulators for DEIRGs included 110 miRNAs and 8 TFs. The abnormal expression of DEIRGs could be regulated by 16 established drugs.
Conclusions
This study identified immune-related biomarkers, regulators, and drugs of DKD. The findings of the present study provide novel insights into immune-related diagnosis and treatment of DKD.
Introduction
Diabetic kidney disease (DKD) is a major complication of diabetes and the leading cause of end-stage renal disease worldwide [1] which places an enormous burden on healthcare systems. An increasing number of studies suggest that renal inflammation and infiltration of immune cells contribute to the development and progression of DKD [2]. Therefore, immune-related genes (IRGs) play critical roles in the activation and recruitment of inflammatory and immune cells. Considering the existing treatments for DKD are not fully satisfactory, more investigations on the novel biomarkers, mechanisms, and therapeutic targets are needed to slow the progression of DKD. Targeting immune-related biomarkers may be a promising novel therapeutic approach for DKD. Due to this, identifying IRGs and elucidating their associated regulators is helpful to develop new therapeutic targets for DKD. Transcription factors (TFs) and microRNAs (miRNAs) are prominent regulators of gene expression [3]. In addition, drugs can bind to target genes and influence the biological process and expression of target genes [4]. Based on these theoretical foundations, if a TF, a miRNA, or a drug oppositely modulated the expression of IRGs involved in DKD, that TF, miRNA, or drug may be used for the treatment of DKD.
Thereby, in the current study, we aimed to identify and validate immune-related biomarkers for DKD. Glomerular podocytes are the main component in maintaining the glomerular normal filtration barrier. Therefore, podocytes are the initial target of cellular damage in the incidence of DKD [5]. And the injury and loss of podocytes are important features of DKD [6]. Therefore, we established the high glucose-induced podocyte injury model and performed reverse transcription-quantitative polymerase chain reaction (RT-qPCR) to further verify the credibility of our in silico findings. And the potential regulators, including TFs, miRNAs, and drugs associated with these experimentally validated biomarkers were explored. The findings of the present study may provide potential novel insights into immune-related diagnosis and treatment of DKD.
Materials and methods
Data source
The Gene Expression Omnibus (GEO) database is an international public repository of high-throughput gene expression datasets [7]. DKD-related glomerular gene expression profiles (GSE1009 and GSE30528) were downloaded from the GEO DataSets (https://www.ncbi.nlm.nih.gov/geo/). The platform annotation file was used to convert the probe expression matrix into a gene expression matrix. The microarray data of GSE1009 was based on the GPL8300 annotation platform, and GSE30528 was based on the GPL571 platform. The Affymetrix data of GSE1009 contained glomerular samples from 3 healthy controls and 3 DKD patients. GSE30528 dataset encompassed glomerular samples from 9 DKD patients and 13 healthy controls.
ImmPort is an open repository of human immunology database [8]. A list of IRGs was downloaded from the ImmPort database (https://www.immport.org/home).
Screening for differentially expressed immune-related genes (DEIRGs)
The Bioconductor R package “limma” [9] was used to obtain differentially expressed genes (DEGs) between DKD and normal controls in GSE1009 and GSE30528 datasets. Genes with |Log fold change (FC)| ≥ 1 and adjusted P value < 0.05 were defined as DEGs. Then the intersection of the DEGs and the IRGs was selected as DEIRGs. The Venn diagram was created with an online tool jvenn [10].
The clinical significance of DEIRGs
The clinical significance of these DEIRGs in DKD patients was validated via the Woroniecka Diabetes Glom dataset in Nephroseq v5 online database (http://v5.nephroseq.org). The comparison of DEIRGs expression levels between DKD and healthy control group was performed by unpaired t-test (data with normal distribution and equal variance) and Mann–Whitney U test for non-normally distributed data. The relationship between DEIRGs and glomerular filtration rate (GFR) was analyzed by Pearson’s correlation analysis. P < 0.05 (two-tailed) was considered as statistically significant.
Cell culture
Mouse Podocyte Clone-5 (MPC5) cell lines were purchased from the American type culture collection cell bank. MPC5 cell lines were routinely cultured in Dulbecco’s modified Eagle’s medium (DMEM) containing 10% fetal bovine serum, 100 U/ml penicillin, 100 μg/ml streptomycin, and 0.25 μg/ml amphotericin B and maintained at 37 °C incubator containing 5% CO2. The culture medium was changed every 2 days. When 80% confluence was reached, the podocytes were divided into the normal control group (NC group, 30 mmol/l mannitol as an osmotic control) and the high-glucose group (HG group, 30 mmol/l glucose). Mannitol and glucose powder were purchased from Solarbio Technology Co.Ltd and dissolved in DMEM for preparing the working stock concentration of 30 mmol/l. The cells were harvested 72 h after treatment.
RT-qPCR
Total RNA was extracted from the cells in NC and HG groups by using SteadyPure Universal RNA Extraction Kit II (Accurate Biotechnology Co. Ltd.) according to the manufacturer’s instructions. The complementary DNA was synthesized using Evo M-MLV RT Premix for qPCR AG (Accurate Biotechnology Co. Ltd.). The reverse transcription was performed on a Bio-Rad S1000TM Thermal Cycler (Bio-Rad Laboratories). RT-qPCR was performed on the Bio-Rad CFX-96 Real-time PCR system (Bio-Rad Laboratories) using the SYBR® Green Premix Pro Taq HS qPCR Kit (Accurate Biotechnology Co. Ltd.). The standard amplification parameters for PCR consisted of an initial denaturation step (95 °C for 30 s), followed by 40 cycles of 95 °C for 5 s and annealing at 60 °C for 30 s. All primers for PCR amplification were purchased from Shanghai Sangon Biotechnology and the primer sequences are shown in Table 1. mRNA expression was normalized to beta-actin expression. Then, the relative expression levels of genes were calculated using the 2−ΔΔCT method. All samples were performed in triplicate and the experimental data were presented as mean ± standard deviation. Between-group comparisons were analyzed by unpaired t test. P < 0.05 (two-tailed) was considered as statistically significant. DEIRGs validated by RT-qPCR were defined as experimentally validated DEIRGs.
Construction of the experimentally validated target DEIRGs-miRNAs-TFs network
miRWalk database is an online platform which can be used to predict the interaction between miRNA and mRNA [11]. TargetScan is an experimentally validated miRNA-mRNA interaction database [12]. miRDB is an online database for the prediction of miRNA targets [13]. The validated DEIRGs were uploaded to miRWalk database. To increase the reliability of predictions, the intersection portion of the TargetScan and miRDB database was selected as predicted miRNAs. NetworkAnalyst [14] is a platform that integrates resources from ChEA [15] and JASPAR [16] databases to make predictions about TFs. TFs simultaneously enriched in ChEA and JASPAR databases were recorded. The interaction network between target genes, miRNAs, and TFs was constructed and visualized by Cytoscape software.
Identification of the drugs target to the experimentally validated DEIRGs
The Comparative Toxicogenomics Database (CTD) [17] is an open online database that can effectively predict the correlation between diseases, drugs, and genes. “Diabetic nephropathies” and “Diabetic kidney diseases” were entered into the CTD database to obtain Chemical-Gene Interactions. Drugs or chemicals that can influence the expression of the experimentally validated DEIRGs were selected to construct the drug-gene interaction network by using Cytoscape software.
Results
Identification of DEIRGs
According to the defined criteria, 1228 and 345 DEGs were extracted from GSE1009 and GSE30528 respectively. Subsequently, the intersection of DEGs and IRGs was performed by the Venn diagram to obtain DEIRGs. Altogether, 16 DKD-associated DEIRGs were identified including 3 common upregulated DEIRGs (CCL19, IL7R, and TRBC1), 9 common downregulated DEIRGs (GHR, FGF1, FYN, VEGFA, F2R, TGFBR3, PTGDS, FGF9, and SEMA5A), and 4 changed genes (Fig. 1). The detailed comparison of upregulated and downregulated DEIRGs in GSE1009 and GSE30528 can be seen in Additional file 1: Fig. S1. The common upregulated and downregulated DEIRGs were selected as key DEIRGs for further analysis.
Validation of the clinical significance of DEIRGs
The expression profiles of TRBC1 could not be retrieved in the Nephroseq database. CCL19, IL7R, GHR, FGF1, FYN, VEGFA, F2R, TGFBR3, PTGDS, FGF9, and SEMA5A could be analyzed in the Nephroseq database. The expression levels of CCL19 and IL7R were significantly upregulated in the glomerular tissues of DKD patients compared with healthy living donors (Fig. 2a, b). The expression levels of GHR, FGF1, FYN, VEGFA, F2R, TGFBR3, PTGDS, FGF9, and SEMA5A were significantly downregulated in the DKD glomerular tissues compared with healthy kidney samples (Fig. 2c–k).
The results of correlation analysis showed that CCL19 and IL7R were negatively correlated with GFR (Fig. 3a, b) while GHR, FGF1, FYN, VEGFA, F2R, TGFBR3, PTGDS, FGF9, and SEMA5A were positively correlated with GFR (Fig. 3c–k). Ultimately, the clinical significance of CCL19, IL7R, GHR, FGF1, FYN, VEGFA, F2R, TGFBR3, PTGDS, FGF9, and SEMA5A were validated.
RT-qPCR validation of differential gene expression
Based on the clinical significance of DEIRGs, CCL19, IL7R, GHR, FGF1, FYN, VEGFA, F2R, TGFBR3, PTGDS, FGF9, and SEMA5A were selected to perform RT-qPCR to further validate their expression in podocytes. The results showed that the relative mRNA expression levels of GHR, FGF1, FYN, TGFBR3, PTGDS, FGF9, and SEMA5A were significantly downregulated in the HG group compared with the NC group (Fig. 4a–k). Therefore, GHR, FGF1, FYN, TGFBR3, PTGDS, FGF9, and SEMA5A were defined as experimentally validated DEIRGs.
Construction of the predicted miRNAs-experimentally validated DEIRGs-TFs network
The enriched miRNAs in TargetScan and miRDB databases were intersected and 110 common miRNAs were recorded. Similarly, the enriched TFs in ChEA and JASPAR databases were intersected and 8 common TFs were recorded. The network among miRNAs, target genes, and TFs can be seen in Fig. 5.
Identification of the drugs target to experimentally validated DEIRGs
The CTD database was used to construct the network between drugs and the experimentally validated DEIRGs associated with DKD. The results showed that various drugs or chemicals such as estradiol, and resveratrol could regulate the abnormal mRNA expression levels of validated DEIRGs (Fig. 6).
Discussion
Accumulating evidence from experimental and clinical studies has demonstrated that immune response is associated with the progression of DKD [18]. In the present study, a total of 12 DEIRGs associated with DKD were identified, including 3 upregulated DEIRGs (CCL19, IL7R, and TRBC1) and 9 downregulated DEIRGs (GHR, FGF1, FYN, VEGFA, F2R, TGFBR3, PTGDS, FGF9, and SEMA5A). The results in the Woroniecka Diabetes Glom dataset of Nephroseq database also proved that the expression levels of CCL19 and IL7R were significantly increased, whereas the expression levels of GHR, FGF1, FYN, VEGFA, F2R, TGFBR3, PTGDS, FGF9, and SEMA5A were significantly decreased in the DKD group. In parallel, the up-regulation of CCL19, IL7R and the down-regulation of GHR, FGF1, FYN, VEGFA, F2R, TGFBR3, PTGDS, FGF9, and SEMA5A were correlated with the deterioration of renal function. Thus, CCL19, IL7R, GHR, FGF1, FYN, VEGFA, F2R, TGFBR3, PTGDS, FGF9, and SEMA5A were hub DEIRGs verified in clinical significance. To further confirm the results from bioinformatics analysis, we performed corresponding in vitro experiments.
And our experimental results confirmed that the relative mRNA expression levels of GHR, FGF1, FYN, TGFBR3, PTGDS, FGF9, and SEMA5A were significantly down-regulated in the HG group. Thus, these experimentally validated DEIRGs could be therapeutic targets for DKD. It has been well established that the activation of the immune inflammatory response can lead to podocyte injury and disrupt the integrity of the glomerular filtration barrier [19] and subsequently induce further deterioration of the renal function and glomerular sclerosis [20]. Consequently, changes in DEIRGs may be useful for timely recognition and treatment of DKD. Therefore, we also explored miRNAs, TFs, and drugs that may regulate the expression levels of DEIRGs.
It has been shown that podocytes express GHR when exposed to growth hormone (GH) [21]. Past research has revealed that under the condition of different chronic kidney diseases, including DKD, the expression level of renal GHR was decreased both in humans and in rodents. And the increase of GHR level was associated with the increase of GFR and the decrease of serum creatinine [22]. However, there are also a few converse opinions. A previous study indicated that chronic exposure to high doses of GH may induce glomerulosclerosis and increased proteinuria [23]. Another study also showed that increased renal expression of the GHR was associated with nephropathy in poorly controlled type 1 diabetes [21]. Therefore, GH/GHR axis involved in the pathogenesis of DKD and maintaining the normal expression level of GHR may be a therapeutic target for DKD.
FGF1 is a well-known mitogen and insulin sensitizer. A previous study has observed a significant reduction in the renal expression of FGF1 in patients with DKD and in mouse models. And FGF1 treatment can prevent the development of DKD by inhibiting the activation of inflammatory signaling cascades such as nuclear factor-κB and c-Jun N-terminal kinases signaling [24]. Meanwhile, FGF1 treatment can also increase insulin sensitization, maintain normoglycemia, and prevent diabetic complications [25]. Similarly, FGF9 can induce ERK and AKT phosphorylation and cAMP-response element binding protein and nuclear factor erythroid-derived 2-like 2 (Nrf2) activation to play the antioxidative function [26] and cardioprotective effect [27]. Our results also proved that the expression levels of FGF1 and FGF9 were decreased in DKD status. Considering our results together with past findings, we can speculate that the down-regulation of FGF1 and FGF9 indicates a weakening antioxidant capacity in DKD status, making the kidney vulnerable to oxidative stress injury. Therefore, FGF1 and FGF9 could be potential therapeutic targets for DKD.
FYN is a proto-oncogene. Previous research has shown that the lipid lowering drug-fenofibrate can protect against DKD by activating Akt2/GSK-3β/FYN/Nrf2 antioxidants [28]. Another study also proved that sulforaphane can ameliorate DKD via GSK-3β/FYN/Nrf2 signaling pathway [29]. Therefore, FYN maybe a promising target in the treatment of DKD. Our results from gene-targeted drugs showed that sulforaphane can also target PTGDS to play a protective role in DKD. PTGDS was downregulated by DNA methylation [30]. Therefore, the downregulation of PTGDS in DKD represents an increased DNA methylation which is closely associated with aging and senescence [31]. Thus, improving the expression level of PTGDS by altering DNA methylation may not only contribute to the control of DKD but also help in delaying or alleviating the consequences of aging.
TGFBR3 is a membrane proteoglycan that functions as a co-receptor with other transforming growth factor receptors. Soluble TGFBR3 may inhibit transforming growth factor-β (TGF-β) signaling [32]. It has been extensively shown that TGF-β is highly activated in glomeruli, tubules, and tubulointerstitium in both patients and animal models with DKD [33]. Hyperglycemia and the formation of advanced glycation end products (AGEs), angiotensin II, and reactive oxidative species induced by high glucose are capable of activating TGF-β [33]. And the activation of TGF-β is closely associated with renal fibrosis and the deterioration of renal function [34]. Thus, the decrease of TGFBR3 proved by the present study may lead to a decrease in its inhibition on TGF-β, which would further increase the expression level of TGF-β and aggravate renal impairment. Therefore, the current study supports the loss of TGFBR3 as one mechanism for the development of DKD.
SEMA5A is a member of the semaphorin family which involved in neuronal and vascular development [35]. SEMA5A can also increase endothelial cell proliferation and decrease cell apoptosis [35]. However, the role of SEMA5A in the occurrence and development of DKD has not been fully studied. We proved that the expression level of SEMA5A was significantly downregulated in DKD glomerulus by bioinformatics analysis. And our subsequent experimental findings confirmed the results predicted by bioinformatics analysis. Considering the effects of SEMA5A in endothelial cell proliferation, migration, apoptosis, neuronal and vascular development, we can speculate that the reduced expression of SEMA5A may involve not only in the progression of diabetic microvascular complications such as DKD, but also in other diabetic complications such as diabetic neuropathy and cognitive dysfunction. Therefore, regulation of SEMA5A may help alleviate diabetic complications.miRNA is an important component of the gene regulatory network, which can regulate the expression of the target gene at the post-transcriptional level [36]. TFs are sequence-specific DNA binding proteins capable of regulating transcription [37]. We identified 110 miRNAs and 8 TFs as regulators of DEIRGs. We could speculate that these regulators participate in the pathological process of DKD and targeting these regulators may help discover potential protective miRNAs or TFs in the future.
Gene-targeted drugs may affect the expression of DEIRGs to retard disease progression. To predict the potential effective drugs for DKD, we applied the CTD database to explore therapeutic agents that might reverse the abnormal expression of DEIRGs in DKD. After exploring the interactions between validated DEIRGs and drugs associated with DKD, we found that 16 drugs could influence the expression of DEIRGs. Among these drugs, estradiol and resveratrol may occupy the most important positions because they can regulate four and three genes respectively. Estradiol is one of the most important sex steroid hormones secreted by the ovary [38]. Many studies have found that premenopausal women have a lower risk for developing DKD [39] because estradiol can attenuate glomerulosclerosis and tubulointerstitial fibrosis by reducing the synthesis of type I and type IV collagen [40,41,42]. Meanwhile, reduced concentrations of plasma estradiol exhibit a greater incidence and progression of DKD [39, 43]. And supplementation of 17β-estradiol markedly inhibited the progression of DKD [40, 44]. Our network pharmacology analysis showed that estradiol can improve the expression of abnormally down-regulated DEIRGs such as GHR, TGFBR3, FYN, and FGF9. Thus, estradiol can inhibit the occurrence or development of DKD by acting on at least 4 gene targets. Our findings also showed that resveratrol can prevent the development of DKD by reversing downregulated DEIRGs such as FGF1, PTGDS, and GHR. In essence, the chemical structure of resveratrol is similar to 17β-estradiol [45]. Consequently, resveratrol is structurally and functionally similar to estrogen [46]. Animal experiments have proved that resveratrol can prevent the progression of DKD by improving gut environment, reducing inflammatory response [47], improving kidney function, relieving proteinuria, ameliorating lipid metabolism, and inhibiting apoptosis [48]. Both resveratrol and 17β-estradiol can significantly reduce blood glucose, and protect pancreatic islet cell [46]. For the above reasons, we support the notion that proper 17β-estradiol intervention is an effective option in the management of DKD. It should also be noted that further experiments are needed to support whether DKD patients can benefit from DEIRGs intervention.
In conclusion, we identified 11 immune-related biomarkers for DKD, including 2 up-regulated DEIRGs (CCL19, IL7R) and 9 down-regulated DEIRGs (GHR, FGF1, FYN, VEGFA, F2R, TGFBR3, PTGDS, FGF9, SEMA5A) by bioinformatics analysis. The up-regulation of CCL19, IL7R and the down-regulation of GHR, FGF1, FYN, VEGFA, F2R, TGFBR3, PTGDS, FGF9, and SEMA5A were correlated with the deterioration of renal function. RT-qPCR showed that the relative mRNA expression levels of GHR, FGF1, FYN, TGFBR3, PTGDS, FGF9, and SEMA5A were significantly down-regulated in the high glucose-induced podocyte injury model. The enriched regulators for DEIRGs included 110 miRNAs and 8 TFs. Estradiol and resveratrol were the most enriched immune-related drugs in the treatment of DKD. The current findings may provide novel potential therapeutic targets for the precise diagnosis and immune-related therapies of DKD.
Availability of data and materials
The datasets analysed during the current study are available in the GEO and Nephroseq repository [Accession numbers: GSE1009 and GSE30528].
Abbreviations
- CCL19 :
-
C-C motif chemokine ligand 19
- IL7R :
-
Interleukin 7 receptor
- TRBC1 :
-
T cell receptor beta constant 1
- GHR :
-
Growth hormone receptor
- FGF1 :
-
Fibroblast growth factor 1
- FYN :
-
FYN proto-oncogene, Src family tyrosine kinase
- VEGFA :
-
Vascular endothelial growth factor A
- F2R :
-
Coagulation factor II thrombin receptor
- TGFBR3 :
-
Transforming growth factor beta receptor 3
- PTGDS :
-
Prostaglandin-H2 D-isomerase
- FGF9 :
-
Fibroblast growth factor 9
- SEMA5A :
-
Semaphorin 5A
References
Zoja C, Xinaris C, Macconi D. Diabetic nephropathy: novel molecular mechanisms and therapeutic targets. Front Pharmacol. 2020;11:586892.
Tang SCW, Yiu WH. Innate immunity in diabetic kidney disease. Nat Rev Nephrol. 2020;16(4):206–22.
Hobert O. Gene regulation by transcription factors and microRNAs. Science. 2008;319(5871):1785–6.
Li S, Cao Y, Li L, Zhang H, Lu X, Bo C, Kong X, Liu Z, Chen L, Liu P, Jiao Y, Wang J, Ning S, Wang L. Building the drug-GO function network to screen significant candidate drugs for myasthenia gravis. PLoS ONE. 2019;14(4):e0214857.
Zhang T, Chi Y, Ren Y, Du C, Shi Y, Li Y. Resveratrol reduces oxidative stress and apoptosis in podocytes via Sir2-related enzymes, sirtuins1 (SIRT1)/peroxisome proliferator-activated receptor γ co-activator 1α (PGC-1α) axis. Med Sci Monit. 2019;25:1220–31.
Hu M, Wang R, Li X, Fan M, Lin J, Zhen J, Chen L, Lv Z. LncRNA MALAT1 is dysregulated in diabetic nephropathy and involved in high glucose-induced podocyte injury via its interplay with β-catenin. J Cell Mol Med. 2017;21(11):2732–47.
Clough E, Barrett T. The gene expression omnibus database. Methods Mol Biol. 2016;1418:93–110.
Bhattacharya S, Dunn P, Thomas CG, Smith B, Schaefer H, Chen J, Hu Z, Zalocusky KA, Shankar RD, Shen-Orr SS, Thomson E, Wiser J, Butte AJ. ImmPort, toward repurposing of open access immunological assay data for translational and clinical research. Sci Data. 2018;5:180015.
Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, Smyth GK. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43(7):e47–e47.
Bardou P, Mariette J, Escudié F, Djemiel C, Klopp C. jvenn: an interactive Venn diagram viewer. BMC Bioinform. 2014;15(1):293.
Dweep H, Sticht C, Pandey P, Gretz N. miRWalk–database: prediction of possible miRNA binding sites by “walking” the genes of three genomes. J Biomed Inform. 2011;44(5):839–47.
Agarwal V, Bell GW, Nam JW, Bartel DP. Predicting effective microRNA target sites in mammalian mRNAs. Elife. 2015;4:e05005.
Chen Y, Wang X. miRDB: an online database for prediction of functional microRNA targets. Nucleic Acids Res. 2020;48(D1):D127-d131.
Zhou G, Soufan O, Ewald J, Hancock REW, Basu N, Xia J. NetworkAnalyst 3.0: a visual analytics platform for comprehensive gene expression profiling and meta-analysis. Nucleic Acids Res. 2019;47(W1):W234–41.
Keenan AB, Torre D, Lachmann A, Leong AK, Wojciechowicz ML, Utti V, Jagodnik KM, Kropiwnicki E, Wang Z. Ma’ayan A: ChEA3: transcription factor enrichment analysis by orthogonal omics integration. Nucleic Acids Res. 2019;47(W1):W212–24.
Castro-Mondragon JA, Riudavets-Puig R, Rauluseviciute I, Lemma RB, Turchi L, Blanc-Mathieu R, Lucas J, Boddie P, Khan A, Manosalva Pérez N, Fornes O, Leung TY, Aguirre A, Hammal F, Schmelter D, Baranasic D, Ballester B, Sandelin A, Lenhard B, Vandepoele K, Wasserman WW, Parcy F, Mathelier A. JASPAR 2022: the 9th release of the open-access database of transcription factor binding profiles. Nucleic Acids Res. 2022;50(D1):D165-d173.
Davis AP, Grondin CJ, Johnson RJ, Sciaky D, Wiegers J, Wiegers TC, Mattingly CJ. Comparative Toxicogenomics Database (CTD): update 2021. Nucleic Acids Res. 2021;49(D1):D1138–43.
Tesch GH. Diabetic nephropathy—is this an immune disorder? Clin Sci (Lond). 2017;131(16):2183–99.
Lazareth H, Henique C, Lenoir O, Puelles VG, Flamant M, Bollée G, Fligny C, Camus M, Guyonnet L, Millien C, Gaillard F, Chipont A, Robin B, Fabrega S, Dhaun N, Camerer E, Kretz O, Grahammer F, Braun F, Huber TB, Nochy D, Mandet C, Bruneval P, Mesnard L, Thervet E, Karras A, Le Naour F, Rubinstein E, Boucheix C, Alexandrou A, Moeller MJ, Bouzigues C, Tharaux P-L. The tetraspanin CD9 controls migration and proliferation of parietal epithelial cells and glomerular disease progression. Nat Commun. 2019;10(1):3303.
Lim AKH, Tesch GH. Inflammation in diabetic nephropathy. Mediat Inflamm. 2012;2012:146154.
Nishad R, Mukhi D, Singh AK, Motrapu M, Chintala K, Tammineni P, Pasupulati AK. Growth hormone induces mitotic catastrophe of glomerular podocytes and contributes to proteinuria. Cell Death Dis. 2021;12(4):342.
Brittain AL, Kopchick JJ. A review of renal GH/IGF1 family gene expression in chronic kidney diseases. Growth Horm IGF Res. 2019;48–49:1–4.
Grunenwald S, Tack I, Chauveau D, Bennet A, Caron P. Impact of growth hormone hypersecretion on the adult human kidney. Ann Endocrinol (Paris). 2011;72(6):485–95.
Liang G, Song L, Chen Z, Qian Y, Xie J, Zhao L, Lin Q, Zhu G, Tan Y, Li X, Mohammadi M, Huang Z. Fibroblast growth factor 1 ameliorates diabetic nephropathy by an anti-inflammatory mechanism. Kidney Int. 2018;93(1):95–109.
Wang D, Zhao T, Zhao Y, Yin Y, Huang Y, Cheng Z, Wang B, Liu S, Pan M, Sun D, Wang Z, Zhu G. PPARγ mediates the anti-epithelial-mesenchymal transition effects of FGF1(ΔHBS) in chronic kidney diseases via inhibition of TGF-β1/SMAD3 signaling. Front Pharmacol. 2021;12:690535.
Chuang JI, Huang JY, Tsai SJ, Sun HS, Yang SH, Chuang PC, Huang BM, Ching CH. FGF9-induced changes in cellular redox status and HO-1 upregulation are FGFR-dependent and proceed through both ERK and AKT to induce CREB and Nrf2 activation. Free Radic Biol Med. 2015;89:274–86.
Sun J, Wang Z, Shi H, Gu L, Wang S, Wang H, Li Y, Wei T, Wang Q, Wang L. LncRNA FAF inhibits fibrosis induced by angiotensinogen II via the TGFβ1-P-Smad2/3 signalling by targeting FGF9 in cardiac fibroblasts. Biochem Biophys Res Commun. 2020;521(3):814–20.
Cheng Y, Zhang X, Ma F, Sun W, Wang W, Yu J, Shi Y, Cai L, Xu Z. The role of Akt2 in the protective effect of fenofibrate against diabetic nephropathy. Int J Biol Sci. 2020;16(4):553–67.
Shang G, Tang X, Gao P, Guo F, Liu H, Zhao Z, Chen Q, Jiang T, Zhang N, Li H. Sulforaphane attenuation of experimental diabetic nephropathy involves GSK-3 beta/Fyn/Nrf2 signaling pathway. J Nutr Biochem. 2015;26(6):596–606.
Xue SS, He JL, Zhang X, Liu YJ, Xue FX, Wang CJ, Ai D, Zhu Y. Metabolomic analysis revealed the role of DNA methylation in the balance of arachidonic acid metabolism and endothelial activation. Biochim Biophys Acta. 2015;1851(10):1317–26.
Reynolds CA, Tan Q, Munoz E, Jylhävä J, Hjelmborg J, Christiansen L, Hägg S, Pedersen NL. A decade of epigenetic change in aging twins: genetic and environmental contributions to longitudinal DNA methylation. Aging Cell. 2020;19(8):e13197.
Álvarez-Silva MC, Yepes S, Torres MM, Barrios AF. Proteins interaction network and modeling of IGVH mutational status in chronic lymphocytic leukemia. Theor Biol Med Model. 2015;12:12.
Wang L, Wang HL, Liu TT, Lan HY. TGF-beta as a master regulator of diabetic nephropathy. Int J Mol Sci. 2021;22(15):7881.
Hong Q, Cai H, Zhang L, Li Z, Zhong F, Ni Z, Cai G, Chen XM, He JC, Lee K. Modulation of transforming growth factor-β-induced kidney fibrosis by leucine-rich ⍺-2 glycoprotein-1. Kidney Int. 2022;101(2):299–314.
Sadanandam A, Rosenbaugh EG, Singh S, Varney M, Singh RK. Semaphorin 5A promotes angiogenesis by increasing endothelial cell proliferation, migration, and decreasing apoptosis. Microvasc Res. 2010;79(1):1–9.
Chen W, Chen H, Zheng D, Zhang H, Deng L, Cui W, Zhang Y, Santos HA, Shen H. Gene-hydrogel microenvironment regulates extracellular matrix metabolism balance in nucleus pulposus. Adv Sci (Weinh). 2019;7(1):1902099–1902099.
Gao Y, Zhu N, Zhu X, Wu M, Jiang C-Z, Grierson D, Luo Y, Shen W, Zhong S, Fu D-Q, Qu G. Diversity and redundancy of the ripening regulatory networks revealed by the fruitENCODE and the new CRISPR/Cas9 and mutants. Hortic Res. 2019;6:39.
Cheng J-C, Fang L, Yan Y, He J, Guo Y, Jia Q, Gao Y, Han X, Sun Y-P. TGF-β1 stimulates aromatase expression and estradiol production through SMAD2 and ERK1/2 signaling pathways in human granulosa-lutein cells. J Cell Physiol. 2021;236(9):6619–29.
Maric C, Sullivan S. Estrogens and the diabetic kidney. Gend Med. 2008;5(Suppl A):S103–13.
Inada A, Inada O, Fujii NL, Nagafuchi S, Katsuta H, Yasunami Y, Matsubara T, Arai H, Fukatsu A, Nabeshima Y-I. Adjusting the 17β-estradiol-to-androgen ratio ameliorates diabetic nephropathy. J Am Soc Nephrol. 2016;27(10):3035–50.
Al-Trad B, Ashankyty IM, Alaraj M. Progesterone ameliorates diabetic nephropathy in streptozotocin-induced diabetic Rats. Diabetol Metab Syndr. 2015;7:97.
Xu H, Wu X, Qin H, Tian W, Chen J, Sun L, Fang M, Xu Y. Myocardin-related transcription factor a epigenetically regulates renal fibrosis in diabetic nephropathy. J Am Soc Nephrol. 2015;26(7):1648–60.
Wells CC, Riazi S, Mankhey RW, Bhatti F, Ecelbarger C, Maric C. Diabetic nephropathy is associated with decreased circulating estradiol levels and imbalance in the expression of renal estrogen receptors. Gend Med. 2005;2(4):227–37.
Catanuto P, Doublier S, Lupia E, Fornoni A, Berho M, Karl M, Striker GE, Xia X, Elliot S. 17 beta-estradiol and tamoxifen upregulate estrogen receptor beta expression and control podocyte signaling pathways in a model of type 2 diabetes. Kidney Int. 2009;75(11):1194–201.
Parazzini F. Resveratrol, tryptophanum, glycine and vitamin E: a nutraceutical approach to sleep disturbance and irritability in peri- and post-menopause. Minerva Ginecol. 2015;67(1):1–5.
Li Y, Huang J, Yan Y, Liang J, Liang Q, Lu Y, Zhao L, Li H. Preventative effects of resveratrol and estradiol on streptozotocin-induced diabetes in ovariectomized mice and the related mechanisms. PLoS ONE. 2018;13(10):e0204499.
Cai T-T, Ye X-L, Li R-R, Chen H, Wang Y-Y, Yong H-J, Pan M-L, Lu W, Tang Y, Miao H, Snijders AM, Mao J-H, Liu X-Y, Lu Y-B, Ding D-F. Resveratrol modulates the gut microbiota and inflammation to protect against diabetic nephropathy in mice. Front Pharmacol. 2020;11:1249.
Zhu H, Zhong S, Yan H, Wang K, Chen L, Zhou M, Li Y. Resveratrol reverts Streptozotocin-induced diabetic nephropathy. Front Biosci (Landmark Ed). 2020;25:699–709.
Acknowledgements
We greatly acknowledge all the researchers and subjects for providing public data. Animals are not involved in our present study.
Funding
This work was supported by National Natural Science Foundation of China [No. 82170755].
Author information
Authors and Affiliations
Contributions
All authors contributed to the study conception and design. SNC and HLJ conceived the idea and organized the study. Data collection and analysis were performed by SNC, BL, and LC. Cell culture and model-building were completed by SNC and BL. RT-qPCR was performed by SNC. The first draft of the manuscript was written by SNC and revised by LC and HLJ. All authors commented on previous versions of the manuscript. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The author has no competing interests in this manuscript.
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.
Differentially expressed genes in GSE1009 and GSE30528.
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.
About this article
Cite this article
Chen, S., Li, B., Chen, L. et al. Identification and validation of immune-related biomarkers and potential regulators and therapeutic targets for diabetic kidney disease. BMC Med Genomics 16, 90 (2023). https://doi.org/10.1186/s12920-023-01519-6
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12920-023-01519-6