We aim to identify the key protein interaction networks and implicated pathways of BK virus nephropathy (BKVN) via bioinformatic methods. The microarray data GSE75693 of 30 patients with stable kidney transplantation and 15 with BKVN were downloaded and analyzed by using the limma package to identify differentially expressed genes (DEGs). Then the gene ontology (GO) functional enrichment analysis and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were done to investigate the molecular function (MF), biological process (BP), cellular components (CC) and pathways of DEGs. Finally, protein-protein interactions (PPIs) were constructed, and the hub proteins were identified. As a result, 249 up-regulated genes and 253 down-regulated genes of BKVN patients were selected based on criteria of P > 0.01 and fold change > 2.0. GO and KEGG showed that DEGs were mainly located in nucleus and cytosol, and were implicated in the immune responses. In the PPI analysis, 26 up-regulated and 8 down-regulated proteins composed the pivotal interaction network. CXCL10, EGF and STAT1 were identified as hub proteins in BKVN. In conclusion, CXCL10, EGF and STAT1 may induce kidney injuries by promoting inflammation and prohibiting reparation of tissue damage in BKVN.