Dynamic viral integration patterns actively participate in the progression of BK polyomavirus-associated diseases after renal transplantation.

The classical lytic infection theory along with large T antigen-mediated oncogenesis cannot explain the BK polyomavirus (BKPyV)-associated tumor secondary to BKPyV-associated nephropathy (BKVAN), viremia/DNAemia, and viruria after renal transplantation. This study performed virome capture sequencing and pathological examination on regularly collected urine sediment and peripheral blood samples, and BKVAN and tumor biopsy tissues of 20 patients with BKPyV-associated diseases of different stages. In the early noncancerous stages, well-amplified integration sites were visualized by in situ polymerase chain reaction, simultaneously with BKPyV inclusion bodies and capsid protein expression. The integration intensity, the proportion of microhomology-mediated end-joining integration, and host PARP-1 and POLQ gene expression levels increased with disease progression. Furthermore, multiomics analysis was performed on BKPyV-associated urothelial carcinoma tissues, identifying tandem-like structures of BKPyV integration using long-read genome sequencing. The carcinogenicity of BKPyV integration was proven to disturb host gene expression and increase viral oncoprotein expression. Fallible DNA double-strand break repair pathways were significantly activated in the parenchyma of BKPyV-associated tumors. Olaparib showed an antitumor activity dose-response effect in the tumor organoids without BRCA1/2 genes mutation. In conclusion, the dynamic viral integration patterns actively participate in the progression of BKPyV-associated diseases and thus could be a potential target for disease monitoring and intervention.