Wastewater-based surveillance (WBS) is an important epidemiological and public health tool for tracking pathogens across the scale of a building, neighbourhood, city, or region. WBS gained widespread adoption globally during the SARS-CoV-2 pandemic for estimating community infection levels by qPCR. Sequencing pathogen genes or genomes from wastewater adds information about pathogen genetic diversity, which can be used to identify viral lineages (including variants of concern) that are circulating in a local population. Capturing the genetic diversity by WBS sequencing is not trivial, as wastewater samples often contain a diverse mixture of viral lineages with real mutations and sequencing errors, which must be deconvoluted computationally from short sequencing reads. In this study we assess nine different computational tools that have recently been developed to address this challenge. We simulated 100 wastewater sequence samples consisting of SARS-CoV-2 BA.1, BA.2, and Delta lineages, in various mixtures, as well as a Delta–Omicron recombinant and a synthetic ‘novel’ lineage. Most tools performed well in identifying the true lineages present and estimating their relative abundances and were generally robust to variation in sequencing depth and read length. While many tools identified lineages present down to 1 % frequency, results were more reliable above a 5 % threshold. The presence of an unknown synthetic lineage, which represents an unclassified SARS-CoV-2 lineage, increases the error in relative abundance estimates of other lineages, but the magnitude of this effect was small for most tools. The tools also varied in how they labelled novel synthetic lineages and recombinants. While our simulated dataset represents just one of many possible use cases for these methods, we hope it helps users understand potential sources of error or bias in wastewater sequencing analysis and to appreciate the commonalities and differences across methods.

Abstract Tomato Brown Rugose Fruit Virus ( ToBRFV ) is a plant pathogen that infects important Solanaceae crop species and can dramatically reduce tomato crop yields. The ToBRFV has rapidly spread around the globe due to its ability to escape detection by antiviral host genes, most notably Tm-2 2 , which are used to confer resistance to other Tobamoviruses in tomato plants. Development of robust and reproducible methods for detecting viruses in the environment aids in the tracking and reduction of pathogen transmission. We detected ToBRFV in municipal wastewater influent (WWI) samples, likely due to its presence in human waste, demonstrating a widespread distribution of ToBRFV in WWI throughout Ontario, Canada. To aid in global ToBRFV surveillance efforts, we developed a tiled-amplicon approach to sequence and track the evolution of ToBRFV genomes in municipal WWI. Our assay recovers 97.5% of the 6393 bp ToBRFV RefSeq genome, omitting the terminal 5’ and 3’ ends. We demonstrate that our sequencing assay is a robust, sensitive, and highly specific method for recovering ToBRFV genomes. Our ToBRFV assay was developed using existing ARTIC Network resources, which includes genome specific primer design, sequencing library prep, and read analysis. Additionally, we adapted our lineage abundance estimation tool, Alcov, to estimate the abundance of ToBRFV clades in samples.

Abstract Wastewater-based surveillance (WBS) is an important epidemiological and public health tool for tracking pathogens across the scale of a building, neighbourhood, city, or region. WBS gained widespread adoption globally during the SARS-CoV-2 pandemic for estimating community infection levels by qPCR. Sequencing pathogen genes or genomes from wastewater adds information about pathogen genetic diversity which can be used to identify viral lineages (including variants of concern) that are circulating in a local population. Capturing the genetic diversity by WBS sequencing is not trivial, as wastewater samples often contain a diverse mixture of viral lineages with real mutations and sequencing errors, which must be deconvoluted computationally from short sequencing reads. In this study we assess nine different computational tools that have recently been developed to address this challenge. We simulated 100 wastewater sequence samples consisting of SARS-CoV-2 BA.1, BA.2, and Delta lineages, in various mixtures, as well as a Delta-Omicron recombinant and a synthetic “novel” lineage. Most tools performed well in identifying the true lineages present and estimating their relative abundances, and were generally robust to variation in sequencing depth and read length. While many tools identified lineages present down to 1% frequency, results were more reliable above a 5% threshold. The presence of an unknown synthetic lineage, which represents an unclassified SARS-CoV-2 lineage, increases the error in relative abundance estimates of other lineages, but the magnitude of this effect was small for most tools. The tools also varied in how they labelled novel synthetic lineages and recombinants. While our simulated dataset represents just one of many possible use cases for these methods, we hope it helps users understand potential sources of noise or bias in wastewater sequencing data and to appreciate the commonalities and differences across methods.

Horizontal gene transfer events between viruses and hosts are widespread across the virosphere. In cyanophage-host systems, such events often involve the transfer of genes involved in photosynthetic processes. The genome of the lytic cyanomyovirus Ma-LMM01 infecting the toxic, bloom-forming, freshwater Microcystis aeruginosa NIES-298 contains a homolog of the non-bleaching A ( nblA ) gene, which was probably acquired from its host. The function of the NblA protein is to disassemble phycobilisomes, cyanobacterial light harvesting complexes that can comprise up to half of the cellular soluble protein content. NblA thus plays an essential dual role in cyanobacteria: it protects the cell from high light intensities and increases the intracellular nitrogen pool under nutrient limitation. NblA has previously been shown to interact with phycocyanin, one of the main components of phycobilisomes. Using structural modeling and protein-protein docking, we show that the NblA dimer of Ma-LMM01 is predicted to have a significantly higher binding affinity for M. aeruginosa NIES-298 phycocyanin (αβ) 6 hexamers, compared to the host homolog. Protein-protein docking suggests that the viral NblA structural model is able to bind deeper into the phycocyanin groove. The main structural difference between the virus and host NblA appears to be an additional α-helix near the N-terminus of the viral NblA, which could be partly responsible for the deeper binding into phycocyanin. This unique helical region, absent in the cellular NblA, would be expected to constitute a viral evolutionary innovation. We propose that a higher binding affinity of NblA to the host phycocyanin may represent a selective advantage for the virus, whose rapid infection cycle requires an increased phycobilisome degradation rate that is not fulfilled by the NblA of the host.

ABSTRACT Through infection and lysis of their coexisting bacterial hosts, viruses impact the biogeochemical cycles sustaining globally significant pelagic oceanic ecosystems. Currently, little is known of the ecological interactions between lytic viruses and their bacterial hosts underlying these biogeochemical impacts at ecosystem scales. This study focused on populations of lytic viruses carrying the B 12 - dependent Class II monomeric ribonucleotide reductase (RNR) gene, ribonucleotide triphosphate reductase (RTPR), documenting seasonal changes in pelagic virioplankton and bacterioplankton using amplicon sequences of RTPR and the 16S rRNA gene, respectively. Amplicon sequence libraries were analyzed using compositional data analysis tools that account for the compositional nature of these data. Both virio- and bacterioplankton communities responded to environmental changes typically seen across seasonal cycles as well as shorter term upwelling–downwelling events. Defining RTPR-carrying viral populations according to major phylogenetic clades proved a more robust means of exploring virioplankton ecology than operational taxonomic units defined by percent sequence homology. Virioplankton RTPR populations showed positive associations with a broad phylogenetic diversity of bacterioplankton including dominant taxa within pelagic oceanic ecosystems such as Prochlorococcus and SAR11. Temporal changes in RTPR-virioplankton, occurring as both free viruses and within infected cells, indicated possible viral–host pairs undergoing sustained infection and lysis cycles throughout the seasonal study. Phylogenetic relationships inferred from RTPR sequences mirrored ecological patterns in virio- and bacterioplankton populations demonstrating possible genome to phenome associations for an essential viral replication gene.

ABSTRACT Toxic algal bloom-forming cyanobacteria are a persistent problem globally for many aquatic environments. Their occurrence is attributed to eutrophication and rising temperatures due to climate change. The result of these blooms is often loss in biodiversity, economic impacts on tourism and fisheries, and risks to human and animal health. Of particular concern is the poorly understood interplay between viruses and toxic species that form blooms because viruses may exacerbate their harmful effects. Concurrently, cyanobacteria are also a source of bioactive compounds other than toxins, which makes them good candidates for drug discovery. We show that virus infection of the cyanobacterium Microcystis aeruginosa , results in as high as a 40-fold increase in the toxin microcystin two days post virus infection (dpi), and predict that microcystin levels may remain high in a body of water up to 7 dpi, long after water discoloration and cell lysis. This implicates viruses as major contributors to toxin release from cyanobacteria and emphasizes the importance of taking them into account in predictive models and in the assessment of water safety. We also show that bioactive compounds of M. aeruginosa inhibit and delay infection of single stranded DNA and single stranded RNA viruses. This highlights the potential of cyanobacteria as an excellent source for the discovery of novel antiviral compounds, and the ease with which screening for cyanobacterial antivirals can be achieved.