There are no standardized protocols for quantifying severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in wastewater to date, especially for population normalization. Here, a pipeline was developed, applied, and assessed to quantify SARS-CoV-2 and key variants of concern (VOCs) RNA in wastewater at Saskatoon, Canada. Normalization approaches using recovery ratio and extraction efficiency, wastewater parameters, or population indicators were assessed by comparing to daily numbers of new cases. Viral load was positively correlated with daily new cases reported in the sewershed. Wastewater surveillance (WS) had a lead time of approximately 7 days, which indicated surges in the number of new cases. WS revealed the variant α and δ driving the third and fourth wave, respectively. The adjustment with the recovery ratio and extraction efficiency improved the correlation between viral load and daily new cases. Normalization of viral concentration to concentrations of the artificial sweetener acesulfame K improved the trend of viral load during the Christmas and New Year holidays when populations were dynamic and variable. Acesulfame K performed better than pepper mild mottle virus, creatinine, and ammonia for population normalization. Hence, quality controls to characterize recovery ratios and extraction efficiencies and population normalization with acesulfame are promising for precise WS programs supporting decision-making in public health.

Monitoring the communal incidence of COVID-19 is important for both government and residents of an area to make informed decisions. However, continuous reliance on one means of monitoring might not be accurate because of biases introduced by government policies or behaviours of residents. Wastewater surveillance was employed to monitor concentrations of SARS-CoV-2 RNA in raw influent wastewater from wastewater treatment plants serving three Canadian Prairie cities with different population sizes. Data obtained from wastewater are not directly influenced by government regulations or behaviours of individuals. The means of three weekly samples collected using 24 h composite auto-samplers were determined. Viral loads were determined by RT-qPCR, and whole-genome sequencing was used to charaterize variants of concern (VOC). The dominant VOCs in the three cities were the same but with different proportions of sub-lineages. Sub-lineages of Delta were AY.12, AY.25, AY.27 and AY.93 in 2021, while the major sub-lineage of Omicron was BA.1 in January 2022, and BA.2 subsequently became a trace-level sub-variant then the predominant VOC. When each VOC was first detected varied among cities; However, Saskatoon, with the largest population, was always the first to present new VOCs. Viral loads varied among cities, but there was no direct correlation with population size, possibly because of differences in flow regimes. Population is one of the factors that affects trends in onset and development of local outbreaks during the pandemic. This might be due to demography or the fact that larger populations had greater potential for inter- and intra-country migration. Hence, wastewater surveillance data from larger cities can typically be used to indicate what to expect in smaller communities.