Abstract Background Diabetic nephropathy (DN) is a prevalent complication of diabetes, often leading to end-stage kidney disease (ESKD). Advanced DN progresses to ESKD rapidly, yet effective diagnostic indicators and treatments are lacking. Methods Two DN-related datasets were obtained from the Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) were identified using the R packages. Pyroptosis-related genes (PRGs) and ferroptosis-related genes (FRGs) were collected from their respective database. Pyroptosis- and ferroptosis-related differentially expressed genes (PFRDEGs) were identified by overlapping DEGs, PRGs, and FRGs for further analysis, including functional enrichment and immune infiltration. Hub genes were identified using a PPI network via MCODE-plugin in Cytoscape. GeneMANIA was utilized to explore intermolecular interactions among hub genes. Based on these hub genes, a diagnostic model was constructed using the receiver operating characteristic curve and potential therapeutic agents were retrieved. Correlation analysis between hub genes and estimated glomerular filtration rate was performed using Nephroseq v5 database, and expression of hub genes was validated in external GEO database, Nephroseq v5 database and DN mice in vivo. Results Four hub genes (CYBB, LCN2, JUN and ADIPOQ) were identified, and three of the four hub genes (CYBB, LCN2 and ADIPOQ) were found to be potential biomarkers for advanced DN. On this basis, three potential therapeutic agents were screened. More importantly, a series of biological experiments confirmed that CYBB and LCN2 were significantly up-regulated in DN mice. Conclusion This study identifies three hub genes as diagnostic biomarkers and mines three potential therapeutic agents for advanced DN, providing new insights into the role of pyroptosis and ferroptosis in advanced DN and laying the foundation for future research. Keywords: diabetic nephropathy, pyroptosis, ferroptosis, bioinformatics, diagnostic biomarkers Introduction As an urgent challenge for public health,[32]^1 diabetic nephropathy (DN), also known as diabetic kidney disease, is a significant chronic microvascular complication of diabetes mellitus.[33]^2 Globally, DN is the leading cause of CKD [including end-stage kidney disease (ESKD)] and an important risk factor for cardiovascular disease and early death in diabetic patients.[34]^3^,[35]^4 DN can be categorized into two groups as early DN and advanced DN according to the levels of urinary albumin-to-creatinine ratio (UACR) and renal function, and advanced DN was defined as UACR > 300mg/g or estimated glomerular filtration rate (eGFR) < 90mL/min.[36]^5 Compared with early DN, patients with advanced DN have a higher risk of declining renal function and progression to ESKD.[37]^6 Nevertheless, there is currently a lack of effective diagnostic means and treatments for advanced DN. In terms of diagnostics, proteinuria and eGFR are nowadays crucial indicators for diagnosing DN, but they have their limitations. For example, among diabetic patients with massive proteinuria, approximately 70% of renal puncture biopsies confirm the presence of DN,[38]^7^,[39]^8 and only 40% of patients with type 2 diabetes and microalbuminuria exhibit the typical pathology of DN.[40]^9 Furthermore, an increasing number of studies have discovered that diabetic patients with normal urinary albumin can still progress to DN, as evidenced by a decrease in eGFR to less than 60mL·min^−1 (1.73 m^2)^−1 and the presence of typical DN in renal puncture pathology.[41]^8^,[42]^10 Above all, there is the absence of effective measures to diagnostic advanced DN. In terms of treatment, although new drugs such as SGLT-2i, DPP-4 inhibitors, and finerenone have provided additional treatment options for DN,[43]^11 many patients with DN still progress to ESKD due to the limited effectiveness of the aforementioned drugs. Consequently, in-depth research on the pathogenesis of DN for identifying effective diagnostic biomarkers and exploring corresponding therapeutic agents will help to address these research gaps. As a complex process involving multiple factors and mechanisms, the occurrence and development of DN is closely linked to the cell death.[44]^12 Ferroptosis, one of the prominent research hotspots in cell death studies, was first proposed by Stockwell et al in 2012.[45]^13 Nowadays increasing evidence suggests that ferroptosis is a key driver of DN.[46]^14–16 Moreover, recent studies have demonstrated the efficacy of various small molecular agents, such as N-acetylcysteine, dapagliflozin, and salusin-β, in improving renal function and alleviating kidney injury by suppressing ferroptosis.[47]^14 On the other hand, pyroptosis is another recently identified mechanism of programmed cell death associated with caspase-1, which triggers a series of inflammatory responses through the release of multiple pro-inflammatory factors such as IL-1β and IL-18.[48]^17 Recent studies also have reported that pyroptosis is related to the occurrence and development of DN, and that renal lesions in DN can be ameliorated by inhibiting the pyroptosis pathway.[49]^17–19 As important indicators of pyroptosis, caspase-11 and gasdermin D (GSDMD) were found to be significantly overexpressed in mice with DN, accompanied by significantly elevated plasma levels of inflammatory factors such as NF-κB, IL-1β, and IL-18, which were downregulated by knockdown of CASPASE-11 or GSDMD.[50]^20 Furthermore, numerous findings have demonstrated that targeting the pertinent pathways of pyroptosis can effectively prevent the occurrence or delay the progression of DN. For instance, Pyrroloquinoline quinone ameliorates renal fibrosis in DN by inhibiting the pyroptosis pathway in C57BL/6 mice and human kidney-2 cells;[51]^18 Inhibition of cellular death through targeted inhibition of caspase-1 significantly attenuates podocyte injury in DN mice;[52]^21 and deletion of Gasdermin D significantly reduces the release of IL-1β and IL-18, thereby attenuating cellular pyroptosis and renal injury in diabetic mice in vivo.[53]^22 More importantly, A growing body of research suggests that multiple layers of interconnections exist between pyroptosis and ferroptosis, including shared triggers, molecular components and protective mechanisms.[54]^23–27 More importantly, many recent studies have found that there is a mechanistic crosstalk between pyroptosis and ferroptosis, which provides new perspectives on the mechanistic study of the disease and new avenues for the treatment of the disease.[55]^28–31 We therefore performed this study to explore potential diagnostic biomarkers linked to advanced DN from the viewpoints of both pyroptosis and ferroptosis by employing bioinformatics techniques. Besides, we also predicted the potential therapeutic agents targeting these biomarkers in an attempt to discover effective treatment measures for advanced DN. It is expected to provide new research perspectives for precise diagnosis and treatment for advanced DN. Material and Methods Data Source For this study, we retrieved [56]GSE142025 from the Gene Expression Omnibus (GEO) database that met the specified conditions using the keywords “diabetic nephropathy” or “diabetic kidney disease”. We then filtered out the dataset [57]GSE30528 using the keywords “diabetic nephropathy” or “diabetic kidney disease” to meet the purpose of the study. The transcriptome datasets [58]GSE142025 and [59]GSE30528 were obtained from the National Center for Biotechnology Information (NCBI) / GEO database. [60]GSE142025 was based on the Illumina HiSeq 4000 platform (Homo sapiens [61]GPL20301).[62]^5^,[63]^32 Used as the primary objects of study, 21 advanced DN samples ([64]GSM4217781–[65]GSM4217801) and 9 control samples ([66]GSM4217808–[67]GSM4217816) were selected from the [68]GSE142025 dataset, which was employed as the training dataset in the study ([69]Supplementary Table S1 presents the data of [70]GSE142025 in detail). Additionally, Based on the [71]GPL571 platform of Affymetrix Human Genome U133A 2.0 Array,[72]^33^,[73]^34 [74]GSE30528 includes 9 DN glomeruli samples ([75]GSM756992–[76]GSM757000) and 13 control glomeruli samples ([77]GSM757001–[78]GSM757013), and it was utilized as an external validation dataset ([79]Supplementary Table S2 presents the data of [80]GSE30528 in detail). Clinical studies associated with [81]GSE142025 was approved by the institutional review board at Shanghai Jiao Tong University Affiliated Sixth People’s Hospital. Clinical studies associated with [82]GSE30528 was approved by the institutional review board of the Albert Einstein College of Medicine and Montefiore Medical Center (2002–202 to K.S). What’s more, pyroptosis-related genes (PRGs) with protein encoding function were retrieved from the Genecards website ([83]https://www.proteinatlas.org/, version 5.19) and the GSEA/Molecular Signatures Database ([84]https://www.gsea-msigdb.org/gsea/msigdb/, version 4.3.3), and [85]Supplementary Table S3 presents the PRGs in detail. Ferroptosis-related genes (FRGs) were downloaded from FerrDb ([86]http://www.zhounan.org/ferrdb/current/) and the GSEA/ Molecular Signatures Database ([87]https://www.gsea-msigdb.org/gsea/msigdb/), [88]Supplementary Table S4 presents the FRGs in detail. Microarray Data Processing The gene probe information was converted into gene symbols using the corresponding annotation profile in [89]GSE30528. In cases where multiple probes corresponded to the same gene, the final gene expression value was determined by calculating the average expression value, and probes that did not map to genes were removed. The download and normalization of the raw matrix for [90]GSE30528 was done with the “GEOquery” (version 2.72.0), “magrittr” (version 2.0.3), and “limma” (version 3.60.3) R packages. The raw matrix for [91]GSE142025 were obtained by manual download and then collated. The standardized plots are provided as [92]Figure S1 in the supplementary materials. Screening the Differentially Expressed Genes (DEGs) and the Pyroptosis- and Ferroptosis-Related Differentially Expressed Genes (PFRDEGs) The [93]GSE142025 transcriptome dataset was analyzed using the “limma” R package (version 3.60.3). DEGs were identified based on adjusted P value < 0.05 and |log fold change(FC)| > 1. Heatmaps and volcano plots were used to visually represent the DEGs. The overlap among DEGs, PRGs and FRGs was visualized through a Venn diagram, displaying PFRDEGs. Gene Enrichment Analysis Gene Ontology (GO) / Kyoto Encyclopedia of Genes and Genomes (KEGG) Pathway Enrichment Analysis for the PFRDEGs were performed by applying the “clusterProfiler” package for enrichment analysis in xiantao online website ([94]https://www.xiantaozi.com/, version 2.0).[95]^35^,[96]^36 GO functional enrichment analysis include analysis of molecular function (MF), biological process (BP), and cellular component (CC) terms. In this study a gene set was regarded as significantly enriched if P.adj < 0.05. Gene Set Enrichment Analysis (GSEA) was performed to explore the potential biological processes of genes involved in advanced DN. GSEA was also performed by applying the “clusterProfiler” R package. P.adj < 0.05 and FDR (q value) < 0.25 was regarded as significantly enriched. Identifying Hub Genes and Exploring Intermolecular Interactions of PFRDEGs To obtain hub genes for more in-depth analysis of the roles of pyroptosis and ferroptosis in advanced DN, STRING database ([97]https://string-db.org/, version 11.5)[98]^37 was utilized to create a protein–protein Interaction (PPI) network for PFRDEGs. MCODE-plugin of Cytoscape (version 3.9.0) was then employed to identify hub genes within the PPI network [Detailed parameters: Degree cutoff: 2; Cluster Finding: Haircut; Node Score Cutoff: 0.2; K-Core: 2; Max Depth: 100.]. GeneMANIA online software ([99]https://genemania.org/) was utilized to delve intermolecular interactions among the four hub genes. Expression of Hub Genes in [100]GSE142025 and Validation in External Dataset Violin plots were used to visually detect differences in hub gene expression between the advanced DN cohort and the control group. In addition, to reduce the potential for false positives, hub gene expression was confirmed using an independent DN dataset from the GEO database. For this validation, [101]GSE30528 was selected, which contains data from 9 DN patients and 13 controls. Hub gene expression data was then acquired and presented using violin plots to facilitate comparative analysis. Analysis of Immune Infiltration and Correlations Between Hub Genes and Immune Cells The CIBERSORTx algorithm serves as a sophisticated computational framework, adept at translating normalized gene expression data into detailed profiles of immune cell infiltration.[102]^38 Leveraging the LM22 reference expression signature, we meticulously computed the relative frequencies of diverse immune cell types, differentiating between advanced DN patients and the control cohort, with the aid of the “CIBERSORT” R package (version 1.03). To deepen our understanding of the interplay between immune cells and pivotal hub genes, we conducted a thorough Pearson correlation analysis. The ensuing data were meticulously processed and rendered visually through the use of both boxplots and heatmaps, facilitated by an online platform ([103]https://www.xiantaozi.com/, version 2.0).[104]^35^,[105]^36 Exploring the Correlation Between Hub Genes and Clinical Indicators To explore the link between hub genes and clinical practice, we performed a data mining analysis between the expression of hub genes and eGFR utilizing the Nephroseq v5 database. The Nephroseq v5 database ([106]http://v5.nephroseq.org, version 4)[107]^39 is a robust information platform designed to evaluate the correlation between gene expression levels and clinical characteristics of kidney diseases. Additionally, we verified the hub genes expression in our study by cross-referencing it with data from this platform. Screening and Validation of Candidate Diagnostic Biomarkers To further explore the clinical significance of the hub genes, receiver operating characteristic (ROC) curves were constructed to identify the potential diagnostic efficacy of the hub genes, which was evaluated based on the area under the curve (AUC). The ROC curves were generated using mRNA expression data from the [108]GSE142025 dataset. The results were validated in the [109]GSE30528 dataset. Statistical significance was set at P < 0.05. Identification of Potential Drug Candidates To obtain potential agents for gene-targeted therapies for advanced DN, three small molecule compounds were retrieved from the ChEMBL database ([110]https://www.ebi.ac.uk/, version 34) based on the hub genes. Information on the chemicals and their biological activities was obtained from the PubChem online website ([111]https://pubchem.ncbi.nlm.nih.gov). The molecular formulae and two-dimensional structures of potential drugs were downloaded from the PubChem website to aid drug discovery. Validation of the Hub Genes in Animal Models for DN Animals and Reagents C57BL/6 mice (SPF, male, 6–8 weeks old, 20 ± 2g) were purchased from the Finoco Biotechnology Co. (Shanghai, China). Mice were provided standard feed and water in a SPF grade animal facility. All animal experimental procedures were approved by the Animal Ethics committee of Finoco Biotechnology Co. (No. AUP-20240731-01) and were conducted in accordance with the National Institutes of Health Guide for the Care and Use of Laboratory Animals. Streptozocin (STZ) was purchased from Sigma Chemical Co. (Louis, USA). FastPure® Cell/Tissue Total RNA Isolation Kit V2, Reverse transcription kit HiScript III RT SuperMix for qPCR (+gDNA wiper, R323-01) and ChamQ Universal SYBR qPCR Master Mix (Q711-03) were purchased from Vazyme Biotech Co. (Nanjing, China). Animal Models The mice were fed a high-fat diet for 4 weeks and then subjected to unilateral nephrectomy. One week after surgery, the mice were followed by intraperitoneal injection of STZ (50 mg/kg dissolved in 0.05 mol/L sterile sodium citrate, pH 4.5) for three consecutive days to induce diabetes. One week after the last STZ injection, animals with plasma glucose levels above 16.7 mmol/mL were considered successful model animals (ie C57BL/6 mouse model for DN) and were recruited into the DN group and were given a high-fat diet for 12 weeks.[112]^20 The control mice were fed with normal diet and were injected intraperitoneally with the same amount of sterile sodium citrate at the same pH as the DN model mice, and were given a normal diet for 12 weeks. After a duration of 12 weeks, euthanasia was performed and kidney tissue was collected from all mice involved in the experiment. Kidneys from 16 mice were harvested in this study, including ten mice in the DN group and six mice in the control group. Randomization was performed using a computer-generated schedule by a third-party. RNA Extraction and RT-qPCR Total RNA was extracted from kidney tissue samples using the FastPure® Cell/Tissue Total RNA Isolation Kit V2 and quantified with a NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific Inc). Subsequently, first-strand complementary DNA (cDNA) was synthesized from the total RNA employing the HiScript III RT SuperMix for qPCR (+gDNA wiper) reverse transcription kit. The cDNA was amplified using the ABI 7500 Fast Real-Time PCR system (Applied Biosystems, Foster City, CA, USA) and ChamQ Universal SYBR qPCR Master Mix. PCR was conducted with a temperature profile consisting of 40 cycles at 95°C for 30s, 95°C for 10s, and 60°C for 30s. The relative expression levels of the target mRNAs were normalized to the expression levels of glyceraldehyde-3-phosphate dehydrogenase (GAPDH) and analyzed utilizing the 2^−ΔΔCt relative quantitative method. The PCR primer sequences are listed in [113]Table 1. Table 1. The PCR Primer Sequences ID Primer Name Primer Sequence (5’–3’) 14433 GAPDH GACATGCCGCCTGGAGAAAC AGCCCAGGATGCCCTTTAGT 16819 LCN2 GCAGGTGGTACGTTGTGGG CTCTTGTAGCTCATAGATGGTGC 11450 ADIPOQ GTTCCCAATGTACCCATTCGC TGTTGCAGTAGAACTTGCCAG 13058 CYBB CCCTTTGGTACAGCCAGTGAAGAT CAATCCCGGCTCCCACTAACATCA [114]Open in a new tab Statistical Analysis Data from external datasets and animal experiments were statistically analyzed and visualized using GraphPad prism 9. The two sets of data from external datasets and animal experiments were initially tested for normality using the Shapiro–Wilk normality test. If the data did not follow a normal distribution, the Mann–Whitney test was performed. Conversely, if the data exhibited a normal distribution, the homogeneity test of variance was conducted. In cases where the variance was considered to be non-uniform, Welch’s t-test was employed, while a group t-test was used when the variance was uniform. Spearman correlation analysis was used for correlation analysis. All statistical analyses were two-tailed with P < 0.05 were regarded statistically significant. Results Overview of the Study Design The overall design scheme of our study is outlined in [115]Figure 1. Firstly, We retrieved [116]GSE142025 and [117]GSE30528 from the GEO database based on the corresponding keywords and normalized the data from the two datasets. Next, we utilized the “limma” R packages to analyze DEGs. We then intersected the DEGs with PRGs and FRGs, resulting in 8 PFRDEGs. We performed GSEA functional enrichment analysis for all the genes in [118]GSE142025 and further analyzed the 8 PFRDEGs using GO/KEGG, PPI study, and Cytoscape software, which led to the identification of hub PFRDEGs. Following we assessed the clinical significance of hub PFRDEGs, evaluated their potential role in diagnosing advanced DN, and validated them using external data and DN model mice. Finally, potential drug prediction for the treatment of advanced DN was performed based on the three hub PFRDEGs. Figure 1. [119]Figure 1 [120]Open in a new tab Flowchart of research design and analyzing process of this study. Abbreviations: DEGs, differentially expressed genes; PRGs, pyroptosis-related genes; FRGs, ferroptosis-related genes; PFRDEGs, Differentially Expressed Pyroptosis-and Ferroptosis- Related Genes; GEO, Gene Expression Omnibus; GO/KEGG, Gene Ontology Enrichment Analysis / Kyoto Encyclopedia of Genes and Genomes Pathway Analysis; PPI, protein–protein interaction; ROC, receiver operating characteristic. Identification of DEGs and the Landscape of the PFRDEGs To identify the pathways and biological mechanisms related to advanced DN, we screened 1855 DEGs from [121]GSE142025 microarray dataset based on the above selection criteria. Among these, 1037 genes were up-regulated and 818 were down-regulated ([122]Figure 2A and [123]B, [124]Supplementary Table S5). Additionally, we screened 490 PRGs and 502 FRGs. By intersecting the DEGs, PRGs and FRGs, we identified 8 PFRDEGs that exhibited significant differences between advanced DN patients and healthy controls ([125]Figure 2C). These PFRDEGs include JUN, CYBB, LCN2, DPEP1, ADIPOQ, ALOX15, PTGS2, and CGAS ([126]Figure 2D). Figure 2. [127]Figure 2 [128]Open in a new tab Screening and identifying DEGs and PFRDEGs. (A) Heatmap of DEGs in [129]GSE142025. The horizontal coordinates of the graph represent individual sample designations and the vertical coordinates represent gene names. (B) Volcano plot of DEGs in [130]GSE142025. (C) Venn diagrams among DEGs, PRGs and FRGs. (D) Expression of the 8 PFRDEGs in the control and advanced DN groups. Blue and red represent advanced-DN group and control group, respectively. Functional Enrichment Analysis To delve into the function of PFRDEGs, GO/KEGG enrichment analyses were performed ([131]Supplementary Table S6 present the results in detail). Firstly, biological features of PFRDEGs were explored through GO annotation: the main biological process (BP) enrichments of PFRDEGs included cellular response to metal ion (GO:0071248), regulation of smooth muscle cell proliferation (GO:0048660), positive regulation of cytokine production (GO:0001819), and so on. In terms of molecular function (MF), the major enrichments in PFRDEGs included R-SMAD binding (GO:0070412), cAMP response element binding (GO:0035497) and glycolipid binding (GO:0051861). Secondly, KEGG pathway analysis was conducted and the results showed that PFRDEGs were remarkably enriched in the IL-17 signaling pathway (hsa04657), AGE-RAGE signaling pathway in diabetic complications (hsa04933), and TNF signaling pathway (hsa04668) ([132]Table 2, [133]Figure 3A and [134]B). Finally, to further investigate the potential pathways in DN, GSEA analyses were performed. The GSEA-enriched pathways mainly participated in oxidative phosphorylation, fcgr activation, fceri mediated NF-κb activation, and so on. Based on the NES value ranking, we selected the top 5 and bottom 5 pathway IDs, and the graphical representation of the results can be viewed in [135]Figure 3C and [136]D, [137]Table 3. [138]Supplementary Table S7 presents the results of GSEA in detail. Table 2. GO/KEGG Enrichment Analysis Results of PFRDEGs Ontology ID Description Gene Ratio P. Adjust q Value BP GO:0071248 Cellular response to metal ion 5/8 0.000 0.000 BP GO:0031953 Negative regulation of protein autophosphorylation 2/8 0.001 0.000 BP GO:0001819 Positive regulation of cytokine production 4/8 0.002 0.001 BP GO:0048660 Regulation of smooth muscle cell proliferation 3/8 0.002 0.001 BP GO:1903555 Regulation of tumor necrosis factor superfamily cytokine production 2/8 0.025 0.010 MF GO:0016175 Superoxide-generating NAD(P)H oxidase activity 1/8 0.023 0.013 MF GO:0035497 cAMP response element binding 1/8 0.030 0.016 MF GO:0043027 Cysteine-type endopeptidase inhibitor activity involved in apoptotic process 1/8 0.035 0.019 MF GO:0070412 R-SMAD binding 1/8 0.036 0.020 MF GO:0051861 Glycolipid binding 1/8 0.042 0.023 KEGG Hsa04657 IL-17 signaling pathway 3/7 0.002 0.002 KEGG Hsa04216 Ferroptosis 2/7 0.015 0.011 KEGG Hsa04933 AGE-RAGE signaling pathway in diabetic complications 2/7 0.038 0.028 KEGG Hsa04625 C-type lectin receptor signaling pathway 2/7 0.038 0.028 KEGG Hsa04668 TNF signaling pathway 2/7 0.038 0.028 [139]Open in a new tab Figure 3. [140]Figure 3 [141]Open in a new tab GO/KEGG enrichment results of PFRDEGs and GSEA results of all the genes in [142]GSE142025. (A) Bubble plot of the results of GO/KEGG enrichment analysis. (B) Histogram of the results of GO/KEGG enrichment analysis. (C) The graphical representation of the last 5 pathways in the GSEA results (The NES values are negative.). (D) The graphical representation of the The first 5 pathways in the GSEA results (The NES values are positive.). Table 3. GSEA Results ID NES P. Adjust q Value Reactome_the_citric_acid_tca_cycle_and_respiratory_electron_transport −4.220 0.000 0.000 Kegg_oxidative_phosphorylation −4.143 0.000 0.000 Wp_electron_transport_chain_oxphos_system_in_mitochondria −4.119 0.000 0.000 Kegg_parkinsons_disease −4.080 0.000 0.000 Reactome_respiratory_electron_transport_atp_synthesis_by_chemiosmotic_c oupling_and_heat_production_by_uncoupling_proteins −4.059 0.000 0.000 Reactome_role_of_lat2_ntal_lab_on_calcium_mobilization 3.749 0.000 0.000 Reactome_cd22_mediated_bcr_regulation 3.739 0.000 0.000 Reactome_role_of_phospholipids_in_phagocytosis 3.692 0.000 0.000 Reactome_fcgr_activation 3.680 0.000 0.000 Reactome_fceri_mediated_nf_kb_activation 3.676 0.000 0.000 [143]Open in a new tab PPI Network Analysis for Identifying Hub Genes and Interactions Among Hub Genes To obtain hub DEGs for a more in-depth analysis of the role of pyroptosis and ferroptosis in DN, PPI network of the PFRDEGs was constructed using the STRING online database. Out of the 8 PFRDEGs, 6 PFRDEGs were found to engage in dialogue with each other, while the remaining 2 PFRDEGs did not show any interaction ([144]Figure 4A). The hub genes were identified using the MCODE-plugin in Cytoscape (Detailed parameters: Degree cutoff: 2; Cluster Finding: Haircut; Node Score Cutoff: 0.2; K-Core: 2; Max Depth: 100), including JUN, LCN2, ADIPOQ and CYBB ([145]Figure 4B). To gain further insights into the intermolecular interaction of these hub genes, the study utilized GeneMANIA online software to search for related genes, resulting in the identification of 20 genes with a total of 425 links ([146]Figure 4C). Figure 4. [147]Figure 4 [148]Open in a new tab PPI network among proteins encoded by PFRDEGs and prediction of intermolecular interactions of the hub genes. (A) PPI network analysis of proteins encoded by PFRDEGs. (B) The four hub genes were identified using the MCODE-plugin of Cytoscape software. (C) Prediction of intermolecular interactions of the hub genes. (The purple line means that the two proteins are co-expressed, the yellow line means that the two proteins are shared protein domains, the pink line means that the two proteins are physical interactions, the blue line means that the two proteins are co-localisation, the Orange line means that the two proteins are predicted.). Immune Infiltration Analysis and Correlations Between Hub Genes and Immune Cells The CIBERSORTx algorithm has been applied to investigate the landscape of immune infiltration in advanced DN. The boxplots of 22 immune cell types between advanced DN group and Control group are shown in [149]Figure 5A, [150]Supplementary Table S8. Specifically, the types of immune cells in advanced DN were B cells naive, T cells CD8, T cells CD4 naive, T cells CD4 memory resting, NK cells activated, Macrophages M0, Macrophages M2 and Mast cells resting. In addition, we analyzed the correlation between each immune cell and the four hub genes. As shown in [151]Figure 5B, [152]Supplementary Table S9, JUN expression was associated with B cells memory, T cells CD4 naive, Macrophages M2, Mast cells resting; CYBB expression was associated with B cells naive, B cells memory, T cells CD8, NK cells resting, NK cells activated, Macrophages M0, Macrophages M2, and Mast cells resting; LCN2 expression was associated with B cells naive, T cells CD8, T cells CD4 memory resting, NK cells activated, Macrophages M0, Macrophages M2 and Mast cells resting. In summary, these results suggest that JUN, LCN2 and CYBB may contribute to the immune microenvironment of advanced DN. Figure 5. [153]Figure 5 [154]Open in a new tab Immune infiltration analysis. (A) The boxplots of 22 immune cells between advanced DN samples and control samples. (B) Correlation between each of the immune cells and four hub genes. Validation of Hub Genes in External Datasets To visually compare the differences in hub genes expression between the advanced DN group and the control group, we used violin plots. As shown in [155]Figure 6A, CYBB and LCN2 in [156]GSE142025 were significantly upregulated, JUN was significantly downregulated in the advanced DN group compared with the control group. To validate the reliability of the four hub genes, we verified the expression of these genes in the [157]GSE30528 dataset. As shown in [158]Figure 6B, CYBB and LCN2 in this dataset were significantly upregulated in DN group compared with the control group. In summary, only CYBB and LCN2 were discrepant in both external datasets. Figure 6. [159]Figure 6 [160]Open in a new tab The expression of hub genes in [161]GSE142025 and validation in the external datasets. (A) The expression of CYBB, LCN2, ADIPOQ and JUN in [162]GSE142025. (B) The expression of CYBB, LCN2, ADIPOQ and JUN in [163]GSE30528 (validation dataset). Hub Genes (CYBB, LCN2 and ADIPOQ) Have Potential Diagnostic Power To evaluate the diagnostic efficacy of the four hub genes identified in the previous analysis, ROC curves were generated and the AUC value was used to evaluate their effectiveness in diagnosing advanced DN. The diagnostic efficacy of three hub genes (CYBB, LCN2 and ADIPOQ) was found to be superior. [164]Figure 7A–C provides an overview of the ROC curves for the three hub genes, showing that the AUC for CYBB was 0.952 (95% CI, 0.859–1.000; P = 0.0003); AUC for LCN2 was 0.968 (95% CI, 0.910–1.000; P = 0.009), and AUC for ADIPOQ was 0.926 (95% CI, 0.836–1.000; P = 0.0034). This diagnostic efficacy of the three genes was further confirmed in the [165]GSE30528 dataset and the AUC was 0.709 (95% CI, 0.481–0.937), 0.821 (95% CI, 0.635–1.000) and 0.744 (95% CI, 0.529–0.958) for CYBB, LCN2 and ADIPOQ, respectively ([166]Figure 7D–F). These findings demonstrate that these hub genes may play a significant role in the pathological mechanisms underlying DN. Figure 7. [167]Figure 7 [168]Open in a new tab Diagnostic value of the hub genes with ROC curves and validation in external dataset. (A-C) ROC curve analysis of the hub genes in the [169]GSE142025. (D-F) Validation of the diagnostic value of the hub genes in the [170]GSE30528. Expression of CYBB and LCN2 Was Significantly Higher in DN Group and Negatively Correlated with eGFR The aforementioned findings indicate that CYBB, LCN2 and ADIPOQ have potential diagnostic value. To further validate the reliability of the data obtained from this bioinformatics analysis, we verified the expression of CYBB and LCN2 by cross-referencing it with data from Nephroseq V5 database. The results obtained by querying Nephroseq V5 database showed that in DN patients, the expression of CYBB and LCN2 was significantly higher than that of healthy living donors [P = 0.0022 ([171]Figure 8A);[172]^40 P < 0.0001 ([173]Figure 8B)[174]^40]; Moreover, the expression of Cybb and Lcn2 was also significantly up-regulated in DN mice [P = 0.0142 ([175]Figure 8C);[176]^41 P < 0.0001 ([177]Figure 8D)[178]^41]. To further explore the clinical significance of the three genes, we performed a correlation analysis between the three genes and eGFR, because eGFR is a crucial impact factor on the prognosis of DN. The results showed that the expression of CYBB was negatively correlated with eGFR [(r = −0.7061, P = 0.0002),[179]^42 [180]Figure 8E], the expression of LCN2 was negatively correlated with eGFR [(r = −0.7588, P = 0.0068),[181]^42 [182]Figure 8F]. These findings suggest that the expression of LCN2 and CYBB can be utilized to a certain extent for predicting renal function in DN. Figure 8. [183]Figure 8 [184]Open in a new tab Expression of CYBB and LCN2 was significantly higher in DN group and negatively correlated with eGFR in external datasets. (A and B) Expression of CYBB and LCN2 was significantly higher in DN patients. (C and D) Expression of Cybb and Lcn2 was significantly higher in DN mice. (E and F) The expression of CYBB and LCN2 was negatively correlated with eGFR in DN patients. The Expression of the Cybb and Lcn2 is Markedly Elevated in C57BL/6 Mice Models for DN Efficacy are key requirement for the translation of any new biomarkers to the clinic. In order to further validate the data obtained in the bioinformatics analysis, a series of biological experiments were conducted, including animal experiments and molecular biology experiments. mRNA expression of the hub genes in the kidney of C57BL/6 mice were measured using RT-qPCR. As depicted in [185]Figure 9, when compared to the control group, the mRNA expression of Cybb and Lcn2 were significantly higher expression in the DN group (Cybb: 1.310±0.417, p=0.0022; Lcn2: 1.921±0.592, p=0.0078). However, there were no statistically differences in the expression of Jun and Adipoq between the two groups. Figure 9. [186]Figure 9 [187]Open in a new tab Validation of the hub genes in C57BL/6 mice models for DN. The mRNA expression of the four hub genes in C57BL/6 mice models for DN were determined by RT-qPCR assay. (control group: n =6, DN group: n=10.). Prediction of Marker Gene-Targeted Therapeutic Agents The another primary aim of this study is to forecast potential therapeutic agents that could be effective in treating advanced DN. To achieve this, the ChEMBL database was thoroughly investigated, and three therapeutic agents—CHEMBL420360, CHEMBL182828, and CHEMBL1823903—were identified as targeting the aforementioned hub genes ([188]Table 4). This discovery offers a robust foundation of both experimental and theoretical evidence to support the development of precision treatments for advanced DN. Table 4. Potential Therapeutic Agents List Compound ID Molecular Formula 2D Structure CHEMBL420360 C19H27N5O4 graphic file with name JIR-17-5763-i0001.jpg CHEMBL182828 C28H25N3O4S graphic file with name JIR-17-5763-i0002.jpg CHEMBL1823903 C24H21Cl2NO3S graphic file with name JIR-17-5763-i0003.jpg [189]Open in a new tab Discussion As a prevalent and recurring microvascular complication of diabetes mellitus, DN has emerged as a significant contributor to ESKD and a public health challenge worldwide.[190]^43 Studies have shown that ferroptosis and pyroptosis can cause renal damage through corresponding pathways, and inhibition of ferroptosis and pyroptosis, respectively, can prevent or delay renal damage and renal fibrosis. Multiple studies have demonstrated the existence of numerous interconnections between pyroptosis and ferroptosis, encompassing shared triggers, molecular components, and protective mechanisms.[191]^23 More importantly, many recent studies have found that there is a mechanistic crosstalk between pyroptosis and ferroptosis, which provides new perspectives on the mechanistic study of the disease and new avenues for the treatment of the disease.[192]^28–31 Consequently, in this study, we conducted the combined study of ferroptosis and pyroptosis in advanced DN by means of bioinformatics. Initially, we identified eight PFRDEGs. Subsequently, We conducted GO/KEGG enrichment analysis on the 8 PFRDEGs and GSEA on all genes in [193]GSE142025. The GO analysis revealed that the PFRDEGs were primarily involved in positive regulation of cytokine production, regulation of smooth muscle cell proliferation, cellular response to metal ion, and so on. The KEGG enrichment analysis showed that the PFRDEGs were enriched in the IL-17 signaling pathway, AGE-RAGE signaling pathway in diabetic complications, and TNF signaling pathway. Additionally, we conducted GSEA, which mainly participated in oxidative phosphorylation, fcgr activation, fceri mediated NF-κb activation, and so on. Subsequently, four hub PFRDEGs (JUN, CYBB, LCN2, and ADIPOQ) were selected, and the three hub PFRDEGs―CYBB, LCN2, and ADIPOQ were up-regulated while JUN was down-regulated. To further explore the clinical significance of hub PFRDEGs, we assessed the diagnostic value of the hub PFRDEGs and did a correlation analysis between these hub PFRDEGs and clinical indicators. We found that CYBB, LCN2, and ADIPOQ exhibited high accuracy (all AUC>0.9) in predicting DN and expression of CYBB and LCN2 was negatively correlated with eGFR. Additionally, we also mined three therapeutic agents targeting CYBB, LCN2, and ADIPOQ. These findings suggest that CYBB, LCN2, and ADIPOQ play a significant role in diagnosing advanced DN, and provide a theoretical basis for the precise treatment for advanced DN. In our study, we found that CYBB could be considered as a potential biomarker for advanced DN. CYBB, also known as Cytochrome b-245 beta chain, encodes the gp91phox protein. Mutations in this gene are associated with X-linked chronic granulomatous disease.[194]^44 Despite the fact that reports on this gene’s involvement in DN are limited, our study found that Cybb expression was significantly overexpressed in renal tissues of DN mice (p=0.0022), which aligns with the only literature results.[195]^45 The role of CYBB in DN remains unclear due to the limited number of studies, but our discovery of the CYBB can provide valuable insights for the diagnostic markers of DN in future research. In addition, ADIPOQ (adiponectin), a protein secreted by adipose tissue, is deemed to be an endogenous insulin sensitizer. ADIPOQ exhibits anti-inflammatory and antiatherogenic effects, and regulates glucose and lipid metabolism as well as insulin action.[196]^46 Abnormal levels of serum adiponectin have been correlated with T2DM, insulin resistance, obesity, cardiovascular diseases, and nephropathy.[197]^47 In our in vivo experiment, we did not observe a significant difference in ADIPOQ expression compared to the control group. The exact mechanism of ADIPOQ’s role in DN has not been extensively studied and requires further research. Furthermore, LCN2, also known as neutrophil gelatinase-associated lipocalin (NGAL), is a transporter protein expressed in a wide range of cells. During acute kidney injury, LCN2 is secreted in high amounts into the urine and blood by tubular cells even before serum creatinine levels increase.[198]^48 In individuals with type 1 diabetes mellitus, LCN2 concentration is elevated in both blood and urine samples prior to the presence of microalbuminuria.[199]^49 Studies have shown a positive correlation between LCN2 and albumin excretion rate as well as hemoglobin A1c levels.[200]^50^,[201]^51 In patients with T2DM, high expression of LCN2 serves as a strong predictor of GFR impairment. We further found the meaningful value of LCN2 in the diagnosis of advanced DN by a joint of bioinformatics analysis and biological experiments based on previous literature. Therefore, LCN2 is expected to serve as a potential biomarker for advanced DN. However, this study has several limitations. The evidence is based on publicly available data, and although we have performed expression validation with another dataset and validated in vivo in mice, further experiments are needed to validate the three diagnostic markers before they can be applied to the clinic. In summary, we applied bioinformatics methods to jointly conduct in-depth data mining of the dataset of advanced DN from two perspectives, namely ferroptosis and pyroptosis, and found that CYBB, LCN2 and ADIPOQ have potential diagnostic and targeted therapeutic value. This study is expected to provide constructive reference and basis for clinical detection and targeted therapy of advanced DN. Conclusions Our study identified CYBB, LCN2, and ADIPOQ as potential biomarkers for advanced DN by utilizing bioinformatics analysis of disease databases. Furthermore, our research validated the diagnostic value of Cybb and Lcn2 through animal experiments. Additionally, we discovered three small molecule drugs via drug screening that show promise for targeted treatment of advanced DN. These findings not only provide new insights into the role of pyroptosis and ferroptosis in advanced DN, but also lay the foundation for future research into non-invasive diagnostic markers and precision therapies for advanced DN. Funding Statement This work is supported by the Senior Integrated Chinese and Western Medicine Talent Training Program in Shanghai, China [ZY (2018–2020)-RCPY-2020]. Author Contributions All authors made a significant contribution to the work reported, whether that is in the conception, study design, execution, acquisition of data, analysis and interpretation, or in all these areas; took part in drafting, revising or critically reviewing the article; gave final approval of the version to be published; have agreed on the journal to which the article has been submitted; and agree to be accountable for all aspects of the work. Disclosure The authors declare no competing interests in this work. References