Characterizing hub biomarkers for post-transplant renal fibrosis and unveiling their immunological functions through RNA sequencing and advanced machine learning techniques

Abstract Background Kidney transplantation stands out as the most effective renal replacement therapy for patients grappling with end-stage renal disease. However, post-transplant renal fibrosis is a prevalent and irreversible consequence, imposing a substantial clinical burden. Unfortunately, the clinical landscape remains devoid of reliable biological markers for diagnosing post-transplant renal interstitial fibrosis. Methods We obtained transcriptome and single-cell sequencing datasets of patients with renal fibrosis from NCBI Gene Expression Omnibus (GEO). Subsequently, we employed Weighted Gene Co-Expression Network Analysis (WGCNA) to identify potential genes by integrating core modules and differential genes. Functional enrichment analysis was conducted to unveil the involvement of potential pathways. To identify key biomarkers for renal fibrosis, we utilized logistic analysis, a LASSO-based tenfold cross-validation approach, and gene topological analysis within Cytoscape. Furthermore, histological staining, Western blotting (WB), and quantitative PCR (qPCR) experiments were performed in a murine model of renal fibrosis to verify the identified hub genes. Moreover, molecular docking and molecular dynamics simulations were conducted to explore possible effective drugs. Results Through WGCNA, the intersection of core modules and differential genes yielded a compendium of 92 potential genes. Logistic analysis, LASSO-based tenfold cross-validation, and gene topological analysis within Cytoscape identified four core genes (CD3G, CORO1A, FCGR2A, and GZMH) associated with renal fibrosis. The expression of these core genes was confirmed through single-cell data analysis and validated using various machine learning methods. Wet experiments also verified the upregulation of these core genes in the murine model of renal fibrosis. A positive correlation was observed between the core genes and immune cells, suggesting their potential role in bolstering immune system activity. Moreover, four potentially effective small molecules (ZINC000003830276-Tessalon, ZINC000003944422-Norvir, ZINC000008214629-Nonoxynol-9, and ZINC000085537014-Cobicistat) were identified through molecular docking and molecular dynamics simulations. Conclusion Four potential hub biomarkers most associated with post-transplant renal fibrosis, as well as four potentially effective small molecules, were identified, providing valuable insights for studying the molecular mechanisms underlying post-transplant renal fibrosis and exploring new targets.