Machine learning-based identification of novel hub genes associated with oxidative stress in lupus nephritis: implications for diagnosis and therapeutic targets

Abstract

Background Lupus nephritis (LN) is a complication of SLE characterised by immune dysfunction and oxidative stress (OS). Limited options exist for LN. We aimed to identify LN-related OS, highlighting the need for non-invasive diagnostic and therapeutic approaches. Methods LN-differentially expressed genes (DEGs) were extracted from Gene Expression Omnibus datasets (GSE32591, GSE112943 and GSE104948) and Molecular Signatures Database for OS-associated DEGs (OSEGs). Functional enrichment analysis was performed for OSEGs related to LN. Weighted gene co-expression network analysis identified hub genes related to OS-LN. These hub OSEGs were refined as biomarker candidates via least absolute shrinkage and selection operator. The predictive value was validated using receiver operating characteristic (ROC) curves and nomogram for LN prognosis. We evaluated LN immune cell infiltration using single-sample gene set enrichment analysis and CIBERSORT. Additionally, gene set enrichment analysis explored the functional enrichment of hub OSEGs in LN. Results The study identified four hub genes, namely STAT1 , PRODH , TXN2 and SETX , associated with OS related to LN. These genes were validated for their diagnostic potential, and their involvement in LN pathogenesis was elucidated through ROC and nomogram. Additionally, alterations in immune cell composition in LN correlated with hub OSEG expression were observed. Immunohistochemical analysis reveals that the hub gene is most correlated with activated B cells and CD8 T cells. Finally, we uncovered that the enriched pathways of OSEGs were mainly involved in the PI3K-Akt pathway and the Janus kinase-signal transducer and activator of transcription pathway. Conclusion These findings contribute to advancing our understanding of the complex interplay between OS, immune dysregulation and molecular pathways in LN, laying a foundation for the identification of potential diagnostic biomarkers and therapeutic targets. Data are available upon reasonable request. The LN cohort’s transcriptional dataset (GSE32591, GSE112943 and GSE104948) was downloaded from the Gene Expression Omnibus (GEO) database ().