Comparison of ultrafiltration and iron chloride flocculation in the preparation of aquatic viromes from contrasting sample types

Accurate prediction of absolute prokaryotic abundance from DNA concentration

Quantification of the absolute microbial abundance in a human stool sample is crucial for a comprehensive understanding of the microbial ecosystem, but this information is lost upon metagenomic sequencing. While several methods exist to measure absolute microbial abundance, they are technically challenging and costly, presenting an opportunity for machine learning. Here, we observe a strong correlation between DNA concentration and the absolute number of 16S ribosomal RNA copies as measured by digital droplet PCR in clinical stool samples from individuals undergoing hematopoietic cell transplantation (BMT CTN 1801). Based on this correlation and additional measurements, we trained an accurate yet simple machine learning model for the prediction of absolute prokaryotic load, which showed exceptional prediction accuracy on an external cohort that includes people living with Parkinson's disease and healthy controls. We propose that, with further validation, this model has the potential to enable accurate absolute abundance estimation based on readily available sample measurements.

Development of a quantitative metagenomic approach to establish quantitative limits and its application to viruses

Quantitative metagenomic methods are maturing but continue to lack clearly-defined analytical limits. Here, we developed a computational tool, QuantMeta, to determine the absolute abundance of targets in metagenomes spiked with synthetic DNA standards. The tool establishes (i) entropy-based detection thresholds to confidently determine the presence of targets, and (ii) an approach to identify and correct read mapping or assembly errors and thus improve the quantification accuracy. Together this allows for an approach to confidently quantify absolute abundance of targets, be they microbial populations, genes, contigs, or metagenome-assembled genomes. We applied the approach to quantify single- and double-stranded DNA viruses in wastewater viral metagenomes, including pathogens and bacteriophages. Concentrations of total DNA viruses in wastewater influent and effluent were >108 copies/ml using QuantMeta. Human-associated DNA viruses were detected and quantifiable with QuantMeta thresholds, including polyomavirus, papillomavirus, and crAss-like phages, at concentrations similar to previous reports that utilized quantitative polymerase chain reaction (PCR)-based assays. Our results highlight the higher detection thresholds of quantitative metagenomics (approximately 500 copies/μl) as compared to PCR-based quantification (approximately 10 copies/μl) despite a sequencing depth of 200 million reads per sample. The QuantMeta approach, applicable to both viral and cellular metagenomes, advances quantitative metagenomics by improving the accuracy of measured target absolute abundances.

Method evaluation for viruses in activated sludge: Concentration, sequencing, and identification

Activated sludge (AS) in wastewater treatment plants is one of the largest artificial microbial ecosystems on earth and it makes enormous contributions to human societies. Viruses are an important component in AS with a high abundance. However, their communities and functionalities have not been as widely explored as those of other microorganisms, such as bacteria. This gap is mainly due to technical challenges in effective viral concentration, extraction, and sequencing. In this study, we compared four kinds of concentration methods, two sequencing approaches, and four identification bioinformatic tools to evaluate the whole analysis workflow for viruses in AS. Results showed flocculation, filtration, and resuspension (FFR) could get the longest DNA lengths and ultracentrifugation obtained the highest DNA yields for viruses in AS. Based on the results of present study, FFR and tangential flow filtration with the membrane pore size of 100 kDa were most recommended to concentrate viruses in AS samples with huge volumes. Besides, different concentration methods could get different viral catalogs and thus multiple methods should be combined to get the whole picture of viruses in the system. In addition, geNomad was the most recommended identification tool for viruses in the present study and the long-read sequencing could improve the assembly statistics of viruses when compared with the short-read sequencing. For the 8192 viral operational taxonomic units in this study, 95.1 % of them were phages and belonged to the same lineage at the order level of Caudovirales. Virulent phages dominated the AS system and Pseudomonadota were the main host. Taken together, this study provides new insights into methods selection for virus research of AS.

Surface microlayer-mediated virome dissemination in the Central Arctic

BACKGROUND

Aquatic viruses act as key players in shaping microbial communities. In polar environments, they face significant challenges such as limited host availability and harsh conditions. However, due to the restricted accessibility of these ecosystems, our understanding of viral diversity, abundance, adaptations, and host interactions remains limited.

RESULTS

To fill this knowledge gap, we studied viruses from atmosphere-close aquatic ecosystems in the Central Arctic and Northern Greenland. Aquatic samples for virus-host analysis were collected from ~60 cm depth and the submillimeter surface microlayer (SML) during the Synoptic Arctic Survey 2021 on icebreaker Oden in the Arctic summer. Water was sampled from a melt pond and open water before undergoing size-fractioned filtration, followed by genome-resolved metagenomic and cultivation investigations. The prokaryotic diversity in the melt pond was considerably lower compared to that of open water. The melt pond was dominated by a Flavobacterium sp. and Aquiluna sp., the latter having a relatively small genome size of 1.2 Mb and the metabolic potential to generate ATP using the phosphate acetyltransferase-acetate kinase pathway. Viral diversity on the host fraction (0.2-5 µm) of the melt pond was strikingly limited compared to that of open water. From the 1154 viral operational taxonomic units (vOTUs), of which two-thirds were predicted bacteriophages, 17.2% encoded for auxiliary metabolic genes (AMGs) with metabolic functions. Some AMGs like glycerol-3-phosphate cytidylyltransferase and ice-binding like proteins might serve to provide cryoprotection for the host. Prophages were often associated with SML genomes, and two active prophages of new viral genera from the Arctic SML strain Leeuwenhoekiella aequorea Arc30 were induced. We found evidence that vOTU abundance in the SML compared to that of ~60 cm depth was more positively correlated with the distribution of a vOTU across five different Arctic stations.

CONCLUSIONS

The results indicate that viruses employ elaborate strategies to endure in extreme, host-limited environments. Moreover, our observations suggest that the immediate air-sea interface serves as a platform for viral distribution in the Central Arctic. Video Abstract.

Effects of parental care on skin microbial community composition in poison frogs

Parent-offspring interactions constitute the first contact of many newborns with their environment, priming community assembly of microbes through priority effects. Early exposure to microbes can have lasting influences on the assembly and functionality of the host's microbiota, leaving a life-long imprint on host health and disease. Studies of the role played by parental care in microbial acquisition have primarily focused on humans and hosts with agricultural relevance. Anuran vertebrates offer the opportunity to examine microbial community composition across life stages as a function of parental investment. In this study, we investigate vertical transmission of microbiota during parental care in a poison frog (Family Dendrobatidae), where fathers transport their offspring piggyback-style from terrestrial clutches to aquatic nurseries. We found that substantial bacterial colonization of the embryo begins after hatching from the vitelline envelope, emphasizing its potential role as microbial barrier during early development. Using a laboratory cross-foster experiment, we demonstrated that poison frogs performing tadpole transport serve as a source of skin microbes for tadpoles on their back. To study how transport impacts the microbial skin communities of tadpoles in an ecologically relevant setting, we sampled frogs and tadpoles of sympatric species that do or do not exhibit tadpole transport in their natural habitat. We found more diverse microbial communities associated with tadpoles of transporting species compared to a non-transporting frog. However, we detected no difference in the degree of similarity between adult and tadpole skin microbiotas, based on whether the frog species exhibits transporting behavior or not. Using a field experiment, we confirmed that tadpole transport can result in the persistent colonization of tadpoles by isolated microbial taxa associated with the caregiver's skin, albeit often at low prevalence. This is the first study to describe vertical transmission of skin microbes in anuran amphibians, showing that offspring transport may serve as a mechanism for transmission of parental skin microbes. Overall, these findings provide a foundation for further research on how vertical transmission in this order impacts host-associated microbiota and physiology.

Large-scale single-virus genomics uncovers hidden diversity of river water viruses and diversified gene profiles

Environmental viruses (primarily bacteriophages) are widely recognized as playing an important role in ecosystem homeostasis through the infection of host cells. However, the majority of environmental viruses are still unknown as their mosaic structure and frequent mutations in their sequences hinder genome construction in current metagenomics. To enable the large-scale acquisition of environmental viral genomes, we developed a new single-viral genome sequencing platform with microfluidic-generated gel beads. Amplification of individual DNA viral genomes in mass-produced gel beads allows high-throughput genome sequencing compared to conventional single-virus genomics. The sequencing analysis of river water samples yielded 1431 diverse viral single-amplified genomes, whereas viral metagenomics recovered 100 viral metagenome-assembled genomes at the comparable sequence depth. The 99.5% of viral single-amplified genomes were determined novel at the species level, most of which could not be recovered by a metagenomic assembly. The large-scale acquisition of diverse viral genomes identified protein clusters commonly detected in different viral strains, allowing the gene transfer to be tracked. Moreover, comparative genomics within the same viral species revealed that the profiles of various methyltransferase subtypes were diverse, suggesting an enhanced escape from host bacterial internal defense mechanisms. Our use of gel bead-based single-virus genomics will contribute to exploring the nature of viruses by accelerating the accumulation of draft genomes of environmental DNA viruses.

Ecogenomics and cultivation reveal distinctive viral-bacterial communities in the surface microlayer of a Baltic Sea slick

Visible surface films, termed slicks, can extensively cover freshwater and marine ecosystems, with coastal regions being particularly susceptible to their presence. The sea-surface microlayer (SML), the upper 1-mm at the air-water interface in slicks (herein slick SML) harbors a distinctive bacterial community, but generally little is known about SML viruses. Using flow cytometry, metagenomics, and cultivation, we characterized viruses and bacteria in a brackish slick SML in comparison to non-slick SML as well as seawater below slick and non-slick areas (subsurface water = SSW). Size-fractionated filtration of all samples distinguished viral attachment to hosts and particles. The slick SML contained higher abundances of virus-like particles, prokaryotic cells, and dissolved organic carbon compared to non-slick SML and SSW. The community of 428 viral operational taxonomic units (vOTUs), 426 predicted as lytic, distinctly differed across all size fractions in the slick SML compared to non-slick SML and SSW. Specific metabolic profiles of bacterial metagenome-assembled genomes and isolates in the slick SML included a prevalence of genes encoding motility and carbohydrate-active enzymes (CAZymes). Several vOTUs were enriched in slick SML, and many virus variants were associated with particles. Nine vOTUs were only found in slick SML, six of them being targeted by slick SML-specific clustered-regularly interspaced short palindromic repeats (CRISPR) spacers likely originating from Gammaproteobacteria. Moreover, isolation of three previously unknown lytic phages for Alishewanella sp. and Pseudoalteromonas tunicata, abundant and actively replicating slick SML bacteria, suggests that viral activity in slicks contributes to biogeochemical cycling in coastal ecosystems.

Robust Performance of SARS-CoV-2 Whole-Genome Sequencing from Wastewater with a Nonselective Virus Concentration Method

The sequencing of human virus genomes from wastewater samples is an efficient method for tracking viral transmission and evolution at the community level. However, this requires the recovery of viral nucleic acids of high quality. We developed a reusable tangential-flow filtration system to concentrate and purify viruses from wastewater for genome sequencing. A pilot study was conducted with 94 wastewater samples from four local sewersheds, from which viral nucleic acids were extracted, and the whole genome of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) was sequenced using the ARTIC V4.0 primers. Our method yielded a high probability (0.9) of recovering complete or near-complete SARS-CoV-2 genomes (>90% coverage at 10× depth) from wastewater when the COVID-19 incidence rate exceeded 33 cases per 100 000 people. The relative abundances of sequenced SARS-CoV-2 variants followed the trends observed from patient-derived samples. We also identified SARS-CoV-2 lineages in wastewater that were underrepresented or not present in the clinical whole-genome sequencing data. The developed tangential-flow filtration system can be easily adopted for the sequencing of other viruses in wastewater, particularly those at low concentrations.

SARS-CoV-2 Variants Detection Strategies in Wastewater Samples Collected in the Bangkok Metropolitan Region

Wastewater surveillance is considered a promising approach for COVID-19 surveillance in communities. In this study, we collected wastewater samples between November 2020 and February 2022 from twenty-three sites in the Bangkok Metropolitan Region to detect the presence of SARS-CoV-2 and its variants for comparison to standard clinical sampling. A total of 215 wastewater samples were collected and tested for SARS-CoV-2 RNA by real-time PCR with three targeted genes (N, E, and ORF1ab); 102 samples were positive (42.5%). The SARS-CoV-2 variants were determined by a multiplex PCR MassARRAY assay to distinguish four SARS-CoV-2 variants, including Alpha, Beta, Delta, and Omicron. Multiple variants of Alpha-Delta and Delta-Omicron were detected in the wastewater samples in July 2021 and January 2022, respectively. These wastewater variant results mirrored the country data from clinical specimens deposited in GISAID. Our results demonstrated that wastewater surveillance using multiple signature mutation sites for SARS-CoV-2 variant detection is an appropriate strategy to monitor the presence of SARS-CoV-2 variants in the community at a low cost and with rapid turn-around time. However, it is essential to note that sequencing surveillance of wastewater samples should be considered complementary to whole genome sequencing of clinical samples to detect novel variants.

In Vitro Viral Recovery Yields under Different Re-Suspension Buffers in Iron Flocculation to Concentrate Viral Hemorrhagic Septicemia Virus Genotype IVa in Seawater

Iron flocculation is widely used to concentrate viruses in water, followed by Fe-virus flocculate formation, collection, and elution. In the elution stage, an oxalic or ascorbic acid re-suspension buffer dissolved iron hydroxide. After the concentration of viral hemorrhagic septicemia virus (VHSV) in seawater (1 × 10¹ to 1 × 10⁵ viral genome copies or plaque-forming unit (PFU)/mL), the recovery yield of the viral genome using quantitative real-time PCR (qRT-PCR) and viral infectivity using the plaque assay were investigated to evaluate the validity of the two re-suspension buffers to concentrate VHSV. The mean viral genome recovery yield with oxalic and ascorbic acid was 71.2 ± 12.3% and 81.4 ± 9.5%, respectively. The mean viral infective recovery yields based on the PFU were significantly different between the two buffers at 23.8 ± 22.7% (oxalic acid) and 4.4 ± 2.7% (ascorbic acid). Notably, although oxalic acid maintains viral infectivity over 60% at a viral concentration above 10⁵ PFU/mL, the infective VHSVs were not sufficiently recovered at a low viral concentration (10² PFU/mL, <10%). To support this result, concentrated VHSV was inoculated in Epithelioma papulosum cyprini (EPC) cells to confirm cell viability, viral gene expression, and extracellular viral titer. All results demonstrated that oxalic acid buffer was superior to ascorbic acid buffer in preserving viral infectivity.

Surfactant Treatment for Efficient Gene Detection of Enteric Viruses and Indicators in Surface Water Concentrated by Ultrafiltration

The hollow fiber ultrafiltration (HFUF)-based microbial concentration method is widely applied for monitoring pathogenic viruses and microbial indicators in environmental water samples. However, the HFUF-based method can co-concentrate substances that interfere with downstream molecular processes-nucleic acid extraction, reverse transcription (RT), and PCR. These inhibitory substances are assumed to be hydrophobic and, therefore, expected to be excluded by a simple surfactant treatment before the silica membrane-based RNA extraction process. In this study, the efficacy and limitations of the sodium deoxycholate (SD) treatment were assessed by quantifying a process control and indigenous viruses using 42 surface water samples concentrated with HFUF. With some exceptions, which tended to be seen in samples with high turbidity (> 4.0 NTU), virus recovery by the ultrafiltration method was sufficiently high (> 10%). RNA extraction-RT-quantitative PCR (RT-qPCR) efficiency of the process control was insufficient (10%) for 30 of the 42 HFUF concentrates without any pretreatments, but it was markedly improved for 21 of the 30 inhibitory concentrates by the SD treatment. Detection rates of indigenous viruses were also improved and no substantial loss of viral RNA was observed. The SD treatment was particularly effective in mitigating RT-qPCR inhibition, although it was not effective in improving RNA extraction efficiency. The methodology is simple and easily applied. These findings indicate that SD treatment can be a good alternative to sample dilution, which is widely applied to mitigate the effect of RT-qPCR inhibition, and can be compatible with other countermeasures.

A rapid, high-throughput, and sensitive PEG-precipitation method for SARS-CoV-2 wastewater surveillance

The effective application of wastewater surveillance is dependent on testing capacity and sensitivity to obtain high spatial resolution testing results for a timely targeted public health response. To achieve this purpose, the development of rapid, high-throughput, and sensitive virus concentration methods is urgently needed. Various protocols have been developed and implemented in wastewater surveillance networks so far, however, most of them lack the ability to scale up testing capacity or cannot achieve sufficient sensitivity for detecting SARS-CoV-2 RNA at low prevalence. In the present study, using positive raw wastewater in Hong Kong, a PEG precipitation-based three-step centrifugation method was developed, including low-speed centrifugation for large particles removal and the recovery of viral nucleic acid, and medium-speed centrifugation for the concentration of viral nucleic acid. This method could process over 100 samples by two persons per day to reach the process limit of detection (PLoD) of 3286 copies/L wastewater. Additionally, it was found that the testing capacity could be further increased by decreasing incubation and centrifugation time without significantly influencing the method sensitivity. The entire procedure uses ubiquitous reagents and instruments found in most laboratories to obtain robust testing results. This high-throughput, cost-effective, and sensitive tool will promote the establishment of nearly real-time wastewater surveillance networks for valuable public health information.

Chloride Enhances DNA Reactivity with Chlorine under Conditions Relevant to Water Treatment

Free available chlorine (FAC) is widely used to inactivate viruses by oxidizing viral components, including genomes. It is commonly assumed that hypochlorous acid (HOCl) is the chlorinating agent responsible for virus inactivation; however, recent studies have underscored that minor constituents of FAC existing in equilibrium with HOCl, such as molecular chlorine (Cl₂), can influence FAC reactivity toward select organic compounds. This study measures the FAC reaction kinetics with dsDNA and ssDNA extracted from representative bacteriophages (T3 and ϕX174) in samples augmented with chloride. Herein, chloride enhances FAC reactivity toward dsDNA and, to a lesser extent, toward ssDNA, especially at pH < 7.5. The enhanced reactivity can be attributed to the formation of Cl₂. Second-order rate constants were determined for reactions of ssDNA and dsDNA with HOCl and Cl₂. DNA chlorination kinetics followed the reactivity-selectivity principle, where the more-reactive nucleophilic species (ssDNA, ∼100× more reactive than dsDNA) reacted less selectively with electrophilic FAC species. The addition of chloride was also shown to enhance the inactivation of bacteriophage T3 (dsDNA genome) by FAC but did not enhance the inactivation of bacteriophage ϕX174 (ssDNA genome). Overall, the results suggest that Cl₂ is an important chlorinating agent of nucleic acids and viruses.

Optimization and evaluation of viral metagenomic amplification and sequencing procedures toward a genome-level resolution of the human fecal DNA virome

INTRODUCTION

Viruses in the human gut have been linked to health and disease. Deciphering the gut virome is dependent on metagenomic sequencing of the virus-like particles (VLPs) purified from the fecal specimens. A major limitation of conventional viral metagenomic sequencing is the low recoverability of viral genomes from the metagenomic dataset.

OBJECTIVES

To develop an optimal method for viral amplification and metagenomic sequencing for maximizing the recovery of viral genomes.

METHODS

We performed parallel virus enrichment and DNA extraction to generate ∼ 30 viral DNA samples from each of 5 fresh fecal specimens and conducted the experiments including 1) optimizing the cycle number for high-fidelity enzyme-based PCR amplification, 2) evaluating the reproducibility of the optimally whole viral metagenomic experimental process, 3) evaluating the reliability of multiple displacement amplification (MDA), 4) testing the capability of long-read sequencing for improving viral metagenomic assembly, and 5) comparing the differences between viral metagenomic and bulk metagenomic approaches.

RESULTS

Our results revealed that the optimal cycle number for PCR amplification is 15. We verified the reliability of MDA and the effectiveness of long-read sequencing. Based on our optimized results, we generated 151 high-quality viruses using the dataset combined from short-read and long-read sequencing. Genomic analysis of these viruses found that most (60.3%) of them were previously unknown and showed a remarkable diversity of viral functions, especially the existence of 206 viral auxiliary metabolic genes. Finally, we uncovered significant differences in the efficiency and coverage of viral identification between viral metagenomic and bulk metagenomic approaches.

CONCLUSIONS

Our study demonstrates the potential of optimized experiment and sequencing strategies in uncovering viral genomes from fecal specimens, which will facilitate future research about the genome-level characterization of complex viral communities.

Evaluation of a virus concentration method based on ultrafiltration and wet foam elution for studying viruses from large-volume water samples

Assessing the presence of viruses in large-volume samples involves cumbersome methods that require specialized training and laboratory equipment. In this study, a large volume concentration (LVC) method, based on dead-end ultrafiltration (DEUF) and Wet Foam Elution™ technology, was evaluated in different type of waters and different microorganisms. Its recovery efficiency was evaluated through different techniques (infectivity assays and molecular detection) by spiking different viral surrogates (bacteriophages PhiX174 and MS2 and Coxsackie virus B5 (CVB5) and Escherichia coli (E. coli). Furthermore, the application of a secondary concentration step was evaluated and compared with skimmed milk flocculation. Viruses present in river water, seawater and groundwater samples were concentrated by applying LVC method and a centrifugal ultrafiltration device (CeUF), as a secondary concentration step and quantified with specific qPCR Human adenoviruses (HAdV) and noroviruses (NoVs). MS2 was used as process control, obtaining a mean viral recovery of 22.0 ± 12.47%. The presence of other viruses was also characterized by applying two different next-generation sequencing approaches. LVC coupled to a secondary concentration step based on CeUF allowed to detect naturally occurring viruses such as HAdV and NoVs in different water matrices. Using HAdV as a human fecal indicator, the highest viral pollution was found in river water samples (100% of positive samples), followed by seawater (83.33%) and groundwater samples (66.67%). The LVC method has also proven to be useful as a virus concentration method in the filed since HAdV and NoVs were detected in the river water and groundwater samples concentrated in the field. All in all, LVC method presents high concentration factor and a low limit of detection and provides viral concentrates useful for subsequent molecular analysis such as PCR and massive sequencing.

A snapshot of the global drinking water virome: Diversity and metabolic potential vary with residual disinfectant use

Viruses are important drivers of microbial community ecology and evolution, influencing microbial mortality, metabolism, and horizontal gene transfer. However, the effects of viruses remain largely unknown in many environments, including in drinking water systems. Drinking water metagenomic studies have offered a whole community perspective of bacterial impacts on water quality, but have not yet considered the influences of viruses. In this study, we address this gap by mining viral DNA sequences from publicly available drinking water metagenomes from distribution systems in six countries around the world. These datasets provide a snapshot of the taxonomic diversity and metabolic potential of the global drinking water virome; and provide an opportunity to investigate the effects of geography, climate, and drinking water treatment practices on viral diversity. Both environmental conditions and differences in sample processing were found to influence the viral composition. Using free chlorine as the residual disinfectant was associated with clear differences in viral taxonomic diversity and metabolic potential, with significantly fewer viral populations and less even viral community structures than observed in distribution systems without residual disinfectant. Additionally, drinking water viruses carry antibiotic resistance genes (ARGs), as well as genes to survive oxidative stress and nitrogen limitation. Through this study, we have demonstrated that viral communities are diverse across drinking water systems and vary with the use of residual disinfectant. Our findings offer directions for future research to develop a more robust understanding of how virus-bacteria interactions in drinking water distribution systems affect water quality.

Detection of SARS-CoV-2 Variants Mu, Beta, Gamma, Lambda, Delta, Alpha, and Omicron in Wastewater Settled Solids Using Mutation-Specific Assays Is Associated with Regional Detection of Variants in Clinical Samples

Changes in the circulation of SARS-CoV-2 variants of concern (VOCs) may require changes in the public health response to the COVID-19 pandemic, as they have the potential to evade vaccines and pharmaceutical interventions and may be more transmissive than other SARS-CoV-2 variants. As such, it is essential to track and prevent their spread in susceptible communities. We developed digital reverse transcription (RT)-PCR assays for mutations characteristic of VOCs and used them to quantify those mutations in samples of wastewater settled solids collected from a publicly owned treatment works (POTW) during different phases of the COVID-19 pandemic. Wastewater concentrations of single mutations characteristic of each VOC, normalized by the concentration of a conserved SARS-CoV-2 N gene, correlate with regional estimates of the proportion of clinical infections caused by each VOC. These results suggest that targeted RT-PCR assays can be used to detect variants circulating in communities and inform the public health response to the pandemic. IMPORTANCE Wastewater represents a pooled biological sample of the contributing community and thus a resource for assessing community health. Here, we show that emergence, spread, and disappearance of SARS-CoV-2 infections caused by variants of concern are reflected in the presence of variant genomic RNA in wastewater settled solids. This work highlights an important public health use case for wastewater.

Detection of SARS-CoV-2 Variants Mu, Beta, Gamma, Lambda, Delta, Alpha, and Omicron in Wastewater Settled Solids Using Mutation-Specific Assays Is Associated with Regional Detection of Variants in Clinical Samples

Changes in the circulation of SARS-CoV-2 variants of concern (VOCs) may require changes in the public health response to the COVID-19 pandemic, as they have the potential to evade vaccines and pharmaceutical interventions and may be more transmissive than other SARS-CoV-2 variants. As such, it is essential to track and prevent their spread in susceptible communities. We developed digital reverse transcription (RT)-PCR assays for mutations characteristic of VOCs and used them to quantify those mutations in samples of wastewater settled solids collected from a publicly owned treatment works (POTW) during different phases of the COVID-19 pandemic. Wastewater concentrations of single mutations characteristic of each VOC, normalized by the concentration of a conserved SARS-CoV-2 N gene, correlate with regional estimates of the proportion of clinical infections caused by each VOC. These results suggest that targeted RT-PCR assays can be used to detect variants circulating in communities and inform the public health response to the pandemic.

IMPORTANCE Wastewater represents a pooled biological sample of the contributing community and thus a resource for assessing community health. Here, we show that emergence, spread, and disappearance of SARS-CoV-2 infections caused by variants of concern are reflected in the presence of variant genomic RNA in wastewater settled solids. This work highlights an important public health use case for wastewater.