Wastewater monitoring outperforms case numbers as a tool to track COVID-19 incidence dynamics when test positivity rates are high

Wastewater-based epidemiology (WBE) has been shown to coincide with, or anticipate, confirmed COVID-19 case numbers. During periods with high test positivity rates, however, case numbers may be underreported, whereas wastewater does not suffer from this limitation. Here we investigated how the dynamics of new COVID-19 infections estimated based on wastewater monitoring or confirmed cases compare to true COVID-19 incidence dynamics. We focused on the first pandemic wave in Switzerland (February to April, 2020), when test positivity ranged up to 26%. SARS-CoV-2 RNA loads were determined 2-4 times per week in three Swiss wastewater treatment plants (Lugano, Lausanne and Zurich). Wastewater and case data were combined with a shedding load distribution and an infection-to-case confirmation delay distribution, respectively, to estimate incidence dynamics. Finally, the estimates were compared to reference incidence dynamics determined by a validated compartmental model. Incidence dynamics estimated based on wastewater data were found to better track the timing and shape of the reference infection peak compared to estimates based on confirmed cases. In contrast, case confirmations provided a better estimate of the subsequent decline in infections. Under a regime of high-test positivity rates, WBE thus provides critical information that is complementary to clinical data to monitor the pandemic trajectory.

Detection of SARS-CoV-2 variants in Switzerland by genomic analysis of wastewater samples

The emergence of SARS-CoV-2 mutants with altered transmissibility, virulence, or immunogenicity emphasizes the need for early detection and epidemiological surveillance of genomic variants. Wastewater samples provide an opportunity to assess circulating viral lineages in the community. We performed genomic sequencing of 122 wastewater samples from three locations in Switzerland to analyze the B.1.1.7, B.1.351, and P.1 variants of SARS-CoV-2 on a population level. We called variant-specific signature mutations and monitored variant prevalence in the local population over time. To enable early detection of emerging variants, we developed a bioinformatics tool that uses read pairs carrying multiple signature mutations as a robust indicator of low-frequency variants. We further devised a statistical approach to estimate the transmission fitness advantage, a key epidemiological parameter indicating the speed at which a variant spreads through the population, and compared the wastewater-based findings to those derived from clinical samples. We found that the local outbreak of the B.1.1.7 variant in two Swiss cities was observable in wastewater up to 8 days before its first detection in clinical samples. We detected a high prevalence of the B.1.1.7 variant in an alpine ski resort popular among British tourists in December 2020, a time when the variant was still very rare in Switzerland. We found no evidence of local spread of the B.1.351 and P.1 variants at the monitored locations until the end of the study (mid February) which is consistent with clinical samples. Estimation of local variant prevalence performs equally well or better for wastewater samples as for a much larger number of clinical samples. We found that the transmission fitness advantage of B.1.1.7, i.e. the relative change of its reproductive number, can be estimated earlier and based on substantially fewer wastewater samples as compared to using clinical samples. Our results show that genomic sequencing of wastewater samples can detect, monitor, and evaluate genetic variants of SARS-CoV-2 on a population level. Our methodology provides a blueprint for rapid, unbiased, and cost-efficient genomic surveillance of SARS-CoV-2 variants.