Droplet Digital RT-PCR to detect SARS-CoV-2 variants of concern in wastewater ============================================================================= * Leo Heijnen * Goffe Elsinga * Miranda de Graaf * Richard Molenkamp * Marion P.G. Koopmans * Gertjan Medema ## Abstract Wastewater surveillance has shown to be a valuable and efficient tool to obtain information about the trends of COVID-19 in the community. Since the recent emergence of new variants, associated with increased transmissibility and/or antibody escape (variants of concern), there is an urgent need for methods that enable specific and timely detection and quantification of the occurrence of these variants in the community. In this study we demonstrate the use of RT-ddPCR on wastewater samples for specific detection of mutation N501Y. This assay enabled simultaneous enumeration of the concentration of variants with the 501Y mutation and Wild Type (WT, containing 501N) SARS-CoV-2 RNA. Detection of N501Y was possible in samples with mixtures of WT with low proportions of lineage B.1.351 (0.5%). The method could accurately determine the proportion of N501Y and WT in mixtures of SARS-CoV-2 RNA. The application to raw sewage samples from the cities of Amsterdam and Utrecht demonstrated that this method can be applied to determine the concentrations and the proportions of WT and N501Y containing SARS-CoV-2 RNA in wastewater samples. The emergence of N501Y in Amsterdam and Utrecht wastewater aligned with the emergence of B.1.1.7 as causative agent of COVID-19 in the Netherlands, indicating that RT-ddPCR of wastewater samples can be used to monitor the emergence of the N501Y mutation in the community. It also indicates that RT-ddPCR could be used for sensitive and accurate monitoring of current (like K417N, E484K) or future mutations present in SARS-CoV-2 variants of concern. Monitoring emergence of these mutations in the community via wastewater is rapid, efficient and valuable in supporting public health decision-making. ## Introduction The first reports of diarrhea caused by SARS-CoV-2 in COVID-19 patients and the detection of SARS-CoV-2 RNA in the feces of 30-60% of infected individuals (Jones, Baluja et al. 2020, Wang, Hu et al. 2020, Wölfel, Corman et al. 2020, Xiao, Tang et al. 2020), sparked the idea that wastewater could be used to monitor the circulation of SARS-CoV-2 at community level. Methods for the detection of SARS-CoV-2 RNA in wastewater were developed and deployed already in the early stage of the pandemic (Lodder and de Roda Husman 2020, Medema, Heijnen et al. 2020). The concentration of SARS-CoV-2 RNA in wastewater has shown to reflect and even precede the trends of the newly reported cases or COVID-19 hospitalizations (Medema, Been et al. 2020, Prado, Fumian et al. 2021). The possibility to obtain objective information about virus circulation in a community via wastewater sampling, independent of diagnostic testing availability, willingness and awareness (i.e. asymptomatic virus carriers), is currently being used to support public health decisions. In situations with low SARS-CoV-2 circulation in the community, wastewater surveillance is being used as an early warning system (Betancourt, Schmitz et al. 2020, Medema, Heijnen et al. 2020, Ahmed, Tscharke et al. 2021) enabling rapid and targeted measures to limit viral spread. The added value of wastewater surveillance has led to worldwide implementation of methods for SARS-CoV-2 monitoring in sewage (Bivins, North et al. 2020). The occurrence of new SARS-CoV-2 variants that have unusual numbers of mutations and are associated with increased transmissibility, antibody escape or both (“Variants of Concern”, VoC), has recently gained attention (Altmann, Boyton et al. 2021, Challen, Brooks-Pollock et al. 2021, Eurosurveillance-editorial-team 2021). Three VoC lineages are currently spreading worldwide and additional variants are under consideration: Lineage B.1.1.7 was recognized first through the large scale genomics initiative in the UK, lineage B.1.351 was first recognized in South Africa and lineage P.1 emerged in late 2020 in Manaus (Brazil). These VoC’s appear to spread at higher efficiency (Davies, Abbott et al. 2021, Faria, Mellan et al. 2021, Volz, Mishra et al. 2021). In addition, lineages B.1.351 and P.1 appear to be partially resistant to neutralizing antibodies from infection, vaccination, or treatment (Garcia-Beltran, Lam et al. 2021, Wang, Liu et al. 2021). Therefore, it is of great importance to monitor the spread of these VoC in the community to understand the dynamics of transmission of VoC and be able to take appropriate public health protection measures. Currently, VoC surveillance is mainly done through whole genome sequencing of patient samples. This approach is relatively expensive, labour intensive and time consuming and may not be achievable in many parts of the world. Variant surveillance via wastewater could be a more efficient and rapid method to monitor the emergence and spread of VoC in a community. However, there are currently no methods available to quantify the presence of different VoC’s in sewage. Sequence analysis of RNA isolated from sewage can specifically detect and identify SARS-CoV-2 RNA sequence variation (Izquierdo-Lara, Elsinga et al. 2020) and sequence analysis has the potential to identify VoC in sewage. However, it is expected that this needs deep sequencing efforts and bioinformatics to find rare mutations in sewage samples which likely contain mixtures of SARS-CoV-2 lineages. The limited sequence differences between VoC and the likely presence of mixtures of sequence variants in sewage makes it difficult to design VoC specific RT-qPCR assays for sewage. In droplet digital RT-PCR (RT-ddPCR) the sample and reagents are partitioned in a large number (~10.000 to 20.000) “water in oil” droplets and PCR-reactions are performed on single molecules in individual closed droplets. This makes it possible to detect rare mutations and to discriminate closely related sequences on the base of probe binding kinetics (Pekin, Skhiri et al. 2011). This makes RT-ddPCR potentially attractive for specific detection presence and proportion of mutations that are specific to VoC’s in wastewater. In this study we evaluated the use of RT-ddPCR to detect the N501Y mutation and to determine the proportion and concentration of VoC containing N501Y in wastewater. The N501Y mutation is present in lineages B.1.1.7, B.1.351 and P.1 and leads to an amino acid change at position 501 in the receptor binding domain of the spike protein. The N501Y mutation appears to influence binding of the spike protein to the human ACE2 receptor (Liu, Zhang et al. 2021, Rynkiewicz, Babbitt et al. 2021). Validation was performed with RNA isolated from wild type virus and a B.1.351 stock and the applicability of the RT-ddPCR was studied on domestic wastewater samples from the Dutch cities of Amsterdam and Utrecht. ## Methods ### Wastewater samples Composite, flow-dependent 24h samples were collected by the operators of the wastewater treatment plants (WWTP) that serve the cities of Amsterdam and Utrecht in the Netherlands. Samples were stored at 4 °C during 24h periods of sampling as previously described (Medema et al 2020). Samples were taken weekly from March 2020 to March 1 2021. RT-qPCR was applied on all samples, RT-ddPCR analysis was applied on the samples from the period November 9 2020 to March 1 2021. ### SARS-CoV-2 viral RNA Reference genomic RNA from WT SARS-CoV-2 virus (WIV04/2019) and variants B.1.351 and B.1.1.7 containing the N501Y mutation was isolated from cell-cultured virus at passage 3 after inoculation of primary patient sample. RNA isolation from 0.2 TCID50 cell-cultured virus was performed by using the MagNA pure 96 (Roche diagnostics) and the total nucleic acid isolation kit (Roche Diagnostics). The sequence of the cell-cultured RNA was confirmed by whole-genome sequencing. ### Virus concentration and nucleic acid extraction Samples were transported to the laboratory and processed as previously described (Medema et al 2020). In short, centrifugation was used to pellet larger particles, viral particles were concentrated from 50 ml supernatant by ultrafiltration through Centricon® Plus-70 centrifugal ultrafilters with a cut-off of 30 or 100 kDa, dependent on supplies (Millipore, Amsterdam, The Netherlands). Approximately 2×104 genomic RNA copies from the murine coronavirus Mouse Hepatitis Virus (MHV)-A59 (obtained from Leiden University Medical Center, The Netherlands) was spiked to each concentrate and co-isolated during the extraction procedure to monitor the possible presence of RT-PCR inhibitors and measure the recovery efficiency of the extraction procedure. Nucleic acid was extracted from the concentrate using the magnetic extraction reagents of the Biomerieux Nuclisens kit (Biomerieux, Amersfoort, the Netherlands) in combination with the semi-automated KingFisher mL (Thermo Scientific, Bleiswijk, The Netherlands) as previously described. Extracted nucleic acid was eluted in a volume of 100 µl. ### RT-qPCR for Coronavirus RNA-quantification in wastewater The N2 assay targeting a fragment of the nucleocapsid gene, as published by US CDC (US-CDC 2020), was used to quantify SARS-CoV-2 RNA in the sewage samples. Reagents and reaction conditions were as previously described (Medema et al. 2020). All RT-PCR’s were run as technical duplicates on 5 µl extracted nucleic acid. RT-qPCR reactions on serial dilutions containing RT-ddPCR calibrated EURM-019 single stranded RNA (provided by the Joint Research Centre) were used to construct calibration curves that subsequently were used to quantify N2 in RNA extracted from the sewage samples. Reactions were considered positive if the cycle threshold was below 40 cycles. Spiked MHV-A59 RNA was detected by performing a MHV-A59 specific RT-qPCR targeting the N-gene using the primers and reaction profile described by Raaben et al. (Raaben, Einerhand et al. 2007). Calibration curves for quantification were generated by performing RT-PCR assays on dilution series of a synthetic quantified gBlock gene fragment containing a partial sequence MHV-A59 N-gene (obtained from IDT, Leuven Belgium). The recovery efficiency was determined by comparing the MHV-A59 RNA concentration in the sewage sample with the concentration MHV-A59 suspension used to spike the samples. This recovery efficiency was used to monitor the performance of the RNA-isolation procedure and the possible presence of RT-PCR inhibitors. The recovery efficiencies are used as process control, but not used to correct the measured SARS-CoV-2 RNA concentrations (Kantor, Nelson et al. 2021). ### PCR for CrAssphage quantification in wastewater A previously described CrAssphage CPQ_064 specific PCR (Stachler, Kelty et al. 2017) was used to quantify this DNA-virus that is ubiquitously present in human intestinal tracts in high concentrations. Assays were performed in duplicate on 5 µl 1:10 diluted extracted nucleic acid. Quantification was performed using PCR assays on dilution series of a synthetic quantified gBlock (obtained from IDT, Leuven, Belgium) containing the CPQ_064 gene fragment. In each wastewater sample, the measured SARS-CoV2 N2-gene concentration was divided by the CrAssphage concentration to normalize for the dilution of the human faecal load by fluctuating human and non-human inputs in the sewer network, such as rain- or groundwater (Medema, Been et al. 2020). ### RT-ddPCR for the N501Y mutation Digital droplet RT-PCR was used to quantify the N501Y and the wild-type (*WIV04/2019*, WT) sequence in one single tube multiplex mutation assay designed by BioRad (Assay ID: dMDS731762551). This assay uses primers that amplify a 80 bp fragment of the Spike gene including the area containing an A to T point mutation that leads to the N501Y amino acid change in the Spike protein (N501Y). Two probes are used to detect PCR-amplification in the droplets: one FAM-labeled probe which perfectly binds to the N501Y mutation and one HEX-labeled probe which perfectly binds to the wild-type SARS-CoV-2 sequence. The ability to perform the PCR-assay in discrete self-contained droplets makes it possible to discriminate between droplets containing SARS-CoV-2 mutant fragments at low frequencies in a background of wild-type fragments. Assays were performed in 20 µl reaction volumes containing the reagents from the One-Step Advance RT-ddPCR for probes: 5 µl RT-ddPCR One-Step Advanced Supermix, 2 µl Reverse Transcriptase, 1 µl DTT (300 mM) supplemented with 1 µl Single tube mutation assay, 6 µl PCR grade and RNAse free water (Applied Biosystems, Fisher Scientific, Landsmeer, The Netherlands) and 5 µl sample-RNA. The BioRad QX200 droplet generator partitioned sample-RNA and reagents in droplets. The temperature profile used for RT-ddPCR was as follows: 60 min. 50°C, 10 min 50°C, 40 cycles with 30 sec. 95°C and 1 min. 55°C followed by 10 min. 98°C, 30 min. 4°C and hold at 12°C. Samples were scanned using the QX200 system (BioRad) and analyzed using the QuantaSoft-Analysis software (BioRad). For each sample, the number of negative and WT or N501Y ddPCR positive droplets were recorded and used to determine the WT or N501Y concentrations. The proportion of Spike gene specific RNA fragments containing the N501Y mutation was calculated by the QuantaSoft-Analysis software as the concentration N501Y in the ddPCR reaction, divided by the sum of WT and N501Y concentrations in the ddPCR reaction. The 95% confidence intervals in the proportion of N501Y were calculated assuming a Poisson distribution of RNA molecules in the droplets. ### Validation experiments Two dilution series were analyzed to evaluate the ability of RT-ddPCR to detect WT and the N501Y variants simultaneously in samples with different concentration ratios of SARS-CoV-2 lineage B.1.351 and WT. The approximate concentration of RNA from WT SARS-CoV-2 virus (Wuhan type) and variant B.1.351 was first quantified using the N2 specific RT-qPCR assay. The first dilution series consisted of a stable concentration of approximately 600 RNA copies of WT, mixed with 2-fold dilutions of lineage B.1.351. The second dilution series contained a stable concentration of approximately 700 RNA copies of variant B.1.351 mixed with 2-fold dilution series of RNA extracted from WT virus. The average concentrations measured with RT-ddPCR in the samples containing stable concentrations of WT or variant B.1.351 respectively were used as values, the dilution factors were subsequently used to calculate the expected proportions of WT and variant B.1.351 in the mixed samples. ## Results ### Method validation To study the ability of RT-ddPCR to differentiate between 501Y and 501N sequences, and to detect low concentrations of SARS-CoV-2 N501Y mutant in the background of WT RNA, two dilutions series were analyzed. For this, we used RNA isolated from cell cultures infected with wild type virus and a B.1.351 strain containing the N501Y mutation. A two-fold dilution series of B.1.351 RNA was made in a stable background of WT RNA (Figure 1) and vice versa (Figure 2), and the concentrations of WT and N501Y was determined by RT-ddPCR. ![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/03/26/2021.03.25.21254324/F1.medium.gif) [Figure 1.](http://medrxiv.org/content/early/2021/03/26/2021.03.25.21254324/F1) Figure 1. Detected copy numbers of WT and N501Y in samples containing ~600 copies WT and 2-fold dilution series of variant lineage B.1.351. ![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/03/26/2021.03.25.21254324/F2.medium.gif) [Figure 2.](http://medrxiv.org/content/early/2021/03/26/2021.03.25.21254324/F2) Figure 2. Detected copy numbers of WT and N501Y in samples containing ~700 copies variant lineage B.1.351WT and 2-fold dilution series of WT. These results demonstrate the ability of the method to simultaneously detect and discriminate between the sequences of WT and the N501Y mutation in lineage B.1.351 from a mixed sample. It also shows that low concentrations of N501Y mutation can be detected in the presence of WT virus RNA and vice versa. Detection of N501Y was possible in a sample dilution containing a theoretical concentration of only 2.7 copies of B.1.351 or WT detection in a sample containing 2.3 copies of WT suggests the feasibility for specific detection of low concentrations of WT and N501Y. Comparing the measured versus expected proportion of the 2-fold dilutions of B.1.351 (Figure 3 A) or the diluted WT (Figure 3B) shows the high linear correlation between the expected proportion and detected proportion of N501Y and WT respectively. ![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/03/26/2021.03.25.21254324/F3.medium.gif) [Figure 3.](http://medrxiv.org/content/early/2021/03/26/2021.03.25.21254324/F3) Figure 3. The expected and detected proportion of N501Y(A) and WT (B) in artificial mixtures of WT and lineage B.1.351 as detected by ddPCR. These results demonstrate that the assay can be used to discriminate between WT SARS-CoV-2 RNA and variants containing the N501Y mutation and simultaneous detection and quantification of these sequences. The specificity of the assay to detect the N501Y mutation was confirmed on RNA from SARS-CoV-2 lineage B.1.1.7 (data not shown). It is expected that the assay can be used to detect low concentrations of WT and N501Y and to obtain reliable insight in the proportion of N501Y in mixtures of viruses, and therefore has the potency to collect information about the spread of variants. ### SARS-CoV-2 VoC surveillance in wastewater of Amsterdam and Utrecht Wastewater samples from the cities of Amsterdam and Utrecht from September 2020 to March 2021 were analyzed to evaluate the applicability of RT-ddPCR to monitor the emergence of VoC’s via wastewater. The concentration SARS-CoV-2 N2 RNA in these samples was first determined using RT-qPCR and compared to newly reported COVID-19 cases in these cities. #### SARS-CoV-2 sewage surveillance in Amsterdam and Utrecht The concentrations of SARS-CoV-2 N2 (normalized for the concentration of CrAssphage in the same sample) in Amsterdam and Utrecht wastewater and the incidence of newly reported cases in Amsterdam and Utrecht are depicted in Figure 4. The results show that the concentrations of SARS-CoV-2 RNA in wastewater reflected the trends in the newly reported positive COVID-19 clinical tests in the population of these cities. Both cities experienced a peak in COVID-19 cases in October and December 2020. ![Figure 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/03/26/2021.03.25.21254324/F4.medium.gif) [Figure 4.](http://medrxiv.org/content/early/2021/03/26/2021.03.25.21254324/F4) Figure 4. The concentration of SARS-CoV-2 N2 gene (normalized for the concentration of Crass-phage) in wastewater of Amsterdam (A) and Utrecht (B) from Aug 26 2020 to March 8 2021 (orange circles) and the newly reported COVID-19 cases in the city population (blue dots) and the 7-day moving average (blue line) COVID-19 data: National Institute of Public Health and the Environment, data source ESRI NL COVID-19 Hub (Esri-Nederland 2021) #### Emergence of the SARS-CoV-2 N501Ymutation in sewage from Amsterdam and Utrecht Analysis of the sewage samples was performed to study the applicability of the N501Y RT-ddPCR assay. The results are summarized in table S1 (Amsterdam) and table S2 (Utrecht). The proportion of N501Y to WT sequences in Amsterdam and Utrecht wastewater was compared with publicly available data from the Dutch National Institute for Public Health and the Environment (RIVM). The available data is based on surveillance of variants by sequencing a random selection of SARS-COV-2 isolates from several hundreds of new cases in the Netherlands every week (Table S3, “Pathogen Genomic surveillance”: [https://www.rivm.nl/coronavirus-covid-19/virus/varianten](https://www.rivm.nl/coronavirus-covid-19/virus/varianten), accessed March 09 2021). The results of these comparisons are shown in Figures 5A (Amsterdam) and 5B (Utrecht). ![Figure 5.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/03/26/2021.03.25.21254324/F5.medium.gif) [Figure 5.](http://medrxiv.org/content/early/2021/03/26/2021.03.25.21254324/F5) Figure 5. Proportion of Spike gene RNA fragments containing the N501Y mutation in wastewater from Amsterdam (A) and Utrecht (B) calculated as relative concentration of N501Ycontaining Spike gene fragments divided by the total concentration of S-gene fragments (WT+ N501Y). Error bars are calculated assuming Poisson distribution of RNA molecules in droplets. The proportion of newly COVID-19 patients infected with N501Y containing variants (Pathogen genomic surveillance) are shown as (the proportional sum of lineages B.1.1.7, B.1.351 and P.1). The earliest detection of N501Y was in the wastewater sample of December 8 2020 in Amsterdam and in the wastewater sample of December 21 2020 in Utrecht. The proportion of N501Y gradually increased to 85% on February 22 2021 in Amsterdam and 81% on February 22 2021 in Utrecht. The total SARS-CoV-2 RNA concentrations were higher in wastewater from Amsterdam compared to Utrecht (Table S1 and S2), yielding a more accurate assessment of the N501Y proportion, as indicated by the error-bars. The first detection of the British variant (lineage B.1.1.7) in patients in the Netherlands was in the first week of December 2020. This corresponds to the first detection of N501Y in wastewater from Amsterdam. The first detection of N501Y in wastewater from Utrecht was two weeks later. After the emergence of variants containing the N501Y mutation, the increase in the proportion of this marker in wastewater follows a trend that is comparable with the increase of N501Y containing variants in patients in the Netherlands. The percentage of N501Y containing variants appears to be higher in Amsterdam wastewater in the complete period suggesting earlier introduction of the British variant in Amsterdam than in Utrecht. The introduction and emergence of SARS-CoV-2 variants containing the N501Y mutation in Utrecht wastewater aligns more in time with the national ‘average’ emergence as obtained through the genomic pathogen surveillance data. #### Comparison RT-qPCR with RT-ddPCR The concentrations SARS-CoV-2 RNA measured in sewage from WWTP Amsterdam and Utrecht with RT-qPCR (N2 amplicon of the N-gene) and RT-ddPCR (sum of WT and N501Y sequences of the Spike gene) are compared in Figure 6. This comparison shows that there is a linear relation between the values measured with both methods. However, concentrations measured with RT-qPCR are consistently higher (61% on average) than concentrations measured with RT-ddPCR. On basis of current knowledge of SARS-CoV-2 sequences in the Dutch population, it is unlikely that there are other sequence variants in the populations that were detected by RT-qPCR but missed with the RT-ddPCR assay. The higher RT-qPCR results could be the result of targeting different gene fragments with potential differences in detected concentrations in wastewater. The different RNA-targets can have different degradation rates in sewage but it can also be the result of higher expression levels of the N-gene compared to the S-gene in coronavirus infected cells and the presence of the N-gene in every subgenomic mRNA (Kim, Lee et al. 2020) making it likely that higher levels of N-sequences will be present in infected intestinal cells and eventually may lead to higher N-sequence concentrations in sewage. ![Figure 6.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/03/26/2021.03.25.21254324/F6.medium.gif) [Figure 6.](http://medrxiv.org/content/early/2021/03/26/2021.03.25.21254324/F6) Figure 6. Comparison of SARS-CoV-2 RNA, measured with RT-qPCR and RT-ddPCR, concentrations in sewage. ## Discussion and conclusions This research describes the first use of RT-ddPCR to specifically detect and quantify SARS-CoV-2 mutations associated with VoC’s in wastewater. This assay can be used for sensitive and simultaneous detection of WT and the N501Y mutation in sewage. Although more extensive research is needed to determine the limit of detection accurately, experiments on SARS-CoV-2 RNA from lineage B.1.351 (with mutation N501Y) and WT in different proportional mixtures demonstrated that proportions as low as 0,5% of lineage B.1.351 can be detected in a background of 700 RNA-copies of WT SARS-CoV-2. These experiments also demonstrated the linearity of the expected proportion and the measured proportion of B.1.351 in the mixtures and vice versa, indicating that this N501Y assay can be used to accurately determine the ratio between SARS-CoV-2 genomes containing N501Y and those containing the WT sequence. The weekly analysis of wastewater from the cities of Amsterdam and Utrecht using RT-qPCR targeting the N-gene (CDC N2 target) clearly demonstrated the relationship between SARS-CoV-2 RNA levels and the Covid-19 incidence in the communities connected to the respective WWTP’s. N501Y containing variants were first detected on December 8 2020 in Amsterdam wastewater and two weeks later in Utrecht wastewater. The proportion of N501Y variants in the weekly samples increased to approximately 80% on March 1, 2021 in both cities. These results demonstrate the applicability of this RT-ddPCR assay to monitor the absolute quantities and ratios of WT viral sequences and sequences containing the N501Y mutation. This makes the use of RT-ddPCR on wastewater samples an attractive, rapid and efficient method to follow the emergence of (mutations associated with) VoC’s in the community, that can generate an early warning of the emergence of VoC’s on the basis of one “population-sample” via domestic wastewater. This proof-of-principle demonstrates the value of RT-ddPCR to detect N501Y, and suggests RT-ddPCR can also be used for the detection of multiple mutations (like f.e. E484K, K417 or P681H) characteristic for various epidemiologically or immunologically relevant VoC’s in sewage to monitor the emergence and spread of VoC in the community. Whole genome sequencing of SARS-CoV-2 from clinical samples from infected individuals is commonly used to identify VoC’s (Tegally, Wilkinson et al. 2021). Although whole genome sequencing has shown to be feasible to give insight in the sequence variations of SARS-CoV-2 genomes in sewage (Izquierdo-Lara, Elsinga et al. 2020) and has been used to detect mutations associated with variants of concern in wastewater(Jahn, Dreifuss et al. 2021), this requires deep sequencing and bioinformatics to determine the presence of mutations associated with VoC’s in wastewater, and is provides less quantitative results than ddPCR. An advantage of whole genome sequencing of wastewater samples is that it can provide more extensive information about the co-occurrence of the range of mutations and deletions associated with VoC. It would be valuable to compare the results of RT-ddPCR with the N501Y (and other) mutations with those of whole genome sequencing on the same wastewater samples and on samples of the COVID-19 cases in the population that contributes to this wastewater. ## Data Availability All data details described in this paper is available and can be provided upon request ## Supplemental material View this table: [Table S1.](http://medrxiv.org/content/early/2021/03/26/2021.03.25.21254324/T1) Table S1. Summary of the RT-ddPCR N501Y and WT results in wastewater samples from WWTP Amsterdam View this table: [Table S2.](http://medrxiv.org/content/early/2021/03/26/2021.03.25.21254324/T2) Table S2. Summary of the RT-ddPCR N501Y and WT results in wastewater samples from WWTP Utrecht View this table: [Table S3.](http://medrxiv.org/content/early/2021/03/26/2021.03.25.21254324/T3) Table S3. N501Y containing SARS-CoV-2 isolates based on sequence analysis of isolates from patients (“ Pathogen Genomic surveillance”: [https://www.rivm.nl/coronavirus-covid-19/virus/varianten](https://www.rivm.nl/coronavirus-covid-19/virus/varianten), accessed March 09 2021) ## Acknowledgements The authors are very grateful for the assistance of the Water Authorities and WWTP operators of Waternet and Hoogheemraadschap de Stichtse Rijnlanden in the Netherlands who organized the sampling and provided the sewage samples. We are also grateful to Meindert de Graaf for sample transport. Eddy van Collenburg (BIORAD) is acknowledged for technical advice and providing the test kits. Coronavirus MHV-A59 was a kind gift from Prof. Eric Snijder and provided by Linda Boomaars (Leids Universitair Medisch Centrum). This study was financed by TKI Health Holland in EUREKA project WASTEWATER4COVID. ## Footnotes * Competing financial interests The authors declare they have no actual or potential competing financial interests. * Received March 25, 2021. * Revision received March 25, 2021. * Accepted March 26, 2021. * © 2021, Posted by Cold Spring Harbor Laboratory The copyright holder for this pre-print is the author. All rights reserved. The material may not be redistributed, re-used or adapted without the author's permission. ## References 1. Ahmed, W., B. Tscharke, P. M. Bertsch, K. Bibby, A. Bivins, P. Choi, L. Clarke, J. Dwyer, J. Edson, T. M. H. Nguyen, J. W. O’Brien, S. L. Simpson, P. Sherman, K. V. Thomas, R. Verhagen, J. Zaugg and J. F. Mueller (2021). “SARS-CoV-2 RNA monitoring in wastewater as a potential early warning system for COVID-19 transmission in the community: A temporal case study.” Sci Total Environ 761: 144216. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.scitotenv.2020.144216&link_type=DOI) 2. Altmann, D. M., R. J. Boyton and R. Beale (2021). “Immunity to SARS-CoV-2 variants of concern.” Science 371(6534): 1103–1104. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEzOiIzNzEvNjUzNC8xMTAzIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDMvMjYvMjAyMS4wMy4yNS4yMTI1NDMyNC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 3. Betancourt, W. W., B. W. Schmitz, G. K. Innes, K. M. Pogreba Brown, S. M. Prasek, E. R. Stark, A. R. Foster, R. S. Sprissler, D. T. Harris, S. P. Sherchan, C. P. Gerba and I. L. Pepper (2020). “Wastewater-based Epidemiology for Averting COVID-19 Outbreaks on The University of Arizona Campus.” medRxiv: 2020.2011.2013.20231340. 4. Bivins, A., D. North, A. Ahmad, W. Ahmed, E. Alm, F. Been, P. Bhattacharya, L. Bijlsma, A. B. Boehm, J. Brown, G. Buttiglieri, V. Calabro, A. Carducci, S. Castiglioni, Z. Cetecioglu Gurol, S. Chakraborty, F. Costa, S. Curcio, F. L. de Los Reyes, 3rd., J. Delgado Vela, K. Farkas, X. Fernandez-Casi, C. Gerba, D. Gerrity, R. Girones, R. Gonzalez, E. Haramoto, A. Harris, P. A. Holden, M. T. Islam, D. L. Jones, B. Kasprzyk-Hordern, M. Kitajima, N. Kotlarz, M. Kumar, K. Kuroda, G. La Rosa, F. Malpei, M. Mautus, S. L. McLellan, G. Medema, J. S. Meschke, J. Mueller, R. J. Newton, D. Nilsson, R. T. Noble, A. van Nuijs, J. Peccia, T. A. Perkins, A. J. Pickering, J. Rose, G. Sanchez, A. Smith, L. Stadler, C. Stauber, K. Thomas, T. van der Voorn, K. Wigginton, K. Zhu and K. Bibby (2020). “Wastewater-Based Epidemiology: Global Collaborative to Maximize Contributions in the Fight Against COVID-19.” Environ Sci Technol 54(13): 7754–7757. 5. Challen, R., E. Brooks-Pollock, J. M. Read, L. Dyson, K. Tsaneva-Atanasova and L. Danon (2021). “Risk of mortality in patients infected with SARS-CoV-2 variant of concern 202012/1: matched cohort study.” BMJ 372: 579. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1136/bmj.n579&link_type=DOI) 6. Davies, N. G., S. Abbott, R. C. Barnard, C. I. Jarvis, A. J. Kucharski, J. D. Munday, C. A. B. Pearson, T. W. Russell, D. C. Tully, A. D. Washburne, T. Wenseleers, A. Gimma, W. Waites, K. L. M. Wong, K. van Zandvoort, J. D. Silverman, K. Diaz-Ordaz, R. Keogh, R. M. Eggo, S. Funk, M. Jit, K. E. Atkins and W. J. Edmunds (2021). “Estimated transmissibility and impact of SARS-CoV-2 lineage B.1.1.7 in England.” Science: eabg3055. 7. Esri-Nederland (2021). “Website: COVID-19 - Historische gegevens RIVM (accessed March 2021).” [https://nlcovid-19-esrinl-content.hub.arcgis.com/](https://nlcovid-19-esrinl-content.hub.arcgis.com/). 8. Eurosurveillance-editorial-team (2021). “Updated rapid risk assessment from ECDC on the risk related to the spread of new SARS-CoV-2 variants of concern in the EU/EEA – first update.” Eurosurveillance 26(3): 2101211. 9. Faria, N. R., T. A. Mellan, C. Whittaker, I. M. Claro, D. D. S. Candido, S. Mishra, M. A. E. Crispim, F. C. Sales, I. Hawryluk, J. T. McCrone, R. J. G. Hulswit, L. A. M. Franco, M. S. Ramundo, J. G. de Jesus, P. S. Andrade, T. M. Coletti, G. M. Ferreira, C. A. M. Silva, E. R. Manuli, R. H. M. Pereira, P. S. Peixoto, M. U. Kraemer, N. Gaburo, C. D. C. Camilo, H. Hoeltgebaum, W. M. Souza, E. C. Rocha, L. M. de Souza, M. C. de Pinho, L. J. T. Araujo, F. S. V. Malta, A. B. de Lima, J. D. P. Silva, D. A. G. Zauli, A. C. S. Ferreira, R. P. Schnekenberg, D. J. Laydon, P. G. T. Walker, H. M. Schlüter, A. L. P. Dos Santos, M. S. Vidal, V. S. Del Caro, R. M. F. Filho, H. M. Dos Santos, R. S. Aguiar, J. L. P. Modena, B. Nelson, J. A. Hay, M. Monod, X. Miscouridou, H. Coupland, R. Sonabend, M. Vollmer, A. Gandy, M. A. Suchard, T. A. Bowden, S. L. K. Pond, C. H. Wu, O. Ratmann, N. M. Ferguson, C. Dye, N. J. Loman, P. Lemey, A. Rambaut, N. A. Fraiji, M. Carvalho, O. G. Pybus, S. Flaxman, S. Bhatt and E. C. Sabino (2021). “Genomics and epidemiology of a novel SARS-CoV-2 lineage in Manaus, Brazil.” medRxiv. 10. Garcia-Beltran, W. F., E. C. Lam, K. S. Denis, A. D. Nitido, Z. H. Garcia, B. M. Hauser, J. Feldman, M. N. Pavlovic, D. J. Gregory, M. C. Poznansky, A. Sigal, A. G. Schmidt, A. J. Iafrate, V. Naranbhai and A. B. Balazs (2021). “Circulating SARS-CoV-2 variants escape neutralization by vaccine-induced humoral immunity.” medRxiv. 11. Izquierdo-Lara, R., G. Elsinga, L. Heijnen, B. B. Oude Munnink, C. M. E. Schapendonk, D. Nieuwenhuijse, M. Kon, L. Lu, F. M. Aarestrup, S. Lycett, G. Medema, M. P. G. Koopmans and M. de Graaf (2020). “Monitoring SARS-CoV-2 circulation and diversity through community wastewater sequencing.” medRxiv: 2020.2009.2021.20198838. 12. Jahn, K., D. Dreifuss, I. Topolsky, A. Kull, P. Ganesanandamoorthy, X. Fernandez-Cassi, C. Bänziger, E. Stachler, L. Fuhrmann, K. P. Jablonski, C. Chen, C. Aquino, T. Stadler, C. Ort, T. Kohn, T. R. Julian and N. Beerenwinkel (2021). “Detection of SARS-CoV-2 variants in Switzerland by genomic analysis of wastewater samples.” medRxiv: 2021.2001.2008.21249379. 13. Jones, D. L., M. Q. Baluja, D. W. Graham, A. Corbishley, J. E. McDonald, S. K. Malham, L. S. Hillary, T. R. Connor, W. H. Gaze, I. B. Moura, M. H. Wilcox and K. Farkas (2020). “Shedding of SARS-CoV-2 in feces and urine and its potential role in person-to-person transmission and the environment-based spread of COVID-19.” Science of The Total Environment 749: 141364. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.scitotenv.2020.141364&link_type=DOI) 14. Kantor, R. S., K. L. Nelson, H. D. Greenwald and L. C. Kennedy (2021). “Challenges in Measuring the Recovery of SARS-CoV-2 from Wastewater.” Environ Sci Technol. 15. Kim, D., J.-Y. Lee, J.-S. Yang, J. W. Kim, V. N. Kim and H. Chang (2020). “The Architecture of SARS-CoV-2 Transcriptome.” Cell 181(4): 914–921.e910. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F26%2F2021.03.25.21254324.atom) 16. Liu, H., Q. Zhang, P. Wei, Z. Chen, K. Aviszus, J. Yang, W. Downing, S. Peterson, C. Jiang, B. Liang, L. Reynoso, G. P. Downey, S. K. Frankel, J. Kappler, P. Marrack and G. Zhang (2021). “The basis of a more contagious 501Y.V1 variant of SARS-COV-2.” bioRxiv. 17. Lodder, W. and A. M. de Roda Husman (2020). “SARS-CoV-2 in wastewater: potential health risk, but also data source.” Lancet Gastroenterol Hepatol 5(6): 533–534. 18. Medema, G., F. Been, L. Heijnen and S. Petterson (2020). “Implementation of environmental surveillance for SARS-CoV-2 virus to support public health decisions: Opportunities and challenges.” Curr Opin Environ Sci Health 17: 49–71. 19. Medema, G., L. Heijnen, G. Elsinga, R. Italiaander and A. Brouwer (2020). “Presence of SARS-Coronavirus-2 in sewage.” medRxiv: 2020.2003.2029.20045880. 20. Medema, G., L. Heijnen, G. Elsinga, R. Italiaander and A. Brouwer (2020). “Presence of SARS-Coronavirus-2 RNA in Sewage and Correlation with Reported COVID-19 Prevalence in the Early Stage of the Epidemic in The Netherlands.” Environmental Science & Technology Letters 7(7): 511–516. 21. Pekin, D., Y. Skhiri, J. C. Baret, D. Le Corre, L. Mazutis, C. B. Salem, F. Millot, A. El Harrak, J. B. Hutchison, J. W. Larson, D. R. Link, P. Laurent-Puig, A. D. Griffiths and V. Taly (2011). “Quantitative and sensitive detection of rare mutations using droplet-based microfluidics.” Lab Chip 11(13): 2156–2166. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1039/c1lc20128j&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21594292&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F26%2F2021.03.25.21254324.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000291555500005&link_type=ISI) 22. Prado, T., T. M. Fumian, C. F. Mannarino, P. C. Resende, F. C. Motta, A. L. F. Eppinghaus, V. H. Chagas do Vale, R. M. S. Braz, J. d. S. R. de Andrade, A.G. Maranhão and M. P. Miagostovich (2021). “Wastewater-based epidemiology as a useful tool to track SARS-CoV-2 and support public health policies at municipal level in Brazil.” Water research 191: 116810–116810. 23. Raaben, M., A. W. C. Einerhand, L. J. A. Taminiau, M. van Houdt, J. Bouma, R. H. Raatgeep, H. A. Büller, C. A. M. de Haan and J. W. A. Rossen (2007). “Cyclooxygenase activity is important for efficient replication of mouse hepatitis virus at an early stage of infection.” Virology journal 4: 55–55. 24. Rynkiewicz, P., G. A. Babbitt, F. Cui, A. O. Hudson and M. L. Lynch (2021). “A comparative survey of Betacoronavirus binding dynamics relevant to the functional evolution of the highly transmissible SARS-CoV-2 variant N501Y.” bioRxiv. 25. Stachler, E., C. Kelty, M. Sivaganesan, X. Li, K. Bibby and O. C. Shanks (2017). “Quantitative CrAssphage PCR Assays for Human Fecal Pollution Measurement.” Environ Sci Technol 51(16): 9146–9154. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1021/acs.est.7b02703&link_type=DOI) 26. Tegally, H., E. Wilkinson, M. Giovanetti, A. Iranzadeh, V. Fonseca, J. Giandhari, D. Doolabh, S. Pillay, E. J. San, N. Msomi, K. Mlisana, A. von Gottberg, S. Walaza, M. Allam, A. Ismail, T. Mohale, A. J. Glass, S. Engelbrecht, G. Van Zyl, W. Preiser, F. Petruccione, A. Sigal, D. Hardie, G. Marais, M. Hsiao, S. Korsman, M. A. Davies, L. Tyers, I. Mudau, D. York, C. Maslo, D. Goedhals, S. Abrahams, O. Laguda-Akingba, A. Alisoltani-Dehkordi, A. Godzik, C. K. Wibmer, B. T. Sewell, J. Lourenço, L. C. J. Alcantara, S. L. Kosakovsky Pond, S. Weaver, D. Martin, R. J. Lessells, J. N. Bhiman, C. Williamson and T. de Oliveira (2021). “Emergence of a SARS-CoV-2 variant of concern with mutations in spike glycoprotein.” Nature. 27. US-CDC (2020). “2019-Novel Coronavirus (2019-nCoV) Real-Time rRT-PCR Panel Primers and Probes.” [https://www.cdc.gov/coronavirus/2019-ncov/lab/rt-pcr-panel-primer-probes.html](https://www.cdc.gov/coronavirus/2019-ncov/lab/rt-pcr-panel-primer-probes.html) (accessed May 19, 2020). 28. Volz, E., S. Mishra, M. Chand, J. C. Barrett, R. Johnson, L. Geidelberg, W. R. Hinsley, D. J. Laydon, G. Dabrera, Á. O’Toole, R. Amato, M. Ragonnet-Cronin, I. Harrison, B. Jackson, C. V. Ariani, O. Boyd, N. J. Loman, J. T. McCrone, S. Gonçalves, D. Jorgensen, R. Myers, V. Hill, D. K. Jackson, K. Gaythorpe, N. Groves, J. Sillitoe, D. P. Kwiatkowski, S. Flaxman, O. Ratmann, S. Bhatt, S. Hopkins, A. Gandy, A. Rambaut and N. M. Ferguson (2021). “Transmission of SARS-CoV-2 Lineage B.1.1.7 in England: Insights from linking epidemiological and genetic data.” medRxiv: 2020.2012.2030.20249034. 29. Wang, D., B. Hu, C. Hu, F. Zhu, X. Liu, J. Zhang, B. Wang, H. Xiang, Z. Cheng, Y. Xiong, Y. Zhao, Y. Li, X. Wang and Z. Peng (2020). “Clinical Characteristics of 138 Hospitalized Patients With 2019 Novel Coronavirus–Infected Pneumonia in Wuhan, China.” JAMA 323(11): 1061–1069. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jama.2020.1585&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F26%2F2021.03.25.21254324.atom) 30. Wang, P., L. Liu, S. Iketani, Y. Luo, Y. Guo, M. Wang, J. Yu, B. Zhang, P. D. Kwong, B. S. Graham, J. R. Mascola, J. Y. Chang, M. T. Yin, M. Sobieszczyk, C. A. Kyratsous, L. Shapiro, Z. Sheng, M. S. Nair, Y. Huang and D. D. Ho (2021). “Increased Resistance of SARS-CoV-2 Variants B.1.351 and B.1.1.7 to Antibody Neutralization.” bioRxiv. 31. Wölfel, R., V. M. Corman, W. Guggemos, M. Seilmaier, S. Zange, M. A. Müller, D. Niemeyer, T. C. Jones, P. Vollmar, C. Rothe, M. Hoelscher, T. Bleicker, S. Brünink, J. Schneider, R. Ehmann, K. Zwirglmaier, C. Drosten and C. Wendtner (2020). “Virological assessment of hospitalized patients with COVID-2019.” Nature 581(7809): 465–469. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41586-020-2196-x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F26%2F2021.03.25.21254324.atom) 32. Xiao, F., M. Tang, X. Zheng, Y. Liu, X. Li and H. Shan (2020). “Evidence for Gastrointestinal Infection of SARS-CoV-2.” Gastroenterology 158(6): 1831–1833.e1833. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1053/j.gastro.2020.02.055&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F26%2F2021.03.25.21254324.atom)