key: cord-0794753-703kcvhf authors: Calderón-Franco, David; Orschler, Laura; Lackner, Susanne; Agrawal, Shelesh; Weissbrodt, David G. title: Monitoring SARS-CoV-2 in sewage: Toward sentinels with analytical accuracy date: 2021-09-10 journal: Sci Total Environ DOI: 10.1016/j.scitotenv.2021.150244 sha: df1e1e53b9ec3728a8046bb6201477e5b4a0caa4 doc_id: 794753 cord_uid: 703kcvhf The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemia has been one of the most difficult challenges humankind has recently faced. Wastewater-based epidemiology has emerged as a tool for surveillance and mitigation of potential viral outbreaks, circumventing biases introduced by clinical patient testing. Due to the situation urgency, protocols followed for isolating viral RNA from sewage were not adapted for such sample matrices. In parallel to their implementation for fast collection of data to sustain surveillance and mitigation decisions, molecular protocols need to be harmonized to deliver accurate, reproducible, and comparable analytical outputs. Here we studied analytical variabilities linked to viral RNA isolation methods from sewage. Three different influent wastewater volumes were used to assess the effects of filtered volumes (50, 100 or 500 mL) for capturing viral particles. Three different concentration strategies were tested: electronegative membranes, polyethersulfone membranes, and anion-exchange diethylaminoethyl cellulose columns. To compare the number of viral particles, different RNA isolation methods (column-based vs. magnetic beads) were compared. The effect of extra RNA purification steps and different RT-qPCR strategies (one step vs. two-step) were also evaluated. Results showed that the combination of 500 mL filtration volume through electronegative membranes and without multiple RNA purification steps (using column-based RNA purification) using two-step RT-qPCR avoided false negatives when basal viral load in sewage are present and yielded more consistent results during the surveillance done during the second-wave in Delft (The Hague area, The Netherlands). By paving the way for standardization of methods for the sampling, concentration and molecular detection of SARS-CoV-2 viruses from sewage, these findings can help water and health surveillance authorities to use and trust results coming from wastewater based epidemiology studies in order to anticipate SARS-CoV-2 outbreaks. In December 2019, China reported the outbreak of a novel coronavirus named SARS-CoV-2 (genus Betacoronavirus, family Coronaviridae), which has spread around the world fast and lethally. The progression of this SARS-CoV-2 pandemic has been monitored primarily by clinically testing symptomatic individuals for presence of viral RNA. However, as asymptomatic persons seem to account for 40 -45% of the infections (Ooi and Low, 2020) , registering solely symptomatic cases does not allow for prediction of the real cumulative incidence of the viral outbreak. In addition, the number of detected cases highly depends on access to diagnostics as well as the threat of isolation and quarantine dissuading people from getting tested. As SARS-CoV-2 RNA can be found in patients' stool, samples of rectal swabs, urine and other bodily fluids, detection of the virus via analysing wastewater has gained interest (Agrawal et al., 2021; Foladori et al., 2020; Kitajima et al., 2020; Medema et al., 2020b; Randazzo et al., 2020a; Wu et al., 2020) . Wastewater-based epidemiology can serve as a tool to trace the circulation of a virus variant within a community or municipality. It does not only protect the anonymity of communitymembers, but it is also unbiased and focused on real-time information, becoming an earlywarning system for viral surveillance (Ahmed et al., 2020a; Larsen and Wigginton, 2020) . It can help make better-informed decisions in public health policies. The urgency of the pandemics pressed the need for a rapid deployment of available analytical methods to track the fate of the virus and its variants. This resulted in a large palette of different protocols for every step. For instance, sampling volumes used vary from a minimum of 2 mL (Rimoldi et al., 2020) , to 45 -60 mL (Hokajärvi et al., 2021; Westhaus et al., 2021) , to 100-200 mL (Medema et al., 2020b; Tanhaei et al., 2021) , up to 1L (Agrawal et al., 2021) ; concentration methods vary from precipitation (Barril et al., 2021; Pérez-Cataluña et al., 2021; Sapula et al., 2021) , to centrifugation (D'Aoust et al., 2021; Kitamura et al., 2021) , ultracentrifugation (Prado et al., methods (Ambrosi et al., 2021; Pérez-Cataluña et al., 2021) and on improving the viral concentration steps (Parra Guardado et al., 2020) . It has been previously mentioned that one of the variables that should be checked is the filtered volume: low sample volume in combination with low recovery yields could lead to false negative or lower measured SARS-CoV-2 concentration (Alygizakis et al., 2021) . Here, we aimed to establish a method to quantify SARS-CoV-2 from influent wastewater. The effects of filtration volumes of influent wastewater (50, 100 and 500 mL), different concentration methods (polyethersulfone membranes, electronegative membranes and anion-exchange chromatography), different RNA isolation methods (column based vs. magnetic beads), and different RT-qPCR strategies (one-step vs. two-step) were assessed. A methodology designed from the selected best parameter testing outputs was used to survey viral gene copies across the second wave of the pandemic in the Delfland water catchment area (The Hague, The Netherlands) and compare it to publicly available clinical data. We deliver an assessment of method variability and biases introduced when basal viral loads are present and establish a methodology implementable in wastewater-based epidemiology. This will support health and water authorities to anticipate, monitor, and control future outbreaks. A summary of the methodology can be found in Figure 1 . In brief, municipal wastewater was collected to test the effect of 3 filtration volumes (50, 100, 500 mL) and 3 virus concentration methods (polyethersulfone membrane, electronegative membrane, anion exchange column) on molecular analyses of 3 genes (spike S, nucleocapsid N, ORF1ab) of SARS-CoV-2 by one-step and two-step reverse transcription (RT) and quantitative PCR (qPCR) after isolation of viral RNA using single (magnetic beads, columns) and multiple purification methods. The optimal method parameters were then applied to monitor a time series over the second wave of the epidemics. ORF1ab) were quantified by RT-qPCR (one-step vs two-step RT-qPCR) after testing the effect of the filtration volume and different concentration methods. Created with BioRender. in The Netherlands (Figure 2) . The 24-composite sample from July 7 2020 was used for method development. All samples were processed in a timeframe of less than 1 day after collection. Influent wastewater samples were heated at 65°C for 30 minutes to deactivate and handle the sewage sample under biosecure conditions (Batéjat et al., 2021; Rabenau et al., 2005) . The effect of used influent wastewater volume and the concentration method was assessed for viral genes quantification. A 24-h composite sample (July 7 2020) was used, representing the beginning of the second wave ( Figure 2) . The goal was to identify the best concentration method when basal viral levels were present in influent wastewater. A volume of 200 µL (50 ng µL -1 ) of an F-specific RNA MS2 phage (Thermo Fisher Scientific, Germany) was spiked as internal standard in the influent wastewater solution before the concentration methods. The group of F specific RNA phages (f2 and MS2 are the best known members of this group) are considered the most attractive subgroup of phages to serve as indicator organisms because their physical structure closely resembles that of enteroviruses in wastewater (Havelaar et al., 1985) . Volumes of 50, 100 and 500 mL of influent wastewater were tested. All the analyses were performed in three biological replicates. These influent wastewater volumes were filtered (membranes) or loaded (chromatographic column) in three technical replicates with three different methodologies: (i) filtration on polyethersulfone (PES) membrane of 0.22 µm pore size (Pall Corporation, USA); (ii) retention on a 1-mL diethylaminoethyl cellulose (DEAE) anion-exchange chromatography (BIA separations, Slovenia) which is used for phages separations (Kattur Venkatachalam et al., 2014; Kramberger et al., 2010) and for isolations of free-floating extracellular DNA from water matrices (Calderón-Franco et al., 2021) ; and (iii) attraction on electronegative membranes of 0.45 µm pore size (Merck Millipore, The Netherlands). Prior to loading on the DEAE chromatography column, the influent wastewater was filtered through the 0.22 µm PES to remove biosolids and particles, which is necessary to protect the DEAE column. Loading and elution in chromatography is more time consuming than direct filtration. However, chromatography is more selective and leads to higher grades of purity of concentrated analytes. Thus, we also considered the DEAE method as a way to identify the recovery of viral RNA when compared to the filtration methods. Extractions of raw viral RNA from membrane filters and DEAE column eluents were performed by using two different commercial kits for RNA extraction. First, we used the MagMax CORE Nucleic Acid Purification kit (Thermo Fisher Scientific, The Netherlands), that employ magnetic beads for nucleic acid purification, together with the MagMax CORE Mechanical Lysis Module (Thermo Fisher Scientific, The Netherlands) by following manufacturer's recommendations. Second, we used the Fast RNA Blue kit (MP Biomedicals), which is based on column-based purification (solid phase extraction), according to the manufacturer's protocol. In addition, the effect of multiple RNA purification steps before quantification of genes by RT-qPCR was assessed by involving the GeneJET NGS Cleanup Kit (Thermo Scientific, USA) according to manufacturer's instructions obtaining the purified RNA sample. The quality and quantity of raw and purified RNA extracts were measured with a NanoDrop Spectrophotometer (ND-1000, USA). Journal Pre-proof On the lowest basal viral load day from our sampling campaign (7th July 2020), we compared one-step and two-step RT-qPCRs in order to assess the effect of quantification of viral RNA genes when the incidence of cases is low. This was assessed to avoid potential false negatives. One-step RT-qPCR was performed using the isolated RNA as a template (raw RNA), which included specific primers for the target genes. Two-step RT-qPCR was performed by synthesizing complementary DNA (cDNA) from the pool of raw RNA isolated per sample using the SuperScript VILO cDNA synthesis kit (Thermo Fisher Scientific, The Netherlands) following manufacturer's recommendations, included random primers, prior to using the cDNA as a template for qPCR. The concentration of the synthesized cDNA was measured with Qubit® dsDNA assays (Thermo Fisher Scientific, USA). We quantified the concentration of SARS-CoV-2 RNA in wastewater samples by measuring marker genes using the TaqPath COVID-19 RT-PCR Kit (Thermo Fisher Scientific, Germany) with a QuantStudio 3 Thermal Cycler, according to a previous study (Agrawal et al., 2021) . The TaqPath COVID-19 RT-PCR Kit includes primer pairs targeting genes that code for the structural transmembrane spike (S) glycoprotein and nucleocapsid (N) protein, and for the non-structural open reading frame of the replicase complex (ORF1ab) and that were used in a multiplex assay ( Table S1 ). The kit includes: (1) TaqPath COVID-19 Assay Multiplex, which contains three primer/probe sets specific to different SARS-CoV-2 genomic regions (i.e. N gene, S gene and Orf1ab gene) and primers/probes for bacteriophage MS2. (2) MS2 Phage Control -RNA control, having a concentration of 10^6 copies per µl, to verify the efficacy of the RNA extraction and the absence of inhibitors in the PCR reaction. (3) TaqPath COVID-19 Control -Positive SARS-CoV-2 RNA control that contains targets specific to the SARS-CoV-2 genomic regions targeted by the assays. The manufacturer (Thermo Fisher Scientific) has not publicly released the primers/probe sets sequences, therefore, we do not have access to the information related to the primers/probe sets sequences and the length of the PCR products. In case of the positive control, we included triplicates of the four different concentration (i.e. (Table S2) . Primer efficiencies were 100.17 ± 5.98% for N, 103.71 ± 11.12% for S, 103.08 ± 5.42% for Orf1ab, and 86.42 ± 6.46% for MS2 phage (n = 3 runs, mean ± SD). The slopes of the standard curves for the quantification were −3.33 ± 0.14 for N, −3.22 ± 0.21 for S, −3.27 ± 0.14 for Orf1ab, and −3.67 ± 0.19 for MS2 phage. The SARS-CoV-2 concentrations detected in the samples are presented without correcting for recovery efficiencies. Reactions were considered positive if the cycle threshold (CT) was below 40 cycles, otherwise negative (i.e., no detection of the SARS-CoV-2 RNA in the sample). The lower limit of detection was 10 gene copies per RT-qPCR reaction from the positive control. Statistical analyses were performed on all molecular datasets with R 3.5.1 (R Foundation for Statistical Computing., 2018) and RStudio (https://www.rstudio.com/). The RT-qPCR abundance data were analyzed in R using ggplot2 (v0.9.3.1). J o u r n a l P r e -p r o o f The effect of different initial influent wastewater volumes (50, 100 and 500 mL), different viral concentration methods, different RNA purification methods, and different one-step/two-step RT-qPCR methods were assessed on a basal viral level timepoint (July 7 2020), corresponding to the beginning of the second wave (Figure 2) and to the lowest day in terms of cumulative incidence of COVID-19 cases. This timepoint served as reference for method standardization as basal viral load levels from influent wastewater samples. The selected methodology was further used for surveilling SARS-CoV-2 gene copies over the second wave (July to December 2020). Among the three different methods tested to concentrate virus particles (electronegative membranes, PES membranes, DEAE column), concentration by electronegative membranes of 0.44 µm pore size was the only one displaying positive results in terms of detection of S, N, ORF1ab genes and MS2 phage control (Figure 3) . When PES membranes were used, SARS-CoV-2 was not detected from the cake, suggesting that there was no retention to the membrane. However, internal control MS2 phage, which was added prior to filtration, was detected in PES membrane. This difference may result from structural viral SARS-CoV-2, among other viruses, has an isoelectric point (pI = ½ (pKa 1 + pKa 2 ) below 7. The pH of the influent wastewater pH ranged between 6.3 -6.8. At pH values below the isoelectric point, the charge of SARS-CoV-2 and MS2 phage could be positive, thus explaining why they J o u r n a l P r e -p r o o f Journal Pre-proof did not bind to the DEAE (positive net charge) column and did bind to the electronegative membranes (Joonaki et al., 2020) . When expressing the RT-qPCR results in concentrations of viral gene copies detected per L, filtering 500 mL influent water through electronegative membranes gave, in average, higher values of gene copies per L (7.2⋅10 4 ± 1.7⋅10 4 viral gene copies L -1 ) when compared to gene copies obtained from 50 mL (7.4⋅10 3 ± 4.1⋅10 3 viral gene copies L -1 ) or 100 mL (3.0⋅10 4 ± 1.9⋅10 4 viral gene copies L -1 ). In Figure 3 , individual gene results were also more consistent when 500 mL were filtered (S gene: 5.6⋅10 4 ± 2.8⋅10 4 viral gene copies L -1 ; N gene: 6.0⋅10 4 ± 2.6⋅10 3 viral gene copies L -1 ; ORF1ab gene: 7.5⋅10 4 ± 2.4⋅10 4 viral gene copies L -1 ) when compared to filtered 50 mL (S gene: 1.3⋅10 4 ± 1.0⋅10 4 viral gene copies L -1 ; N gene: 4.2⋅10 3 ± 2.3⋅10 3 viral gene copies L -1 ; ORF1ab gene: 6.9⋅10 3 ± 2.9⋅10 3 viral gene copies L -1 ) and 100 mL (S gene: 5.2⋅10 4 ± 1.6⋅10 4 viral gene copies L -1 ; N gene: 9.7⋅10 3 ± 6.7⋅10 3 viral gene copies L -1 ; ORF1ab gene: 2.9⋅10 4 ± 4.7⋅10 3 viral gene copies L -1 ). Here we demonstrated that filtration of larger volumes (i.e., 500 mL) in combination of electronegative membranes allowed, in the lowest day of COVID-19 incidence of the second wave, to adsorb and concentrate enough SARS-CoV-2 for highest and most consistent (uniform results among gene tested) detection and quantification by two-step RT-qPCR. Journal Pre-proof One of the major analytical concerns related to the wastewater-based epidemiology of SARS-CoV-2 is the quality and purity of the isolated RNA, especially when low viral loads are present. This can potentially compromise the accuracy of the results obtained through false negatives. The effect of multiple RNA purification steps and one-step/two-step RT-qPCR strategies were assessed. The best concentration method obtained during the first phase of this study (i.e., 500 mL influent wastewater filtered on electronegative membranes) was utilized. gene assays in gene copies per 500 mL, using one-step (raw RNA and purified RNA) and two-step (cDNA-raw RNA and cDNA-purified RNA) RT-qPCR. In Figure 4 , it can be observed that for the three genes tested (N, ORF1ab and S), the two-step RT-qPCR performed from raw and purified RNA displayed, on average, higher values (1.6⋅10 5 ± 1.0⋅10 4 viral gene copies 500 mL -1 ) than one-step RT-qPCR (no positive results from tested genes). No detection was obtained from one-step RT-qPCR both from raw RNA and purified al., 2012) . This can affect the sensitivity of the assay and result in false-negative results (Ahmed et al., 2020b) . The three genes were consistently detected by two-step RT-qPCRfrom both raw RNA and purified RNA. However, RNA purification induced higher variability on the average of all the tested genes (1.7⋅10 5 ± 1.1⋅10 5 gene copies in 500 mL, 65% relative standard deviation (RSD)) when compared to the lower standard deviation obtained with two-step RT-qPCR (1.5⋅10 5 ± 3.7⋅10 4 gene copies in 500 mL, 25% RSD). When individual genes were tested, the N gene exhibited the highest variability (2-fold) between raw (1.2⋅10 5 ± 8.0⋅10 4 gene copies in 500 mL) and purified (7.3⋅10 4 ± 6.0⋅10 4 gene copies in 500 mL) RNA. Thus, multiple RNA purification stages before RT-qPCR are not recommended: RNA from viral particles were lost during the purification process even if detected with two-step RT-qPCR. The same trend was observed during the time course of the 6-months sampling campaign, where two-step RT-qPCR measurements displayed higher and more consistent results (i.e., uniform values between the three tested genes) ( Figure S1 ). Our observation about how relevant it is to consider one-step vs. two-step RT-qPCR when viral load in wastewater samples is very low is in consensus with a previous study (Chik et al., 2021) . This was of high importance as it was the day, from the sampling campaign, with the lowest viral (2) and run a two-step RT-qPCR analysis, which allow enhancing the resolution. This is because cDNA is firstly synthesized from total viral RNA using random primers followed by specific primers for qPCR analysis, being able to detect low-levels of SARS-CoV-2 RNA copies in sewage. Journal Pre-proof The range of commercial options to isolate RNA and quantify SARS-CoV-2 viral particles is diverse and complex. Two different RNA isolation methods were chosen to survey the number of viral particles over the sampling campaign ( Figure 5) : one column-based microbial RNA extraction kit (Fast RNA) and one magnetic beads-based kit (MagMax Core kit). Comparative analysis between RNA extraction methods are displayed in Figure 5 . Both methods showed similar results in terms of gene copies per liter (no significant difference; p>0.05) during the whole sampling campaign with the exception of October 20, 2020 (p<0.05), which happened to be the peak of the second wave. The column-based method resulted in less variation between replicates when compared to the magnetic beads method. Our results can be compared with the results from a previous measuring campaign performed in The Hague (Medema et al., 2020a) , where results were reported as concentrations of N gene copies per mL on the first wave of the pandemics from February 2020 until August 2020. The results of July and August were in the range of 10 4 -10 5 viral N gene copies L -1 . The concentrations obtained with our optimized methodology are 10-fold higher (10 5 -10 6 gene copies L -1 ) taking into account the average of these three genes (N, S and ORF1ab) than those obtained by Medema et al. (2020a) . It allowed a more sensitive quantification of viral gene copies at the beginning of the pandemic second-wave. Early detection of SARS-CoV-2 and other future possible viral agents from wastewater requires accurate and adapted-to-matrix methodologies. Thus, being able to quantify low viral levels with non-clinical samples is an important factor to consider in wastewater-based epidemiology, as early communication to the competent authorities is necessary to have enough reaction time to prevent viral outbreaks in urban settlements. The second peak started to increase in the beginning of September 2020 (relative increase of 76.8% when compared the numbers from July to September 2020) while significant viral copies could already be quantified by August 18, 2020 (1.73⋅10 6 viral gene copies L -1 ). This is in accordance with Larsen and Wigginton (2020), who have claimed that wastewater-based epidemiology remediates clinical biases that result from low access to isolation diagnostics and quarantine threat dissuading citizens from getting tested. Randazzo et al. (2020b) have highlighted the importance of wastewater-based epidemiology as an early indicator, revealing that members of the community were shedding SARS-CoV-2 RNA particles in their stool even before the first cases had been reported by the authorities in Murcia (Spain). Overall, a global effort is necessary to generate, harmonize, and use standardized analytical protocols for an accurate, reproducible, and comparable detection and quantification of SARS-CoV-2 (Bivins et al., 2020) . On a wastewater-based epidemiology basis, crucial considerations are needed from sampling to preparation and measurements in order to define (i) the sampling points, (ii) the sample size and frequency, (iii) the method to efficiently concentrate viral particles, (iv) the quantification methods used, (v) the primer set(s), (vi) the controls, and (vii) the normalization methods. Some initiatives on protocol standardization have been taken such as the NORMAN-SCORE joint initiative to facilitate data comparison between "SARS-CoV-2 in sewages'' studies (Lundy et al., 2021) . Knowing the inherent biases associated with qPCR analysis of wastewater samples due to the variability in the matrix of the samples (Agrawal et al., 2021; Chik et al., 2021) and the urgency of time, we believe that it is currently important to focus on the robustness and reproducibility of the qPCR method used in each respective lab, especially during low incidence values. Analytical robustness should help translate risk in the form of viral concentrations per volume filtered to react fast enough to minimize undesirable infective and spreading consequences. Molecular sentinels require analytical accuracy for inter-comparison and information delivery. From this work, we conclude that: 1. The higher the influent wastewater sample volume was (500 mL instead of 100 or 50 mL), the more uniform the concentrations of the three viral genes (N, S, ORF1ab) were, especially when low incidence were reported 2. Employing a two-step RT-qPCR method would help to enhance the resolution for samples with low viral loads. 3. Electronegative membranes were optimal and cheap for the concentration of SARS-CoV-2 from pre-settled influent wastewater. 4. Multiple RNA purification steps are not recommended, due to plausible loss of SARS-CoV-2 RNA and thereby, affecting qPCR analysis outcomes. While comparing the column-based vs. magnetic beads method for the RNA purification, we found that the column-based purification method displayed less variability amongst the biological triplicates throughout the sampling campaign than magnetic beads. The authors declare no conflict of interest. DCF, SA and DGW designed the study. DCF, SA and LO performed the experimental investigations. DCF wrote the manuscript with direct contribution, edits, and critical feedback by all authors. Research Council (NWO). The manuscript has been released as a Pre-Print at medRxiv: https://doi.org/10. 1101/2021.06.30.21259751 J o u r n a l P r e -p r o o f Declaration of interests ☒ The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper. J o u r n a l P r e -p r o o f Long-term monitoring of SARS-CoV-2 RNA in wastewater of the Frankfurt metropolitan area in Southern Germany 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 Surveillance of SARS-CoV-2 RNA in wastewater: Methods optimization and quality control are crucial for generating reliable public health information Analytical methodologies for the detection of SARS-CoV-2 in wastewater: Protocols and future perspectives SARS-CoV-2: Comparative analysis of different RNA extraction methods Evaluation of viral concentration methods for SARS-CoV-2 recovery from wastewaters Heat inactivation of the severe acute respiratory syndrome coronavirus 2 Comparison of approaches to quantify SARS-CoV-2 in wastewater using RT-qPCR: Results and implications from a collaborative inter-laboratory study in Canada Quantitative analysis of SARS-CoV-2 RNA from wastewater solids in communities with low COVID-19 incidence and prevalence SARS-CoV-2 from faeces to wastewater treatment: What do we know? A review SARS-CoV-2 RNA F specific RNA bacteriophages in sewage: Methodology and occurrence The detection and stability of the SARS-CoV-2 RNA biomarkers in wastewater influent in Helsinki Surface Chemistry Can Unlock Drivers of Surface Stability of SARS-CoV-2 in a Variety of Environmental Conditions Concentration and purification of enterovirus 71 using a weak anion-exchange monolithic column Making Waves: Collaboration in the time of SARS-CoV-2 -rapid development of an international cooperation and wastewater surveillance database to support public health decision-making 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 Asymptomatic SARS-CoV-2 infection Development and optimization of a new method for direct extraction of SARS-CoV-2 RNA from municipal wastewater using magnetic beads Comparing analytical methods to detect SARS-CoV-2 in wastewater Preliminary results of SARS-CoV-2 detection in sewerage system in niterói municipality Stability and inactivation of SARS coronavirus Metropolitan wastewater analysis for COVID-19 epidemiological surveillance SARS-CoV-2 RNA in wastewater anticipated COVID-19 occurrence in a low prevalence area Presence and infectivity of SARS-CoV-2 virus in wastewaters and rivers Humidity and deposition solution play a critical role in virus inactivation by heat treatment on N95 respirators An optimized and robust PEG precipitation method for detection of SARS-CoV-2 in wastewater PCR inhibitors -occurrence, properties and removal The first detection of SARS-CoV-2 RNA in the wastewater of Tehran Detection of SARS-CoV-2 in raw and treated wastewater in Germany -Suitability for COVID-19 surveillance and potential transmission risks Sewage, Salt, Silica, and SARS-CoV-2 (4S): An Economical Kit-Free Method for Direct Capture of SARS-CoV-2 RNA from Wastewater SARS-CoV-2 Titers in Wastewater Are Higher than Expected from Clinically Confirmed Cases Evaluation of lockdown effect on SARS-CoV-2 dynamics through viral genome quantification in waste water We are grateful to Mariska Ronteltap from the Delfland Water Authority and Abdelmaoula Idrissi from Delfluent Services and WWTP Harnashpolder, The Netherlands, for helping us arrange this