Identification of Hub Genes and Potential ceRNA Networks of Diabetic Nephropathy by Weighted Gene Co-Expression Network Analysis

Abstract

Diabetic nephropathy (DN) is one of the most common microvascular complications in diabetic patients, and is the main cause of end-stage renal disease. The exact molecular mechanism of DN is not fully understood. The aim of this study was to identify novel biomarkers and mechanisms for DN disease progression by weighted gene co-expression network analysis (WGCNA). From the GSE142153 dataset based on the peripheral blood monouclear cells (PBMC) of DN, we identified 234 genes through WGCNA and differential expression analysis. Gene Ontology (GO) annotations mainly included inflammatory response, leukocyte cell-cell adhesion, and positive regulation of proteolysis. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways mostly included IL-17 signaling pathway, MAPK signaling pathway, and PPAR signaling pathway in DN. A total of four hub genes (IL6, CXCL8, MMP9 and ATF3) were identified by cytoscape, and the relative expression levels of hub genes were also confirmed by RT-qPCR. ROC curve analysis determined that the expression of the four genes could distinguish DN from controls (the area under the curve is all greater than 0.8), and Pearson correlation coefficient analysis suggested that the expression of the four genes was related to estimated glomerular filtration rate (eGFR) of DN. Finally, through database prediction and literature screening, we constructed lncRNA-miRNA-mRNA network. We propose that NEAT1/XIST/KCNQ1T1-let-7b-5p-IL6, NEAT1/XIST-miR-93-5p-CXCL8 and NEAT1/XIST/KCNQ1T1-miR-27a-3p/miR-16-5p-ATF3 might be potential RNA regulatory pathways to regulate the disease progression of early DN. In conclusion, we identified four hub genes, namely, IL6, CXCL8, MMP9, and ATF3, as markers for early diagnosis of DN, and provided insight into the mechanisms of disease development in DN at the transcriptome level.