key: cord-0759291-uy0z62s6
authors: Carcereny, Albert; Martínez-Velázquez, Adán; Bosch, Albert; Allende, Ana; Truchado, Pilar; Cascales, Jenifer; Romalde, Jesús L; Lois, Marta; Polo, David; Sánchez, Gloria; Pérez-Cataluña, Alba; Díaz-Reolid, Azahara; Antón, Andrés; Gregori, Josep; Garcia-Cehic, Damir; Quer, Josep; Palau, Margarita; Ruano, Cristina González; Pintó, Rosa M; Guix, Susana
title: Monitoring Emergence of the SARS-CoV-2 B.1.1.7 Variant through the Spanish National SARS-CoV-2 Wastewater Surveillance System (VATar COVID-19)
date: 2021-08-16
journal: Environ Sci Technol
DOI: 10.1021/acs.est.1c03589
sha: 7d50954daf20e9bdd4b19b8e88b0b311d1ca9d8d
doc_id: 759291
cord_uid: uy0z62s6

[Image: see text] Since its first identification in the United Kingdom in late 2020, the highly transmissible B.1.1.7 variant of SARS-CoV-2 has become dominant in several countries raising great concern. We developed a duplex real-time RT-qPCR assay to detect, discriminate, and quantitate SARS-CoV-2 variants containing one of its mutation signatures, the ΔHV69/70 deletion, and used it to trace the community circulation of the B.1.1.7 variant in Spain through the Spanish National SARS-CoV-2 Wastewater Surveillance System (VATar COVID-19). The B.1.1.7 variant was detected earlier than clinical epidemiological reporting by the local authorities, first in the southern city of Málaga (Andalucía) in week 20_52 (year_week), and multiple introductions during Christmas holidays were inferred in different parts of the country. Wastewater-based B.1.1.7 tracking showed a good correlation with clinical data and provided information at the local level. Data from wastewater treatment plants, which reached B.1.1.7 prevalences higher than 90% for ≥2 consecutive weeks showed that 8.1 ± 2.0 weeks were required for B.1.1.7 to become dominant. The study highlights the applicability of RT-qPCR-based strategies to track specific mutations of variants of concern as soon as they are identified by clinical sequencing and their integration into existing wastewater surveillance programs, as a cost-effective approach to complement clinical testing during the COVID-19 pandemic.

Environmental surveillance of specimens contaminated by human feces is used to monitor enteric virus disease transmission in the population, and several countries have implemented SARS-CoV-2 wastewater monitoring networks to aid decision making during the COVID-19 pandemic. 1−3 In Spain, a nationwide COVID-19 wastewater surveillance project (VATar COVID-19) was launched in June 2020 (https:// www.miteco.gob.es/es/agua/temas/concesiones-y-autorizaciones/vertidos-de-aguas-residuales/alerta-temprana-covid19/ default.aspx) and has weekly monitored SARS-CoV-2 levels in untreated wastewater from initially 32 wastewater treatment plants (WWTPs) since then. On March 2021, the European Commission adopted a recommendation on a common approach to establish and make greater use of systematic wastewater surveillance of SARS-CoV-2 as a new source of independent information on the spread of the virus and its variants in the European Union. 4 In situations with low or absent SARS-CoV-2 circulation in the community, wastewater surveillance has proven to be a useful tool as an early warning system, 5−9 and several studies have also tried to infer disease incidence in a community, independent of diagnostic testing availability based on SARS-CoV-2 wastewater concentrations, with considerable uncertainties. 10−12 Despite titanic efforts based on confinement measures and mass-vaccination programs, the emergence of novel variants of concern (VOCs), mainly B.1.1.7, B.1.351, B.1.1.28.1, and recently B1.617.2, so far, suggests that continued surveillance is required to control the COVID-19 pandemic in the long run. Since January 2021, countries within and outside Europe have observed a substantial increase in the number and proportion of SARS-CoV-2 cases of the B.1.1.7 variant, first reported in the United Kingdom. 13, 14 Since the B.1.1.7 variant has been shown to be more transmissible than the previously predominant circulating variants and since its infections may be more severe, 15 countries where the variant has spread and become dominant are concerned on whether the occurrence of the variant will result in increases in total COVID-19 incidences, hospitalizations, and excess mortality due to overstretched health systems.

The emergence of SARS-CoV-2 variants that may increase transmissibility and/or immune escape points to an imperative need for the implementation of targeted surveillance methods. While sequencing should be the gold standard for variant characterization, cost-effective molecular assays, which could be rapidly established and scaled up, may offer several advantages and provide valuable quantitative information without delay.

This study included the development and validation of a onetube duplex quantitative real-time RT-PCR (RT-qPCR) assay to detect, discriminate, and quantitate SARS-CoV-2 variants containing the ΔHV69/70 deletion from variants lacking it, using allelic discrimination probes. Confirmatory sequencing of a subset of samples was performed to be able to ascertain the validity of these assays to trace the community circulation of the B.1.1.7 variant. The RT-qPCR-based assay improved the current variant tracking capability and could be easily implemented for monitoring the emergence of ΔHV69/70-containing SARS-CoV-2 variants (mainly B.1.1.7) in Spain through the nationwide wastewater surveillance network.

■ METHODS Wastewater Sampling. Influent water grab samples were weekly collected from 32 WWTPs (one weekly sample per site) located in 15 different autonomous communities in Spain, from the middle of December 2020 (week 20_51, year_week number) to the end of March 2021 (week 21_13) (last week of 2020 was not sampled). All samples were transported on ice to one of the four participating laboratories of analysis (A, B, C, and D), stored at 4°C and processed within 1−2 days upon arrival. The time between sample collection and arrival to the laboratory ranged between 3−24 h.

Sample Concentration, Nucleic Acid Extraction, and Process Control. Wastewater samples (200 mL) were concentrated by the aluminum hydroxide adsorption-precipitation method, as previously described. 7, 16 Briefly, samples were adjusted to pH 6.0, a 1:100 v:v of 0.9 N AlCl 3 solution was added, and samples were gently mixed for 15 min at room temperature. The precipitate was collected by centrifugation at 1700× g for 20 min, and the pellet was resuspended in 10 mL of 3% beef extract (pH 7.4). After a 10 min shaking at 150 rpm, samples were centrifuged at 1900× g for 30 min, and concentrates were resuspended in 1−2 mL of phosphate buffered saline (PBS). All samples were spiked with a known amount of an animal coronavirus used as a process control virus. Animal coronaviruses differed between participant laboratories and included the attenuated PUR46-MAD strain of transmissible gastroenteritis enteric virus (TGEV), 17 porcine epidemic diarrhea virus (PEDV) strain CV777 (kindly provided by Prof. A. Carvajal from the University of Leon), and murine hepatitis virus (MHV) strain ATCC VR-764. 18 Depending on each laboratory, between 10 and 100 μL of the animal coronavirus stock were added to 200 mL of sample to obtain final concentrations of 2.5 × 10 4 −2.5 × 10 5 copies/mL (PEDV and MHV) or 6.9 × 10 3 TCID50/mL (TGEV). Nucleic acid extraction from concentrates was performed from 300 μL of sample using the Maxwell RSC PureFood GMO and Authentication Kit (Promega Corporation, Madison, US) or from 150 μL of sample using the NucleoSpin RNA Virus kit (Macherey-Nagel GmbH & Co., Duren, Germany), following the manufacturer's instructions. Each extraction included a negative control and a process virus control used to estimate the virus recovery efficiency. RT-qPCR for process control viruses was performed as previously described. 19−21 Parallel to ISO 15216-1:2017 22 for the determination of norovirus and hepatitis A virus in the food chain, samples with a virus recovery ≥1% were considered acceptable.

SARS-CoV-2 RT-qPCR Assays. The N1 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, using the PrimeScript One-Step RT-PCR Kit (Takara Bio, USA) and 2019-nCoV RUO qPCR Probe Assay primer/probe mix (IDT, Integrated DNA Technologies, Leuven, Belgium). Different instruments were used by different participating labs, including CFX96 BioRad, LightCycler 480 (Roche Diagnostics, Germany), Stratagene Mx3005P (Applied Biosystems, USA), and QuantStudio 5 (Applied Biosystems, USA).

The S gene was analyzed by a duplex gene allelic discrimination TaqMan RT-qPCR assay, using 400 nM of the following primers targeting the S gene (For-S21708 5′ATT-CAACTCAGGACTTGTTCTTACCTT3′ and Rev-S21796 5′TAAATGGTAGGACAGGGTTATCAAAC3′) and 200 nM of the following probes (S_Probe6970in 5′FAM-TCCATGC-TATACATGTCTCTGGGACCAATG BHQ1−3′ and S_Pro-be6970del 5′HEX-TTCCATGCTATCTCTGGGAC-CAATGGTACT BHQ1−3′). The RT-qPCR mastermix was prepared using the PrimeScript One-Step RT-PCR Kit (Takara Bio, USA), and the temperature program was 10 min at 50°C, 3 min at 95°C, and 45 cycles of 3 s at 95°C and 30 s at 60°C.

RT-qPCR analysis for each target included the analysis of duplicate wells containing undiluted RNA and duplicate wells containing a 10-fold dilution to monitor the presence of inhibitors. Every RT-qPCR assay included four wells corresponding to negative controls (two nuclease-free water and two negative extraction controls). Commercially available Twist Synthetic SARS-CoV-2 RNA Controls (Control 2, MN908947.3 and Control 14, EPI_ISL_710528) were used to prepare standard curves for genome quantitation. Both synthetic RNA controls were quantified by droplet-based digital PCR using the One-Step RT-ddPCR Advanced Kit for probes in a QX200 System (BioRad), to estimate the exact concentration of genome copies (GC)/μl, prior to construction of RT-qPCR standard curves. The limit of detection (LOD) and limit of quantification (LOQ) were determined for each specific target by running a series of dilutions of the target with 4−10 replicates per dilution. Parameters of all standard curves and estimated LOD and LOQ for the four participating laboratories are summarized in Supporting Information Table S1 .

RT-qPCR Data Analysis and Interpretation. The following criteria were used to estimate SARS-CoV-2 gene viral titers. For each specific target, Cq values ≤ 40 were converted into GC/L using the corresponding standard curve and volumes tested. Occurrence of inhibition was estimated by comparing average viral titers obtained from duplicate wells tested on undiluted RNA with duplicate wells tested on 10-fold diluted RNA. Inhibition was ascertained when difference in Environmental Science & Technology pubs.acs.org/est Article average viral titers was higher than 0.5 log 10 , and if that occurred, viral titers were inferred from the 10-fold RNA dilution. The percentage of SARS-CoV-2 genomes containing the ΔHV69/70 deletion within the S gene was calculated using the following formula:

In cases with one of the concentrations <LOQ, the percentage was calculated using the corresponding LOQ of the assay. Data with both concentrations <LOQ were not considered.

S Gene Sequencing and Single Nucleotide Polymorphism (SNP) Identification. Full-length S gene sequencing of wastewater samples was performed following the ARTIC Network protocol (https://artic.network/ncov-2019, with minor modifications) with selected v3 primers (Integrated DNA Technologies) for genome amplification and KAPA HyperPrep Kit (Roche Applied Science) for library preparation. 23 Libraries were loaded in MiSeq Reagent Kit 600v3 cartridges and sequenced on a MiSeq platform (Illumina). The raw sequenced reads were cleaned from low-quality segments and mapped against the Wuhan-Hu-1 reference genome sequence to find out variant-specific signature mutations (mutations and indels).

Duplex SARS-CoV-2 S Gene Allelic Discrimination RT-qPCR Assay Validation. To confirm the ability of the RT-qPCR assay to discriminate and estimate the proportion of targets containing the ΔHV69/70 deletion in different scenarios, nine preparations containing 90:10, 50:50, or 10:90 proportions of B.1.1.7 and Wuhan-Hu-1 synthetic control RNAs at three different total concentration levels (1 × 10 4 , 1 × 10 3 , and 1 × 10 2 GC/rxn) were made and analyzed ( Figure 1 ). Assays performed using B.1.1.7 and Wuhan-Hu-1 synthetic control RNAs did not show cross-reactivity (data not shown). Results showed the ability of the method to detect and quantitatively discriminate sequences containing the ΔHV69/70 deletion from wildtype sequences in mixed samples through a wide range of total RNA concentrations often observed in extracted wastewater samples.

Additionally, to confirm that sequences detected in natural samples collected during the study period containing the ΔHV69/70 deletion corresponded to the B.1.1.7 variant, a subset of eight samples (four from week 21_11 and four from week 21_12), with B.1.1.7 proportion estimates ranging between 35 and 100% were analyzed by NGS sequencing of (Figure 2 ). Seven additional amino acid substitutions/deletions, which were not specific for the B.1.1.7 variant, were also detected: G142V, A222V, G257V, S375P, ΔT376, F377L, and K537E. Of these substitutions, as of April 28, G142V and G257V had already been reported in 1128 and 259 sequences published at the GISAID database, respectively, but the others had only been reported at low frequencies.

Substitutions S375P/ΔT376/F377L affecting three consecutive residues within the receptor-binding domain (RBD) and identified in 35% of sequences from WWTP-31 in Madrid have been individually reported less than 15 times in several countries, including Spain. However, our data report the occurrence of the three mutations within the same variant. Of note is that mutations in these residues have been related to antigenic drift. 24 1.7 in the territory. Total levels of SARS-CoV-2 RNA were determined by RT-qPCR using N1 target as well as the S discriminatory RT-qPCR, without normalization by the population number (Figure 3) . A moderate correlation was observed between both SARS-CoV-2 genome concentration measures between N1 and total S gene titers (R 2 = 0.303, data not shown). From all samples analyzed for different SARS-CoV- Environmental Science & Technology pubs.acs.org/est Article regions, etc.) and were associated to the nationwide state of alarm, in place since October 2020. As most European countries, at the clinical level, a peak in COVID-19 cases occurred between the end of December 2020 (week 20_52) and early February 2021 (week 21_05). As measured through the N1 target, a peak in SARS-CoV-2 genome levels in wastewater was observed in week 21_01 in nine regions including Zaragoza ( Figure 3B ; WWTP-15), Las Palmas in Canary Islands ( Figure 3D ; WWTP-18), three cities in Catalonia ( Figure 3H ; WWTP-26, -27, and -28), and Madrid ( Figure 3I ; WWTP-07, -08, -31, and -32); in week 21_02 in regions including Coŕdoba and Granada in Andalucía ( Figure 3A ; Figure 3N ; WWTP-13), and Oviedo ( Figure 3O ; WWTP-16), and in week 21_05 in Sevilla ( Figure 3A ; WWTP-10) and Vitoria ( Figure 3N ; WWTP-12). SARS-CoV-2 genome titers in Tenerife (Canary Islands) ( Figure 3D ; WWTP-19) were already at peak titers in the first week of the study (Figure 3 ). First detection of the B.1.1.7 variant in wastewater samples occurred in the southern city of Maĺaga (Andalucía) in week 20_52 (Figure 4 ). In the first week of January 2021, it could also be detected in the two largest cities (Madrid and Barcelona), in two northern cities (Santander and Vitoria), another city in Andalucía (Coŕdoba), and in Tenerife (Canary Islands), suggesting that multiple introductions occurred during Christmas holidays in different parts of the Spanish territory and representing 20% of all sampled WWPTs. B.1.1.7 levels detected in week 21_01 were lower than 10% in most WWTPs, with the exception of WWTP-32 in Madrid, which was 22.9%. Percentages of WWTPs with B.1.1.7 detection increased progressively, up to 56% by week 21_04, 91% by week 21_08, and 97% by week 21_13. Our data also showed that while the total level of SARS-CoV-2 RNA genomes measured by the N1 target showed a slight decline for several weeks after the first peak observed during early 2021, this negative trend was slowed or reversed in most WWTPs from the time when the proportion of B.1.1.7 became more abundant (Figure 3 ). In six cities, a significant increase of N1 RNA titers higher than 1 log 10 with respect to the preceding week, as adopted by the VATar COVID-19 Spanish Reporting System (available at https:// www.miteco.gob.es/es/agua/temas/concesiones-y-autorizaciones/nota-tecnica-vatar-miterd_tcm30-517518.pdf), was observed in the last week of the study ( Figure 3E WWTP-20, Figure 3F WWTP-29, Figure 3G WWTP-23, and Figure 3I WWTP-07, -08, and -30).

The relative proportion of the B.1.1.7 variant in wastewater could be estimated for 91% of positive samples. Figure 4 shows the heatmap of the evolution of B.1.1.7 prevalence in wastewater over time. Predominance of the B.1.1.7 variant with prevalences ≥ 90% was reached in all WWTPs except in WWTP-08, -09, and -32 in Madrid, and WWTP-15 in Zaragoza. In Madrid, B.1.1.7 reached 48 and 86% and in Zaragoza 51% at the end of the study period, although it kept progressively increasing thereafter (data not shown). When considering data from nine WWTPs, which showed at least 2 consecutive weeks with B.1.1.7 percentages near fixation (90−100%) as a confirmation of predominance (Table S2 ), we could estimate that approximately 8.1 ± 2.0 weeks' time was required for the B.1.1.7 variant to become predominant in sewage. This would correspond to an average increase rate of 11.7% (9.6−15.1%) per week.

The proportion of B.1.1.7 was compared to the prevalence detected at the clinical level. The abundance of B.1.1.7, as a fraction of all sequenced clinical specimens by the local authorities in each autonomous community, was obtained from update reports of the epidemiological situation of the variants of SARS-CoV-2 of importance published by the Spanish Ministry of Health. 27 When comparing the proportion of B.1.1.7 estimated from wastewater with the proportion estimated from sequencing of clinical isolates, a good correlation was observed (R 2 = 0.5012) ( Figure 5A ). Analysis comparing the B.1.1.7 prevalence estimate with clinical specimens with prevalence estimated from wastewater 1 or 2 weeks before did not increase the correlation determinant (data not shown), suggesting that wastewater analysis did not allow us to anticipate the increase in B.1.1.7 proportion. On a geographical/temporal analysis, wastewater testing allowed the confirmation of circulation of the B.1.1.7 variant before its identification in clinical specimens, especially in part due to a noticeable clinical undertesting in most regions. Data for 3 selected weeks are shown in Figure 5B . By week 21_04 at the end of January, at the clinical level, the variant was only detected on clinical samples in Galicia and País Vasco, mainly due to the low number of sequenced isolates in most regions, while it was detected in 18/32 (56%) WWTPs. At the end of the study, all autonomous communities' public health departments reported percentages higher than 60%, but while some WWTPs showed very high percentages, some others showed relatively lower proportions, indicating differences at the local level.

As a cost-effective approach to screen thousands of inhabitants, wastewater-based epidemiology (WBE) is a valuable tool to anticipate the circulation of specific pathogens in a community and to closely track their incidence evolution through space and time. 28, 29 This surveillance has proven to be extremely useful in monitoring the circulation of total SARS-CoV-2 in different parts of the world 1−3 and may be effective for tracking novel SARS-CoV-2 VOCs. Despite NGS allowing the definitive identification of specific variants and has already been applied on sewage samples, 30, 31 it is resource and time-consuming, thus limiting the number of samples that can be processed and the number of labs which can implement this approach on a regular basis. In addition, most deep sequencing protocols allow the identification of signature mutations within short individual reads, which when detected in wastewater samples containing a mixture of different isolates may not be proof of co-occurrence of such mutations within the same genome, thus providing only an indirect evidence of the presence of a certain variant. 32 Finally, the use of NGS is also challenged by the presence of inhibitors in samples, limiting its success rate and depth coverage. Given its higher tolerance to inhibitors, droplet digital RT-PCR has been acknowledged as a suitable approach to simultaneously enumerate the concentration of variants with the N501Y mutation and wildtype in wastewater, 33 but the widespread use of droplet digital RT-PCR may be limited nowadays because of the high economic investment in instrumentation. On the other hand, the use of RT-qPCR methods offers the advantage of rapid turnaround time, lower Environmental Science & Technology pubs.acs.org/est Article cost, and immediate availability in most public health laboratories. In the current study, we validated a duplex RT-qPCR assay to discriminate and enumerate SARS-CoV-2 variants containing the ΔHV69/70 deletion from variants lacking it. Among molecular markers specific for the B.1.1.7 variant, the 6-nucleotide deletion corresponding to residues 69/ 70 was chosen because it offers the possibility to design highly specific robust probes to be used in wastewater samples, which, unlike clinical specimens, will contain mixed sequences in most cases. Similar to the TaqPath COVID-19 assay (Thermo Fisher Scientific) and other RT-qPCR protocols designed for clinical diagnosis, 34 the novel duplex RT-qPCR assay developed in this study proved highly specific and discriminatory. The ΔHV69/70 deletion is located within the N-terminal domain of the S glycoprotein and has been described to be located at a recurrent deletion region (RDR), and phylogenetic studies showed that it has arisen independently at least 13 times. 35 In addition to being a signature mutation of the highly transmissible B. During the study period, weekly wastewater estimates of the proportion of B.1.1.7, representing a larger and more comprehensive proportion of typed cases including both symptomatic and asymptomatic cases, well reflected the trends in the reported sequenced clinical cases in most regions. Although the number of clinical specimens sequenced from public health laboratories was not high during the study period and showed strong geographic differences, a correlation was observed between the proportion of B.1.1.7 cases observed at the clinical level and the data estimated from sewage when using samples from the same week ( Figure 5A ). Of note is that this association was not more robust when using data from 1−2 previous weeks (data not shown). The lack of anticipation ability could be due to several unknown factors, including differences in shedding and kinetic levels between variants, differences in the proportion of asymptomatic infections, and differences in environmental stability. Sewage surveillance allowed the identification of B.1.1.7 circulation in the Spanish territory in the southern city of Maĺaga before it was confirmed at the clinical level by national public health authorities and allowed us to infer multiple simultaneous introductions during Christmas and New Year's holidays in distant parts of the country (Madrid, Barcelona, Santander, Vitoria, Coŕdoba, and Tenerife). By the end of January 2021, only 13% (2/15) autonomous communities had reported B.1.1.7 clinical cases, while circulation in sewage had been confirmed in 67% (10/15) of them, confirming its use as an early warning approach. Data from 9 WWTPs, which reached B.1.1.7 near fixation rates, defined as higher than 90% for ≥2 consecutive weeks, showed that 8.1 ± 2.0 weeks were required to reach B.1.1.7 predominance, which would be a slightly shorter time than what has been locally observed at the clinical level. A research publication reported first detection of imported B.1.1.7 clinical cases in Madrid in week 20_52 (December 2020) and a proportion of 62% of total newly diagnosed COVID-19 cases 10 weeks later. 40 Data from other studies are also in the UK reported by ECDC show that B.1.1.7 cases went from less than 5% of all positive cases to more than 60% in less than 6 weeks during November to mid-December 2020, 41 and Davies et al. demonstrated that it became dominant throughout the country. 15 Estimates from the US indicate that B.1.1.7 would become dominant in most states 4 months after its first identification in late November 2020. 14 Our data also showed that predominance of the B.1.1.7 variant appeared to correspond to a slowdown in the negative trend of total SARS-CoV-2 wastewater levels, which had been observed from early 2021 in most cities, probably due to lockdown measures and the mass-vaccination campaign, which was initiated in the last week of 2020. Even in some cities, including Santander ( Figure 3E ; WWTP-20), Cuenca ( Figure  3F ; WWTP-29), Valladolid ( Figure 3G ; WWTP-23), and Madrid ( Figure 3I ; WWTP-07, -08, and -30), total SARS-CoV-2 levels showed a positive trend at the end of the study. These results suggest that the emergence of B.1.1.7 cases could have produced a higher transmission rate and a slight increase in COVID-19 incidence, as confirmed by clinical epidemiological data reported by the Spanish Ministry of Health, reporting an incidence peak between the end of March and April 2020. 42 Despite this positive trend markedly observed in some regions, the smooth running of the mass-vaccination campaign starting on December 27, 2020, in addition to nonpharmaceutical interventions, likely contributed to minimizing the impact of B.1.1.7 emergence. As of the end of March 2021, the percentage of the Spanish population who had been partially immunized or totally vaccinated were of 13.2 and 6.8%, respectively.

Finally, despite the state of alarm decreed by the Spanish government, as a measure to unify confinement and restriction measures across the country, was maintained throughout the study period, predominance of the B.1.1.7 variant was not homogeneous, and dynamics were variable among cities across the country. For instance, a rapid predominance of the B.1.1.7 variant was observed in Granada (WWTP-4) and Caćeres (WWTP-36), while in other cities, B.1.1.7 reached prevalences higher than 90% by week 3−6 after first positive detection and decreased thereafter for 2−3 weeks. Several reasons could explain these waves, including differences in regional social distancing behaviors, repetitive B.1.1.7 case imports, introduction of additional variants, climatic effect on sewage composition, size of the WWTP, or variability related to the use of grab samples instead of composite samples.

This study highlights the use of WBE as a cost-effective, noninvasive, and unbiased approach, which may complement clinical testing during the COVID-19 pandemic, and demonstrates the applicability of duplex RT-qPCR assays on sewage surveillance as a rapid, attractive, and resourceful method to Environmental Science & Technology pubs.acs.org/est Article track the early circulation and emergence of the known VOC in a population, especially at times when clinical typing is insufficient and when signature mutations can be unequivocally assigned to a specific VOC. The current strategy could be readily adaptable to track specific mutations of other VOCs as soon as they are identified by clinical genomic sequencing in the future and integrated into existing wastewater surveillance programs.

■ ASSOCIATED CONTENT * sı Supporting Information

The Supporting Information is available free of charge at https://pubs.acs.org/doi/10.1021/acs.est.1c03589.

Additional data on parameters defining standard curves, LOD and LOQ for RT-qPCR assays used in the study (Table S1) , and data used to estimate the time required to reach B.1.1.7 prevalence >90% in wastewater ( 

Implementation of Environmental Surveillance for SARS-CoV-2 Virus to Support Public Health Decisions: Opportunities and Challenges

Wastewater-Based Epidemiology as a Useful Tool to Track SARS-CoV-2 and Support Public Health Policies at Municipal Level in Brazil

COMMISSION RECOMMENDATION of 17.3.2021 on a Common Approach to Establish a Systematic Surveillance of SARS-CoV-2 and Its Variants in Wastewaters in the EU

First Confirmed Detection of SARS-CoV-2 in Untreated Wastewater in Australia: A Proof of Concept for the Wastewater Surveillance of COVID-19 in the Community

First Detection of SARS-CoV-2 in Untreated Wastewaters in Italy

SARS-CoV-2 RNA in Wastewater Anticipated COVID-19 Occurrence in a Low Prevalence Area

Time Evolution of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) in Wastewater during the First Pandemic Wave of COVID-19 in the Metropolitan Area of Barcelona

Evaluation of Lockdown Effect on SARS-CoV-2 Dynamics through Viral Genome Quantification in Waste Water

Making Waves: Wastewater-Based Epidemiology for COVID-19 -Approaches and Challenges for Surveillance and Prediction

Tracking COVID-19 via Sewage

COVID-19 (SARS-CoV-2) Outbreak Monitoring Using Wastewater-Based Epidemiology in Qatar

Investigation of Novel SARS-COV-2 Variant Variant of Concern 202012/01 (PHE),-Public Health England

Detection Of Enteric Viruses

Sola, I. Identification of a Coronavirus Transcription Enhancer

A Hepatitis Virus Complicating Studies With Mouse Leukemia

Cyclooxygenase Activity Is Important for Efficient Replication of Mouse Hepatitis Virus at an Early Stage of Infection

A Real-Time TaqMan RT-PCR Assay with an Internal Amplification Control for Rapid Detection of Transmissible Gastroenteritis Virus in Swine Fecal Samples

Comparison and Evaluation of Conventional RT-PCR, SYBR Green I and TaqMan Real-Time RT-PCR Assays for the Detection of Porcine Epidemic Diarrhea Virus

Microbiology of the Food Chain -Horizontal Method for Determination of Hepatitis A Virus and Norovirus Using Real-Time RT-PCR -Part 1: Method for Quantification

Cleavage Site in the Viral Quasispecies of COVID19 Patients

Disease and Diplomacy: GISAID's Innovative Contribution to Global Health. Global challenges

Complete Mapping of Mutations to the SARS-CoV-2 Spike Receptor-Binding Domain That Escape Antibody Recognition

Hepatitis E in Italy: 5 Years of National Epidemiological, Virological and Environmental Surveillance

Monitoring Human Enteric Viruses in Wastewater and Relevance to Infections Encountered in the Clinical Setting: A One-Year Experiment in Central France

Monitoring SARS-CoV-2 Circulation and Diversity through Community Wastewater Sequencing, the Netherlands and Belgium

Detection of SARS-CoV-2 Variants in Switzerland by Genomic Analysis of Wastewater Samples

Droplet Digital RT-PCR to Detect SARS-CoV-2 Variants of Concern in Wastewater

Network for Genomic Surveillance in South Africa

Multiplex QPCR Discriminates Variants of Concern to Enhance Global Surveillance of SARS-CoV-2

Recurrent Deletions in the SARS-CoV-2 Spike Glycoprotein Drive Antibody Escape

The Circulating SARS-CoV-2 Spike Variant N439K Maintains Fitness While Evading Antibody-Mediated Immunity

Josset, L.; the COVID-Diagnosis HCL Study Group Two-Step Strategy for the Identification of SARS-CoV-2 Variant of Concern 202012/01 and Other Variants with Spike Deletion H69− V70

SARS-CoV-2 Mutations and Variants of

SARS-CoV-2 -Increased Circulation of Variants of Concern and Vaccine Rollout in the EU/EEA, 14th Update −

The authors declare no competing financial interest.Environmental Science & Technology pubs.acs.org/est Article