Over the course of the Corona Virus Disease-19 (COVID-19) pandemic in 2020–2022, monitoring of the severe acute respiratory syndrome coronavirus 2 ribonucleic acid (SARS-CoV-2 RNA) in wastewater has rapidly evolved into a supplementary surveillance instrument for public health. Short term trends (2 weeks) are used as a basis for policy and decision making on measures for dealing with the pandemic. Normalisation is required to account for the dilution rate of the domestic wastewater that can strongly vary due to time- and location-dependent sewer inflow of runoff, industrial discharges and extraneous waters. The standard approach in sewage surveillance is normalisation using flow measurements, although flow based normalisation is not effective in case the wastewater volume sampled does not match the wastewater volume produced. In this paper, two alternative normalisation methods, using electrical conductivity and crAssphage have been studied and compared with the standard approach using flow measurements. For this, a total of 1116 24-h flow-proportional samples have been collected between September 2020 and August 2021 at nine monitoring locations. In addition, 221 stool samples have been analysed to determine the daily crAssphage load per person. Results show that, although crAssphage shedding rates per person vary greatly, on a population-level crAssphage loads per person per day were constant over time and similar for all catchments. Consequently, crAssphage can be used as a quantitative biomarker for populations above 5595 persons. Electrical conductivity is particularly suitable to determine dilution rates relative to dry weather flow concentrations. The overall conclusion is that flow normalisation is necessary to reliably determine short-term trends in virus circulation, and can be enhanced using crAssphage and/or electrical conductivity measurement as a quality check.

Monitoring of SARS-CoV-2 in wastewater (WW) is a promising tool for epidemiological surveillance, correlating not only viral RNA levels with the infection dynamics within the population, but also to viral diversity. However, the complex mixture of viral lineages in WW samples makes tracking of specific variants or lineages circulating in the population a challenging task. We sequenced sewage samples of 9 WW-catchment areas within the city of Rotterdam, used specific signature mutations from individual SARS-CoV-2 lineages to estimate their relative abundances in WW and compared them against those observed in clinical genomic surveillance of infected individuals between September 2020 and December 2021. We showed that especially for dominant lineages, the median of the frequencies of signature mutations coincides with the occurrence of those lineages in Rotterdam's clinical genomic surveillance. This, along with digital droplet RT-PCR targeting signature mutations of specific variants of concern (VOCs), showed that several VOCs emerged, became dominant and were replaced by the next VOC in Rotterdam at different time points during the study. In addition, single nucleotide variant (SNV) analysis provided evidence that spatio-temporal clusters can also be discerned from WW samples. We were able to detect specific SNVs in sewage, including one resulting in the Q183H amino acid change in the Spike gene, that was not captured by clinical genomic surveillance. Our results highlight the potential use of WW samples for genomic surveillance, increasing the set of epidemiological tools to monitor SARS-CoV-2 diversity.