Abstract
Wastewater-based surveillance (WBS) has become an effective tool around the globe for indirect monitoring of COVID-19 in communities. Quantities of viral fragments of SARS-CoV-2 in wastewater are related to numbers of clinical cases of COVID-19 reported within the corresponding sewershed. Variants of Concern (VOCs) have been detected in wastewater by use of reverse transcription quantitative polymerase chain reaction (RT-qPCR) or sequencing. A multiplex RT-qPCR assay to detect and estimate the prevalence of multiple VOCs, including Omicron/Alpha, Beta, Gamma, and Delta, in wastewater RNA extracts was developed and validated. The probe-based multiplex assay, named “N200” focuses on amino acids 199-202, a region of the N gene that contains several mutations that are associated with variants of SARS- CoV-2 within a single amplicon. Each of the probes in the N200 assay are specific to the targeted mutations and worked equally well in single- and multi-plex modes. To estimate prevalence of each VOC, the abundance of the targeted mutation was compared with a non- mutated region within the same amplified region. The N200 assay was applied to monitor frequencies of VOCs in wastewater extracts from six sewersheds in Ontario, Canada collected between December 1, 2021, and January 4, 2022. Using the N200 assay, the replacement of the Delta variant along with the introduction and rapid dominance of the Omicron variant were monitored in near real-time, as they occurred nearly simultaneously at all six locations. The N200 assay is robust and efficient for wastewater surveillance can be adopted into VOC monitoring programs or replace more laborious assays currently being used to monitor SARS- CoV-2 and its VOCs.
Introduction
Communities around the world have quickly adopted wastewater-based surveillance (WBS) of SARS-CoV-2 viral fragments to detect and infer prevalence of COVID-19 disease burden in communities (D’Aoust et al., 2021; Kitajima et al., 2020; Medema et al., 2020). Because it does not require testing of individual sensitivity and includes all individuals, symptomatic or asymptomatic, that contribute to the sewershed, WBS has several advantages over more costly and invasive clinical testing. WBS utilizes reverse transcription quantitative real-time polymerase chain reaction (RT-qPCR) based techniques targeting various regions of the SARS-CoV-2 genome, including the popular US-Centers for Disease Control (CDC) N1, N2 (nucleocapsid) and E (envelope), to test for the presence of SARS-CoV-2 in wastewater (Hamouda et al., 2021; Kitajima et al., 2020; Kumblathan et al., 2021).
Continuous emergence of Variants of Concern (VOCs) and Variants of Interest (VOIs) has presented additional challenges to the surveillance of SARS-CoV-2 and has been an area where WBS can provide more timely information quicker than clinical testing. In addition to clinical tracking of abundance of SARS-CoV-2 infections, RT-qPCR methods have been used to detect and monitor prevalence of emerging VOCs in wastewater (Carcereny et al., 2021; Graber et al., 2021; Lee et al., 2021; Peterson et al., 2022; Wurtzer et al., 2021). This information provides an integrated estimation of the relative abundance of variants within populations. There have been several examples of RT-qPCR assays for identification and quantifications of VOCs wastewater. The mutations N501Y and Δ69-70 del in the spike protein or the presence of mutation D3L in the N gene have been targeted by RT-qPCR assays to monitor for the presence of VOCs in wastewaters across Canada (Peterson et al., 2022). Using an allele specific RT-qPCR, with an artificial mismatch, three mutations, HV69/70del, Y144del and A570D, were targeted to discriminate and quantify the Alpha from VOC relative to the wild type in wastewater of 19 communities in the USA (Lee et al, 2021). Another assay targets mutations Δ69-70 in the S gene to detect the Alpha and Beta variants respectively (Yaniv et al., 2021). An assay based on mutation D3L of the gene, was applied for detection and quantification of relative proportion of the Alpha VOC (Graber et al., 2021). Those authors were able to effectively trace the increase of the proportion of Alpha variant in wastewater of Ottawa, Canada, in almost real time. The trends in incidence of Alpha were similar to that of the number of clinical cases in the community. These various studies demonstrate that WBS can rapidly provide information on incidence and prevalence of SARS-CoV-2 to authorities with sensitivity and lineage specificity. RT-qPCR assays are specific and can be deployed rapidly across existing WBS networks. The challenges with using RT-qPCR assays in surveillance of VOCs are the amount of time needed to develop and validate new assays each time a new VOC emerges as well as increasing the number assays being conducted to continually monitor for various VOCs and VOIs.
The objective of this study was to develop a multiplex RT-qPCR assay that could be used to simultaneously screen for several VOCs, as well as a marker of the total SARS-CoV-2 present in the sample so that an estimation of the relative proportion of VOCs can be reliably determined. The N gene region was selected for several reasons, including that this region is targeted by the N1 and N2 assays, developed by the US CDC, which are commonly used for monitoring wastewater (Lu et al., 2020). There is evidence to suggest that assays developed for the N-region of the SARS-CoV-2 genome have greater sensitivity for measuring various VOCs of SARS-CoV-2 in wastewater than assays targeting than does the S-(Yaniv et al., 2021) or the envelope regions (Pérez-Cataluña et al., 2021). The N-gene was selected as the target region as the target of the assay since it is also rich in mutations unique to VOCs (Kiryanov et al., 2022). Specifically, the 121-basepair region in the N-gene open reading frame (ORF) (nucleotides 28,837 to 28,958) that includes single or multiple nucleotide variants for each VOC was selected (Table 1). According to GISAID data summarized by nextrain.org (Hadfield et al., 2018), the nucleotide 28880- 28882 (AA N:203) have one of the highest entropy levels in the SARS-CoV-2 genome, at 0.972, with nucleotides 28883-28885 (AA N:204) also having high entropy level of 0.711 (nextstrain.org; accessed April 1, 2022). The non-synonymous functional polymorphisms present in this region include R203M, present in Delta, which increases the spread of the virus (Syed et al., 2021); R203K and G204R, (both present in Alpha, Gamma Lambda, and Omicron, which are related to increased infectivity, fitness and virulence (Johnson et al., 2021; Lee et al., 2021; Wu et al., 2021); and T205I, present in Beta and Mu variants that were predicted to have reduced antigenicity and greater affinity for HLA-I (Antonio et al., 2021). The presence of multiple mutations in a single region allowed for the design of a single amplicon, multiple probe assay, that can distinguish among several VOCs. A Universal probe, that targets an area of the amplicon without mutations allows for the quantification of the total SARS-CoV-2 signal, regardless of which variants are present in wastewater. This allows for the measurement of the total SARS-CoV-2 signal as well as providing estimations of the relative abundance of individual VOCs.
Amino acid mutations and corresponding nucleotide mutations present in the targeted region of the SARS-CoV-2 nucleocapisd genome (AA 199-214).
Prior to being incorporated into a wastewater surveillance program, the N200 assay was validated using several synthetic standards as well as samples of wastewaters. As an example of the utility of this assay, between November and January 2021, samples from six wastewater treatment facilities in Ontario, Canada, were examined with the N200 assay. During this period the N200 assay was utilized to test for the presence and increased frequency of the R203K/G204R mutation (present in the Omicron VOC), while comparing directly to the decrease in frequency of the R203M mutation (present in the Delta VOC), relative to a VOC-agnostic Universal marker, all within a single RT-qPCR reaction.
Materials and Methods
Target selection
To select an appropriate RT-qPCR target for testing wastewater, mutations that were indicative of each VOC, encompassing amino acids N:202-205 were identified. Mutations selected for each VOC and their known prevalence are presented in Table 1. The template sequence used for designs were retrieved from accession numbers provided by TWIST Bioscience (South San Francisco, CA, USA) for synthetic controls 14 (Alpha), 16 (Beta), 17 (Gamma) and 23 (Delta) and sequences were aligned using MAFFT (Katoh et al., 2002).
Assay design
Candidate forward and reverse primers that encompassed an amplicon 121 base pairs long containing all the desired mutations were designed with a target predicted Tm around 61°C and 40-60% GC content. Candidate primers were screened for hairpins, homodimers and heterodimers in silico using OligoAnalyzer by integrated DNA technologies (IDT; RRID:SCR_001363) with the “qPCR” parameter selected. Candidate primers that passed were then screened against NCBI’s non-redundant nucleotide database using BLAST (Altschul et al., 1990) and Primer-BLAST (Ye et al., 2012) to maximize mismatches with as much non-target genetic material as possible (human, microbes, food, etc.).
Five probes were designed to be selective for each of the mutations (Table 1, 2), as well as a universal probe for the “Universal” probe complementary to a nearby, highly conserved (bp 28891-28905, mutation rates <0.87% according to covidcg.org) region of the amplicon to detect total SARS-CoV-2 (Table 2, Figure 1). Probes were designed with the mutated base(s) at their center using PrimerExpress v3.0.1 with a GC content of 40-60%. BHQplus technology (LCG, Biosearch Technologies, Inc., California, USA) was used to increase the Tm of each probe, to have a predicted Tm of 69°C, which enabled the design of short probes approximately 15 nucleotides long. Probes were then tested for hairpins, homodimers, and heterodimers. Candidate probes were also screened for % identity to non-target sequences using NCBI’s BLAST tool (Agarwala et al., 2016).
Forward (F) and Reverse (R) primer and probe designs for the N200 VOC assay. The reporter and quencher moieties, and final working concentrations are included for each oligonucleotide.
N200 assay amplicon with the location of the N200 forward and reverse primers and all probes are shown on a portion of the N gene of the SARS-CoV-2 genome. The sequence, nucleotide numbers, and ammino acid displayed are based on the SARS-CoV-2 reference genome (Wuhan; NC_045512.2; figure generated using Benchling.com).
Assay Validation
Initial tests of the N200 assay were performed on synthetic templates. TWIST synthetic RNA controls 2, 14, 16, 17, and 23 (TWIST Bioscience) were employed to generate standard curves (sequences representative in lieu of Wuhan-Hu-1 (Genbank MN908957.3), as well as the VOCs, Alpha, Beta, Gamma, and Delta VOCs, respectively). The N200 standard curves was assessed for linearity and efficiency by use of each of the probes in singleplex. After validation with synthetic templates assays were tested with wastewater RNA extracts to ensure successful amplification. Temperature gradients and primer/probe concentration gradients were performed to determine the optimal assay conditions in singleplex.
Cross-reactivity of probes was assessed by use of synthetic RNA controls. Each probe, in singleplex, was tested for reactivity with TWIST RNA controls 2, 14, 16, 17 and 23 at a concentration of 500 copies per reaction. Specificity was also determined in a triplex assay iteration (Universal, R203M, and R203K probes) by use of digital PCR (dPCR) at multiple concentrations of TWIST RNA controls (see dPCR section for additional details).
The N200 assay was validated in multiplex by comparing the efficiency, R2 and y- intercept of standard curves, as well as the estimated concentration and proportion of variants tested in wastewater relative to singleplex. Sensitivity of the N200 assay with the Universal, R203K, and R203M probes were determined by use of a 12-point standard curve (TWIST synthetic controls) with 15-replicates each. The limit of detection (LOD) was defined as the standard concentration at which >95% of the replicates were detected. This was performed by two independent labs (at the University of Waterloo and the Canadian National Microbiology Laboratory – Winnipeg) using the same suppliers of primers and probes and standards.
Wastewater sample collection, concentration, and extraction
Wastewater from the Regions of Peel (wastewater treatment plant influents GE Booth and Clarkson), York (access points Humber Air Management Facility (AMF) and Warden) and Waterloo (wastewater treatment plant influents Kitchener and Waterloo) Ontario, Canada. At each site wastewater was collected by plant operators using 24-hour cooled composite samplers (combined sample of three grabs at Warden) 3-5 times a week between November 28th, 2021, and January 4th, 2022. Samples were aliquoted into 250 mL HPDE bottles (Systems Plus, Baden, Canada) and transported on ice to the University of Waterloo.
To concentrate RNA, 250 mL samples were well-mixed by inversion, 40 mL aliquots poured into 50 mL conical tubes with 4 g PEG-8000 and 0.9 g NaCl, then mixed on an orbital shaker at low speed at 4°C for 2 h before being stored at 4°C overnight (Wu et al., 2020). Samples were then centrifuged at 12,000 x g at 4°C for 90 minutes with no brake. The supernatant was discarded, and the remaining sample was centrifuged again at 12,000 x g at 4°C for 5 min with no brake. The remaining supernatant was decanted and pipetted out and the wet weight of the resulting pellet recorded. Up to 250 mg of the pellet was used for extraction of RNA by use of the RNeasy PowerMicrobiome Kit (Qiagen, Germantown, MD) as per manufacturers instructions and included the addition of 100 µL of TRIzol (Thermo Fisher, Mississauga, Canada) to the pellet before bead-beating. Total RNA was eluted in 100 µL of RNase-free water.
RT-qPCR
One-step RT-qPCR was performed using TaqPath 1-Step Master Mix (Thermo Fisher, Mississauga, Canada). The primer and probe sequences and final working concentrations are presented in Table 2. RT-qPCR reactions were run in triplicate using 5 µL RNA template and a final reaction volume of 20 µL. Cycling was performed on an OPUS or CFX 96 Touch qPCR thermocycler (BioRad, Hercules, CA) as follows: RT at 50°C for 15 min, 95°C for 2 min, 45 cycles of 95°C for 3 sec followed by 57°C for 30 sec. All qPCR plates were run with controls including no template controls (NTCs), a positive control (wastewater sample with previously determined amount of Delta variant) and a standard curve based on a homologous synthetic RNA (i.e., EDX, TWIST).
Digital PCR (dPCR)
For precise quantification, commercial standards were analyzed via singleplex one-step RT-dPCR (QIAcuity, Qiagen, Hilden, Germany). QIAquity One-Step Viral RT-PCR Kit (Qiagen) together with CDC_N1 probe and primers (Lu et al., 2020) were used to determine absolute copy number (cp) of RNA in the following SARS-CoV-2 genomic RNA standards: TWIST Bioscience Controls 2, 14, 16, 17, 23. Testing for specificity of probes was also conducted using the N200 assay in triplex (Universal, R203M, and R203K probes) with two concentrations of TWIST Bioscience Controls 14, and 23. Reactions consisted of 10 µL of RNA template, forward and reverse primers (working concentration of 500 nM each), probe (working concentration of 125 nM for N1, or 100, 400 and 400 nM for the Universal, R203M and R203K probes respectively), 10 µL of supermix, 0.4 µL of reverse transcriptase, and 14.6 µL of PCR-grade water to make up the balance of the 40 µL reaction. All reactions were performed in triplicate using 24 well x 26,000 channel plates. Cycling was performed on the QIAcuity as follows: RT at 50°C for 30 min, 95°C for 2 min, with 40 cycles at 95°C for 3 sec then annealing at 55°C or 57°C (for N1 and N200 assays respectively) for 30 sec. The N200 assay was also tested on RT-dPCR in triplex (Universal, R203M, and R203K probes) using the TWIST Bioscience Controls 23 and 48 at multiple concentrations to determine the specificity of the assay.
Data Analysis
Samples were quantified (copies/well) by use of a 7-point standard curve with the relevant TWIST standard. Samples were then corrected for elution volume and volume of wastewater extracted and reported as copies/mL wastewater (Quantity). The proportion of total Quantity associated with each VOC in a sample was determined by use of the equation mutation specific for the VOC (Equation 1).
Results
N200 assay conditions and performance
The N200 assay performed well with all probes in singleplex. When tested with a serial dilution of TWIST standards, the singleplex standard curves were linear, had an efficiency near 100%, and the slopes and y-intercepts were comparable between each of the probes (Table 3). The N200 assay performed similarly when tested in singleplex, duplex and triplex (Table 3). All probes targeting mutations were found to be highly specific with little to no cross-reactivity for non-target templates using qPCR (in singleplex; Table 4). Using RT-qPCR, this was demonstrated by the lack of amplification on all VOC TWIST templates without the target mutation (Table 4). Some cross-reactivity was observed between the T205I probe and the Wuhan variant control (TWIST Control 2; Table 4), but this accounted for less than 5% of the signal observed on the homologous template. Positive amplification with similar sensitivity was seen on of all TWIST templates tested was observed with the Universal probe (Table 4). High target specificity of each probe was confirmed by RT-dPCR, where in the presence of a high concentration of non- target template, almost no amplification was detected (Table 5). Copy numbers reported by the mutation probes (R203M, R203K) were comparable to the copy numbers reported by the Universal probe in the RT-dPCR reactions at multiple concentrations, which further suggests that each probe has similar sensitivities and efficiencies (Table 5). The LOD (95% detection) was determined to be 4-6 copies per reaction for the Universal, R203K, and R203M probes respectively when tested in two separate laboratories.
Performance metrics of the N200 VOC assay for each of the probes in singleplex, duplex or triplex.
Cycle Threshold (CT) value (mean) of N200 assay performed in singleplex with probes targeting a mutation (R203K; T205I; S202S; R203M) or a common region (Universal) of the SARS-CoV-2 genome. Individual probes were performed on each of five synthetic templates synthesized by TWIST Bioscience (San Francisco, USA), controls for the original Wuhan strain and for each VOC. ND = not detected
Gene copies/µL of RNA (mean ± SD) detected by the Universal, R203M, or R203K probes in the N200 assay (triplexed) when tested against TWIST controls for variants of SARS-CoV-2 using RT-dPCR.
Performance of N200 assay with samples of wastewater
The triplex N200 assay was applied to wastewater samples throughout December 2021 and January 2022 to monitor for community transmission of the Omicron VOC, which would be signified by decrease in the proportion of the R203M (Delta) mutation and an increase in the R203K mutation (Omicron). The total SARS-CoV-2 (Universal) signal in samples during this period started out low, but rapidly increased over the sampling period (Figure 2). During this time, we observed a rapid transition in relative abundance of mutations from being dominated by R203M (Delta) to being dominated by R203K (Omicron) at all monitoring locations (Figure 3).
SARS-CoV-2 signal measured in municipal wastewater using the N200 assay with the Universal probe in a multiplexed rection (Log 2 scale – copies/mL). Samples from the Peel (GE Booth and Clarkson), York (Humber AMF and Warden) and Waterloo (Waterloo and Kitchener) Regions were collected during the period encompassing both the introduction of Omicron into these communities and the concomitant reduction of the endemic Delta VOC.
Application of the N200 triplex assay to monitor the relative abundance (%) of the mutations R203M (black circles- presumed Delta) or R203K/G204R (red squares – presumed Omicron) in municipal wastewater effluent. Effluent was collected from the Region of Peel (A, B), York Region (C, D) and Waterloo Region (E, F) using composite (A-B, D-F) and grab (C) sampling. Lines drawn through points are a 5-point sigmoidal curve calculated using GraphPad PRISM (v 9.3.1).
Discussion
This study presents an RT-qPCR approach for detection of SARS-CoV-2 and VOCs Omicron/Alpha, Beta, Gamma and Delta in a multiplexed assay with a single amplicon. The specificity and sensitivity of the N200 assay allow for simultaneous detection and quantification of mutations associated with VOCs. These individual assays were multiplexed and validated in duplex and triplex reactions, which allowed for rapid, cost-effective detection and quantitation of mutations associated with VOCs. Relative amounts of these mutations are then compared to a non-mutated region of the same amplicon to estimate the prevalence of current VOCs in a wastewater sample. This assay was able to track the initial presence and then the transition from the R203M mutation (i.e., Delta) as it was replaced by the R203K/G204R mutation (i.e., Omicron/Alpha) in multiple communities in south-central Ontario in December 2021. This is representative of the transition from the Delta variant to the Omicron variant observed in clinical cases at this time in these communities.
In addition to R203M (Delta) and R203K/G204R (Omicron/Alpha) targets, probes targeting the T205I (Beta) mutation and the silent mutation in the amino acid S202 (Gamma), allow for surveillance of known VOCs using a single assay. Having multiplexed, variant specific probes that are compared to a Universal reference within the same amplicon is beneficial for improved consistency and accuracy of measuring relative prevalence of variants. As the gene targets are all in the same 121 bp amplicon, there is minimal concern about the dynamics of different regions of the SARS-CoV-2 genome in either the wastewater system itself, or in the extraction process. Because the reactions occur together in a multiplexed assay there is reduced variability and replication compared to running assays separately. This approach saves resources and time, as fewer reactions are run per sample. Critically, this assay also saves sample volume, which is often limited, allowing for additional assays to be run on the same extract or reducing the need for multiple extractions to be performed. Ultimately it makes the process more efficient, reduces turnaround time and improves quality.
The combined use of Universal and VOC-specific probes in the N200 assay to estimate prevalence of VOCs in wastewater is advantageous for several reasons. First, since when multiplexed, the mutation and comparator are being assessed within the same amplicon as well the same reaction, this method provides increased precision of VOCs prevalence estimates. Additional confidence in prevalence of VOCs is provided by the fact that the mutation and Universal signals can be quantified using the same standard. In comparison, other VOC qPCR assays estimate prevalence by comparing a mutant allele to a genetic marker of SARS-COV-2 in another amplicon, /genetic locus such as the CDC_N1 (Vogels et al., 2021; Wolfe et al., 2021), or by comparing the frequencies of the mutant and “wild-type” (e.g., ancestral/endemic allele) at the same locus (Chan et al., 2022; Graber et al., 2021; Lee et al., 2021; Peterson et al., 2022). While these assays can produce accurate data (from both analytical and clinical perspectives), there are additional controls needed to account for variations in efficiencies of variation (inter- assay variability and matching standard quantities), that are eliminated by the N200 assay. Estimations of VOC prevalence are important during wastewater surveillance SARS-CoV-2 in that they can provide useful information to public health units in a timely manner.
During analyses of wastewaters through late 2021 and early 2022, a rapid increase in the SARS-CoV-2 signal (Universal probe) as well as a quick change over of the dominant variant from Delta (R203M) to Omicron (R203K/G204R) was observed at all sites monitored. The assay first detected the R203K mutation in a sample collected from GE Booth (Peel Region) on December 3rd, 2021. While it is possible that this early detection of the R203K/G204R mutation could have been a variant other than Omicron, within a week of this detection, the R203K/G204R mutation was consistently detected in the wastewater at this site while R203M decreased in proportion. The use of two mutation targets provided confidence in the data as it was being collected and reported because both a decrease in the proportion of the R203M mutation (Delta) and an increase in the proportion of the R203K/G204R mutations (Omicron) occurred simultaneously. This data was reported to public health partners within a week of samples being collected and this was immediately disseminated to the public by the Public Health Unit as part of a weekly Epi-Report (https://data.peelregion.ca/documents/profile-of- covid-19-cases/explore). Rapid emergence of Omicron presented a major challenge for clinical testing and sequencing capacity. This was also confounded by the holiday season in Ontario. The wastewater data in general, quickly became important to Public Health Units as the clinical testing became less reliable because of changes to testing eligibility in Ontario. They were able to confirm the rapid replacement of Delta with the Omicron variant using wastewater that was independent of clinical testing (Arts et al., 2022).
The N200 multiplex qPCR was developed to improve estimates of the proportion of various variants in wastewater, while also making the assay more efficient to use in terms of reduced use of extracts, reagents, and time. The assay was validated using synthetic templates and was successfully applied to monitor for the emergence of the Omicron variant in several communities in Ontario. Although, this assay is one among several now available for detection of VOCs in wastewater, the unique design, with a single amplicon, multiple target mutations, and universal comparator targeting a highly mutable area of the N-gene make it very useful in wastewater-based surveillance of SARS-CoV-2.
Data Availability
All data produced in the present work are available upon reasonable request to the authors.
Acknowledgements
Funding for this project was provided by the Ontario Ministry of the Environment, Conservation and Parks, and NSERC Discovery and Canada Research Chairs. The research published in this paper is part of the project titled “Next generation solutions to ensure healthy water resources for future generations” funded by the Global Water Futures program, Canada First Research Excellence Fund. Additional information is available at www.globalwaterfutures.ca. The Municipal and Public Health Unit staff contributed to sample collection and provided valuable advice. Codey Dueck (PHAC) for laboratory analysis supporting LOD determinations. Yash Badlani, Carly Sing-Judge, Emily Dodsworth and Samantha Maerten assisted with sample preparation and analysis.
References
