BKTyper: Free Online Tool for Polyoma BK Virus VP1 and NCCR Typing

Seasonality of respiratory, enteric, and urinary viruses revealed by wastewater genomic surveillance

Wastewater surveillance can reveal population-level infectious disease burden and emergent public health threats can be reliably assessed through wastewater surveillance. While molecular methods for wastewater monitoring of microorganisms have traditionally relied on PCR-based approaches, next-generation sequencing (NGS) can provide deeper insights via genomic analyses of multiple diverse pathogens. We conducted a year-long sequencing surveillance of 1,408 composite wastewater samples collected from 12 neighborhood-level access points in the greater Tempe area, Arizona, USA, and show that variation in wastewater viruses is driven by seasonal time and location. The temporal dynamics of viruses in wastewater were influenced cyclically, with the most dissimilarity between samples 23 weeks apart (i.e., winter vs summer, spring vs fall). We identified diverse urinary and enteric viruses including polyomaviruses, astroviruses, and noroviruses, and showed that their genotypes/subtypes shifted across seasons. We show that while wastewater data of certain respiratory viruses like severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) strongly correlate with clinical case rates, laboratory-reported case incidences were discordant with surges of high viral load in wastewater for other viruses like human coronavirus 229E. These results demonstrate the utility of wastewater sequencing for informing decision-making in public health.IMPORTANCEWastewater surveillance can provide insights into the spread of pathogens in communities. Advances in next-generation sequencing (NGS) methodologies allow for more precise detection of viruses in wastewater. Long-term wastewater surveillance of viruses is an important tool for public health preparedness. This system can act as a public health observatory that gives real-time early warning for infectious disease outbreaks and improved response times.

Time-dependent variations in BK polyomavirus genome from kidney transplant recipients with persistent viremia

BK polyomavirus (BKPyV) is a human DNA virus that resides latent in the host's renal tissue. Reactivation occurs occasionally and in case of kidney transplantation, it can lead to polyomavirus-associated nephropathy (PVAN). Due to the lack of specific antivirals for BKPyV and despite the risk of allograft rejection, reduction of immunosuppression remains the main approach for treating PVAN. Current data suggests that mutations can accumulate over time in the major capsid protein VP1 and can lead to neutralization escape in kidney transplant recipients. Herein, we show that mutations occur throughout the entire BKPyV genome, including in VP1. Changes were identified by per-patient comparison of viral genome sequences obtained in samples from 32 kidney recipients with persistent viremia collected at different post-transplant time-points. Amino acid changes were observed in both earlier and later post-transplant samples, although some of them were only found in later samples. Changes in VP1 mainly consisted in the introduction of a new amino acid. A switch back to the conservative amino acid was also observed. This should be considered in future approaches for treating BKPyV infection in kidney transplant recipients.

Evidence for virus-mediated oncogenesis in bladder cancers arising in solid organ transplant recipients

A small percentage of bladder cancers in the general population have been found to harbor DNA viruses. In contrast, up to 25% of tumors of solid organ transplant recipients, who are at an increased risk of developing bladder cancer and have an overall poorer outcomes, harbor BK polyomavirus (BKPyV). To better understand the biology of the tumors and the mechanisms of carcinogenesis from potential oncoviruses, we performed whole genome and transcriptome sequencing on bladder cancer specimens from 43 transplant patients. Nearly half of the tumors from this patient population contained viral sequences. The most common were from BKPyV (N=9, 21%), JC polyomavirus (N=7, 16%), carcinogenic human papillomaviruses (N=3, 7%), and torque teno viruses (N=5, 12%). Immunohistochemistry revealed variable Large T antigen expression in BKPyV-positive tumors ranging from 100% positive staining of tumor tissue to less than 1%. In most cases of BKPyV-positive tumors, the viral genome appeared to be clonally integrated into the host chromosome consistent with microhomology-mediated end joining and coincided with focal amplifications of the tumor genome similar to other virus-mediated cancers. Significant changes in host gene expression consistent with the functions of BKPyV Large T antigen were also observed in these tumors. Lastly, we identified four mutation signatures in our cases, with those attributable to APOBEC3 and SBS5 being the most abundant. Mutation signatures associated with an antiviral drug, ganciclovir, and aristolochic acid, a nephrotoxic compound found in some herbal medicines, were also observed. The results suggest multiple pathways to carcinogenesis in solid organ transplant recipients with a large fraction being virus-associated.

Polyomavirus BK Genome Comparison Shows High Genetic Diversity in Kidney Transplant Recipients Three Months after Transplantation

BK polyomavirus (BKPyV) is a human DNA virus generally divided into twelve subgroups based on the genetic diversity of Viral Protein 1 (VP1). BKPyV can cause polyomavirus-associated nephropathy (PVAN) after kidney transplantation. Detection of BKPyV DNA in blood (viremia) is a source of concern and increase in plasma viral load is associated with a higher risk of developing PVAN. In this work, we looked for possible associations of specific BKPyV genetic features with higher plasma viral load in kidney transplant patients. We analyzed BKPyV complete genome in three-month samples from kidney recipients who developed viremia during their follow-up period. BKPyV sequences were obtained by next-generation sequencing and were de novo assembled using the new BKAnaLite pipeline. Based on the data from 72 patients, we identified 24 viral groups with unique amino acid sequences: three in the VP1 subgroup IVc2, six in Ib1, ten in Ib2, one in Ia, and four in II. In none of the groups did the mean plasma viral load reach a statistically significant difference from the overall mean observed at three months after transplantation. Further investigation is needed to better understand the link between the newly described BKPyV genetic variants and pathogenicity in kidney transplant recipients.

Rearrangement in the Hypervariable Region of JC Polyomavirus Genomes Isolated from Patient Samples and Impact on Transcription Factor-Binding Sites and Disease Outcomes

JC polyomavirus (JCPyV) is the causative agent of the fatal, incurable, neurological disease, progressive multifocal leukoencephalopathy (PML). The virus is present in most of the adult population as a persistent, asymptotic infection in the kidneys. During immunosuppression, JCPyV reactivates and invades the central nervous system. A main predictor of disease outcome is determined by mutations within the hypervariable region of the viral genome. In patients with PML, JCPyV undergoes genetic rearrangements in the noncoding control region (NCCR). The outcome of these rearrangements influences transcription factor binding to the NCCR, orchestrating viral gene transcription. This study examines 989 NCCR sequences from patient isolates deposited in GenBank to determine the frequency of mutations based on patient isolation site and disease status. The transcription factor binding sites (TFBS) were also analyzed to understand how these rearrangements could influence viral transcription. It was determined that the number of TFBS was significantly higher in PML samples compared to non-PML samples. Additionally, TFBS that could promote JCPyV infection were more prevalent in samples isolated from the cerebrospinal fluid compared to other locations. Collectively, this research describes the extent of mutations in the NCCR that alter TFBS and how they correlate with disease outcome.

In Vivo Generation of BK and JC Polyomavirus Defective Viral Genomes in Human Urine Samples Associated with Higher Viral Loads

Defective viral genomes (DVGs) are parasitic viral sequences containing point mutations, deletions, or duplications that might interfere with replication. DVGs are often associated with viral passage at high multiplicities of infection in culture systems but have been increasingly reported in clinical specimens. To date however, only RNA viruses have been shown to contain DVGs in clinical specimens. Here, using direct deep sequencing with multiple library preparation strategies and confirmatory digital droplet PCR (ddPCR) of urine samples taken from immunosuppressed individuals, we show that clinical BK polyomavirus (BKPyV) and JC polyomavirus (JCPyV) strains contain widespread genomic rearrangements across multiple loci that likely interfere with viral replication. BKPyV DVGs were derived from BKPyV genotypes Ia, Ib-1, and Ic. The presence of DVGs was associated with specimens containing higher viral loads but never reached clonality, consistent with a model of parasitized replication. These DVGs persisted during clinical infection as evidenced in two separate pairs of samples containing BK virus collected from the same individual up to 302 days apart. In a separate individual, we observed the generation of DVGs after a 57.5-fold increase in viral load. In summary, by extending the presence of DVGs in clinical specimens to DNA viruses, we demonstrate the ubiquity of DVGs in clinical virology.

IMPORTANCE Defective viral genomes (DVGs) can have a significant impact on the production of infectious virus particles. DVGs have only been identified in cultured viruses passaged at high multiplicities of infection and RNA viruses collected from clinical specimens; no DNA virus in the wild has been shown to contain DVGs. Here, we identified BK and JC polyomavirus DVGs in clinical urine specimens and demonstrated that these DVGs are more frequently identified in samples with higher viral loads. The strains containing DVGs had rearrangements throughout their genomes, with the majority affecting genes required for viral replication. Longitudinal analysis showed that these DVGs can persist during an infection but do not reach clonality within the chronically infected host. Our identification of polyomavirus DVGs suggests that these parasitic sequences exist across the many classes of viruses capable of causing human disease.

BK Polyomavirus Genotypes in Two Patients after Hematopoietic Cell Transplant

BK polyomavirus (BKPyV) infection can lead to nephropathy and hemorrhagic cystitis (HC). We evaluated BKPyV genotypes in two individuals after hematopoietic cell transplant (HCT). The first case developed HC and was infected with genotype Ib-2, while the second did not develop HC and was infected with genotype Ia.

BK polyomavirus (BKPyV) infection can lead to nephropathy and hemorrhagic cystitis (HC). We evaluated BKPyV genotypes in two individuals after hematopoietic cell transplant (HCT). The first case developed HC and was infected with genotype Ib-2, while the second did not develop HC and was infected with genotype Ia.