The continuing emergence of SARS-CoV-2 variants of concern and variants of interest emphasizes the need for early detection and epidemiological surveillance of novel variants. We used genomic sequencing of 122 wastewater samples from three locations in Switzerland to monitor the local spread of B.1.1.7 (Alpha), B.1.351 (Beta) and P.1 (Gamma) variants of SARS-CoV-2 at a population level. We devised a bioinformatics method named COJAC (Co-Occurrence adJusted Analysis and Calling) that uses read pairs carrying multiple variant-specific signature mutations as a robust indicator of low-frequency variants. Application of COJAC revealed that a local outbreak of the Alpha variant in two Swiss cities was observable in wastewater up to 13 d before being first reported in clinical samples. We further confirmed the ability of COJAC to detect emerging variants early for the Delta variant by analysing an additional 1,339 wastewater samples. While sequencing data of single wastewater samples provide limited precision for the quantification of relative prevalence of a variant, we show that replicate and close-meshed longitudinal sequencing allow for robust estimation not only of the local prevalence but also of the transmission fitness advantage of any variant. We conclude that genomic sequencing and our computational analysis can provide population-level estimates of prevalence and fitness of emerging variants from wastewater samples earlier and on the basis of substantially fewer samples than from clinical samples. Our framework is being routinely used in large national projects in Switzerland and the UK.

In December 2020, the United Kingdom (UK) reported a SARS-CoV-2 Variant of Concern (VoC) which is now named B.1.1.7. Based on initial data from the UK and later data from other countries, this variant was estimated to have a transmission fitness advantage of around 40–80 % (Volz et al., 2021; Leung et al., 2021; Davies et al., 2021). This study aims to estimate the transmission fitness advantage and the effective reproductive number of B.1.1.7 through time based on data from Switzerland. We generated whole genome sequences from 11.8 % of all confirmed SARS-CoV-2 cases in Switzerland between 14 December 2020 and 11 March 2021. Based on these data, we determine the daily frequency of the B.1.1.7 variant and quantify the variant’s transmission fitness advantage on a national and a regional scale. We estimate B.1.1.7 had a transmission fitness advantage of 43–52 % compared to the other variants circulating in Switzerland during the study period. Further, we estimate B.1.1.7 had a reproductive number above 1 from 01 January 2021 until the end of the study period, compared to below 1 for the other variants. Specifically, we estimate the reproductive number for B.1.1.7 was 1.24 [1.07–1.41] from 01 January until 17 January 2021 and 1.18 [1.06–1.30] from 18 January until 01 March 2021 based on the whole genome sequencing data. From 10 March to 16 March 2021, once B.1.1.7 was dominant, we estimate the reproductive number was 1.14 [1.00–1.26] based on all confirmed cases. For reference, Switzerland applied more non-pharmaceutical interventions to combat SARS-CoV-2 on 18 January 2021 and lifted some measures again on 01 March 2021. The observed increase in B.1.1.7 frequency in Switzerland during the study period is as expected based on observations in the UK. In absolute numbers, B.1.1.7 increased exponentially with an estimated doubling time of around 2–3.5 weeks. To monitor the ongoing spread of B.1.1.7, our plots are available online.

Background: Throughout the COVID-19 pandemic, SARS-CoV-2 genetic variants of concern (VOCs) have repeatedly and independently arisen. VOCs are characterised by increased transmissibility, increased virulence or reduced neutralisation by antibodies obtained from prior infection or vaccination. Tracking the introduction and transmission of VOCs relies on sequencing, typically whole genome sequencing of clinical samples. Wastewater surveillance is increasingly used to track the introduction and spread of SARS-CoV-2 variants through sequencing approaches. Aim: Here, we adapt and apply a rapid, high-throughput method for detection and quantification of the relative frequency of two deletions characteristic of the Alpha, Beta, and Gamma VOCs in wastewater. Methods: We developed drop-off RT-dPCR assays and an associated statistical approach implemented in the R package WWdPCR to analyse temporal dynamics of SARS-CoV-2 signature mutations (spike Δ69–70 and ORF1a Δ3675–3677) in wastewater and quantify transmission fitness advantage of the Alpha VOC. Results: Based on analysis of Zurich wastewater samples, the estimated transmission fitness advantage of SARS-CoV-2 Alpha based on the spike Δ69–70 was 0.34 (95% confidence interval (CI): 0.30–0.39) and based on ORF1a Δ3675–3677 was 0.53 (95% CI: 0.49–0.57), aligning with the transmission fitness advantage of Alpha estimated by clinical sample sequencing in the surrounding canton of 0.49 (95% CI: 0.38–0.61). Conclusion: Digital PCR assays targeting signature mutations in wastewater offer near real-time monitoring of SARS-CoV-2 VOCs and potentially earlier detection and inference on transmission fitness advantage than clinical sequencing.

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.

The large amount and diversity of viral genomic datasets generated by next-generation sequencing technologies poses a set of challenges for computational data analysis workflows, including rigorous quality control, scaling to large sample sizes, and tailored steps for specific applications. Here, we present V-pipe 3.0, a computational pipeline designed for analyzing next-generation sequencing data of short viral genomes. It is developed to enable reproducible, scalable, adaptable, and transparent inference of genetic diversity of viral samples. By presenting 2 large-scale data analysis projects, we demonstrate the effectiveness of V-pipe 3.0 in supporting sustainable viral genomic data science.

Wastewater-based genomic surveillance (WWGS) has proven effective for monitoring SARS-CoV-2 and other viruses within communities. It enables rapid detectio