Integrated Analysis of Multiple Microarray Studies to Identify Core Gene-Expression Signatures Involved in Tubulointerstitial Injury in Diabetic Nephropathy

Diabetic nephropathy is a leading cause of end-stage renal disease in both developed and developing countries. It is lack of specific diagnosis, and the pathogenesis remains unclarified in diabetic nephropathy, following the unsatisfactory effects of existing treatments. Therefore, it is very meaningful to find biomarkers with high specificity and potential targets. Two datasets, GSE30529 and GSE47184 from GEO based on diabetic nephropathy tubular samples, were downloaded and merged after batch effect removal. A total of 545 different expression genes screened with log2FC > 0.5 were weighted gene coexpression correlation network analysis, and green module and blue module were identified. The results of KEGG analyses both in green module and GSEA analysis showed the same two enriched pathway, focal adhesion and viral myocarditis. Based on the intersection among WGCNA focal adhesion/Viral myocarditis, GSEA focal adhesion/viral myocarditis, and PPI network, 17 core genes, ACTN1, CAV1, PRKCB, PDGFRA, COL1A2, COL6A3, RHOA, VWF, FN1, HLA-F, HLA-DPB1, ITGB2, HLA-DRA, HLA-DMA, HLA-DPA1, HLA-B, and HLA-DMB, were identified as potential biomarkers in diabetic tubulointerstitial injury and were further validated externally for expression at GSE99325 and GSE104954 and clinical feature at nephroseq V5 online platform. CMap analysis suggested that two compounds, LY-294002 and bufexamac, may be new insights for therapeutics of diabetic tubulointerstitial injury. Conclusively, it was raised that a series of core genes may be as potential biomarkers for diagnosis and two prospective compounds.