Tracing the origin of SARS-CoV-2 omicron-like spike sequences detected in an urban sewershed: a targeted, longitudinal surveillance study of a cryptic wastewater lineage

BACKGROUND

The origin of novel SARS-CoV-2 spike sequences found in wastewater, without corresponding detection in clinical specimens, remains unclear. We sought to determine the origin of one such cryptic wastewater lineage by tracking and characterising its persistence and genomic evolution over time.

METHODS

We first detected a cryptic lineage, WI-CL-001, in municipal wastewater in Wisconsin, USA, in January, 2022. To determine the source of WI-CL-001, we systematically sampled wastewater from targeted sub-sewershed lines and maintenance holes using compositing autosamplers. Viral concentrations in wastewater samples over time were measured by RT digital PCR. In addition to using metagenomic 12s rRNA sequencing to determine the virus's host species, we also sequenced SARS-CoV-2 spike receptor binding domains, and, where possible, whole viral genomes to identify and characterise the evolution of this lineage.

FINDINGS

We traced WI-CL-001 to its source at a single commercial building. There we detected the cryptic lineage at concentrations as high as 2·7 × 109 genome copies per L. The majority of 12s rRNA sequences detected in wastewater leaving the identified source building were human. Additionally, we generated over 100 viral receptor binding domain and whole-genome sequences from wastewater samples containing the cryptic lineage collected over the 13 consecutive months this virus was detectable (January, 2022, to January, 2023). These sequences contained a combination of fixed nucleotide substitutions characteristic of Pango lineage B.1.234, which circulated in humans in Wisconsin at low levels from October, 2020, to February, 2021. Despite this, mutations in the spike gene and elsewhere resembled those subsequently found in omicron variants.

INTERPRETATION

We propose that prolonged detection of WI-CL-001 in wastewater indicates persistent shedding of SARS-CoV-2 from a single human initially infected by an ancestral B.1.234 virus. The accumulation of convergent omicron-like mutations in WI-CL-001's ancestral B.1.234 genome probably reflects persistent infection and extensive within-host evolution. People who shed cryptic lineages could be an important source of highly divergent viruses that sporadically emerge and spread.

FUNDING

The Rockefeller Foundation, Wisconsin Department of Health Services, Centers for Disease Control and Prevention, National Institute on Drug Abuse, and the Center for Research on Influenza Pathogenesis and Transmission.

Wolbachia-mediated resistance to Zika virus infection in Aedes aegypti is dominated by diverse transcriptional regulation and weak evolutionary pressures

A promising candidate for arbovirus control and prevention relies on replacing arbovirus-susceptible Aedes aegypti populations with mosquitoes that have been colonized by the intracellular bacterium Wolbachia and thus have a reduced capacity to transmit arboviruses. This reduced capacity to transmit arboviruses is mediated through a phenomenon referred to as pathogen blocking. Pathogen blocking has primarily been proposed as a tool to control dengue virus (DENV) transmission, however it works against a range of viruses, including Zika virus (ZIKV). Despite years of research, the molecular mechanisms underlying pathogen blocking still need to be better understood. Here, we used RNA-seq to characterize mosquito gene transcription dynamics in Ae. aegypti infected with the wMel strain of Wolbachia that are being released by the World Mosquito Program in Medellín, Colombia. Comparative analyses using ZIKV-infected, uninfected tissues, and mosquitoes without Wolbachia revealed that the influence of wMel on mosquito gene transcription is multifactorial. Importantly, because Wolbachia limits, but does not completely prevent, replication of ZIKV and other viruses in coinfected mosquitoes, there is a possibility that these viruses could evolve resistance to pathogen blocking. Therefore, to understand the influence of Wolbachia on within-host ZIKV evolution, we characterized the genetic diversity of molecularly barcoded ZIKV virus populations replicating in Wolbachia-infected mosquitoes and found that within-host ZIKV evolution was subject to weak purifying selection and, unexpectedly, loose anatomical bottlenecks in the presence and absence of Wolbachia. Together, these findings suggest that there is no clear transcriptional profile associated with Wolbachia-mediated ZIKV restriction, and that there is no evidence for ZIKV escape from this restriction in our system.

Wolbachia -mediated resistance to Zika virus infection in Aedes aegypti is dominated by diverse transcriptional regulation and weak evolutionary pressures

A promising candidate for arbovirus control and prevention relies on replacing arbovirus-susceptible Aedes aegypti populations with mosquitoes that have been colonized by the intracellular bacterium Wolbachia and thus have a reduced capacity to transmit arboviruses. This reduced capacity to transmit arboviruses is mediated through a phenomenon referred to as pathogen blocking. Pathogen blocking has primarily been proposed as a tool to control dengue virus (DENV) transmission, however it works against a range of viruses, including Zika virus (ZIKV). Despite years of research, the molecular mechanisms underlying pathogen blocking still need to be better understood. Here, we used RNA-seq to characterize mosquito gene transcription dynamics in Ae. aegypti infected with the w Mel strain of Wolbachia that are being released by the World Mosquito Program in Medellín, Colombia. Comparative analyses using ZIKV-infected, uninfected tissues, and mosquitoes without Wolbachia revealed that the influence of w Mel on mosquito gene transcription is multifactorial. Importantly, because Wolbachia limits, but does not completely prevent, replication of ZIKV and other viruses in coinfected mosquitoes, there is a possibility that these viruses could evolve resistance to pathogen blocking. Therefore, to understand the influence of Wolbachia on within-host ZIKV evolution, we characterized the genetic diversity of molecularly barcoded ZIKV virus populations replicating in Wolbachia -infected mosquitoes and found that within-host ZIKV evolution was subject to weak purifying selection and, unexpectedly, loose anatomical bottlenecks in the presence and absence of Wolbachia . Together, these findings suggest that there is no clear transcriptional profile associated with Wolbachia -mediated ZIKV restriction, and that there is no evidence for ZIKV escape from this restriction in our system.

Author Summary

When Wolbachia bacteria infect Aedes aegypti mosquitoes, they dramatically reduce the mosquitoes' susceptibility to infection with a range of arthropod-borne viruses, including Zika virus (ZIKV). Although this pathogen-blocking effect has been widely recognized, its mechanisms remain unclear. Furthermore, because Wolbachia limits, but does not completely prevent, replication of ZIKV and other viruses in coinfected mosquitoes, there is a possibility that these viruses could evolve resistance to Wolbachia -mediated blocking. Here, we use host transcriptomics and viral genome sequencing to examine the mechanisms of ZIKV pathogen blocking by Wolbachia and viral evolutionary dynamics in Ae. aegypti mosquitoes. We find complex transcriptome patterns that do not suggest a single clear mechanism for pathogen blocking. We also find no evidence that Wolbachia exerts detectable selective pressures on ZIKV in coinfected mosquitoes. Together our data suggest that it may be difficult for ZIKV to evolve Wolbachia resistance, perhaps due to the complexity of the pathogen blockade mechanism.