Quantitative SARS-CoV-2 tracking of variants Delta, Delta plus, Kappa and Beta in wastewater by allele-specific RT-qPCR ======================================================================================================================= * Wei Lin Lee * Xiaoqiong Gu * Federica Armas * Franciscus Chandra * Hongjie Chen * Fuqing Wu * Mats Leifels * Amy Xiao * Feng Jun Desmond Chua * Germaine WC Kwok * Shreshtha Jolly * Claire YJ Lim * Janelle Thompson * Eric J Alm ## ABSTRACT The Delta (B.1.617.2) variant has caused major devastation in India and other countries around the world. First detected in October 2020, it has now spread to more than 100 countries, prompting WHO to declare it as a global variant of concern (VOC). The Delta (B.1.617.2), Delta plus (B.1.617.2.1) and Kappa (B.1.617.1) variants are all sub-lineages of the original B.1.617 variant. Prior to the inception of B.1.617, vaccine rollout, safe-distancing and timely lockdowns greatly reduced COVID-19 hospitalizations and deaths. However, the Delta variant, allegedly more infectious and for which existing vaccines seemed less effective, has catalyzed the resurgence of cases. Therefore, there is an imperative need for increased surveillance of the B.1.617 variants. While the Beta variant is increasingly outpaced by the Delta variant, the spread of the Beta variant remains of concern due to its vaccine resistance. Efforts have been made to utilize wastewater-based surveillance for community-based tracking of SARS-CoV-2 variants, however wastewater with its low SARS-CoV-2 viral titers and mixtures of viral variants, requires assays to be variant-specific yet accurately quantitative for meaningful interpretation. Following on the design principles of our previous assays for the Alpha variant, here we report allele-specific and multiplex-compatible RT-qPCR assays targeting mutations T19R, D80A, K417N, T478K and E484Q, for quantitative detection and discrimination of the Delta, Delta plus, Kappa and Beta variants in wastewater. This method is open-sourced and can be implemented using commercially available RT-qPCR protocols, and would be an important tool for tracking the spread of B.1.617 and the Beta variants in communities. Key Words * SARS-CoV-2 * COVID-19 * Variant * RT-qPCR * Wastewater * WBE * WBS ## INTRODUCTION The coronavirus disease-2019 (COVID-19) pandemic, caused by the etiological agent SARS-CoV-2, has triggered major public health, economic and psycho-social consequences worldwide (Zoumpourlis et al., 2020). Over the course of the pandemic, new viral genetic variants of SARS-CoV-2 have emerged, yielding variants of concern (VOCs) and variants of interest (VOIs) (Lauring and Hodcroft, 2021). Variants may possess increased transmissibility, more severe disease manifestation, or decrease effectiveness of available medical countermeasures like vaccinations and warrants close monitoring (CDC, 2021). Kappa (B.1.617.1) and Delta (B.1.617.2) variants first emerged in India in October, 2020 (PHE, 2021). The Delta variant has been predicted to be 40-60 % more transmissible than the Alpha variant (Mahase, 2021; PHE, 2021). In comparison to the Alpha variant, the Delta variant is more resistant to vaccine- and infection-induced immunity (Lustig et al., 2021; Bernal et al., 2020), with reduced sensitivity to antibody neutralization (Planas et al., 2021a). The Delta variant has spread rapidly to 104 countries, catalysing waves of COVID-19 infections worldwide (WHO, 2021). In mid March 2021, sequences of the Delta variant with a mutation K417N in the spike protein were detected. This was dubbed the Delta plus (B.1.617.2.1) variant. This K417N mutation is also present in the VOC Beta (B.1.351), which was first discovered in South Africa in 2020 and now reported in 123 countries. While its global prevalence has fallen in the presence of the Delta variant, there remains vigilance on the spread of the Beta variant due to it being more vaccine-resistant than others in circulation (Planas et al., 2021b, Charmet et al., 2021). The emergence of VOCs and VOIs means that tracking their introduction and spread in populations becomes essential to manage the pandemic. This has mainly been performed via genomic sequencing of clinical samples (Behrmann and Spiegel, 2020). However genomic sequencing of individual clinical samples is expensive and requires specialized infrastructure (Gwinn et al., 2019). Furthermore, for effective surveillance, a significant fraction of infected individuals need to be tested. A companion to clinical surveillance is wastewater-based surveillance, which has been shown to be effective at capturing temporal trends in viral circulation during this COVID-19 pandemic (Hassard et al., 2021; Thompson et al., 2020; Wu et al., 2020a). Current mainstream methods for variant tracking in wastewater rely on enriching and sequencing of the environmental SARS-CoV-2 genome (Napit et al., 2021, Crits-Christoph et al., 2021; Fontenele et al., 2021). This method, constrained by the same bottlenecks of cost and infrastructure requirement as clinical sequencing, is further hampered by challenges in detection of low-frequency variants and is poorly quantitative due to the lack of robust modelling to quantify variant titers (Fuqua et al., 2021; Van Poelvoorde et al., 2021). RT-qPCR-based methods have been developed for variant identification in clinical samples (Clark et al., 2021; Vogels et al., 2021; Wang et al., 2021) but are not yet validated for quantification of variant mixtures, which are expected in wastewater samples. This motivates the development of specialized methods for quantification of variant mixtures in environmental samples such as wastewater (Graber et al., 2021; Yaniv et al., 2021a, 2021b; Wurtzer et al., 2021). While RT-qPCR methods cannot discover new variants or identify variants beyond what they are designed for, they are more sensitive and quantitative, providing a readout from sample to data, in just hours. In our previous work, we developed and validated Allele-Specific RT-qPCR (AS RT-qPCR) assays for quantitative detection of Alpha variant B.1.1.7 in wastewater, tracking its occurrence over time in 19 communities across the United States (Lee et al., 2021). Here, we develop a similar AS RT-qPCR-based assay for the tracking of the Delta, Delta plus, Kappa and Beta variants by targeting five loci - T19R, D80A, K417N, T478K and E484Q that would identify and differentiate these variants. These assays could be performed as individual reactions or for increased throughput, as triplexes. The AS RT-qPCR assay for variant tracking is easily implementable in the conventional RT-qPCR based surveillance workflow commonly used worldwide. ## RESULTS AND DISCUSSION Following the methodology established in our previous study demonstrating single nucleotide discrimination for mutations associated with the SARS-CoV-2 Alpha variant (Lee et al., 2021), we developed an Allele-Specific RT-qPCR (AS RT-qPCR) panel for tracking mutations indicative of the SARS-CoV-2 variants Delta (B.1.617.2), Delta plus (B.1.617.2.1), Kappa (B.1.617.1) and Beta (B.1.351) and validated this approach for synthetic mixtures of Beta and Kappa VOC RNA in a wastewater RNA matrix (Delta assay validation is currently in progress, subject to availability of commercial RNA standards). To design allele specific primers we screened a panel of primers targeting mutations characteristic of four variants - Delta, Delta plus, Kappa and Beta, and identified primers targeting five loci - T19R, D80A, K417N, T478K and E484Q as having optimal sensitivity and specificity. These five targets are highly predictive for detection and discrimination of Delta, Delta plus, Kappa and Beta variants among currently circulating strains (**Figure 1**). The Kappa variant can be determined using positive tests for E484Q and the Delta variant by T19R and/or T478K. The Delta plus variant is indicated by additional detection of K417N on top of the ones for the Delta variant. K417N is also present in the Beta variant, which can be indicated by positive tests for D80A. ![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/08/06/2021.08.03.21261298/F1.medium.gif) [Figure 1.](http://medrxiv.org/content/early/2021/08/06/2021.08.03.21261298/F1) Figure 1. Schematic for deciding which assays to run for identification of Delta, Delta plus, Kappa and Beta variants. **Left**,Venn diagram showing targeted mutations that are representative of each variant. **Right**, Single assays for identification and discrimination of Delta, Delta plus, Kappa and Beta variants. E484Q is unique to the Kappa variant. D80A is unique to the Beta variant. Either T19R or T478K is representative of Delta (inclusive of Delta plus). Delta plus (without Delta) is indicated by subtracting the abundance of K417N (present in both Delta and Beta) by the abundance of D80A (unique to Beta). To exclusively indicate the abundance of the Delta variant, without that of the Delta plus variant, would require subtracting the abundance of the Delta variant by that of the Delta plus. Primers and probes for all assays shown are described in **Table 2**. We tabulated the presence of these mutations across SARS-CoV-2 variants **(Table 1**) to ensure the absence of cross reactivity with other variants. Both T19R and T478K are unique to the Delta and the Delta plus variant. K417N is present in both Beta and Delta plus. E484Q is unique to the Kappa variant, though variants Beta, Gamma and Eta possess a different mutation, E484K at the same site. View this table: [Table 1.](http://medrxiv.org/content/early/2021/08/06/2021.08.03.21261298/T1) Table 1. Characteristic mutations at target loci in Wild Type SARS-CoV-2 and its variants. Mutations to the WT sequence are shown in red (CDC, 2021). (-) denotes no change from wild type SARS-CoV-2. View this table: [Table 2.](http://medrxiv.org/content/early/2021/08/06/2021.08.03.21261298/T2) Table 2. AS RT-qPCR primer sequences developed in this paper for SARS-CoV-2 variants. Allele-specific nucleotides are marked bold. Synthetic mismatches are underlined. Assays can be run as individual singleplex assays (**Figure 1**) or for increased throughput, as triplex assays depending on the variants to be tracked (**Figure 2**). Triplex setup 1 enables identification of Kappa (E484Q), Delta (including Delta plus) (T19R) and Beta (D80A) variants within the same reaction. This triplex does not enable discrimination between Delta and Delta plus variants. Triplex setup 2 enables identification of Delta (T19R), Delta plus (T19R, K417N) and Beta (D80A, K417N) variants, using D80A to differentiate Delta plus from the Beta variant. K417/N and E484/Q cannot be combined in the same multiplex due to primer interference. ![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/08/06/2021.08.03.21261298/F2.medium.gif) [Figure 2.](http://medrxiv.org/content/early/2021/08/06/2021.08.03.21261298/F2) Figure 2. Triplexes for increased throughput of variant detection and quantitation. **Left**, Triplex setup 1 enables identification within the same reaction, the Kappa variant through E484Q, the Delta variant (including Delta plus) through T19R and the Beta variant through D80A. **Right**, Triplex setup 2 enables identification and differentiation of the Delta, Delta plus and Beta variants. The Delta plus variant contains K417N on top of the T19R present in both Delta and Delta plus variants. The Beta variant contains D80A and shares K417N with the Delta plus variant. As such Delta plus is indicated by excluding K417N by D80A. Primers targeting K417/N and E484/Q cannot be combined in the same multiplex due to primer interference. ### Validation of assays that indicate the Beta and Kappa variants Here we present analytical validation of the sensitivity and specificity of these five pairs of primer-probe sets. Setups are shown in **Table 3**. To improve readability, primers targeting WT are indicated as the amino acid, followed by the position of the amino acid in the protein sequence i.e. T19, D80, K417, T478 and E484 while primers targeting the variants are indicated by a suffix designating the mutant residue i.e. T19R, D80A, K417N, T478K and E484Q. Loci are respectively named T19/R, D80/A, K417/N, T478/K and E484/Q to refer to the loci on both WT and mutant sequences. View this table: [Table 3.](http://medrxiv.org/content/early/2021/08/06/2021.08.03.21261298/T3) Table 3. The AS RT-qPCR panel for identification and discimination of Delta, Delta plus, Kappa and Beta variants. We examined the specificity of the AS RT-qPCR assays for Beta and Kappa variant loci D80/A, K417/N and E484/Q by screening them against their respective WT and mutant genome targets in the SARS-CoV-2 RNA (**Figure 3**). Specificity data for T19/R and T478/K is shown in **Figure 4**. Synthetic RNA of the Beta variant was used for D80A and K417N, and Kappa RNA for E484Q. The amplification efficiency of each primer and probe set were between 89.3 to 104% for the correct RNA (i.e. Kappa assay for Kappa RNA, Beta assay for Beta RNA and WT assay for WT RNA). All three assays for mutant sequences and two of three assays for WT sequences were highly specific with cross reactivity only observed at or above 103 copies of non-target sequence per PCR. WT assay K417 had lower specificity indicating cross reactivity at 102 RNA copies of the mutant genotype. All assays were sufficiently specific to support quantification of mixtures at concentrations of SARS-CoV-2 RNA commonly observed in wastewater RNA preparations (Wu et al., 2020a, b, 2021). ![](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/08/06/2021.08.03.21261298/F3/graphic-7.medium.gif) [](http://medrxiv.org/content/early/2021/08/06/2021.08.03.21261298/F3/graphic-7) ![](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/08/06/2021.08.03.21261298/F3/graphic-8.medium.gif) [](http://medrxiv.org/content/early/2021/08/06/2021.08.03.21261298/F3/graphic-8) Figure 3. Specificity and cross-reactivity of the three sets of primers against full length synthetic WT, Beta and Kappa variant RNA. Full length synthetic RNA for the Delta and Delta plus variants were not yet commercially available at the time of this work. Primers targeting loci D80/A and K417/N are screened against WT and Beta variant RNA. Primers targeting loci E484/Q are screened against WT and Kappa variant RNA. Exact concentrations of synthetic RNA used were obtained by RT-ddPCR. Colored symbols represent tests against the matching genotype (WT-specific primers to WT RNA and mutant specific primers to mutant RNA) and grey symbols denote tests against RNA of the opposite genotype. Diamonds and squares denote tests using primers designed to target WT and variant RNA respectively. RT-qPCR efficiency and y-intercept cycle threshold (Ct) values were calculated for the primers against their respective RNA target sequences and shown in the table. ![Figure 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/08/06/2021.08.03.21261298/F4.medium.gif) [Figure 4.](http://medrxiv.org/content/early/2021/08/06/2021.08.03.21261298/F4) Figure 4. Specificity and cross-reactivity of the two sets of primers against synthetic DNA containing WT or Delta variant sequences. As full length synthetic Delta variant RNA was not yet commercially available at the time of this work, primers targeting loci T19 and T478 are screened against **a**) WT and mutant synthetic DNA and **b**) full length synthetic WT RNA. Exact concentrations of synthetic DNA and RNA were obtained by ddPCR and RT-ddPCR respectively. Data shown reflect two sets of independent measurements taken on different days. Colored symbols represent tests against the matching genotype (WT-specific primers to WT DNA and mutant specific primers to mutant DNA) and grey symbols denote tests against DNA of the opposite genotype. Diamonds and squares denote tests using primers designed to target WT and variant templates respectively. qPCR and RT-qPCR efficiencies and y-intercept cycle threshold (Ct) values were calculated for the primers against their respective RNA target sequences and shown in the table. ### Validation of assays T19R and T478K that indicate the Delta variant As commercial Delta variant RNA was not available in the laboratory at the time of this work, primers targeting Delta variant loci T19/R and T478/K were validated against synthetic DNA containing WT and variant sequences (**Figure 4a**) (DNA sequences shown in **Table 5**) and full length synthetic WT RNA (**Figure 4b**). The assays targeting T478/K appeared more specific than those targeting T19/R, though all remain discriminant up to at least 102 copies of DNA of the opposite genotype. Screening against WT RNA (**Figure 4b**), the amplification efficiencies for the T19 and T478 primers were between 74.0-98.2%. Mutant primers targeting T19R and T478K did not cross react with up to 103 copies of WT RNA. ### Validation of the singleplex and triplex assays in wastewater RNA The validation so far suggests most primer-probe sets to be discriminant, with minimal cross-reactivity against at least 102 copies of the opposite genotype. To improve throughput of these assays, we explored which of these primer-probe sets could be multiplexed by combining primers targeting WT loci in the same reaction and primers targeting mutant loci in the other reaction. We found that K417/N and E484/Q cannot be combined in the same multiplex due to primer interference. Further, while primers targeting T478/K were more specific than those for T19/R for identifying the Delta variant, only the latter was incorporated into the triplex assays given overlapping primers, with T478/K being very close in proximity to E484Q. As such, we developed two triplexes (**Figure 2, Table 4**). The first triplex (Triplex 1) can identify Kappa (E484Q), Delta (including Delta plus) (T19R) and Beta (D80A) variants within the same reaction. Triplex 2 can identify and differentiate Delta (T19R), Delta plus (T19R, K417N) and Beta (D80A, K417N) variants, using D80A to differentiate Delta plus from the Beta variant. Amplification efficiencies between the singleplex and the two triplexes appeared similar for most of the loci (**Figure 5 and 6**). View this table: [Table 4.](http://medrxiv.org/content/early/2021/08/06/2021.08.03.21261298/T4) Table 4. Setups for AS RT-qPCR triplexes. All primers and probes are used at a final concentration of 500 nM and 200 nM respectively. ![](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/08/06/2021.08.03.21261298/F5/graphic-12.medium.gif) [](http://medrxiv.org/content/early/2021/08/06/2021.08.03.21261298/F5/graphic-12) ![](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/08/06/2021.08.03.21261298/F5/graphic-13.medium.gif) [](http://medrxiv.org/content/early/2021/08/06/2021.08.03.21261298/F5/graphic-13) Figure 5. Ct value comparison of singleplex to multiplex AS RT-qPCR using 10-fold serial dilutions of full length synthetic RNA for loci T19, D80/A, K417/N and E484/K. **Left panel**, WT primers are assayed against WT RNA. **Right panel**, primers targeting D80A and K417N were assayed against Beta RNA, and primers targeting E484Q were assayed against Kappa RNA. Diamonds and squares denote tests using primers designed to target WT and variant templates respectively. T19R is not tested against the RNA of the Delta variant and only screened using synthetic DNA templates containing WT and Delta variant sequences, shown in Figure 6. Singleplex and multiplex reactions were performed in the presence of 500 nM of primers and 200 nM of probes. RT-qPCR efficiencies and y-intercept cycle threshold (Ct) values were calculated for the primers against their respective RNA target sequences and shown in the tables. Replicate values represent that of the triplex 1 and 2 respectively. ![Figure 6.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/08/06/2021.08.03.21261298/F6.medium.gif) [Figure 6.](http://medrxiv.org/content/early/2021/08/06/2021.08.03.21261298/F6) Figure 6. Ct value comparison of multiplex to singleplex AS RT-qPCR using 10-fold serial dilutions of synthetic DNA for loci T19/R. **Left**, WT T19 primers are assayed against synthetic DNA with WT sequences. **Right**, Delta variant targeting T19R primers were assayed against synthetic DNA with Delta variant sequences. ![Figure 7.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/08/06/2021.08.03.21261298/F7.medium.gif) [Figure 7.](http://medrxiv.org/content/early/2021/08/06/2021.08.03.21261298/F7) Figure 7. Partial validation of singleplex and triplex assays in a mixture of WT, Beta and Kappa RNA in the presence and absence of wastewater RNA. (Note: Delta variant was not assessed as the synthetic RNA is not yet available in the laboratory). Kappa, Beta and WT RNA were added at 25:25:950, 100:100:800 and 400:400:200 copies per reaction respectively to make up 2.5%, 10% and 40% of Beta and Kappa RNA against a total of 1000 copies of SARS-CoV-2 RNA per reaction. WW refers to assays performed in the presence of wastewater RNA which contained PMMoV at an average of Ct 32–34. Non-specific amplification of delta variant loci was not observed for triplex 1 or 2. We tested the singleplex and triplex assays in a mixture of WT, Beta and Kappa RNA in the presence and absence of wastewater RNA derived from SARS-CoV-2 negative wastewater. We confirm no amplification with the singleplex and triplex assays targeting Beta, Kappa and Delta variants against wastewater RNA without adding synthetic SARS-CoV-2 RNA. Both the assays (singleplex and triplex) were able to robustly detect synthetic SARS-CoV-2 variant RNA, down to an abundance of 2.5% against a backdrop of 1000 SARS-CoV-2 RNA copies. Similar Ct values were acquired in both assays, when conducted in the presence or absence of wastewater RNA. On average, SARS-CoV-2 RNA derived across wastewater samples tend to be below 103 viral copies per ml, which gives rise to less than 103 viral RNA per RT-qPCR reaction (Wu et al., 2020a, b, 2021). This confirms that the singleplex and multiplex assays are specific and quantitative in wastewater RNA. ## DISCUSSION The present study presents five assays that could be applied to detecting and differentiating SARS-CoV-2 variants Delta (B.1.617.2), Delta plus (B.1.617.2.1), Kappa (B.1.617.1) and Beta (B.1.351) in wastewater. We followed the screening and validation strategy developed in our previous work on the alpha variant to yield reliable assays that enable discrimination of single nucleotide mutations (Lee et al., 2021). Primer-probe sets were developed for targeting both WT and mutant sequences on loci T19/R, D80/A, K417/N, T478/K and E484/Q. To increase throughput, we propose primer-probe sets that could be combined into triplexes for parallel interrogation of more loci within the same reaction. One of these triplexes enables detection of the Delta, Kappa and Beta variants while the other detects and discriminates Delta and Beta from Delta plus variants. We confirm low cross-reactivity for all the primers to at least 102 copies of RNA of the opposite genotype and similar quantitative performance of singleplex and triplex assays. Assays for T19/R and T478/K were only validated on synthetic DNA carrying WT and mutant sequences and synthetic WT RNA. Amplification efficiencies were only derived for T19/R and T478/K against synthetic WT RNA and not Delta variant RNA. However WT and mutant assays differ by only a single nucleotide and hybridization kinetics should be similar. Nonetheless detailed analytical validation against the Delta variant RNA will be performed and reported as an update to this preprint. In this work, we developed primer-probe sets targeting both WT and mutant loci, though variant detection and quantitation would only require performing the reactions targeting the mutant loci. As shown in our previous manuscript (Lee et al., 2021), quantification of both WT and mutant loci in wastewater enables determination of the proportion of WT to variant sequence at each target loci. Validation of the assays reported herein in SARS-CoV-2 positive wastewater is commencing and will be reported when available. This work expands on the utility of AS RT-qPCR for the quantitative detection of variants in wastewater for population-based tracking of SARS-CoV-2 variants. ## MATERIALS AND METHODS ### Assay design We designed AS RT-qPCR reactions to detect five non-synonymous single nucleotide variants in the spike gene. Primers and probes were designed following our previous work (Lee et al., 2021) and using the Integrated DNA Technologies (IDT)’s PrimerQuest Tool. Target mutations were placed at the 3’ end of either the forward or reverse primer. All primers were designed to have a melting temperature in the range of 59–65°C and the probes in the range of 64–72°C. Probes were designed to anneal to the same strand as the AS primer, with the probe as close to the 3′-end of the AS primers as possible. Guanines are avoided at the 5′-end of the probe. AS primers were designed to include an artificial mismatch near the 3’ terminal nucleotide to improve discrimination. All primers and probes were purchased from IDT (**Table 2**). ### DNA standards and their quantification by ddPCR DNA standards were synthesised by IDT. Sequences (shown in **Table 5**) were designed to span the primer and probe sets of each assay and then concatenated into a single sequence. This DNA were quantified using ddPCR Supermix for Probes #1863026 (Bio-Rad) following manufacturer’s recommendations. Reaction mixtures consisted of 10 μl of 2× Supermix, 900nM primers, 250 nM probe, 1 μl of DNA, topped up with molecular grade water to a final volume of 20 μl. A Bio-Rad QX200 ddPCR droplet generator (Bio-Rad, USA) was used to convert the reaction mix into droplets. Thermal cycling was performed on a Bio-Rad CFX96 as follows: 10 min at 95°C, followed by 40 cycles of 30 s at 94°C and 1 min at 60°C (ramp rate of 2°C/sec), followed by enzyme inactivation at 98°C for 10 min and holding at 4°C. View this table: [Table 5.](http://medrxiv.org/content/early/2021/08/06/2021.08.03.21261298/T5) Table 5. Sequences of DNA standards used for validation of primer sets. ### Analysis of assays against DNA standards by qPCR qPCR reactions were performed using TaqMan Fast advanced Master Mix (ThermoFisher #4444556) at a final volume of 10 μL in duplicates, according to the manufacturer’s recommendations with a final primer concentration of 0.5 μM and probe of 0.2 μM with 1 μL of DNA template, at 20 s at 95 °C, 40 cycles of 3 s at 95 °C and 30 s at 60 °C) A single reverse or forward primer and probe was used with each set of allele-specific forward or reverse primers (**Table 3**) ### RNA Standards and their quantification by ddPCR Twist synthetic SARS-CoV-2 RNA Controls, control 2 (Wuhan-Hu-1), control 16 (South Africa/KRISP-EC-K005299/2020), control 18 (India/CT-ILSGS00361/2021), were used as RNA standards, representing WT, Beta (B.1.351) and Kappa (B.1.617.1) variants respectively. RNA standards were prepared as single-use aliquots. Controls 2, 16, and 18 were quantified by digital droplet PCR (ddPCR) to be 4.11 × 105, 8.08 × 105, and 3.24 × 105 copies/μL, respectively. Quantification was performed using One-Step RT-ddPCR Advanced Kit for Probes #1864022 (Bio-Rad) following manufacturer’s recommendations. ### Analysis of assays against RNA standards by RT-qPCR AS RT-qPCR was performed using the Taqman Virus 1-Step master mix (Thermofisher #4444434) with technical duplicates, at a final volume of 10 µL, according to the manufacturer’s recommendations. For singleplex reactions, a single reverse or forward primer and probe was used with each set of allele-specific forward or reverse primers (**Table 3**). The final concentration of the AS RT-qPCR primers were 500 nM, probe at 200 nM, with 1 µL of template. Multiplex reactions were set up with wild-type targeting primer-probe sets in the same reaction, and mutation targeting primer-probe sets in the other reaction, as shown in **Table 4**, with final concentration of each AS RT-qPCR primer as 500 nM and probe at 200 nM. No template controls were included for each assay and none of them amplified. The reactions are setup using electronic pipettes (Eppendorf) and performed on a Bio-Rad CFX384 real-time PCR instrument under the following conditions, 5 min at 50 °C and 20 s at 95 °C, followed by 45 cycles of 3 s at 95 °C and 30 s at 60 °C. ## Data Availability Source data will be made available upon request. ## Data analysis Data was analysed using Microsoft Excel and Graphpad prism. Graphs were presented using Graphpad Prism. ## Data availability Source data will be made available upon request. ## Declaration of competing interests EJA is an advisor to Biobot Analytics and holds shares in the company. ## Funding Statement This research is supported by the National Research Foundation, Prime Minister’s Office, Singapore, under its Campus for Research Excellence and Technological Enterprise (CREATE) program funding to the Singapore-MIT Alliance for Research and Technology (SMART) Antimicrobial Resistance Interdisciplinary Research Group (AMR IRG) and the Intra-CREATE Thematic Grant (Cities) grant NRF2019-THE001-0003a to JT and EJA and funding from the Singapore Ministry of Education and National Research Foundation through an RCE award to Singapore Centre for Environmental Life Sciences Engineering (SCELSE). ## Contributions EJA and JT conceptualized the project. WLL designed the experiments. WLL, XG, FA, FC, HC, FW, ML, AX, SJ and CYJL analyzed the data. WLL, FJDC, SJ and CYJL performed experiments. All authors contributed to writing the manuscript. WLL, JT and EJA supervised the project. All authors read and approved the manuscript. ## Acknowledgements We thank members of the Biobot Analytics, Inc team for helpful discussions. * Received August 3, 2021. * Revision received August 3, 2021. * Accepted August 6, 2021. * © 2021, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution 4.0 International), CC BY 4.0, as described at [http://creativecommons.org/licenses/by/4.0/](http://creativecommons.org/licenses/by/4.0/) ## REFERENCES 1. 1.Behrmann, O., Spiegel, M., 2020. COVID-19: from rapid genome sequencing to fast decisions. Lancet Infect. Dis. 20, 1218. [https://doi.org/10.1016/S1473-3099(20)30580-6](https://doi.org/10.1016/S1473-3099(20)30580-6) 2. 2.Bernal, J.L., Andrews, N., Gower, C., Gallagher, E., Simmons, R., Thelwall, S., Stowe, J., Tessier, E., Groves, N., Dabrera, G., Myers, R., Campbell, C., Amirthalingam, G., Edmunds, M., Zambon, M., Brown, K., Hopkins, S., Chand, M., Ramsay, M., 2021. Effectiveness of COVID-19 vaccines against the B.1.617.2 variant. medRxiv 2021.05.22.21257658. [https://doi.org/10.1101/2021.05.22.21257658](https://doi.org/10.1101/2021.05.22.21257658) 3. 3.CDC, 2021. SARS-CoV-2 Variant Classifications and Definitions [WWW Document]. URL [https://www.cdc.gov/coronavirus/2019-ncov/variants/variant-info.html](https://www.cdc.gov/coronavirus/2019-ncov/variants/variant-info.html) (accessed 7.9.21). 4. 4.Clark, A.E., Wang, Z., Cantarel, B., Kanchwala, M., Xing, C., Chen, L., Irwin, P., Xu, Y., Oliver, D., Lee, F., Gagan, J.R., Filkins, L., Muthukumar, A., Park, J.Y., Sarode, R., SoRelle, J.A., 2021. Multiplex Fragment Analysis Identifies SARS-CoV-2 Variants. medRxiv 2021.04.15.21253747. [https://doi.org/10.1101/2021.04.15.21253747](https://doi.org/10.1101/2021.04.15.21253747) 5. 5.Crits-Christoph, A., Kantor, R.S., Olm, M.R., Whitney, O.N., Al-Shayeb, B., Lou, Y.C., Flamholz, A., Kennedy, L.C., Greenwald, H., Hinkle, A., Hetzel, J., Spitzer, S., Koble, J., Tan, A., Hyde, F., Schroth, G., Kuersten, S., Banfield, J.F., Nelson, K.L., 2021. Genome Sequencing of Sewage Detects Regionally Prevalent SARS-CoV-2 Variants. MBio 12, e02703–20. [https://doi.org/10.1128/mBio.02703-20](https://doi.org/10.1128/mBio.02703-20) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1128/mBio.02703-20&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=33468686&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F06%2F2021.08.03.21261298.atom) 6. 6.Fontenele, R.S., Kraberger, S., Hadfield, J., Driver, E.M., Bowes, D., Holland, L.A., Faleye, T.O.C., Adhikari, S., Kumar, R., Inchausti, R., Holmes, W.K., Deitrick, S., Brown, P., Duty, D., Smith, T., Bhatnagar, A., Yeager, R.A., Holm, R.H., von Reitzenstein, N.H., Wheeler, E., Dixon, K., Constantine, T., Wilson, M.A., Lim, E.S., Jiang, X., Halden, R.U., Scotch, M., Varsani, A., 2021. High-throughput sequencing of SARS-CoV-2 in wastewater provides insights into circulating variants. medRxiv 2021.01.22.21250320. [https://doi.org/10.1101/2021.01.22.21250320](https://doi.org/10.1101/2021.01.22.21250320) 7. 7.Fuqua, J.L., Rouchka, E.C., Waigel, S., Sokoloski, K., Chung, D., Zacharias, W., Zhang, M., Chariker, J., Talley, D., Santisteban, I., Varsani, A., Moyer, S., Holm, R.H., Yeager, R.A., Smith, T., Bhatnagar, A., 2021. A rapid assessment of wastewater for genomic surveillance of SARS-CoV-2 variants at sewershed scale in Louisville, KY. medRxiv 2021.03.18.21253604. [https://doi.org/10.1101/2021.03.18.21253604](https://doi.org/10.1101/2021.03.18.21253604) 8. 8.Graber, T.E., Mercier, É., D’Aoust, P.M., Hoang, H.-D., Tian, X., Tasneem, S., Bhatnagar, K., Delatolla, R., 2021. An allele-specific primer extension assay to quantify the proportion of B.1.1.7-specific SARS-CoV-2 RNA in wastewater. medRxiv 2021.02.22.21252041. [https://doi.org/10.1101/2021.02.22.21252041](https://doi.org/10.1101/2021.02.22.21252041) 9. 9.Gwinn, M., MacCannell, D., Armstrong, G.L., 2019. Next-Generation Sequencing of Infectious Pathogens. JAMA 321, 893–894. [https://doi.org/10.1001/jama.2018.21669](https://doi.org/10.1001/jama.2018.21669) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jama.2018.21669&link_type=DOI) 10. 10.Hassard, F., Lundy, L., Singer, A.C., Grimsley, J., Di Cesare, M., 2021. Innovation in wastewater near-source tracking for rapid identification of COVID-19 in schools. The Lancet Microbe 2, e4–e5. [https://doi.org/10.1016/S2666-5247(20)30193-2](https://doi.org/10.1016/S2666-5247(20)30193-2) 11. 11.Lauring, A.S., Hodcroft, E.B., 2021. Genetic Variants of SARS-CoV-2—What Do They Mean? JAMA 325, 529–531. [https://doi.org/10.1001/jama.2020.27124](https://doi.org/10.1001/jama.2020.27124) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F06%2F2021.08.03.21261298.atom) 12. 12.Lee, W.L., Imakaev, M., Armas, F., McElroy, K.A., Gu, X., Duvallet, C., Chandra, F., Chen, H., Leifels, M., Mendola, S., Floyd-O’Sullivan, R., Powell, M.M., Wilson, S.T., Berge, K.L.J., Lim, C.Y.J., Wu, F., Xiao, A., Moniz, K., Ghaeli, N., Matus, M., Thompson, J., Alm, E.J., 2021. Quantitative SARS-CoV-2 Alpha Variant B.1.1.7 Tracking in Wastewater by Allele-Specific RT-qPCR. Environ. Sci. Technol. Lett. [https://doi.org/10.1021/acs.estlett.1c00375](https://doi.org/10.1021/acs.estlett.1c00375) 13. 13.Lustig, Y., Zuckerman, N., Nemet, I., Atari, N., Kliker, L., Regev-Yochay, G., Sapir, E., Mor, O., Alroy-Preis, S., Mendelson, E., Mandelboim, M., 2021. Neutralising capacity against Delta (B.1.617.2) and other variants of concern following Comirnaty (BNT162b2, BioNTech/Pfizer) vaccination in health care workers, Israel. Euro Surveill. Bull. Eur. sur les Mal. Transm. = Eur. Commun. Dis. Bull. 26. [https://doi.org/10.2807/1560-7917.ES.2021.26.26.2100557](https://doi.org/10.2807/1560-7917.ES.2021.26.26.2100557) 14. 14.Mahase, E., 2021. Delta variant: What is happening with transmission, hospital admissions, and restrictions? BMJ 373, 1513. [https://doi.org/10.1136/bmj.n1513](https://doi.org/10.1136/bmj.n1513) 15. 15.Napit, R., Manandhar, P., Chaudhary, A., Shrestha, B., Poudel, A., Raut, R., Pradhan, S., Raut, S., Mathema, S., Rajbhandari, R., Dixit, S., Schwind, J.S., Johnson, C.K., Mazet, J.K., Karmacharya, D., 2021. Rapid genomic surveillance of SARS-CoV-2 in a dense urban community using environmental (sewage) samples. medRxiv 2021.03.29.21254053. [https://doi.org/10.1101/2021.03.29.21254053](https://doi.org/10.1101/2021.03.29.21254053) 16. 16.PHE, 2021. SARS-CoV-2 variants of concern and variants under investigation in England [WWW Document]. URL [https://assets.publishing.service.gov.uk/government/uploads/system/uploads/attachment\_data/file/997414/Variants\_of\_Concern\_VOC\_Technical\_Briefing\_16.pdf](https://assets.publishing.service.gov.uk/government/uploads/system/uploads/attachment\_data/file/997414/Variants\_of\_Concern\_VOC_Technical_Briefing_16.pdf) 17. 17.Planas, D., Bruel, T., Grzelak, L., Guivel-Benhassine, F., Staropoli, I., Porrot, F., Planchais, C., Buchrieser, J., Rajah, M.M., Bishop, E., Albert, M., Donati, F., Prot, M., Behillil, S., Enouf, V., Maquart, M., Smati-Lafarge, M., Varon, E., Schortgen, F., Yahyaoui, L., Gonzalez, M., De Sèze, J., Péré, H., Veyer, D., Sève, A., Simon-Lorière, E., Fafi-Kremer, S., Stefic, K., Mouquet, H., Hocqueloux, L., van der Werf, S., Prazuck, T., Schwartz, O., 2021a. Sensitivity of infectious SARS-CoV-2 B.1.1.7 and B.1.351 variants to neutralizing antibodies. Nat. Med. 27, 917–924. [https://doi.org/10.1038/s41591-021-01318-5](https://doi.org/10.1038/s41591-021-01318-5) 18. 18.Planas, D., Veyer, D., Baidaliuk, A., Staropoli, I., Guivel-Benhassine, F., Rajah, M.M., Planchais, C., Porrot, F., Robillard, N., Puech, J., Prot, M., Gallais, F., Gantner, P., Velay, A., Le Guen, J., Kassis-Chikhani, N., Edriss, D., Belec, L., Seve, A., Courtellemont, L., Péré, H., Hocqueloux, L., Fafi-Kremer, S., Prazuck, T., Mouquet, H., Bruel, T., Simon-Lorière, E., Rey, F.A., Schwartz, O., 2021b. Reduced sensitivity of SARS-CoV-2 variant Delta to antibody neutralization. Nature. [https://doi.org/10.1038/s41586-021-03777-9](https://doi.org/10.1038/s41586-021-03777-9) 19. 19.Thompson, J.R., Nancharaiah, Y. V, Gu, X., Lee, W.L., Rajal, V.B., Haines, M.B., Girones, R., Ng, L.C., Alm, E.J., Wuertz, S., 2020. Making waves: Wastewater surveillance of SARS-CoV-2 for population-based health management. Water Res. 184, 116181. [https://doi.org/10.1016/j.watres.2020.116181](https://doi.org/10.1016/j.watres.2020.116181) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.watres.2020.116181&link_type=DOI) 20. 20.Van Poelvoorde, L.A.E., Delcourt, T., Coucke, W., Herman, P., De Keersmaecker, S.C.J., Saelens, X., Roosens, N., Vanneste, K., 2021. Strategy and performance evaluation of low-frequency variant calling for SARS-CoV-2 in wastewater using targeted deep Illumina sequencing. medRxiv 2021.07.02.21259923. [https://doi.org/10.1101/2021.07.02.21259923](https://doi.org/10.1101/2021.07.02.21259923) 21. 21.Vogels, C.B.F., Breban, M.I., Ott, I.M., Alpert, T., Petrone, M.E., Watkins, A.E., Kalinich, C.C., Earnest, R., Rothman, J.E., Goes de Jesus, J., Morales Claro, I., Magalhães Ferreira, G., Crispim, M.A.E., Network, B.-U.C.G., Singh, L., Tegally, H., Anyaneji, U.J., Africa, N. for G.S. in S., Hodcroft, E.B., Mason, C.E., Khullar, G., Metti, J., Dudley, J.T., MacKay, M.J., Nash, M., Wang, J., Liu, C., Hui, P., Murphy, S., Neal, C., Laszlo, E., Landry, M.L., Muyombwe, A., Downing, R., Razeq, J., de Oliveira, T., Faria, N.R., Sabino, E.C., Neher, R.A., Fauver, J.R., Grubaugh, N.D., 2021. Multiplex qPCR discriminates variants of concern to enhance global surveillance of SARS-CoV-2. PLOS Biol. 19, e3001236. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pbio.3001236&link_type=DOI) 22. 22.Wang, H., Miller, J.A., Verghese, M., Sibai, M., Solis, D., Mfuh, K.O., Jiang, B., Iwai, N., Mar, M., Huang, C., Yamamoto, F., Sahoo, M.K., Zehnder, J., Pinsky, B.A., 2021. Multiplex SARS-CoV-2 Genotyping PCR for Population-Level Variant Screening and Epidemiologic Surveillance. J. Clin. Microbiol. 59. [https://doi.org/10.1101/2021.04.20.21255480](https://doi.org/10.1101/2021.04.20.21255480) 23. 23.WHO, 2021. Weekly epidemiological update on COVID-19 - 6 July 2021 [WWW Document]. Emerg. Situational Update. URL [https://www.who.int/publications/m/item/weekly-epidemiological-update-on-covid-196-july-2021](https://www.who.int/publications/m/item/weekly-epidemiological-update-on-covid-196-july-2021) (accessed 7.13.21). 24. 24.Wu, F., Xiao, A., Zhang, J., Moniz, K., Endo, N., Armas, F., Bonneau, R., Brown, M.A., Bushman, M., Chai, P.R., Duvallet, C., Erickson, T.B., Foppe, K., Ghaeli, N., Gu, X., Hanage, W.P., Huang, K.H., Lee, W.L., Matus, M., McElroy, K.A., Nagler, J., Rhode, S.F., Santillana, M., Tucker, J.A., Wuertz, S., Zhao, S., Thompson, J., Alm, E.J., 2020a. SARS-CoV-2 titers in wastewater foreshadow dynamics and clinical presentation of new COVID-19 cases. medRxiv Prepr. Serv. Heal. Sci. [https://doi.org/10.1101/2020.06.15.20117747](https://doi.org/10.1101/2020.06.15.20117747) 25. 25.Wu, F., Zhang, J., Xiao, A., Gu, X., Lee, W.L., Armas, F., Kauffman, K., Hanage, W., Matus, M., Ghaeli, N., Endo, N., Duvallet, C., Poyet, M., Moniz, K., Washburne, A.D., Erickson, T.B., Chai, P.R., Thompson, J., Alm, E.J., 2020b. SARS-CoV-2 Titers in Wastewater Are Higher than Expected from Clinically Confirmed Cases. mSystems 5, e00614–20. [https://doi.org/10.1128/mSystems.00614-20](https://doi.org/10.1128/mSystems.00614-20) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1128/mSystems.00614-20&link_type=DOI) 26. 26.Wu, F., Xiao, A., Zhang, J., Moniz, K., Endo, N., Armas, F., Bushman, M., Chai, P.R., Duvallet, C., Erickson, T.B., Foppe, K., Ghaeli, N., Gu, X., Hanage, W.P., Huang, K.H., Lee, W.L., McElroy, K.A., Rhode, S.F., Matus, M., Wuertz, S., Thompson, J., Alm, E.J., 2021b. Wastewater surveillance of SARS-CoV-2 across 40 U.S. states from February to June 2020. Water Res. 202, 117400. [https://doi.org/10.1016/j.watres.2021.117400](https://doi.org/10.1016/j.watres.2021.117400) 27. 27.Wurtzer, S., Waldman, P., Levert, M., Mouchel, J.M., Gorgé, O., Boni, M., Maday, Y., Marechal, V., Moulin, L., 2021. Monitoring the propagation of SARS CoV2 variants by tracking identified mutation in wastewater using specific RT-qPCR. medRxiv 2021.03.10.21253291. [https://doi.org/10.1101/2021.03.10.21253291](https://doi.org/10.1101/2021.03.10.21253291) 28. 28.Yaniv, K., Ozer, E., Kushmaro, A., 2021a. SARS-CoV-2 variants of concern, Gamma (P.1) and Delta (B.1.617), sensitive detection and quantification in wastewater employing direct RT-qPCR. medRxiv 2021.07.14.21260495. [https://doi.org/10.1101/2021.07.14.21260495](https://doi.org/10.1101/2021.07.14.21260495) 29. 29.Yaniv, K., Ozer, E., Plotkin, N., Bhandarkar, N.S., Kushmaro, A., 2021b. RT-qPCR assay for detection of British (B.1.1.7) and South Africa (B.1.351) variants of SARS-CoV-2. medRxiv 2021.02.25.21252454. [https://doi.org/10.1101/2021.02.25.21252454](https://doi.org/10.1101/2021.02.25.21252454) 30. 30.Zoumpourlis, V., Goulielmaki, M., Rizos, E., Baliou, S., Spandidos, D.A., 2020. [Comment] The COVID-19 pandemic as a scientific and social challenge in the 21st century. Mol. Med. Rep. 22, 3035–3048. [https://doi.org/10.3892/mmr.2020.11393](https://doi.org/10.3892/mmr.2020.11393)