Detecting SARS-CoV-2 lineages and mutational load in municipal wastewater and a use-case in the metropolitan area of Thessaloniki, Greece

A gold standard dataset and evaluation of methods for lineage abundance estimation from wastewater

During the SARS-CoV-2 pandemic, genome-based wastewater surveillance sequencing has been a powerful tool for public health to monitor circulating and emerging viral variants. As a medium, wastewater is very complex because of its mixed matrix nature, which makes the deconvolution of wastewater samples more difficult. Here we introduce a gold standard dataset constructed from synthetic viral control mixtures of known composition, spiked into a wastewater RNA matrix and sequenced on the Oxford Nanopore Technologies platform. We compare the performance of eight of the most commonly used deconvolution tools in identifying SARS-CoV-2 variants present in these mixtures. The software evaluated was primarily chosen for its relevance to the CDC wastewater surveillance reporting protocol, which until recently employed a pipeline that incorporates results from four deconvolution methods: Freyja, kallisto, Kraken 2/Bracken, and LCS. We also tested Lollipop, a deconvolution method used by the Swiss SARS-CoV-2 Sequencing Consortium, and three additional methods not used in the C-WAP pipeline: lineagespot, Alcov, and VaQuERo. We found that the commonly used software Freyja outperformed the other CDC pipeline tools in correct identification of lineages present in the control mixtures, and that the VaQuERo method was similarly accurate, with minor differences in the ability of the two methods to avoid false negatives and suppress false positives. Our results also provide insight into the effect of the tiling primer scheme and wastewater RNA extract matrix on viral sequencing and data deconvolution outcomes.

Online dashboards for SARS-CoV-2 wastewater-based epidemiology

Aim: Wastewater-based epidemiology (WBE) is increasingly used to monitor pandemics. In this manuscript, we review methods and limitations of WBE, as well as their online dashboards. Materials & methods: Online dashboards were retrieved using PubMed and search engines, and annotated for timeliness, availability of English version, details on SARS-CoV-2 sublineages, normalization by population and PPMoV load, availability of case/hospitalization count charts and of raw data for export. Results: We retrieved 51 web portals, half of them from Europe. Africa is represented from South Africa only, and only seven portals are available from Asia. Conclusion: WBS provides near-real-time cost-effective monitoring of analytes across space and time in populations. However, tremendous heterogeneity still persists in the SARS-CoV-2 WBE literature.

Tracking SARS-CoV-2 variants of concern in wastewater: an assessment of nine computational tools using simulated genomic data

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.

Detection of SARS-COV-2 variants and their proportions in wastewater samples using next-generation sequencing in Finland

Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) variants may have different characteristics, e.g., in transmission, mortality, and the effectiveness of vaccines, indicating the importance of variant detection at the population level. Wastewater-based surveillance of SARS-CoV-2 RNA fragments has been shown to be an effective way to monitor the COVID-19 pandemic at the population level. Wastewater is a complex sample matrix affected by environmental factors and PCR inhibitors, causing insufficient coverage in sequencing, for example. Subsequently, results where part of the genome does not have sufficient coverage are not uncommon. To identify variants and their proportions in wastewater over time, we utilized next-generation sequencing with the ARTIC Network's primer set and bioinformatics pipeline to evaluate the presence of variants in partial genome data. Based on the wastewater data from November 2021 to February 2022, the Delta variant was dominant until mid-December in Helsinki, Finland's capital, and thereafter in late December 2022 Omicron became the most common variant. At the same time, the Omicron variant of SARS-CoV-2 outcompeted the previous Delta variant in Finland in new COVID-19 cases. The SARS-CoV-2 variant findings from wastewater are in agreement with the variant information obtained from the patient samples when visually comparing trends in the sewerage network area. This indicates that the sequencing of wastewater is an effective way to monitor temporal and spatial trends of SARS-CoV-2 variants at the population level.

Emergence and spread of SARS-CoV-2 variants of concern in Canada: a retrospective analysis from clinical and wastewater data

BACKGROUND

The spread of SARS-CoV-2 has been studied at unprecedented levels worldwide. In jurisdictions where molecular analysis was performed on large scales, the emergence and competition of numerous SARS-CoV-2lineages have been observed in near real-time. Lineage identification, traditionally performed from clinical samples, can also be determined by sampling wastewater from sewersheds serving populations of interest. Variants of concern (VOCs) and SARS-CoV-2 lineages associated with increased transmissibility and/or severity are of particular interest.

METHOD

Here, we consider clinical and wastewater data sources to assess the emergence and spread of VOCs in Canada retrospectively.

RESULTS

We show that, overall, wastewater-based VOC identification provides similar insights to the surveillance based on clinical samples. Based on clinical data, we observed synchrony in VOC introduction as well as similar emergence speeds across most Canadian provinces despite the large geographical size of the country and differences in provincial public health measures.

CONCLUSION

In particular, it took approximately four months for VOC Alpha and Delta to contribute to half of the incidence. In contrast, VOC Omicron achieved the same contribution in less than one month. This study provides significant benchmarks to enhance planning for future VOCs, and to some extent for future pandemics caused by other pathogens, by quantifying the rate of SARS-CoV-2 VOCs invasion in Canada.

Impact of reference design on estimating SARS-CoV-2 lineage abundances from wastewater sequencing data

BACKGROUND

Sequencing of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) RNA from wastewater samples has emerged as a valuable tool for detecting the presence and relative abundances of SARS-CoV-2 variants in a community. By analyzing the viral genetic material present in wastewater, researchers and public health authorities can gain early insights into the spread of virus lineages and emerging mutations. Constructing reference datasets from known SARS-CoV-2 lineages and their mutation profiles has become state-of-the-art for assigning viral lineages and their relative abundances from wastewater sequencing data. However, selecting reference sequences or mutations directly affects the predictive power.

RESULTS

Here, we show the impact of a mutation- and sequence-based reference reconstruction for SARS-CoV-2 abundance estimation. We benchmark 3 datasets: (i) synthetic "spike-in"' mixtures; (ii) German wastewater samples from early 2021, mainly comprising Alpha; and (iii) samples obtained from wastewater at an international airport in Germany from the end of 2021, including first signals of Omicron. The 2 approaches differ in sublineage detection, with the marker mutation-based method, in particular, being challenged by the increasing number of mutations and lineages. However, the estimations of both approaches depend on selecting representative references and optimized parameter settings. By performing parameter escalation experiments, we demonstrate the effects of reference size and alternative allele frequency cutoffs for abundance estimation. We show how different parameter settings can lead to different results for our test datasets and illustrate the effects of virus lineage composition of wastewater samples and references.

CONCLUSIONS

Our study highlights current computational challenges, focusing on the general reference design, which directly impacts abundance allocations. We illustrate advantages and disadvantages that may be relevant for further developments in the wastewater community and in the context of defining robust quality metrics.

Effectiveness of environmental surveillance of SARS-CoV-2 as an early-warning system: Update of a systematic review during the second year of the pandemic

The aim of this updated systematic review was to offer an overview of the effectiveness of environmental surveillance (ES) of SARS-CoV-2 as a potential early-warning system (EWS) for COVID-19 and new variants of concerns (VOCs) during the second year of the pandemic. An updated literature search was conducted to evaluate the added value of ES of SARS-CoV-2 for public health decisions. The search for studies published between June 2021 and July 2022 resulted in 1,588 publications, identifying 331 articles for full-text screening. A total of 151 publications met our inclusion criteria for the assessment of the effectiveness of ES as an EWS and early detection of SARS-CoV-2 variants. We identified a further 30 publications among the grey literature. ES confirms its usefulness as an EWS for detecting new waves of SARS-CoV-2 infection with an average lead time of 1-2 weeks for most of the publication. ES could function as an EWS for new VOCs in areas with no registered cases or limited clinical capacity. Challenges in data harmonization and variant detection require standardized approaches and innovations for improved public health decision-making. ES confirms its potential to support public health decision-making and resource allocation in future outbreaks.

Long-term monitoring of COVID-19 prevalence in raw and treated wastewater in Salvador, the largest capital of the Brazilian Northeast

Wastewater-based epidemiology (WBE) becomes an interesting epidemiological approach to monitoring the prevalence of SARS-CoV-2 broadly and non-invasively. Herein, we employ for the first time WBE, associated or not with the PEG 8000 precipitation method, for the detection of SARS-CoV-2 in samples of raw or treated wastewater from 22 municipal wastewater treatment stations (WWTPs) located in Salvador, the fourth most populous city in Brazil. Our results demonstrate the success of the application of WBE for detecting SARS-CoV-2 in both types of evaluated samples, regardless of the usage of PEG 8000 concentration procedure. Further, an increase in SARS-CoV-2 positivity rate was observed in samples collected in months that presented the highest number of confirmed COVID-19 cases (May/2021, June/2021 and January/2022). While PEG 8000 concentration step was found to significantly increase the positivity rate in treated wastewater samples (p < 0.005), a strong positive correlation (r: 0.84; p < 0.002) between non-concentrated raw wastewater samples with the number of new cases of COVID-19 (April/2021-February/2022) was observed. In general, the present results reinforce the efficiency of WBE approach to monitoring the presence of SARS-CoV-2 in either low- or high-capacity WWTPs. The successful usage of WBE even in raw wastewater samples makes it an interesting low-cost tool for epidemiological surveillance.

Wastewater surveillance uncovers regional diversity and dynamics of SARS-CoV-2 variants across nine states in the USA

Wastewater-based epidemiology (WBE) is a non-invasive and cost-effective approach for monitoring the spread of a pathogen within a community. WBE has been adopted as one of the methods to monitor the spread and population dynamics of the SARS-CoV-2 virus, but significant challenges remain in the bioinformatic analysis of WBE-derived data. Here, we have developed a new distance metric, CoVdist, and an associated analysis tool that facilitates the application of ordination analysis to WBE data and the identification of viral population changes based on nucleotide variants. We applied these new approaches to a large-scale dataset from 18 cities in nine states of the USA using wastewater collected from July 2021 to June 2022. We found that the trends in the shift between the Delta and Omicron SARS-CoV-2 lineages were largely consistent with what was seen in clinical data, but that wastewater analysis offered the added benefit of revealing significant differences in viral population dynamics at the state, city, and even neighborhood scales. We also were able to observe the early spread of variants of concern and the presence of recombinant lineages during the transitions between variants, both of which are challenging to analyze based on clinically-derived viral genomes. The methods outlined here will be beneficial for future applications of WBE to monitor SARS-CoV-2, particularly as clinical monitoring becomes less prevalent. Additionally, these approaches are generalizable, allowing them to be applied for the monitoring and analysis of future viral outbreaks.

Wastewater early warning system for SARS-CoV-2 outbreaks and variants in a Coruña, Spain

Wastewater-based epidemiology has been widely used as a cost-effective method for tracking the COVID-19 pandemic at the community level. Here we describe COVIDBENS, a wastewater surveillance program running from June 2020 to March 2022 in the wastewater treatment plant of Bens in A Coruña (Spain). The main goal of this work was to provide an effective early warning tool based in wastewater epidemiology to help in decision-making at both the social and public health levels. RT-qPCR procedures and Illumina sequencing were used to weekly monitor the viral load and to detect SARS-CoV-2 mutations in wastewater, respectively. In addition, own statistical models were applied to estimate the real number of infected people and the frequency of each emerging variant circulating in the community, which considerable improved the surveillance strategy. Our analysis detected 6 viral load waves in A Coruña with concentrations between 10³ and 10⁶ SARS-CoV-2 RNA copies/L. Our system was able to anticipate community outbreaks during the pandemic with 8-36 days in advance with respect to clinical reports and, to detect the emergence of new SARS-CoV-2 variants in A Coruña such as Alpha (B.1.1.7), Delta (B.1.617.2), and Omicron (B.1.1.529 and BA.2) in wastewater with 42, 30, and 27 days, respectively, before the health system did. Data generated here helped local authorities and health managers to give a faster and more efficient response to the pandemic situation, and also allowed important industrial companies to adapt their production to each situation. The wastewater-based epidemiology program developed in our metropolitan area of A Coruña (Spain) during the SARS-CoV-2 pandemic served as a powerful early warning system combining statistical models with mutations and viral load monitoring in wastewater over time.

VirPool: model-based estimation of SARS-CoV-2 variant proportions in wastewater samples

BACKGROUND

The genomes of SARS-CoV-2 are classified into variants, some of which are monitored as variants of concern (e.g. the Delta variant B.1.617.2 or Omicron variant B.1.1.529). Proportions of these variants circulating in a human population are typically estimated by large-scale sequencing of individual patient samples. Sequencing a mixture of SARS-CoV-2 RNA molecules from wastewater provides a cost-effective alternative, but requires methods for estimating variant proportions in a mixed sample.

RESULTS

We propose a new method based on a probabilistic model of sequencing reads, capturing sequence diversity present within individual variants, as well as sequencing errors. The algorithm is implemented in an open source Python program called VirPool. We evaluate the accuracy of VirPool on several simulated and real sequencing data sets from both Illumina and nanopore sequencing platforms, including wastewater samples from Austria and France monitoring the onset of the Alpha variant.

CONCLUSIONS

VirPool is a versatile tool for wastewater and other mixed-sample analysis that can handle both short- and long-read sequencing data. Our approach does not require pre-selection of characteristic mutations for variant profiles, it is able to use the entire length of reads instead of just the most informative positions, and can also capture haplotype dependencies within a single read.

Viral variant-resolved wastewater surveillance of SARS-CoV-2 at national scale

SARS-CoV-2 surveillance by wastewater-based epidemiology is poised to provide a complementary approach to sequencing individual cases. However, robust quantification of variants and de novo detection of emerging variants remains challenging for existing strategies. We deep sequenced 3,413 wastewater samples representing 94 municipal catchments, covering >59% of the population of Austria, from December 2020 to February 2022. Our system of variant quantification in sewage pipeline designed for robustness (termed VaQuERo) enabled us to deduce the spatiotemporal abundance of predefined variants from complex wastewater samples. These results were validated against epidemiological records of >311,000 individual cases. Furthermore, we describe elevated viral genetic diversity during the Delta variant period, provide a framework to predict emerging variants and measure the reproductive advantage of variants of concern by calculating variant-specific reproduction numbers from wastewater. Together, this study demonstrates the power of national-scale WBE to support public health and promises particular value for countries without extensive individual monitoring.

SARS-CoV-2 wastewater monitoring using a novel PCR-based method rapidly captured the Delta-to-Omicron ΒΑ.1 transition patterns in the absence of conventional surveillance evidence

Conventional SARS-CoV-2 surveillance based on genotyping of clinical samples is characterized by challenges related to the available sequencing capacity, population sampling methodologies, and is time, labor, and resource-demanding. Wastewater-based variant surveillance constitutes a valuable supplementary practice, since it does not require extensive sampling, and provides information on virus prevalence in a timely and cost-effective manner. Consequently, we developed a sensitive real-time RT-PCR-based approach that exclusively amplifies and quantifies SARS-CoV-2 genomic regions carrying the S:Δ69/70 deletion, indicative of the Omicron BA.1 variant, in wastewater. The method was incorporated in the analysis of composite daily samples taken from the main Wastewater Treatment Plant of Thessaloniki, Greece, from 1 December 2021. The applicability of the methodology is dependent on the epidemiological situation. During Omicron BA.1 global emergence, Thessaloniki was experiencing a massive epidemic wave attributed solely to the Delta variant, according to genomic surveillance data. Since Delta does not possess the S:Δ69/70, the emergence of Omicron BA.1 could be monitored via the described methodology. Omicron BA.1 was detected in sewage samples on 19 December 2021 and a rapid increase of its viral load was observed in the following 10-day period, with an estimated early doubling time of 1.86 days. The proportion of the total SARS-CoV-2 load attributed to BA.1 reached 91.09 % on 7 January, revealing a fast Delta-to-Omicron transition pattern. The detection of Omicron BA.1 subclade in wastewater preceded the outburst of reported (presumable) Omicron cases in the city by approximately 7 days. The proposed wastewater surveillance approach based on selective PCR amplification of a genomic region carrying a deletion signature enabled rapid, real-time data acquisition on Omicron BA.1 prevalence and dynamics during the slow remission of the Delta wave. Timely provision of these results to State authorities readily influences the decision-making process for targeted public health interventions, including control measures, awareness, and preparedness.

SARS-CoV-2 wastewater monitoring using a novel PCR-based method rapidly captured the Delta-to-Omicron ΒΑ.1 transition patterns in the absence of conventional surveillance evidence

Conventional SARS-CoV-2 surveillance based on genotyping of clinical samples is characterized by challenges related to the available sequencing capacity, population sampling methodologies, and is time, labor, and resource-demanding. Wastewater-based variant surveillance constitutes a valuable supplementary practice, since it does not require extensive sampling, and provides information on virus prevalence in a timely and cost-effective manner. Consequently, we developed a sensitive real-time RT-PCR-based approach that exclusively amplifies and quantifies SARS-CoV-2 genomic regions carrying the S:Δ69/70 deletion, indicative of the Omicron BA.1 variant, in wastewater. The method was incorporated in the analysis of composite daily samples taken from the main Wastewater Treatment Plant of Thessaloniki, Greece, from 1 December 2021. The applicability of the methodology is dependent on the epidemiological situation. During Omicron BA.1 global emergence, Thessaloniki was experiencing a massive epidemic wave attributed solely to the Delta variant, according to genomic surveillance data. Since Delta does not possess the S:Δ69/70, the emergence of Omicron BA.1 could be monitored via the described methodology. Omicron BA.1 was detected in sewage samples on 19 December 2021 and a rapid increase of its viral load was observed in the following 10-day period, with an estimated early doubling time of 1.86 days. The proportion of the total SARS-CoV-2 load attributed to BA.1 reached 91.09 % on 7 January, revealing a fast Delta-to-Omicron transition pattern. The detection of Omicron BA.1 subclade in wastewater preceded the outburst of reported (presumable) Omicron cases in the city by approximately 7 days. The proposed wastewater surveillance approach based on selective PCR amplification of a genomic region carrying a deletion signature enabled rapid, real-time data acquisition on Omicron BA.1 prevalence and dynamics during the slow remission of the Delta wave. Timely provision of these results to State authorities readily influences the decision-making process for targeted public health interventions, including control measures, awareness, and preparedness.

Wastewater Sequencing-An Innovative Method for Variant Monitoring of SARS-CoV-2 in Populations

The SARS-CoV-2 outbreak has already affected more than 555 million people, and 6.3 million people have died. Due to its high infectivity, it is crucial to track SARS-CoV-2 outbreaks early to prevent the spread of infection. Wastewater monitoring appears to be a powerful and effective tool for managing epidemiological situations. Due to emerging mutations of SARS-CoV-2, there is a need to monitor mutations in order to control the pandemic. Since the sequencing of randomly chosen individuals is time-consuming and expensive, sequencing of wastewater plays an important role in revealing the dynamics of infection in a population. The sampling method used is a crucial factor and significantly impacts the results. Wastewater can be collected as a grab sample or as a 24 h composite sample. Another essential factor is the sample volume, as is the method of transport used. This review discusses different pretreatment procedures and RNA extraction, which may be performed using various methods, such as column-based extraction, TRIzol, or magnetic extraction. Each of the methods has its advantages and disadvantages, which are described accordingly. RT-qPCR is a procedure that confirms the presence of SARS-CoV-2 genes before sequencing. This review provides an overview of currently used methods for preparing wastewater samples, from sampling to sequencing.

SARS-CoV-2 Amino Acid Mutations Detection in Greek Patients Infected in the First Wave of the Pandemic

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), a novel virus that belongs to the Coronoviridae family, emerged in December 2019, causing the COVID-19 pandemic in March 2020. Unlike previous SARS and Middle East respiratory syndrome (MERS) outbreaks, this virus has a higher transmissibility rate, albeit a lower case fatality rate, which results in accumulation of a significant number of mutations and a faster evolution rate. Genomic studies on the mutation rate of the virus, as well as the identification of mutations that prevail and their impact on disease severity, are of great importance for pandemic surveillance and vaccine and drug development. Here, we aim to identify mutations on the SARS-CoV-2 viral genome and their effect on the proteins they are located in, in Greek patients infected in the first wave of the pandemic. To this end, we perform SARS-CoV-2 amplicon-based NGS sequencing on nasopharyngeal swab samples from Greek patients and bioinformatic analysis of the results. Although SARS-CoV-2 is considered genetically stable, we discover a variety of mutations on the viral genome. In detail, 18 mutations are detected in total on 10 SARS-CoV-2 isolates. The mutations are located on ORF1ab, S protein, M protein, ORF3a and ORF7a. Sixteen are also detected in patients from other regions around the world, and two are identified for the first time in the present study. Most of them result in amino acid substitutions. These substitutions are analyzed using computational tools, and the results indicate minor or major impact on the proteins’ structural stability, which could probably affect viral transmissibility and pathogenesis. The correlation of these variations with the viral load levels is examined, and their implication for disease severity and the biology of the virus are discussed.