Abstract
Wastewater monitoring of SARS-CoV-2 allows for early detection and monitoring of COVID-19 burden in communities and can track specific variants of concern. Targeted assays enabled relative proportions of SARS-CoV-2 Omicron and Delta variants to be determined across 30 municipalities covering >75% of the province of Alberta (pop. 4.5M) in Canada, from November 2021 to January 2022. Larger cities like Calgary and Edmonton exhibited a more rapid emergence of Omicron relative to smaller and more remote municipalities. Notable exceptions were Banff, a small international resort town, and Fort McMurray, a more remote northern city with a large fly-in worker population. The integrated wastewater signal revealed that the Omicron variant represented close to 100% of SARS-CoV-2 burden prior to the observed increase in newly diagnosed clinical cases throughout Alberta, which peaked two weeks later. These findings demonstrate that wastewater monitoring offers early and reliable population-level results for establishing the extent and spread of emerging pathogens including SARS-CoV-2 variants.
Introduction
The COVID-19 pandemic has led to rapid scientific progress in wastewater-based surveillance of community infections. Measuring levels of RNA from SARS-CoV-2 in sewage samples began being used as a complementary surveillance tool early in the pandemic, resulting in hundreds of wastewater COVID-19 monitoring groups and online dashboards around the world (Naughton et al. 2021), including in Alberta (https://covid-tracker.chi-csm.ca). This strategy is premised on the fecal shedding of SARS-CoV-2 by infected individuals (Cevik et al. 2021; Yuan et al. 2021) and enabled by modifying RT-qPCR workflows used for diagnosing patients to quantify viral RNA in sewage sampled at wastewater treatment plants (WWTPs) or other nodes within the sewer network (Acosta et al. 2021a, b; Qiu et al. 2022) at regular intervals. Teams in Alberta and elsewhere demonstrated during pandemic waves that wastewater is a leading indicator of COVID-19, with results typically preceding clinical diagnosed cases by 4-6 days (e.g., D’Aoust et al. 2021, Medema et al. 2020; Nemudryi et al. 2020; Randazzo et al. 2020). Sampling, testing and rapidly reporting wastewater virus RNA levels provides early warning of the population-wide disease burden to policy makers, health officials and the public, enabling evidence-based decision making for preparedness and disease control.
On November 24, 2021, South Africa first reported the emergence of a novel SARS-CoV-2 variant associated with rapid community transmission in the Gauteng province (WHO, 2021). By November 26th the World Health Organization had labelled Omicron as a new variant of concern (VOC). Omicron was subsequently rapidly identified in countries around the world, including in Canada where on November 28th cases were detected in inbound international travellers. Alberta’s first case of Omicron from clinical specimen testing was confirmed on November 30. By December and into January this virus had spread rapidly throughout large and smaller communities, prompting re-introduction of public health restrictions.
Wastewater testing can also differentiate changes in disease burden caused by different VOCs in communities (Lee et al. 2021). As soon as viral genomes of VOCs become available within the international scientific community, e.g., via GISAID (Elbe & Buckland-Merrett, 2017), variant-specific PCR primers and probes can be developed and deployed on regularly collected wastewater samples to understand the dynamics of community disease burden caused by VOCs (Peterson et al. 2022). While sequencing viral genomes from wastewater is technically feasible, either via targeted amplicon tiling protocols (Rios et al. 2021; Lin et al. 2021) or shotgun metagenomics (Rothman et al. 2021; Pérez-Cataluña et al. 2022), these comprehensive approaches are significantly more costly and time consuming than targeted RT-qPCR screening of RNA extracted from wastewater that can provide accurate data on VOCs at a fraction of the cost, and in near-real time. In Alberta wastewater is sampled, processed and analyzed in university laboratories in Calgary and Edmonton, and reported to health officials and online to the public two days after sampling. In this study, variant-specific PCR assays were employed to assess the emergence and temporal change in prevalence of the Omicron and Delta variants in Alberta by monitoring wastewater in 30 municipalities ranging from small towns (pop. <10,000) to large cities (pop. >1M) up to three times per week. Wastewater surveillance demonstrated changes in COVID-19 burden associated with emergence of the new Omicron variant for >75% of Alberta’s population of 4.5M between late November 2021 and mid-January 2022.
Methods
Wastewater was collected from municipal WWTPs across the province as 24-h composite samples up to 3 times per week. RNA was isolated from wastewater using either affinity binding columns that purify nucleic acids directly (Whitney et al. 2021) or ultrafiltration followed by RNA extraction (Qiu et al. 2022). These two approaches were used to process 233 and 209 WWTP samples, respectively (Figure 1), with the same method applied to a given sampling site over time throughout the entire study period. Wastewater samples from three geographically disparate WWTPs in Calgary, Fort McMurray and Lethbridge, comprising 11% of all samples in the study, were processed using both methods for comparison and revealed no difference (Mann-Whitney test P = 0.46, 0.39 and 0.59, respectively; Fig. S1).
Wastewater sample processing using affinity columns and ultrafiltration was compared by testing samples from WWTPs in Calgary-1 (n=14), Fort McMurray (n=18) and Lethbridge (n=15) using both methods. The proportion of Omicron obtained following either processing method was compared using Mann-Whitney tests revealing no significant difference in median values between the two sample processing methods. Median and interquartile ranges are indicated as the middle, top, and bottom lines of each box. Ends of the whiskers mark the lowest and highest ratios determined in each sample series.
Proportion of Omicron relative to Delta SARS-CoV-2 variants in community wastewater samples assessed using RT-qPCR assays for specific variants following sample processing using ultrafiltration (A) or affinity columns (B). Lines of best fit plotted with second order smoothing are shown for different WWTPs, including 3 that had samples processed using both ultrafiltration (A) and affinity columns (B) for comparison (see dashed lines for Calgary-1, Fort McMurray and Lethbridge; Fig. S1). Monitoring lasted for 53 days beginning on November 28th (plotted as consecutive days on the x-axes). The grey shaded area on the right side (A, B) shows the 7-day rolling average of new clinical cases reported in Alberta (right y-axis), which increased after the Omicron variant was predominant in municipal wastewater from 30 communities sampled at 21 WWTPs throughout the province (C). Calgary and Edmonton are served by 3 and 2 WWTPs, respectively (A, B), and some individual WWTPs also serve several municipalities (e.g., Edmonton-2 serves 6 others; Red Deer serves 3 others; Calgary’s WWTPs serve 3 others). The timing (in days) of the Omicron-to-Delta ratio passing 10%, 50% and 90% of community COVID-19 burden (D-F) reveals general trends of decreasing population size (bubble diameter and y-axis) and distance from the nearest airport in Calgary, Edmonton or Fort McMurray (bubble shading). Bubble plots only include data from Calgary-1 and Edmonton-1 WWTPs (the largest WWTP from each city), scaled to the population of the corresponding sewershed sub-catchment in those cities.
RNA quantification by RT-qPCR incorporated a newly designed set of assays that selectively amplify the B.1.1.529 Omicron variant or the Delta variant by targeting mutations in amino acids 203 and 204 of the nucleocapsid gene. Omicron and Delta were the only two variants detected in Alberta by clinical screening during the study period (Alberta Health, 2022). Total SARS-CoV-2 levels were quantified with widely used universal assays targeting the N1 & N2 regions of the nucleocapsid gene in the wild type virus (Acosta et al. 2021a, b) and all other VOCs identified to date. Omicron, Delta and total SARS-CoV-2 assays were triplexed together enabling an Omicron-to-Delta ratio to be estimated in each wastewater sample using the Omicron signal (R203K-G204R assay) and the Delta signal (R203M assay). This allowed the emergence and prevalence of Omicron to be tracked at the population level throughout the province.
The daily number of new cases of COVID-19 clinically diagnosed across the province were collected from Data Analytics of Alberta Health Services, and are reported using a 7-day rolling average.
Results
Wastewater separation, identification and quantification of SARS-CoV-2 is intrinsically more complicated than conducting the same PCR strategy on clinical samples (i.e., nasopharyngeal swabs). It is not normally recommended that results be directly compared between different WWTPs, due to intrinsic heterogeneities, e.g., physiochemical differences manifesting different PCR inhibition potential, different proportions of urban, industrial and agricultural inputs, and different flow rates and distances impacting signal degradation (Pecson et al. 2021). While these limitations apply to total SARS-CoV-2 quantification, they are mitigated when determining Omicron-to-Delta ratios within the same multiplex RT-qPCR reaction, since RNA genomes derived from either variant are expected to react similarly to the factors mentioned above.
Figure 1 shows the emergence of the Omicron variant (corresponding to the displacement of the Delta variant) in different Alberta communities. Omicron was first detected in Alberta community wastewater during late November and early December (Fig. 1; Table S1). In Calgary, four consecutive samples collected during December 5-9 revealed the sustained presence of 3-9% Omicron (compared to >90% Delta) among infected individuals contributing to the sewershed in this cosmopolitan city of 1.3 million people. Omicron was first detected in wastewater in the capital city of Edmonton (pop. 1.1 million) on December 10th (15% Omicron; 85% Delta). The rate of increase of Omicron in the international resort town of Banff was higher than in larger cities such as Calgary and Edmonton (Fig. 1) and surpassed 80% in thrice weekly samples taken during December 20-23. By this time, Calgary and Edmonton had just surpassed 50%, and the proportion of Omicron infections was growing in smaller bedroom communities adjacent to these two large urban centres (e.g., Okotoks, High River, Strathmore and especially Airdrie, which are all <70 km away from Calgary; Fig. 1, Fig. S2). Communities that experienced the most delayed emergence of Omicron were smaller and more remote, with Brooks (pop. 14,451; 190 km from Calgary) and Taber (pop. 19,070; 263 km from Calgary) not reaching high proportions of Omicron until December 29th.
Proportion of Omicron in Alberta municipalities plotted as a function of (A) the timing of its emergence, (B) the population served by the municipal WWTP, and (C) the distance of the municipality from Calgary. Timing in panel A corresponds to the midpoint between sampling dates with values below and above 10%. Panels B and C highlight communities that don’t follow the general trends of Omicron emerging earlier in larger municipalities (e.g., Banff has the smallest resident population but increases early), and emerging later in locations farther away Calgary’s major international airport (e.g., Fort McMurray is the farthest away from Calgary but increases early), respectively. Edmonton and Calgary results are limited to the largest WWTPs in the two cities (i.e., Edmonton-1 and Calgary-1 from Figure 1).
Interpretation
The expected general trend demonstrated by this analysis of objective wastewater evidence is that large cities encounter the emergence of a new virus before the smaller centres farther away from the cities, but with notable exceptions. Banff is 127 km west of Calgary and experienced a more rapid onset of Omicron infection than anywhere else in the province despite the resident population in Banff (pop. 13,427) being <1% of Calgary. Banff is an international resort community in Banff National Park – Canada’s busiest national park – which attracts >4 million visitors annually from around the world (Banff, 2022). Early detection of Omicron in Banff may correspond to attracting tourists at the onset of the ski season in November and December. Interestingly the nearby and slightly larger mountain town of Canmore (pop. 27,664) located 105 km east of Calgary (22 km east of Banff, and outside the national park), experienced a much later onset of Omicron infection. This is likely related to Canmore hosting fewer international tourists than Banff, and featuring much less high-density dormitory-style living among the worker population in Banff that supports the tourism industry.
More remote communities located a greater distance away from Alberta’s large international airports exhibited later emergence of the Omicron variant (Fig. 1; Fig. S2). The Calgary International Airport serves 16 million travellers per year with direct flights arriving from 15 countries (YYC, 2022), compared to 8 million travellers and 6 international connections for the Edmonton International Airport (YEG, 2022). Plotting Omicron dynamics in Alberta municipalities as a function of distance from Calgary (Fig. S2) suggests a link to international travel. A notable exception to this trend is Omicron emergence in Fort McMurray. Despite being a remote, relatively small (pop. 79,205) northern city that is farther from Calgary than any other municipality sampled, Fort McMurray exhibited an Omicron emergence comparable to Calgary’s rapid onset. Fort McMurray has one of Canada’s busiest airports to accommodate shift workers commuting from across Canada to work in the oil sands industry (YMM, 2022). This high level of contact with other parts of the country is likely to facilitate rapid introduction of an emerging virus such as the Omicron variant. Workers travelling to Fort McMurray from other provinces, or Alberta’s major urban centres, likely contributed to accelerated Omicron emergence relative to other smaller and/or remotely situated Alberta municipalities.
Wastewater results also demonstrate that the emergence of Omicron was the driver of clinical cases increasing in December and January during Alberta’s 5th wave (grey shaded area in Fig. 1A, B). During this time COVID-19 public health surveillance shifted to focus PCR testing on patients at risk for severe illness and eligible for early treatment, patients presenting to emergency wards with more serious illness, and essential workers (AHS, 2022). This resulted in PCR testing dramatically underestimating total disease burden in the population as a whole relative to earlier waves. Clinical cases that were reported still show a steep increase after the emergence and propagation of Omicron revealed by wastewater testing. These dynamics mirror Alberta’s shift from Delta to Omicron infections confirmed by screening sub-sets of clinical samples using PCR assays for VOCs and associated viral genome sequencing, which revealed Omicron to be >50% by December 16th and >95% by December 28th (Alberta Health, 2022). This demonstrates that wastewater surveillance reliably provides important information needed and acted upon by public health officials.
VOC information derived from viral genome sequencing of clinical samples is typically non-random, e.g, being intentionally biased towards clinical cases of interest (e.g., outbreaks; hospitalizations) or incoming international travellers (e.g., Williams et al. 2021). Similarly, clinical PCR testing is susceptible to changes in testing policies, capacity limitations or individuals not getting tested, e.g., by personal choice or when infections are asymptomatic (Green et al. 2021). Wastewater testing on the other hand offers an unbiased representation of disease prevalence, capturing all individuals and groups contributing to the sewershed for a tiny fraction of the cost of clinical testing on a per-capita basis. In large cities like Calgary and Edmonton with over one million residents (Fig. 1D-F; Fig. S1) monitoring wastewater for COVID-19 community burden costs only a few cents per person per year (based on three times weekly testing in Alberta) and can provide objective information about community infection to public health authorities, policy makers and the public in near real time. With COVID-19 clinical testing strategies and resources become more targeted, wastewater testing offers an important, early, objective population-based metric of disease burden surveillance that can be easily adapted for emerging pathogens throughout large jurisdictions like Alberta.
Data Availability
All data produced in the present study are available upon reasonable request to the authors
Chronology of detection of the Omicron variant in Alberta community wastewater
Acknowledgements
This work was funded by the Government of Alberta. The authors wish to thank Dr. Norma Ruecker and Dr. Rhonda Clark for help with sampling and logistics. This work would not have been possible without collaboration from provincial and municipal leaders and WWTP operators in Calgary, Edmonton, Fort McMurray, Grande Prairie, Cold Lake, Edson, Lacombe, Red Deer, Banff, Canmore, Drumheller, Strathmore, Okotoks, High River, Brooks, Medicine Hat, Taber and Lethbridge.
