Identification of VCAN as Hub Gene for Diabetic Kidney Disease Immune Injury Using Integrated Bioinformatics Analysis

Background: Diabetic kidney disease (DKD) is a leading cause of chronic kidney disease in China. Tubular injury contributes to the progression of DKD. Our study was conducted to explore the differential gene expression profiles between kidneys from patients with DKD and kidney living donors (LDs). Methods: In total, seven DKD and eighteen LD gene expression profiles from the GSE104954 dataset were downloaded from the Gene Expression Omnibus database. Differentially expressed genes (DEGs) were analyzed in R with the limma package. DEGs were uploaded to the g:Profiler online database to explore the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Ingenuity pathway analysis (IPA) was carried out using online IPA software. Weighted gene co-expression network analysis (WGCNA) was performed using the WGCNA R package. By integrating DEGs and genes from the top 1 phenotype-gene associated module, we determined the hub gene. We next tested the hub gene, VCAN, in the GSE30122 dataset. We also validated the versican levels in human kidney tissues, explored immune cell type enrichment using an online database xCell, and investigated the correlation between cell types and VCAN expression. Results: A total of 563 DEGs was identified. A large number of pathways were involved in the immune response process according to the results of GO, KEGG, and IPA. Using WGCNA, we selected the lightcyan module in which genes showed the strongest correlation with the phenotype and smallest P-value. We also identified VCAN as a hub gene by integrating DEG analysis and WGCNA. Versican expression was upregulated in human diabetic kidney tissue. Moreover, versican was speculated to play a role in immune injury according to the enrichment of functions and signaling pathways. VCAN transcript levels correlate with the assembly of immune cells in the kidney. Conclusion: Immune processes played an essential role in DKD tubulointerstitium injury. The hub gene VCAN contributed to this process.