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.

An outbreak of SARS-CoV-2 Alpha variant (Pango lineage B.1.1.7) was detected at a primary school (School X) in Lansingerland, the Netherlands, in December 2020. The outbreak was studied retrospectively, and population-based screening was used to assess the extent of virus circulation and decelerate transmission. Cases were SARS-CoV-2 laboratory confirmed and were residents of Lansingerland (November 16th 2020 until February 22th 2021), or had an epidemiological link with School X or neighbouring schools. The SARS-CoV-2 variant was determined using variant PCR or whole genome sequencing. A questionnaire primarily assessed clinical symptoms. A total of 77 Alpha variant cases were found with an epidemiological link to School X, 16 Alpha variant cases linked to the neighbouring schools, and 146 Alpha variant cases among residents of Lansingerland without a link to the schools. The mean number of self-reported symptoms was not significantly different among Alpha variant infected individuals compared to non-Alpha infected individuals. The secondary attack rate (SAR) among Alpha variant exposed individuals in households was 52% higher compared to non-Alpha variant exposed individuals (p = 0.010), with the mean household age, and mean number of children and adults per household as confounders. Sequence analysis of 60 Alpha variant sequences obtained from cases confirmed virus transmission between School X and neighbouring schools, and showed that multiple introductions of the Alpha variant had already taken place in Lansingerland at the time of the study. The alpha variant caused a large outbreak at both locations of School X, and subsequently spread to neighbouring schools, and households. Population-based screening (together with other public health measures) nearly stopped transmission of the outbreak strain, but did not prevent variant replacement in the Lansingerland municipality.