Exploring the mechanism of BK polyomavirus-associated nephropathy through consensus gene network approach

BK polyomavirus-associated nephropathy occurs in kidney transplant recipients under immunosuppressive treatment. BK polyomavirus is implicated in cancer development and invasion, and case reports of renal cell carcinoma and urothelial carcinoma possibly associated with BK polyomavirus has been reported. Further, it has been suggested that the immune responses of KT-related diseases could play a role in the pathogenesis and progression of renal cell carcinoma. Thus, we thought to examine the relationship between BK polyomavirus-associated nephropathy and renal cell carcinoma in terms of gene expression. To identify the common and specific immune responses involved in kidney transplantation-related diseases with a specific focus on BK polyomavirus-associated nephropathy, we performed consensus weighted gene co-expression network analysis on gene profile datasets of renal biopsy samples from different institutions. After the identification of gene modules and validation of the obtained network by immunohistochemistry of the marker across kidney transplantation-related diseases, the relationship between prognosis of renal cell carcinoma and modules was assessed. We included the data from 248 patients and identified the 14 gene clusters across the datasets. We revealed that one cluster related to the translation regulating process and DNA damage response was specifically upregulated in BK polyomavirus-associated nephropathy. There was a significant association between the expression value of hub genes of the identified cluster including those related to cGAS-STING pathway and DNA damage response, and the prognosis of renal cell carcinoma. The study suggested the potential link between kidney transplantation-related diseases, especially specific transcriptomic signature of BK polyomavirus associated nephropathy and renal cell carcinoma.