Exploring Canine Picornavirus Diversity in the USA Using Wastewater Surveillance: From High-Throughput Genomic Sequencing to Immuno-Informatics and Capsid Structure Modeling

The SARS-CoV-2 pandemic resulted in a scale-up of viral genomic surveillance globally. However, the wet lab constraints (economic, infrastructural, and personnel) of translating novel virus variant sequence information to meaningful immunological and structural insights that are valuable for the development of broadly acting countermeasures (especially for emerging and re-emerging viruses) remain a challenge in many resource-limited settings. Here, we describe a workflow that couples wastewater surveillance, high-throughput sequencing, phylogenetics, immuno-informatics, and virus capsid structure modeling for the genotype-to-serotype characterization of uncultivated picornavirus sequences identified in wastewater. Specifically, we analyzed canine picornaviruses (CanPVs), which are uncultivated and yet-to-be-assigned members of the family Picornaviridae that cause systemic infections in canines. We analyzed 118 archived (stored at -20 °C) wastewater (WW) samples representing a population of ~700,000 persons in southwest USA between October 2019 to March 2020 and October 2020 to March 2021. Samples were pooled into 12 two-liter volumes by month, partitioned (into filter-trapped solids [FTSs] and filtrates) using 450 nm membrane filters, and subsequently concentrated to 2 mL (1000×) using 10,000 Da MW cutoff centrifugal filters. The 24 concentrates were subjected to RNA extraction, CanPV complete capsid single-contig RT-PCR, Illumina sequencing, phylogenetics, immuno-informatics, and structure prediction. We detected CanPVs in 58.3% (14/24) of the samples generated 13,824,046 trimmed Illumina reads and 27 CanPV contigs. Phylogenetic and pairwise identity analyses showed eight CanPV genotypes (intragenotype divergence <14%) belonging to four clusters, with intracluster divergence of <20%. Similarity analysis, immuno-informatics, and virus protomer and capsid structure prediction suggested that the four clusters were likely distinct serological types, with predicted cluster-distinguishing B-cell epitopes clustered in the northern and southern rims of the canyon surrounding the 5-fold axis of symmetry. Our approach allows forgenotype-to-serotype characterization of uncultivated picornavirus sequences by coupling phylogenetics, immuno-informatics, and virus capsid structure prediction. This consequently bypasses a major wet lab-associated bottleneck, thereby allowing resource-limited settings to leapfrog from wastewater-sourced genomic data to valuable immunological insights necessary for the development of prophylaxis and other mitigation measures.