Background: Tubulointerstitial injury (TIL) is common in chronic kidney disease (CKD), which in turn, leads to loss of renal function. The aims of present study were to screen critical genes with tubulointerstitial lesion in CKD by weighted gene correlation network analysis (WGCNA). Methods: GSE104954 gene expression data downloaded from GEO database were used for analysis for differential expression genes(DEGs); Meanwhile, 90 expression data of tubulointerstitial samples in GSE47185 combined with clinic information were applied to WGCNA. According to the enrichment analysis and relationship with estimated glomerular filtration rate (eGFR), the eGFR-associated modules were defined, and the hub gene is selected according to the average intramodular connectivity (Kwithin); Clinical data from Nephroseq were obtained to further validate the relationship between CKD and hub genes. Results: Totally 294 DEGs were screened. Using the WGCNA, we identified 15 co-expressed modules, the blue, brown and yellow modules exhibited strongly association with eGFR, and were significantly enriched in several signaling pathways that have been reported involved in pathogenesis of CKD. Furthermore, it was found that the 4 genes (PLG, ITGB2, CTSS and CCL5) was one of the DEGs which also be identified as hub genes according to Kwithin. Finally, the Nephroseq online tool showed that the tubulointerstitial expression levels of PLG significantly positively correlated with the eGFR, while ITGB2, CTSS and CCL5 connected negatively to the eGFR. Conclusion: WGCNA is an efficient approach to system biology. By this procedure, the present study improved the understanding of the transcriptome status of renal tubulointerstitial injury in CKD.