RESUMO
Little is known about the urinary virome and how it interacts with the host, particularly in renal transplant diseases. Using metagenomic sequencing, we characterized the urinary virome of 23 kidney transplant recipients longitudinally (11 BKV+ patients and 12 BKV- patients). We applied linear mixed effects models, PERMANOVA, k-means clustering, and MaAsLin2 algorithms to determine virome signatures associated with post-transplant time, BK viremia status, and patient sex. We found that the richness and alpha diversity of urinary virome were significantly different in renal transplant recipients with BKV+ over time in comparison to BKV- (richness P = 0.012, alpha P < 0.0001). Female BKV- patients had significantly higher virome richness than males (P = 0.0063). Virome beta diversity was significantly different between patients by BKV status (P < 0.001). Additionally, we identified underlying interactions between patient sex and BKV status, in terms of virome beta diversity (P = 0.008). BK polyomavirus infections were primarily of subtypes IA, IB1, and IB2. The non-BK dominant samples clustered into six urinary virome community states. BKV- samples had more anelloviruses than BKV+ samples though this difference was not statistically significant. Lastly, we identified specific viruses, associated with BKV+ and time in our samples. Our results indicate that dynamic alterations in the urinary virome over the post-transplant period in kidney transplant recipients can be shaped by BK viremia and patient sex. These findings advance our fundamental understanding of the urinary virome and support a new line of investigation in renal disease and transplantation. IMPORTANCE: The urinary microbiome is increasingly implicated in renal health and disease. While most research focuses on bacteria communities of the microbiome, factors that influence the urinary virome are not understood. Here, we investigated the urinary virome of 23 adult kidney transplant recipients longitudinally over 14 weeks post-transplant. We show that alterations in the urinary virome are associated with kidney transplant recipients with BK polyomavirus viremia that can lead to BK nephropathy and allograft rejection. By modeling the temporal dynamics post-transplant, we delineated specific profiles of the urinary virome associated with patient sex and urinary community states. These findings reveal fundamental aspects of the urinary virome that can be leveraged to better manage kidney diseases.