Viral integration in BK polyomavirus-associated urothelial carcinoma in renal transplant recipients: multistage carcinogenesis revealed by next-generation virome capture sequencing

BK polyomavirus (BKPyV)-associated cancer after transplantation has gained increasing attention. However, the role of BKPyV integration on oncogenesis is still unclear. In this study, next-generation virome capture sequencing of primary and metastatic tumors were performed in three patients with BKPyV-associated urothelial carcinoma after renal transplantation. As a result, a total of 332 viral integration sites were identified in the six tumors. Integration of BKPyV in both primary and metastatic tumors followed the mechanism of microhomology-mediated end joining mostly, since microhomologies between human and BKPyV genomes were significantly enriched in flanking regions of 84% of the integration sites. Viral DNA breakpoints were nonrandom and tended to assemble in large T gene, small T gene and viral protein 2 gene. There were three, one and one consensus integration sites between the primary and metastatic tumors, which affected LINC01924, eIF3c, and NEIL2 genes in the three cases respectively. Thus, we concluded that integration of BKPyV was a continuous process occurring in both primary and metastatic tumors, generating heterogenous tumor cell populations. Through this ongoing process, certain cell populations might have gained growth advantage or metastatic potential, as a result of viral integration either affecting the cellular genes where the viral DNA integrated to or altering the expression or function of the viral genes.